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.

Unbiased phenotype and genotype matching maximizes gene discovery and diagnostic yield.

PURPOSE: Widespread application of next-generation sequencing, combined with data exchange platforms, has provided molecular diagnoses for countless families. To maximize diagnostic yield, we implemented an unbiased semi-automated genematching algorithm based on genotype and phenotype matching. METHODS: Rare homozygous variants identified in 2 or more affected individuals, but not in healthy individuals, were extracted from our local database of ∼12,000 exomes. Phenotype similarity scores (PSS), based on human phenotype ontology terms, were assigned to each pair of individuals matched at the genotype level using HPOsim. RESULTS: 33,792 genotype-matched pairs were discovered, representing variants in 7567 unique genes. There was an enrichment of PSS ≥0.1 among pathogenic/likely pathogenic variant-level pairs (94.3% in pathogenic/likely pathogenic variant-level matches vs 34.75% in all matches). We highlighted founder or region-specific variants as an internal positive control and proceeded to identify candidate disease genes. Variant-level matches were particularly helpful in cases involving inframe indels and splice region variants beyond the canonical splice sites, which may otherwise have been disregarded, allowing for detection of candidate disease genes, such as KAT2A, RPAIN, and LAMP3. CONCLUSION: Semi-automated genotype matching combined with PSS is a powerful tool to resolve variants of uncertain significance and to identify candidate disease genes.

Both cell autonomous and non-autonomous processes modulate the association between replication timing and mutation rate.

Cancer somatic mutations are the product of multiple mutational and repair processes, some of which are tightly associated with DNA replication. Mutation rates (MR) are known to be higher in late replication timing (RT) regions, but different processes can affect this association. Systematic analysis of the mutational landscape of 2787 tumors from 32 tumor types revealed that approximately one third of the tumor samples show weak association between replication timing and mutation rate. Further analyses revealed that those samples have unique mutational signatures and are enriched with mutations in genes involved in DNA replication, DNA repair and chromatin structure. Surprisingly, analysis of differentially expressed genes between weak and strong RT-MR association groups revealed that tumors with weak association are enriched with genes associated with cell-cell communication and the immune system, suggesting a non-autonomous response to DNA damage.

Characterization of Unidirectional Replication Forks in the Mouse Genome.

Origins of replication are genomic regions in which replication initiates in a bidirectional manner. Recently, a new methodology (origin-derived single-stranded DNA sequencing; ori-SSDS) was developed that allows the detection of replication initiation in a strand-specific manner. Reanalysis of the strand-specific data revealed that 18-33% of the peaks are non-symmetrical, suggesting a single direction of replication. Analysis of replication fork direction data revealed that these are origins of replication in which the replication is paused in one of the directions, probably due to the existence of a replication fork barrier. Analysis of the unidirectional origins revealed a preference of G4 quadruplexes for the blocked leading strand. Taken together, our analysis identified hundreds of genomic locations in which the replication initiates only in one direction, and suggests that G4 quadruplexes may serve as replication fork barriers in such places.

Mutational signatures association with replication timing in normal cells reveals similarities and differences with matched cancer tissues.

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

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.

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.

RIF1-ASF1-mediated high-order chromatin structure safeguards genome integrity.

The 53BP1-RIF1 pathway antagonizes resection of DNA broken ends and confers PARP inhibitor sensitivity on BRCA1-mutated tumors. However, it is unclear how this pathway suppresses initiation of resection. Here, we identify ASF1 as a partner of RIF1 via an interacting manner similar to its interactions with histone chaperones CAF-1 and HIRA. ASF1 is recruited to distal chromatin flanking DNA breaks by 53BP1-RIF1 and promotes non-homologous end joining (NHEJ) using its histone chaperone activity. Epistasis analysis shows that ASF1 acts in the same NHEJ pathway as RIF1, but via a parallel pathway with the shieldin complex, which suppresses resection after initiation. Moreover, defects in end resection and homologous recombination (HR) in BRCA1-deficient cells are largely suppressed by ASF1 deficiency. Mechanistically, ASF1 compacts adjacent chromatin by heterochromatinization to protect broken DNA ends from BRCA1-mediated resection. Taken together, our findings identify a RIF1-ASF1 histone chaperone complex that promotes changes in high-order chromatin structure to stimulate the NHEJ pathway for DSB repair.

Asynchronous Replication Timing: A Mechanism for Monoallelic Choice During Development.

Developmental programming is carried out by a sequence of molecular choices that epigenetically mark the genome to generate the stable cell types which make up the total organism. A number of important processes, such as genomic imprinting, selection of immune or olfactory receptors, and X-chromosome inactivation in females are dependent on the ability to stably choose one single allele in each cell. In this perspective, we propose that asynchronous replication timing (ASRT) serves as the basis for a sophisticated universal mechanism for mediating and maintaining these decisions.

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.

Chromosomal coordination and differential structure of asynchronous replicating regions.

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Study of mitotic chromatin supports a model of bookmarking by histone modifications and reveals nucleosome deposition patterns.

Mitosis encompasses key molecular changes including chromatin condensation, nuclear envelope breakdown, and reduced transcription levels. Immediately after mitosis, the interphase chromatin structure is reestablished and transcription resumes. The reestablishment of the interphase chromatin is probably achieved by "bookmarking," i.e., the retention of at least partial information during mitosis. To gain a deeper understanding of the contribution of histone modifications to the mitotic bookmarking process, we merged proteomics, immunofluorescence, and ChIP-seq approaches. We focused on key histone modifications and employed HeLa-S3 cells as a model system. Generally, in spite of the general hypoacetylation observed during mitosis, we observed a global concordance between the genomic organization of histone modifications in interphase and mitosis, suggesting that the epigenomic landscape may serve as a component of the mitotic bookmarking process. Next, we investigated the nucleosome that enters nucleosome depleted regions (NDRs) during mitosis. We observed that in ∼60% of the NDRs, the entering nucleosome is distinct from the surrounding highly acetylated nucleosomes and appears to have either low levels of acetylation or high levels of phosphorylation in adjacent residues (since adjacent phosphorylation may interfere with the ability to detect acetylation). Inhibition of histone deacetylases (HDACs) by the small molecule TSA reverts this pattern, suggesting that these nucleosomes are specifically deacetylated during mitosis. Altogether, by merging multiple approaches, our study provides evidence to support a model where histone modifications may play a role in mitotic bookmarking and uncovers new insights into the deposition of nucleosomes during mitosis.

Germline DNA replication timing shapes mammalian genome composition.

Mammalian DNA replication is a highly organized and regulated process. Large, Mb-sized regions are replicated at defined times along S-phase. Replication Timing (RT) is thought to play a role in shaping the mammalian genome by affecting mutation rates. Previous analyses relied on somatic RT profiles. However, only germline mutations are passed on to offspring and affect genomic composition. Therefore, germ cell RT information is necessary to evaluate the influences of RT on the mammalian genome. We adapted the RT mapping technique for limited amounts of cells, and measured RT from two stages in the mouse germline - primordial germ cells (PGCs) and spermatogonial stem cells (SSCs). RT in germline cells exhibited stronger correlations to both mutation rate and recombination hotspots density than those of RT in somatic tissues, emphasizing the importance of using correct tissues-of-origin for RT profiling. Germline RT maps exhibited stronger correlations to additional genetic features including GC-content, transposable elements (SINEs and LINEs), and gene density. GC content stratification and multiple regression analysis revealed independent contributions of RT to SINE, gene, mutation, and recombination hotspot densities. Together, our results establish a central role for RT in shaping multiple levels of mammalian genome composition.

Perturbations in the Replication Program Contribute to Genomic Instability in Cancer.

Cancer and genomic instability are highly impacted by the deoxyribonucleic acid (DNA) replication program. Inaccuracies in DNA replication lead to the increased acquisition of mutations and structural variations. These inaccuracies mainly stem from loss of DNA fidelity due to replication stress or due to aberrations in the temporal organization of the replication process. Here we review the mechanisms and impact of these major sources of error to the replication program.

Genome-wide Determination of Mammalian Replication Timing by DNA Content Measurement.

Replication of the genome occurs during S phase of the cell cycle in a highly regulated process that ensures the fidelity of DNA duplication. Each genomic region is replicated at a distinct time during S phase through the simultaneous activation of multiple origins of replication. Time of replication (ToR) correlates with many genomic and epigenetic features and is linked to mutation rates and cancer. Comprehending the full genomic view of the replication program, in health and disease is a major future goal and challenge. This article describes in detail the "Copy Number Ratio of S/G1 for mapping genomic Time of Replication" method (herein called: CNR-ToR), a simple approach to map the genome wide ToR of mammalian cells. The method is based on the copy number differences between S phase cells and G1 phase cells. The CNR-ToR method is performed in 6 steps: 1. Preparation of cells and staining with propidium iodide (PI); 2. Sorting G1 and S phase cells using fluorescence-activated cell sorting (FACS); 3. DNA purification; 4. Sonication; 5. Library preparation and sequencing; and 6. Bioinformatic analysis. The CNR-ToR method is a fast and easy approach that results in detailed replication maps.

Biases in the SMART-DNA library preparation method associated with genomic poly dA/dT sequences.

Avoiding biases in next generation sequencing (NGS) library preparation is crucial for obtaining reliable sequencing data. Recently, a new library preparation method has been introduced which has eliminated the need for the ligation step. This method, termed SMART (switching mechanism at the 5' end of the RNA transcript), is based on template switching reverse transcription. To date, there has been no systematic analysis of the additional biases introduced by this method. We analysed the genomic distribution of sequenced reads prepared from genomic DNA using the SMART methodology and found a strong bias toward long (≥12bp) poly dA/dT containing genomic loci. This bias is unique to the SMART-based library preparation and does not appear when libraries are prepared with conventional ligation based methods. Although this bias is obvious only when performing paired end sequencing, it affects single end sequenced samples as well. Our analysis demonstrates that sequenced reads originating from SMART-DNA libraries are heavily skewed toward genomic poly dA/dT tracts. This bias needs to be considered when deciding to use SMART based technology for library preparation.

Persistence to anti-cancer treatments in the stationary to proliferating transition.

The heterogeneous responses of clonal cancer cells to treatment is understood to be caused by several factors, including stochasticity, cell-cycle dynamics, and different micro-environments. In a tumor, cancer cells may encounter fluctuating conditions and transit from a stationary culture to a proliferating state, for example this may occur following treatment. Here, we undertake a quantitative evaluation of the response of single cancerous lymphoblasts (L1210 cells) to various treatments administered during this transition. Additionally, we developed an experimental system, a "Mammalian Mother Machine," that tracks the fate of thousands of mammalian cells over several generations under transient exposure to chemotherapeutic drugs. Using our developed system, we were able to follow the same cell under repeated treatments and continuously track many generations. We found that the dynamics of the transition between stationary and proliferative states are highly variable and affect the response to drug treatment. Using cell-cycle markers, we were able to isolate a subpopulation of persister cells with distinctly higher than average survival probability. The higher survival rate encountered with cell-cycle phase specific drugs was associated with a significantly longer time-till-division, and was reduced by a non cell-cycle specific drug. Our results suggest that the variability of transition times from the stationary to the proliferating state may be an obstacle hampering the effectiveness of drugs and should be taken into account when designing treatment regimens.

The mutation spectrum in genomic late replication domains shapes mammalian GC content.

Genome sequence compositions and epigenetic organizations are correlated extensively across multiple length scales. Replication dynamics, in particular, is highly correlated with GC content. We combine genome-wide time of replication (ToR) data, topological domains maps and detailed functional epigenetic annotations to study the correlations between replication timing and GC content at multiple scales. We find that the decrease in genomic GC content at large scale late replicating regions can be explained by mutation bias favoring A/T nucleotide, without selection or biased gene conversion. Quantification of the free dNTP pool during the cell cycle is consistent with a mechanism involving replication-coupled mutation spectrum that favors AT nucleotides at late S-phase. We suggest that mammalian GC content composition is shaped by independent forces, globally modulating mutation bias and locally selecting on functional element. Deconvoluting these forces and analyzing them on their native scales is important for proper characterization of complex genomic correlations.

Distinguishing between stochasticity and determinism: Examples from cell cycle duration variability.

We describe a recent approach for distinguishing between stochastic and deterministic sources of variability, focusing on the mammalian cell cycle. Variability between cells is often attributed to stochastic noise, although it may be generated by deterministic components. Interestingly, lineage information can be used to distinguish between variability and determinism. Analysis of correlations within a lineage of the mammalian cell cycle duration revealed its deterministic nature. Here, we discuss the sources of such variability and the possibility that the underlying deterministic process is due to the circadian clock. Finally, we discuss the "kicked cell cycle" model and its implication on the study of the cell cycle in healthy and cancerous tissues.

Lineage correlations of single cell division time as a probe of cell-cycle dynamics.

Stochastic processes in cells are associated with fluctuations in mRNA, protein production and degradation, noisy partition of cellular components at division, and other cell processes. Variability within a clonal population of cells originates from such stochastic processes, which may be amplified or reduced by deterministic factors. Cell-to-cell variability, such as that seen in the heterogeneous response of bacteria to antibiotics, or of cancer cells to treatment, is understood as the inevitable consequence of stochasticity. Variability in cell-cycle duration was observed long ago; however, its sources are still unknown. A central question is whether the variance of the observed distribution originates from stochastic processes, or whether it arises mostly from a deterministic process that only appears to be random. A surprising feature of cell-cycle-duration inheritance is that it seems to be lost within one generation but to be still present in the next generation, generating poor correlation between mother and daughter cells but high correlation between cousin cells. This observation suggests the existence of underlying deterministic factors that determine the main part of cell-to-cell variability. We developed an experimental system that precisely measures the cell-cycle duration of thousands of mammalian cells along several generations and a mathematical framework that allows discrimination between stochastic and deterministic processes in lineages of cells. We show that the inter- and intra-generation correlations reveal complex inheritance of the cell-cycle duration. Finally, we build a deterministic nonlinear toy model for cell-cycle inheritance that reproduces the main features of our data. Our approach constitutes a general method to identify deterministic variability in lineages of cells or organisms, which may help to predict and, eventually, reduce cell-to-cell heterogeneity in various systems, such as cancer cells under treatment.