Atlas of Circadian Metabolism Reveals System-wide Coordination and Communication between Clocks.

Metabolic diseases are often characterized by circadian misalignment in different tissues, yet how altered coordination and communication among tissue clocks relate to specific pathogenic mechanisms remains largely unknown. Applying an integrated systems biology approach, we performed 24-hr metabolomics profiling of eight mouse tissues simultaneously. We present a temporal and spatial atlas of circadian metabolism in the context of systemic energy balance and under chronic nutrient stress (high-fat diet [HFD]). Comparative analysis reveals how the repertoires of tissue metabolism are linked and gated to specific temporal windows and how this highly specialized communication and coherence among tissue clocks is rewired by nutrient challenge. Overall, we illustrate how dynamic metabolic relationships can be reconstructed across time and space and how integration of circadian metabolomics data from multiple tissues can improve our understanding of health and disease.

Partitioning circadian transcription by SIRT6 leads to segregated control of cellular metabolism.

Circadian rhythms are intimately linked to cellular metabolism. Specifically, the NAD(+)-dependent deacetylase SIRT1, the founding member of the sirtuin family, contributes to clock function. Whereas SIRT1 exhibits diversity in deacetylation targets and subcellular localization, SIRT6 is the only constitutively chromatin-associated sirtuin and is prominently present at transcriptionally active genomic loci. Comparison of the hepatic circadian transcriptomes reveals that SIRT6 and SIRT1 separately control transcriptional specificity and therefore define distinctly partitioned classes of circadian genes. SIRT6 interacts with CLOCK:BMAL1 and, differently from SIRT1, governs their chromatin recruitment to circadian gene promoters. Moreover, SIRT6 controls circadian chromatin recruitment of SREBP-1, resulting in the cyclic regulation of genes implicated in fatty acid and cholesterol metabolism. This mechanism parallels a phenotypic disruption in fatty acid metabolism in SIRT6 null mice as revealed by circadian metabolome analyses. Thus, genomic partitioning by two independent sirtuins contributes to differential control of circadian metabolism.

Personalized RNA neoantigen vaccines stimulate T cells in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is lethal in 88% of patients1, yet harbours mutation-derived T cell neoantigens that are suitable for vaccines 2,3. Here in a phase I trial of adjuvant autogene cevumeran, an individualized neoantigen vaccine based on uridine mRNA-lipoplex nanoparticles, we synthesized mRNA neoantigen vaccines in real time from surgically resected PDAC tumours. After surgery, we sequentially administered atezolizumab (an anti-PD-L1 immunotherapy), autogene cevumeran (a maximum of 20 neoantigens per patient) and a modified version of a four-drug chemotherapy regimen (mFOLFIRINOX, comprising folinic acid, fluorouracil, irinotecan and oxaliplatin). The end points included vaccine-induced neoantigen-specific T cells by high-threshold assays, 18-month recurrence-free survival and oncologic feasibility. We treated 16 patients with atezolizumab and autogene cevumeran, then 15 patients with mFOLFIRINOX. Autogene cevumeran was administered within 3 days of benchmarked times, was tolerable and induced de novo high-magnitude neoantigen-specific T cells in 8 out of 16 patients, with half targeting more than one vaccine neoantigen. Using a new mathematical strategy to track T cell clones (CloneTrack) and functional assays, we found that vaccine-expanded T cells comprised up to 10% of all blood T cells, re-expanded with a vaccine booster and included long-lived polyfunctional neoantigen-specific effector CD8+ T cells. At 18-month median follow-up, patients with vaccine-expanded T cells (responders) had a longer median recurrence-free survival (not reached) compared with patients without vaccine-expanded T cells (non-responders; 13.4 months, P = 0.003). Differences in the immune fitness of the patients did not confound this correlation, as responders and non-responders mounted equivalent immunity to a concurrent unrelated mRNA vaccine against SARS-CoV-2. Thus, adjuvant atezolizumab, autogene cevumeran and mFOLFIRINOX induces substantial T cell activity that may correlate with delayed PDAC recurrence.

Reprogramming of the circadian clock by nutritional challenge.

Circadian rhythms and cellular metabolism are intimately linked. Here, we reveal that a high-fat diet (HFD) generates a profound reorganization of specific metabolic pathways, leading to widespread remodeling of the liver clock. Strikingly, in addition to disrupting the normal circadian cycle, HFD causes an unexpectedly large-scale genesis of de novo oscillating transcripts, resulting in reorganization of the coordinated oscillations between coherent transcripts and metabolites. The mechanisms underlying this reprogramming involve both the impairment of CLOCK:BMAL1 chromatin recruitment and a pronounced cyclic activation of surrogate pathways through the transcriptional regulator PPARγ. Finally, we demonstrate that it is specifically the nutritional challenge, and not the development of obesity, that causes the reprogramming of the clock and that the effects of the diet on the clock are reversible.

Single-cell genomic variation induced by mutational processes in cancer.

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Ovarian cancer mutational processes drive site-specific immune evasion.

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability1-4 patterned by distinct mutational processes5,6, tumour heterogeneity7-9 and intraperitoneal spread7,8,10. Immunotherapies have had limited efficacy in HGSOC11-13, highlighting an unmet need to assess how mutational processes and the anatomical sites of tumour foci determine the immunological states of the tumour microenvironment. Here we carried out an integrative analysis of whole-genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumour sites from 42 treatment-naive patients with HGSOC. Homologous recombination-deficient HRD-Dup (BRCA1 mutant-like) and HRD-Del (BRCA2 mutant-like) tumours harboured inflammatory signalling and ongoing immunoediting, reflected in loss of HLA diversity and tumour infiltration with highly differentiated dysfunctional CD8+ T cells. By contrast, foldback-inversion-bearing tumours exhibited elevated immunosuppressive TGFß signalling and immune exclusion, with predominantly naive/stem-like and memory T cells. Phenotypic state associations were specific to anatomical sites, highlighting compositional, topological and functional differences between adnexal tumours and distal peritoneal foci. Our findings implicate anatomical sites and mutational processes as determinants of evolutionary phenotypic divergence and immune resistance mechanisms in HGSOC. Our study provides a multi-omic cellular phenotype data substrate from which to develop and interpret future personalized immunotherapeutic approaches and early detection research.

Clonal fitness inferred from time-series modelling of single-cell cancer genomes.

Progress in defining genomic fitness landscapes in cancer, especially those defined by copy number alterations (CNAs), has been impeded by lack of time-series single-cell sampling of polyclonal populations and temporal statistical models1-7. Here we generated 42,000 genomes from multi-year time-series single-cell whole-genome sequencing of breast epithelium and primary triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), revealing the nature of CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy. Using a new Wright-Fisher population genetics model8,9 to infer clonal fitness, we found that TP53 mutation alters the fitness landscape, reproducibly distributing fitness over a larger number of clones associated with distinct CNAs. Furthermore, in TNBC PDX models with mutated TP53, inferred fitness coefficients from CNA-based genotypes accurately forecast experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDXs resulted in cisplatin-resistant clones emerging from low-fitness phylogenetic lineages in the untreated setting. Conversely, high-fitness clones from treatment-naive controls were eradicated, signalling an inversion of the fitness landscape. Finally, upon release of drug, selection pressure dynamics were reversed, indicating a fitness cost of treatment resistance. Together, our findings define clonal fitness linked to both CNA and therapeutic resistance in polyclonal tumours.

Probabilistic cell-type assignment of single-cell RNA-seq for tumor microenvironment profiling.

Single-cell RNA sequencing has enabled the decomposition of complex tissues into functionally distinct cell types. Often, investigators wish to assign cells to cell types through unsupervised clustering followed by manual annotation or via 'mapping' to existing data. However, manual interpretation scales poorly to large datasets, mapping approaches require purified or pre-annotated data and both are prone to batch effects. To overcome these issues, we present CellAssign, a probabilistic model that leverages prior knowledge of cell-type marker genes to annotate single-cell RNA sequencing data into predefined or de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while controlling for batch and sample effects. We demonstrate the advantages of CellAssign through extensive simulations and analysis of tumor microenvironment composition in high-grade serous ovarian cancer and follicular lymphoma.

Distinct metabolic adaptation of liver circadian pathways to acute and chronic patterns of alcohol intake.

Binge drinking and chronic exposure to ethanol contribute to alcoholic liver diseases (ALDs). A potential link between ALDs and circadian disruption has been observed, though how different patterns of alcohol consumption differentially impact hepatic circadian metabolism remains virtually unexplored. Using acute versus chronic ethanol feeding, we reveal differential reprogramming of the circadian transcriptome in the liver. Specifically, rewiring of diurnal SREBP transcriptional pathway leads to distinct hepatic signatures in acetyl-CoA metabolism that are translated into the subcellular patterns of protein acetylation. Thus, distinct drinking patterns of alcohol dictate differential adaptation of hepatic circadian metabolism.

CircadiOmics: circadian omic web portal.

Circadian rhythms play a fundamental role at all levels of biological organization. Understanding the mechanisms and implications of circadian oscillations continues to be the focus of intense research. However, there has been no comprehensive and integrated way for accessing and mining all circadian omic datasets. The latest release of CircadiOmics (http://circadiomics.ics.uci.edu) fills this gap for providing the most comprehensive web server for studying circadian data. The newly updated version contains high-throughput 227 omic datasets corresponding to over 74 million measurements sampled over 24 h cycles. Users can visualize and compare oscillatory trajectories across species, tissues and conditions. Periodicity statistics (e.g. period, amplitude, phase, P-value, q-value etc.) obtained from BIO_CYCLE and other methods are provided for all samples in the repository and can easily be downloaded in the form of publication-ready figures and tables. New features and substantial improvements in performance and data volume make CircadiOmics a powerful web portal for integrated analysis of circadian omic data.

Comparative Circadian Metabolomics Reveal Differential Effects of Nutritional Challenge in the Serum and Liver.

Diagnosis and therapeutic interventions in pathological conditions rely upon clinical monitoring of key metabolites in the serum. Recent studies show that a wide range of metabolic pathways are controlled by circadian rhythms whose oscillation is affected by nutritional challenges, underscoring the importance of assessing a temporal window for clinical testing and thereby questioning the accuracy of the reading of critical pathological markers in circulation. We have been interested in studying the communication between peripheral tissues under metabolic homeostasis perturbation. Here we present a comparative circadian metabolomic analysis on serum and liver in mice under high fat diet. Our data reveal that the nutritional challenge induces a loss of serum metabolite rhythmicity compared with liver, indicating a circadian misalignment between the tissues analyzed. Importantly, our results show that the levels of serum metabolites do not reflect the circadian liver metabolic signature or the effect of nutritional challenge. This notion reveals the possibility that misleading reads of metabolites in circulation may result in misdiagnosis and improper treatments. Our findings also demonstrate a tissue-specific and time-dependent disruption of metabolic homeostasis in response to altered nutrition.

What time is it? Deep learning approaches for circadian rhythms.

MOTIVATION: Circadian rhythms date back to the origins of life, are found in virtually every species and every cell, and play fundamental roles in functions ranging from metabolism to cognition. Modern high-throughput technologies allow the measurement of concentrations of transcripts, metabolites and other species along the circadian cycle creating novel computational challenges and opportunities, including the problems of inferring whether a given species oscillate in circadian fashion or not, and inferring the time at which a set of measurements was taken. RESULTS: We first curate several large synthetic and biological time series datasets containing labels for both periodic and aperiodic signals. We then use deep learning methods to develop and train BIO_CYCLE, a system to robustly estimate which signals are periodic in high-throughput circadian experiments, producing estimates of amplitudes, periods, phases, as well as several statistical significance measures. Using the curated data, BIO_CYCLE is compared to other approaches and shown to achieve state-of-the-art performance across multiple metrics. We then use deep learning methods to develop and train BIO_CLOCK to robustly estimate the time at which a particular single-time-point transcriptomic experiment was carried. In most cases, BIO_CLOCK can reliably predict time, within approximately 1 h, using the expression levels of only a small number of core clock genes. BIO_CLOCK is shown to work reasonably well across tissue types, and often with only small degradation across conditions. BIO_CLOCK is used to annotate most mouse experiments found in the GEO database with an inferred time stamp. AVAILABILITY AND IMPLEMENTATION: All data and software are publicly available on the CircadiOmics web portal: circadiomics.igb.uci.edu/ CONTACTS: fagostin@uci.edu or pfbaldi@uci.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The pervasiveness and plasticity of circadian oscillations: the coupled circadian-oscillators framework.

MOTIVATION: Circadian oscillations have been observed in animals, plants, fungi and cyanobacteria and play a fundamental role in coordinating the homeostasis and behavior of biological systems. Genetically encoded molecular clocks found in nearly every cell, based on negative transcription/translation feedback loops and involving only a dozen genes, play a central role in maintaining these oscillations. However, high-throughput gene expression experiments reveal that in a typical tissue, a much larger fraction ([Formula: see text]) of all transcripts oscillate with the day-night cycle and the oscillating species vary with tissue type suggesting that perhaps a much larger fraction of all transcripts, and perhaps also other molecular species, may bear the potential for circadian oscillations. RESULTS: To better quantify the pervasiveness and plasticity of circadian oscillations, we conduct the first large-scale analysis aggregating the results of 18 circadian transcriptomic studies and 10 circadian metabolomic studies conducted in mice using different tissues and under different conditions. We find that over half of protein coding genes in the cell can produce transcripts that are circadian in at least one set of conditions and similarly for measured metabolites. Genetic or environmental perturbations can disrupt existing oscillations by changing their amplitudes and phases, suppressing them or giving rise to novel circadian oscillations. The oscillating species and their oscillations provide a characteristic signature of the physiological state of the corresponding cell/tissue. Molecular networks comprise many oscillator loops that have been sculpted by evolution over two trillion day-night cycles to have intrinsic circadian frequency. These oscillating loops are coupled by shared nodes in a large network of coupled circadian oscillators where the clock genes form a major hub. Cells can program and re-program their circadian repertoire through epigenetic and other mechanisms. AVAILABILITY AND IMPLEMENTATION: High-resolution and tissue/condition specific circadian data and networks available at http://circadiomics.igb.uci.edu. CONTACT: pfbaldi@ics.uci.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Circadian control of fatty acid elongation by SIRT1 protein-mediated deacetylation of acetyl-coenzyme A synthetase 1.

The circadian clock regulates a wide range of physiological and metabolic processes, and its disruption leads to metabolic disorders such as diabetes and obesity. Accumulating evidence reveals that the circadian clock regulates levels of metabolites that, in turn, may regulate the clock. Here we demonstrate that the circadian clock regulates the intracellular levels of acetyl-CoA by modulating the enzymatic activity of acetyl-CoA Synthetase 1 (AceCS1). Acetylation of AceCS1 controls the activity of the enzyme. We show that acetylation of AceCS1 is cyclic and that its rhythmicity requires a functional circadian clock and the NAD(+)-dependent deacetylase SIRT1. Cyclic acetylation of AceCS1 contributes to the rhythmicity of acetyl-CoA levels both in vivo and in cultured cells. Down-regulation of AceCS1 causes a significant decrease in the cellular acetyl-CoA pool, leading to reduction in circadian changes in fatty acid elongation. Thus, a nontranscriptional, enzymatic loop is governed by the circadian clock to control acetyl-CoA levels and fatty acid synthesis.

Ongoing genome doubling promotes evolvability and immune dysregulation in ovarian cancer.

Whole-genome doubling (WGD) is a critical driver of tumor development and is linked to drug resistance and metastasis in solid malignancies. Here, we demonstrate that WGD is an ongoing mutational process in tumor evolution. Using single-cell whole-genome sequencing, we measured and modeled how WGD events are distributed across cellular populations within tumors and associated WGD dynamics with properties of genome diversification and phenotypic consequences of innate immunity. We studied WGD evolution in 65 high-grade serous ovarian cancer (HGSOC) tissue samples from 40 patients, yielding 29,481 tumor cell genomes. We found near-ubiquitous evidence of WGD as an ongoing mutational process promoting cell-cell diversity, high rates of chromosomal missegregation, and consequent micronucleation. Using a novel mutation-based WGD timing method, doubleTime , we delineated specific modes by which WGD can drive tumor evolution: (i) unitary evolutionary origin followed by significant diversification, (ii) independent WGD events on a pre-existing background of copy number diversity, and (iii) evolutionarily late clonal expansions of WGD populations. Additionally, through integrated single-cell RNA sequencing and high-resolution immunofluorescence microscopy, we found that inflammatory signaling and cGAS-STING pathway activation result from ongoing chromosomal instability and are restricted to tumors that remain predominantly diploid. This contrasted with predominantly WGD tumors, which exhibited significant quiescent and immunosuppressive phenotypic states. Together, these findings establish WGD as an evolutionarily 'active' mutational process that promotes evolvability and dysregulated immunity in late stage ovarian cancer.

Identification of transcriptional programs using dense vector representations defined by mutual information with GeneVector.

Deciphering individual cell phenotypes from cell-specific transcriptional processes requires high dimensional single cell RNA sequencing. However, current dimensionality reduction methods aggregate sparse gene information across cells, without directly measuring the relationships that exist between genes. By performing dimensionality reduction with respect to gene co-expression, low-dimensional features can model these gene-specific relationships and leverage shared signal to overcome sparsity. We describe GeneVector, a scalable framework for dimensionality reduction implemented as a vector space model using mutual information between gene expression. Unlike other methods, including principal component analysis and variational autoencoders, GeneVector uses latent space arithmetic in a lower dimensional gene embedding to identify transcriptional programs and classify cell types. In this work, we show in four single cell RNA-seq datasets that GeneVector was able to capture phenotype-specific pathways, perform batch effect correction, interactively annotate cell types, and identify pathway variation with treatment over time.

Immune profiling after allogeneic hematopoietic cell transplantation in pediatric acute myeloid leukemia.

Although allogeneic hematopoietic cell transplant (allo-HCT) is curative for high-risk pediatric acute myeloid leukemia (AML), disease relapse remains the primary cause of posttransplant mortality. To identify pressures imposed by allo-HCT on AML cells that escape the graft-versus-leukemia effect, we evaluated immune signatures at diagnosis and posttransplant relapse in bone marrow samples from 4 pediatric patients using a multimodal single-cell proteogenomic approach. Downregulation of major histocompatibility complex class II expression was most profound in progenitor-like blasts and accompanied by correlative changes in transcriptional regulation. Dysfunction of activated natural killer cells and CD8+ T-cell subsets at relapse was evidenced by the loss of response to interferon gamma, tumor necrosis factor α signaling via NF-κB, and interleukin-2/STAT5 signaling. Clonotype analysis of posttransplant relapse samples revealed an expansion of dysfunctional T cells and enrichment of T-regulatory and T-helper cells. Using novel computational methods, our results illustrate a diverse immune-related transcriptional signature in posttransplant relapses not previously reported in pediatric AML.

Determinants of Survival with Combined HER2 and PD-1 Blockade in Metastatic Esophagogastric Cancer.

PURPOSE: We report updated clinical outcomes from a phase II study of pembrolizumab, trastuzumab, and chemotherapy (PTC) in metastatic esophagogastric cancer in conjunction with outcomes from an independent Memorial Sloan Kettering (MSK) cohort. PATIENTS AND METHODS: The significance of pretreatment 89Zr-trastuzumab PET, plasma circulating tumor DNA (ctDNA) dynamics, and tumor HER2 expression and whole exome sequencing was evaluated to identify prognostic biomarkers and mechanisms of resistance in patients treated on-protocol with PTC. Additional prognostic features were evaluated using a multivariable Cox regression model of trastuzumab-treated MSK patients (n = 226). Single-cell RNA sequencing (scRNA-seq) data from MSK and Samsung were evaluated for mechanisms of therapy resistance. RESULTS: 89Zr-trastuzumab PET, scRNA-seq, and serial ctDNA with CT imaging identified how pre-treatment intrapatient genomic heterogeneity contributes to inferior progression-free survival (PFS). We demonstrated that the presence of intensely avid lesions by 89Zr-trastuzumab PET declines in tumor-matched ctDNA by 3 weeks, and clearance of tumor-matched ctDNA by 9 weeks were minimally invasive biomarkers of durable PFS. Paired pre- and on-treatment scRNA-seq identified rapid clearance of HER2-expressing tumor clones with expansion of clones expressing a transcriptional resistance program, which was associated with MT1H, MT1E, MT2A, and MSMB expression. Among trastuzumab-treated patients at MSK, ERBB2 amplification was associated with improved PFS, while alterations in MYC and CDKN2A/B were associated with inferior PFS. CONCLUSIONS: These findings highlight the clinical relevance of identifying baseline intrapatient heterogeneity and serial ctDNA monitoring of HER2-positive esophagogastric cancer patients to identify early evidence of treatment resistance, which could guide proactive therapy escalation or deescalation.

Accurate determination of CRISPR-mediated gene fitness in transplantable tumours.

Assessing tumour gene fitness in physiologically-relevant model systems is challenging due to biological features of in vivo tumour regeneration, including extreme variations in single cell lineage progeny. Here we develop a reproducible, quantitative approach to pooled genetic perturbation in patient-derived xenografts (PDXs), by encoding single cell output from transplanted CRISPR-transduced cells in combination with a Bayesian hierarchical model. We apply this to 181 PDX transplants from 21 breast cancer patients. We show that uncertainty in fitness estimates depends critically on the number of transplant cell clones and the variability in clone sizes. We use a pathway-directed allelic series to characterize Notch signaling, and quantify TP53 / MDM2 drug-gene conditional fitness in outlier patients. We show that fitness outlier identification can be mirrored by pharmacological perturbation. Overall, we demonstrate that the gene fitness landscape in breast PDXs is dominated by inter-patient differences.