
==== Front
Nucleic Acids Res
Nucleic Acids Res
nar
Nucleic Acids Research
0305-1048
1362-4962
Oxford University Press

38994563
10.1093/nar/gkae600
gkae600
AcademicSubjects/SCI00010
Gene regulation, Chromatin and Epigenetics
SMCHD1 activates the expression of genes required for the expansion of human myoblasts
Wong Matthew Man-Kin Sprott Center for Stem Cell Research, Regenerative Medicine Program, Ottawa Hospital Research Institute; Ottawa, ON K1H 8L6, Canada
Department of Cellular and Molecular Medicine, University of Ottawa; Ottawa, ON K1H 8L6, Canada

Hachmer Sarah Department of Cell and Regenerative Biology, University of Wisconsin; Madison, WI 53705, USA

Gardner Ed Sprott Center for Stem Cell Research, Regenerative Medicine Program, Ottawa Hospital Research Institute; Ottawa, ON K1H 8L6, Canada
Department of Cellular and Molecular Medicine, University of Ottawa; Ottawa, ON K1H 8L6, Canada

Runfola Valeria Division of Genetics and Cell Biology, IRCCS San Raffaele Scientific Institute, Milano 20132, Italy

Arezza Eric Sprott Center for Stem Cell Research, Regenerative Medicine Program, Ottawa Hospital Research Institute; Ottawa, ON K1H 8L6, Canada

https://orcid.org/0000-0002-6824-8569
Megeney Lynn A Sprott Center for Stem Cell Research, Regenerative Medicine Program, Ottawa Hospital Research Institute; Ottawa, ON K1H 8L6, Canada
Department of Cellular and Molecular Medicine, University of Ottawa; Ottawa, ON K1H 8L6, Canada

Emerson Charles P Jr Wellstone Muscular Dystrophy Program, Department of Neurology, University of Massachusetts Chan Medical School, Worcester, MA 01655, USA

https://orcid.org/0000-0002-3811-4121
Gabellini Davide Division of Genetics and Cell Biology, IRCCS San Raffaele Scientific Institute, Milano 20132, Italy

https://orcid.org/0000-0002-6265-8102
Dilworth F Jeffrey Sprott Center for Stem Cell Research, Regenerative Medicine Program, Ottawa Hospital Research Institute; Ottawa, ON K1H 8L6, Canada
Department of Cellular and Molecular Medicine, University of Ottawa; Ottawa, ON K1H 8L6, Canada
Department of Cell and Regenerative Biology, University of Wisconsin; Madison, WI 53705, USA

To whom correspondence should be addressed. Tel: +1 608 265 3758; Email: fdilworth@wisc.edu
09 9 2024
12 7 2024
12 7 2024
52 16 94509462
27 6 2024
1 6 2024
22 11 2023
© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

SMCHD1 is an epigenetic regulatory protein known to modulate the targeted repression of large chromatin domains. Diminished SMCHD1 function in muscle fibers causes Facioscapulohumeral Muscular Dystrophy (FSHD2) through derepression of the D4Z4 chromatin domain, an event which permits the aberrant expression of the disease-causing gene DUX4. Given that SMCHD1 plays a broader role in establishing the cellular epigenome, we examined whether loss of SMCHD1 function might affect muscle homeostasis through additional mechanisms. Here we show that acute depletion of SMCHD1 results in a DUX4-independent defect in myoblast proliferation. Genomic and transcriptomic experiments determined that SMCHD1 associates with enhancers of genes controlling cell cycle to activate their expression. Amongst these cell cycle regulatory genes, we identified LAP2 as a key target of SMCHD1 required for the expansion of myoblasts, where the ectopic expression of LAP2 rescues the proliferation defect of SMCHD1-depleted cells. Thus, the epigenetic regulator SMCHD1 can play the role of a transcriptional co-activator for maintaining the expression of genes required for muscle progenitor expansion. This DUX4-independent role for SMCHD1 in myoblasts suggests that the pathology of FSHD2 may be a consequence of defective muscle regeneration in addition to the muscle wasting caused by spurious DUX4 expression.

Graphical Abstract

Graphical Abstract

ERA-Net for Research on Rare Diseases Canadian Institutes of Health Research 10.13039/501100000024 FDN-143330
==== Body
pmcIntroduction

The structural maintenance of chromosomes flexible hinge domain containing 1 (SMCHD1) gene encodes an epigenetic modifier first identified in a screen for genes involved in transgene array silencing (1). SMCHD1 is best known for its role in X chromosome inactivation (XCI) where it plays an essential role in merging the intermediate compartments to ensure proper gene silencing (2,3). SMCHD1 contributes to the establishment of heterochromatin by regulating DNA methylation, where it is required for Dnmt3b-mediated methylation of CpG islands on the inactive X chromosome (4), while also inhibiting the activity of Tet DNA demethylases (5). SMCHD1 can also actively repress gene transcription by limiting long-range chromatin interactions where it acts to deplete CTCF and cohesin from the topologically associated domain (TAD) boundaries at specific loci, including the Hox and protocadherin gene clusters (3,6–8).

Loss of function mutations in the SMCHD1 gene are observed in a subgroup of patients with the muscle-wasting condition Facioscapulohumeral Muscular Dystrophy (FSHD), one of the most frequent neuromuscular diseases. However, the predominant cause of FSHD is a deletion of the subtelomeric D4Z4 tandem DNA array on chromosome 4. This contraction of the D4Z4 locus (10 or less repeats) induces an epigenetic de-repression of the locus, which results in the ectopic expression of an RNA encoding double homeobox 4 (DUX4), a transcription factor whose presence in muscle fibers induces apoptosis (9–11). Based on this mechanism of action, it was later demonstrated that mutations reducing SMCHD1 function in muscle could lead to the de-repression of a slightly longer D4Z4 tandem array (up to 16 repeats) to permit the expression of DUX4 in patients with a variant of the disease called FSHD2 (12–14). Confirming the importance of Dux4 expression as causal in FSHD, the onset of the muscle-wasting disease requires a patient to have a specific 4qA allele, which introduces a polyA signal downstream of the last copy of the D4Z4 array to permit the expression of stable DUX4 mRNA (15,16). Indeed, individuals that have a 4qB D4Z4 allele will not develop FSHD despite possessing a D4Z4 array contraction or an SMCHD1 mutation that opens chromatin at the Dux4 gene. These facts clearly establish that regardless of the mechanism of D4Z4 locus de-repression, spurious accumulation of the DUX4 transcription factor protein is causal in FSHD. Nevertheless, patients who have mutations of SMCHD1 in addition to a contraction of the D4Z4 locus often display a more severe muscle-wasting phenotype (17,18). This raises the question of the potential for the general epigenetic silencer SMCHD1 to have additional roles in muscle cells that might exasperate the disease phenotype.

Here, we set out to investigate the function of SMCHD1 in human myoblasts that might affect the health of muscle independently of DUX4 expression. For this, we took advantage of the fact that the susceptibility of individuals to FSHD depends upon a specific polymorphism within the D4Z4 locus since the 4qA chromosome variant creates a polyadenylation signal which is necessary to stabilize DUX4 transcripts (16,19). We thus depleted SMCHD1 in cells derived from a healthy individual with a 4qB polymorphism at the D4Z4 locus on both chromosomes that cannot give rise to stable DUX4 transcripts (20). Using this approach, we uncovered a previously unappreciated DUX4-independent role for SMCHD1 as a co-activator of genes involved in cell cycle progression during myogenic progenitor cell expansion.

Materials and methods

Reagents and software

The names, sources, and catalogue numbers of the reagents used in the experiments are listed in Supplementary Table S1. All the software and databases used in the data analyses are listed in Supplementary Table S2.

Cell culture

Immortalized non-FSHD myoblast cell line WS234 (donor 15Vbic) (20) was obtained from the Wellstone Muscular Dystrophy Cooperative Research Center for FSHD, University of Massachusetts Medical School, USA. This cell line was previously tested to possess a 4qB D4Z4 chromosome variant on both chromosomes and have no detectable DUX4 mRNA expression (21). In proliferating condition, cells were cultured in LHCN medium [4:1 DMEM : Medium 199 supplemented with 15% BGS (v/v), 0.03 ug/ml Zinc Sulfate, 1.4 ug/ml vitamin B12, 0.055 ug/ml dexamethasone, 2.5 ng/ml HGF, 10 ng/ml FGF, 0.6× penicillin/streptomycin, 20 mM HEPES (pH 7.4)] on plates coated with 0.1% gelatin (w/v) at 37°C in 5% CO2 incubator. HEK293T cells were cultured in DMEM growth medium [DMEM supplemented with 10% BGS and 1× penicillin/streptomycin], at 37°C in 5% CO2 incubator.

Plasmid generation

The plasmid vector pDONR223 containing the open reading frame of human LAP2β was purchased from the genome editing and molecular biology facility of the University of Ottawa. The open reading frame of LAP2β was amplified by PCR using the forward primer 5′-catggaggatccgccaccatgccggagttcctggaag-3′ and the reverse primer 5′-cgtatggtacctcattagttggattttctagggtca-3′. The resulting PCR product was cloned in the lentiviral vector pLenti-EF1a-Blank and the sequence of LAP2β was confirmed by Sanger sequencing.

Preparation of lentivirus

HEK293T cells were transfected with the pMD2.G, psPAX2 and shRNA vector plasmids (Supplementary Table S3) using polyethyleneimine as previously described (22). Lentivirus-containing medium supernatant was centrifuged at 126 086 × g for 2.5 h at 4°C using Optima L-100 XP ultracentrifuge (Beckman Coulter), and the pellet was resuspended in LHCN medium. To transduce the cells, purified lentivirus and 8 μg/ml of polybrene was added to the culture medium. The cells were incubated for 16 hours at 37°C and then fresh LHCN medium was replaced.

Growth curve, EdU proliferation essay

The cell proliferation rate was determined by EdU incorporation in newly synthesized DNA as well as cell counting. Three days after lentiviral transduction, an equivalent number of cells were seeded on each well of a gelatin-coated six-well plate. Fresh LHCN medium was replaced every 2 days. Five days after lentiviral transduction, 10 μM EdU was added to the culture medium for 90 min at 37°C. Cells were fixed to perform immunostaining according to the manufacturer's protocol of the EdU staining kit. For cell counting experiments, cells were trypsinized and resuspended LHCN medium. Trypan blue solution was added and the mixture was loaded to a hemocytometer for cell counting. The cell number from each well was calculated from the average of 6 independent hemocytometer counting.

RNA-Seq and data analysis

Five days after SMCHD1 or non-silencing lentiviral transduction, RNA was collected using a PureLink RNA mini kit according to the manufacturer's protocol. The depletion of rRNA was performed using the Ribo-Zero Magnetic Gold Kit, and the RNA library was prepared using KAPA stranded RNA-Seq kit according to the manufacturer's protocol. Sequencing was performed on an Illumina HiSeq4000 platform. RNA-Seq reads from fastq files were aligned to human hg38 genome using the STAR (23), counts were summarized using featureCounts function (24) in Rsubread and used as input for DESeq2 (25) to identify the differentially expressed genes. RNA-Seq datasets (26) downloaded from GEO were analyzed using the same settings as the SMCHD1 depletion RNA-Seq dataset. Volcano plots were made using EnhancedVolcano (https://github.com/kevinblighe/EnhancedVolcano), dot plots were made using ggplot2 (27) and heatmaps using pheatmap (https://github.com/raivokolde/pheatmap) in R. Gene ontology analysis was performed using DAVID (28), and values were used as inputs for dot plots.

Gene set enrichment analysis

The down-regulated or up-regulated genes after SMCHD1 knockdown (adjusted P <= 0.05) were used as two separate gene set for Gene Set Enrichment Analysis (GSEA) (29). The analysis was performed using the GSEA software according to the user guide with the number of permutations set to 1000.

Binding and expression target analysis (BETA)

The analysis was done following the user manual (30). Peaks within 150 kb of TSS and differentially expressed genes with adjusted P value ⇐ 0.05 were selected for the analysis, without using CTCF as a boundary to filter peaks around genes.

ChIP-Seq and data analysis

The SMCHD1 ChIP-Seq was performed as previously described (31) with some modifications. Briefly, proliferating WS234 myoblasts were fixed in 1% (w/v) formaldehyde for 20 min at room temperature and quenched with 125 mM glycine for 10 min at room temperature. Cells were rinsed with PBS, collected using cell scrapers, centrifuged, and resuspended in swelling buffer [50 mM Tris base (pH 8.0), 300 mM sucrose, 10 mM NaCl, 0.5% NP-40 (v/v)]. Cells were incubated for 30 min in swelling buffer at 4°C, and sonicated using Bioruptor (Diagenode) for 2 cycles (HI power, cycles of 30 s on/1 min off). Nuclei were spun down at 1000 xg for 5 min at 4°C, resuspended in sonication buffer [50 mM HEPES (pH 7.9), 140 mM NaCl, 1 mM EDTA, 1% Triton X-100 (v/v), 0.1% Na-deoxycholate (w/v),1% SDS (w/v)], and sonicated for 13 cycles using the same settings as the previous step. Next, 7.5 ug of SMCHD1 antibody or normal rabbit IgG was added to 250 ug of precleared chromatin and incubated for 3 h at room temperature. The chromatin was then washed and purified as described previously (31). Library was prepared using Kapa Hyperprep kit and paired-end sequencing was performed on an Illumina HiSeq4000 platform. The reads were aligned to the human hg38 genome using bowtie2 (32) and peaks were called by MACS2 (33) using a q value cutoff = 0.05. Meta profile plots of ChIP-Seq signal around SMCHD1 peaks were generated using EnrichedHeatmap and ggplot2 in R.

CUT&Tag and data analysis

Analysis of histone modification enrichment across the genome was performed using CUT&Tag as previously described (34). Antibodies against either H3K4me3, H3K4me1, H3K9me3, H3K27ac or H3K27me3 (or control of rabbit IgG) were incubated with digitonin permeabilized WS234 myoblasts. After washing away the unbound antibodies, permeabilized cells were then incubated with the Tn5-PA protein. After a subsequent wash, cells were incubated with MgCl2 to induce tagmentation of the chromatin at sites enriched for the specific histone marks of interest. Tagmented DNA was amplified by PCR to generate a library for sequencing on an Illumina NovaSeq 6000 with 100 bp paired-end reads to a depth of 10M reads per library. Reads were aligned to the human hg38 genome using Bowtie2 (32) and peaks were called by MACS2 (33) using a q value cutoff = 0.05.

Motif identification using HOMER (Hypergeometric Optimization of Motif EnRichment)

The bed files generated by MACS2 were converted to fasta files using ‘getfasta’ function in BEDTools (35). The converted fasta files were used as input files for HOMER analysis (36) using the ‘findMotifs.pl’ function.

Statistical analyses

For RNA-Seq, differentially expressed genes were defined as genes with a fold change >0.5 and an adjusted P-value <= 0.05, unless otherwise specified. For the identification of the ChIP-Seq peak identification using MACS2 (33), the cutoff value q value cutoff = 0.05. For overlapping ChIP-Seq peaks between two datasets, P values and Z scores were calculated by permutation test using regioneR (37) with number of permutations = 1000. Two-way ANOVA (for cell counting experiments) and Student's two-tailed t-test (for other experiments) were performed using the statistical graphing software GraphPad Prism v6.0, and *P< 0.05, **P< 0.01, ***P< 0.001, ****P< 0.0001, ns = not significant.

Results

SMCHD1 is required for the efficient proliferation of human myoblasts

To determine the roles of SMCHD1 that are independent of DUX4 in adult myogenesis, we performed an acute depletion of SMCHD1 in immortalized human WS234 myoblasts (20). WS234 cells were chosen as they possess the non-pathogenic 4qB D4Z4 allele on both copies of chromosome 4, preventing the expression of a stable DUX4 mRNA transcript (21). We opted to use a shRNA-mediated approach for depletion to prevent the selection of clones that had adapted to the complete loss of SMCHD1. Using three different shRNAs targeting distinct regions of the SMCHD1 transcript, we observed a reduction in mRNA and protein accumulation that varied between 50 and 90% when examined 5 days after lentivirus infection (Supplementary Figure S1A and S1B). These reduced levels of SMCHD1 would be consistent with the haploinsufficiency in FSHD2 patients possessing a heterozygous mutation of the gene. Phenotypically, the reduction in SMCHD1 levels in cells caused the muscle progenitor cells to adopt a less elongated morphology compared to the controls (Figure 1A). This change in cell appearance suggests that SMCHD1 may be important for myoblast function.

Figure 1. Loss of SMCHD1 leads to FSHD-like repression of PAX7 target genes in myoblasts without DUX4 target genes activation. (A) Nine days after lentiviral transduction, cells were observed under the fluorescence microscope. (White arrows) Transduced mCherry+SMCHD1 shRNA expressing cells. (Black arrows) Cells that do not express shRNA expression. (B) Volcano plot showing the DEGs after SMCHD1 knockdown. Genes with adjusted p-value <= 0.05 and log2 fold change > 0.5 were labeled red. (C) DUX4 score of samples transduced with non-silencing scrambled shRNA (Control) or SMCHD1 targeting shRNA (shSMCHD1) based on DUX4 targets identified by Choi et al. (39). (D) PAX7 scores of samples transduced with non-silencing scrambled shRNA (Control) or SMCHD1 targeting shRNA (shSMCHD1). (E) Venn diagram of the target genes activated by PAX7 and differentially expressed upon SMCHD1 depletion. (F) Number of PAX7-activated genes activated by PAX7 that were down-regulated (red) and up-regulated (blue) upon SMCHD1 depletion. (G) UCSC genome browser tracts showing RNA-Seq reads of samples transduced with non-silencing scrambled shRNA (Control) or SMCHD1 targeting shRNA (shSMCHD1) near the MYF5 gene. (H) Heatmap showing the expression of DUX4 target genes identified by Yao et. al. (41) upon SMCHD1 depletion. (I) UCSC genome browser tracts showing RNA-Seq reads of samples transduced with non-silencing scrambled shRNA (Control) or SMCHD1 targeting shRNA (shSMCHD1) near the DUX4 gene.

Loss of SMCHD1 leads to aberrant expression of PAX7 target genes in myoblasts

To understand how the loss of SMCHD1 impairs myoblast function, we performed RNA-sequencing analysis on myoblasts depleted of SMCHD1 using two distinct shRNAs. Analysis of gene expression changes in myoblasts depleted of SMCHD1 identified 1631 differentially expressed genes (932 down-regulated and 699 up-regulated) (Figure 1B and Supplementary File S1). Though we were using myoblasts that possess only the 4qB D4Z4 allele, we wanted to confirm the absence of DUX4 in our myoblasts. Indeed, no transcripts of DUX4 were detected in our deep-sequencing of RNA from WS234 cells under any conditions tested. Due to the difficulty of detecting DUX4 transcripts in muscle cells, measurement of DUX4-target gene expression is often used as a surrogate for detecting the presence of DUX4 (38). Measuring the previously described ‘DUX4 score’ (38), we observed that SMCHD1 knockdown did not significantly change the expression of previously established DUX4 biomarkers (39–41) (Figure 1C and Supplementary Figure S1C). Reasoning that the loss of SMCHD1 may only change the expression of the most robust DUX4-FSHD biomarkers, we re-assessed a DUX4 score using only the top biomarkers identified by Yao et al. (41). Again we found that the depletion of SMCHD1 did not significantly change the expression of the four core sets of DUX4-FSHD biomarkers (Figure 1H). Thus, the proliferation defects observed in the absence of SMCHD1 are not the result of induced DUX4 expression.

A hallmark of myoblasts from FSHD patients is the repression of Pax7 target genes (42). As such, we next examined whether the expression of known Pax7 target genes was altered upon depletion of SMCHD1. For this purpose, we calculated a PAX7 score (38) based on changes in transcript abundance upon SMCHD1 depletion (Figure 1D). We observed that ∼50% of the genes normally activated by Pax7 show reduced expression in myoblasts depleted of SMCHD1 (Figure 1E–G), while ∼11% of genes repressed by PAX7 expression show increased expression upon SMCHD1 depletion (Supplementary Figure S1D, E). Therefore, our data showed that SMCHD1 depletion reproduces some of the altered gene expression as observed in myoblasts-derived from FSHD patients despite the fact that DUX4 itself is not expressed (Figure 1I).

To better understand how the loss of SMCHD1 was affecting myoblast function, we performed gene ontology (GO) enrichment analysis. Consistent with results observed in FSHD patients with SMCHD1 mutations (43), loss of the epigenetic modifier in myoblasts results in the upregulation of genes associated with the extracellular matrix (Figure 2A and Supplementary Figure S2A). Though SMCHD1 is known as an epigenetic silencer, we were intrigued by the large number of genes involved in cell cycle (DNA replication, nucleosome assembly and G1/S transition) that were down-regulated upon loss of SMCHD1 (Figure 2B). Using RT-qPCR, we confirmed that SMCHD1-depletion leads to the downregulation of the genes coding for cyclins and cyclin-dependent kinases (Supplementary Figure S2B), centromeres (Supplementary Figure S2C) and core histone proteins (Supplementary Figure S2D). These findings suggest that loss of SMCHD1 reduces the expression of genes associated with cell cycle progression and may affect the proliferation of myoblasts. To explore this possibility, we examined the expansion of SMCHD1-depleted myoblasts (Figure 2C-D). Measuring the rate of EdU incorporation into newly synthesized DNA, we observed that loss of SMCHD1 resulted in fewer myoblasts passing through S-phase of the cell cycle (Figure 2C). This decreased passage through cell cycle resulted in a reduced expansion of the progenitor population as determined by direct counting of cell numbers over time (Figure 2D). Taken together this data indicates that the loss of SMCHD1 in human myoblasts leads to a decreased ability of the progenitor population to expand their numbers.

Figure 2. Loss of SMCHD1 in myoblasts leads to a decrease in proliferation. (A) GO terms of up-regulated genes after SMCHD1 depletion. (B) GO terms of down-regulated genes after SMCHD1 depletion. (C) Five days after lentiviral transduction, an EdU cell proliferation assay was performed, and the percentage of EdU+ cells was quantified. Error bars represent standard deviation from three independent experiments. (D) An equal number of myoblasts was seeded 3 days after lentiviral infection, and cell number was counted from day 5 to day 8 after lentiviral transduction. Error bars represent standard deviation from 3 independent experiments.

Having observed that loss of SMCHD1 results in FSHD-like changes in PAX7 target gene expression, we next examined whether the down-regulation of cell cycle genes observed upon depletion of SMCHD1 in myoblasts might be an unappreciated defect in myoblasts from FSHD patients with SMCHD1 mutations. Published datasets from FSHD patients with D4Z4 contractions (42,44) or SMCHD1 mutations (26) were analyzed in combination with RNA-Seq data from our SMCHD1-depleted myoblast using Gene Set Enrichment Analysis (GSEA). Since these studies looked at differentiation of primary myoblasts, we limited our comparative analysis to the 0h time point (ie proliferation) of these published datasets (26,42,44). Interestingly, we observed that genes that change their expression (both up- or down-regulated) in myoblasts upon loss of SMCHD1 overlapped significantly with genes that change their expression in FSHD myoblasts compared to healthy controls (Figure 3A and Supplementary Figure S3A). This strong correlation between FSHD patient myoblasts and myoblasts depleted of SMCHD1 suggested these cells may share similar defects, and led us to examine whether the expression of cell cycle genes was also dysregulated in FSHD patients. GO analysis showed that FSHD2 myoblasts from patients have altered expression of genes involved in cell division, and DNA replication (Supplementary Figure S3B), as well as negative regulation of cell proliferation (Supplementary Figure S3C). Unexpectedly, we observed that the changes in gene expression related to cell cycle and cell proliferation were specific to FSHD2 (Figure 3B), but not in FSHD1 myoblasts (Supplementary Figure S3D-F). The defect in cell cycle gene expression could also be observed (Figure 3C) in biopsies taken from FSHD2 patients (43,45). Thus, down-regulation of genes involved in cell cycle progression is a characteristic shared between myoblasts depleted of SMCHD1 (myoblasts with shRNA-mediated depletion of SMCDH1 or primary FSHD2 patient myoblasts), but not myoblasts from FSHD patients that retain SMCHD1 (FSHD1). The similar changes in cell cycle gene expression between SMCHD1-depleted myoblasts and in FSHD2 patient samples was intriguing given that previous studies have demonstrated that the doubling time of healthy and FSHD1 myoblasts were similar (46,47).

Figure 3. Comparison of changes in gene expression profiles upon SMCHD1 depletion and under FSHD2 condition. (A) Gene set enrichment analysis (GSEA) of DEGs in FSHD2 myoblasts compared to down-regulated (left panel) or up-regulated (right panel) genes upon SMCHD1 depletion. (B) (Left panel) Venn diagram of commonly down-regulated genes and (Right panel) dot plot of the top GO terms of the genes commonly down-regulated in SMCHD1-depleted myoblasts and in FSHD2 myoblasts. (C) (Left panel) Upset plot showing the commonly down-regulated genes upon SMCHD1 depletion and in FSHD2 biopsies from Wang et al. (45) and Wong et al. (43) studies. (Right panel) GO terms of commonly down-regulated genes upon SMCHD1 depletion and in FSHD2 biopsies.

SMCHD1 binds at transcriptional regulatory regions of cell cycle genes in human myoblasts

SMCHD1 is known as a transcriptional silencer, yet the loss of SMDHD1 results in the down-regulation of cell cycle genes. This led us to examine whether genes down-regulated by the loss of SMCHD1 might be directly targeted by the epigenetic regulatory protein. Using chromatin immunoprecipitation sequencing (ChIP-Seq), we probed the genome-wide binding of SMCHD1 in WS234 cells. We identified 15 560 sites (q value <= 0.05) of binding of SMCHD1 in proliferating myoblasts (Supplementary File S2). Knowing that SMCHD1 is involved in X chromosome inactivation (2), we first examined the distribution of SMCHD1 across the different chromosomes. As expected, we found that SMCHD1 binding was enriched on the X chromosome, with ∼4 fold more SMCHD1 binding sites relative to individual autosomes (15.1 sites per Mb on X chromosome versus 3.6 sites per Mb on autosomes) (Supplementary Figure S4A). Similarly, we observed a strong enrichment of SMCHD1 (Supplementary Figure S4B-S4D) at sites previously known to be targeted by SMCHD1, including telomeres (48) as well as the 5S rRNA and tRNA clusters on chromosome 1 (49). This established that SMCHD1 in human myoblasts localizes to many of the regions previously described for other cell types and provides confidence in the specificity of our ChIP-Seq.

We next wanted to associate SMCHD1 binding at specific genes to changes in gene expression in order to identify direct target genes. For this purpose, we performed Binding and Expression Target Analysis (BETA) which directly integrates data from ChIP-Seq and RNA-Seq (30). Using this approach, we identified 358 genes bound by SMCHD1 that changed their expression upon SMCHD1 depletion (Supplementary File S2). Though SMCHD1 has known silencing activity, we observed that many target genes are repressed upon loss of SMCHD1 (181 down-regulated and 177 up-regulated). We then asked whether SMCHD1’s binding motifs are different in those sites associated with down-regulated versus up-regulated targets. We found that the majority of the enriched motifs were similar between the two conditions (Supplementary Figure S5) with some exceptions. For example, the ZBTB3 motif of CCACTVYAC was specifically enriched at the peaks associated with SMCHD1 co-activated genes (Supplementary Figure S5A), while the ZNF416 motif of CCCAGCT was specifically enriched at peaks associated with SMCHD1 co-repressed genes (Supplementary Figure S5B). To identify epigenetic characteristics associated with SMCHD1 binding sites, we performed CUT&Tag analysis for various histone marks (H3K27ac, H3K9me3, H3K4me3, H3K4me1 and H3K27me3) in proliferating WS234 myoblasts. Globally we observed that SMCHD1 preferentially bound at heterochromatic regions enriched for H3K9me3 marks. This finding was expected and is consistent with the binding of SMCHD1 to the telomers and the inactive X chromosome, and work showing that the binding of SMCHD1 at the D4Z4 locus is H3K9me3 dependent (50). However, when analysis was restricted to SMCHD1 peaks associated with genes that show reduced expression upon SMCHD1-depletion, we observed an enrichment for the transcriptionally permissive H3K4me3 mark but no enrichment of the repressive marks H3K9me3 or H3K27me3 (Figure 4A). This finding led us to examine whether SMCHD1 might be associated with transcriptional regulatory regions. Using the human GeneHancer database (51), our analysis determined that 38% of SMCHD1 direct target genes (Z score = 22.7) showed SMCHD1 binding at previously identified promoters or enhancers (Figure 4B). Consistent with previous findings (52) we observed the co-binding of CTCF at many of these SMCHD1-bound enhancers (Supplementary Figure S5C). This indicates a role for SMCHD1 in helping establish active enhancers at genes involved in cell cycle regulation.

Figure 4. Binding regions of SMCHD1 in the human myoblast genome. (A) The enrichment of H3K9me3, H3K4me1, H3K4me3, H3K27me3 and H3K27ac around SMCHD1 peak centers at SMCHD1 co-activated genes. (B) Percentage of direct targets-associated peaks that were located at promoters and enhancers. Z score and P-value were calculated by permutation test using regioneR. (C) UCSC genome browser tracts showing ENCODE cis-regulatory elements, GeneHancer annotated enhancer/promoter regions (¶) and their regulatory interactions, as well as ChIP-Seq (SMCHD1), CUT&Tag (H3K4me3 and H3K27ac) and RNA-Seq reads (control: proliferating myoblasts expressing non-silencing scrambled shRNA, shSMCHD1: proliferating myoblasts expressing shRNA targeting SMCHD1) near the LAP2 gene.

SMCHD1 regulates the proliferation of myoblasts via its direct target gene LAP2

Having identified a list of genes involved in cell cycle regulation as direct targets of SMCHD1, we next sought to identify genes that contribute to the cell proliferation defect in SMCHD1-depleted cells. Amongst these cell cycle regulators directly regulated by SMCHD1, we focused on Lamina-associated polypeptide 2 (LAP2, also known as Thymopoietin or TMPO) (Figure 4C), since our analysis of published RNA-Seq data (26) showed the gene was down-regulated in proliferating primary myoblasts isolated from FSHD2 patients (Figure 5A) while being unaltered in myoblasts from FSHD1 patients (Supplementary Figure S6A). Among the different isoforms of the LAP2 gene (Figure 5B), LAP2α is a lamin A/C binding protein that has been shown to be essential for the proliferation of glioblastoma cells (53), while LAP2β is essential for nuclear envelope reassembly after mitosis in xenopus (54). We therefore tested whether the loss of LAP2 would also lead to a proliferation defect in human myoblasts. We were able to obtain a strong depletion of LAP2 expression in myoblasts using two different shRNA targeting LAP2 (Figure 5B, and Supplementary Figure S6B). Examining cell morphology at 9 days after shRNA treatment, we observed that myoblasts depleted of LAP2 became more flattened and spread out compared to control cells (expressing non-silencing shRNA) which continued to maintain a normal myoblast morphology (Figure 5C). Counting the cell numbers as the population proliferated over 10 days, myoblasts depleted in LAP2 showed a reduced rate of proliferation (Figure 5D). This indicates that the effect of losing SMCHD1 on the proliferation and morphology of myoblasts can be recapitulated by the loss of LAP2. Examination of exon inclusion from RNA-seq data indicated that LAP2β is the predominant isoform of the gene expressed in proliferating human myoblasts (Figure 4C). To determine whether LAP2 could rescue the SMCHD1 depletion-mediated proliferation defect, we overexpressed LAP2β in myoblasts depleted of SMCHD1. Excitingly, ectopic expression of LAP2β was sufficient to significantly improve the proliferation of SMCHD1-depleted myoblasts (Figure 5E). In summary, our data suggest that SMCHD1 helps ensure the efficient proliferation of human myoblasts by activating the expression of its target genes involved in cell cycle progression, including LAP2.

Figure 5. LAP2 depletion slows myoblast proliferation. (A) Heatmap showing the RNA expression of LAP2 in FSHD2 and control non-FSHD myoblasts (GEO accession number GSE143493). (B) Different LAP2 isoforms and target regions of different LAP2 shRNAs. (C) Nine days after LAP2 shRNA transduction, cells were observed under a fluorescence microscope for GFP+ shRNA-expressing cells. (D) Three days after LAP2 depletion, an equal number of myoblasts were seeded on 6 well plates. Cell number was counted on day 4, 6, 8 and 10 after lentiviral transduction. Error bars represent standard deviation from 3 independent experiments. (E) Lentivirus was added to the culture medium of proliferating myoblasts as indicated (control: non-silencing scrambled shRNA, shSMCHD1: shRNA targeting SMCHD1, vector: plenti-GIII vector alone, LAP2β: plenti-GIII vector expressing LAP2β transcript). Three days after lentiviral transduction, equal number of myoblasts were seeded on six-well plates. Cell number was counted on days 4, 6 and 10 after lentiviral transduction. Error bars represent standard deviation from three independent experiments.

Discussion

The detrimental effects of reduced SMCHD1 levels in mature muscle myofibers have been shown to be at the heart of the FSHD disease pathology in FSHD2 patients. However, the importance of the ubiquitously expressed SMCHD1 protein in maintaining and expanding the myoblast population has yet to be fully established. Here, we show that SMCHD1 is required for the efficient expansion of muscle progenitor cells. We found that myoblasts depleted of SMCHD1 show a reduced rate of proliferation that correlates with reduced expression of cell cycle-related genes at the transcription level. While these changes in gene expression occur in the absence of DUX4 expression, it is notable that the same dysregulation of cell cycle gene expression can be observed in muscle progenitor cells obtained from FSHD2 patients with SMCHD1 mutations. This is consistent with FSHD being a proposed satellite cell-opathy (55,56) and suggests that impaired expansion of myoblasts during regeneration may contribute to the disease pathology in FSHD2 patients.

The application of a knockdown approach to diminish SMCHD1 levels has allowed us to delineate a role for the chromatin-modifying enzyme in regulating the proliferation of human muscle progenitor cells. The reduced expansion of myoblasts upon reduction of SMCHD1 levels is associated with the down-regulation of genes involved in DNA replication and cell cycle progression. Consistent with the role of SMCHD1 in promoting cell proliferation, it has been shown that SMCHD1 is required for the expression of cell cycle genes in the early mouse embryo and that depletion of maternal SMCHD1 in the oocyte reduces the number of cells present in the blastocyst (57). Similarly, deletion of SMCHD1 in the hematopoietic system resulted in hematopoietic stem cells (HSCs) that were inefficient at competing with wild-type HSCs for repopulating the bone marrow niche (58). This suggests that SMCHD1 plays a role in the expansion of various stem cell populations.

We identified LAP2 as a key SMCHD1-target gene that ensures efficient proliferation of the myoblasts. Indeed, SMCHD1 binds to an enhancer for LAP2, and is required for efficient expression of this protein that links the nuclear envelope to chromatin. The LAP proteins have been shown to play an important role in cell proliferation through the regulation of Rb/E2F1 activity as well as through the re-organization of nuclear envelope-associated proteins after cell division (59). Importantly, we demonstrated that exogenous expression of the LAP2β isoform was able to partially rescue the proliferation defect observed in conditions of diminished SMCHD1 levels. Our findings are particularly relevant for patients with FSHD2 since the diminished expression of LAP2 was also observed in individuals carrying SMCHD1 mutations. Although it is generally accepted that FSHD2 patients do not present unique clinical phenotypes (60), we observed that muscle progenitor cells from FSHD2 patients show reduced expression of many of the cell cycle-related genes that were identified upon SMCHD1 depletion. As such, patients with mutations in SMCHD1 might expect to show diminished regeneration efficiency in addition to the loss of muscle due to spurious DUX4 expression. Though patients with FSHD1 have been shown to have defective regeneration due to impaired muscle differentiation (56), our identification of this DUX4-independent role for SMCHD1 in muscle regeneration (Figure 6), provides an additional explanation for previous findings that SMCHD1 mutations can act as a disease modifier that increases the disease severity FSHD1 patients with D4Z4 locus contractions that cause muscle wasting (17). In these patients, not only would the muscle fibers undergo wasting but the muscle would also have a reduced ability to regenerate due to poor expansion of the muscle progenitors. Thus, we have identified a key role for SMCHD1 in maintaining the expression of genes involved in the proliferation of muscle progenitor cells.

Figure 6. Working model of how the loss of SMCHD1 function results in FSHD2-associated defects. (Left panel) Loss of SMCHD1 binding on the promoter and enhancer regions of its target genes in myoblasts may result in additional problems to myoblasts growth, on top of those caused by spurious DUX4 expression. (Right panel) Heterozygous SMCHD1 mutations in FSHD2 patients lead to a reduced SMCHD1 binding on the D4Z4 locus and the expression of DUX4 in muscle, which results in the degeneration of myofiber.

A surprising finding from our study was the determination that SMCHD1 plays the role of a transcriptional activator in certain contexts. SMCHD1 is classically known for its role in silencing large chromatin domains (1). During the onset of X-chromosome inactivation, recruitment of SMCHD1 to the intermediate compartments along the chromosome allows their fusion to form repressive megadomains that ensure gene silencing (3). As a non-classical SMC protein, it is thought that SMCHD1 creates these repressive megadomains by entrapping DNA in dense loops of chromatin (61). However, we and others have observed that depletion of SMCHD1 results in a large number of genes being down-regulated, a result consistent with a transcriptional co-activator function (this study; 57–58). As expected, genome-wide analysis of SMCHD1 binding showed an accumulation of the epigenetic modifier at sites-specific repressive domains in myoblasts. Unexpectedly, we also observed that SMCHD1 binds at transcriptional enhancers and/or promoters for many of the genes that are down-regulated upon removal of the epigenetic modifier. This is consistent with SMCHD1 acting as a transcriptional co-activator in certain genomic contexts. Of note, we found that ZBTB3 motifs are specifically enriched at SMCHD1 binding sites associated with co-activated genes. This is interesting since ZBTB3 is known as a structural regulator of enhancer-promoter interactions (62) where it co-activates genes to modulate the self-renewal of embryonic stem cells (63) and the proliferation of cancer cells (64). The enrichment of ZBTB3 motifs at these co-activated genes suggests that this site-specific DNA-binding factor could be responsible for recruiting SMCHD1 to selected genes, where SMCHD1 then acts as the ‘motor’ to mediate the loop extrusion that brings the enhancers in contact with the promoters to facilitate transcriptional elongation. Future studies will be required to examine the relationship between SMCHD1 and ZBTB3 in co-activating cell cycle genes in proliferating myoblasts. Importantly, SMCHD1 is not the first example in which an epigenetic regulator has been found to have a dual role in regulating gene expression. While the polycomb repressive complexes PRC1 and PRC2 play an important role in maintaining transcriptional repression at the D4Z4 locus (65–67), it has been shown that the binding of PRC1 (68–70) and PRC2 (71,72) at enhancers can also promote the activation of developmental gene expression by establishing long-range chromatin interactions. Thus, it appears that multiple different transcriptional repressors that mediate long-range interactions can be co-opted to facilitate the activation of gene expression in certain contexts.

In summary, we have identified a role for SMCHD1 in modulating the expansion of muscle progenitor cells through the activation of genes that control cell cycle progression. These findings suggest that impaired muscle regeneration may contribute to the disease pathology in FSHD2 patients with mutations that induce SMCHD1 haploinsufficiency.

Supplementary Material

gkae600_Supplemental_Files

Acknowledgements

We thank Odile Neyret for ChIP-Seq and RNA-Seq library sequencing; Sarah Boissel for CUT&Tag library sequencing. Hina Bandukwala for help with TelomereHunter analysis, and Lifang Li for technical help.

Data availability

The RNA-Seq and ChIP-Seq data sets generated during this study are available in Gene Expression Omnibus as GSE190830, while CUT&Tag data sets are available as GSE267354.

Supplementary data

Supplementary Data are available at NAR Online.

Funding

ERA-Net for Research on Rare Diseases [to D.G., F.J.D.]; Canadian Institutes of Health Research [FDN-143330 to F.J.D.]. Funding for open access charge: Canadian Institutes for Health Research.

Conflict of interest statement. None declared.

Notes

Present address: Matthew Man-Kin Wong, Department of Genetics, Stanford University, Palo Alto, CA 94304, USA.
==== Refs
References

1. Blewitt M.E. , VickaryousN.K., HemleyS.J., AsheA., BruxnerT.J., PreisJ.I., ArkellR., WhitelawE. An N-ethyl-N-nitrosourea screen for genes involved in variegation in the mouse. Proc. Natl. Acad. Sci. U.S.A. 2005; 102 :7629–7634.15890782
2. Blewitt M.E. , GendrelA.V., PangZ., SparrowD.B., WhitelawN., CraigJ., ApedaileA., HiltonD., DunwoodieS., BrockdorffN.et al . SMCHD1, containing a structural-maintenance-of-chromosomes hinge domain, has a critical role in X inactivation. Nat. Genet. 2008; 40 :663–669.18425126
3. Wang C.Y. , JéguT., ChuH.P., OhH.J., LeeJ.T. SMCHD1 merges chromosome compartments and assists formation of super-structures on the inactive X. Cell. 2018; 174 :406–421.29887375
4. Gendrel A.V. , ApedaileA., CokerH., TermanisA., ZvetkovaI. SMCHD1-dependent and -independent pathways determine developmental dynamics of CpG island methylation on the inactive X chromosome. Dev. Cell. 2012; 23 :265–279.22841499
5. Huang Z. , YuJ., CuiW., JohnsonB.K., KimK., PfeiferG.P. The chromosomal protein SMCHD1 regulates DNA methylation and the 2c-like state of embryonic stem cells by antagonizing TET proteins. Sci. Adv. 2021; 7 :eabb9149.33523915
6. Jansz N. , KeniryA., TrussartM., BildsoeH., BeckT., TonksI.D., MouldA., HickeyP., BreslinK., IminitoffM.et al . SMCHD1 regulates long-range chromatin interactions on the inactive X chromosome and at Hox clusters. Nat. Struct. Mol. Biol. 2018; 25 :766–777.30127357
7. Chen K. , HuJ., MooreD.L., LiuR., KessansS.A., BreslinK., LucetI.S., KeniryA., LeongH.S., ParishC.L.et al . Genome-wide binding and mechanistic analyses of SMCHD1-mediated epigenetic regulation. Proc. Natl. Acad. Sci. U.S.A. 2015; 112 :E3535–E3544.26091879
8. Gdula M.R. , NesterovaT.B., PintacudaG., GodwinJ., ZhanY., OzadamH., McClellanM., MoralliD., KruegerF., GreenC.M.et al . The non-canonical SMC protein SMCHD1 antagonises TAD formation and compartmentalisation on the inactive X chromosome. Nat. Commun. 2019; 10 :30.30604745
9. Kowaljow V. , MarcowyczA., AnsseauE., CondeC.B., SauvageS., MattéottiC., AriasC., CoronaE.D., NuñezN.G., LeoO.et al . The DUX4 gene at the FSHD1A locus encodes a pro-apoptotic protein. Neuromuscul. Disord. 2007; 17 :611–623.17588759
10. Wallace L.M. , GarwickS.E., MeiW., BelayewA., CoppeeF., LadnerK.J., GuttridgeD., YangJ., HarperS.Q. DUX4, a candidate gene for facioscapulohumeral muscular dystrophy, causes p53-dependent myopathy in vivo. Ann. Neurol. 2011; 69 :540–552.21446026
11. Bosnakovski D. , XuZ., GangE.J., GalindoC.L., LiuM., SimsekT., GarnerH.R., Agha-MohammadiS., TassinA., CoppéeF.et al . An isogenetic myoblast expression screen identifies DUX4-mediated FSHD-associated molecular pathologies. EMBO J. 2008; 27 :2766–2779.18833193
12. Lemmers R.J. , TawilR., PetekL.M., BalogJ., BlockG.J., SantenG.W.E., AmellA.M., van der VlietP.J., AlmomaniR., StraasheijmK.R.et al . Digenic inheritance of an SMCHD1 mutation and an FSHD-permissive D4Z4 allele causes facioscapulohumeral muscular dystrophy type 2. Nat. Genet. 2012; 44 :1370–1374.23143600
13. Lemmers R.J. , GoemanJ.J., van der VlietP.J., van NieuwenhuizenM.P., BalogJ., Vos-VersteegM., CamanoP., ArroyoM.A.R., JericoI., RogersM.T.et al . Inter-individual differences in CpG methylation at D4Z4 correlate with clinical variability in FSHD1 and FSHD2. Hum. Mol. Genet. 2015; 24 :659–669.25256356
14. Sacconi S. , Briand-SuleauA., GrosM., BaudoinC., LemmersR.J.L.F., RondeauS., LaghaN., NigumannP., CambieriC., PumaA.et al . FSHD1 and FSHD2 form a disease continuum. Neurology. 2019; 92 :e2273–e2285.30979860
15. Lemmers R.J. , de KievitP., SandkuijlL., PadbergG.W., van OmmenG.J., FrantsR.R., van der MaarelS.M. Facioscapulohumeral muscular dystrophy is uniquely associated with one of the two variants of the 4q subtelomere. Nat. Genet. 2002; 32 :235–236.12355084
16. Lemmers R.J. , WohlgemuthM., van der GaagK.J., van der VlietP.J., van TeijlingenC.M., de KnijffP., PadbergG.W., FrantsR.R., van der MaarelS.M. Specific sequence variations within the 4q35 region are associated with facioscapulohumeral muscular dystrophy. Am. J. Hum. Genet. 2007; 81 :884–894.17924332
17. Sacconi S. , LemmersR.J., BalogJ., van der VlietP.J., LahautP., van NieuwenhuizenM.P., StraasheijmK.R., DebipersadR.D., Vos-VersteegM., SalviatiL.et al . The FSHD2 gene SMCHD1 is a modifier of disease severity in families affected by FSHD1. Am. J. Hum. Genet. 2013; 93 :744–751.24075187
18. Larsen M. , RostS., El HajjN., FerbertA., DeschauerM., WalterM.C., SchoserB., TacikP., KressW., MüllerC.R. Diagnostic approach for FSHD revisited: SMCHD1 mutations cause FSHD2 and act as modifiers of disease severity in FSHD1. Eur. J. Hum. Genet. 2015; 23 :808–816.25370034
19. Lemmers R.J. , van der VlietP.J., KloosterR., SacconiS., CamañoP., DauwerseJ.G., SniderL., StraasheijmK.R., van OmmenG.J., PadbergG.W.et al . A unifying genetic model for facioscapulohumeral muscular dystrophy. Science. 2010; 329 :1650–1653.20724583
20. Homma S. , ChenJ.C., RahimovF., BeermannM.L., HangerK., BibatG.M., WagnerK.R., KunkelL.M., EmersonC.P.J., MillerJ.B A unique library of myogenic cells from facioscapulohumeral muscular dystrophy subjects and unaffected relatives: family, disease and cell function. Eur. J. Hum. Genet. 2012; 20 :404–410.22108603
21. Jones T.I. , ChenJ.C., RahimovF., HommaS., ArashiroP., BeermannM.L., KingO.D., MillerJ.B., KunkelL.M., EmersonC.P.J.et al . Facioscapulohumeral muscular dystrophy family studies of DUX4 expression: evidence for disease modifiers and a quantitative model of pathogenesis. Hum. Mol. Genet. 2012; 21 :4419–4430.22798623
22. Longo P.A. , KavranJ.M., KimM.S., LeahyD.J. Transient mammalian cell transfection with polyethylenimine (PEI). Methods Enzymol. 2013; 529 :227–240.24011049
23. Dobin A. , DavisC.A., SchlesingerF., DrenkowJ., ZaleskiC., JhaS., BatutP., ChaissonM., GingerasT.R. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29 :15–21.23104886
24. Liao Y. , SmythG.K., ShiW. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019; 47 :e47.30783653
25. Love M.I. , HuberW., AndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15 :550.25516281
26. Jiang S. , WilliamsK., KongX., ZengW., NguyenN.V., MaX., TawilR., YokomoriK., MortazaviA Single-nucleus RNA-seq identifies divergent populations of FSHD2 myotube nuclei. PLoS Genet. 2020; 16 :e1008754.32365093
27. Wickham H. ggplot2: Elegant Graphics for Data Analysis. 2016; NY Springer-Verlag.
28. Huang D. , ShermanB.T., LempickiR.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009; 4 :44–57.19131956
29. Subramanian A. , TamayoP., MoothaV.K., MukherjeeS., EbertB.L., GilletteM.A., PaulovichA., PomeroyS.L., GolubT.R., LanderE.S.et al . Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 2005; 102 :15545–15550.16199517
30. Wang S. , SunH., MaJ., ZangC., WangC., WangJ., TangQ., MeyerC.A., ZhangY., LiuX.S. Target analysis by integration of transcriptome and ChIP-seq data with BETA. Nat. Protoc. 2013; 8 :2502–2515.24263090
31. Singh K. , CassanoM., PlanetE., SebastianS., JangS.M., SohiG., FaralliH., ChoiJ., YounH.D., DilworthF.J.et al . A KAP1 phosphorylation switch controls MyoD function during skeletal muscle differentiation. Genes Dev. 2015; 29 :513–525.25737281
32. Langmead B. , SalzbergS.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods. 2012; 9 :357–359.22388286
33. Zhang Y. , LiuT., MeyerC.A., EeckhouteJ., JohnsonD.S., BernsteinB.E., NusbaumC., MyersR.M., BrownM., LiW.et al . Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008; 9 :R137.18798982
34. Li Y. , NakkaK., OlenderT., Gingras-GelinasP., WongM.M., RobinsonD.C.L., BandukwalaH., PaliiC.G., NeyretO., BrandM.et al . Chromatin and transcription factor profiling in rare stem cell populations using CUT&Tag. STAR Protoc. 2021; 2 :100751.34467227
35. Quinlan A.R. , HallI.M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26 :841–842.20110278
36. Heinz S. , BennerC., SpannN., BertolinoE., LinY.C., LasloP., ChengJ.X., MurreC., SinghH., GlassC.K.et al . Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell. 2010; 38 :576–589.20513432
37. Gel B. , Díez-VillanuevaA., SerraE., BuschbeckM., PeinadoM.A., MalinverniR. regioneR: an R/Bioconductor package for the association analysis of genomic regions based on permutation tests. Bioinformatics. 2016; 32 :289–291.26424858
38. Banerji C. , ZammitP.S. PAX7 target gene repression is a superior FSHD biomarker than DUX4 target gene activation, associating with pathological severity and identifying FSHD at the single-cell level. Hum. Mol. Genet. 2019; 28 :2224–2236.31067297
39. Choi S.H. , GearhartM.D., CuiZ., BosnakovskiD., KimM., SchennumN., KybaM. DUX4 recruits p300/CBP through its C-terminus and induces global H3K27 acetylation changes. Nucleic Acids Res. 2016; 44 :5161–5173.26951377
40. Geng L.N. , YaoZ., SniderL., FongA.P., CechJ.N., YoungJ.M., van der MaarelS.M., RuzzoW.L., GentlemanR.C., TawilR.et al . DUX4 activates germline genes, retroelements, and immune mediators: implications for facioscapulohumeral dystrophy. Dev. Cell. 2012; 22 :38–51.22209328
41. Yao Z. , SniderL., BalogJ., LemmersR.J., Van Der MaarelS.M., TawilR., TapscottS.J. DUX4-induced gene expression is the major molecular signature in FSHD skeletal muscle. Hum. Mol. Genet. 2014; 23 :5342–5352.24861551
42. Banerji C. , PanamarovaM., HebaishiH., WhiteR.B., RelaixF., SeveriniS., ZammitP.S. PAX7 target genes are globally repressed in facioscapulohumeral muscular dystrophy skeletal muscle. Nat. Commun. 2017; 8 :2152.29255294
43. Wong C.J. , WangL.H., FriedmanS.D., ShawD., CampbellA.E., BudechC.B., LewisL.M., LemmersR.J.F.L., StatlandJ.M., van der MaarelS.M.et al . Longitudinal measures of RNA expression and disease activity in FSHD muscle biopsies. Hum. Mol. Genet. 2020; 29 :1030–1043.32083293
44. Krom Y.D. , DumonceauxJ., MamchaouiK., den HamerB., MariotV., NegroniE., GengL.N., MartinN., TawilR., TapscottS.J.et al . Generation of isogenic D4Z4 contracted and noncontracted immortal muscle cell clones from a mosaic patient: a cellular model for FSHD. Am. J. Pathol. 2012; 181 :1387–1401.22871573
45. Wang L.H. , FriedmanS.D., ShawD., SniderL., WongC.J., BudechC.B., PoliachikS.L., GoveN.E., LewisL.M., CampbellA.E. MRI-informed muscle biopsies correlate MRI with pathology and DUX4 target gene expression in FSHD. Hum. Mol. Genet. 2019; 28 :476–486.30312408
46. Barro M. , CarnacG., FlavierS., MercierJ., VassetzkyY., Laoudj-ChenivesseD. Myoblasts from affected and non-affected FSHD muscles exhibit morphological differentiation defects. J. Cell. Mol. Med. 2010; 14 :275–289.18505476
47. Vilquin J.T. , MarolleauJ.P., SacconiS., GarcinI., LacassagneM.N., RobertI., TernauxB., BouazzaB., LargheroJ., DesnuelleC. Normal growth and regenerating ability of myoblasts from unaffected muscles of facioscapulohumeral muscular dystrophy patients. Gene Ther. 2005; 12 :1651–1662.15973444
48. Grolimund L. , AebyE., HamelinR., ArmandF., ChiappeD., MoniatteM., LingnerJ. A quantitative telomeric chromatin isolation protocol identifies different telomeric states. Nat. Commun. 2013; 4 :2848.24270157
49. Mason A.G. , SliekerR.C., BalogJ., LemmersR., WongC.J., YaoZ., LimJ.W., FilippovaG.N., NeE., TawilR.et al . SMCHD1 regulates a limited set of gene clusters on autosomal chromosomes. Skelet Muscle. 2017; 7 :12.28587678
50. Zeng W. , ChenY.Y., NewkirkD.A., WuB., BalogJ., KongX., BallA.R.Jr, ZanottiS., TawilR., HashimotoN.et al . Genetic and epigenetic characteristics of FSHD-associated 4q and 10q D4Z4 that are distinct from non-4q/10q D4Z4 homologs. Hum. Mutat. 2014; 35 :998–1010.24838473
51. Fishilevich S. , NudelR., RappaportN., HadarR., PlaschkesI., SteinT.I., RosenN., KohnA., TwikM., SafranM.et al . GeneHancer: genome-wide integration of enhancers and target genes in GeneCards. Database (Oxford). 2017; 2017 :bax028.28605766
52. Laberthonnière C. , DelourmeM., ChevalierR., DionC., GanneB., HirstD., CaronL., PerrinP., AdélaïdeJ., ChaffanetM.et al . In skeletal muscle and neural crest cells, SMCHD1 regulates biological pathways relevant for Bosma syndrome and facioscapulohumeral dystrophy phenotype. Nucleic Acids Res. 2023; 51 :7269–7287.37334829
53. Zhang L. , WangG., ChenS., DingJ., JuS., CaoH., TianH. Depletion of thymopoietin inhibits proliferation and induces cell cycle arrest/apoptosis in glioblastoma cells. World J. Surg. Oncol. 2016; 14 :267.27756319
54. Dubińska-Magiera M. , ChmielewskaM., KoziołK., MachowskaM., HutchisonC.J., GoldbergM.W., RzepeckiR. Xenopus LAP2β protein knockdown affects location of lamin B and nucleoporins and has effect on assembly of cell nucleus and cell viability. Protoplasma. 2016; 253 :943–956.26209045
55. Ganassi M. , ZammitP.S. Involvement of muscle satellite cell dysfunction in neuromuscular disorders: Expanding the portfolio of satellite cell-opathies. Eur. J. Transl. Myol. 2022; 32 :10064.35302338
56. Banerji C.R.S. , HendersonD., TawilR.N., ZammitP.S. Skeletal muscle regeneration in facioscapulohumeral muscular dystrophy is correlated with pathological severity. Hum. Mol. Genet. 2020; 29 :2746–2760.32744322
57. Midic U. , VincentK.A., WangK., LokkenA., SeveranceA.L., RalstonA., KnottJ.G., LathamK.E. Novel key roles for structural maintenance of chromosome flexible domain containing 1 (SMCHD1) during preimplantation mouse development. Mol. Reprod. Dev. 2018; 85 :635–648.29900695
58. Kinkel S.A. , LiuJ., BeckT., BreslinK.A., IminitoffM., HickeyP., BlewittM.E. Epigenetic modifier SMCHD1 maintains a normal pool of long-term hematopoietic stem cells. iScience. 2022; 25 :104684.35856023
59. Barton L.J. , SoshnevA.A., GeyerP.K. Networking in the nucleus: A spotlight on LEM-domain proteins. Curr. Opin. Cell Biol. 2015; 34 :1–8.25863918
60. de Greef J.C. , LemmersR.J., CamañoP., DayJ.W., SacconiS., DunandM., van EngelenB.G.M., Kiuru-EnariS., PadbergG.W., RosaA.L.et al . Clinical features of facioscapulohumeral muscular dystrophy 2. Neurology. 2010; 75 :1548–1554.20975055
61. Gurzau A.D. , BlewittM.E., CzabotarP.E., MurphyJ.M., BirkinshawR.W. Relating SMCHD1 structure to its function in epigenetic silencing. Biochem. Soc. Trans. 2020; 48 :1751–1763.32779700
62. Wang Z. , LiangQ., QianX., HuB., ZhengZ., WangJ., HuY., BaoZ., ZhaoK., ZhouY.et al . An autoimmune pleiotropic SNP modulates IRF5 alternative promoter usage through ZBTB3-mediated chromatin looping. Nat. Commun. 2023; 14 :1208.36869052
63. Ye B. , LiuB., YangL., ZhuX., ZhangD., WuW.x, ZhuP., WangY., WangS., XiaP.et al . LncKdm2b controls self-renewal of embryonic stem cells via activating expression of transcription factor Zbtb3. EMBO J. 2018; 37 :e97174.29535137
64. Lim J.H. Zinc finger and BTB domain-containing protein 3 is essential for the growth of cancer cells. BMB Rep. 2014; 47 :405–410.24856827
65. Haynes P. , BomsztykK., MillerD.G. Sporadic DUX4 expression in FSHD myocytes is associated with incomplete repression by the PRC2 complex and gain of H3K9 acetylation on the contracted D4Z4 allele. Epigenetics Chromatin. 2018; 11 :47.30122154
66. Casa V. , RunfolaV., MicheloniS., AzizA., DilworthF.J., GabelliniD. Polycomb repressive complex 1 provides a molecular explanation for repeat copy number dependency in FSHD muscular dystrophy. Hum. Mol. Genet. 2017; 26 :753–767.28040729
67. Cabianca D.S. , CasaV., BodegaB., XynosA., GinelliE., TanakaY., TanakaY., GabelliniD. A long ncRNA links copy number variation to a polycomb/trithorax epigenetic switch in FSHD muscular dystrophy. Cell. 2012; 149 :819–831.22541069
68. Morey L. , SantanachA., BlancoE., AloiaL., NoraE.P., BruneauB.G., Di CroceL. Polycomb regulates mesoderm cell fate-specification in embryonic stem cells through activation and repression mechanisms. Cell Stem Cell. 2015; 17 :300–315.26340528
69. Cohen I. , ZhaoD., BarC., ValdesV.J., Dauber-DeckerK.L., NguyenM.B., NakayamaM., RendlM., BickmoreW.A., KosekiH.et al . PRC1 fine-tunes gene repression and activation to safeguard skin development and stem cell specification. Cell Stem Cell. 2018; 22 :726–739.29727681
70. Loubiere V. , PapadopoulosG.L., SzaboQ., MartinezA.M., CavalliG. Widespread activation of developmental gene expression characterized by PRC1-dependent chromatin looping. Sci. Adv. 2020; 6 :eaax4001.31950077
71. Gentile C. , BerlivetS., MayranA., PaquetteD., Guerard-MilletF., BajonE., DostieJ., KmitaM. PRC2-associated chromatin contacts in the developing limb reveal a possible mechanism for the atypical role of PRC2 in HoxA gene expression. Dev. Cell. 2019; 50 :184–196.31204170
72. Gentile C. , KmitaM. Polycomb repressive complexes in Hox gene regulation: silencing and beyond: the functional dynamics of Polycomb repressive complexes in Hox gene regulation. Bioessays. 2020; 42 :e1900249.32743818
