Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

FOXP1 and KLF2 reciprocally regulate checkpoints of stem-like to effector transition in CAR T cells.

In cancer and infections, self-renewing stem-like CD8+ T cells mediate the response of immunotherapies and replenish terminally exhausted T cells and effector-like T cells. However, the programs governing the lineage choice in chimeric antigen receptor (CAR) T cells are unclear. Here, by simultaneously profiling single-cell chromatin accessibility and transcriptome in the same CAR T cells, we identified heterogeneous chromatin states within CD8+ T cell subsets that foreshadowed transcriptional changes and were primed for regulation by distinct transcription factors. Transcription factors that controlled each CD8+ T cell subset were regulated by high numbers of enhancers and positioned as hubs of gene networks. FOXP1, a hub in the stem-like network, promoted expansion and stemness of CAR T cells and limited excessive effector differentiation. In the effector network, KLF2 enhanced effector CD8+ T cell differentiation and prevented terminal exhaustion. Thus, we identified gene networks and hub transcription factors that controlled the differentiation of stem-like CD8+ CAR T cells into effector or exhausted CD8+ CAR T cells.

Kinome-wide siRNA screen identifies a DCLK2-TBK1 oncogenic signaling axis in clear cell renal cell carcinoma.

TANK-binding kinase 1 (TBK1) is a potential therapeutic target in multiple cancers, including clear cell renal cell carcinoma (ccRCC). However, targeting TBK1 in clinical practice is challenging. One approach to overcome this challenge would be to identify an upstream TBK1 regulator that could be targeted therapeutically in cancer specifically. In this study, we perform a kinome-wide small interfering RNA (siRNA) screen and identify doublecortin-like kinase 2 (DCLK2) as a TBK1 regulator in ccRCC. DCLK2 binds to and directly phosphorylates TBK1 on Ser172. Depletion of DCLK2 inhibits anchorage-independent colony growth and kidney tumorigenesis in orthotopic xenograft models. Conversely, overexpression of DCLK2203, a short isoform that predominates in ccRCC, promotes ccRCC cell growth and tumorigenesis in vivo. Mechanistically, DCLK2203 elicits its oncogenic signaling via TBK1 phosphorylation and activation. Taken together, these results suggest that DCLK2 is a TBK1 activator and potential therapeutic target for ccRCC.

Paracrine Induction of HIF by Glutamate in Breast Cancer: EglN1 Senses Cysteine.

The HIF transcription factor promotes adaptation to hypoxia and stimulates the growth of certain cancers, including triple-negative breast cancer (TNBC). The HIFα subunit is usually prolyl-hydroxylated by EglN family members under normoxic conditions, causing its rapid degradation. We confirmed that TNBC cells secrete glutamate, which we found is both necessary and sufficient for the paracrine induction of HIF1α in such cells under normoxic conditions. Glutamate inhibits the xCT glutamate-cystine antiporter, leading to intracellular cysteine depletion. EglN1, the main HIFα prolyl-hydroxylase, undergoes oxidative self-inactivation in the absence of cysteine both in biochemical assays and in cells, resulting in HIF1α accumulation. Therefore, EglN1 senses both oxygen and cysteine.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer.

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

Blockade of the checkpoint receptor TIGIT prevents NK cell exhaustion and elicits potent anti-tumor immunity.

Checkpoint blockade enhances effector T cell function and has elicited long-term remission in a subset of patients with a broad spectrum of cancers. TIGIT is a checkpoint receptor thought to be involved in mediating T cell exhaustion in tumors; however, the relevance of TIGIT to the dysfunction of natural killer (NK) cells remains poorly understood. Here we found that TIGIT, but not the other checkpoint molecules CTLA-4 and PD-1, was associated with NK cell exhaustion in tumor-bearing mice and patients with colon cancer. Blockade of TIGIT prevented NK cell exhaustion and promoted NK cell-dependent tumor immunity in several tumor-bearing mouse models. Furthermore, blockade of TIGIT resulted in potent tumor-specific T cell immunity in an NK cell-dependent manner, enhanced therapy with antibody to the PD-1 ligand PD-L1 and sustained memory immunity in tumor re-challenge models. This work demonstrates that TIGIT constitutes a previously unappreciated checkpoint in NK cells and that targeting TIGIT alone or in combination with other checkpoint receptors is a promising anti-cancer therapeutic strategy.

Regioselective hydroformylation of propene catalysed by rhodium-zeolite.

Hydroformylation is an industrial process for the production of aldehydes from alkenes1,2. Regioselective hydroformylation of propene to high-value n-butanal is particularly important, owing to a wide range of bulk applications of n-butanal in the manufacture of various necessities in human daily life3. Supported rhodium (Rh) hydroformylation catalysts, which often excel in catalyst recyclability, ease of separation and adaptability for continuous-flow processes, have been greatly exploited4. Nonetheless, they usually consist of rotationally flexible and sterically unconstrained Rh hydride dicarbonyl centres, only affording limited regioselectivity to n-butanal5-8. Here we show that proper encapsulation of Rh species comprising Rh(I)-gem-dicarbonyl centres within a MEL zeolite framework allows the breaking of the above model. The optimized catalyst exhibits more than 99% regioselectivity to n-butanal and more than 99% selectivity to aldehydes at a product formation turnover frequency (TOF) of 6,500 h-1, surpassing the performance of all heterogeneous and most homogeneous catalysts developed so far. Our comprehensive studies show that the zeolite framework can act as a scaffold to steer the reaction pathway of the intermediates confined in the space between the zeolite framework and Rh centres towards the exclusive formation of n-butanal.

An oncogenic JMJD6-DGAT1 axis tunes the epigenetic regulation of lipid droplet formation in clear cell renal cell carcinoma.

Characterized by intracellular lipid droplet accumulation, clear cell renal cell carcinoma (ccRCC) is resistant to cytotoxic chemotherapy and is a lethal disease. Through an unbiased siRNA screen of 2-oxoglutarate (2-OG)-dependent enzymes, which play a critical role in tumorigenesis, we identified Jumonji domain-containing 6 (JMJD6) as an essential gene for ccRCC tumor development. The downregulation of JMJD6 abolished ccRCC colony formation in vitro and inhibited orthotopic tumor growth in vivo. Integrated ChIP-seq and RNA-seq analyses uncovered diacylglycerol O-acyltransferase 1 (DGAT1) as a critical JMJD6 effector. Mechanistically, JMJD6 interacted with RBM39 and co-occupied DGAT1 gene promoter with H3K4me3 to induce DGAT1 expression. JMJD6 silencing reduced DGAT1, leading to decreased lipid droplet formation and tumorigenesis. The pharmacological inhibition (or depletion) of DGAT1 inhibited lipid droplet formation in vitro and ccRCC tumorigenesis in vivo. Thus, the JMJD6-DGAT1 axis represents a potential new therapeutic target for ccRCC.

Genome-wide Screening Identifies SFMBT1 as an Oncogenic Driver in Cancer with VHL Loss.

von Hippel-Lindau (VHL) is a critical tumor suppressor in clear cell renal cell carcinomas (ccRCCs). It is important to identify additional therapeutic targets in ccRCC downstream of VHL loss besides hypoxia-inducible factor 2α (HIF2α). By performing a genome-wide screen, we identified Scm-like with four malignant brain tumor domains 1 (SFMBT1) as a candidate pVHL target. SFMBT1 was considered to be a transcriptional repressor but its role in cancer remains unclear. ccRCC patients with VHL loss-of-function mutations displayed elevated SFMBT1 protein levels. SFMBT1 hydroxylation on Proline residue 651 by EglN1 mediated its ubiquitination and degradation governed by pVHL. Depletion of SFMBT1 abolished ccRCC cell proliferation in vitro and inhibited orthotopic tumor growth in vivo. Integrated analyses of ChIP-seq, RNA-seq, and patient prognosis identified sphingosine kinase 1 (SPHK1) as a key SFMBT1 target gene contributing to its oncogenic phenotype. Therefore, the pVHL-SFMBT1-SPHK1 axis serves as a potential therapeutic avenue for ccRCC.

Interface nano-optics with van der Waals polaritons.

Polaritons are hybrid excitations of matter and photons. In recent years, polaritons in van der Waals nanomaterials-known as van der Waals polaritons-have shown great promise to guide the flow of light at the nanoscale over spectral regions ranging from the visible to the terahertz. A vibrant research field based on manipulating strong light-matter interactions in the form of polaritons, supported by these atomically thin van der Waals nanomaterials, is emerging for advanced nanophotonic and opto-electronic applications. Here we provide an overview of the state of the art of exploiting interface optics-such as refractive optics, meta-optics and moiré engineering-for the control of van der Waals polaritons. This enhanced control over van der Waals polaritons at the nanoscale has not only unveiled many new phenomena, but has also inspired valuable applications-including new avenues for nano-imaging, sensing, on-chip optical circuitry, and potentially many others in the years to come.

Topological polaritons and photonic magic angles in twisted α-MoO3 bilayers.

Twisted two-dimensional bilayer materials exhibit many exotic electronic phenomena. Manipulating the 'twist angle' between the two layers enables fine control of the electronic band structure, resulting in magic-angle flat-band superconductivity1,2, the formation of moiré excitons3-8 and interlayer magnetism9. However, there are limited demonstrations of such concepts for photons. Here we show how analogous principles, combined with extreme anisotropy, enable control and manipulation of the photonic dispersion of phonon polaritons in van der Waals bilayers. We experimentally observe tunable topological transitions from open (hyperbolic) to closed (elliptical) dispersion contours in bilayers of α-phase molybdenum trioxide (α-MoO3), arising when the rotation between the layers is at a photonic magic twist angle. These transitions are induced by polariton hybridization and are controlled by a topological quantity. At the transitions the bilayer dispersion flattens, exhibiting low-loss tunable polariton canalization and diffractionless propagation with a resolution of less than λ0/40, where λ0 is the free-space wavelength. Our findings extend twistronics10 and moiré physics to nanophotonics and polaritonics, with potential applications in nanoimaging, nanoscale light propagation, energy transfer and quantum physics.

METTL3 is essential for normal progesterone signaling during embryo implantation via m6A-mediated translation control of progesterone receptor.

Embryo implantation, a crucial step in human reproduction, is tightly controlled by estrogen and progesterone (P4) via estrogen receptor alpha and progesterone receptor (PGR), respectively. Here, we report that N6-methyladenosine (m6A), the most abundant mRNA modification in eukaryotes, plays an essential role in embryo implantation through the maintenance of P4 signaling. Conditional deletion of methyltransferase-like 3 (Mettl3), encoding the m6A writer METTL3, in the female reproductive tract using a Cre mouse line with Pgr promoter (Pgr-Cre) resulted in complete implantation failure due to pre-implantation embryo loss and defective uterine receptivity. Moreover, the uterus of Mettl3 null mice failed to respond to artificial decidualization. We further found that Mettl3 deletion was accompanied by a marked decrease in PGR protein expression. Mechanistically, we found that Pgr mRNA is a direct target for METTL3-mediated m6A modification. A luciferase assay revealed that the m6A modification in the 5' untranslated region (5'-UTR) of Pgr mRNA enhances PGR protein translation efficiency in a YTHDF1-dependent manner. Finally, we demonstrated that METTL3 is required for human endometrial stromal cell decidualization in vitro and that the METTL3-PGR axis is conserved between mice and humans. In summary, this study provides evidence that METTL3 is essential for normal P4 signaling during embryo implantation via m6A-mediated translation control of Pgr mRNA.

Chromatin accessibility mediated transcriptome changes contribute to flavor substance alterations and jasmonic acid hyperaccumulation during oolong tea withering process.

During the oolong tea withering process, abiotic stresses induce significant changes in the content of various flavor substances and jasmonic acid (JA). However, the changes in chromatin accessibility during withering and their potential impact remain poorly understood. By integrating ATAC-seq, RNA-seq, metabolite, and hormone assays, we characterized the withering treatment-induced changes in chromatin accessibility, gene expression levels, important metabolite contents, and JA and JA-ILE contents. Additionally, we analyzed the effects of chromatin accessibility alterations on gene expression changes, content changes of important flavor substances, and JA hyperaccumulation. Our analysis identified a total of 3451 open- and 13 426 close-differentially accessible chromatin regions (DACRs) under withering treatment. Our findings indicate that close-DACRs-mediated down-regulated differentially expressed genes (DEGs) resulted in the reduced accumulation of multiple catechins during withering, whereas open-DACRs-mediated up-regulated DEGs contributed to the increased accumulation of important terpenoids, JA, JA-ILE and short-chain C5/C6 volatiles. We further highlighted important DACRs-mediated DEGs associated with the synthesis of catechins, terpenoids, JA and JA and short-chain C5/C6 volatiles and confirmed the broad effect of close-DACRs on catechin synthesis involving almost all enzymes in the pathway during withering. Importantly, we identified a novel MYB transcription factor (CsMYB83) regulating catechin synthesis and verified the binding of CsMYB83 in the promoter-DACRs regions of key catechin synthesis genes using DAP-seq. Overall, our results not only revealed a landscape of chromatin alters-mediated transcription, flavor substance and hormone changes under oolong tea withering, but also provided target genes for flavor improvement breeding in tea plant.

A social theory-enhanced graph representation learning framework for multitask prediction of drug-drug interactions.

Current machine learning-based methods have achieved inspiring predictions in the scenarios of mono-type and multi-type drug-drug interactions (DDIs), but they all ignore enhancive and depressive pharmacological changes triggered by DDIs. In addition, these pharmacological changes are asymmetric since the roles of two drugs in an interaction are different. More importantly, these pharmacological changes imply significant topological patterns among DDIs. To address the above issues, we first leverage Balance theory and Status theory in social networks to reveal the topological patterns among directed pharmacological DDIs, which are modeled as a signed and directed network. Then, we design a novel graph representation learning model named SGRL-DDI (social theory-enhanced graph representation learning for DDI) to realize the multitask prediction of DDIs. SGRL-DDI model can capture the task-joint information by integrating relation graph convolutional networks with Balance and Status patterns. Moreover, we utilize task-specific deep neural networks to perform two tasks, including the prediction of enhancive/depressive DDIs and the prediction of directed DDIs. Based on DDI entries collected from DrugBank, the superiority of our model is demonstrated by the comparison with other state-of-the-art methods. Furthermore, the ablation study verifies that Balance and Status patterns help characterize directed pharmacological DDIs, and that the joint of two tasks provides better DDI representations than individual tasks. Last, we demonstrate the practical effectiveness of our model by a version-dependent test, where 88.47 and 81.38% DDI out of newly added entries provided by the latest release of DrugBank are validated in two predicting tasks respectively.

Membrane-dependent actin polymerization mediated by the Legionella pneumophila effector protein MavH.

L. pneumophila propagates in eukaryotic cells within a specialized niche, the Legionella-containing vacuole (LCV). The infection process is controlled by over 330 effector proteins delivered through the type IV secretion system. In this study, we report that the Legionella MavH effector localizes to endosomes and remodels host actin cytoskeleton in a phosphatidylinositol 3-phosphate (PI(3)P) dependent manner when ectopically expressed. We show that MavH recruits host actin capping protein (CP) and actin to the endosome via its CP-interacting (CPI) motif and WH2-like actin-binding domain, respectively. In vitro assays revealed that MavH stimulates actin assembly on PI(3)P-containing liposomes causing their tubulation. In addition, the recruitment of CP by MavH negatively regulates F-actin density at the membrane. We further show that, in L. pneumophila-infected cells, MavH appears around the LCV at the very early stage of infection and facilitates bacterium entry into the host. Together, our results reveal a novel mechanism of membrane tubulation induced by membrane-dependent actin polymerization catalyzed by MavH that contributes to the early stage of L. pneumophila infection by regulating host actin dynamics.

SOX2 downregulation of PML increases HCMV gene expression and growth of glioma cells.

The presence of human cytomegalovirus (HCMV) in glioblastoma (GBM) and improved outcomes of GBM patients receiving therapies targeting the virus have implicated HCMV in GBM progression. However, a unifying mechanism that accounts for the contribution of HCMV to the malignant phenotype of GBM remains incompletely defined. Here we have identified SOX2, a marker of glioma stem cells (GSCs), as a key determinant of HCMV gene expression in gliomas. Our studies demonstrated that SOX2 downregulated promyelocytic leukemia (PML) and Sp100 and consequently facilitated viral gene expression by decreasing the amount of PML nuclear bodies in HCMV-infected glioma cells. Conversely, the expression of PML antagonized the effects of SOX2 on HCMV gene expression. Furthermore, this regulation of SOX2 on HCMV infection was demonstrated in a neurosphere assay of GSCs and in a murine xenograft model utilizing xenografts from patient-derived glioma tissue. In both cases, SOX2 overexpression facilitated the growth of neurospheres and xenografts implanted in immunodeficient mice. Lastly, the expression of SOX2 and HCMV immediate early 1 (IE1) protein could be correlated in tissues from glioma patients, and interestingly, elevated levels of SOX2 and IE1 were predictive of a worse clinical outcome. These studies argue that HCMV gene expression in gliomas is regulated by SOX2 through its regulation of PML expression and that targeting molecules in this SOX2-PML pathway could identify therapies for glioma treatment.

Pan-cancer analysis of CDKN2A alterations identifies a subset of gastric cancer with a cold tumor immune microenvironment.

BACKGROUND: Although CDKN2A alteration has been explored as a favorable factor for tumorigenesis in pan-cancers, the association between CDKN2A point mutation (MUT) and intragenic deletion (DEL) and response to immune checkpoint inhibitors (ICIs) is still disputed. This study aims to determine the associations of CDKN2A MUT and DEL with overall survival (OS) and response to immune checkpoint inhibitors treatment (ICIs) among pan-cancers and the clinical features of CDKN2A-altered gastric cancer. METHODS: This study included 45,000 tumor patients that underwent tumor sequencing across 33 cancer types from four cohorts, the MSK-MetTropism, MSK-IMPACT, OrigiMed2020 and TCGA cohorts. Clinical outcomes and genomic factors associated with response to ICIs, including tumor mutational burden, copy number alteration, neoantigen load, microsatellite instability, tumor immune microenvironment and immune-related gene signatures, were collected in pan-cancer. Clinicopathologic features and outcomes were assessed in gastric cancer. Patients were grouped based on the presence of CDKN2A wild type (WT), CDKN2A MUT, CDKN2A DEL and CDKN2A other alteration (ALT). RESULTS: Our research showed that CDKN2A-MUT patients had shorter survival times than CDKN2A-WT patients in the MSK MetTropism and TCGA cohorts, but longer OS in the MSK-IMPACT cohort with ICIs treatment, particularly in patients having metastatic disease. Similar results were observed among pan-cancer patients with CDKN2A DEL and other ALT. Notably, CDKN2A ALT frequency was positively related to tumor-specific objective response rates to ICIs in MSK MetTropism and OrigiMed 2020. Additionally, individuals with esophageal carcinoma or stomach adenocarcinoma who had CDKN2A MUT had poorer OS than patients from the MSK-IMPACT group, but not those with adenocarcinoma. We also found reduced levels of activated NK cells, T cells CD8 and M2 macrophages in tumor tissue from CDKN2A-MUT or DEL pan-cancer patients compared to CDKN2A-WT patients in TCGA cohort. Gastric cancer scRNA-seq data also showed that CDKN2A-ALT cancer contained less CD8 T cells but more exhausted T cells than CDKN2A-WT cancer. A crucial finding of the pathway analysis was the inhibition of three immune-related pathways in the CDKN2A ALT gastric cancer patients, including the interferon alpha response, inflammatory response, and interferon gamma response. CONCLUSIONS: This study illustrates the CDKN2A MUT and DEL were associated with a poor outcome across cancers. CDKN2A ALT, on the other hand, have the potential to be used as a biomarker for choosing patients for ICI treatment, notably in esophageal carcinoma and stomach adenocarcinoma.

Impact of Bmal1 Rescue and Time-Restricted Feeding on Liver and Muscle Proteomes During the Active Phase in Mice.

Molecular clocks and daily feeding cycles support metabolism in peripheral tissues. Although the roles of local clocks and feeding are well defined at the transcriptional level, their impact on governing protein abundance in peripheral tissues is unclear. Here, we determine the relative contributions of local molecular clocks and daily feeding cycles on liver and muscle proteomes during the active phase in mice. LC-MS/MS was performed on liver and gastrocnemius muscle harvested 4 h into the dark phase from WT, Bmal1 KO, and dual liver- and muscle-Bmal1-rescued mice under either ad libitum feeding or time-restricted feeding during the dark phase. Feeding-fasting cycles had only minimal effects on levels of liver proteins and few, if any, on the muscle proteome. In contrast, Bmal1 KO altered the abundance of 674 proteins in liver and 80 proteins in muscle. Local rescue of liver and muscle Bmal1 restored ∼50% of proteins in liver and ∼25% in muscle. These included proteins involved in fatty acid oxidation in liver and carbohydrate metabolism in muscle. For liver, proteins involved in de novo lipogenesis were largely dependent on Bmal1 function in other tissues (i.e., the wider clock system). Proteins regulated by BMAL1 in liver and muscle were enriched for secreted proteins. We found that the abundance of fibroblast growth factor 1, a liver secreted protein, requires BMAL1 and that autocrine fibroblast growth factor 1 signaling modulates mitochondrial respiration in hepatocytes. In liver and muscle, BMAL1 is a more potent regulator of dark phase proteomes than daily feeding cycles, highlighting the need to assess protein levels in addition to mRNA when investigating clock mechanisms. The proteome is more extensively regulated by BMAL1 in liver than in muscle, and many metabolic pathways in peripheral tissues are reliant on the function of the clock system as a whole.

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice.

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.