3D genome organization contributes to genome instability at fragile sites

A three-dimensional genomics view for speciation research

Genomic information is folded in a three-dimensional (3D) structure, a rarely explored evolutionary driver of speciation. Technological advances now enable the study of 3D genome structures (3DGSs) across the Tree of Life. At the onset of 3D speciation genomics, we discuss the putative roles of 3DGSs in speciation.

Linear interaction between replication and transcription shapes DNA break dynamics at recurrent DNA break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unravel the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA reveals co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity uncovers that head-to-head interactions between replication and transcription machinery result in 60% DSB contribution to the head-on compared to 40% for co-directional. Taken together we reveal TTR as a fragile class and show how the linear interaction between transcription and replication impacts genome stability.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

A Unifying Hypothesis for the Genome Dynamics Proposed to Underlie Neuropsychiatric Phenotypes

The sheer number of gene variants and the extent of the observed clinical and molecular heterogeneity recorded in neuropsychiatric disorders (NPDs) could be due to the magnified downstream effects initiated by a smaller group of genomic higher-order alterations in response to endogenous or environmental stress. Chromosomal common fragile sites (CFS) are functionally linked with microRNAs, gene copy number variants (CNVs), sub-microscopic deletions and duplications of DNA, rare single-nucleotide variants (SNVs/SNPs), and small insertions/deletions (indels), as well as chromosomal translocations, gene duplications, altered methylation, microRNA and L1 transposon activity, and 3-D chromosomal topology characteristics. These genomic structural features have been linked with various NPDs in mostly isolated reports and have usually only been viewed as areas harboring potential candidate genes of interest. The suggestion to use a higher level entry point (the 'fragilome' and associated features) activated by a central mechanism ('stress') for studying NPD genetics has the potential to unify the existing vast number of different observations in this field. This approach may explain the continuum of gene findings distributed between affected and unaffected individuals, the clustering of NPD phenotypes and overlapping comorbidities, the extensive clinical and molecular heterogeneity, and the association with certain other medical disorders.

Time is ticking faster for long genes in aging

Recent studies of aging organisms have identified a systematic phenomenon, characterized by a negative correlation between gene length and their expression in various cell types, species, and diseases. We term this phenomenon gene-length-dependent transcription decline (GLTD) and suggest that it may represent a bottleneck in the transcription machinery and thereby significantly contribute to aging as an etiological factor. We review potential links between GLTD and key aging processes such as DNA damage and explore their potential in identifying disease modification targets. Notably, in Alzheimer's disease, GLTD spotlights extremely long synaptic genes at chromosomal fragile sites (CFSs) and their vulnerability to postmitotic DNA damage. We suggest that GLTD is an integral element of biological aging.

Cell cycle gene alterations associate with a redistribution of mutation risk across chromosomal domains in human cancers

Mutations in human cells exhibit increased burden in heterochromatic, late DNA replication time (RT) chromosomal domains, with variation in mutation rates between tissues mirroring variation in heterochromatin and RT. We observed that regional mutation risk further varies between individual tumors in a manner independent of cell type, identifying three signatures of domain-scale mutagenesis in >4,000 tumor genomes. The major signature reflects remodeling of heterochromatin and of the RT program domains seen across tumors, tissues and cultured cells, and is robustly linked with higher expression of cell proliferation genes. Regional mutagenesis is associated with loss of activity of the tumor-suppressor genes RB1 and TP53, consistent with their roles in cell cycle control, with distinct mutational patterns generated by the two genes. Loss of regional heterogeneity in mutagenesis is associated with deficiencies in various DNA repair pathways. These mutation risk redistribution processes modify the mutation supply towards important genes, diverting the course of somatic evolution.

Insights into common fragile site instability: DNA replication challenges at DNA repeat sequences

Common fragile sites (CFS) are specific genomic regions prone to chromosomal instability under conditions of DNA replication stress. CFSs manifest as breaks, gaps, and constrictions on metaphase chromosomes under mild replication stress. These replication-sensitive CFS regions are preferentially unstable during cancer development, as reflected by their association with copy number variants (CNVs) frequently arise in most tumor types. Over the years, it became clear that a combination of different characteristics underlies the enhanced sensitivity of CFSs to replication stress. As of today, there is a strong evidence that the core fragility regions along CFSs overlap with actively transcribed large genes with delayed replication timing upon replication stress. Recently, the mechanistic basis for CFS instability was further extended to regions which span topologically associated domain (TAD) boundaries, generating a fragility signature composed of replication, transcription and genome organization. The presence of difficult-to-replicate AT-rich repeats was one of the early features suggested to characterize a subgroup of CFSs. These long stretches of AT-dinucleotide have the potential to fold into stable secondary structures which may impede replication fork progression, leaving the region under-replicated. Here, we focus on the molecular mechanisms underlying repeat instability at CFSs and on the proteins involved in the resolution of secondary structure impediments arising along repetitive sequence elements which are essential for the maintenance of genome stability.

Identification of Genomic Instability in Cows Infected with BVD Virus

An important factor for dairy cattle farmers is the profitability of cattle rearing, which is influenced by the animals' health and reproductive parameters, as well as their genomic stability and integrity. Bovine viral diarrhea (BVD) negatively affects the health of dairy cattle and causes reproductive problems. The aim of the study was to identify genomic instability in cows with reproductive disorders following infection with the BVD virus. The material for analysis was peripheral blood from Holstein-Friesian cows with reproductive problems, which had tested positive for BVD, and from healthy cows with no reproductive problems, which had tested negative for BVD. Three cytogenetic tests were used: the sister chromatid exchange assay, fragile sites assay, and comet assay. Statistically significant differences were noted between the groups and between the individual cows in the average frequency of damage. The assays were good biomarkers of genomic stability and enabled the identification of individuals with an increased frequency of damage to genetic material that posed a negative impact on their health. The assays can be used to prevent disease during its course and evaluate the genetic resistance of animals. This is especially important for the breeder, both for economic and breeding reasons. Of the three assays, the comet assay proved to be the most sensitive for identifying DNA damage in the animals.

Replication stress in mammalian embryo development, differentiation, and reprogramming

Duplicating a genome of 3 billion nucleotides is challenged by a variety of obstacles that can cause replication stress and affect the integrity of the genome. Recent studies show that replication fork slowing and stalling is prevalent in early mammalian development, resulting in genome instability and aneuploidy, and constituting a barrier to development in human reproduction. Genome instability resulting from DNA replication stress is a barrier to the cloning of animals and to the reprogramming of differentiated cells to induced pluripotent stem cells, as well as a barrier to cell transformation. Remarkably, the regions most impacted by replication stress are shared in these different cellular contexts, affecting long genes and flanking intergenic areas. In this review we integrate our knowledge of DNA replication stress in mammalian embryos, in programming, and in reprogramming, and we discuss a potential role for fragile sites in sensing replication stress and restricting cell cycle progression in health and disease.

Mistimed origin licensing and activation stabilize common fragile sites under tight DNA-replication checkpoint activation

Genome integrity requires replication to be completed before chromosome segregation. The DNA-replication checkpoint (DRC) contributes to this coordination by inhibiting CDK1, which delays mitotic onset. Under-replication of common fragile sites (CFSs), however, escapes surveillance, resulting in mitotic chromosome breaks. Here we asked whether loose DRC activation induced by modest stresses commonly used to destabilize CFSs could explain this leakage. We found that tightening DRC activation or CDK1 inhibition stabilizes CFSs in human cells. Repli-Seq and molecular combing analyses showed a burst of replication initiations implemented in mid S-phase across a subset of late-replicating sequences, including CFSs, while the bulk genome was unaffected. CFS rescue and extra-initiations required CDC6 and CDT1 availability in S-phase, implying that CDK1 inhibition permits mistimed origin licensing and firing. In addition to delaying mitotic onset, tight DRC activation therefore supports replication completion of late origin-poor domains at risk of under-replication, two complementary roles preserving genome stability.

Genomic rearrangements and evolutionary changes in 3D chromatin topologies in the cotton tribe (Gossypieae)

BACKGROUND

Analysis of the relationship between chromosomal structural variation (synteny breaks) and 3D-chromatin architectural changes among closely related species has the potential to reveal causes and correlates between chromosomal change and chromatin remodeling. Of note, contrary to extensive studies in animal species, the pace and pattern of chromatin architectural changes following the speciation of plants remain unexplored; moreover, there is little exploration of the occurrence of synteny breaks in the context of multiple genome topological hierarchies within the same model species.

RESULTS

Here we used Hi-C and epigenomic analyses to characterize and compare the profiles of hierarchical chromatin architectural features in representative species of the cotton tribe (Gossypieae), including Gossypium arboreum, Gossypium raimondii, and Gossypioides kirkii, which differ with respect to chromosome rearrangements. We found that (i) overall chromatin architectural territories were preserved in Gossypioides and Gossypium, which was reflected in their similar intra-chromosomal contact patterns and spatial chromosomal distributions; (ii) the non-random preferential occurrence of synteny breaks in A compartment significantly associate with the B-to-A compartment switch in syntenic blocks flanking synteny breaks; (iii) synteny changes co-localize with open-chromatin boundaries of topologically associating domains, while TAD stabilization has a greater influence on regulating orthologous expression divergence than do rearrangements; and (iv) rearranged chromosome segments largely maintain ancestral in-cis interactions.

CONCLUSIONS

Our findings provide insights into the non-random occurrence of epigenomic remodeling relative to the genomic landscape and its evolutionary and functional connections to alterations of hierarchical chromatin architecture, on a known evolutionary timescale.

Effect of Selected Micro- and Macroelements and Vitamins on the Genome Stability of Bovine Embryo Transfer Recipients following In Vitro Fertilization

Genome instability can lead to a wide variety of diseases. Many endogenous and exogenous factors influence the level of damage to genetic material. Genome integrity depends on factors such as the fidelity of DNA replication, normal DNA organization in the chromosomes, and repair mechanisms. Genome stability influences fertility, embryonic development, and the maintenance of pregnancy. In the case of in vitro fertilization, it can be an important factor determining the success of the procedure. The aim of the study was to assess the stability of the genomes of recipient cows following in vitro fertilization using cytogenetic tests and to analyze the effects of selected vitamins and micro- and macroelements on genome integrity. Genome stability was analyzed using the sister chromatid exchange, fragile site, and comet assays. The material for analysis was peripheral blood from 20 Holstein-Friesian heifers that were embryo transfer recipients. The effect of selected micro- and macroelements and vitamins on the genome stability of the cows was analyzed. Folic acid was shown to significantly influence the level of damage identified using the SCE, FS, and SCGE assays, while iron affected SCE and SCGE results, and zinc affected FS.

GEN1 promotes common fragile site expression

Our genomes harbor conserved DNA sequences, known as common fragile sites (CFSs), that are difficult to replicate and correspond to regions of genome instability. Following replication stress, CFS loci give rise to breaks or gaps (termed CFS expression) where under-replicated DNA subsequently undergoes mitotic DNA synthesis (MiDAS). We show that loss of the structure-selective endonuclease GEN1 reduces CFS expression, leading to defects in MiDAS, ultrafine anaphase bridge formation, and DNA damage in the ensuing cell cycle due to aberrant chromosome segregation. GEN1 knockout cells also exhibit an elevated frequency of bichromatid constrictions consistent with the presence of unresolved regions of under-replicated DNA. Previously, the role of GEN1 was thought to be restricted to the nucleolytic resolution of recombination intermediates. However, its ability to cleave under-replicated DNA at CFS loci indicates that GEN1 plays a dual role resolving both DNA replication and recombination intermediates before chromosome segregation.

Potential Role of the Fragile Histidine Triad in Cancer Evo-Dev

Cancer development follows an evolutionary pattern of "mutation-selection-adaptation" detailed by Cancer Evolution and Development (Cancer Evo-Dev), a theory that represents a process of accumulating somatic mutations due to the imbalance between the mutation-promoting force and the mutation-repairing force and retro-differentiation of the mutant cells to cancer initiation cells in a chronic inflammatory microenvironment. The fragile histidine triad (FHIT) gene is a tumor suppressor gene whose expression is often reduced or inactivated in precancerous lesions during chronic inflammation or virus-induced replicative stress. Here, we summarize evidence regarding the mechanisms by which the FHIT is inactivated in cancer, including the loss of heterozygosity and the promoter methylation, and characterizes the role of the FHIT in bridging macroevolution and microevolution and in facilitating retro-differentiation during cancer evolution and development. It is suggested that decreased FHIT expression is involved in several critical steps of Cancer Evo-Dev. Future research needs to focus on the role and mechanisms of the FHIT in promoting the transformation of pre-cancerous lesions into cancer.

Computational workflow for integrative analyses of DNA replication timing, epigenomic, and transcriptomic data

Temporal profiling of DNA replication timing (RT) in combination with chromatin modifications, chromatin accessibility, and gene expression provides new insights into the causal relationships between chromatin and RT during cell cycle. Here, we describe a protocol for in-depth integrative computational analyses of Repli-seq, ATAC-seq, RNA-seq, and ChIP-seq or CUT&RUN data for multiple marks at various time points across cell cycle and changes in their interrelationships upon an experimental perturbation (e.g., knockdown or overexpression of a regulatory protein). For complete details on the use and execution of this protocol, please refer to Van Rechem et al. (2021).

3D chromatin connectivity underlies replication origin efficiency in mouse embryonic stem cells

In mammalian cells, chromosomal replication starts at thousands of origins at which replisomes are assembled. Replicative stress triggers additional initiation events from 'dormant' origins whose genomic distribution and regulation are not well understood. In this study, we have analyzed origin activity in mouse embryonic stem cells in the absence or presence of mild replicative stress induced by aphidicolin, a DNA polymerase inhibitor, or by deregulation of origin licensing factor CDC6. In both cases, we observe that the majority of stress-responsive origins are also active in a small fraction of the cell population in a normal S phase, and stress increases their frequency of activation. In a search for the molecular determinants of origin efficiency, we compared the genetic and epigenetic features of origins displaying different levels of activation, and integrated their genomic positions in three-dimensional chromatin interaction networks derived from high-depth Hi-C and promoter-capture Hi-C data. We report that origin efficiency is directly proportional to the proximity to transcriptional start sites and to the number of contacts established between origin-containing chromatin fragments, supporting the organization of origins in higher-level DNA replication factories.

Replication stress generates distinctive landscapes of DNA copy number alterations and chromosome scale losses

BACKGROUND

A major driver of cancer chromosomal instability is replication stress, the slowing or stalling of DNA replication. How replication stress and genomic instability are connected is not known. Aphidicolin-induced replication stress induces breakages at common fragile sites, but the exact causes of fragility are debated, and acute genomic consequences of replication stress are not fully explored.

RESULTS

We characterize DNA copy number alterations (CNAs) in single, diploid non-transformed cells, caused by one cell cycle in the presence of either aphidicolin or hydroxyurea. Multiple types of CNAs are generated, associated with different genomic regions and features, and observed copy number landscapes are distinct between aphidicolin and hydroxyurea-induced replication stress. Coupling cell type-specific analysis of CNAs to gene expression and single-cell replication timing analyses pinpointed the causative large genes of the most recurrent chromosome-scale CNAs in aphidicolin. These are clustered on chromosome 7 in RPE1 epithelial cells but chromosome 1 in BJ fibroblasts. Chromosome arm level CNAs also generate acentric lagging chromatin and micronuclei containing these chromosomes.

CONCLUSIONS

Chromosomal instability driven by replication stress occurs via focal CNAs and chromosome arm scale changes, with the latter confined to a very small subset of chromosome regions, potentially heavily skewing cancer genome evolution. Different inducers of replication stress lead to distinctive CNA landscapes providing the opportunity to derive copy number signatures of specific replication stress mechanisms. Single-cell CNA analysis thus reveals the impact of replication stress on the genome, providing insights into the molecular mechanisms which fuel chromosomal instability in cancer.

Structural variations in cancer and the 3D genome

Structural variations (SVs) affect more of the cancer genome than any other type of somatic genetic alteration but difficulties in detecting and interpreting them have limited our understanding. Clinical cancer sequencing also increasingly aims to detect SVs, leading to a widespread necessity to interpret their biological and clinical relevance. Recently, analyses of large whole-genome sequencing data sets revealed features that impact rates of SVs across the genome in different cancers. A striking feature has been the extent to which, in both their generation and their influence on the selective fitness of cancer cells, SVs are more specific to individual cancer types than other genetic alterations such as single-nucleotide variants. This Perspective discusses how the folding of the 3D genome, and differences in its folding across cell types, affect observed SV rates in different cancer types as well as how SVs can impact cancer cell fitness.

Replication stress impairs chromosome segregation and preimplantation development in human embryos

Human cleavage-stage embryos frequently acquire chromosomal aneuploidies during mitosis due to unknown mechanisms. Here, we show that S phase at the 1-cell stage shows replication fork stalling, low fork speed, and DNA synthesis extending into G2 phase. DNA damage foci consistent with collapsed replication forks, DSBs, and incomplete replication form in G2 in an ATR- and MRE11-dependent manner, followed by spontaneous chromosome breakage and segmental aneuploidies. Entry into mitosis with incomplete replication results in chromosome breakage, whole and segmental chromosome errors, micronucleation, chromosome fragmentation, and poor embryo quality. Sites of spontaneous chromosome breakage are concordant with sites of DNA synthesis in G2 phase, locating to gene-poor regions with long neural genes, which are transcriptionally silent at this stage of development. Thus, DNA replication stress in mammalian preimplantation embryos predisposes gene-poor regions to fragility, and in particular in the human embryo, to the formation of aneuploidies, impairing developmental potential.

Transcription and genome integrity

Transcription can cause genome instability by promoting R-loop formation but also act as a mutation-suppressing machinery by sensing of DNA lesions leading to the activation of DNA damage signaling and transcription-coupled repair. Recovery of RNA synthesis following the resolution of repair of transcription-blocking lesions is critical to avoid apoptosis and several new factors involved in this process have recently been identified. Some DNA repair proteins are recruited to initiating RNA polymerases and this may expediate the recruitment of other factors that participate in the repair of transcription-blocking DNA lesions. Recent studies have shown that transcription of protein-coding genes does not always give rise to spliced transcripts, opening the possibility that cells may use the transcription machinery in a splicing-uncoupled manner for other purposes including surveillance of the transcribed genome.

New Era of Mapping and Understanding Common Fragile Sites: An Updated Review on Origin of Chromosome Fragility

Common fragile sites (CFSs) are specific genomic loci prone to forming gaps or breakages upon replication perturbation, which correlate well with chromosomal rearrangement and copy number variation. CFSs have been actively studied due to their important pathophysiological relevance in different diseases such as cancer and neurological disorders. The genetic locations and sequences of CFSs are crucial to understanding the origin of such unstable sites, which require reliable mapping and characterizing approaches. In this review, we will inspect the evolving techniques for CFSs mapping, especially genome-wide mapping and sequencing of CFSs based on current knowledge of CFSs. We will also revisit the well-established hypotheses on the origin of CFSs fragility, incorporating novel findings from the comprehensive analysis of finely mapped CFSs regarding their locations, sequences, and replication/transcription, etc. This review will present the most up-to-date picture of CFSs and, potentially, a new framework for future research of CFSs.

Breaking the paradigm: early insights from mammalian DNA breakomes

DNA double-strand breaks (DSBs) can result from both exogenous and endogenous sources and are potentially toxic lesions to the human genome. If improperly repaired, DSBs can threaten genome integrity and contribute to premature ageing, neurodegenerative disorders and carcinogenesis. Through decades of work on genome stability, it has become evident that certain regions of the genome are inherently more prone to breakage than others, known as genome instability hotspots. Recent advancements in sequencing-based technologies now enable the profiling of genome-wide distributions of DSBs, also known as breakomes, to systematically map these instability hotspots. Here, we review the application of these technologies and their implications for our current understanding of the genomic regions most likely to drive genome instability. These breakomes ultimately highlight both new and established breakage hotspots including actively transcribed regions, loop boundaries and early-replicating regions of the genome. Further, these breakomes challenge the paradigm that DNA breakage primarily occurs in hard-to-replicate regions. With these advancements, we begin to gain insights into the biological mechanisms both invoking and protecting against genome instability.

Human topoisomerases and their roles in genome stability and organization

Human topoisomerases comprise a family of six enzymes: two type IB (TOP1 and mitochondrial TOP1 (TOP1MT), two type IIA (TOP2A and TOP2B) and two type IA (TOP3A and TOP3B) topoisomerases. In this Review, we discuss their biochemistry and their roles in transcription, DNA replication and chromatin remodelling, and highlight the recent progress made in understanding TOP3A and TOP3B. Because of recent advances in elucidating the high-order organization of the genome through chromatin loops and topologically associating domains (TADs), we integrate the functions of topoisomerases with genome organization. We also discuss the physiological and pathological formation of irreversible topoisomerase cleavage complexes (TOPccs) as they generate topoisomerase DNA–protein crosslinks (TOP-DPCs) coupled with DNA breaks. We discuss the expanding number of redundant pathways that repair TOP-DPCs, and the defects in those pathways, which are increasingly recognized as source of genomic damage leading to neurological diseases and cancer.

Topoisomerases have essential roles in transcription, DNA replication, chromatin remodelling and, as recently revealed, 3D genome organization. However, topoisomerases also generate DNA–protein crosslinks coupled with DNA breaks, which are increasingly recognized as a source of disease-causing genomic damage.

Characterization of Chromosomal Instability in Glioblastoma

Glioblastoma multiforme (GBM) is a malignant tumor of the central nervous system (CNS). The poor prognosis of GBM due to resistance to therapy has been associated with high chromosomal instability (CIN). Replication stress is a major cause of CIN that manifests as chromosome rearrangements, fragility, and breaks, including those cytologically expressed within specific chromosome regions named common fragile sites (CFSs). In this work, we characterized the expression of human CFSs in the glioblastoma U-251 MG cell line upon treatment with the inhibitor of DNA polymerase alpha aphidicolin (APH). We observed 52 gaps/breaks located within previously characterized CFSs. We found 17 to be CFSs in GBM cells upon treatment with APH, showing a frequency equal to at least 1% of the total gaps/breaks. We report that two CFSs localized to regions FRA2E (2p13/p12) and FRA2F (2q22) were only found in U-251 MG cells, but not lymphocytes or fibroblasts, after APH treatment. Notably, these glioblastoma-specific CFSs had a relatively high expression compared to the other CFSs with breakage frequency between ∼7 and 9%. Presence of long genes, incomplete replication, and delayed DNA synthesis during mitosis (MiDAS) after APH treatment suggest that an impaired replication process may contribute to this loci-specific fragility in U-251 MG cells. Altogether, our work offers a characterization of common fragile site expression in glioblastoma U-251 MG cells that may be further exploited for cytogenetic and clinical studies to advance our understanding of this incurable cancer.

Deciphering the Role of 3D Genome Organization in Breast Cancer Susceptibility

Cancer risk by environmental exposure is modulated by an individual's genetics and age at exposure. This age-specific period of susceptibility is referred to as the "Window of Susceptibility" (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

Epigenomic signatures associated with spontaneous and replication stress-induced DNA double strand breaks

Common fragile sites (CFSs) are specific regions of all individuals' genome that are predisposed to DNA double strand breaks (DSBs) and undergo subsequent rearrangements. CFS formation can be induced in vitro by mild level of DNA replication stress, such as DNA polymerase inhibition or nucleotide pool disturbance. The mechanisms of CFS formation have been linked to DNA replication timing control, transcription activities, as well as chromatin organization. However, it is unclear what specific cis- or trans-factors regulate the interplay between replication and transcription that determine CFS formation. We recently reported genome-wide mapping of DNA DSBs under replication stress induced by aphidicolin in human lymphoblastoids for the first time. Here, we systematically compared these DSBs with regards to nearby epigenomic features mapped in the same cell line from published studies. We demonstrate that aphidicolin-induced DSBs are strongly correlated with histone 3 lysine 36 trimethylation, a marker for active transcription. We further demonstrate that this DSB signature is a composite effect by the dual treatment of aphidicolin and its solvent, dimethylsulfoxide, the latter of which potently induces transcription on its own. We also present complementing evidence for the association between DSBs and 3D chromosome architectural domains with high density gene cluster and active transcription. Additionally, we show that while DSBs were detected at all but one of the fourteen finely mapped CFSs, they were not enriched in the CFS core sequences and rather demarcated the CFS core region. Related to this point, DSB density was not higher in large genes of greater than 300 kb, contrary to reported enrichment of CFS sites at these large genes. Finally, replication timing analyses demonstrate that the CFS core region contain initiation events, suggesting that altered replication dynamics are responsible for CFS formation in relatively higher level of replication stress.

OUP accepted manuscript

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription–replication interactions. Here, we report a novel method to identify genomic loci prone to transcription–replication interactions termed transcription–replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription–replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Camptothecin-Induced Replication Stress Affects DNA Replication Profiling by E/L Repli-Seq

E/L Repli-seq is a powerful tool for detecting cell type-specific replication landscapes in mammalian cells, but its potential to monitor DNA replication under replication stress awaits better understanding. Here, we used E/L Repli-seq to examine the temporal order of DNA replication in human retinal pigment epithelium cells treated with the topoisomerase I inhibitor camptothecin. We found that the replication profiles by E/L Repli-seq exhibit characteristic patterns after replication-stress induction, including the loss of specific initiation zones within individual early replication timing domains. We also observed global disappearance of the replication timing domain structures in the profiles, which can be explained by checkpoint-dependent suppression of replication initiation. Thus, our results demonstrate the effectiveness of E/L Repli-seq at identifying cells with replication-stress-induced altered DNA replication programs.

Machine learning meets omics: applications and perspectives

The innovation of biotechnologies has allowed the accumulation of omics data at an alarming rate, thus introducing the era of 'big data'. Extracting inherent valuable knowledge from various omics data remains a daunting problem in bioinformatics. Better solutions often need some kind of more innovative methods for efficient handlings and effective results. Recent advancements in integrated analysis and computational modeling of multi-omics data helped address such needs in an increasingly harmonious manner. The development and application of machine learning have largely advanced our insights into biology and biomedicine and greatly promoted the development of therapeutic strategies, especially for precision medicine. Here, we propose a comprehensive survey and discussion on what happened, is happening and will happen when machine learning meets omics. Specifically, we describe how artificial intelligence can be applied to omics studies and review recent advancements at the interface between machine learning and the ever-widest range of omics including genomics, transcriptomics, proteomics, metabolomics, radiomics, as well as those at the single-cell resolution. We also discuss and provide a synthesis of ideas, new insights, current challenges and perspectives of machine learning in omics.

Encounters in Three Dimensions: How Nuclear Topology Shapes Genome Integrity

Almost 25 years ago, the phosphorylation of a chromatin component, histone H2AX, was discovered as an integral part of the DNA damage response in eukaryotes. Much has been learned since then about the control of DNA repair in the context of chromatin. Recent technical and computational advances in imaging, biophysics and deep sequencing have led to unprecedented insight into nuclear organization, highlighting the impact of three-dimensional (3D) chromatin structure and nuclear topology on DNA repair. In this review, we will describe how DNA repair processes have adjusted to and in many cases adopted these organizational features to ensure accurate lesion repair. We focus on new findings that highlight the importance of chromatin context, topologically associated domains, phase separation and DNA break mobility for the establishment of repair-conducive nuclear environments. Finally, we address the consequences of aberrant 3D genome maintenance for genome instability and disease.

Programmed DNA elimination: silencing genes and repetitive sequences in somatic cells

In a multicellular organism, the genomes of all cells are in general the same. Programmed DNA elimination is a notable exception to this genome constancy rule. DNA elimination removes genes and repetitive elements in the germline genome to form a reduced somatic genome in various organisms. The process of DNA elimination within an organism is highly accurate and reproducible; it typically occurs during early embryogenesis, coincident with germline-soma differentiation. DNA elimination provides a mechanism to silence selected genes and repeats in somatic cells. Recent studies in nematodes suggest that DNA elimination removes all chromosome ends, resolves sex chromosome fusions, and may also promote the birth of novel genes. Programmed DNA elimination processes are diverse among species, suggesting DNA elimination likely has evolved multiple times in different taxa. The growing list of organisms that undergo DNA elimination indicates that DNA elimination may be more widespread than previously appreciated. These various organisms will serve as complementary and comparative models to study the function, mechanism, and evolution of programmed DNA elimination in metazoans.

Collective regulation of chromatin modifications predicts replication timing during cell cycle

Replication timing (RT) associates with genome architecture, while having a mixed relationship to histone marks. By profiling replication at high resolution and assessing broad histone marks across the cell cycle at the resolution of RT with and without genetic perturbation, we address the causal relationship between histone marks and RT. Four primary chromatin states, including an uncharacterized H3K36me2 state, emerge and define 97% of the mappable genome. RT and local replication patterns (e.g., initiation zones) quantitatively associate with chromatin states, histone mark dynamics, and spatial chromatin structure. Manipulation of broad histone marks and enhancer elements by overexpressing the histone H3 lysine 9/36 tri-demethylase KDM4A impacts RT across 11% of the genome. Broad histone modification changes were strong predictors of the observed RT alterations. Lastly, replication within H3K36me2-enriched neighborhoods is sensitive to KDM4A overexpression and is controlled at a megabase scale. These studies establish a role for collective chromatin mark regulation in modulating RT.

Common Threads: Aphidicolin-Inducible and Folate-Sensitive Fragile Sites in the Human Genome

The human genome has many chromosomal regions that are fragile, demonstrating chromatin breaks, gaps, or constrictions on exposure to replication stress. Common fragile sites (CFSs) are found widely distributed in the population, with the largest subset of these sites being induced by aphidicolin (APH). Other fragile sites are only found in a subset of the population. One group of these so-called rare fragile sites (RFSs) is induced by folate stress. APH-inducible CFSs are generally located in large transcriptionally active genes that are A + T rich and often enriched for tracts of AT-dinucleotide repeats. In contrast, all the folate-sensitive sites mapped to date consist of transcriptionally silenced CGG microsatellites. Thus, all the folate-sensitive fragile sites may have a very similar molecular basis that differs in key ways from that of the APH CFSs. The folate-sensitive FSs include FRAXA that is associated with Fragile X syndrome (FXS), the most common heritable form of intellectual disability. Both CFSs and RFSs can cause chromosomal abnormalities. Recent work suggests that both APH-inducible fragile sites and FRAXA undergo Mitotic DNA synthesis (MiDAS) when exposed to APH or folate stress, respectively. Interestingly, blocking MiDAS in both cases prevents chromosome fragility but increases the risk of chromosome mis-segregation. MiDAS of both APH-inducible and FRAXA involves conservative DNA replication and POLD3, an accessory subunit of the replicative polymerase Pol δ that is essential for break-induced replication (BIR). Thus, MiDAS is thought to proceed via some form of BIR-like process. This review will discuss the recent work that highlights the similarities and differences between these two groups of fragile sites and the growing evidence for the presence of many more novel fragile sites in the human genome.

Chromosomal Location of Genes Differentially Expressed in Tumor Cells Surviving High-Dose X-ray Irradiation: A Preliminary Study on Radio-Fragile Sites

Altered gene expression is a common feature of tumor cells after irradiation. Our previous study showed that this phenomenon is not only an acute response to cytotoxic stress, instead, it was persistently detected in tumor cells that survived 10 Gy irradiation (IR cells). The current understanding is that DNA double-strand breaks (DSBs) are recognized by the phosphorylation of histone H2AX (H2AX) and triggers the ataxia-telangiectasia mutated (ATM) protein or the ATM- and Rad3-related (ATR) pathway, which activate or inactivate the DNA repair or apoptotic or senescence related molecules and causes the expression of genes in many instances. However, because changes in gene expression persist after passaging in IR cells, it may be due to the different pathways from these transient intracellular signaling pathways caused by DSBs. We performed microarray analysis of 30,000 genes in radiation-surviving cells (H1299-IR and MCF7-IR) and found an interesting relation between altered genes and their chromosomal loci. These loci formed a cluster on the chromosome, especially on 1q21 and 6p21-p22 in both irradiated cell lines. These chromosome sites might be regarded as “radio-fragile” sites.

Transcription-Replication Collisions-A Series of Unfortunate Events

Transcription-replication interactions occur when DNA replication encounters genomic regions undergoing transcription. Both replication and transcription are essential for life and use the same DNA template making conflicts unavoidable. R-loops, DNA supercoiling, DNA secondary structure, and chromatin-binding proteins are all potential obstacles for processive replication or transcription and pose an even more potent threat to genome integrity when these processes co-occur. It is critical to maintaining high fidelity and processivity of transcription and replication while navigating through a complex chromatin environment, highlighting the importance of defining cellular pathways regulating transcription-replication interaction formation, evasion, and resolution. Here we discuss how transcription influences replication fork stability, and the safeguards that have evolved to navigate transcription-replication interactions and maintain genome integrity in mammalian cells.

Locus-specific transcription silencing at the FHIT gene suppresses replication stress-induced copy number variant formation and associated replication delay

Impaired replication progression leads to de novo copy number variant (CNV) formation at common fragile sites (CFSs). We previously showed that these hotspots for genome instability reside in late-replicating domains associated with large transcribed genes and provided indirect evidence that transcription is a factor in their instability. Here, we compared aphidicolin (APH)-induced CNV and CFS frequency between wild-type and isogenic cells in which FHIT gene transcription was ablated by promoter deletion. Two promoter-deletion cell lines showed reduced or absent CNV formation and CFS expression at FHIT despite continued instability at the NLGN1 control locus. APH treatment led to critical replication delays that remained unresolved in G2/M in the body of many, but not all, large transcribed genes, an effect that was reversed at FHIT by the promoter deletion. Altering RNase H1 expression did not change CNV induction frequency and DRIP-seq showed a paucity of R-loop formation in the central regions of large genes, suggesting that R-loops are not the primary mediator of the transcription effect. These results demonstrate that large gene transcription is a determining factor in replication stress-induced genomic instability and support models that CNV hotspots mainly result from the transcription-dependent passage of unreplicated DNA into mitosis.

Mind the replication gap

Unlike bacteria, mammalian cells need to complete DNA replication before segregating their chromosomes for the maintenance of genome integrity. Thus, cells have evolved efficient pathways to restore stalled and/or collapsed replication forks during S-phase, and when necessary, also to delay cell cycle progression to ensure replication completion. However, strong evidence shows that cells can proceed to mitosis with incompletely replicated DNA when under mild replication stress (RS) conditions. Consequently, the incompletely replicated genomic gaps form, predominantly at common fragile site regions, where the converging fork-like DNA structures accumulate. These branched structures pose a severe threat to the faithful disjunction of chromosomes as they physically interlink the partially duplicated sister chromatids. In this review, we provide an overview discussing how cells respond and deal with the under-replicated DNA structures that escape from the S/G2 surveillance system. We also focus on recent research of a mitotic break-induced replication pathway (also known as mitotic DNA repair synthesis), which has been proposed to operate during prophase in an attempt to finish DNA synthesis at the under-replicated genomic regions. Finally, we discuss recent data on how mild RS may cause chromosome instability and mutations that accelerate cancer genome evolution.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Cytolethal Distending Toxin Promotes Replicative Stress Leading to Genetic Instability Transmitted to Daughter Cells

The cytolethal distending toxin (CDT) is produced by several Gram-negative pathogenic bacteria. In addition to inflammation, experimental evidences are in favor of a protumoral role of CDT-harboring bacteria such as Escherichia coli, Campylobacter jejuni, or Helicobacter hepaticus. CDT may contribute to cell transformation in vitro and carcinogenesis in mice models, through the genotoxic action of CdtB catalytic subunit. Here, we investigate the mechanism of action by which CDT leads to genetic instability in human cell lines and colorectal organoids from healthy patients’ biopsies. We demonstrate that CDT holotoxin induces a replicative stress dependent on CdtB. The slowing down of DNA replication occurs mainly in late S phase, resulting in the expression of fragile sites and important chromosomic aberrations. These DNA abnormalities induced after CDT treatment are responsible for anaphase bridge formation in mitosis and interphase DNA bridge between daughter cells in G1 phase. Moreover, CDT-genotoxic potential preferentially affects human cycling cells compared to quiescent cells. Finally, the toxin induces nuclear distension associated to DNA damage in proliferating cells of human colorectal organoids, resulting in decreased growth. Our findings thus identify CDT as a bacterial virulence factor targeting proliferating cells, such as human colorectal progenitors or stem cells, inducing replicative stress and genetic instability transmitted to daughter cells that may therefore contribute to carcinogenesis. As some CDT-carrying bacterial strains were detected in patients with colorectal cancer, targeting these bacteria could be a promising therapeutic strategy.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Dealing with DNA lesions: When one cell cycle is not enough

Subversion of genome integrity fuels cellular adaptation and is a prerequisite for organismal evolution, yet genomic lesions are also the harmful driving force of cancer and other age-related human diseases. Genome integrity maintenance is inherently linked to genome organization and nuclear architecture, which are substantially remodeled during the cell cycle. Here we discuss recent findings on how actively dividing cells cope with endogenous genomic lesions that occur frequently at repetitive, heterochromatic, and late replicating regions as byproducts of genome duplication. We discuss how such lesions, rather than being resolved immediately when they occur, are dealt with in subsequent cell cycle phases, and even after mitotic cell division, and how this in turn affects genome organization, stability, and function.

Hotspots of Human Mutation

Mutation of the human genome results in three classes of genomic variation: single nucleotide variants; short insertions or deletions; and large structural variants (SVs). Some mutations occur during normal processes, such as meiotic recombination or B cell development, and others result from DNA replication or aberrant repair of breaks in sequence-specific contexts. Regardless of mechanism, mutations are subject to selection, and some hotspots can manifest in disease. Here, we discuss genomic regions prone to mutation, mechanisms contributing to mutation susceptibility, and the processes leading to their accumulation in normal and somatic genomes. With further, more accurate human genome sequencing, additional mutation hotspots, mechanistic details of their formation, and the relevance of hotspots to evolution and disease are likely to be discovered.

Chromosome structural variation in tumorigenesis: mechanisms of formation and carcinogenesis

With the rapid development of next-generation sequencing technology, chromosome structural variation has gradually gained increased clinical significance in tumorigenesis. However, the molecular mechanism(s) underlying this structural variation remain poorly understood. A search of the literature shows that a three-dimensional chromatin state plays a vital role in inducing structural variation and in the gene expression profiles in tumorigenesis. Structural variants may result in changes in copy number or deletions of coding sequences, as well as the perturbation of structural chromatin features, especially topological domains, and disruption of interactions between genes and their regulatory elements. This review focuses recent work aiming at elucidating how structural variations develop and misregulate oncogenes and tumor suppressors, to provide general insights into tumor formation mechanisms and to provide potential targets for future anticancer therapies.