A genome-wide analysis of common fragile sites: what features determine chromosomal instability in the human genome?

Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression

BACKGROUND

DNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression.

RESULTS

To characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density.

CONCLUSIONS

Altogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

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.

Association between the copy number variation of CCSER1 gene and growth traits in Chinese Capra hircus (goat) populations

Recently, Coiled-coil serine-rich protein 1 (CCSER1) gene is reported to be related to economic traits in livestock, and become a hotspot. In our study, we detected CCSER1 gene CNV in 693 goats from six breeds (GZB, GZW, AN, BH, HG, TH) by quantitative real-time PCR (qPCR) and the association analysis between the types of CNV and growth traits. Then, CCSER1 gene expression pattern was discovered in seven tissues from NB goats. Our results showed that the CCSER1 gene copy numbers were distributed differently in the aforementioned six breeds. The type of CCSER1 gene CNV was significantly associated with body weight and heart girth traits in GZW goat, in which individuals with deletion type were dominant in body weight trait (P < 0.05), while the normal type individuals were more advantageous in heart girth trait (P < 0.01); and there was a significant association with heart girth in TH goat (P < 0.05), which normal type was the dominant one. The expression profile revealed that CCSER1 gene has the highest level in the lung, followed by the small intestine and heart. In conclusion, our result is dedicated to an in-depth study of the novel CCSER1 gene CNV site and to provide essential information for Chinese goats molecular selective breeding in the future.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Completing genome replication outside of S phase

Mitotic DNA synthesis (MiDAS) is an unusual form of DNA replication that occurs during mitosis. Initially, MiDAS was characterized as a process associated with intrinsically unstable loci known as common fragile sites that occurs after cells experience DNA replication stress (RS). However, it is now believed to be a more widespread "salvage" mechanism that is called upon to complete the duplication of any under-replicated genomic region. Emerging data suggest that MiDAS is a DNA repair process potentially involving two or more pathways working in parallel or sequentially. In this review, we introduce the causes of RS, regions of the human genome known to be especially vulnerable to RS, and the strategies used to complete DNA replication outside of S phase. Additionally, because MiDAS is a prominent feature of aneuploid cancer cells, we will discuss how targeting MiDAS might potentially lead to improvements in cancer therapy.

TRAIP resolves DNA replication-transcription conflicts during the S-phase of unperturbed cells

Cell division is the basis for the propagation of life and requires accurate duplication of all genetic information. DNA damage created during replication (replication stress) is a major cause of cancer, premature aging and a spectrum of other human disorders. Over the years, TRAIP E3 ubiquitin ligase has been shown to play a role in various cellular processes that govern genome integrity and faultless segregation. TRAIP is essential for cell viability, and mutations in TRAIP ubiquitin ligase activity lead to primordial dwarfism in patients. Here, we have determined the mechanism of inhibition of cell proliferation in TRAIP-depleted cells. We have taken advantage of the auxin induced degron system to rapidly degrade TRAIP within cells and to dissect the importance of various functions of TRAIP in different stages of the cell cycle. We conclude that upon rapid TRAIP degradation, specifically in S-phase, cells cease to proliferate, arrest in G2 stage of the cell cycle and undergo senescence. Our findings reveal that TRAIP works in S-phase to prevent DNA damage at transcription start sites, caused by replication-transcription conflicts.

Dicentric chromosome breakage in Drosophila melanogaster is influenced by pericentric heterochromatin and occurs in nonconserved hotspots

Chromosome breakage plays an important role in the evolution of karyotypes and can produce deleterious effects within a single individual, such as aneuploidy or cancer. Forces that influence how and where chromosomes break are not fully understood. In humans, breakage tends to occur in conserved hotspots called common fragile sites (CFS), especially during replication stress. By following the fate of dicentric chromosomes in Drosophila melanogaster, we find that breakage under tension also tends to occur in specific hotspots. Our experimental approach was to induce sister chromatid exchange in a ring chromosome to generate a dicentric chromosome with a double chromatid bridge. In the following cell division, the dicentric bridges may break. We analyzed the breakage patterns of 3 different ring-X chromosomes. These chromosomes differ by the amount and quality of heterochromatin they carry as well as their genealogical history. For all 3 chromosomes, breakage occurs preferentially in several hotspots. Surprisingly, we found that the hotspot locations are not conserved between the 3 chromosomes: each displays a unique array of breakage hotspots. The lack of hotspot conservation, along with a lack of response to aphidicolin, suggests that these breakage sites are not entirely analogous to CFS and may reveal new mechanisms of chromosome fragility. Additionally, the frequency of dicentric breakage and the durability of each chromosome's spindle attachment vary significantly between the 3 chromosomes and are correlated with the origin of the centromere and the amount of pericentric heterochromatin. We suggest that different centromere strengths could account for this.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage^1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes^4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas⁸ (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus⁹ (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Different classes of genomic inserts contribute to human antibody diversity

The enormous diversity of antibodies is a key element to combat infections. Antibodies containing pathogen receptors were a surprising discovery that contrasted antibody diversification through classic recombination events. However, such insert-containing antibodies were thus far exclusively detected in African individuals exposed to malaria parasites and were identified as screening byproducts or through hypothesis-driven search. The prevalence and complexity of insertion events remained elusive. In this study, we devise an unbiased, systematic approach to identify inserts in the human antibody repertoire. We show that inserts from distant genomic regions occur in the majority of donors and are independent of Plasmodium falciparum preexposure. Our findings suggest that four distinct classes of insertion events contribute diversity to the human antibody repertoire.

Recombination of antibody genes in B cells can involve distant genomic loci and contribute a foreign antigen-binding element to form hybrid antibodies with broad reactivity for Plasmodium falciparum. So far, antibodies containing the extracellular domain of the LAIR1 and LILRB1 receptors represent unique examples of cross-chromosomal antibody diversification. Here, we devise a technique to profile non-VDJ elements from distant genes in antibody transcripts. Independent of the preexposure of donors to malaria parasites, non-VDJ inserts were detected in 80% of individuals at frequencies of 1 in 10⁴ to 10⁵ B cells. We detected insertions in heavy, but not in light chain or T cell receptor transcripts. We classify the insertions into four types depending on the insert origin and destination: 1) mitochondrial and 2) nuclear DNA inserts integrated at VDJ junctions; 3) inserts originating from telomere proximal genes; and 4) fragile sites incorporated between J-to-constant junctions. The latter class of inserts was exclusively found in memory and in in vitro activated B cells, while all other classes were already detected in naïve B cells. More than 10% of inserts preserved the reading frame, including transcripts with signs of antigen-driven affinity maturation. Collectively, our study unravels a mechanism of antibody diversification that is layered on the classical V(D)J and switch recombination.

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.

Extrachromosomal DNA in Cancer

In cancer, complex genome rearrangements and other structural alterations, including the amplification of oncogenes on circular extrachromosomal DNA (ecDNA) elements, drive the formation and progression of tumors. ecDNA is a particularly challenging structural alteration. By untethering oncogenes from chromosomal constraints, it elevates oncogene copy number, drives intratumoral genetic heterogeneity, promotes rapid tumor evolution, and results in treatment resistance. The profound changes in DNA shape and nuclear architecture generated by ecDNA alter the transcriptional landscape of tumors by catalyzing new types of regulatory interactions that do not occur on chromosomes. The current suite of tools for interrogating cancer genomes is well suited for deciphering sequence but has limited ability to resolve the complex changes in DNA structure and dynamics that ecDNA generates. Here, we review the challenges of resolving ecDNA form and function and discuss the emerging tool kit for deciphering ecDNA architecture and spatial organization, including what has been learned to date about how this dramatic change in shape alters tumor development, progression, and drug resistance.

Prognostic impact of chromosomal changes at relapse after allogeneic hematopoietic cell transplantation for acute myeloid leukemia or myelodysplastic syndrome

Chromosome analysis is a powerful prognostic tool in myeloid malignancies. Recipients who experience relapse after allogeneic hematopoietic cell transplantation (allo-HCT) often show chromosomal changes between diagnosis and relapse. However, the clinical impact of chromosomal changes and the efficacy of post-relapse treatment according to chromosomal changes have not been fully investigated. We retrospectively analyzed 72 recipients who had experienced relapse after allo-HCT for acute myeloid leukemia or myelodysplastic syndrome. We categorized them into two groups: with or without clonal chromosomal changes at relapse after allo-HCT. Post-relapse survival was shorter in the clonal chromosomal change group (median 117 days vs 275 days, P = 0.019). Moreover, acquisition of chromosome 7 abnormality or complex changes tended to be associated with inferior survival in a univariate analysis (median 92 days vs median 173 days, P = 0.043), and this adverse impact was confirmed in a multivariate analysis (hazard ratio 2.07, P = 0.024). The patterns of chromosomal changes from diagnosis to relapse after allo-HCT were heterogenous, and further investigations are required to clarify the effect of individual chromosomal changes.

First Comprehensive Characterization of Phayre’s Leaf-Monkey (Trachypithecus phayrei) Karyotype

The chromosomal homologies of human (Homo sapiens—HSA) and Trachypithecus phayrei (TPH—Phayre’s leaf-monkey, family Cercopithecidae) have previously been studied by using classical chromosome staining/banding and fluorescence in situ hybridization (FISH) from the 1970s to 1990s. In this study, we carried out molecular cytogenetics applying human multicolor banding (MCB), locus-specific, and human heterochromatin-specific probes to establish the first detailed chromosomal map of TPH, which was not available until now. Accordingly, it was possible to precisely determine evolutionary-conserved breakpoints (ECBs) and the orientation of evolutionary-conserved segments compared to HSA. It could be shown that five chromosomes remained completely unchanged between these two species, and 16 chromosomes underwent only intrachromosomal changes. In addition, 50 ECBs that failed to be resolved in previous reports were exactly identified and characterized in this study. It could also be shown that 43.5% of TPH centromere positions were conserved and 56.5% were altered compared to HSA. Interestingly, 82% ECBs in TPH corresponded to human fragile sites. Overall, this study is an essential contribution to future studies and reviews on chromosomal evolution in Cercopithecidae.

Integrated genome and transcriptome analyses reveal the mechanism of genome instability in ataxia with oculomotor apraxia 2

Mutations in the SETX gene, which encodes Senataxin, are associated with the progressive neurodegenerative diseases ataxia with oculomotor apraxia 2 (AOA2) and amyotrophic lateral sclerosis 4 (ALS4). To identify the causal defect in AOA2, patient-derived cells and SETX knockouts (human and mouse) were analyzed using integrated genomic and transcriptomic approaches. A genome-wide increase in chromosome instability (gains and losses) within genes and at chromosome fragile sites was observed, resulting in changes to gene-expression profiles. Transcription stress near promoters correlated with high GCskew and the accumulation of R-loops at promoter-proximal regions, which localized with chromosomal regions where gains and losses were observed. In the absence of Senataxin, the Cockayne syndrome protein CSB was required for the recruitment of the transcription-coupled repair endonucleases (XPG and XPF) and RAD52 recombination protein to target and resolve transcription bubbles containing R-loops, leading to genomic instability. These results show that transcription stress is an important contributor to SETX mutation-associated chromosome fragility and AOA2.

FusionAI: Predicting fusion breakpoint from DNA sequence with deep learning

Identifying the molecular mechanisms related to genomic breakage is an important goal of cancer mechanism studies. Among diverse locations of structural variants, fusion genes, which have the breakpoints in the gene bodies and are typically identified from the split reads of RNA-seq data, can provide a highlighted structural variant resource for studying the genomic breakages with expression and potential pathogenic impacts. In this study, we developed FusionAI, which utilizes deep learning to predict gene fusion breakpoints based on DNA sequence and let us identify fusion breakage code and genomic context. FusionAI leverages the known fusion breakpoints to provide a prediction model of the fusion genes from the primary genomic sequences via deep learning, thereby helping researchers a more accurate selection of fusion genes and better understand genomic breakage.

•

FusionAI predicts fusion gene breakpoints from a DNA sequence

•

FusonAI reduce the effort for validating fusion genes with other tools

•

High feature importance regions were apart 100nt from the exon junction BPs

•

High feature importance regions were overlapped with 44 human genomic features

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.

Genomic Mosaicism Formed by Somatic Variation in the Aging and Diseased Brain

Over the past 20 years, analyses of single brain cell genomes have revealed that the brain is composed of cells with myriad distinct genomes: the brain is a genomic mosaic, generated by a host of DNA sequence-altering processes that occur somatically and do not affect the germline. As such, these sequence changes are not heritable. Some processes appear to occur during neurogenesis, when cells are mitotic, whereas others may also function in post-mitotic cells. Here, we review multiple forms of DNA sequence alterations that have now been documented: aneuploidies and aneusomies, smaller copy number variations (CNVs), somatic repeat expansions, retrotransposons, genomic cDNAs (gencDNAs) associated with somatic gene recombination (SGR), and single nucleotide variations (SNVs). A catch-all term of DNA content variation (DCV) has also been used to describe the overall phenomenon, which can include multiple forms within a single cell’s genome. A requisite step in the analyses of genomic mosaicism is ongoing technology development, which is also discussed. Genomic mosaicism alters one of the most stable biological molecules, DNA, which may have many repercussions, ranging from normal functions including effects of aging, to creating dysfunction that occurs in neurodegenerative and other brain diseases, most of which show sporadic presentation, unlinked to causal, heritable genes.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Common Fragile Sites Are Characterized by Faulty Condensin Loading after Replication Stress

Cells coordinate interphase-to-mitosis transition, but recurrent cytogenetic lesions appear at common fragile sites (CFSs), termed CFS expression, in a tissue-specific manner after replication stress, marking regions of instability in cancer. Despite such a distinct defect, no model fully provides a molecular explanation for CFSs. We show that CFSs are characterized by impaired chromatin folding, manifesting as disrupted mitotic structures visible with molecular fluorescence in situ hybridization (FISH) probes in the presence and absence of replication stress. Chromosome condensation assays reveal that compaction-resistant chromatin lesions persist at CFSs throughout the cell cycle and mitosis. Cytogenetic and molecular lesions are marked by faulty condensin loading at CFSs, a defect in condensin-I-mediated compaction, and are coincident with mitotic DNA synthesis (MIDAS). This model suggests that, in conditions of exogenous replication stress, aberrant condensin loading leads to molecular defects and CFS expression, concomitantly providing an environment for MIDAS, which, if not resolved, results in chromosome instability.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Implications of Endogenous Retroelements in the Etiopathogenesis of Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. While its etiology remains elusive, current understanding suggests a multifactorial process with contributions by genetic, immunologic, hormonal, and environmental factors. A hypothesis that combines several of these factors proposes that genomic elements, the L1 retrotransposons, are instrumental in SLE pathogenesis. L1 retroelements are transcriptionally activated in SLE and produce two proteins, ORF1p and ORF2p, which are immunogenic and can drive type I interferon (IFN) production by producing DNA species that activate cytosolic DNA sensors. In addition, these two proteins reside in RNA-rich macromolecular assemblies that also contain well-known SLE autoantigens like Ro60. We surmise that cells expressing L1 will exhibit all the hallmarks of cells infected by a virus, resulting in a cellular and humoral immune response similar to those in chronic viral infections. However, unlike exogenous viruses, L1 retroelements cannot be eliminated from the host genome. Hence, dysregulated L1 will cause a chronic, but perhaps episodic, challenge for the immune system. The clinical and immunological features of SLE can be at least partly explained by this model. Here we review the support for, and the gaps in, this hypothesis of SLE and its potential for new diagnostic, prognostic, and therapeutic options in SLE.

SMC5/6 is required for replication fork stability and faithful chromosome segregation during neurogenesis

Mutations of SMC5/6 components cause developmental defects, including primary microcephaly. To model neurodevelopmental defects, we engineered a mouse wherein Smc5 is conditionally knocked out (cKO) in the developing neocortex. Smc5 cKO mice exhibited neurodevelopmental defects due to neural progenitor cell (NPC) apoptosis, which led to reduction in cortical layer neurons. Smc5 cKO NPCs formed DNA bridges during mitosis and underwent chromosome missegregation. SMC5/6 depletion triggers a CHEK2-p53 DNA damage response, as concomitant deletion of the Trp53 tumor suppressor or Chek2 DNA damage checkpoint kinase rescued Smc5 cKO neurodevelopmental defects. Further assessment using Smc5 cKO and auxin-inducible degron systems demonstrated that absence of SMC5/6 leads to DNA replication stress at late-replicating regions such as pericentromeric heterochromatin. In summary, SMC5/6 is important for completion of DNA replication prior to entering mitosis, which ensures accurate chromosome segregation. Thus, SMC5/6 functions are critical in highly proliferative stem cells during organism development.

Human papilloma virus (HPV) integration signature in Cervical Cancer: identification of MACROD2 gene as HPV hot spot integration site

Background

Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality with infection by human papilloma virus (HPV) being the most important risk factor. We analysed the association between different viral integration signatures, clinical parameters and outcome in pre-treated CCs.

Methods

Different integration signatures were identified using HPV double capture followed by next-generation sequencing (NGS) in 272 CC patients from the BioRAIDs study [NCT02428842]. Correlations between HPV integration signatures and clinical, biological and molecular features were assessed.

Results

Episomal HPV was much less frequent in CC as compared to anal carcinoma (p < 0.0001). We identified >300 different HPV-chromosomal junctions (inter- or intra-genic). The most frequent integration site in CC was in MACROD2 gene followed by MIPOL1/TTC6 and TP63. HPV integration signatures were not associated with histological subtype, FIGO staging, treatment or PFS. HPVs were more frequently episomal in PIK3CA mutated tumours (p = 0.023). Viral integration type was dependent on HPV genotype (p < 0.0001); HPV18 and HPV45 being always integrated. High HPV copy number was associated with longer PFS (p = 0.011).

Conclusions

This is to our knowledge the first study assessing the prognostic value of HPV integration in a prospectively annotated CC cohort, which detects a hotspot of HPV integration at MACROD2; involved in impaired PARP1 activity and chromosome instability.

Pervasive chromosomal instability and karyotype order in tumour evolution

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Unveiling Human Non-Random Genome Editing Mechanisms Activated in Response to Chronic Environmental Changes: I. Where Might These Mechanisms Come from and What Might They Have Led To?

This article challenges the notion of the randomness of mutations in eukaryotic cells by unveiling stress-induced human non-random genome editing mechanisms. To account for the existence of such mechanisms, I have developed molecular concepts of the cell environment and cell environmental stressors and, making use of a large quantity of published data, hypothesised the origin of some crucial biological leaps along the evolutionary path of life on Earth under the pressure of natural selection, in particular, (1) virus-cell mating as a primordial form of sexual recombination and symbiosis; (2) Lamarckian CRISPR-Cas systems; (3) eukaryotic gene development; (4) antiviral activity of retrotransposon-guided mutagenic enzymes; and finally, (5) the exaptation of antiviral mutagenic mechanisms to stress-induced genome editing mechanisms directed at "hyper-transcribed" endogenous genes. Genes transcribed at their maximum rate (hyper-transcribed), yet still unable to meet new chronic environmental demands generated by "pollution", are inadequate and generate more and more intronic retrotransposon transcripts. In this scenario, RNA-guided mutagenic enzymes (e.g., Apolipoprotein B mRNA editing catalytic polypeptide-like enzymes, APOBECs), which have been shown to bind to retrotransposon RNA-repetitive sequences, would be surgically targeted by intronic retrotransposons on opened chromatin regions of the same "hyper-transcribed" genes. RNA-guided mutagenic enzymes may therefore "Lamarkianly" generate single nucleotide polymorphisms (SNP) and gene copy number variations (CNV), as well as transposon transposition and chromosomal translocations in the restricted areas of hyper-functional and inadequate genes, leaving intact the rest of the genome. CNV and SNP of hyper-transcribed genes may allow cells to surgically explore a new fitness scenario, which increases their adaptability to stressful environmental conditions. Like the mechanisms of immunoglobulin somatic hypermutation, non-random genome editing mechanisms may generate several cell mutants, and those codifying for the most environmentally adequate proteins would have a survival advantage and would therefore be Darwinianly selected. Non-random genome editing mechanisms represent tools of evolvability leading to organismal adaptation including transgenerational non-Mendelian gene transmission or to death of environmentally inadequate genomes. They are a link between environmental changes and biological novelty and plasticity, finally providing a molecular basis to reconcile gene-centred and "ecological" views of evolution.

Genome-wide high-resolution mapping of mitotic DNA synthesis sites and common fragile sites by direct sequencing

Common fragile sites (CFSs) are genomic loci prone to the formation of breaks or gaps on metaphase chromosomes. They are hotspots for chromosome rearrangements and structural variations, which have been extensively implicated in carcinogenesis, aging, and other pathological processes. Although many CFSs were identified decades ago, a consensus is still lacking for why they are particularly unstable and sensitive to replication perturbations. This is in part due to the lack of high-resolution mapping data for the vast majority of the CFSs, which has hindered mechanistic interrogations. Here, we seek to map human CFSs with high resolution on a genome-wide scale by sequencing the sites of mitotic DNA synthesis (MiDASeq) that are specific for CFSs. We generated a nucleotide-resolution atlas of MiDAS sites (MDSs) that covered most of the known CFSs, and comprehensively analyzed their sequence characteristics and genomic features. Our data on MDSs tallied well with long-standing hypotheses to explain CFS fragility while highlighting the contributions of late replication timing and large transcription units. Notably, the MDSs also encompassed most of the recurrent double-strand break clusters previously identified in mouse neural stem/progenitor cells, thus bridging evolutionarily conserved break points across species. Moreover, MiDAseq provides an important resource that can stimulate future research on CFSs to further unravel the mechanisms and biological relevance underlying these labile genomic regions.

RTEL1 suppresses G-quadruplex-associated R-loops at difficult-to-replicate loci in the human genome

Oncogene activation during tumorigenesis generates DNA replication stress, a known driver of genome rearrangements. In response to replication stress, certain loci, such as common fragile sites and telomeres, remain under-replicated during interphase and subsequently complete locus duplication in mitosis in a process known as 'MiDAS'. Here, we demonstrate that RTEL1 (regulator of telomere elongation helicase 1) has a genome-wide role in MiDAS at loci prone to form G-quadruplex-associated R-loops, in a process that is dependent on its helicase function. We reveal that SLX4 is required for the timely recruitment of RTEL1 to the affected loci, which in turn facilitates recruitment of other proteins required for MiDAS, including RAD52 and POLD3. Our findings demonstrate that RTEL1 is required for MiDAS and suggest that RTEL1 maintains genome stability by resolving conflicts that can arise between the replication and transcription machineries.

Comprehensive analysis of indels in whole-genome microsatellite regions and microsatellite instability across 21 cancer types

Microsatellites are repeats of 1- to 6-bp units, and approximately 10 million microsatellites have been identified across the human genome. Microsatellites are vulnerable to DNA mismatch errors and have thus been used to detect cancers with mismatch repair deficiency. To reveal the mutational landscape of microsatellite repeat regions at the genome level, we analyzed approximately 20.1 billion microsatellites in 2717 whole genomes of pan-cancer samples across 21 tissue types. First, we developed a new insertion and deletion caller (MIMcall) that takes into consideration the error patterns of different types of microsatellites. Among the 2717 pan-cancer samples, our analysis identified 31 samples, including colorectal, uterus, and stomach cancers, with a higher proportion of mutated microsatellite (≥0.03), which we defined as microsatellite instability (MSI) cancers of genome-wide level. Next, we found 20 highly mutated microsatellites that can be used to detect MSI cancers with high sensitivity. Third, we found that replication timing and DNA shape were significantly associated with mutation rates of microsatellites. Last, analysis of mutations in mismatch repair genes showed that somatic SNVs and short indels had larger functional impacts than germline mutations and structural variations. Our analysis provides a comprehensive picture of mutations in the microsatellite regions and reveals possible causes of mutations, as well as provides a useful marker set for MSI detection.

The Landscape of Somatic Copy Number Alterations in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide. Somatic copy number alterations (CNAs) play a significant role in the development of this lethal cancer. In this study, we present a meta-analysis of CNAs for a total of 1,395 HNSCC samples. Publicly available R packages and in-house scripts were used for genomic array data processing, including normalization, segmentation and CNA calling. We detected 125 regions of significant gains or losses using GISTIC algorithm and found several potential driver genes in these regions. The incidence of chromothripsis in HNSCC was estimated to be 6%, and the chromosome pulverization hotspot regions were detected. We determined 323 genomic locations significantly enriched for breakpoints, which indicate HNSCC-specific genomic instability regions. Unsupervised clustering of genome-wide CNA data revealed a sub-cluster predominantly composed of nasopharynx tumors and presented a large proportion of HPV-positive samples. These results will facilitate the discovery of therapeutic candidates and extend our molecular understanding of HNSCC.

Immuno-genomic landscape of osteosarcoma

Limited clinical activity has been seen in osteosarcoma (OS) patients treated with immune checkpoint inhibitors (ICI). To gain insights into the immunogenic potential of these tumors, we conducted whole genome, RNA, and T-cell receptor sequencing, immunohistochemistry and reverse phase protein array profiling (RPPA) on OS specimens from 48 pediatric and adult patients with primary, relapsed, and metastatic OS. Median immune infiltrate level was lower than in other tumor types where ICI are effective, with concomitant low T-cell receptor clonalities. Neoantigen expression in OS was lacking and significantly associated with high levels of nonsense-mediated decay (NMD). Samples with low immune infiltrate had higher number of deleted genes while those with high immune infiltrate expressed higher levels of adaptive resistance pathways. PARP2 expression levels were significantly negatively associated with the immune infiltrate. Together, these data reveal multiple immunosuppressive features of OS and suggest immunotherapeutic opportunities in OS patients.

The efficacy of immune checkpoint inhibitors (ICI) in osteosarcoma has been limited. Here, the authors investigate the immunogenomic landscape of osteosarcoma, and integrated analyses highlight features related to a suppressed immune microenvironment.

BRCA1 and S phase DNA repair pathways restrict LINE-1 retrotransposition in human cells

Long interspersed element-1 (LINE-1 or L1) is the only autonomous retrotransposon active in human cells. Different host factors have been shown to influence L1 mobility however, systematic analyses of these factors are limited. Here, we developed a high-throughput microscopy-based retrotransposition assay that identified the Double-Stranded Break (DSB) repair and Fanconi Anemia factors active in the S/G2 phase as potent inhibitors and regulators of L1 activity. In particular BRCA1, an E3 ubiquitin ligase with a key role in several DNA repair pathways, directly affects L1 retrotransposition frequency and structure and also plays a distinct role in controlling L1 ORF2 protein translation through L1 mRNA binding. These results suggest the existence of a “battleground” at the DNA replication fork between HR factors and L1 retrotransposons, and revealing a potential role for L1 in the genotypic evolution of tumors characterized by BRCA1 and HR repair deficiencies.

Transcriptional effects of copy number alterations in a large set of human cancers

Copy number alterations (CNAs) can promote tumor progression by altering gene expression levels. Due to transcriptional adaptive mechanisms, however, CNAs do not always translate proportionally into altered expression levels. By reanalyzing >34,000 gene expression profiles, we reveal the degree of transcriptional adaptation to CNAs in a genome-wide fashion, which strongly associate with distinct biological processes. We then develop a platform-independent method-transcriptional adaptation to CNA profiling (TACNA profiling)-that extracts the transcriptional effects of CNAs from gene expression profiles without requiring paired CNA profiles. By applying TACNA profiling to >28,000 patient-derived tumor samples we define the landscape of transcriptional effects of CNAs. The utility of this landscape is demonstrated by the identification of four genes that are predicted to be involved in tumor immune evasion when transcriptionally affected by CNAs. In conclusion, we provide a novel tool to gain insight into how CNAs drive tumor behavior via altered expression levels.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

Autophagy gene haploinsufficiency drives chromosome instability, increases migration, and promotes early ovarian tumors

Autophagy, particularly with BECN1, has paradoxically been highlighted as tumor promoting in Ras-driven cancers, but potentially tumor suppressing in breast and ovarian cancers. However, studying the specific role of BECN1 at the genetic level is complicated due to its genomic proximity to BRCA1 on both human (chromosome 17) and murine (chromosome 11) genomes. In human breast and ovarian cancers, the monoallelic deletion of these genes is often co-occurring. To investigate the potential tumor suppressor roles of two of the most commonly deleted autophagy genes in ovarian cancer, BECN1 and MAP1LC3B were knocked-down in atypical (BECN1+/+ and MAP1LC3B+/+) ovarian cancer cells. Ultra-performance liquid chromatography mass-spectrometry metabolomics revealed reduced levels of acetyl-CoA which corresponded with elevated levels of glycerophospholipids and sphingolipids. Migration rates of ovarian cancer cells were increased upon autophagy gene knockdown. Genomic instability was increased, resulting in copy-number alteration patterns which mimicked high grade serous ovarian cancer. We further investigated the causal role of Becn1 haploinsufficiency for oncogenesis in a MISIIR SV40 large T antigen driven spontaneous ovarian cancer mouse model. Tumors were evident earlier among the Becn1+/- mice, and this correlated with an increase in copy-number alterations per chromosome in the Becn1+/- tumors. The results support monoallelic loss of BECN1 as permissive for tumor initiation and potentiating for genomic instability in ovarian cancer.

Investigation of somatic CNVs in brains of synucleinopathy cases using targeted SNCA analysis and single cell sequencing

Synucleinopathies are mostly sporadic neurodegenerative disorders of partly unexplained aetiology, and include Parkinson's disease (PD) and multiple system atrophy (MSA). We have further investigated our recent finding of somatic SNCA (α-synuclein) copy number variants (CNVs, specifically gains) in synucleinopathies, using Fluorescent in-situ Hybridisation for SNCA, and single-cell whole genome sequencing for the first time in a synucleinopathy. In the cingulate cortex, mosaicism levels for SNCA gains were higher in MSA and PD than controls in neurons (> 2% in both diseases), and for MSA also in non-neurons. In MSA substantia nigra (SN), we noted SNCA gains in > 3% of dopaminergic (DA) neurons (identified by neuromelanin) and neuromelanin-negative cells, including olig2-positive oligodendroglia. Cells with CNVs were more likely to have α-synuclein inclusions, in a pattern corresponding to cell categories mostly relevant to the disease: DA neurons in Lewy-body cases, and other cells in the striatonigral degeneration-dominant MSA variant (MSA-SND). Higher mosaicism levels in SN neuromelanin-negative cells may correlate with younger onset in typical MSA-SND, and in cingulate neurons with younger death in PD. Larger sample sizes will, however, be required to confirm these putative findings. We obtained genome-wide somatic CNV profiles from 169 cells from the substantia nigra of two MSA cases, and pons and putamen of one. These showed somatic CNVs in ~ 30% of cells, with clonality and origins in segmental duplications for some. CNVs had distinct profiles based on cell type, with neurons having a mix of gains and losses, and other cells having almost exclusively gains, although control data sets will be required to determine possible disease relevance. We propose that somatic SNCA CNVs may contribute to the aetiology and pathogenesis of synucleinopathies, and that genome-wide somatic CNVs in MSA brain merit further study.

Comprehensive identification and characterization of somatic copy number alterations in triple‑negative breast cancer

Triple‑negative breast cancer (TNBC) accounts for ~15% of all breast cancer diagnoses each year. Patients with TNBC tend to have a higher risk for early relapse and a worse prognosis. TNBC is characterized by extensive somatic copy number alterations (CNAs). However, the DNA CNA profile of TNBC remains to be extensively investigated. The present study assessed the genomic profile of CNAs in 201 TNBC samples, aiming to identify recurrent CNAs that may drive the pathogenesis of TNBC. In total, 123 regions of significant amplification and deletion were detected using the Genomic Identification of Significant Targets in Cancer algorithm, and potential driver genes for TNBC were identified. A total of 31 samples exhibited signs of chromothripsis and revealed chromosome pulverization hotspot regions. The present study further determined 199 genomic locations that were significantly enriched for breakpoints, which indicated TNBC‑specific genomic instability regions. Unsupervised hierarchical clustering of tumors resulted in three main subgroups that exhibited distinct CNA profiles, which may reveal the heterogeneity of molecular mechanisms in TNBC subgroups. These results will extend the molecular understanding of TNBC and will facilitate the discovery of therapeutic and diagnostic target candidates.

Higher Order Organization of the mtDNA: Beyond Mitochondrial Transcription Factor A

The higher order organization of eukaryotic and prokaryotic genomes is pivotal in the regulation of gene expression. Specifically, chromatin accessibility in eukaryotes and nucleoid accessibility in bacteria are regulated by a cohort of proteins to alter gene expression in response to diverse physiological conditions. By contrast, prior studies have suggested that the mitochondrial genome (mtDNA) is coated solely by mitochondrial transcription factor A (TFAM), whose increased cellular concentration was proposed to be the major determinant of mtDNA packaging in the mitochondrial nucleoid. Nevertheless, recent analysis of DNase-seq and ATAC-seq experiments from multiple human and mouse samples suggest gradual increase in mtDNA occupancy during the course of embryonic development to generate a conserved footprinting pattern which correlate with sites that have low TFAM occupancy in vivo (ChIP-seq) and tend to adopt G-quadruplex structures. These findings, along with recent identification of mtDNA binding by known modulators of chromatin accessibility such as MOF, suggest that mtDNA higher order organization is generated by cross talk with the nuclear regulatory system, may have a role in mtDNA regulation, and is more complex than once thought.

Transcription-mediated organization of the replication initiation program across large genes sets common fragile sites genome-wide

Common fragile sites (CFSs) are chromosome regions prone to breakage upon replication stress known to drive chromosome rearrangements during oncogenesis. Most CFSs nest in large expressed genes, suggesting that transcription could elicit their instability; however, the underlying mechanisms remain elusive. Genome-wide replication timing analyses here show that stress-induced delayed/under-replication is the hallmark of CFSs. Extensive genome-wide analyses of nascent transcripts, replication origin positioning and fork directionality reveal that 80% of CFSs nest in large transcribed domains poor in initiation events, replicated by long-travelling forks. Forks that travel long in late S phase explains CFS replication features, whereas formation of sequence-dependent fork barriers or head-on transcription–replication conflicts do not. We further show that transcription inhibition during S phase, which suppresses transcription–replication encounters and prevents origin resetting, could not rescue CFS stability. Altogether, our results show that transcription-dependent suppression of initiation events delays replication of large gene bodies, committing them to instability.

Common Fragile Sites (CFSs) are chromosome regions prone to breakage upon replication stress known to drive chromosome rearrangements during oncogenesis. Here the authors use genome-wide and single cell techniques to assess how replication timing and transcriptional activity correlate with genome stability.

RNase H2-RED carpets the path to eukaryotic RNase H2 functions

Eukaryotic RNases H2 have dual functions in initiating the removal of ribonucleoside monophosphates (rNMPs) incorporated by DNA polymerases during DNA synthesis and in cleaving the RNA moiety of RNA/DNA hybrids formed during transcription and retrotransposition. The other major cellular RNase H, RNase H1, shares the hybrid processing activity, but not all substrates. After RNase H2 incision at the rNMPs in DNA the Ribonucleotide Excision Repair (RER) pathway completes the removal, restoring dsDNA. The development of the RNase H2-RED (Ribonucleotide Excision Defective) mutant enzyme, which can process RNA/DNA hybrids but is unable to cleave rNMPs embedded in DNA has unlinked the two activities and illuminated the roles of RNase H2 in cellular metabolism. Studies mostly in Saccharomyces cerevisiae, have shown both activities of RNase H2 are necessary to maintain genome integrity and that RNase H1 and H2 have overlapping as well as distinct RNA/DNA hybrid substrates. In mouse RNase H2-RED confirmed that rNMPs in DNA during embryogenesis induce lethality in a p53-dependent DNA damage response. In mammalian cell cultures, RNase H2-RED helped identifying DNA lesions produced by Top1 cleavage at rNMPs and led to determine that RNase H2 participates in the retrotransposition of LINE-1 elements. In this review, we summarize the studies and conclusions reached by utilization of RNase H2-RED enzyme in different model systems.

Impact of AAV2 and Hepatitis B Virus Integration Into Genome on Development of Hepatocellular Carcinoma in Patients with Prior Hepatitis B Virus Infection

PURPOSE

Hepatitis B viral (HBV) DNA is frequently integrated into the genomes of hepatocellular carcinoma (HCC) in patients with chronic HBV infection (chronic HBV, hereafter), whereas the frequency of HBV integration in patients after the disappearance of HBV (prior HBV, hereafter) has yet to be determined. This study aimed to detect integration of HBV and adeno-associated virus type 2 (AAV2) into the human genome as a possible oncogenic event.

EXPERIMENTAL DESIGN

Virome capture sequencing was performed, using HCC and liver samples obtained from 243 patients, including 73 with prior HBV without hepatitis C viral (HCV) infection and 81 with chronic HBV.

RESULTS

Clonal HBV integration events were identified in 11 (15.0%) cases of prior HBV without HCV and 61 (75.3%) cases of chronic HBV (P < 0.001). Several driver genes were commonly targeted by HBV, leading to transcriptional activation of these genes; TERT [four (5.4%) vs. 15 (18.5%)], KMT2B [two (2.7%) vs. five (6.1%)], CCNE1 [zero vs. one (1.2%)], CCNA2 [zero vs. one (1.2%)]. Conversely, CCNE1 and CCNA2 were, respectively, targeted by AAV2 only in prior HBV. In liver samples, HBV genome recurrently integrated into fibrosis-related genes FN1, HS6ST3, KNG1, and ROCK1 in chronic HBV. There was not history of alcohol abuse and 3 patients with a history of nucleoside analogue treatment for HBV in 8 prior HBV with driver gene integration.

CONCLUSIONS

Despite the seroclearance of hepatitis B surface antigen, HBV or AAV2 integration in prior HBV was not rare; therefore, such patients are at risk of developing HCC.

Radiation-Induced Chromosomal Breaks may be DNA Repair Fragile Sites with Larger-scale Correlations to Eight Double-Strand-Break Related Data Sets over the Human Genome

In this work, we compared the genomic distribution of common radiation-induced chromosomal breaks to eight different data sets covering the whole human genome. Sites with a high probability of chromatid breakage after exposure to low and high ionization density radiations were often located inside common and rare fragile sites, indicating that they may be a new and more local type of DNA repair-related fragility. Breaks in specific chromosome bands after acute exposure to oil and benzene also showed strong correlation with these sites and fragile sites. In addition, close correlation was found with cytologically detected chiasma and MLH1 immunofluorescence sites and with the HapMap recombination density distributions. Also, of interest, copy number changes occurred predominantly at radiation-induced breaks and fragile sites, at least for breast cancers with poor prognosis, and they decreased weakly but significantly in regions with increasing recombination and CpG density. An increased CpG density is linked to regions of high gene density to secure high-fidelity reproduction and survival. To minimize cancer induction, cancer-related genes are often located in regions of decreased recombination density and/or higher-than-average CpG density. It is compelling that all these data sets were influenced by the cells' handling of double-strand breaks and, more generally, DNA damage on its genome. In fact, the DNA repair genes systematically avoid regions with a high recombination density, as they need to be intact to accurately handle repairable DNA lesions.

Genome-wide mapping and profiling of γH2AX binding hotspots in response to different replication stress inducers

BACKGROUND

Replication stress (RS) gives rise to DNA damage that threatens genome stability. RS can originate from different sources that stall replication by diverse mechanisms. However, the mechanism underlying how different types of RS contribute to genome instability is unclear, in part due to the poor understanding of the distribution and characteristics of damage sites induced by different RS mechanisms.

RESULTS

We use ChIP-seq to map γH2AX binding sites genome-wide caused by aphidicolin (APH), hydroxyurea (HU), and methyl methanesulfonate (MMS) treatments in human lymphocyte cells. Mapping of γH2AX ChIP-seq reveals that APH, HU, and MMS treatments induce non-random γH2AX chromatin binding at discrete regions, suggesting that there are γH2AX binding hotspots in the genome. Characterization of the distribution and sequence/epigenetic features of γH2AX binding sites reveals that the three treatments induce γH2AX binding at largely non-overlapping regions, suggesting that RS may cause damage at specific genomic loci in a manner dependent on the fork stalling mechanism. Nonetheless, γH2AX binding sites induced by the three treatments share common features including compact chromatin, coinciding with larger-than-average genes, and depletion of CpG islands and transcription start sites. Moreover, we observe significant enrichment of SINEs in γH2AX sites in all treatments, indicating that SINEs may be a common barrier for replication polymerases.

CONCLUSIONS

Our results identify the location and common features of genome instability hotspots induced by different types of RS, and help in deciphering the mechanisms underlying RS-induced genetic diseases and carcinogenesis.

Hemizygosity Enables a Mutational Transition Governing Fungal Virulence and Commensalism

Candida albicans is a commensal fungus of human gastrointestinal and reproductive tracts, but also causes life-threatening systemic infections. The balance between colonization and pathogenesis is associated with phenotypic plasticity, with alternative cell states producing different outcomes in a mammalian host. Here, we reveal that gene dosage of a master transcription factor regulates cell differentiation in diploid C. albicans cells, as EFG1 hemizygous cells undergo a phenotypic transition inaccessible to "wild-type" cells with two functional EFG1 alleles. Notably, clinical isolates are often EFG1 hemizygous and thus licensed to undergo this transition. Phenotypic change corresponds to high-frequency loss of the functional EFG1 allele via de novo mutation or gene conversion events. This phenomenon also occurs during passaging in the gastrointestinal tract with the resulting cell type being hypercompetitive for commensal and systemic infections. A "two-hit" genetic model therefore underlies a key phenotypic transition in C. albicans that enables adaptation to host niches.

Modeling double strand break susceptibility to interrogate structural variation in cancer

BACKGROUND

Structural variants (SVs) are known to play important roles in a variety of cancers, but their origins and functional consequences are still poorly understood. Many SVs are thought to emerge from errors in the repair processes following DNA double strand breaks (DSBs).

RESULTS

We used experimentally quantified DSB frequencies in cell lines with matched chromatin and sequence features to derive the first quantitative genome-wide models of DSB susceptibility. These models are accurate and provide novel insights into the mutational mechanisms generating DSBs. Models trained in one cell type can be successfully applied to others, but a substantial proportion of DSBs appear to reflect cell type-specific processes. Using model predictions as a proxy for susceptibility to DSBs in tumors, many SV-enriched regions appear to be poorly explained by selectively neutral mutational bias alone. A substantial number of these regions show unexpectedly high SV breakpoint frequencies given their predicted susceptibility to mutation and are therefore credible targets of positive selection in tumors. These putatively positively selected SV hotspots are enriched for genes previously shown to be oncogenic. In contrast, several hundred regions across the genome show unexpectedly low levels of SVs, given their relatively high susceptibility to mutation. These novel coldspot regions appear to be subject to purifying selection in tumors and are enriched for active promoters and enhancers.

CONCLUSIONS

We conclude that models of DSB susceptibility offer a rigorous approach to the inference of SVs putatively subject to selection in tumors.

Genome-wide profiling of Epstein-Barr virus integration by targeted sequencing in Epstein-Barr virus associated malignancies

Rationale: Epstein-Barr virus (EBV) is associated with multiple malignancies with expression of viral oncogenic proteins and chronic inflammation as major mechanisms contributing to tumor development. A less well-studied mechanism is the integration of EBV into the human genome possibly at sites which may disrupt gene expression or genome stability.

Methods: We sequenced tumor DNA to profile the EBV sequences by hybridization-based enrichment. Bioinformatic analysis was used to detect the breakpoints of EBV integrations in the genome of cancer cells.

Results: We identified 197 breakpoints in nasopharyngeal carcinomas and other EBV-associated malignancies. EBV integrations were enriched at vulnerable regions of the human genome and were close to tumor suppressor and inflammation-related genes. We found that EBV integrations into the introns could decrease the expression of the inflammation-related genes, TNFAIP3, PARK2, and CDK15, in NPC tumors. In the EBV genome, the breakpoints were frequently at oriP or terminal repeats. These breakpoints were surrounded by microhomology sequences, consistent with a mechanism for integration involving viral genome replication and microhomology-mediated recombination.

Conclusion: Our finding provides insight into the potential of EBV integration as an additional mechanism mediating tumorigenesis in EBV associated malignancies.

The role of fork stalling and DNA structures in causing chromosome fragility

Alternative non-B form DNA structures, also called secondary structures, can form in certain DNA sequences under conditions that produce single-stranded DNA, such as during replication, transcription, and repair. Direct links between secondary structure formation, replication fork stalling, and genomic instability have been found for many repeated DNA sequences that cause disease when they expand. Common fragile sites (CFSs) are known to be AT-rich and break under replication stress, yet the molecular basis for their fragility is still being investigated. Over the past several years, new evidence has linked both the formation of secondary structures and transcription to fork stalling and fragility of CFSs. How these two events may synergize to cause fragility and the role of nuclease cleavage at secondary structures in rare and CFSs are discussed here. We also highlight evidence for a new hypothesis that secondary structures at CFSs not only initiate fragility but also inhibit healing, resulting in their characteristic appearance.

Endogenous DNA Double-Strand Breaks during DNA Transactions: Emerging Insights and Methods for Genome-Wide Profiling

DNA double-strand breaks (DSBs) jeopardize genome integrity and can-when repaired unfaithfully-give rise to structural rearrangements associated with cancer. Exogenous agents such as ionizing radiation or chemotherapy can invoke DSBs, but a vast amount of breakage arises during vital endogenous DNA transactions, such as replication and transcription. Additionally, chromatin looping involved in 3D genome organization and gene regulation is increasingly recognized as a possible contributor to DSB events. In this review, we first discuss insights into the mechanisms of endogenous DSB formation, showcasing the trade-off between essential DNA transactions and the intrinsic challenges that these processes impose on genomic integrity. In the second part, we highlight emerging methods for genome-wide profiling of DSBs, and discuss future directions of research that will help advance our understanding of genome-wide DSB formation and repair.

The Protective Role of Dormant Origins in Response to Replicative Stress

Genome stability requires tight regulation of DNA replication to ensure that the entire genome of the cell is duplicated once and only once per cell cycle. In mammalian cells, origin activation is controlled in space and time by a cell-specific and robust program called replication timing. About 100,000 potential replication origins form on the chromatin in the gap 1 (G1) phase but only 20⁻30% of them are active during the DNA replication of a given cell in the synthesis (S) phase. When the progress of replication forks is slowed by exogenous or endogenous impediments, the cell must activate some of the inactive or "dormant" origins to complete replication on time. Thus, the many origins that may be activated are probably key to protect the genome against replication stress. This review aims to discuss the role of these dormant origins as safeguards of the human genome during replicative stress.