Project description:Numerous anti-cancer drugs perturb thymidylate biosynthesis and lead to genomic uracil incorporation contributing to their antiproliferative effect. Still, it is not yet characterized if uracil incorporations have any positional preference. Here, we aimed to uncover genome-wide alterations in uracil pattern upon drug treatments in human cancer cell line models derived from HCT116. We developed a straightforward U-DNA sequencing method (U-DNA-Seq) that was combined with in situ super-resolution imaging. Using a novel robust analysis pipeline, we found broad regions with elevated probability of uracil occurrence both in treated and non-treated cells. Correlation with chromatin markers and other genomic features shows that non-treated cells possess uracil in the late replicating constitutive heterochromatic regions, while drug treatment induced a shift of incorporated uracil towards segments that are normally more active/functional. Data were corroborated by colocalization studies via dSTORM microscopy. This approach can be applied to study the dynamic spatio-temporal nature of genomic uracil.
Project description:We report a new immunoprecipitation-coupled sequencing (DIP-Seq) application termed U-DNA-Seq, where a tailored and catalytically inactive uracil-DNA glycosylase (UNG) was used as uracil-DNA sensor to immunoprecipitate uracil containing genomic DNA fragments. Genomic uracil was profiled in drug-treated (5-fluoro-2'-deoxyuridine (5FdUR) or raltitrexed (RTX)) or non-treated (NT) HCT116 cells expressing the UNG inhibitor (UGI). The same experiments were also performed in the mismatch repair proficient version of the HCT116 cells (HCT116MMR), where chromosome 3 is reinserted to restore functional MMR (PMID: 8044777). Moreover, wild-type HCT116 or K562 cells were also measured. We found that regions of uracil enrichment in this assay were rather broad as compared to the sharp peaks typical in ChIP-seq. Therefore, we applied an approach alternative to the conventional peak calling. Namely, we calculated genome scaled coverage tracks and log2 ratio tracks of the enriched versus the input samples using deepTools package (bamCoverage and bigwigCompare tools, respectively) to provide a more appropriate description of uracil-enriched genomic regions. Interval (bed) files were also derived from these log2 ratio tracks to be able to screen large datasets for colocalizing features with them. For wider context of the study, see the related publication.
