Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions

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.

Transcription elongation defects link oncogenic SF3B1 mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human diseases remains unexplored. Using isogenic cell lines, patient samples, and a mutant mouse model, we investigated how cancer-associated mutations in SF3B1 alter transcription. We found that these mutations reduce the elongation rate of RNA polymerase II (RNAPII) along gene bodies and its density at promoters. The elongation defect results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1. The decreased promoter-proximal RNAPII density reduces both chromatin accessibility and H3K4me3 marks at promoters. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC/H3K4me pathway, which, when modulated, reverse both transcription and chromatin changes. Our findings reveal how splicing factor mutant states behave functionally as epigenetic disorders through impaired transcription-related changes to the chromatin landscape. We also present a rationale for targeting the Sin3/HDAC complex as a therapeutic strategy.

Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

Fork coupling directs DNA replication elongation and termination

DNA replication is initiated at multiple loci to ensure timely duplication of eukaryotic genomes. Sister replication forks progress bidirectionally, and replication terminates when two convergent forks encounter one another. To investigate the coordination of replication forks, we developed a replication-associated in situ HiC method to capture chromatin interactions involving nascent DNA. We identify more than 2000 fountain-like structures of chromatin contacts in human and mouse genomes, indicative of coupling of DNA replication forks. Replication fork interaction not only occurs between sister forks but also involves forks from two distinct origins to predetermine replication termination. Termination-associated chromatin fountains are sensitive to replication stress and lead to coupled forks-associated genomic deletions in cancers. These findings reveal the spatial organization of DNA replication forks within the chromatin context.

Integrative analysis of DNA replication origins and ORC/MCM binding sites in human cells reveals a lack of overlap

Based on experimentally determined average inter-origin distances of ∼100 kb, DNA replication initiates from ∼50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the Origin Recognition Complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and 5 ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ∼7.5 million union origins identified by all datasets, only 0.27% were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques (20,250 shared origins), suggesting extensive variability in origin usage and identification. 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF binding sites, G-quadruplex sites and activating histone marks, these overlaps are comparable or less than that of known Transcription Start Sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ∼13,000 reproducible ORC binding sites in human cancer cells, and only 4.5% were within 1 kb of the ∼11,000 union MCM2-7 binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, S. cerevisiae . Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

The chromatin network helps prevent cancer-associated mutagenesis at transcription-replication conflicts

Genome instability is a feature of cancer cells, transcription being an important source of DNA damage. This is in large part associated with R-loops, which hamper replication, especially at head-on transcription-replication conflicts (TRCs). Here we show that TRCs trigger a DNA Damage Response (DDR) involving the chromatin network to prevent genome instability. Depletion of the key chromatin factors INO80, SMARCA5 and MTA2 results in TRCs, fork stalling and R-loop-mediated DNA damage which mostly accumulates at S/G2, while histone H3 Ser10 phosphorylation, a mark of chromatin compaction, is enriched at TRCs. Strikingly, TRC regions show increased mutagenesis in cancer cells with signatures of homologous recombination deficiency, transcription-coupled nucleotide excision repair (TC-NER) and of the AID/APOBEC cytidine deaminases, being predominant at head-on collisions. Thus, our results support that the chromatin network prevents R-loops and TRCs from genomic instability and mutagenic signatures frequently associated with cancer.

Cohesin maintains replication timing to suppress DNA damage on cancer genes

Cohesin loss-of-function mutations are frequently observed in tumors, but the mechanism underlying its role in tumorigenesis is unclear. Here, we found that depletion of RAD21, a core subunit of cohesin, leads to massive genome-wide DNA breaks and 147 translocation hotspot genes, co-mutated with cohesin in multiple cancers. Increased DNA damages are independent of RAD21-loss-induced transcription alteration and loop anchor disruption. However, damage-induced chromosomal translocations coincide with the asymmetrically distributed Okazaki fragments of DNA replication, suggesting that RAD21 depletion causes replication stresses evidenced by the slower replication speed and increased stalled forks. Mechanistically, approximately 30% of the human genome exhibits an earlier replication timing after RAD21 depletion, caused by the early initiation of >900 extra dormant origins. Correspondingly, most translocation hotspot genes lie in timing-altered regions. Therefore, we conclude that cohesin dysfunction causes replication stresses induced by excessive DNA replication initiation, resulting in gross DNA damages that may promote tumorigenesis.

OKseqHMM: a genome-wide replication fork directionality analysis toolkit

During each cell division, tens of thousands of DNA replication origins are co-ordinately activated to ensure the complete duplication of the human genome. However, replication fork progression can be challenged by many factors, including co-directional and head-on transcription-replication conflicts (TRC). Head-on TRCs are more dangerous for genome integrity. To study the direction of replication fork movement and TRCs, we developed a bioinformatics toolkit called OKseqHMM (https://github.com/CL-CHEN-Lab/OK-Seq, https://doi.org/10.5281/zenodo.7428883). Then, we used OKseqHMM to analyse a large number of datasets obtained by Okazaki fragment sequencing to directly measure the genome-wide replication fork directionality (RFD) and to accurately predict replication initiation and termination at a fine resolution in organisms including yeast, mouse and human. We also successfully applied our analysis to other genome-wide sequencing techniques that also contain RFD information (e.g. eSPAN, TrAEL-seq). Our toolkit can be used to predict replication initiation and fork progression direction genome-wide in a wide range of cell models and growth conditions. Comparing the replication and transcription directions allows identifying loci at risk of TRCs, particularly head-on TRCs, and investigating their role in genome instability by checking DNA damage data, which is of prime importance for human health.

Transcription elongation defects link oncogenic splicing factor mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human disease remains unexplored. Here, we investigated the impact of non-synonymous mutations in SF3B1 and U2AF1, two commonly mutated splicing factors in cancer, on transcription. We find that the mutations impair RNA Polymerase II (RNAPII) transcription elongation along gene bodies leading to transcription-replication conflicts, replication stress and altered chromatin organization. This elongation defect is linked to disrupted pre-spliceosome assembly due to impaired association of HTATSF1 with mutant SF3B1. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC complex, which, when modulated, normalize transcription defects and their downstream effects. Our findings shed light on the mechanisms by which oncogenic mutant spliceosomes impact chromatin organization through their effects on RNAPII transcription elongation and present a rationale for targeting the Sin3/HDAC complex as a potential therapeutic strategy.

GRAPHICAL ABSTRACT

HIGHLIGHTS

Oncogenic mutations of SF3B1 and U2AF1 cause a gene-body RNAPII elongation defectRNAPII transcription elongation defect leads to transcription replication conflicts, DNA damage response, and changes to chromatin organization and H3K4me3 marksThe transcription elongation defect is linked to disruption of the early spliceosome formation through impaired interaction of HTATSF1 with mutant SF3B1.Changes to chromatin organization reveal potential therapeutic strategies by targeting the Sin3/HDAC pathway.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

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

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

Effects of replication domains on genome-wide UV-induced DNA damage and repair

Nucleotide excision repair is the primary repair mechanism that removes UV-induced DNA lesions in placentals. Unrepaired UV-induced lesions could result in mutations during DNA replication. Although the mutagenesis of pyrimidine dimers is reasonably well understood, the direct effects of replication fork progression on nucleotide excision repair are yet to be clarified. Here, we applied Damage-seq and XR-seq techniques and generated replication maps in synchronized UV-treated HeLa cells. The results suggest that ongoing replication stimulates local repair in both early and late replication domains. Additionally, it was revealed that lesions on lagging strand templates are repaired slower in late replication domains, which is probably due to the imbalanced sequence context. Asymmetric relative repair is in line with the strand bias of melanoma mutations, suggesting a role of exogenous damage, repair, and replication in mutational strand asymmetry.

Determination of human DNA replication origin position and efficiency reveals principles of initiation zone organisation

Replication of the human genome initiates within broad zones of ∼150 kb. The extent to which firing of individual DNA replication origins within initiation zones is spatially stochastic or localised at defined sites remains a matter of debate. A thorough characterisation of the dynamic activation of origins within initiation zones is hampered by the lack of a high-resolution map of both their position and efficiency. To address this shortcoming, we describe a modification of initiation site sequencing (ini-seq), based on density substitution. Newly replicated DNA is rendered 'heavy-light' (HL) by incorporation of BrdUTP while unreplicated DNA remains 'light-light' (LL). Replicated HL-DNA is separated from unreplicated LL-DNA by equilibrium density gradient centrifugation, then both fractions are subjected to massive parallel sequencing. This allows precise mapping of 23,905 replication origins simultaneously with an assignment of a replication initiation efficiency score to each. We show that origin firing within early initiation zones is not randomly distributed. Rather, origins are arranged hierarchically with a set of very highly efficient origins marking zone boundaries. We propose that these origins explain much of the early firing activity arising within initiation zones, helping to unify the concept of replication initiation zones with the identification of discrete replication origin sites.

The SWI/SNF chromatin remodeling complex helps resolve R-loop-mediated transcription-replication conflicts

ATP-dependent chromatin remodelers are commonly mutated in human cancer. Mammalian SWI/SNF complexes comprise three conserved multisubunit chromatin remodelers (cBAF, ncBAF and PBAF) that share the BRG1 (also known as SMARCA4) subunit responsible for the main ATPase activity. BRG1 is the most frequently mutated Snf2-like ATPase in cancer. In the present study, we have investigated the role of SWI/SNF in genome instability, a hallmark of cancer cells, given its role in transcription, DNA replication and DNA-damage repair. We show that depletion of BRG1 increases R-loops and R-loop-dependent DNA breaks, as well as transcription-replication (T-R) conflicts. BRG1 colocalizes with R-loops and replication fork blocks, as determined by FANCD2 foci, with BRG1 depletion being epistatic to FANCD2 silencing. Our study, extended to other components of SWI/SNF, uncovers a key role of the SWI/SNF complex, in particular cBAF, in helping resolve R-loop-mediated T-R conflicts, thus, unveiling a new mechanism by which chromatin remodeling protects genome integrity.

Transcription shapes DNA replication initiation to preserve genome integrity

BACKGROUND

Early DNA replication occurs within actively transcribed chromatin compartments in mammalian cells, raising the immediate question of how early DNA replication coordinates with transcription to avoid collisions and DNA damage.

RESULTS

We develop a high-throughput nucleoside analog incorporation sequencing assay and identify thousands of early replication initiation zones in both mouse and human cells. The identified early replication initiation zones fall in open chromatin compartments and are mutually exclusive with transcription elongation. Of note, early replication initiation zones are mainly located in non-transcribed regions adjacent to transcribed regions. Mechanistically, we find that RNA polymerase II actively redistributes the chromatin-bound mini-chromosome maintenance complex (MCM), but not the origin recognition complex (ORC), to actively restrict early DNA replication initiation outside of transcribed regions. In support of this finding, we detect apparent MCM accumulation and DNA replication initiation in transcribed regions due to anchoring of nuclease-dead Cas9 at transcribed genes, which stalls RNA polymerase II. Finally, we find that the orchestration of early DNA replication initiation by transcription efficiently prevents gross DNA damage.

CONCLUSION

RNA polymerase II redistributes MCM complexes, but not the ORC, to prevent early DNA replication from initiating within transcribed regions. This RNA polymerase II-driven MCM redistribution spatially separates transcription and early DNA replication events and avoids the transcription-replication initiation collision, thereby providing a critical regulatory mechanism to preserve genome stability.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

Monitoring genome-wide replication fork directionality by Okazaki fragment sequencing in mammalian cells

The ability to monitor DNA replication fork directionality at the genome-wide scale is paramount for a greater understanding of how genetic and environmental perturbations can impact replication dynamics in human cells. Here we describe a detailed protocol for isolating and sequencing Okazaki fragments from asynchronously growing mammalian cells, termed Okazaki fragment sequencing (Ok-seq), for the purpose of quantitatively determining replication initiation and termination frequencies around specific genomic loci by meta-analyses. Briefly, cells are pulsed with 5-ethynyl-2'-deoxyuridine (EdU) to label newly synthesized DNA, and collected for DNA extraction. After size fractionation on a sucrose gradient, Okazaki fragments are concentrated and purified before click chemistry is used to tag the EdU label with a biotin conjugate that is cleavable under mild conditions. Biotinylated Okazaki fragments are then captured on streptavidin beads and ligated to Illumina adapters before library preparation for Illumina sequencing. The use of Ok-seq to interrogate genome-wide replication fork initiation and termination efficiencies can be applied to all unperturbed, asynchronously growing mammalian cells or under conditions of replication stress, and the assay can be performed in less than 2 weeks.

Discovery and Protein Modeling Studies of Novel Compound Mutations Causing Resistance to Multiple Tyrosine Kinase Inhibitors in Chronic Myeloid Leukemia

OBJECTIVE

BCR-ABL fusion oncogene is the hallmark of chronic myeloid leukemia (CML), causing genomic instability which leads to accumulation of mutations in BCR-ABL as well as other genes. BCR-ABL mutations are the cause of tyrosine kinase inhibitors (TKIs) resistance in CML. Recently, compound BCR-ABL mutations have been reported to resist all FDA approved TKIs. Therefore, finding novel compound BCR-ABL mutations can help and clinically manage CML. Therefore, our objective was to find out novel drug-resistant compound BCR-ABL mutations in CML and carry out their protein modelling studies.

METHODOLOGY

Peripheral blood samples were collected from ten imatinib resistant CML patients receiving nilotinib treatment. BCR-ABL transcript mutations were investigated by employing capillary sequencing. Patient follow-up was carried out using European LeukemiaNet guidelines. Protein modeling  studies were carried out for new compound mutations using PyMol to see the effects of mutations at structural level.

RESULTS

A novel compound mutation (K245N mutation along with G250W mutation) and previously known T351I utation was detected in two of the nilotinib resistance CML patients respectively while in the rest of 8 nilotinib responders, no resistant mutations were detected. Protein modelling studies indicated changes in BCR-ABL mutant protein which may have negatively impacted its binding with nilotinib leading to drug resistance.

CONCLUSION

We report a novel nilotinib resistant BCR-ABL compound mutation (K245N along with G250W mutation) which impacts structural modification in BCR-ABL mutant protein leading to drug resistance. As compound mutations pose a new threat by causing resistance to all FDA approved tyrosine kinase inhibitors in BCR-ABL+ leukemias, our study opens a new direction for in vitro characterization of novel BCR-ABL compound mutations and their resistant to second  generation and third generation TKIs.

Identification of new driver and passenger mutations within APOBEC-induced hotspot mutations in bladder cancer

BACKGROUND

APOBEC-driven mutagenesis and functional positive selection of mutated genes may synergistically drive the higher frequency of some hotspot driver mutations compared to other mutations within the same gene, as we reported for FGFR3 S249C. Only a few APOBEC-associated driver hotspot mutations have been identified in bladder cancer (BCa). Here, we systematically looked for and characterised APOBEC-associated hotspots in BCa.

METHODS

We analysed 602 published exome-sequenced BCas, for part of which gene expression data were also available. APOBEC-associated hotspots were identified by motif-mapping, mutation signature fitting and APOBEC-mediated mutagenesis comparison. Joint analysis of DNA hairpin stability and gene expression was performed to predict driver or passenger hotspots. Aryl hydrocarbon receptor (AhR) activity was calculated based on its target genes expression. Effects of AhR knockout/inhibition on BCa cell viability were analysed.

RESULTS

We established a panel of 44 APOBEC-associated hotspot mutations in BCa, which accounted for about half of the hotspot mutations. Fourteen of them overlapped with the hotspots found in other cancer types with high APOBEC activity. They mostly occurred in the DNA lagging-strand templates and the loop of DNA hairpins. APOBEC-associated hotspots presented systematically a higher prevalence than the other mutations within each APOBEC-target gene, independently of their functional impact. A combined analysis of DNA loop stability and gene expression allowed to distinguish known passenger from known driver hotspot mutations in BCa, including loss-of-function mutations affecting tumour suppressor genes, and to predict new candidate drivers, such as AHR Q383H. We further characterised AHR Q383H as an activating driver mutation associated with high AhR activity in luminal tumours. High AhR activity was also found in tumours presenting amplifications of AHR and its co-receptor ARNT. We finally showed that BCa cells presenting those different genetic alterations were sensitive to AhR inhibition.

CONCLUSIONS

Our study identified novel potential drivers within APOBEC-associated hotspot mutations in BCa reinforcing the importance of APOBEC mutagenesis in BCa. It could allow a better understanding of BCa biology and aetiology and have clinical implications such as AhR as a potential therapeutic target. Our results also challenge the dogma that all hotspot mutations are drivers and mostly gain-of-function mutations affecting oncogenes.

Targeting BMP signaling in the bone marrow microenvironment of myeloid leukemia

The bone morphogenetic protein (BMP) pathway regulates the fate and proliferation of normal hematopoietic stem cells (HSC) as well as interactions with their niche. While BMP2 and BMP4 promote HSC differentiation, only BMP4 maintains HSC pool and favors interactions with their niche. In myeloid leukemia, we have identified intrinsic and extrinsic dysregulations of the BMP pathway in Chronic Myeloid Leukemia (CML) and Acute Myeloid leukemia (AML) responsible for leukemic stem cells (LSC) survival. In AML, BMP pathway alterations sustain and promote resistant immature-like leukemic cells by activating a new signaling cascade. Binding of BMP4 to BMPR1A leads to ΔNp73 expression, which in turn induces NANOG, altogether associated with a poor patient's prognosis. Despite efficient targeted therapies, like Tyrosine Kinase Inhibitors (TKI) in CML, many patients retain LSCs. Our laboratory demonstrated that the BMP pathway sustains a permanent pool of LSCs expressing high levels of BMPR1B receptor, that evolve upon treatment to progressively implement a BMP4 autocrine loop, leading to TKI-resistant cells. Single cell RNA-Seq analysis of TKI-persisting LSCs showed a co-enrichment of BMP with Jak2-signaling, quiescence and stem cell (SC) signatures. Using a new model of persisting LSCs, we recently demonstrated that BMPR1B+ cells display co-activated Smad1/5/8 and Stat3 pathways and could be targeted by blocking BMPR1B/Jak2 signal. Lastly, a specific BMPR1B inhibitor impaired BMP4-mediated LSC protection against TKIs. Altogether, data based on various studies including ours, indicate that BMP targeting could eliminate leukemic cells within a protective bone marrow microenvironment to efficiently impact residual resistance or persistence of LSCs in myeloid leukemia.

Genomic methods for measuring DNA replication dynamics

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 analog 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.

APOBEC-mediated Mutagenesis as a Likely Cause of FGFR3 S249C Mutation Over-representation in Bladder Cancer

FGFR3 is one of the most frequently mutated genes in bladder cancer and a driver of an oncogenic dependency. Here we report that only the most common recurrent FGFR3 mutation, S249C (TCC→TGC), represents an APOBEC-type motif and is probably caused by the APOBEC-mediated mutagenic process, accounting for its over-representation. We observed significant enrichment of the APOBEC mutational signature and overexpression of AID/APOBEC gene family members in bladder tumors with S249C compared to tumors with other recurrent FGFR3 mutations. Analysis of replication fork directionality suggests that the coding strand of FGFR3 is predominantly replicated as a lagging strand template that could favor the formation of hairpin structures, facilitating mutagenic activity of APOBEC enzymes. In vitro APOBEC deamination assays confirmed S249 as an APOBEC target. We also found that the FGFR3 S249C mutation was common in three other cancer types with an APOBEC mutational signature, but rare in urothelial tumors without APOBEC mutagenesis and in two diseases probably related to aging. PATIENT SUMMARY: We propose that APOBEC-mediated mutagenesis can generate clinically relevant driver mutations even within suboptimal motifs, such as in the case of FGFR3 S249C, one of the most common mutations in bladder cancer. Knowledge about the etiology of this mutation will improve our understanding of the molecular mechanisms of bladder cancer.