Search

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1–3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4–7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (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 papillomavirus9 (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.

Supplementary Figures 1-11, Supplementary Tables 1-2

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

BackgroundCancer somatic mutations are the product of multiple mutational and repair processes, which are tightly associated with DNA replication. Distinctive patterns of somatic mutations accumulation in tumors, termed mutational signatures, are indicative of processes the tumors underwent. While tumor mutational load is correlated with late replicating regions and spatial genome organization, much is unknown about the association of many different mutational processes and replication timing, and the interplay with chromatin structure remains an open question.MethodsWe systematically analyzed the mutational landscape of 2,787 WGS tumors from 32 different tumor types separately for early and late replicating regions. We used sequence context normalization and chromatin data to account for sequence and chromatin accessibility differences between early and late replicating regions. Moreover, we expanded the signature analyses to doublet base substitutions and small insertions and deletions by developing an artificial genomes-based approach to account for sequence differences between various genomic regions.ResultsWe revealed the replication timing (RT) association of single base, doublet base and small insertions and deletions mutational signatures. The association is signature specific: some are associated with early or late replication (such as UV-exposure signatures SBS7b and SBS7a, respectively) and others have no association. Most associations exist even after normalizing for genome accessibility. We further developed a focused mutational signature identification approach, which uses RT information to improve signature identification, and found that SBS16, which is biased towards early replication, is strongly associated with better survival rates in liver cancer.ConclusionsOur comprehensive analyses enabled a more robust classification of RT association of single base, doublet base and indels signatures. By doing so, we demonstrated a variation in the association with RT, as many mutational processes biased towards either early or late replication timing, and others have an equal RT distribution. These associations were independent from chromatin accessibility in most cases. This work highlights that restricting signatures analyses to concise genomic regions improves identification of signatures, such as SBS16, and demonstrates its clinically relevance as a predictor of improved survival of liver cancer patients.

Funder: National Institute of Health (NIH) DP2-GM123495, DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome���s replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) ��� sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Mutations do not occur uniformly across genomes but instead show biased associations with various genomic features, most notably late replication timing. However, it remains contested which mutation types in human cells relate to DNA replication dynamics and to what extents. Previous studies have been limited by the absence of cell-type-specific replication timing profiles and lack of consideration of inter-individual variation. To overcome these limitations, we performed high-resolution comparisons of mutational landscapes between and within lymphoblastoid cell lines from 1662 individuals, 151 chronic lymphocytic leukemia patients, and three colon adenocarcinoma cell lines including two with mismatch repair deficiency. Using cell type-matched replication timing profiles, we demonstrate how mutational pathways can exhibit heterogeneous replication timing associations. We further identified global mutation load as a novel, pervasive determinant of mutational landscape heterogeneity across individuals. Specifically, elevated mutation load corresponded to increased late replication timing bias as well as replicative strand asymmetries of clock-like mutations and off-target somatic hypermutation. The association of somatic hypermutation with DNA replication timing was further influenced by mutational clustering. Considering these multivariate factors, and by incorporating mutation phasing at an unprecedented scale, we identified a unique mutational landscape on the inactive X-chromosome. Overall, we report underappreciated complexity of mutational pathways and their relationship to replication timing and identify specific factors underlying differential mutation landscapes among cell types and individuals.

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species’ phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited inter-individual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and impacts regulatory evolution.

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.

Alteration of normal ploidy (aneuploidy) can have a number of opposing effects, such as unbalancing protein abundances and inhibiting cell growth but also accelerating genetic diversification and rapid adaptation. The interplay of these detrimental and beneficial effects remains puzzling. Here, to understand how cells develop tolerance to aneuploidy, we subject disomic (i.e. with an extra chromosome copy) strains of yeast to long-term experimental evolution under strong selection, by forcing disomy maintenance and daily population dilution. We characterize mutations, karyotype alterations and gene expression changes, and dissect the associated molecular strategies. Cells with different extra chromosomes accumulated mutations at distinct rates and displayed diverse adaptive events. They tended to evolve towards normal ploidy through chromosomal DNA loss and gene expression changes. We identify genes with recurrent mutations and altered expression in multiple lines, revealing a variant that improves growth under genotoxic stresses. These findings support rapid evolvability of disomic strains that can be used to characterize fitness effects of mutations under different stress conditions., Aneuploidy (abnormal chromosome number) can enable rapid adaptation to stress conditions, but it also entails fitness costs from gene imbalance. Here, the authors experimentally evolve yeast while forcing maintenance of aneuploidy to identify the mechanisms that promote tolerance of aneuploidy.

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

Genome instability (GIN) is a main contributing factor to congenital and somatic diseases, but its sporadic occurrence in individual cell cycles makes it difficult to study mechanistically. One profound manifestation of GIN is the formation of micronuclei (MN), the engulfment of chromosomes or chromosome fragments in their own nuclear structures separate from the main nucleus. Here, we developed MN-seq, an approach for sequencing the DNA contained within micronuclei. We applied MN-seq to mice with mutations in Mcm4 and Rad9a, which disrupt DNA replication, repair, and damage responses. Data analysis and simulations show that centromere presence, fragment length, and a heterogenous landscape of chromosomal fragility all contribute to the patterns of DNA present within MN. In particular, we show that long genes, but also gene-poor regions, are associated with chromosome breaks that lead to the enrichment of particular genomic sequences in MN, in a genetic background-specific manner. Finally, we introduce single-cell micronucleus sequencing (scMN-seq), an approach to sequence the DNA present in MN of individual cells. Together, sequencing micronuclei provides a systematic approach for studying GIN and reveals novel molecular associations with chromosome breakage and segregation.

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. MCM10 mutant cells demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in MCM10-mutant cells was predominantly comprised of replication initiation defects. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel modulator of DNA replication timing.

Regulation of DNA replication and copy number is necessary to promote genome stability and maintain cell and tissue function. DNA replication is regulated temporally in a process known as replication timing (RT). Rap1-interacting factor 1 (Rif1) is a key regulator of RT and has a critical function in copy number control in polyploid cells. Previously, we demonstrated that Rif1 functions with SUUR to inhibit replication fork progression and promote underreplication (UR) of specific genomic regions. How Rif1-dependent control of RT factors into its ability to promote UR is unknown. By applying a computational approach to measure RT in Drosophila polyploid cells, we show that SUUR and Rif1 have differential roles in controlling UR and RT. Our findings reveal that Rif1 acts to promote late replication, which is necessary for SUUR-dependent underreplication. Our work provides new insight into the process of UR and its links to RT., In polyploid cells, copy number is not uniform throughout the genome. Underreplication causes pericentric heterochromatin and defined regions of euchromatin to have reduced copy number relative to overall ploidy. SUUR and Rif1 are key regulators of underreplication, and Rif1 is a known regulator of replication timing (RT). Das et al. take a computational approach to measure replication timing in polyploid cells. Their results demonstrate that, while Rif1 and SUUR both promote underreplication, they differentially affect RT.

Fanconi anemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink (ICL) repair resulting in chromosome breakage1–3. The FA repair pathway comprises at least 22 FANC proteins including BRCA1 and BRCA24–6, and protects against carcinogenic endogenous and exogenous aldehydes7–10. Individuals with FA are hundreds to thousands-fold more likely to develop head and neck (HNSCC), esophageal and anogenital squamous cell carcinomas (SCCs) with a median onset age of 31 years11. The aggressive nature of these tumors and poor patient tolerance of platinum and radiation-based therapy have been associated with short survival in FA11–16. Molecular studies of SCCs from individuals with FA (FA SCCs) have been limited, and it is unclear how they relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or human papillomavirus (HPV) infection17. Here, by sequencing FA SCCs, we demonstrate that the primary genomic signature of FA-deficiency is the presence of a high number of structural variants (SVs). SVs are enriched for small deletions, unbalanced translocations, and fold-back inversions that arise in the context of TP53 loss. The SV breakpoints preferentially localize to early replicating regions, common fragile sites, tandem repeats, and SINE elements. SVs are often connected forming complex rearrangements. Resultant genomic instability underlies elevated copy number alteration (CNA) rates of key HNSCC-associated genes, including PIK3CA, MYC, CSMD1, PTPRD, YAP1, MXD4, and EGFR. In contrast to sporadic HNSCC, we find no evidence of HPV infection in FA HNSCC, although positive cases were identified in gynecologic tumors. A murine allograft model of FA pathway-deficient SCC was enriched in SVs, exhibited dramatic tumor growth advantage, more rapid epithelial-to-mesenchymal transition (EMT), and enhanced autonomous inflammatory signaling when compared to an FA pathway-proficient model. In light of the protective role of the FA pathway against SV formation uncovered here, and recent findings of FA pathway insufficiency in the setting of increased formaldehyde load resulting in hematopoietic stem cell failure and carcinogenesis18–20, we propose that high copy-number instability in sporadic HNSCC may result from functional overload of the FA pathway by endogenous and exogenous DNA crosslinking agents. Our work lays the foundation for improved FA patient treatment and demonstrates that FA SCC is a powerful model to study tumorigenesis resulting from DNA crosslinking damage.

Haploid human embryonic stem cells (ESCs) provide a powerful genetic system but diploidize at high rates. We hypothesized that diploidization results from aberrant DNA replication. To test this, we profiled DNA replication timing in isogenic haploid and diploid ESCs. The greatest difference was the earlier replication of the X Chromosome in haploids, consistent with the lack of X-Chromosome inactivation. We also identified 21 autosomal regions that had delayed replication in haploids, extending beyond the normal S phase and into G2/M. Haploid-delays comprised a unique set of quiescent genomic regions that are also underreplicated in polyploid placental cells. The same delays were observed in female ESCs with two active X Chromosomes, suggesting that increased X-Chromosome dosage may cause delayed autosomal replication. We propose that incomplete replication at the onset of mitosis could prevent cell division and result in re-entry into the cell cycle and whole genome duplication.

The composition of a plant virome may represent spatiotemporal patterns of plant virus abundance. Using next generation sequencing, we investigated the viromes of watermelon fruits grown in two adjacent open fields located in Eastern Israel: Kalia and Mitzpe-Shalem. The two viromes were comprised of distinct virus species and genera. Studying spatial and temporal effects on virus occurrence, we detected the crinivirus Cucurbit yellow stunting disorder virus in Kalia and not in Mitzpe-Shalem, irrespective of collection time. A spatial effect was also observed regarding the occurrence of the begomovirus Watermelon chlorotic stunt virus, which was not detected in watermelons from Kalia, but rather in watermelon fruit in Northern Israel in 2018. A temporal effect was observed on the appearance of the ipomovirus Cucumber vein yellowing virus, which was detected in watermelons from Kalia in 2017 and was absent in the 2016 virome analysis. Importantly, regardless of temporal and spatial effects, the crinivirus Cucurbit chlorotic yellows virus, which was new to Israeli landscape, was detected in all the tested watermelon plants. These changes in viral abundance were associated with an early whitefly (Bemisia tabaci) occurrence, which might be the cause for the severe disease spread in watermelons.

Fanconi anemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink (ICL) repair resulting in chromosome breakage. The FA repair pathway comprises at least 22 FANC proteins including BRCA1 and BRCA2 and protects against carcinogenic endogenous and exogenous aldehydes. Individuals with FA are hundreds to thousands-fold more likely to develop head and neck (HNSCC), esophageal and anogenital squamous cell carcinomas (SCCs) with a median onset age of 31 years. The aggressive nature of these tumors and poor patient tolerance of platinum and radiation-based therapy have been associated with short survival in FA. Molecular studies of SCCs from individuals with FA (FA SCCs) have been limited, and it is unclear how they relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or human papillomavirus (HPV) infection. Here, by sequencing FA SCCs, we demonstrate that the primary genomic signature of FA-deficiency is the presence of a high number of structural variants (SVs). SVs are enriched for small deletions, unbalanced translocations, and fold-back inversions that arise in the context of TP53 loss. The SV breakpoints preferentially localize to early replicating regions, common fragile sites, tandem repeats, and SINE elements. SVs are often connected forming complex rearrangements. Resultant genomic instability underlies elevated copy number alteration (CNA) rates of key HNSCC-associated genes, including PIK3CA, MYC, CSMD1, PTPRD, YAP1, MXD4, and EGFR. In contrast to sporadic HNSCC, we find no evidence of HPV infection in FA HNSCC, although positive cases were identified in gynecologic tumors. A murine allograft model of FA pathway-deficient SCC was enriched in SVs, exhibited dramatic tumor growth advantage, more rapid epithelial-to-mesenchymal transition, and enhanced autonomous inflammatory signaling when compared to an FA pathway-proficient model. In light of the protective role of the FA pathway against SV formation uncovered here, and recent findings of FA pathway insufficiency in the setting of increased formaldehyde load resulting in hematopoietic stem cell failure and carcinogenesis, we propose that high copy-number instability in sporadic HNSCC may result from functional overload of the FA pathway by endogenous and exogenous DNA crosslinking agents. Our work lays the foundation for improved FA patient treatment and demonstrates that FA SCC is a powerful model to study tumorigenesis resulting from DNA crosslinking damage. Citation Format: Andrew L. Webster, Mathijs A. Sanders, Krupa Patel, Ralf Dietrich, Raymond J. Noonan, Francis P. Lach, Ryan R. White, Audrey M. Goldfarb, Kevin Hadi, Matthew M. Edwards, Frank X. Donovan, Moonjung Jung, Sunandini Sridhar, Olivier Fedrigo, Huasong Tian, Joel Rosiene, Thomas Heineman, Jennifer Kennedy, Lorenzo Bean, Rasim O. Rosti, Rebecca Tryon, Ashlyn-Maree Gonzalez, Allana Rosenberg, Ji-Dung Luo, Thomas Carrol, Eunike Velleuer, Jeff C. Rastatter, Susanne I. Wells, Jordi Surrallés, Grover Bagby, Margaret L. MacMillan, John E. Wagner, Maria Cancio, Farid Boulad, Theresa Scognamiglio, Roger Vaughan, Amnon Koren, Marcin Imielinski, Settara Chandrasekharappa, Arleen D. Auerbach, Bhuvanesh Singh, David Kutler, Peter J. Campbell, Agata Smogorzewska. Fanconi anemia pathway deficiency drives copy number variation in squamous cell carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6196.

DNA replication initiates from replication origins firing at different times throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We developed an approach to profile DNA replication from whole-genome sequencing of thousands of single cells. We describe "in silico flow cytometry", a method for discriminating replicating cells with superior accuracy to FACS and staging them across S phase. Using two microfluidic platforms, we analyzed up to 2,428 individual replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that the vast majority are in confined genomic regions. While initiation occurs in a remarkably similar order across cells, we unexpectedly identified a subset of initiation regions that constitutively fire in late S phase, and another subset firing randomly throughout S phase. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression.

Significance Ataxia with oculomotor apraxia (AOA) is a progressive neurodegenerative disease characterized by early‐onset autosomal recessive cerebellar ataxia with oculomotor apraxia, peripheral axonal neuropathy, and impaired motor functions. The AOA-2 subgroup results from mutations in an RNA/DNA helicase, Senataxin, which is encoded by the SETX gene. Here, we carried out integrated genome and transcriptome analyses of cell lines derived from individuals with AOA2, as well as CRISPR/Cas9 generated SETX knockouts, and observed genome-wide chromosome fragility. Genome instability was caused by increased transcription stress and the accumulation of RNA/DNA hybrids near gene promotors, resulting in aberrant DNA repair that led to changes in gene-expression profiles. The results indicate that SETX-defective cells exhibit transcription stress that leads to chromosome fragility., 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.

Human cleavage stage embryos frequently acquire chromosomal aneuploidies during mitosis. We recently showed that double strand breaks (DSBs) introduced by Cas9 result in frequent chromosomal alterations. Here we show that spontaneous DSBs arise de novo during embryonic DNA replication and show as chromosomal breaks that occur primarily, in centromeric, and gene-poor areas, with the latter harboring long genes implicated in neurodevelopmental disorders, including CNTN5, IMMP2L, LRP1B, RYR2, CSMD1. While transcription-replication conflicts are not required for chromosome fragility, fragile, gene-poor regions contain long-traveling replication forks, and complete replication late in G2 phase, hours after most of the genome has been duplicated. Incomplete replication at these sites in both human and mouse zygotes, results in DSBs, whole and segmental chromosome errors, micronucleation, chromosome fragmentation and poor embryo quality. Thus, failure to complete DNA replication is a mechanism of developmental failure, and a source of detrimental genetic alterations in humans.

Motivation Genomic DNA replicates according to a reproducible spatiotemporal program, with some loci replicating early in S phase while others replicate late. Despite being a central cellular process, DNA replication timing studies have been limited in scale due to technical challenges. Results We present TIGER (Timing Inferred from Genome Replication), a computational approach for extracting DNA replication timing information from whole genome sequence data obtained from proliferating cell samples. The presence of replicating cells in a biological specimen leads to non-uniform representation of genomic DNA that depends on the timing of replication of different genomic loci. Replication dynamics can hence be observed in genome sequence data by analyzing DNA copy number along chromosomes while accounting for other sources of sequence coverage variation. TIGER is applicable to any species with a contiguous genome assembly and rivals the quality of experimental measurements of DNA replication timing. It provides a straightforward approach for measuring replication timing and can readily be applied at scale. Availability and implementation TIGER is available at https://github.com/TheKorenLab/TIGER. Supplementary information Supplementary data are available at Bioinformatics online.

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure and gene regulation. However, the genetic basis of the mammalian DNA replication timing program is poorly understood1–3. To systematically identify genetic regulators of DNA replication timing, we exploited inter-individual variation in 457 human pluripotent stem cell lines from 349 individuals. We show that the human genome’s replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs4) – base-pair-resolution sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, to affect replication timing. Analysis of rtQTL locations reveals a histone code for replication initiation, composed of bivalent histone H3 trimethylation marks on a background of histone hyperacetylation. The H3 trimethylation marks are individually repressive yet synergize to promote early replication. We further identify novel positive and negative regulators of DNA replication timing, the former comprised of pluripotency-related transcription factors while the latter involve boundary elements. Human replication timing is controlled by a multi-layered mechanism that operates on target DNA sequences, is composed of dozens of effectors working combinatorially, and follows principles analogous to transcription regulation: a histone code, activators and repressors, and a promoter-enhancer logic.

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon (POLE) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations.

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

During October 2014, unfamiliar mild mosaic and mottling symptoms were identified on leaves of pepper (Capsicum chinense cv. Habanero) seedlings grown in the Arava valley in Israel 2–3 weeks post planting. Symptomatic plants were tested positive by ELISA using laboratory-produced antisera for tobamovirus species. Typical tobamovirus rod-shaped morphology was observed by transmission electron microscopy (TEM) analysis of purified virion preparation that was used for mechanical inoculation of laboratory test plants for the completion of Koch’s postulates. The complete viral genome was sequenced from small interfering RNA purified from symptomatic pepper leaves and fruits by next-generation sequencing (NGS) using Illumina MiSeq platform. The contigs generated by the assembly covered 80% of the viral genome. RT-PCR amplification and Sanger sequencing were employed in order to validate the sequence generated by NGS technology. The nucleotide sequence of the complete viral genome was 99% identical to the complete genome of Paprika mild mottle virus isolate from Japan (PaMMV-J), and the deduced amino acid sequence was 99% identical to PaMMV-J protein. Amplicons from seed RNA showed 100% identity to the viral isolate from the collected symptomatic pepper plants. Partial host range analysis revealed a slow development of systemic infection in inoculated tomato plants (Lycopersicon esculentum). Interestingly, double inoculation of susceptible wild-type tomato plants and Tm-22-resistant tomato plants with the PaMMV-IL and Tomato brown rugose fruit virus (ToBRFV) resulted in accelerated viral expression in the plants.

Microsatellites (MSs) are tracts of variable-length repeats of short DNA motifs that exhibit high rates of mutation in the form of insertions or deletions (indels) of the repeated motif. Despite their prevalence, the contribution of somatic MS indels to cancer has been largely unexplored, owing to difficulties in detecting them in short-read sequencing data. Here we present two tools: MSMuTect, for accurate detection of somatic MS indels, and MSMutSig, for identification of genes containing MS indels at a higher frequency than expected by chance. Applying MSMuTect to whole-exome data from 6, 747 human tumors representing 20 tumor types, we identified >1, 000 previously undescribed MS indels in cancer genes. Additionally, we demonstrate that the number and pattern of MS indels can accurately distinguish microsatellite-stable tumors from tumors with microsatellite instability, thus potentially improving classification of clinically relevant subgroups. Finally, we identified seven MS indel driver hotspots: four in known cancer genes (ACVR2A, RNF43, JAK1, and MSH3) and three in genes not previously implicated as cancer drivers (ESRP1, PRDM2, and DOCK3).

Sensitive detection of somatic copy number alterations (SCNA) in cancer genomes is confounded by “waviness” in read depth data. We present dryclean, a signal processing algorithm to optimize SCNA detection in whole genome (WGS) and targeted sequencing platforms through foreground detection and background subtraction of read depth data. Application of dryclean to WGS demonstrates that WGS waviness is driven by replication timing. Re-analysis of thousands of tumor profiles reveals that dryclean provides superior detection of biologically relevant SCNAs relative to state-of-the-art algorithms. Applied to in silico tumor dilutions, dryclean improves the sensitivity of relapse detection 10-fold relative to current standards. dryclean is available as an R package in the GitHub repository https://github.com/mskilab/dryclean

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analogue incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

Some diseases are caused by coinfection of several pathogens in the same plant. However, studies on the complexity of these coinfection events under different environmental conditions are scarce. Our ongoing research involves late wilting disease of cucumber caused by coinfection of Cucumber green mottle mosaic virus (CGMMV) and Pythium spp. We specifically investigated the role of various temperatures (18, 25, 32°C) on the coinfection by CGMMV and two predominant Pythium species occurring in cucumber greenhouses under Middle Eastern climatic conditions. During the summer months, Pythium aphanidermatum was most common, whereas P. spinosum predominated during the winter–spring period. P. aphanidermatum preferred higher temperatures while P. spinosum preferred low temperatures and caused very low levels of disease at 32°C when the 6-day-old seedlings were infected with P. spinosum alone. Nevertheless, after applying a later coinfection with CGMMV on the 14-day-old plants, a synergistic effect was detected for both Pythium species at optimal and suboptimal temperatures, with P. spinosum causing high mortality incidence even at 32°C. The symptoms caused by CGMMV infection appeared earlier as the temperature increased. However, within each temperature, no significant influence of the combined infection was detected. Our results demonstrate the complexity of coinfection in changing environmental conditions and indicate its involvement in disease development and severity as compared with infection by each of the pathogens alone.

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.