<figure><img src="image_1.png" /></figure>

HAL
open science

# Does DNA Methylation Matter in FSHD?

Valentina Salsi, Frédérique Magdinier, Rossella Tupler

► To cite this version:

Valentina Salsi, Frédérique Magdinier, Rossella Tupler. Does DNA Methylation Matter in FSHD?. Genes, 2020, 11 (3), pp.258. 10.3390/genes11030258 . hal-02503971

HAL Id: hal-02503971

https://amu.hal.science/hal-02503971

Submitted on 16 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Does DNA Methylation Matter in FSHD?

Valentina Salsi<fnref n="1" />, Frédérique Magdinier<fnref n="2" /> ID and Rossella Tupler<fnref n="3" />,4,5,6,* ID

Received: 23 January 2020; Accepted: 25 February 2020; Published: 28 February 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Facioscapulohumeral muscular dystrophy (FSHD) has been associated with the genetic and epigenetic molecular features of the CpG-rich D4Z4 repeat tandem array at 4q35. Reduced DNA methylation of D4Z4 repeats is considered part of the FSHD mechanism and has been proposed as a reliable marker in the FSHD diagnostic procedure. We considered the assessment of D4Z4 DNA methylation status conducted on distinct cohorts using different methodologies. On the basis of the reported results we conclude that the percentage of DNA methylation detected at D4Z4 does not correlate with the disease status. Overall, data suggest that in the case of FSHD1, D4Z4 hypomethylation is a consequence of the chromatin structure present in the contracted allele, rather than a proxy of its function. Besides, CpG methylation at D4Z4 DNA is reduced in patients presenting diseases unrelated to muscle progressive wasting, like Bosma Arhinia and Microphthalmia syndrome, a developmental disorder, as well as ICF syndrome. Consistent with these observations, the analysis of epigenetic reprogramming at the D4Z4 locus in human embryonic and induced pluripotent stem cells indicate that other mechanisms, independent from the repeat number, are involved in the control of the epigenetic structure at D4Z4.

**Keywords:** FSHD1; FSHD2; D4Z4 macrosatellite; DNA methylation; epigenetics

## 1. Introduction

Facioscapulohumeral muscular dystrophy (FSHD, OMIM#158900) is the third most common myopathy with a reported prevalence ranging between 1 in 8333 [1] and 1 in 20,000 [2]. FSHD is characterized by insidious onset and progressive wasting of a highly selective set of muscle groups, the facial, limb girdle, and foot extensor muscles, and by a great variability of clinical expression among patients and within families [3,4]. The disease appears significantly earlier in males [5-7] and this determines male patients to be in higher number and more severely affected than females [5,6,8-10]. In some families, individuals affected by FSHD can be found only in one generation [5,7,11,12]. These differences are striking in discordant monozygotic twins [13-15].

Genetically, FSHD has been considered a Mendelian disease with autosomal dominant inheritance twins [16-18]. On this basis, the FSHD genetic locus was mapped on chromosome 4q35 by genetic linkage analysis [16,18,19] and associated with rearrangements of an array of tandemly repeated 3.3 kb

<footnote n="1">¹ Department of Life Sciences, University of Modena and Reggio Emilia, 4, 41121 Modena, Italy; valentina.salsi@unimore.it</footnote>

<footnote n="2">² Aix Marseille Univ, INSERM, MMG, U 1251, 13005 Marseille, France; Frederique.MAGDINIER@univ-amu.fr</footnote>

<footnote n="3">³ Department of Biomedical, Metabolic and Neural Sciences, University of Modena and Reggio Emilia, 4, 41121 Modena, Italy</footnote>

<footnote n="4">⁴ Center for Neuroscience and Neurotechnology, University of Modena and Reggio Emilia, 4, 41121 Modena, Italy</footnote>

<footnote n="5">⁵ Department of Molecular Cell and Cancer Biology, University of Massachusetts Medical School, Worcester, MA 01003, USA</footnote>

<footnote n="6">⁶ Li Weibo Institute for Rare Diseases Research at the University of Massachusetts Medical School, Worcester, MA 01003, USA</footnote>

<footnote>* Correspondence: rossella.tupler@unimore.it; Tel.: +3905-9205-5414</footnote>
segments (D4Z4) [20]. FSHD is the only human disease causally linked to copy number variation of macrosatellite DNA elements [21]. The number of D4Z4 repeats varies from 11 to 100 in the general population, whereas 10 repeats or fewer are usually found in sporadic and familial FSHD patients. As a general rule, alleles composed by 11–100 copies of D4Z4 repeats constitute the normal size range of D4Z4 alleles, whereas alleles with 8 or fewer D4Z4 repeats are considered diagnostic for the disease [22,23].

The routine DNA molecular testing, based on the identification of D4Z4 arrays with less
than 10 units at 4q35, has been considered highly sensitive and specific [23,24]. However, several
genotype-phenotype studies have shown that D4Z4 reduction (4–8 reduced units, RU) is also present
in 3% of the general population [25–28] and in some cases is associated with distinct myopathic
phenotypes not reminiscent of FSHD [29,30]. Moreover, 10% of FSHD patients carry D4Z4 alleles
of borderline size (9–10 RU), which are also found in healthy people or in subjects with a different
myopathy [31] (Ricci et al., submitted). Finally, 5–10% of FSHD patients carry D4Z4 arrays of size
within the range of the general population (11 RU or more) on both chromosomes 4q [32,33]. These
subjects represent a second form of disease, FSHD2. Therefore, the number of D4Z4 RU at 4q35 does
not per se characterize FSHD (Figure 1).

As a matter of fact, genotype-phenotype studies have shown a large spectrum of clinical
phenotypes in myopathic subjects [29–31], carrying a D4Z4 reduced allele as well as reduced penetrance
among relatives carrying the same D4Z4 reduced allele [7,9,10,12,25,34,35]. All this has relevant effects
on clinical practice, complicating diagnosis, prognosis, and genetic counseling.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Molecular complexity in facioscapulohumeral muscular dystrophy (FSHD). D4Z4 contractions are associated with a permissive 4qA genotype that involves the aberrant expression of the *DUX4* retrogene and is responsible for FSHD1, but also occurs in 1.3–2% of the normal population. 4qA is also found in cases presenting complex phenotypes and as well as in 30–50% of healthy relatives. Mutations in *SMCHD1* or *DNMT3B* genes have been associated with FSHD2 and are responsible for ICF and BAMS syndromes. Parents of ICF patients are heterozygous for *DNMT3B* pathogenic variants but do not show any sign of muscular dystrophy. Other myopathic patients carry *SMCHD1* mutations. Grey dot = methyl CpG; white dot = unmethylated CpG.</figcaption></figure>

2. Molecular Features and the Epigenetic Model for FSHD

The D4Z4 repeat is part of a family of 3.3 kb sequences frequently found in heterochromatic regions, such as the short arms of the acrocentric chromosomes [36], and a nearly identical and equally polymorphic D4Z4 array reside at the sub-telomere of chromosome 10q [37,38].

D4Z4 is a CpG-rich (73%) macrosatellite DNA element encompassing more than 16 nucleosomes and containing multiple repeat sequences normally associated with heterochromatin [36]. The D4Z4 repeat unit contains the open reading frame of the retrogene *DUX4*.

Studies showed that D4Z4 plays a role in control of gene silencing at 4q35 through the recruitment
of a multi-protein repressor complex, the D4Z4 Recognition Complex (DRC), composed of YY1,
HMGB2, and Nucleolin, binding a 27bp DNA element within the D4Z4 sequence (D4Z4 Binding Element, DBE) [39]. The factor YY1 is a component of multiple chromatin regulatory complexes, including the Polycomb Repressive Complex 2 (PRC2), which includes the H3K27 methyltransferase EZH2 [40]. YY1 also interacts with histone deacetylase 1 (HDAC1), HDAC2, and HDAC3 [41], as well as with PARP1 [42]. A long non coding RNA, DBE-T, encoded within the 4q35 locus, is selectively transcribed in FSHD samples and participates in the transcriptional and epigenetic regulation of the 4q35 genes [43].

Several clinical features, such as penetrance variability, gender bias in severity [6], asymmetric
muscle wasting, and discordance in monozygotic twins [13–15], suggest that FSHD development
involves epigenetic factors which might influence gene expression through local modification of
chromatin structure.

Indeed, the great clinical heterogeneity and molecular uncertainty has pointed to the need for additional
markers to support FSHD diagnosis, genetic counseling, and patient stratification for clinical trials.

In the rarer cases of FSHD2 (OMIM#158901), more than 80% of subjects carry heterozygous mutations in the *SMCHD1* (structural maintenance of chromosomes flexible hinge domain containing 1) gene, which encodes a chromatin remodeling protein required for normal DNA methylation levels and transcriptional repression at certain loci, including the inactive X chromosome, imprinted genomic regions, and the D4Z4 arrays [44–46]. A few FSHD2 cases carry heterozygous mutations in the *DNMT3B* (DNA methyltransferase 3 beta) gene [47] which is responsible for the establishment of the proper de novo cytosine methylation profile during development. FSHD1 and, at a greater extent, FSHD2 patients present D4Z4 DNA hypomethylation [33,48,49].

All these findings supported the hypothesis that significant alteration of the 4q35 epigenetic
landscape is central to FSHD, and that this occurs in FSHD1 by removal of a significant number of D4Z4
heterochromatic elements, and in FSHD2, via SMCHD1 or DNMT3B haploinsufficiency in presence of
a specific telomeric polymorphism, 4qA, which allows the expression of the most distal copy of the
DUX4 gene. More specifically, the PAS present on 4qA alleles allows for DUX4 mRNA polyadenylation,
stabilization, and translation [50–52]. In turn, all these alterations have brought forward a unifying
model for FSHD pathogenesis, involving the loss of epigenetic silencing and the consequent aberrant
expression of the *DUX4* retrogene [52] (Figure 1). At present, these molecular characteristics are used
to explain FSHD pathogenesis and to develop therapeutic approaches.

3. The Debated Role of DNA Methylation in FSHD: Clinical and Families Studies

DNA methylation, a covalent post-synthetic modification of cytosines engaged in CpG dinucleotides, is a heritable epigenetic mark. In healthy individuals, 70% to 90% of the CpGs are methylated in somatic tissues, representing between 0.75% and 1% of the total number of bases in the diploid human genome [53]. In human cells, dispersed CpGs are methylated whereas CpG clusters (CpG islands) are mostly unmethylated [54]. Repetitive DNA sequences, enriched in CpGs, are usually packed as condensed and repressed chromatin by dense methylation, a protective mechanism that inhibits the invasion of the genome by the reactivation of transposable elements [55,56]. DNA methylation plays important roles in a number of physiological processes, such as development [57,58] and ageing [59–61]. Aberrant DNA methylation plays a causal role in a variety of diseases, including cancer [62,63]. In general, cancer cells exhibit DNA hypermethylation of promoter regions of tumor suppressor genes and global hypomethylation of repetitive DNA sequences accompanied with an increased genomic instability or loss of heterozygosity. Methylation changes also contribute to a number of diseases [64,65]. In particular, methylation have been investigated in diseases associated with short tandem repeats DNA [66] and CpG methylation has been considered a possible disease marker or modifier [67]. In fact, DNA hypermethylation is a critical feature of expanded short repeat arrays, as reported for Fragile X syndrome [66] or Friedreich's ataxia [68].

In FSHD, reduction of epigenetic silencing and aberrant gene expression seems to constitute the
underlying mechanism leading to disease. It has thus been thought that reduced methylation at the
D4Z4 locus might represent a proxy indicator of the reduced transcriptional silencing associated with FSHD. As a consequence, low CpG methylation of the D4Z4 sequence has been proposed as a reliable marker in the FSHD diagnostic procedure [33,48,49].

To investigate this hypothesis on a large-scale, different works have been conducted on distinct clinical cohorts. These studies assessed the methylation status at D4Z4, trying to correlate the percentage of CpG methylation at the D4Z4 array with one or more FSHD features [69–73].

In early works, DNA methylation has been investigated by Southern blotting after digestion with methylation-sensitive restriction enzymes (MSRE assay) [33,48] and hybridization with the p13E11 probe, which detects a unique region upstream of the first repeat of the D4Z4 array (Figure 2). The percentage of CpG methylation of the D4Z4 was estimated based on the density of the hybridization signals of restriction fragments obtained prior or after restriction digestion with MSRE. The most proximal *FseI* site within the D4Z4 array was considered as the most sensitive for this assay with an estimated DNA methylation of 70–80% in healthy individuals and a significant decrease in FSHD1 patients. The MSRE assay possesses a restricted range of investigation since it evaluates only the few CpGs detected by the various of MSREs without distinction between the D4Z4 repeat array on the 4q and 10q chromosome [69], thus offering a narrow view of the D4Z4 methylation status.

To overcome these limitations, the bisulfite treatment of DNA followed by DNA sequencing (BSS) has become the technique of choice to assess the D4Z4 methylation, virtually allowing the detection of the methylation status of any CpGs within the array at once.

With this technical improvement, many works succeeded in reporting a more exhaustive analysis of many CpGs within all D4Z4 units, even if without distinction between the 4q and 10q alleles [70,71]. Recently, Jones et al. [69] developed a 4q- and 10q-specific protocol for BSS analyses, which is concomitantly able to capture the epigenetic status of the 3' end of the most distal D4Z4 repeat and the abutting A-type locus, which has been correlated with the pathogenic expression of *DUX4* on short D4Z4 arrays. A further advance in the methodology used in methylation studies was obtained by Roche et al. [74] consisting in the assessment of DNA methylation at D4Z4 by using a custom deep sequencing method after sodium bisulfite conversion of genomic DNA.

Nevertheless, the different attempts to use BSS analysis to assess D4Z4 epigenetic status suffer from several limitations, both from a conceptual and a technical point of view.

As shown in Table 1, we compared the results of the main reports in the field. The comprehensive analysis of these studies highlights some problematic points: (1) the lack of a clear description of the clinical status of individuals within each cohort [69–72,74]; (2) the ambiguous definition of asymptomatic or healthy carriers of a D4Z4 allele of reduced size, since they are included either in the group of people with FSHD even if their clinical score is equal to zero, or in the control group despite their liability to develop the disease [69,72]; (3) the constant attempt to correlate the methylation status with the disease severity [70,72]; and (4) the effort to uniquely correlate hypomethylation at D4Z4 with *DUX4* expression [69,70,72].

From a technical point of view, the scenario is even more complicated by the fact that each work used different primers sets in the BSS analyses to evaluate CpGs within different regions of the D4Z4 array (Figure 2) and applied different statistical tools. First, for a proper analysis of DNA methylation by bisulfite sequencing, (1) primers should not contain CpG sites within their sequence to ensure unbiased amplification of both methylated and unmethylated DNA; and (2) primers should be designed in a region containing enough number of non-CpG cytosines to avoid amplification of incompletely modified genomic DNA [75]. Despite these recommendations in some of the works considered here, the primer designed for BSS analysis included CpGs [69] or just one or few non-CpG cytosines [73], which might affect the amplification of the converted DNA. Secondly, the selection of different regions within the array led to controversial results about which is/are the region/s more representative for the methylation status at D4Z4.
<table><caption>**Table 1.** Systematic analysis and comparison of the major literature reports on D4Z4 methylation in FSHD.</caption><thead><tr><th>REFs</th><th>Method</th><th>Assayed Region (s)</th><th>4q-Specific</th><th>Nr. of Subjects</th><th>Clinical Status (CCEF/ACSS)</th><th>Hypomethylation in FSHD1 /FSHD2</th><th>Position within the Locus in which Hypomethylation was Found to be Discriminant</th><th>Correlation with D4Z4 Size</th><th>Correlation with Disease Severity</th><th>Matters</th></tr></thead><tbody><tr><td>[69]</td><td>BSS</td><td>4qA; 56 CpGs or 4qAL; 30 CpGs DUX4 5'; 59CpGs</td><td>Yes No</td><td>Families 12 Controls 10 Asymptomatic 22 FSHD1</td><td>Not reoprted</td><td>Yes; Asymptomatyc carriers with intermediate methylation levels.</td><td>4qA: pLAM only No DUX4 5'</td><td>No</td><td>Not tested</td><td>.Different results using 5' primers set. .Small differences between pathogenic contracted) and non-pathogenic allele. .No correlation with allele- specific % of DUX4 expression. .Unclear interpretation of data due to a global estimation of methylation as a function of D4Z4 repeat lengths. .Non-penetrant carriers of a DRA included in analysis.</td></tr><tr><td>[70]</td><td>MRSE1</td><td>FseI site</td><td>No</td><td>254 controls ? 25 Asymptomatic 186 FSHD1 individuals 74 FSHD2-(just SMCHD1 mut) 20 Controls 29</td><td>ASCC</td><td>Yes; Asymptomatyc carriers not clearly defined</td><td>FseI site</td><td>Yes for 1-6 repeats carriers No for 7-10 repeats carriers</td><td>No for FSHD1; Yes for FSHD2</td><td></td></tr><tr><td>[71]</td><td>meDIP BSS</td><td>D4Z4-5P; 22 CpGs D4Z4-MID; 31 CpGs D4Z4-3P; 14 CpGs</td><td>No</td><td>Asymptomatic 37 FSHD1 9 FSHD2 8 FSHD1 4FSHD2 1Asymptomatic 7 Controls 51 Controls 2Asymptomatic 44FSHD1 17 FSHD2 (just SMCHD1 mut)</td><td>Not reported</td><td>Yes; Asymptomatyc carriers with methylation levels not different than controls</td><td>D4Z4-5' only No MID or 3'</td><td>No</td><td>No</td><td>.Limited BSS analysis. No correlation with the global number of residual repeats.</td></tr><tr><td>[72]</td><td>BSS MRSE1</td><td>DUX4-PAS FseI site</td><td>Yes No</td><td></td><td>ASCC</td><td>Yes; Asymptomatyc carriers not clearly defined No</td><td>DUX4-PAS: pLAM only No FseI</td><td>Yes</td><td>Yes</td><td>.Arbitrary selection of just one (number 6) CpG for statistical analysis. .Non-penetrant carriers of a DRA included in analysis.</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>REFs</th><th>Method</th><th>Assayed Region (s)</th><th>4q-Specific</th><th>Nr. of Subjects</th><th>Clinical Status (CCEF/ACSS)</th><th>Hypomethylation in FSHD1 /FSHD2</th><th>Position within the Locus in which Hypomethylation was Found to be Discriminant</th><th>Correlation with D4Z4 Size</th><th>Correlation with Disease Severity</th><th>Matters</th></tr></thead><tbody><tr><td>[73]</td><td>MRSE1+ MRSE2 BSS+ MRSE1+ MRSE2</td><td>FseI site BsaAI site 4qA; 56 CpGs or 4qAL; 30 CpGs</td><td>No Yes Yes</td><td>122P+110: 88 FSHD1 P + 47FSHD1R 34 FSHD2P+ 2 FSHD 2 R 61 Asymptomatic R 1 Control 11FSHD1 8 Asymptomatic</td><td>CCEF</td><td>Assumed: beyond the purpose of the paper. Asymptomatic carriers with methylation levels not different than index cases</td><td>FseI site 4qA: pLAM</td><td>No</td><td>No</td><td>-Limited BSS analysis</td></tr><tr><td>[74]</td><td>BSS +NGS</td><td>D4Z4-5P; 22 CpGs D4Z4-MID; 31 CpGs D4Z4-3P; 14CpGs DR1; 31CpGs 4aAint;14 CpGs 4qAext;42 CpGs</td><td>No Yes</td><td>10 Controls 29 FSHD1 10 FSHD2-(just SMCHD1 mut)</td><td>Not reported</td><td>Yes Asymptomatic carriers not included in analysis</td><td>D4Z4-5' only</td><td>No</td><td>No</td><td>.Non-penetrant carriers of a DRA not included in analysis.</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Overview of the D4Z4 repeats and of the region selected in methylation analyses. The D4Z4 array at 4q35 with an enlarged schematic representation of the D4Z4 distal most repeat (from position 1 to 3303 given relative to the two flanking *Kpn*I sites). The different regions within D4Z4 are indicated: LSau repeat (position 1–340), Region A (position 869–1071), *hhspm3* (position 1313–1780), *DUX4 ORF* (position 1792–3063), plus the 3' pLAM region. The figure highlights the proximal *BsaAI* and *FseI* methylation-sensitive restriction sites analyzed by Southern blotting, and the location of bisulfite (BS) PCR products used in the selected literature reports, each represented with a differently colored line. The position of each region within the array, starting from the first *Kpn*I site, is indicated together with the assayed number of CpGs. The white lines within these regions indicate the presence of CpGs in PCR primers.</figcaption></figure>

As anticipated, we also observed that statistical analyses performed in each individual study could not be compared with each other. For instance, in the work by Calandra et al. [72], the assessment of the correlation between methylation levels and disease was performed focusing on a single CpG out of ten within a region where D4Z4 hypomethylation is more evident.

In the study by Lemmers et al. [70] the global methylation levels were estimated by means of a novel statistical model considering the D4Z4 array as a linear string of mathematic units: the methylation level was calculated as a linear function of D4Z4 repeats number. In biological terms, this arbitrary approach has an intrinsic bias, since the repeat length of each allele affects the 3D chromatin structure of the individual D4Z4 locus, and most likely its methylation status. Longer alleles generate a condensed chromatin structure and beyond a certain number of repetitive elements, the D4Z4 region displays a high percentage of CpG methylation irrespective of the discrete D4Z4 repeat number of each individual long allele.

Moreover, when considering alleles with fewer D4Z4 repeats, the effect on DNA methylation differs in each D4Z4 allele and depends on the individual size. In fact, in families in which it was possible to follow the segregation of D4Z4 alleles analyzed by bisulfite sequencing, it was possible to discriminate between the two D4Z4 alleles carrying different repeat numbers. In these cases, the degree of methylation roughly correlates with number of alleles with 10 or fewer repeats, but it was not associated with the clinical condition [73].

Regardless of the methodologies used, data indicates that D4Z4 methylation is highly variable also in healthy individuals and that hypomethylation concerns only a limited number of CpGs within
the D4Z4 sequence. Hypomethylated CpGs are mostly clustered in the proximal part of the repeat but not in the *DUX4* promoter region per se or in the distal part of the coding sequence and 3' UTR [69].

Overall, the observed hypomethylation of D4Z4 in FSHD patients shows a correlation with the reduction of D4Z4 repeat units, even though the accurate review of the published data reveals a lack of association between D4Z4 methylation level and FSHD1 clinical status.

In conclusion, a detailed phenotypic characterization should be recommended, together with a common study design and a consensus approach in order to reach a more accurate and unbiased evaluation of the global methylation pattern and to draw any conclusion on its clinical significance.

## 4. DNA Methylation in FSHD2

The SMCHD1 chromatin-associated factor has been implicated in FSHD on the basis of Whole Exome Sequencing of 16 FSHD2 families presenting residual D4Z4 methylation below a threshold of 25% at the most proximal *FseI* site [76]. It is now recognized that approximately 80% of FSHD2 patients are carriers of one *SMCHD1* variant [77]. In a subsequent study, the threshold that defines hypomethylation was estimated to be 30% of residual methylation at this proximal *FseI* site [70]. However, there seems to be no strict correlation between *SMCHD1* variants, D4Z4 hypomethylation, and the appearance of clinical signs. This is exemplified by the non-affected individuals who display a low DNA methylation level <30%; 26% for Rf854-case 2398; 30% for Rf854-case 2434; and 29% for Rf300-Cases 3 and 4, presented in the original work of Lemmers et al. [76]. More strikingly, considering the clinical phenotype, the marked D4Z4 hypomethylation detected in the majority of FSHD2 patients carrying a mutation in *SMCHD1* is not accompanied either with a more severe phenotype or earlier disease onset. Moreover, D4Z4 hypomethylation, in the presence of the haplotype 4qA/PAS distal to the D4Z4 array, which is considered permissive of *DUX4* expression, is not associated with a muscle phenotype in patients with Bosma Arhinia and Microphthalmia syndrome (BAMS), patients affected with ICF (Immunodeficiency, Centromeric Instability, and Facial anomalies) homozygotes for *DNMT3B* mutation and their heterozygote parents [78], or patients carrying a deletion of the 18p chromosome containing the *SMCHD1* gene [79]. Besides, in the few families with FSHD index cases carrying a mutation in *DNMT3B* [47] or an 18p deletion [80], the presence of a DNA mutation and D4Z4 hypomethylation does not segregate with clinical signs of the disease. Furthermore, as highlighted by molecular combing [81], a significant proportion of individuals affected by FSHD display an atypical genotype, with complex distal rearrangements [82], presence of proximal deletions, additional 10q alleles, in the absence of D4Z4 array shortening, *SMCHD1* variants, or D4Z4 hypomethylation [83].

So far, further investigations are required to answer a few important questions, such as the role of epigenetic alterations in FSHD, whether hypomethylation leads to *DUX4* activation and, above all, how this translates into a specific muscle phenotype.

## 5. Changes in DNA Methylation at D4Z4 Upon Reprogramming

Pluripotent stem cells represent an important tool to model human genetic disorders, as they can be used to analyze the effect of specific genomic alterations in early development or during differentiation [83–86]. Cell reprogramming always involves profound epigenetic changes and the acquisition of an epigenetic pattern similar to that of pluripotent embryonic stem cells. The pattern of DNA methylation is not fully erased after reprogramming. Hence, reprogrammed cells keep memory of the pattern of the parental somatic cells they were derived from, but also acquire DNA methylation profiles specific of induced pluripotent stem cells (iPSCs). iPSCs obtained by reprogramming of patients’ cells have been recently used to clarify the molecular basis of several “epigenetic” diseases. For instance, in human hiPSCs from Fragile X (FXS) patients, the *FMR1* promoter region aside the 5’UTR with expanded CGG triplets (n > 200) remains in most cases hypermethylated, suggesting that once established in FXS fibroblasts, these epigenetic marks are stably maintained after reprogramming, whereas in other cases, the methylation pattern of the disease-associated locus does not reflect the profile of the donor cells [87–90].
Concerning FSHD1, the investigation of DNA methylation at D4Z4 in iPSCs from affected subjects and healthy donors revealed that these cells do not retain the methylation pattern inherited from the donor cells, but acquire a new methylation profile. This argues for the presence of two mechanisms acting at D4Z4 upon reprogramming: an active "erasure" of the cell-of-origin epigenetic profile and a "rewriting" of a de novo methylation pattern at the array. Moreover, reports demonstrate that the D4Z4 methylation was identical from clone to clone from both FSHD1 patients and controls, indicating that remethylation of D4Z4 upon epigenetic reprogramming does not depend on the residual number of D4Z4 repeats [79].

More in details, the comparison of the methylation profile between FSHD iPSCs and human Embryonic Stem Cells (hESCs) with a short or long D4Z4 array, revealed a trend toward a high methylation level, without any significant difference between FSHD1 and control cells and between iPSCs and ESCs.

The higher methylation level found at D4Z4 in pluripotent cells further suggests that D4Z4 methylation status does not correlate with the number of repeats but is a feature of stemness, which highlights once more the complex but yet unknown regulatory mechanisms of this locus [79].

Interestingly, recent findings revealed a key role for *DUX4* at very early stages of human development and the activation of embryonic genes at the zygotic genome activation stage (ZGA) [91–94], a stage characterized by profound methylation changes. A high level of methylation was observed in hiPSCs and hESCs. These cells represent a later developmental stage compared to ZGA since they are derived from the blastocyst inner cell mass (ICM), a developmental stage characterized by remethylation of the zygotic genome after post-fertilization waves of demethylation up to the morula stage. These findings rule out the connection between *DUX4* expression and D4Z4 hypomethylation.

Additional features, such as the size of the abutting telomere, might be also implicated in the regulation of D4Z4 methylation. Interestingly, the highest level of methylation is observed in cells where telomerase is reactivated, such as hiPSCs and ES cells [95–97]. These observations suggest the possible link between D4Z4 methylation and telomeres and also draw attention to the complexity of the epigenetic regulation of this macrosatellite element during development.

## 6. Trans-Acting Factors

The possibility to use the percentage values of methylated CpGs at D4Z4 as a diagnostic marker for FSHD is further weakened by the fact that DNA methylation could not be considered per se. It is now well established that DNA methylation is only one among numerous indicators of chromatin structure at D4Z4 and several research groups reported that alterations of chromatin structure is reflected by post-translational histone modifications [48,98–100], by the binding of various non-histone proteins and RNAs on the repeat array [39,43,101], and by higher order chromatin structures formation [102–104].

Epigenetic alteration of chromatin relies on DNA methylation and on histone tails modifications in nucleosomes, which together determine the so called "epigenetic signature" of a specific region. It has been reported that the D4Z4 array contains histone H3 Lysine 9 trimethylation (H3K9me3), a repressive mark associated with heterochromatin formation, together with H3K27me3, a repressive chromatin mark associated with Polycomb silencing [33,43,98]. A specific loss of H3K9me3 followed by the loss of HP1γ and cohesin binding at D4Z4, suggestive of a more relaxed chromatin structure [98], was reported both for FSHD1, in the presence of a reduced number of repeats, and also in FSHD2 subjects. The loss of repressive histone marks is not a direct consequence of DNA hypomethylation since H3K9me3 enrichment at D4Z4 was unaltered in ICF syndrome, which displays severe DNA hypomethylation at the 4q35 locus [49,98]. Besides, H3K27 presence at D4Z4 did not show any changes between the control and FSHD muscle cells [43,102].

Other observations argue for a role for the Polycomb proteins in the epigenetic regulation at D4Z4: The methyl-transferase EZH2 (Enhancer of zeste homolog 2), a component of Polycomb repressive complex 2 (PRC2), was shown to be reduced in FSHD muscle cells, while ASH1L (ASH1 Like Histone Lysine Methyltransferase), a member of the Trithorax complex associated with transcriptionally
active chromatin, was found to be specifically recruited to D4Z4 in FSHD cells [43]. Notably, ASH1L recruitment has been shown to be dependent on the expression of the non-coding RNA DBE-T, transcribed by to the most proximal D4Z4 unit. In particular, DBE-T-mediated ASH1L recruitment to D4Z4 in FSHD cells is reported to mediate histone H3 lysine 36 dimethylation (H3K36me2), a major histone mark associated with transcriptional activation, leading to chromatin remodeling and 4q35 gene transcription [43].

Among transcription factors, two proteins have mainly been associated with D4Z4, YY1 Yin Yang 1) and CTCF (CCCTC-binding factor). Reports demonstrated the association of a repressor protein complex, named DRC (D4Z4 Recognition Complex) composed of YY1, a known transcriptional repressor, HMGB2 (high-mobility group protein 2), an architectural protein, and nucleolin, with the 27 bp DBE (D4Z4 binding element) sequence contained within each D4Z4 [39]. YY1 was reported to be the major factor bound to the DBE, although other DRC components contributed to DNA binding affinity and specificity. It is noteworthy that YY1 binds DNA only when it is not methylated [105]. The region encompassing the DBE element is hypermethylated in FSHD and control subjects, except from the CpGs neighboring the YY1 consensus site, which show variable levels of methylation (Tupler, personal observation and [106]). Conversely, CTCF binding is disrupted by CpG methylation [107]. Two CTCF binding sites are located within the 5' region of D4Z4 repeat, which is reported to display hypomethylation in the presence of a DRA [74,101,106,108]. CTCF binding to 4q35 displays an inverse correlation to D4Z4 copy number and function as a Lamin A/C-dependent chromatin insulator, protecting D4Z4 from epigenetic silencing by surrounding heterochromatic regions, therefore keeping the D4Z4 chromatin open primarily in a DRA context [102]. Since CTCF has also been shown to mediate chromatin loop formation and to generate TADs (topologically associated domain) [108–111] the increased CTCF binding to D4Z4 in FSHD may result in altered nuclear and chromatin organization. In particular, using an integrated genome wide approach (4C-seq) it has been reported that 4q-D4Z4 interactome is altered in FSHD1, leading to a chromatin switch toward an active state (mediated by enhancer-promoter interactions), which in turn results in the transcriptional activation of genes involved in muscular atrophy [103].

All these observations further complicate the comprehension of D4Z4 activity regulation and confirm that DNA methylation and its consequence on DUX4 activation are only a part of the fine tuning controlling the transition between healthy and disease status.

## 7. The Significance of D4Z4 Hypomethylation in FSHD and Its Implication for Clinical Counseling

The main current theory about the FSHD pathogenic model involves the aberrant expression of *DUX4* retrogene caused by epigenetic changes in the D4Z4 region [93]. This should be due to reduced methylation, either for D4Z4 repeats reduction in FSHD1 or for the effect of the *SMCHD1* gene in FSHD2 [52,76]. However, the literature reports clearly reveals that hypomethylation correlates with the reduction of the arrayed repeats, but not with the disease status [69–72]. It is thus conceivable that in FSHD1, hypomethylation is a consequence of the chromatin structure present in the contracted allele, instead of an indicator of its function.

This also holds true for FSHD2, especially for patients who carry mutations in *SMCHD1*, deletions of the 18p locus or mutations in *DNMT3B*, and especially considering the molecular overlaps (D4Z4 hypomethylation) with other unrelated diseases, such as BAMS, ICF, or 18p deletion syndromes (Figure 3).

All these observations indicate that accurate diagnosis should not completely rely on molecular findings. Large studies based on clinical evaluation demonstrate that (1) 2% of healthy individuals from the general population bear one FSHD1-sized D4Z4 allele with the possibility of expressing the DUX4 transcript [27]; (2) the clinical severity of FSHD does not clearly associate with the size of the deletion of the D4Z4 allele [7,12,112]; (3) the penetrance of FSHD is not complete, with 20–50% of relatives carrying FSHD1-sized D4Z4 alleles being healthy [7,9,10,25,34,35]; (4) the disease penetrance decreases among
second-fifth degree relatives in FSHD families [7]; (5) there are families showing affected subjects in one generation, families with a single affected subject, and also subjects with atypical phenotypes in which additional mutations in genes causing other neuromuscular diseases have been found [7,30,34]; in addition, independent studies showed that (6) *SMCHD1* mutations cause the Bosma developmental syndrome (OMIM #603457) and D4Z4 hypomethylation has been found in these patients with no documented signs of FSHD [113,114]; (7) *SMCHD1* mutations have been found in subjects with Limb Girdle Muscular Dystrophies (LGMDs) (Tupler, unpublished NGS data); (8) homozygous or compound heterozygous mutations in *DNMT3B* cause immunodeficiency-centromeric instability-facial anomalies syndrome-1 (ICF1) (OMIM #242860), a syndrome characterized by hypo- or agammaglobulinemia and frequent infections that determine poor life expectancy; and (9) some ICF1 cases showed D4Z4 hypomethylation with no signs of FSHD [78] (Figure 3).

Thus, *SMCHD1* and *DNMT3B* mutations have been associated with congenital diseases, whereas FSHD is an adult progressive disease with no congenital presentation, and very importantly, no FSHD patients carrying *SMCHD1* mutations have signs of Bosma syndrome even when they carry overlapping missense mutations or mutations in the same regions of the coding sequence. These data are consistent with the frequent observation that genes play roles in multiple pathways, and that different mutations within the same gene can contribute to distinct phenotypes with different mechanisms of action (Figure 3).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** The lack of a unique molecular signature in FSHD complicates genotype-phenotype correlation in clinical practice. The figure illustrates the complex molecular scenario observed in individuals with clinical features of FSHD. Hypomethylation is observed in a large subset of them, including in patients who carry a mutation in *SMCHD1* or *DNMT3B*. This scheme also highlights that hypomethylation of D4Z4 is observed in other rare diseases linked to mutations in *SMCHD1* (BAMS), *DNMT3B* (ICF) or carrying a deletion of the 18p locus that contains the *SMCHD1* gene.</figcaption></figure>

**8. Conclusions**

Unraveling the mechanism leading to FSHD reveals to be a prickly issue. One of the major challenges we face in clinical practice is that there is no unique molecular signature discriminating FSHD patients. We are still trying to promote a linear common diagnostic approach to FSHD, an ambitious objective that should involve all clinical practitioners in the field, as displayed by the use of the FSHD Comprehensive Clinical Evaluation Form [29,115]. At present, current available molecular markers are not sufficient for a proper diagnosis or prognosis, nor for genetic counseling, but patients are still looking for support and assistance on a daily basis. Because of the large phenotypic variations observed among affected subjects, the clinical assessment should be performed carefully considering the family context, with the attempt to find a pattern of inheritance and shared features
Our hypothesis is that FSHD clinical variability may result from the combinatorial effects of variant alleles in genes exerting a detrimental effect on muscle function, together with epigenetic changes influencing their expression. The genetic background of each individual, including harmful gene variants or variants in regulatory elements, might have a crucial role in disease penetrance, high clinical variability and unpredictable progression. In view of the results of genotype-phenotype studies on FSHD cohorts and of the recent findings on the effects of *SMCHD1* mutations, it is advisable to apply a stringent analysis of clinical phenotypes to decipher molecular data. An oversimplification may cause misdiagnosis and biased interpretation of molecular findings with severe effects not only for research, following incorrect hints, but also for clinical practice. The systematic collection of precise phenotypic data of patients and families should become an ineludible tool to support research in human molecular genetics and to translate molecular findings into clinical practice.

Besides, the D4Z4 array and the regulation of its chromatin dynamics has become more and more intriguing over the recent years because of the identification of numerous factors able to regulate its chromatin state and topology within the genome. This element has often been considered separately and not necessarily in the context of the 4q35 locus itself, which also contains several microsatellites of unknown function upstream of D4Z4, and repetitive beta satellite elements distal to the D4Z4, not mentioning the nearby telomere. Variations in macrosatellite number and associated features suggest that not only D4Z4, but also other elements of the large family of repetitive sequences, might play important regulatory functions either locally or at long distances on gene expression. It is speculated that copy number variations of repetitive elements might be associated with variations in human phenotype and they could explain different susceptibility to the disease and missing heritability. In FSHD, the molecular features of the 4q35 sub-telomeric region could be an example of this molecular and clinical variegation.

In conclusion, the study of CpG methylation, performed by several groups and in distinct ways, failed to reveal a clear and unifying theory about the "epigenetic" basis of FSHD. Indeed, the D4Z4 methylation status does not mirror the clinical expression of the disease. The measurement of this epigenetic mark must be interpreted with caution and may support the characterization of individual FSHD families in clinical practice rather than be assumed as a marker of D4Z4 function in FSHD.

**Author Contributions:** The concept and idea and of this review was done by R.T., V.S. and F.M.; all authors wrote the draft of the manuscript; R.T. revised the entire text; and All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by Telethon Italy GUP13012, by Association Francaise contre les Myopathies, TRIM-RD, by Regione Emilia Romagna progetto RARER.

**Acknowledgments:** We thank Antonio Vallarola for the in-depth critique of the manuscript and Maria Francesca Di Feo for editing the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Deenen, J.C.W.; Arnts, H.; Van Der Maarel, S.M.; Padberg, G.W.; Verschuuren, J.J.G.M.; Bakker, E.; Weinreich, S.S.; Verbeek, A.L.M.; Van Engelen, B.G.M. EPopulation-based incidence and prevalence of facioscapulohumeral dystrophy. *Neurology* **2014**, *83*, 1056–1059. [CrossRef] [PubMed]

2. Mostacciuolo, M.L.; Pastorello, E.; Vazza, G.; Miorin, M.; Angelini, C.; Tomelleri, G.; Galluzzi, G.; Trevisan, C.P. Facioscapulohumeral muscular dystrophy: Epidemiological and molecular study in a north-east Italian population sample. *Clin. Genet.* **2009**, *75*, 550–555. [CrossRef]

3. Padberg, G.W.A.M.; LUMC. Facioscapulohumeral Disease. Ph.D. Thesis, Leiden University, Leiden, The Netherlands, 1982.

4. Flanigan, K.M. *Facioscapulohumeral Muscular Dystrophy and Scapulohumeral Syndrome*; McGraw-Hill: New York, NY, USA, 2004.

5. Zatz, M.; Marie, S.K.; Passos-Bueno, M.R.; Vainzof, M.; Campiotto, S.; Cerqueira, A.; Wijmenga, C.; Padberg, G.; Frants, R. High proportion of new mutations and possible anticipation in Brazilian facioscapulohumeral muscular dystrophy families. *Am. J. Hum. Genet.* **1995**, *56*, 99–105. [PubMed]
6. Tonini, M.M.O.; Passos-Bueno, M.R.; Cerqueira, A.; Matioli, S.R.; Pavanello, R.; Zatz, M. Asymptomatic carriers and gender differences in facioscapulohumeral muscular dystrophy (FSHD). *Neuromuscul. Disord.* **2004**, *14*, 33–38. [CrossRef] [PubMed]

7. Ricci, G.; Scionti, I.; Sera, F.; Govi, M.; D'Amico, R.; Frambolli, I.; Mele, F.; Filosto, M.; Vercelli, L.; Ruggiero, L.; et al. Large scale genotype–phenotype analyses indicate that novel prognostic tools are required for families with facioscapulohumeral muscular dystrophy. *Brain* **2013**, *136*, 3408–3417. [CrossRef]

8. Park, H.J.; Hong, J.M.; Lee, J.H.; Lee, H.S.; Shin, H.Y.; Kim, S.M.; Ki, C.S.; Lee, J.H.; Choi, Y.C. Low D4Z4 copy number and gender difference in Korean patients with facioscapulohumeral muscular dystrophy type 1. *Neuromuscul. Disord.* **2015**, *25*, 859–864. [CrossRef]

9. Wu, Z.Y.; Wang, Z.Q.; Murong, S.X.; Wang, N. FSHD in Chinese population: Characteristics of translocation and genotype-phenotype correlation. *Neurology* **2004**, *63*, 581–583. [CrossRef]

10. Lin, F.; Wang, Z.-Q.; Lin, M.-T.; Murong, S.-X.; Wang, N. New Insights into Genotype-phenotype Correlations in Chinese Facioscapulohumeral Muscular Dystrophy. *Chin. Med. J. (Engl)*. **2015**, *128*, 1707–1713. [CrossRef]

11. Tawil, R.; Kissel, J.T.; Heatwole, C.; Pandya, S.; Gronseth, G.; Benatar, M.; Guideline Development, Dissemination, and Implementation Subcommittee of the American Academy of Neurology; Practice Issues Review Panel of the American Association of Neuromuscular & Electrodiagnostic Medicine. Evidence-based guideline summary: Evaluation, diagnosis, and management of facioscapulohumeral muscular dystrophy. *Neurology* **2015**, *85*, 357–364. [CrossRef]

12. Goto, K.; Nishino, I.; Hayashi, Y.K. Very low penetrance in 85 Japanese families with facioscapulohumeral muscular dystrophy 1A. *J. Med. Genet.* **2004**, *41*, e12. [CrossRef]

13. Tawil, R.; Storvick, D.; Feasby, T.E.; Weiffenbach, B.; Griggs, R.C. Extreme variability of expression in monozygotic twins with FSH muscular dystrophy. *Neurology* **1993**, *43*, 345–348. [CrossRef] [PubMed]

14. Griggs, R.C.; Tawil, R.; McDermott, M.; Forrester, J.; Figlewicz, D.; Weiffenbach, B. Monozygotic twins with facioscapulohumeral dystrophy (FSHD): Implications for genotype/phenotype correlation. *Muscle Nerve* **1995**, *18*, S50–S55. [CrossRef]

15. Tupler, R.; Barbierato, L.; Memmi, M.; Sewry, C.A.; De Grandis, D.; Maraschio, P.; Tiepolo, L.; Ferlini, A. Identical de novo mutation at the D4F104S1 locus in monozygotic male twins affected by facioscapulohumeral muscular dystrophy (FSHD) with different clinical expression. *J. Med. Genet.* **1998**, *35*, 778–783. [CrossRef] [PubMed]

16. Wijmenga, C.; Frants, R.R.; Brouwer, O.F.; van der Klift, H.M.; Khan, P.M.; Padberg, G.W. Facioscapulohumeral muscular dystrophy gene in Dutch families is not linked to markers for familial adenomatous polyposis on the long arm of chromosome 5. *J. Neurol. Sci.* **1990**, *95*, 225–229. [CrossRef]

17. Padberg, G.W.; Lunt, P.W.; Koch, M.; Fardeau, M. Diagnostic criteria for facioscapulohumeral muscular dystrophy. *Neuromuscul. Disord.* **1991**, *1*, 231–234. [CrossRef]

18. Weiffenbach, B.; Bagley, R.G.; Falls, K.; Dubois, J.; Hyser, C.; Storvick, D.; Schultz, P.; Mendell, J.R.; Milner, E.C.; Jacobsen, S.J. Framework multipoint map of the long arm of human chromosome 4 and telomeric localization of the gene for FSHD. *Mamm. Genome* **1992**, *3*, 143–150. [CrossRef]

19. Sarfarazi, M.; Wijmenga, C.; Upadhyaya, M.; Weiffenbach, B.; Hyser, C.; Mathews, K.; Murray, J.; Gilbert, J.; Pericak-Vance, M.; Lunt, P. Regional mapping of facioscapulohumeral muscular dystrophy gene on 4q35: Combined analysis of an international consortium. *Am. J. Hum. Genet.* **1992**, *51*, 396–403.

20. Wijmenga, C.; Hewitt, J.E.; Sandkuijl, L.A.; Clark, L.N.; Wright, T.J.; Dauwerse, H.G.; Gruter, A.-M.; Hofker, M.H.; Moerer, P.; Williamson, R.; et al. Chromosome 4q DNA rearrangements associated with facioscapulohumeral muscular dystrophy. *Nat. Genet.* **1992**, *2*, 26–30. [CrossRef]

21. Tremblay, D.C.; Alexander, G.; Moseley, S.; Chadwick, B.P. Expression, tandem repeat copy number variation and stability of four macrosatellite arrays in the human genome. *BMC Genom.* **2010**, *11*, 632. [CrossRef]

22. Lunt, P.W.; Jardine, P.E.; Koch, M.; Maynard, J.; Osborn, M.; Williams, M.; Harper, P.S.; Upadhyaya, M. Phenotypic-genotypic correlation will assist genetic counseling in 4q35-facioscapulohumeral muscular dystrophy. *Muscle Nerve. Suppl.* **1995**, *2*, S103–S109. [CrossRef]

23. Upadhyaya, M.; Maynard, J.; Rogers, M.T.; Lunt, P.W.; Jardine, P.; Ravine, D.; Harper, P.S. Improved molecular diagnosis of facioscapulohumeral muscular dystrophy (FSHD): Validation of the differential double digestion for FSHD. *J. Med. Genet.* **1997**, *34*, 476–479. [CrossRef] [PubMed]
24. Tawil, R.; Van Der Maarel, S.M. Facioscapulohumeral muscular dystrophy. *Muscle Nerve* **2006**, *34*, 1-15.
[CrossRef] [PubMed]

25. NAKAGAWA, M.; MATSUZAKI, T.; HIGUCHI, I.; FUKUNAGA, H.; INUI, T.; NAGAMITSU, S.; YAMADA, H.; ARIMURA, K.; OSAME, M. Facioscapulohumeral Muscular Dystrophy: Clinical Diversity and Genetic Abnormalities in Japanese Patients. *Intern. Med.* **1997**, *36*, 333-339. [CrossRef] [PubMed]

26. Wohlgemuth, M.; Lemmers, R.J.; van der Kooi, E.L.; van der Wielen, M.J.; van Overveld, P.G.; Dauwerse, H.; Bakker, E.; Frants, R.R.; Padberg, G.W.; van der Maarel, S.M. Possible phenotypic dosage effect in patients compound heterozygous for FSHD-sized 4q35 alleles. *Neurology* **2003**, *61*, 909-913. [CrossRef] [PubMed]

27. Scionti, I.; Greco, F.; Ricci, G.; Govi, M.; Arashiro, P.; Vercelli, L.; Berardinelli, A.; Angelini, C.; Antonini, G.; Cao, M.; et al. Large-scale population analysis challenges the current criteria for the molecular diagnosis of fascioscapulohumeral muscular dystrophy. *Am. J. Hum. Genet.* **2012**, *90*, 628-635. [CrossRef] [PubMed]

28. Overveld, P.G.M.v.; Lemmers, R.J.; Deidda, G.; Sandkuijl, L.; Padberg, G.W.; Frants, R.R.; van der Maarel, S.M. Interchromosomal repeat array interactions between chromosomes 4 and 10: A model for subtelomeric plasticity. *Hum. Mol. Genet.* **2000**, *9*, 2879-2884. [CrossRef]

29. Ricci, G.; Ruggiero, L.; Vercelli, L.; Sera, F.; Nikolic, A.; Govi, M.; Mele, F.; Daolio, J.; Angelini, C.; Antonini, G.; et al. A novel clinical tool to classify facioscapulohumeral muscular dystrophy phenotypes. *J. Neurol.* **2016**, *263*, 1204-1214. [CrossRef]

30. Ricci, G.; Zatz, M.; Tupler, R. Facioscapulohumeral Muscular Dystrophy: More Complex than it Appears. *Curr. Mol. Med.* **2014**, *14*, 1052-1068. [CrossRef]

31. Butz, M.; Koch, M.C.; Muller-Felber, W.; Lemmers, R.J.L.F.; van der Maarel, S.M.; Schreiber, H. Facioscapulohumeral muscular dystrophy. *J. Neurol.* **2003**, *250*, 932-937. [CrossRef]

32. de Greef, J.C.; Lemmers, R.J.L.F.; Camano, P.; Day, J.W.; Sacconi, S.; Dunand, M.; van Engelen, B.G.M.; Kiuru-Enari, S.; Padberg, G.W.; Rosa, A.L.; et al. Clinical features of facioscapulohumeral muscular dystrophy 2. *Neurology* **2010**, *75*, 1548-1554. [CrossRef]

33. van Overveld, P.G.M.; Lemmers, R.J.F.L.; Sandkuijl, L.A.; Enthoven, L.; Winokur, S.T.; Bakels, F.; Padberg, G.W.; van Ommen, G.-J.B.; Frants, R.R.; van der Maarel, S.M. Hypomethylation of D4Z4 in 4q-linked and non-4q-linked facioscapulohumeral muscular dystrophy. *Nat. Genet.* **2003**, *35*, 315-317. [CrossRef] [PubMed]

34. Salort-Campana, E.; Nguyen, K.; Bernard, R.; Jouve, E.; Solé, G.; Nadaj-Pakleza, A.; Niederhauser, J.; Charles, E.; Ollagnon, E.; Bouhour, F.; et al. Low penetrance in facioscapulohumeral muscular dystrophy type 1 with large pathological D4Z4 alleles: A cross-sectional multicenter study. *Orphanet J. Rare Dis.* **2015**, *10*, 2. [CrossRef] [PubMed]

35. Sakellariou, P.; Kekou, K.; Fryssira, H.; Sofocleous, C.; Manta, P.; Panousopoulou, A.; Gounaris, K.; Kanavakis, E. Mutation spectrum and phenotypic manifestation in FSHD Greek patients. *Neuromuscul. Disord.* **2012**, *22*, 339-349. [CrossRef] [PubMed]

36. Hewitt, J.E.; Lyle, R.; Clark, L.N.; Valleley, E.M.; Wright, T.J.; Wijmenga, C.; van Deutekom, J.C.T.; Francis, F.; Sharpe, P.T.; Hofker, M.; et al. Analysis of the tandem repeat locus D4Z4 associated with facioscapulohumeral muscular dystropothy. *Hum. Mol. Genet.* **1994**, *3*, 1287-1295. [CrossRef]

37. Bakker, E.; Wijmenga, C.; Vossen, R.H.; Padberg, G.W.; Hewitt, J.; van der Wielen, M.; Rasmussen, K.; Frants, R.R. The FSHD-linked locus D4F104S1 (p13E-11) on 4q35 has a homologue on 10qter. *Muscle Nerve. Suppl.* **1995**, *2*, S39-S44. [CrossRef]

38. Deidda, G.; Cacurri, S.; Grisanti, P.; Vigneti, E.; Piazzo, N.; Felicetti, L. Physical mapping evidence for a duplicated region on chromosome 10qter showing high homology with the facioscapulohumeral muscular dystrophy locus on chromosome 4qter. *Eur. J. Hum. Genet.* **1995**, *3*, 155-167. [CrossRef]

39. Gabellini, D.; Green, M.R.; Tupler, R. Inappropriate gene activation in FSHD: A repressor complex binds a chromosomal repeat deleted in dystrophic muscle. *Cell* **2002**, *110*, 339-348. [CrossRef]

40. Yang, W.-M.; Yao, Y.L.; Seto, E. The FK506-binding protein 25 functionally associates with histone deacetylases and with transcription factor YY1. *EMBO J.* **2001**, *20*, 4814-4825. [CrossRef]

41. Yao, Y.L.; Yang, W.M.; Seto, E. Regulation of transcription factor YY1 by acetylation and deacetylation. *Mol. Cell. Biol.* **2001**, *21*, 5979-5991. [CrossRef]
42. Oei, S.L.; Shi, Y. Transcription Factor Yin Yang 1 Stimulates Poly(ADP-Ribosylation) and DNA Repair. *Biochem. Biophys. Res. Commun.* **2001**, *284*, 450–454. [CrossRef]

43. Cabianca, D.S.; Casa, V.; Bodega, B.; Xynos, A.; Ginelli, E.; Tanaka, Y.; Gabellini, D. A long ncRNA links copy number variation to a polycomb/trithorax epigenetic switch in FSHD muscular dystrophy. *Cell* **2012**, *149*, 819–831. [CrossRef] [PubMed]

44. Massah, S.; Hollebakken, R.; Labrecque, M.P.; Kolybaba, A.M.; Beischlag, T.V.; Prefontaine, G.G. Epigenetic characterization of the growth hormone gene identifies SmcHD1 as a regulator of autosomal gene clusters. *PLoS ONE* **2014**, *9*, e97535. [CrossRef] [PubMed]

45. Chen, K.; Hu, J.; Moore, D.L.; Liu, R.; Kessans, S.A.; Breslin, K.; Lucet, I.S.; Keniry, A.; Leong, H.S.; Parish, C.L.; et al. Genome-wide binding and mechanistic analyses of Smchd1-mediated epigenetic regulation. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, E3535–E3544. [CrossRef]

46. Mould, A.W.; Pang, Z.; Pakusch, M.; Tonks, I.D.; Stark, M.; Carrie, D.; Mukhopadhyay, P.; Seidel, A.; Ellis, J.J.; Deakin, J.; et al. Smchd1 regulates a subset of autosomal genes subject to monoallelic expression in addition to being critical for X inactivation. *Epigenetics Chromatin* **2013**, *6*, 19. [CrossRef] [PubMed]

47. Van Den Boogaard, M.L.; Lemmers, R.J.L.F.; Balog, J.; Wohlgemuth, M.; Auranen, M.; Mitsuhashi, S.; Van Der Vliet, P.J.; Straasheijm, K.R.; Van Den Akker, R.F.P.; Kriek, M.; et al. Mutations in DNMT3B Modify Epigenetic Repression of the D4Z4 Repeat and the Penetrance of Facioscapulohumeral Dystrophy. *Am. J. Hum. Genet.* **2016**, *98*, 1020–1029. [CrossRef] [PubMed]

48. de Greef, J.C.; Lemmers, R.J.L.F.; van Engelen, B.G.M.; Sacconi, S.; Venance, S.L.; Frants, R.R.; Tawil, R.; van der Maarel, S.M. Common epigenetic changes of D4Z4 in contraction-dependent and contraction-independent FSHD. *Hum. Mutat.* **2009**, *30*, 1449–1459. [CrossRef]

49. de Greef, J.C.; Wohlgemuth, M.; Chan, O.A.; Hansson, K.B.; Smeets, D.; Frants, R.R.; Weemaes, C.M.; Padberg, G.W.; van der Maarel, S.M. Hypomethylation is restricted to the D4Z4 repeat array in phenotypic FSHD. *Neurology* **2007**, *69*, 1018–1026. [CrossRef]

50. Dixit, M.; Ansseau, E.; Tassin, A.; Winokur, S.; Shi, R.; Qian, H.; Sauvage, S.; Mattéotti, C.; Van Acker, A.M.; Leo, O.; et al. DUX4, a candidate gene of facioscapulohumeral muscular dystrophy, encodes a transcriptional activator of PITX1. *Proc. Natl. Acad. Sci. USA* **2007**, *104*, 18157–18162. [CrossRef]

51. Kowaljow, V.; Marcowycz, A.; Ansseau, E.; Conde, C.B.; Sauvage, S.; Mattéotti, C.; Arias, C.; Corona, E.D.; Nuñez, N.G.; Leo, O.; et al. The DUX4 gene at the FSHD1A locus encodes a pro-apoptotic protein. *Neuromuscul. Disord.* **2007**, *17*, 611–623. [CrossRef]

52. Lemmers, R.J.L.F.; Van Der Vliet, P.J.; Klooster, R.; Sacconi, S.; Camaño, P.; Dauwerse, J.G.; Snider, L.; Straasheijm, K.R.; Van Ommen, G.J.; Padberg, G.W.; et al. A unifying genetic model for facioscapulohumeral muscular dystrophy. *Science* **2010**, *329*, 1650–1653. [CrossRef]

53. Ehrlich, M.; Gama-Sosa, M.A.; Huang, L.H.; Midgett, R.M.; Kuo, K.C.; Mccune, R.A.; Gehrke, C. Amount and distribution of 5-methylcytosine in human DNA from different types of tissues or cells. *Nucleic Acids Res.* **1982**, *10*, 2709–2721. [CrossRef] [PubMed]

54. Horard, B.; Eymery, A.; Fourel, G.; Vassetzky, N.; Puechberty, J.; Roizes, G.; Lebrigand, K.; Barbry, P.; Laugraud, A.; Gautier, C.; et al. Global analysis of DNA methylation and transcription of human repetitive sequences. *Epigenetics* **2009**, *4*, 66–77. [CrossRef] [PubMed]

55. Lander, E.S.; Linton, L.M.; Birren, B.; Nusbaum, C.; Zody, M.C.; Baldwin, J.; Devon, K.; Dewar, K.; Doyle, M.; Fitzhugh, W.; et al. Initial sequencing and analysis of the human genome. *Nature* **2001**, *409*, 860–921. [PubMed]

56. Francastel, C.; Magdinier, F. DNA methylation in satellite repeats disorders. *Essays Biochem.* **2019**, *63*, 757–771. [CrossRef] [PubMed]

57. Cantone, I.; Fisher, A.G. Epigenetic programming and reprogramming during development. *Nat. Struct. Mol. Biol.* **2013**, *20*, 282–289. [CrossRef] [PubMed]

58. Smith, Z.D.; Meissner, A. DNA methylation: Roles in mammalian development. *Nat. Rev. Genet.* **2013**, *14*, 204–220. [CrossRef]

59. Toyota, M.; Issa, J.P.J. The role of DNA hypermethylation in human neoplasia. *Electrophoresis* **2000**, *21*, 329–333. [CrossRef]

60. Bocker, M.T.; Hellwig, I.; Breiling, A.; Eckstein, V.; Ho, A.D.; Lyko, F. Genome-wide promoter DNA methylation dynamics of human hematopoietic progenitor cells during differentiation and aging. *Blood* **2011**, *117*, e182–e189. [CrossRef]
61. Martino, D.J.; Tulic, M.K.; Gordon, L.; Hodder, M.; Richman, T.; Metcalfe, J.; Prescott, S.L.; Saffery, R. Evidence for age-related and individual-specific changes in DNA methylation profile of mononuclear cells during early immune development in humans. *Epigenetics* **2011**, *6*, 1. [CrossRef]

62. Jones, P.A.; Baylin, S.B. The fundamental role of epigenetic events in cancer. *Nat. Rev. Genet.* **2002**, *3*, 415–428.
[CrossRef]

63. Feinberg, A.P. A genetic approach to cancer epigenetics. *Cold Spring Harb. Symp. Quant. Biol.* **2005**, *70*, 335–341. [CrossRef]

64. Wilson, A.S.; Power, B.E.; Molloy, P.L. DNA hypomethylation and human diseases. *Biochim. Biophys. Acta* **2007**, *1775*, 138–162. [CrossRef]

65. Bergman, Y.; Cedar, H. DNA methylation dynamics in health and disease. *Nat. Struct. Mol. Biol.* **2013**, *20*, 274–281. [CrossRef] [PubMed]

66. Gatchel, J.R.; Zoghbi, H.Y. Diseases of unstable repeat expansion: Mechanisms and common principles. *Nat. Rev. Genet.* **2005**, *6*, 743–755. [CrossRef] [PubMed]

67. Gijselinck, I.; Van Mossevelde, S.; Van Der Zee, J.; Sieben, A.; Engelborghs, S.; De Bleecker, J.; Ivanoiu, A.; Deryck, O.; Edbauer, D.; Zhang, M.; et al. The C9orf72 repeat size correlates with onset age of disease, DNA methylation and transcriptional downregulation of the promoter. *Mol. Psychiatry* **2016**, *21*, 1112–1124. [CrossRef] [PubMed]

68. Evans-Galea, M.V.; Carrodus, N.; Rowley, S.M.; Corben, L.A.; Tai, G.; Saffery, R.; Galati, J.C.; Wong, N.C.; Craig, J.M.; Lynch, D.R.; et al. FXN methylation predicts expression and clinical outcome in Friedreich ataxia. *Ann. Neurol.* **2012**, *71*, 487–497. [CrossRef] [PubMed]

69. Jones, T.I.; King, O.D.; Himeda, C.L.; Homma, S.; Chen, J.C.J.J.; Beermann, M.L.; Yan, C.; Emerson, C.P.; Miller, J.B.; Wagner, K.R.; et al. Individual epigenetic status of the pathogenic D4Z4 macrosatellite correlates with disease in facioscapulohumeral muscular dystrophy. *Clin. Epigenetics* **2015**, *7*, 37. [CrossRef]

70. Lemmers, R.J.L.F.; Goeman, J.J.; van der Vliet, P.J.; van Nieuwenhuizen, M.P.; Balog, J.; Vos-Versteeg, M.; Camano, P.; Ramos Arroyo, M.A.; Jerico, I.; Rogers, M.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. [CrossRef]

71. Gaillard, M.-C.; Roche, S.; Dion, C.; Tasmadjian, A.; Bouget, G.; Salort-Campana, E.; Vovan, C.; Chaix, C.; Broucqsault, N.; Morere, J.; et al. Differential DNA methylation of the D4Z4 repeat in patients with FSHD and asymptomatic carriers. *Neurology* **2014**, *83*, 733–742. [CrossRef] [PubMed]

72. Calandra, P.; Cascino, I.; Lemmers, R.J.L.F.L.F.; Galluzzi, G.; Teveroni, E.; Monforte, M.; Tasca, G.; Ricci, E.; Moretti, F.; Van Der Maarel, S.M.; et al. Allele-specific DNA hypomethylation characterises FSHD1 and FSHD2. *J. Med. Genet.* **2016**, *53*, 348–355. [CrossRef]

73. Nikolic, A.; Jones T.I.; Govi, M.; Mele, F.; Maranda, L.; Sera, F.; Ricci, G.; Ruggiero, L.; Vercelli, L.; Portaro, S.; et al. Interpretation of the epigenetic signature of facioscapulohumeral muscular dystrophy in light of genotype-phenotype studies. *IJMS* under review.

74. Roche, S.; Dion, C.; Broucqsault, N.; ere, C.; Gaillard, M.-C.; Robin, J.D.; Lagarde, A.; Puppo, F.; Vovan, C.; Chaix, C.; et al. Methylation hotspots evidenced by deep sequencing in patients with facioscapulohumeral dystrophy and mosaicism. *Neurol. Genet.* **2019**, *5*, e372. [CrossRef] [PubMed]

75. Li, L.C.; Dahiya, R. MethPrimer: Designing primers for methylation PCRs. *Bioinformatics* **2002**, *18*, 1427–1431.
[CrossRef]

76. Lemmers, R.J.L.F.; Tawil, R.; Petek, L.M.; Balog, J.; Block, G.J.; Santen, G.W.E.; Amell, A.M.; van der Vliet, P.J.; Almomani, R.; Straasheijm, K.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.
[CrossRef] [PubMed]

77. Gaillard, M.-C.; Puppo, F.; Roche, S.; Dion, C.; Campana, E.S.; Mariot, V.; Chaix, C.; Vovan, C.; Mazaleyrat, K.; Tasmadjian, A.; et al. Segregation between SMCHD1 mutation, D4Z4 hypomethylation and Facio-Scapulo-Humeral Dystrophy: A case report. *BMC Med. Genet.* **2016**, *17*, 66. [CrossRef] [PubMed]

78. Kondo, T. Whole-genome methylation scan in ICF syndrome: Hypomethylation of non-satellite DNA repeats D4Z4 and NBL2. *Hum. Mol. Genet.* **2000**, *9*, 597–604. [CrossRef]

79. Dion, C.; Roche, S.; Laberthonnière, C.; Broucqsault, N.; Mariot, V.; Xue, S.; Gurzau, A.D.; Nowak, A.; Gordon, C.T.; Gaillard, M.-C.; et al. SMCHD1 is involved in *de novo* methylation of the *DUX4* -encoding D4Z4 macrosatellite. *Nucleic Acids Res.* **2019**, *47*, 2822–2839. [CrossRef] [PubMed]
80. Balog, J.; Goossens, R.; Lemmers, R.J.L.F.; Straasheijm, K.R.; Van Der Vliet, P.J.; Van Den Heuvel, A.; Cambieri, C.; Capet, N.; Feasson, L.; Manel, V.; et al. Monosomy 18p is a risk factor for facioscapulohumeral dystrophy. *J. Med. Genet.* **2018**, *55*, 469–478. [CrossRef]

81. Nguyen, K.; Walrafen, P.; Bernard, R.; Attarian, S.; Chaix, C.; Vovan, C.; Renard, E.; Dufrane, N.; Pouget, J.; Vannier, A.; et al. Molecular combing reveals allelic combinations in facioscapulohumeral dystrophy. *Ann. Neurol.* **2011**, *70*, 627–633. [CrossRef]

82. Nguyen, K.; Puppo, F.; Roche, S.; Gaillard, M.-C.; Chaix, C.; Lagarde, A.; Pierret, M.; Vovan, C.; Olschwang, S.; Salort-Campana, E.; et al. Molecular combing reveals complex 4q35 rearrangements in Facioscapulohumeral dystrophy. *Hum. Mutat.* **2017**, *38*, 1432–1441. [CrossRef]

83. Nguyen, K.; Broucqsault, N.; Chaix, C.; Roche, S.; Robin, J.D.; Vovan, C.; Gerard, L.; Mégarbané, A.; Urtizberea, J.A.; Bellance, R.; et al. Deciphering the complexity of the 4q and 10q subtelomeres by molecular combing in healthy individuals and patients with facioscapulohumeral dystrophy. *J. Med. Genet.* **2019**, *56*, 590–601. [CrossRef]

84. Lengerke, C.; Daley, G.Q. Disease models from pluripotent stem cells: Turning back time in disease pathogenesis. *Ann. N. Y. Acad. Sci.* **2009**, *1176*, 191–196. [CrossRef] [PubMed]

85. Urbach, A. Modeling for Lesch-Nyhan Disease by Gene Targeting in Human Embryonic Stem Cells. *Stem Cells* **2004**, *22*, 635–641. [CrossRef] [PubMed]

86. Lowry, W.E.; Richter, L.; Yachechko, R.; Pyle, A.D.; Tchieu, J.; Sridharan, R.; Clark, A.T.; Plath, K. Generation of human induced pluripotent stem cells from dermal fibroblasts. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 2883–2888. [CrossRef] [PubMed]

87. Eiges, R.; Urbach, A.; Malcov, M.; Frumkin, T.; Schwartz, T.; Amit, A.; Yaron, Y.; Eden, A.; Yanuka, O.; Benvenisty, N.; et al. Developmental Study of Fragile X Syndrome Using Human Embryonic Stem Cells Derived from Preimplantation Genetically Diagnosed Embryos. *Cell Stem Cell* **2007**, *1*, 568–577. [CrossRef] [PubMed]

88. Okano, M.; Bell, D.W.; Haber, D.A.; Li, E. DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. *Cell* **1999**, *99*, 247–257. [CrossRef]

89. Huang, K.; Wu, Z.; Liu, Z.; Hu, G.; Yu, J.; Chang, K.H.; Kim, K.P.; Le, T.; Faull, K.F.; Rao, N.; et al. Selective demethylation and altered gene expression are associated with ICF syndrome in human-induced pluripotent stem cells and mesenchymal stem cells. *Hum. Mol. Genet.* **2014**, *23*, 6448–6457. [CrossRef]

90. Avitzour, M.; Mor-Shaked, H.; Yanovsky-Dagan, S.; Aharoni, S.; Altarescu, G.; Renbaum, P.; Eldar-Geva, T.; Schonberger, O.; Levy-Lahad, E.; Epsztejn-Litman, S.; et al. FMR1 epigenetic silencing commonly occurs in undifferentiated fragile X-affected embryonic stem cells. *Stem Cell Rep.* **2014**, *3*, 699–706. [CrossRef]

91. De Iaco, A.; Planet, E.; Coluccio, A.; Verp, S.; Duc, J.; Trono, D. DUX-family transcription factors regulate zygotic genome activation in placental mammals. *Nat. Genet.* **2017**, *49*, 941–945. [CrossRef]

92. Hendrickson, P.G.; Doráis, J.A.; Grow, E.J.; Whiddon, J.L.; Lim, J.W.; Wike, C.L.; Weaver, B.D.; Pflueger, C.; Emery, B.R.; Wilcox, A.L.; et al. Conserved roles of mouse DUX and human DUX4 in activating cleavage-stage genes and MERVL/HERVL retrotransposons. *Nat. Genet.* **2017**, *49*, 925–934. [CrossRef]

93. Geng, L.N.; Yao, Z.; Snider, L.; Fong, A.P.; Cech, J.N.; Young, J.M.; vanderMaarel, S.M.; Ruzzo, W.L.; Gentleman, R.C.; Tawil, R.; et al. DUX4 Activates Germline Genes, Retroelements, and Immune Mediators: Implications for Facioscapulohumeral Dystrophy. *Dev. Cell* **2012**, *22*, 38–51. [CrossRef]

94. Snider, L.; Geng, L.N.; Lemmers, R.J.L.F.; Kyba, M.; Ware, C.B.; Nelson, A.M.; Tawil, R.; Filippova, G.N.; van der Maarel, S.M.; Tapscott, S.J.; et al. Facioscapulohumeral Dystrophy: Incomplete Suppression of a Retrotransposed Gene. *PLoS Genet.* **2010**, *6*, e1001181. [CrossRef] [PubMed]

95. Zechner, U.; Nolte, J.; Wolf, M.; Shirneshan, K.; Hajj, N.E.; Weise, D.; Kaltwasser, B.; Zovoilis, A.; Haaf, T.; Engel, W. Comparative methylation profiles and telomerase biology of mouse multipotent adult germline stem cells and embryonic stem cells. *MHR Basic Sci. Reprod. Med.* **2009**, *15*, 345–353. [CrossRef] [PubMed]

96. Marion, R.M.; Strati, K.; Li, H.; Tejera, A.; Schoeftner, S.; Ortega, S.; Serrano, M.; Blasco, M.A. Telomeres Acquire Embryonic Stem Cell Characteristics in Induced Pluripotent Stem Cells. *Cell Stem Cell* **2009**, *4*, 141–154. [CrossRef] [PubMed]

97. Marion, R.M.; Blasco, M.A. Telomeres and telomerase in adult stem cells and pluripotent embryonic stem cells. *Adv. Exp. Med. Biol.* **2010**, *695*, 118–131. [PubMed]
98. Zeng, W.; de Greef, J.C.; Chen, Y.-Y.; Chien, R.; Kong, X.; Gregson, H.C.; Winokur, S.T.; Pyle, A.; Robertson, K.D.; Schmiesing, J.A.; et al. Specific loss of histone H3 lysine 9 trimethylation and HP1gamma/cohesin binding at D4Z4 repeats is associated with facioscapulohumeral dystrophy (FSHD). *PLoS Genet.* **2009**, *5*, e1000559.
[CrossRef]

99. Zeng, W.; Chen, Y.-Y.; Newkirk, D.A.; Wu, B.; Balog, J.; Kong, X.; Ball, A.R.; Zanotti, S.; Tawil, R.; Hashimoto, N.; 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. [CrossRef]

100. Jiang, G.; Yang, F.; van Overveld, P.G.M.; Vedanarayanan, V.; van der Maarel, S.; Ehrlich, M. Testing the position-effect variegation hypothesis for facioscapulohumeral muscular dystrophy by analysis of histone modification and gene expression in subtelomeric 4q. *Hum. Mol. Genet.* **2003**, *12*, 2909–2921. [CrossRef]

101. Ottaviani, A.; Rival-Gervier, S.; Boussouar, A.; Foerster, A.M.; Rondier, D.; Sacconi, S.; Desnuelle, C.; Gilson, E.; Magdinier, F. The D4Z4 macrosatellite repeat acts as a CTCF and A-type lamins-dependent insulator in Facio-Scapulo-Humeral dystrophy. *PLoS Genet.* **2009**, *5*, e1000394. [CrossRef]

102. Bodega, B.; Ramirez, G.D.C.; Grasser, F.; Cheli, S.; Brunelli, S.; Mora, M.; Meneveri, R.; Marozzi, A.; Mueller, S.; Battaglioli, E.; et al. Remodeling of the chromatin structure of the facioscapulohumeral muscular dystrophy (FSHD) locus and upregulation of FSHD-related gene 1 (FRG1) expression during human myogenic differentiation. *BMC Biol.* **2009**, *7*, 41. [CrossRef]

103. Cortesi, A.; Pesant, M.; Sinha, S.; Marasca, F.; Sala, E.; Gregoretti, F.; Antonelli, L.; Oliva, G.; Chiereghin, C.; Soldà, G.; et al. 4q-D4Z4 chromatin architecture regulates the transcription of muscle atrophic genes in facioscapulohumeral muscular dystrophy. *Genome Res.* **2019**, *29*, 883–895. [CrossRef]

104. Petrov, A.; Pirozhkova, I.; Carnac, G.; Laoudj, D.; Lipinski, M.; Vassetzky, Y.S. Chromatin loop domain organization within the 4q35 locus in facioscapulohumeral dystropny patients versus normal human myoblasts. *Proc. Natl. Acad. Sci. USA* **2006**, *103*, 6982–6987. [CrossRef] [PubMed]

105. Kim, J.; Kollhoff, A.; Bergmann, A.; Stubbs, L. Methylation-sensitive binding of transcription factor YY1 to an insulator sequence within the paternally expressed imprinted gene, Peg3. *Hum. Mol. Genet.* **2003**, *12*, 233–245. [CrossRef] [PubMed]

106. Huichalaf, C.; Micheloni, S.; Ferri, G.; Caccia, R.; Gabellini, D. DNA methylation analysis of the macrosatellite repeat associated with FSHD muscular dystrophy at single nucleotide level. *PLoS ONE* **2014**, *9*, e115278. [CrossRef] [PubMed]

107. Wang, H.; Maurano, M.T.; Qu, H.; Varley, K.E.; Gertz, J.; Pauli, F.; Lee, K.; Canfield, T.; Weaver, M.; Sandstrom, R.; et al. Widespread plasticity in CTCF occupancy linked to DNA methylation. *Genome Res.* **2012**, *22*, 1680–1688. [CrossRef]

108. Hartweck, L.M.; Anderson, L.J.; Lemmers, R.J.; Dandapat, A.; Toso, E.A.; Dalton, J.C.; Tawil, R.; Day, J.W.; Van Der Maarel, S.M.; Kyba, M. A focal domain of extreme demethylation within D4Z4 in FSHD2. *Neurology* **2013**, *80*, 392–399. [CrossRef]

109. Petrov, A.; Allinne, J.; Pirozhkova, I.; Laoudj, D.; Lipinski, M.; Vassetzky, Y.S. A nuclear matrix attachment site in the 4q35 locus has an enhancer-blocking activity in vivo: Implications for the facio-scapulo-humeral dystrophy. *Genome Res.* **2008**, *18*, 39–45. [CrossRef]

110. Gaszner, M.; Felsenfeld, G. Insulators: Exploiting transcriptional and epigenetic mechanisms. *Nat. Rev. Genet.* **2006**, *7*, 703–713. [CrossRef]

111. Ghirlando, R.; Felsenfeld, G. CTCF: Making the right connections. *Genes Dev.* **2016**, *30*, 881–891. [CrossRef]

112. Nikolic, A.; Ricci, G.; Sera, F.; Bucci, E.; Govi, M.; Mele, F.; Rossi, M.; Ruggiero, L.; Vercelli, L.; Ravaglia, S.; et al. Clinical expression of facioscapulohumeral muscular dystrophy in carriers of 1-3 D4Z4 reduced alleles: Experience of the FSHD Italian National Registry. *BMJ Open* **2016**, *6*, 1–10. [CrossRef]

113. Shaw, N.D.; Brand, H.; Kupchinsky, Z.A.; Bengani, H.; Plummer, L.; Jones, T.I.; Erdin, S.; Williamson, K.A.; Rainger, J.; Stortchevoi, A.; et al. SMCHD1 mutations associated with a rare muscular dystrophy can also cause isolated arhinia and Bosma arhinia microphthalmia syndrome. *Nat. Genet.* **2017**, *49*, 238–248. [CrossRef]
114. Gordon, C.T.; Xue, S.; Yigit, G.; Filali, H.; Chen, K.; Rosin, N.; Yoshiura, K.I.; Oufadem, M.; Beck, T.J.; McGowan, R.; et al. De novo mutations in SMCHD1 cause Bosma arhinia microphthalmia syndrome and abrogate nasal development. *Nat. Genet.* **2017**, *49*, 249–255. [CrossRef] [PubMed]

115. Ricci, G.; Cammish, P.; Siciliano, G.; Tupler, R.; Lochmuller, H.; Evangelista, T. Phenotype may predict the clinical course of facioscapolohumeral muscular dystrophy. *Muscle Nerve* **2019**, *59*, 711–713. [CrossRef] [PubMed]

<figure><img src="image_7.png" /></figure>

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).