Single-cell spatial transcriptomics reveals a dystrophic trajectory following a developmental bifurcation of myoblast cell fates in facioscapulohumeral muscular dystrophy.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in "FSHD-hi" and "FSHD-lo" states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

Single-nucleus RNA-seq identifies divergent populations of FSHD2 myotube nuclei.

FSHD is characterized by the misexpression of DUX4 in skeletal muscle. Although DUX4 upregulation is thought to be the pathogenic cause of FSHD, DUX4 is lowly expressed in patient samples, and analysis of the consequences of DUX4 expression has largely relied on artificial overexpression. To better understand the native expression profile of DUX4 and its targets, we performed bulk RNA-seq on a 6-day differentiation time-course in primary FSHD2 patient myoblasts. We identify a set of 54 genes upregulated in FSHD2 cells, termed FSHD-induced genes. Using single-cell and single-nucleus RNA-seq on myoblasts and differentiated myotubes, respectively, we captured, for the first time, DUX4 expressed at the single-nucleus level in a native state. We identified two populations of FSHD myotube nuclei based on low or high enrichment of DUX4 and FSHD-induced genes ("FSHD-Lo" and "FSHD Hi", respectively). FSHD-Hi myotube nuclei coexpress multiple DUX4 target genes including DUXA, LEUTX and ZSCAN4, and also upregulate cell cycle-related genes with significant enrichment of E2F target genes and p53 signaling activation. We found more FSHD-Hi nuclei than DUX4-positive nuclei, and confirmed with in situ RNA/protein detection that DUX4 transcribed in only one or two nuclei is sufficient for DUX4 protein to activate target genes across multiple nuclei within the same myotube. DUXA (the DUX4 paralog) is more widely expressed than DUX4, and depletion of DUXA suppressed the expression of LEUTX and ZSCAN4 in late, but not early, differentiation. The results suggest that the DUXA can take over the role of DUX4 to maintain target gene expression. These results provide a possible explanation as to why it is easier to detect DUX4 target genes than DUX4 itself in patient cells and raise the possibility of a self-sustaining network of gene dysregulation triggered by the limited DUX4 expression.

Relationship of DUX4 and target gene expression in FSHD myocytes.

Facioscapulohumeral dystrophy (FSHD) is associated with the upregulation of the DUX4 transcription factor and its target genes. However, low-frequency DUX4 upregulation in patient myocytes is difficult to detect and examining the relationship and dynamics of DUX4 and target gene expression has been challenging. Using RNAScope in situ hybridization with highly specific probes, we detect the endogenous DUX4 and target gene transcripts in situ in patient skeletal myotubes during 13-day differentiation in vitro. We found that the endogenous DUX4 transcripts primarily localize as foci in one or two nuclei as compared with the accumulation of the recombinant DUX4 transcripts in the cytoplasm. We also found the continuous increase of DUX4 and target gene-positive myotubes after Day 3, arguing against its expected immediate cytotoxicity. Interestingly, DUX4 and target gene expression become discordant later in differentiation with the increase of DUX4-positive/target gene-negative as well as DUX4-negative/target gene-positive myotubes. Depletion of DUX4-activated transcription factors, DUXA and LEUTX, specifically repressed a DUX4-target gene, KDM4E, later in differentiation, suggesting that after the initial activation by DUX4, target genes themselves contribute to the maintenance of downstream gene expression. Together, the study provides important new insights into the dynamics of the DUX4 transcriptional network in FSHD patient myocytes.

Striated myocyte structural integrity: Automated analysis of sarcomeric z-discs.

As sarcomeres produce the force necessary for contraction, assessment of sarcomere order is paramount in evaluation of cardiac and skeletal myocytes. The uniaxial force produced by sarcomeres is ideally perpendicular to their z-lines, which couple parallel myofibrils and give cardiac and skeletal myocytes their distinct striated appearance. Accordingly, sarcomere structure is often evaluated by staining for z-line proteins such as α-actinin. However, due to limitations of current analysis methods, which require manual or semi-manual handling of images, the mechanism by which sarcomere and by extension z-line architecture can impact contraction and which characteristics of z-line architecture should be used to assess striated myocytes has not been fully explored. Challenges such as isolating z-lines from regions of off-target staining that occur along immature stress fibers and cell boundaries and choosing metrics to summarize overall z-line architecture have gone largely unaddressed in previous work. While an expert can qualitatively appraise tissues, these challenges leave researchers without robust, repeatable tools to assess z-line architecture across different labs and experiments. Additionally, the criteria used by experts to evaluate sarcomeric architecture have not been well-defined. We address these challenges by providing metrics that summarize different aspects of z-line architecture that correspond to expert tissue quality assessment and demonstrate their efficacy through an examination of engineered tissues and single cells. In doing so, we have elucidated a mechanism by which highly elongated cardiomyocytes become inefficient at producing force. Unlike previous manual or semi-manual methods, characterization of z-line architecture using the metrics discussed and implemented in this work can quantitatively evaluate engineered tissues and contribute to a robust understanding of the development and mechanics of striated muscles.

Biphasic recruitment of TRF2 to DNA damage sites promotes non-sister chromatid homologous recombination repair.

TRF2 (TERF2) binds to telomeric repeats and is critical for telomere integrity. Evidence suggests that it also localizes to non-telomeric DNA damage sites. However, this recruitment appears to be precarious and functionally controversial. We find that TRF2 recruitment to damage sites occurs by a two-step mechanism: the initial rapid recruitment (phase I), and stable and prolonged association with damage sites (phase II). Phase I is poly(ADP-ribose) polymerase (PARP)-dependent and requires the N-terminal basic domain. The phase II recruitment requires the C-terminal MYB/SANT domain and the iDDR region in the hinge domain, which is mediated by the MRE11 complex and is stimulated by TERT. PARP-dependent recruitment of intrinsically disordered proteins contributes to transient displacement of TRF2 that separates two phases. TRF2 binds to I-PpoI-induced DNA double-strand break sites, which is enhanced by the presence of complex damage and is dependent on PARP and the MRE11 complex. TRF2 depletion affects non-sister chromatid homologous recombination repair, but not homologous recombination between sister chromatids or non-homologous end-joining pathways. Our results demonstrate a unique recruitment mechanism and function of TRF2 at non-telomeric DNA damage sites.

Scc1 sumoylation by Mms21 promotes sister chromatid recombination through counteracting Wapl.

DNA double-strand breaks (DSBs) fuel cancer-driving chromosome translocations. Two related structural maintenance of chromosomes (Smc) complexes, cohesin and Smc5/6, promote DSB repair through sister chromatid homologous recombination (SCR). Here we show that the Smc5/6 subunit Mms21 sumoylates multiple lysines of the cohesin subunit Scc1. Mms21 promotes cohesin-dependent small ubiquitin-like modifier (SUMO) accumulation at laser-induced DNA damage sites in S/G2 human cells. Cells expressing the nonsumoylatable Scc1 mutant (15KR) maintain sister chromatid cohesion during mitosis but are defective in SCR and sensitive to ionizing radiation (IR). Scc1 15KR is recruited to DNA damage sites. Depletion of Wapl, a negative cohesin regulator, rescues SCR defects of Mms21-deficient or Scc1 15KR-expressing cells. Expression of the acetylation-mimicking Smc3 mutant does not bypass the requirement for Mms21 in SCR. We propose that Scc1 sumoylation by Mms21 promotes SCR by antagonizing Wapl at a step after cohesin loading at DSBs and in a way not solely dependent on Smc3 acetylation.

Femtosecond near-infrared laser microirradiation reveals a crucial role for PARP signaling on factor assemblies at DNA damage sites.

Laser microirradiation is a powerful tool for real-time single-cell analysis of the DNA damage response (DDR). It is often found, however, that factor recruitment or modification profiles vary depending on the laser system employed. This is likely due to an incomplete understanding of how laser conditions/dosages affect the amounts and types of damage and the DDR. We compared different irradiation conditions using a femtosecond near-infrared laser and found distinct damage site recruitment thresholds for 53BP1 and TRF2 correlating with the dose-dependent increase of strand breaks and damage complexity. Low input-power microirradiation that induces relatively simple strand breaks led to robust recruitment of 53BP1 but not TRF2. In contrast, increased strand breaks with complex damage including crosslinking and base damage generated by high input-power microirradiation resulted in TRF2 recruitment to damage sites with no 53BP1 clustering. We found that poly(ADP-ribose) polymerase (PARP) activation distinguishes between the two damage states and that PARP activation is essential for rapid TRF2 recruitment while suppressing 53BP1 accumulation at damage sites. Thus, our results reveal that careful titration of laser irradiation conditions allows induction of varying amounts and complexities of DNA damage that are gauged by differential PARP activation regulating protein assembly at the damage site.

Single-nucleus RNA-seq of differentiating human myoblasts reveals the extent of fate heterogeneity.

Myoblasts are precursor skeletal muscle cells that differentiate into fused, multinucleated myotubes. Current single-cell microfluidic methods are not optimized for capturing very large, multinucleated cells such as myotubes. To circumvent the problem, we performed single-nucleus transcriptome analysis. Using immortalized human myoblasts, we performed RNA-seq analysis of single cells (scRNA-seq) and single nuclei (snRNA-seq) and found them comparable, with a distinct enrichment for long non-coding RNAs (lncRNAs) in snRNA-seq. We then compared snRNA-seq of myoblasts before and after differentiation. We observed the presence of mononucleated cells (MNCs) that remained unfused and analyzed separately from multi-nucleated myotubes. We found that while the transcriptome profiles of myoblast and myotube nuclei are relatively homogeneous, MNC nuclei exhibited significant heterogeneity, with the majority of them adopting a distinct mesenchymal state. Primary transcripts for microRNAs (miRNAs) that participate in skeletal muscle differentiation were among the most differentially expressed lncRNAs, which we validated using NanoString. Our study demonstrates that snRNA-seq provides reliable transcriptome quantification for cells that are otherwise not amenable to current single-cell platforms. Our results further indicate that snRNA-seq has unique advantage in capturing nucleus-enriched lncRNAs and miRNA precursors that are useful in mapping and monitoring differential miRNA expression during cellular differentiation.

Dynamic regulation of the PR-Set7 histone methyltransferase is required for normal cell cycle progression.

Although the PR-Set7/Set8/KMT5a histone H4 Lys 20 monomethyltransferase (H4K20me1) plays an essential role in mammalian cell cycle progression, especially during G2/M, it remained unknown how PR-Set7 itself was regulated. In this study, we discovered the mechanisms that govern the dynamic regulation of PR-Set7 during mitosis, and that perturbation of these pathways results in defective mitotic progression. First, we found that PR-Set7 is phosphorylated at Ser 29 (S29) specifically by the cyclin-dependent kinase 1 (cdk1)/cyclinB complex, primarily from prophase through early anaphase, subsequent to global accumulation of H4K20me1. While S29 phosphorylation did not affect PR-Set7 methyltransferase activity, this event resulted in the removal of PR-Set7 from mitotic chromosomes. S29 phosphorylation also functions to stabilize PR-Set7 by directly inhibiting its interaction with the anaphase-promoting complex (APC), an E3 ubiquitin ligase. The dephosphorylation of S29 during late mitosis by the Cdc14 phosphatases was required for APC(cdh1)-mediated ubiquitination of PR-Set7 and subsequent proteolysis. This event is important for proper mitotic progression, as constitutive phosphorylation of PR-Set7 resulted in a substantial delay between metaphase and anaphase. Collectively, we elucidated the molecular mechanisms that control PR-Set7 protein levels during mitosis, and demonstrated that its orchestrated regulation is important for normal mitotic progression.

The annealing helicase HARP is recruited to DNA repair sites via an interaction with RPA.

HepA-related protein (HARP) (also known as SMARCAL1) is an ATP-driven annealing helicase that catalyzes the formation of dsDNA from complementary Replication protein A (RPA)-bound ssDNA. Here we find that HARP contains a conserved N-terminal motif that is necessary and sufficient for binding to RPA. This RPA-binding motif is not required for annealing helicase activity, but is essential for the recruitment of HARP to sites of laser-induced DNA damage. These findings suggest that the interaction of HARP with RPA increases the concentration of annealing helicase activity in the vicinity of ssDNA regions to facilitate processes such as DNA repair.

Chromatin dynamics during DNA repair revealed by pair correlation analysis of molecular flow in the nucleus.

Chromatin dynamics modulate DNA repair factor accessibility throughout the DNA damage response. The spatiotemporal scale upon which these dynamics occur render them invisible to live cell imaging. Here we present a believed novel assay to monitor the in vivo structural rearrangements of chromatin during DNA repair. By pair correlation analysis of EGFP molecular flow into chromatin before and after damage, this assay measures millisecond variations in chromatin compaction with submicron resolution. Combined with laser microirradiation we employ this assay to monitor the real-time accessibility of DNA at the damage site. We find from comparison of EGFP molecular flow with a molecule that has an affinity toward double-strand breaks (Ku-EGFP) that DNA damage induces a transient decrease in chromatin compaction at the damage site and an increase in compaction to adjacent regions, which together facilitate DNA repair factor recruitment to the lesion with high spatiotemporal control.

Genetic and epigenetic characteristics of FSHD-associated 4q and 10q D4Z4 that are distinct from non-4q/10q D4Z4 homologs.

Facioscapulohumeral dystrophy (FSHD) is one of the most prevalent muscular dystrophies. The majority of FSHD cases are linked to a decreased copy number of D4Z4 macrosatellite repeats on chromosome 4q (FSHD1). Less than 5% of FSHD cases have no repeat contraction (FSHD2), most of which are associated with mutations of SMCHD1. FSHD is associated with the transcriptional derepression of DUX4 encoded within the D4Z4 repeat, and SMCHD1 contributes to its regulation. We previously found that the loss of heterochromatin mark (i.e., histone H3 lysine 9 tri-methylation (H3K9me3)) at D4Z4 is a hallmark of both FSHD1 and FSHD2. However, whether this loss contributes to DUX4 expression was unknown. Furthermore, additional D4Z4 homologs exist on multiple chromosomes, but they are largely uncharacterized and their relationship to 4q/10q D4Z4 was undetermined. We found that the suppression of H3K9me3 results in displacement of SMCHD1 at D4Z4 and increases DUX4 expression in myoblasts. The DUX4 open reading frame (ORF) is disrupted in D4Z4 homologs and their heterochromatin is unchanged in FSHD. The results indicate the significance of D4Z4 heterochromatin in DUX4 gene regulation and reveal the genetic and epigenetic distinction between 4q/10q D4Z4 and the non-4q/10q homologs, highlighting the special role of the 4q/10q D4Z4 chromatin and the DUX4 ORF in FSHD.

Engineered FSHD mutations results in D4Z4 heterochromatin disruption and feedforward DUX4 network activation.

Facioscapulohumeral dystrophy (FSHD) is linked to contraction of D4Z4 repeats on chromosome 4q with SMCHD1 mutations acting as a disease modifier. D4Z4 heterochromatin disruption and abnormal upregulation of the transcription factor DUX4, encoded in the D4Z4 repeat, are the hallmarks of FSHD. However, defining the precise effect of D4Z4 contraction has been difficult because D4Z4 repeats are primate-specific and DUX4 expression is very rare in highly heterogeneous patient myocytes. We generated isogenic mutant cell lines harboring D4Z4 and/or SMCHD1 mutations in a healthy human skeletal myoblast line. We found that the mutations affect D4Z4 heterochromatin differently, and that SMCHD1 mutation or disruption of DNA methylation stabilizes otherwise variegated DUX4 target activation in D4Z4 contraction mutant cells, demonstrating the critical role of modifiers. Our study revealed amplification of the DUX4 signal through downstream targets, H3.X/Y and LEUTX. Our results provide important insights into how rare DUX4 expression leads to FSHD pathogenesis.

Cohesin mediates chromatin interactions that regulate mammalian ß-globin expression.

The ß-globin locus undergoes dynamic chromatin interaction changes in differentiating erythroid cells that are thought to be important for proper globin gene expression. However, the underlying mechanisms are unclear. The CCCTC-binding factor, CTCF, binds to the insulator elements at the 5' and 3' boundaries of the locus, but these sites were shown to be dispensable for globin gene activation. We found that, upon induction of differentiation, cohesin and the cohesin loading factor Nipped-B-like (Nipbl) bind to the locus control region (LCR) at the CTCF insulator and distal enhancer regions as well as at the specific target globin gene that undergoes activation upon differentiation. Nipbl-dependent cohesin binding is critical for long-range chromatin interactions, both between the CTCF insulator elements and between the LCR distal enhancer and the target gene. We show that the latter interaction is important for globin gene expression in vivo and in vitro. Furthermore, the results indicate that such cohesin-mediated chromatin interactions associated with gene regulation are sensitive to the partial reduction of Nipbl caused by heterozygous mutation. This provides the first direct evidence that Nipbl haploinsufficiency affects cohesin-mediated chromatin interactions and gene expression. Our results reveal that dynamic Nipbl/cohesin binding is critical for developmental chromatin organization and the gene activation function of the LCR in mammalian cells.

The acetyltransferase activity of San stabilizes the mitotic cohesin at the centromeres in a shugoshin-independent manner.

Proper sister chromatid cohesion is critical for maintaining genetic stability. San is a putative acetyltransferase that is important for sister chromatid cohesion in Drosophila melanogaster, but not in budding yeast. We showed that San is critical for sister chromatid cohesion in HeLa cells, suggesting that this mechanism may be conserved in metazoans. Furthermore, although a small fraction of San interacts with the NatA complex, San appears to mediate cohesion independently. San exhibits acetyltransferase activity in vitro, and its activity is required for sister chromatid cohesion in vivo. In the absence of San, Sgo1 localizes correctly throughout the cell cycle. However, cohesin is no longer detected at the mitotic centromeres. Furthermore, San localizes to the cytoplasm in interphase cells; thus, it may not gain access to chromosomes until mitosis. Moreover, in San-depleted cells, further depletion of Plk1 rescues the cohesion along the chromosome arms, but not at the centromeres. Collectively, San may be specifically required for the maintenance of the centromeric cohesion in mitosis.

Specific loss of histone H3 lysine 9 trimethylation and HP1gamma/cohesin binding at D4Z4 repeats is associated with facioscapulohumeral dystrophy (FSHD).

Facioscapulohumeral dystrophy (FSHD) is an autosomal dominant muscular dystrophy in which no mutation of pathogenic gene(s) has been identified. Instead, the disease is, in most cases, genetically linked to a contraction in the number of 3.3 kb D4Z4 repeats on chromosome 4q. How contraction of the 4qter D4Z4 repeats causes muscular dystrophy is not understood. In addition, a smaller group of FSHD cases are not associated with D4Z4 repeat contraction (termed "phenotypic" FSHD), and their etiology remains undefined. We carried out chromatin immunoprecipitation analysis using D4Z4-specific PCR primers to examine the D4Z4 chromatin structure in normal and patient cells as well as in small interfering RNA (siRNA)-treated cells. We found that SUV39H1-mediated H3K9 trimethylation at D4Z4 seen in normal cells is lost in FSHD. Furthermore, the loss of this histone modification occurs not only at the contracted 4q D4Z4 allele, but also at the genetically intact D4Z4 alleles on both chromosomes 4q and 10q, providing the first evidence that the genetic change (contraction) of one 4qD4Z4 allele spreads its effect to other genomic regions. Importantly, this epigenetic change was also observed in the phenotypic FSHD cases with no D4Z4 contraction, but not in other types of muscular dystrophies tested. We found that HP1gamma and cohesin are co-recruited to D4Z4 in an H3K9me3-dependent and cell type-specific manner, which is disrupted in FSHD. The results indicate that cohesin plays an active role in HP1 recruitment and is involved in cell type-specific D4Z4 chromatin regulation. Taken together, we identified the loss of both histone H3K9 trimethylation and HP1gamma/cohesin binding at D4Z4 to be a faithful marker for the FSHD phenotype. Based on these results, we propose a new model in which the epigenetic change initiated at 4q D4Z4 spreads its effect to other genomic regions, which compromises muscle-specific gene regulation leading to FSHD pathogenesis.

Comparative analysis of different laser systems to study cellular responses to DNA damage in mammalian cells.

Proper recognition and repair of DNA damage is critical for the cell to protect its genomic integrity. Laser microirradiation ranging in wavelength from ultraviolet A (UVA) to near-infrared (NIR) can be used to induce damage in a defined region in the cell nucleus, representing an innovative technology to effectively analyze the in vivo DNA double-strand break (DSB) damage recognition process in mammalian cells. However, the damage-inducing characteristics of the different laser systems have not been fully investigated. Here we compare the nanosecond nitrogen 337 nm UVA laser with and without bromodeoxyuridine (BrdU), the nanosecond and picosecond 532 nm green second-harmonic Nd:YAG, and the femtosecond NIR 800 nm Ti:sapphire laser with regard to the type(s) of damage and corresponding cellular responses. Crosslinking damage (without significant nucleotide excision repair factor recruitment) and single-strand breaks (with corresponding repair factor recruitment) were common among all three wavelengths. Interestingly, UVA without BrdU uniquely produced base damage and aberrant DSB responses. Furthermore, the total energy required for the threshold H2AX phosphorylation induction was found to vary between the individual laser systems. The results indicate the involvement of different damage mechanisms dictated by wavelength and pulse duration. The advantages and disadvantages of each system are discussed.

DNA damage induced during mitosis undergoes DNA repair synthesis.

Understanding the mitotic DNA damage response (DDR) is critical to our comprehension of cancer, premature aging and developmental disorders which are marked by DNA repair deficiencies. In this study we use a micro-focused laser to induce DNA damage in selected mitotic chromosomes to study the subsequent repair response. Our findings demonstrate that (1) mitotic cells are capable of DNA repair as evidenced by DNA synthesis at damage sites, (2) Repair is attenuated when DNA-PKcs and ATM are simultaneously compromised, (3) Laser damage may permit the observation of previously undetected DDR proteins when damage is elicited by other methods in mitosis, and (4) Twenty five percent of mitotic DNA-damaged cells undergo a subsequent mitosis. Together these findings suggest that mitotic DDR is more complex than previously thought and may involve factors from multiple repair pathways that are better understood in interphase.

Lysosome-targeting turn-on red/NIR BODIPY probes for imaging hypoxic cells.

A small series of fluorescent lysosome-targeting probes based on the BODIPY fluorophore and containing morpholine and nitro groups were rationally designed. These probes emitted light from green to NIR wavelengths, and provided specificity for imaging the lysosomes of hypoxic cells. The electron-withdrawing nitrophenyl group at the meso position was found to lead to highly efficient nonradiative decay of the S1 state, and hence a recovery of fluorescence when reduction of the nitro group occurred under hypoxic conditions.

NAD+ consumption by PARP1 in response to DNA damage triggers metabolic shift critical for damaged cell survival.

DNA damage signaling is critical for the maintenance of genome integrity and cell fate decision. Poly(ADP-ribose) polymerase 1 (PARP1) is a DNA damage sensor rapidly activated in a damage dose- and complexity-dependent manner playing a critical role in the initial chromatin organization and DNA repair pathway choice at damage sites. However, our understanding of a cell-wide consequence of its activation in damaged cells is still limited. Using the phasor approach to fluorescence lifetime imaging microscopy and fluorescence-based biosensors in combination with laser microirradiation, we found a rapid cell-wide increase of the bound NADH fraction in response to nuclear DNA damage, which is triggered by PARP-dependent NAD+ depletion. This change is linked to the metabolic balance shift to oxidative phosphorylation (oxphos) over glycolysis. Inhibition of oxphos, but not glycolysis, resulted in parthanatos due to rapid PARP-dependent ATP deprivation, indicating that oxphos becomes critical for damaged cell survival. The results reveal the novel prosurvival response to PARP activation through a change in cellular metabolism and demonstrate how unique applications of advanced fluorescence imaging and laser microirradiation-induced DNA damage can be a powerful tool to interrogate damage-induced metabolic changes at high spatiotemporal resolution in a live cell.