A Cyclic Permutation Approach to Removing Spatial Dependency between Clustered Gene Ontology Terms.

Traditional gene set enrichment analysis falters when applied to large genomic domains, where neighboring genes often share functions. This spatial dependency creates misleading enrichments, mistaking mere physical proximity for genuine biological connections. Here we present Spatial Adjusted Gene Ontology (SAGO), a novel cyclic permutation-based approach, to tackle this challenge. SAGO separates enrichments due to spatial proximity from genuine biological links by incorporating the genes' spatial arrangement into the analysis. We applied SAGO to various datasets in which the identified genomic intervals are large, including replication timing domains, large H3K9me3 and H3K27me3 domains, HiC compartments and lamina-associated domains (LADs). Intriguingly, applying SAGO to prostate cancer samples with large copy number alteration (CNA) domains eliminated most of the enriched GO terms, thus helping to accurately identify biologically relevant gene sets linked to oncogenic processes, free from spatial bias.

S Phase Duration Is Determined by Local Rate and Global Organization of Replication.

The duration of the cell cycle has been extensively studied and a wide degree of variability exists between cells, tissues and organisms. However, the duration of S phase has often been neglected, due to the false assumption that S phase duration is relatively constant. In this paper, we describe the methodologies to measure S phase duration, summarize the existing knowledge about its variability and discuss the key factors that control it. The local rate of replication (LRR), which is a combination of fork rate (FR) and inter-origin distance (IOD), has a limited influence on S phase duration, partially due to the compensation between FR and IOD. On the other hand, the organization of the replication program, specifically the amount of replication domains that fire simultaneously and the degree of overlap between the firing of distinct replication timing domains, is the main determinant of S phase duration. We use these principles to explain the variation in S phase length in different tissues and conditions.

Mitotic H3K9ac is controlled by phase-specific activity of HDAC2, HDAC3, and SIRT1.

Histone acetylation levels are reduced during mitosis. To study the mitotic regulation of H3K9ac, we used an array of inhibitors targeting specific histone deacetylases. We evaluated the involvement of the targeted enzymes in regulating H3K9ac during all mitotic stages by immunofluorescence and immunoblots. We identified HDAC2, HDAC3, and SIRT1 as modulators of H3K9ac mitotic levels. HDAC2 inhibition increased H3K9ac levels in prophase, whereas HDAC3 or SIRT1 inhibition increased H3K9ac levels in metaphase. Next, we performed ChIP-seq on mitotic-arrested cells following targeted inhibition of these histone deacetylases. We found that both HDAC2 and HDAC3 have a similar impact on H3K9ac, and inhibiting either of these two HDACs substantially increases the levels of this histone acetylation in promoters, enhancers, and insulators. Altogether, our results support a model in which H3K9 deacetylation is a stepwise process-at prophase, HDAC2 modulates most transcription-associated H3K9ac-marked loci, and at metaphase, HDAC3 maintains the reduced acetylation, whereas SIRT1 potentially regulates H3K9ac by impacting HAT activity.

Cancer Mutational Processes Vary in Their Association with Replication Timing and Chromatin Accessibility.

Cancer somatic mutations are the product of multiple mutational and repair processes, both of which are tightly associated with DNA replication. Distinctive patterns of somatic mutation accumulation, termed mutational signatures, are indicative of processes sustained within tumors. However, the association of various mutational processes with replication timing (RT) remains an open question. In this study, we systematically analyzed the mutational landscape of 2,787 tumors from 32 tumor types separately for early and late replicating regions using sequence context normalization and chromatin data to account for sequence and chromatin accessibility differences. To account for sequence differences between various genomic regions, an artificial genome-based approach was developed to expand the signature analyses to doublet base substitutions and small insertions and deletions. The association of mutational processes and RT was signature specific: Some signatures were associated with early or late replication (such as SBS7b and SBS7a, respectively), and others had no association. Most associations existed even after normalizing for genome accessibility. A focused mutational signature identification approach was also developed that uses RT information to improve signature identification; this approach found that SBS16, which is biased toward early replication, is strongly associated with better survival rates in liver cancer. Overall, this novel and comprehensive approach provides a better understanding of the etiology of mutational signatures, which may lead to improved cancer prevention, diagnosis, and treatment. SIGNIFICANCE: Many mutational processes associate with early or late replication timing regions independently of chromatin accessibility, enabling development of a focused identification approach to improve mutational signature detection.

Cardiac leptin overexpression in the context of acute MI and reperfusion potentiates myocardial remodeling and left ventricular dysfunction.

BACKGROUND: Acute MI induces leptin expression in the heart, however the role of myocardial leptin in cardiac ischemia and reperfusion (IR) remains unknown. To shed light on the effects of elevated levels of leptin in the myocardium, we overexpressed cardiac leptin and assessed local remodeling and myocardial function in this context. METHODS AND RESULTS: Cardiac leptin overexpression was stimulated in mice undergoing IR by a single intraperitoneal injection of leptin antagonist (LepA). All mice exhibited a normal pattern of body weight gain. A rapid, long-term upregulation of leptin mRNA was demonstrated in the heart, adipose, and liver tissues in IR/LepA-treated mice. Overexpressed cardiac leptin mRNA extended beyond postoperative day (POD) 30. Plasma leptin peaked 7.5 hours postoperatively, especially in IR/LepA-treated mice, subsiding to normal levels by 24 hours. On POD-30 IR/LepA-treated mice demonstrated cardiomyocyte hypertrophy and perivascular fibrosis compared to IR/saline controls. Echocardiography on POD-30 demonstrated eccentric hypertrophy and systolic dysfunction in IR/LepA. We recorded reductions in Ejection Fraction (p<0.001), Fraction Shortening (p<0.01), and Endocardial Fraction Area Change (p<0.01), and an increase in Endocardial Area Change (p<0.01). Myocardial remodeling in the context of IR and cardiac leptin overexpression was associated with increased cardiac TGFß ligand expression, activated Smad2, and downregulation of STAT3 activity. CONCLUSIONS: Cardiac IR coinciding with increased myocardial leptin synthesis promotes cardiomyocyte hypertrophy and fibrosis and potentiates myocardial dysfunction. Plasma leptin levels do not reflect cardiac leptin synthesis, and may not predict leptin-related cardiovascular morbidity. Targeting cardiac leptin is a potential treatment for cardiac IR damage.