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.

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.

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.

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.

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.

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.

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.

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.

Laser Microirradiation to Study In Vivo Cellular Responses to Simple and Complex DNA Damage.

DNA damage induces specific signaling and repair responses in the cell, which is critical for protection of genome integrity. Laser microirradiation became a valuable experimental tool to investigate the DNA damage response (DDR) in vivo. It allows real-time high-resolution single-cell analysis of macromolecular dynamics in response to laser-induced damage confined to a submicrometer region in the cell nucleus. However, various laser conditions have been used without appreciation of differences in the types of damage induced. As a result, the nature of the damage is often not well characterized or controlled, causing apparent inconsistencies in the recruitment or modification profiles. We demonstrated that different irradiation conditions (i.e., different wavelengths as well as different input powers (irradiances) of a femtosecond (fs) near-infrared (NIR) laser) induced distinct DDR and repair protein assemblies. This reflects the type of DNA damage produced. This protocol describes how titration of laser input power allows induction of different amounts and complexities of DNA damage, which can easily be monitored by detection of base and crosslinking damages, differential poly (ADP-ribose) (PAR) signaling, and pathway-specific repair factor assemblies at damage sites. Once the damage conditions are determined, it is possible to investigate the effects of different damage complexity and differential damage signaling as well as depletion of upstream factor(s) on any factor of interest.

The Use of Laser Microirradiation to Investigate the Roles of Cohesins in DNA Repair.

In addition to their mitotic and transcriptional functions, cohesin plays critical roles in DNA damage response (DDR) and repair. Specifically, cohesin promotes homologous recombination (HR) repair of DNA double-strand breaks (DSBs), which is conserved from yeast to humans, and is a critical effector of ATM/ATR DDR kinase-mediated checkpoint control in mammalian cells. Optical laser microirradiation has been instrumental in revealing the damage site-specific functions of cohesin and, more recently, uncovering the unique role of cohesin-SA2, one of the two cohesin complexes uniquely present in higher eukaryotes, in DNA repair in human cells. In this review, we briefly describe what we know about cohesin function and regulation in response to DNA damage, and discuss the optimized laser microirradiation conditions used to analyze cohesin responses to DNA damage in vivo.

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.

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.

Combination of low level light therapy and nitrosyl-cobinamide accelerates wound healing.

Low level light therapy (LLLT) has numerous therapeutic benefits, including improving wound healing, but the precise mechanisms involved are not well established; in particular, the underlying role of cytochrome C oxidase (C-ox) as the primary photoacceptor and the associated biochemical mechanisms still require further investigation. We previously showed the nitric oxide (NO) donating drug nitrosyl-cobinamide (NO-Cbi) enhances wound healing through a cGMP/cGMP-dependent protein kinase/ERK1/2 mechanism. Here, we show that the combination of LLLT and NO-Cbi markedly improves wound healing compared to either treatment alone. LLLT-enhanced wound healing proceeded through an electron transport chain-C-ox-dependent mechanism with a reduction of reactive oxygen species and increased adenosine triphosphate production. C-ox was validated as the primary photoacceptor by three observations: increased oxygen consumption, reduced wound healing in the presence of sodium azide, and disassociation of cyanide, a known C-ox ligand, following LLLT. We conclude that LLLT and NO-Cbi accelerate wound healing through two independent mechanisms, the electron transport chain-C-ox pathway and cGMP signaling, respectively, with both resulting in ERK1/2 activation.

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.

Distinct functions of human cohesin-SA1 and cohesin-SA2 in double-strand break repair.

Cohesin is an essential multiprotein complex that mediates sister chromatid cohesion critical for proper segregation of chromosomes during cell division. Cohesin is also involved in DNA double-strand break (DSB) repair. In mammalian cells, cohesin is involved in both DSB repair and the damage checkpoint response, although the relationship between these two functions is unclear. Two cohesins differing by one subunit (SA1 or SA2) are present in somatic cells, but their functional specificities with regard to DNA repair remain enigmatic. We found that cohesin-SA2 is the main complex corecruited with the cohesin-loading factor NIPBL to DNA damage sites in an S/G(2)-phase-specific manner. Replacing the diverged C-terminal region of SA1 with the corresponding region of SA2 confers this activity on SA1. Depletion of SA2 but not SA1 decreased sister chromatid homologous recombination repair and affected repair pathway choice, indicating that DNA repair activity is specifically associated with cohesin recruited to damage sites. In contrast, both cohesin complexes function in the intra-S checkpoint, indicating that cell cycle-specific damage site accumulation is not a prerequisite for cohesin's intra-S checkpoint function. Our findings reveal the unique ways in which cohesin-SA1 and cohesin-SA2 participate in the DNA damage response, coordinately protecting genome integrity in human cells.

Nitrosyl-cobinamide (NO-Cbi), a new nitric oxide donor, improves wound healing through cGMP/cGMP-dependent protein kinase.

Nitric oxide (NO) donors have been shown to improve wound healing, but the mechanism is not well defined. Here we show that the novel NO donor nitrosyl-cobinamide (NO-Cbi) improved in vitro wound healing in several cell types, including an established line of lung epithelial cells and primary human lung fibroblasts. On a molar basis, NO-Cbi was more effective than two other NO donors, with the effective NO-Cbi concentration ranging from 3 to 10µM, depending on the cell type. Improved wound healing was secondary to increased cell migration and not cell proliferation. The wound healing effect of NO-Cbi was mediated by cGMP, mainly through cGMP-dependent protein kinase type I (PKGI), as determined using pharmacological inhibitors and activators, and siRNAs targeting PKG type I and II. Moreover, we found that Src and ERK were two downstream mediators of NO-Cbi's effect. We conclude that NO-Cbi is a potent inducer of cell migration and wound closure, acting via cGMP, PKG, Src, and extracellular signal regulated kinase (ERK).

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.

Condensin I recruitment to base damage-enriched DNA lesions is modulated by PARP1.

Condensin I is important for chromosome organization and segregation in mitosis. We previously showed that condensin I also interacts with PARP1 in response to DNA damage and plays a role in single-strand break repair. However, whether condensin I physically associates with DNA damage sites and how PARP1 may contribute to this process were unclear. We found that condensin I is preferentially recruited to DNA damage sites enriched for base damage. This process is dictated by PARP1 through its interaction with the chromosome-targeting domain of the hCAP-D2 subunit of condensin I.