Estradiol differentially regulates DUX4, β-catenin and PAX3/PAX7 in primary myoblasts of facioscapulohumeral muscular dystrophy patients

Ceren Hangul; Esin Guvenir Celik; Hacer Kaya; Onur Eroglu; Hilmi Uysal; Sibel Berker Karauzum

doi:10.1515/tjb-2020-0351

Article Open Access

Estradiol differentially regulates DUX4, β-catenin and PAX3/PAX7 in primary myoblasts of facioscapulohumeral muscular dystrophy patients

Ceren Hangul , Esin Guvenir Celik , Hacer Kaya , Onur Eroglu , Hilmi Uysal and Sibel Berker Karauzum

Published/Copyright: January 5, 2021

Published by

Become an author with De Gruyter Brill

Submit Manuscript Author Information Explore this Subject

From the journal Turkish Journal of Biochemistry Volume 46 Issue 4

Abstract

Objectives

There is a clinical variability and heterogeneity among Facioscapulohumeral Muscular Dystrophy (FSHD) patients. Escalation after menopause in women, early onset in men suggests that estrogen might be a protective factor on the course of FSHD. In spite of few molecular studies supporting the protective role of estrogen in FSHD in vitro, there is no study revealing the effect of estradiol on the protein levels of DUX4, β-catenin and PAX3/PAX7. In present study, we investigated the effect of estradiol treatment on the expressions of DUX4, β-catenin and PAX3/PAX7 protein levels.

Materials and Methods

Primary myoblasts of 63 and 71 years old (63yM/71yM) males; 47 years old (47yF) female FSHD patients were used. Cells were processed under these conditions; (i) untreated, (ii) 10 nM-30 min estradiol and (iii) 10 nM-4 h estradiol treated. The expression of DUX4, PAX3/PAX7 and β-catenin were examined by western-blotting.

Results

Expression of DUX4 significantly downregulated after 4 h treatment of estradiol while PAX3/PAX7 56 kDa variant expression upregulated in 71yM cells. β-catenin and PAX3 expression was variable among the samples.

Conclusion

Our results suggest that estrogen might be a palliative treatment option via downregulation of DUX4 protein in DUX4 expressing FSHD patients.

ÖZ

Amaç

Fasioskapulohumeral Musküler Distrofi (FSHD) hastaları arasında klinik değişkenlik ve heterojenite gözlenmektedir. Erkeklerde hastalığın daha erken başlaması, kadınlarda ise menopoz ile bulguların şiddetlenmesi, östrojenin FSHD’nin seyri üzerinde koruyucu bir faktör olabileceğini düşündürmüştür. FSHD’deki koruyucu rolü üzerine birkaç in vitro çalışma mevcut olmakla beraber, östradiolün FSHD’ de DUX4, β-katenin ve PAX3/PAX7 proteinlerinin düzeyleri üzerine olan etkisini ortaya koyan bir çalışma bulunmamaktadır. Mevcut çalışmamızda, östradiol muamelesinin DUX4, β-katenin ve PAX3/PAX7 protein düzeylerine etkisi araştırılmıştır.

Gereç ve Yöntem

63 ve 71 yaşında erkek(63yE/71yE) ve 47 yaşında kadın (47yK) FSHD hastalarına ait primer miyoblast hücreleri; (i) kontrol (hiçbir muamele yapılmamış), (ii) 10 nM–30 dakika estradiol ve (iii) 10 nM-4-saat östradiol muamelesi sonrasında; DUX4, β-katenin ve PAX3/PAX7 ekspresyonunu belirlemek üzere western blot metodu ile incelenmiştir.

Bulgular

Östradiol ile 4 saat muamele sonucunda 71yE hastasına ait miyoblast hücrelerinde DUX4 ekspresyonunun down-regülasyonu, PAX3/PAX7 56 kDa varyantının up-regülasyonu gözlenmiştir. β-katenin ve PAX3 proteinleri, miyoblast örnekleri arasında değişken ekspresyon yanıtları sergilemiştir.

Sonuç

Bu çalışmanın sonuçları, DUX4 eksprese eden FSHD hastalarında östradiolün DUX4 down-regülasyonu yaparak palyatif bir tedavi seçeneği sunabileceğine işaret etmektedir.

Keywords: β-catenin; DUX4; estradiol; facioscapulohumeral muscular dystrophy; FSHD; PAX3; PAX7; skeletal muscle

Anahtar Kelimeler: β-katenin; DUX4; Fasioskapulohumeral Musküler Distrofi; FSHD; iskelet kası; östradiol; PAX3; PAX7

Introduction

Facioscapulohumeral Muscular Dystrophy (FSHD) (OMIM#158900), is an autosomal dominant genetic disease. Deletion in D4Z4 macrosatellite repeat located at chromosome 4q35 region is responsible from the disease. In normal population, D4Z4 contains 11 to 100 repeats and each unit consists of 3,300 bases involving DUX4 gene. In FSHD1 (95% of the patients) this repeat number is between 1 and 10 resulting from a deletion in this area [1]. This deletion causes chromatin relaxation ending up with the activation of DUX4 which is normally inactive in adult tissues except testis in men. If there is qA allele at the end of this repeat, this activation leads to the stabilization of DUX4 mRNA and the change of transcriptional network in FSHD muscle [2].

DUX4 is a double homeobox protein, especially expressed in embryological development [3], [4], its structure is similar to PAX protein family with their homeodomains. In one study, it was revealed that PAX3 and PAX7 are the only members that can compete with DUX4, while other homeodomain proteins as PAX6, Pitx2c, OTX1, Rax, Hesx1, MIXL1 and Tbx1 do not have this property [5]. Supporting this, overexpression of PAX3 and PAX7 in FSHD cell lines resulted in a decrease of DUX4 revealing a converse relation between DUX4 and PAX3/PAX7 [6]. DUX4 gene activation is a hallmark of FSHD. In addition to this, PAX7 target gene repression had been revealed as a biomarker in FSHD [7] just recently. PAX3 and PAX7 have splice variants; PAX3 have eight splice variants as 44.8/45.2/22.7/53.5/24.1/53/56/53.3 kDa and PAX7 have three splice variants as 56.9/55.1/56.7 kDa. Because it has not been clearly identified which variant or variants take a role in the pathology of FSHD, common antibody had been used.

DUX4 is a transcription factor and has many target proteins [8]. It was revealed that DUX4 overexpression resulted in either upregulation or downregulation of wide range of genes [9]. Until now, it has been difficult to attribute one or few responsible gene(s) for explaining the pathophysiology. In 2015, Banerji et al. used a new network approach and suggested β-catenin as a central protein of this network [10] supporting previous studies that indicated impaired Wnt/β-catenin pathway in FSHD [11], [12].

FSHD is a progressive disease with an onset around 20–30s [13]. As named, Facioscapulohumeral Muscular Dystrophy affects primarily facial, scapular, humeral muscles and tibialis anterior muscles. There are clear variations in onset and clinical severity even in same family members carrying same repeat unit number. This indicates to the existence of additional modifier factors affecting severity. In men, onset is earlier and phenotype is more severe compared to women [14]. In women, it is observed that after childbirth [15], [16] and especially after menopause, there is clear aggravation of symptoms arousing a suspicion for the effect of estrogen on FSHD severity. Mul et al. identified no relation of clinical severity with lifetime estrogen exposure [17]. However, in our previous clinical study ratio of estradiol and progesterone to testosterone had been related to clinical severity independent of D4Z4 repeat unit [18]. Supportively, in one clinical study it is observed that anti-estradiol cancer therapy caused an aggravation of FSHD symptoms in contrast to other chemotherapeutic agents [19]. Later on, Teveroni et al. revealed the protective molecular effect of estradiol on FSHD cell lines by showing hypermethylation of DUX4’s binding sites on DNA. Although direct regulation of DUX4 couldn’t be shown, they revealed the downregulation of DUX4 target genes via epigenetic modulation after estrogen treatment [20]. Independent from DUX4, estrogen had also been shown to contribute satellite cell maintenance and muscle regeneration in females [21].

Based on the clinical difference between men and women in FSHD and in the light of aforementioned studies in literature, the purpose of this study was to reveal whether estrogen had an effect on the pathophysiology of FSHD at a molecular level. For this purpose, four different transcription factors which were indicated as key factors in literature had been investigated (Figure 1). These key transcription factors were; (i) DUX4 protein which locates in each D4Z4 repeat at chromosome 4q35 region; (ii) β-catenin which had been demonstrated as a central protein in the network of the FSHD, (iii) PAX3 and PAX7 which are known to compete with DUX4 for DNA binding and also essential in embryologic development of skeletal muscles. Since the best reflecting in vitro study would be to evaluate the effect of estradiol on FSHD primary myoblasts, patient derived myoblasts had been used for investigation.

Figure 1:

Representative schematic of the study.

Materials and methods

This study had been approved by Akdeniz University Faculty of Medicine Clinical Research Ethical Committee; Decision Number: 284, Date: 28.10.2015.

Human myoblast cell culture

FSHD myoblast cell lines were obtained from Massachusetts University in the United States, from Professor Rabi Tawil and Don Henderson team. FSHD cell lines belonging to 63 and 71 years old males; 47 years old female who were diagnosed with FSHD1 had been used in the experiments. These cell lines are DUX4 positive cell lines in the level of DUX4 mRNA (Table 1).

Table 1:

Cell lines and their detailed information.

Cell line code	Age at Bx*	Gender	Onset	CSS**	Repeat length & allele	Bx muscle
71yM	71	M	40	8	31 kb A	Deltoid
63yM	63	M	15	6	23 kb A	Vastus lateralis
47yF	47	F	47	1	19 kb A	Vastus lateralis

*Bx, Biopsy; **CSS, Clinic Severity Score. A scoring system for severity of the disease and could be reached on https://www.urmc.rochester.edu/MediaLibraries/URMCMedia/fields-center/documents/ClinicalSeverityScoring.pdf. Reference: [22], [23].

Selected cell lines reflect; differently sized contracted repeat units, variable clinical severity, both male and female genders and different muscles where muscle biopsy had been sampled (Table 1). Instead of performing all procedures in one cell line, three cell lines had been used, to be able to reflect the clinical variability in FSHD patients [13], FSHD myoblast cells were grown F10 Nutrient Medium (Gibco11550-043) containing 20% heat-inactivated fetal bovine serum(Gibco26140-079), 1% penicillin/streptomycin (Gibco15070-063 10 ng/mL), Fibroblast Growth Factor (bFGF) (PromegaG5071), 1 µM Dexamethazone (Sigma water soluble D2915), at 37 °C containing 5% CO₂ incubator.

FSHD myoblast cells were plated in 75 cm² flasks. Cells were allowed to attach to flasks overnight. Media was replenished for a period of 2–3 days after adhesion. As the cell number of the flasks reached 60–80% confluency, medium was removed for cell passaging; the cells were washed with 2 mL of Phosphate Buffered Saline (PBS), incubated for 3 min with 3 mL of trypsin-EDTA (Ethylene-Diamine Tetraacetic Acid) (0.1%). At the end of this period, 3 mL of F10 medium was added to the flasks to stop trypsin activity. Cells were collected and centrifuged at 1,800 rpm for 7 min, then the supernatant was discarded, and 1 mL of fresh medium was added onto cell pellet. Homogenized cell suspension was counted in the hemocytometer. 1.5 × 10⁶ cells were plated to 75 cm² flasks. After two or three passages, when reached 60–80% confluency, adhesive cells in different plates were treated with 10 nM estradiol (E2-Sigma Aldrich-E2758) for 30 min, 4 h and non-treated (control). All cells were incubated at 37 °C in 5% CO₂.

Treatments

Estradiol had been dissolved in ethanol, diluted in medium to reach 10 nM concentration. 10 nM concentration had been identified from estradiol-FSHD cell study by Teveroni et al. 2017 and other two studies with myoblasts cultures in literature [20], [24], [25]. Estradiol treatment time had been identified as 72 h from Teveroni et al. 2017 [20]. However, after observing morphological destruction and high amounts of cell death in first experiments, literature re-scanning had been performed and time period was identified as two stages: 30 min and 4 h [26] to reveal time dependent expression of investigated proteins.

Western blotting

Cells were collected by 1XPBS and centrifuged for 2 min at 13,200 rpm. The supernatant was discarded and radioimmuno-precipitation assay (RIPA) solution was added as an elution buffer on the cells. Phenyl-methyl-sulfonyl-fluoride (PMSF) was added to prevent degradation of proteins with protease activity just before addition to the RIPA solution. The samples were incubated at room temperature for 20 min on a shaker and centrifuged at +4 °C for 20 min at 13,200 rpm. Bradford method was used to determine the amount of protein [27]. Lysates were separated by SDS-PAGE gel electrophoresis using 12% polyacrilamide gels and then transferred to nitrocellulose membranes with semi-dry transfer system (iBlot2). Protein transfer was controlled by Ponceau Red staining.

Nitrocellulose membranes were blocked for 1 h with 5% milk powder without fat in Tris Buffered Saline (TBS) solution. After blocking, membranes were incubated with primary antibodies at 4 °C overnight in the same conditions. Primary antibodies used in 1:1,000 dilution as follows: anti DUX4 rabbit mab (ABCAM-E5-5), anti PAX3/PAX7 mouse mab (Santa Cruz-365843), anti β-catenin mab (Cell Signaling 9582S), β-actin (Cell Signaling) and GAPDH (Cell Signaling-2118).

After primary antibody labeling, membranes were washed with Tris Buffered Saline with Tween 20 (TBS-T) for 5 min of three times and then incubated with secondary antibodies attached with Horse Radish Peroxidase (HRP). For secondary antibody incubation 1:5,000 dilution was used for 1 h at 4 °C. Proteins were screened using ECL (ThermoFisher Scientific) in Syngene G:Box Chemi XRQ system without film.

For the proteins which have similar molecular weight and band pattern, stripping and re-staining had been performed. 0.5 M buffer with 6.8 pH (containing 20 mL SDS 10%, 12.5 mL Tris HCl, 0.8 mL β-mercaptoethanol, 67.5 mL distilled water) (abcam) had been used for harsh stripping protocol. Buffer solution was pre-heated to 50 °C; nitrocellulose membranes had been covered by buffer and incubated in 50 °C for 30–45 min with some agitation. Stripping solution was disposed, and membrane rinsed with water for 2 min. Then membrane was washed with TBS-T for 5 min. After this step, membrane was blocked with 5% milk powder and aforementioned steps were repeated for antibody labeling.

Statistical analysis

“ImageJ application” program is used for the density measurement of the monitored bands. Each protein band had been measured three times and mean value of three measurement had been used. For the calculation of relative expression levels; each value had been divided to values of β-actin. Graphs with error bars including p values (determined by one-way Anova test) had been prepared in GraphPad Prism.

Results

Three cell lines belonging to 71yM, 63yM, 47yF were investigated as following: (i) control (no estradiol treatment), (ii)30 min 10 nM 17-β estradiol (E2) treatment, (iii) 4 h 10 nM 17-β estradiol (E2) treatment.

The effect of E2 on DUX4 expression

DUX4 expression decreased to zero at fourth hour after E2 treatment in 71yM (p<0.0001) (Figure 2A). In the myoblasts of 63yM and 47yF there were no DUX4 protein expression neither in control myoblasts nor E2 treated myoblasts (Figure 2B, C).

Figure 2:

Demonstration of time-dependent DUX4 expression by western blotting in the presence of E2. (A) in 71yM, **** p<0.0001, (B) in 63yM, (C) in 47yF. Control (untreated); 30 min (exposure to 10 nM E2 for 30 min); 4 h (exposure to 10 nM E2 for 4 h). In western blotting images; above DUX4, below β-actin had been demonstrated.

The effect of E2 on β-catenin expression

In 71yM, β-catenin exhibited significant increase at 30 min (p=0.0084) and 4 h (p=0.0376) after estradiol treatment compared to control myoblasts (Figure 3A). Similar results were observed in 47yF (β-catenin increased at 30 min [p=0.0269] and 4 h treatment [p=0.0022] compared to control myoblasts [Figure 3C]) while in 63yM, after 4 h treatment, β-catenin expression was found to be lowest level (p<0.0001) compared to control myoblasts (Figure 3B). These results revealed divergent β-catenin responses to E2 treatment in FSHD myoblasts.

Figure 3:

Demonstration of time-dependent β-catenin expression by western blotting in the presence of E2. (A) in 71yM, ** p=0.0084, * p=0.0376, (B) in 63yM **** p<0.0001, (C) in 47yF * p=0.0269, ** p=0.0022. Control (Untreated); 30 min (exposure to 10 nM E2 for 30 min); 4 h (exposure to 10 nM E2 for 4 h). In western blotting images; above β-catenin, below β-actin had been demonstrated.

The effect of E2 on PAX3/PAX7 expression

PAX3/PAX7 labeled variant which expressed in all myoblasts belonging to three FSHD patients was 45 kDa. The expression of this variant was also present in control myoblasts where no E2 treatment, and the level of this protein exhibited divergent responses to E2 treatment similar to β-catenin while these responses were just in an opposite direction that β-catenin exhibited. In 71yM its expression gradually decreased in time dependent manner (p<0.0001), after 4 h of E2 treatment there were no expression (Figure 4A). Similar to that, in 47yM expression was found to be the lowest level at 4 h after E2 (p<0.0001) (Figure 4B). In 63yM, in contrast, its expression reached the highest level at 4 h after E2 (p<0.0001) (Figure 4C).

Figure 4:

Demonstration of time-dependent 45 kDa PAX3 expression by western blotting in the presence of E2. (A) 71yM **** p<0.0001, (B) 45 kDa PAX3 variant in 63yM **** p<0.0001, (C) 45 kDa PAX3 variant in 47yF **** p<0.0001. Control (Untreated); 30 min (exposure to 10 nM E2 for 30 min); 4 h (exposure to 10 nM E2 for 4 h). In western blotting images; above PAX3, below β-actin had been demonstrated. Because PAX3 and β-actin have similar weight, stripping protocol had been applied.

Additionally, in 71yM, expression of another PAX3/PAX7 variant which was 56 kDa had been observed after 4 h of E2 (p<0.0001), while there were no expression of this protein neither in non-E2 treated myoblasts nor 30 min E2 treated myoblasts (Figure 5A). 56 kDa variant was observed to be expressed only in 71yM myoblasts similar to DUX4 of whose expression zeroized at same time point.

Figure 5:

(A) Demonstration of time-dependent 56 kDa PAX3/PAX7 expression by western blotting in the presence of E2 in 71yM **** p<0.0001, (B) Demonstration of time-dependent GAPDH expression by western blotting in the presence of E2 without housekeeping protein, (C) Demonstration of time-dependent GAPDH expression **** p<0.0001.

The effect of E2 on GAPDH

GAPDH was used as housekeeping protein at first experiments. However, GAPDH protein bands were unexpectedly observed to be dependent to E2 treatment (Figure 5B). Repetition of GAPDH experiment via using β-actin as housekeeping protein, revealed that GAPDH expression had clearly been affected from E2. After 4 h treatment GAPDH expression dramatically decreased in time dependent manner compared to control (p<0.0001) (Figure 5C). Therefore, β-actin had been replaced as housekeeping protein in continuing experiments.

Discussion

In present study, the effect of E2 on FSHD myoblast had been researched through DUX4, β-catenin and PAX3/PAX7 protein levels. Teveroni et al. (2017) observed no change of DUX4 mRNA levels. On the other hand, after E2 treatment, downregulation of DUX4 target proteins (such as TRIM43 and ZSCAN4) had been revealed in differentiating FSHD myoblasts [20]. In our study, however, DUX4 protein level had shown to be affected by E2 treatment. After 4 h of E2, the expression of DUX4 was zeroised in myoblasts of 71yM. Inconsistency in the results of our and the previous study might arise from (i) using different cell types as myoblast cells vs. differentiating cells (ii) having different estrogen exposure times as 4 vs. 72 h and (iii) identifying different DUX4 structures as mRNA vs. protein. Despite divergent results, both of these studies underpin the protective role of estrogen in FSHD via DUX4.

In muscles of FSHD patients; DUX4 is expressed only in 1 of every 1,000 nuclei [28], [29] and usually not detectable in every FSHD patient cell line due to its low levels of expression [5]. In our study, DUX4 protein had been detected in untreated primary myoblasts of solely 71yM but not in myoblasts of 63yM or 47yF. The reason why DUX4 expression is seen only in this patient may be due to the following differences: 71yM’s Clinical Severity Score (CSS) was eight with highest score in between cell lines, on the other hand 63yM’s score was 6 and 47yF’s score was one, at the time of the biopsy. Also, variance in estrogen response between women and men might be conceivable additional factor effecting DUX4 expression in myoblasts. In 47yF, due to gender factor, relative higher estrogen response may suppress DUX4 expression. Another critical point to discriminate is the muscle group where the biopsy had been sampled. Muscle involvement of FSHD usually starts from the face and upper extremity; then proceeds to the lower extremity [30]. In 71yM, biopsy was taken from the deltoid muscle, while in 63yM and 47yF biopsies were taken from vastus lateralis muscle which is later involved at the progression of FSHD (Table 1). Because deltoid muscle is expected to be more severely affected than vastus lateralis, 71yM is expected to have higher DUX4 levels in which sampling was from deltoid muscle.

Temporal specific action of Wnt/β-catenin signaling has multiple functions in myogenesis including myoblast differentiation, muscle regeneration and satellite cell proliferation [31]. Role of β-catenin in FSHD pathophysiology is not completely clear yet, however, there are studies that investigated its function. In one study, opposite interplay had been suggested for DUX4 and β-catenin. It was concluded that DUX4 transcription was suppressed when Wnt/β-catenin signaling was activated in FSHD1 and FSHD2 primary human myoblasts [11]. Supporting this result, Banerji et al. revealed that DUX4 disrupted Wnt/β-catenin signaling pathway in C2C12 mice myoblasts [10]. In a study with human peripheral blood mononuclear cells of FSHD patients, a significant decrease in both cytosolic and nuclear fraction of β-catenin was observed compared to the control group [12]. In present study we observed a decrease in DUX4 and an increase in β-catenin expression in 71yM via E2 treatment. This reverse relation of DUX4 and β-catenin is compatible with previous researches in literature and this response in 71yM with E2 treatment supported the protective role of estrogen. However, non-DUX4 expressing myoblasts (63yM, 47yF) exhibited divergent responses to E2. These divergent responses put β-catenin into a more complex position in FSHD pathophysiology.

There is not any expression of neither PAX3 nor PAX7 protein from the advanced stages of skeletal muscle differentiation. Expression of PAX3 and PAX7 is limited to undifferentiated muscle cells [32], that’s why in our study, both PAX3 and PAX7 proteins were expected to be expressed in myoblasts. We observed two differently sized protein bands for PAX3 and PAX7 in western blotting as; 56 and 45 kDa in 71yM.

Since both PAX3 and PAX7 have 56 kDa splice forms, this band might belong to any of these forms. Strikingly, 56 kDa variant of PAX3/PAX7 expression had been stimulated 4 h after estradiol treatment, while DUX4 expression had been downregulated with estradiol treatment and there was no expression at the time of 4 h. 56 kDa PAX3/PAX7 response support the competition with DUX4 that was revealed by Bosnakovski et al. [5]. They also revealed when PAX3/PAX7 upregulated DUX4 is no longer toxic. Altogether, previously published studies together with our results support the protective role of estrogen by upregulating 56 kDA variant of PAX3/PAX7 expression and downregulating DUX4 expression.

Since PAX7 does not have a 45 kDa form, 45 kDa protein band corresponds to 44.8 kDa variant or 45.2 kDa variant of PAX3. This 45 kDa variant of PAX3 observed in all myoblasts (71yM, 63yM, 47yF) and exhibited divergent responses to estradiol treatment, interestingly just in an opposite direction to β-catenin responses. The divergent responses of PAX3 to estradiol might be explained with different variant expression which was not possible show technically. Not only this, but also there are differences both in the formation and development of facial and body muscles embryonically [33]. Therefore, different muscle groups may show different gene expression patterns. Good proof for that is a study with rats that investigated PAX3 and PAX7 levels in between extensor digitorum longus and in soleus muscles. As a result, PAX3 mRNA levels were found to be similar to each other in different muscle groups, while PAX7 mRNA levels were not similar to each other [32]. Because of that, independent from the severity, sampling from different muscle origin from 71yM and 63yM/47yF might also explain the variance in the expression of DUX4 protein levels. However, β-catenin exhibited also divergent responses among cell lines. An additional possibility is the variations in the estrogen response element (ERE) sequence that affect the binding and the transcriptional activation [34] of PAX3 and β-catenin leading to divergent responses. Because each myoblast belongs to one individual’s unique genome, genetic or epigenetic polymorphisms located in ERE sequences might potentially affect the responses.

GAPDH have six different spliced variants; 36.1 kDa (three variants), 31.5 kDa (two variants), 27.9 kDa (one variant). In our study, we observed a band profile at 32 kDa level. This result indicates that 31.5 kDa isoform(s) is expressed in FSHD myoblasts. Expression of 31.5 kDa variant of GAPDH protein isoform in FSHD myoblasts is a new knowledge for the literature. Additionally, this variant shows to be affected by estradiol treatment. Joe et al. found that GAPDH gene had the target region for 17-β estradiol binding on rat brain due to the rapid and significant increases in GAPDH activity by addition of 10 nM estradiol [35]. This study might explain the response in FSHD myoblast cells. To reveal the effect clearly, by using β-actin and GAPDH together, we showed that in FSHD myoblasts, estradiol treatment clearly changed GAPDH protein levels. Because of this, following experiments had been performed by using β-actin as internal control instead of GAPDH. For future FSHD studies concerning estradiol treatment, we recommend not to use GAPDH as housekeeping protein.

Limitations of the study and conclusion

With this novel study, it was revealed that estradiol could protect myoblasts via regulating DUX4, and estradiol can significantly modify β-catenin and PAX3/PAX7 levels in FSHD derived primary myoblasts. There were some limitations of present study: (i) primary myoblast cells are quite difficult to handle, especially when affection of a muscle disease would be considered. Additionally, there are limited period of passages to perform experiments before the differentiation of myoblasts. (ii) the levels of DUX4 protein is too low to detect especially in protein structure; (iii) neither β-catenin protein nor PAX3/PAX7 protein variants have not been studied with estradiol treatment in FSHD myoblasts, and there were not any literature finding to compare our results during experimental period.

Even though an uphill struggle, we thought the best way to reflect our hypothesis would be using FSHD myoblasts and identified key transcription factors. We observed different responses at 30 min and 4 h with estradiol, because of this we concluded that the temporal timing is very important to design future studies including multiple time sections rather than only one temporally design. Otherwise, it would be impossible to capture zero level of DUX4 and stimulation of 56 kDa PAX3/PAX7 expression.

Level changes of all these transcription factors from different patients’ cell lines at 30 min and 4 h after estradiol had been summarized in Table 2. In this table, each patient’s myoblast exhibited completely divergent DUX4, β-catenin and PAX3/PAX7 responses with estradiol treatment. Knowing that each of these transcription factors have multiple targets, marks of distinct molecular signatures among FSHD patients will be highlighted better. These divergent molecular responses bear considerable resemblance to clinical variability of FSHD patients and that indicates a unique molecular underlying structure for each patient. Supporting this, in literature, studies on FSHD revealed that methylation status of D4Z4 is shown to be different among patients [36] in addition to D4Z4 contraction. There are also multiple other pathways are corrupted as mitochondrial disfunction, oxidative stress [37] which can be driven by multiple other genetic and epigenetic factors. Conclusion here is that, it is necessary to evaluate FSHD clinical features and pathophysiology in combination of genetic, individual and environmental factors. As information is growing up, FSHD shifts to multifactorial disease rather than being a sole genetic disease.

Table 2:

Summary of the affection (increase: ▲/decrease: ▼/no change: (–)/no expression: 0) of selected proteins in all cell lines after estradiol treatment.

	DUX4 (45 kDa)			β-catenin (92 kDa)			PAX3/PAX7 (56 kDa)			PAX3 (45 kDa)			GAPDH (31.5 kDa)
	71yM	63yM	47yF	71yM	63yM	47yF	71yM	63yM	47yF	71yM	63yM	47yF	71yM
C-30 min	–	0	0	▲	–	▲	–	0	0	▼	–	–	–
C-4 h	▼	0	0	▲	▼	▲	▲	0	0	▼	▲	▼	▼

c, control; 30 min, 30 minutes; 4 h, 4 hours.

Our major purpose in investigating for the role of estrogen in female and male FSHD patients was to contribute to the FSHD treatment for future generations. Estrogen is a candidate treatment option in FSHD. However, known potential risks as cancer and thromboembolism make estrogen treatment questionable. In our previous clinical study on FSHD, we revealed approximately 12.5% of our FSHD patients had high fibrinogen levels [38]. That is why, it would be more logical to use local therapies or softer estrogen analogs. Supporting this Banerji et al. recently revealed that phytoestrogens as genistein rescued the impaired muscle phenotype in FSHD [39] , which could be one of the potential therapeutic options for alleviating FSHD soon.

Corresponding author: Sibel Berker Karauzum, PhD, Department of Medical Biology and Genetics, Faculty of Medicine, Akdeniz University, Antalya, 07058, Turkey, Phone: +0090 532 727 89 41, E-mail: sibelberkerkarauzum@gmail.com

Funding source: Akdeniz University Research Foundation

Award Identifier / Grant number: TDK-2016-1460

Acknowledgments

Thanks to Professor Rabi Tawil and Don Henderson for sending cell lines for our doctorate research project only for open a road to FSHD research. Thanks to Judit Balog and Professor Silvere van der Maarel for helping to contact with Rabi Tawil. Thanks to Z. Cetin for his support in reading and editing support before publication. This study has been presented as a poster: at 26th Annual FSHD International Research Congress, Marseille, France (19–20 July 2019). Orally presented at: 13th National Medical Genetic Congress (November 7th–11th 2018 Antalya) and 54th National Neurology Congress (November 30th 2018–December 6th 2019, Antalya).

Research funding: This study was funded by Akdeniz University Research Foundation, Project number: TDK-2016-1460.
Author contributions: Doctorate thesis of CH. Conceptualization CH, EGC, HK, OE, HU, SBK; methodology, CH, EGC, HK, OE, SBK; software, CH, EGC.; validation, CH, OE, HU, SBK; investigation, CH, EGC, HK, OE, SBK; resources CH, SBK.; data curation, CH, SBK, OE; writing—original draft preparation, CH, SBK.; writing—review and editing, CH, EGC, HK, SBK, OE, HU.; supervision, SBK, OE; project administration, SBK.; funding acquisition, SBK. All authors have read and agreed to the published version of the manuscript.
Competing interests: The authors declare no conflict of interest.
Informed consent: Human biopsy samples were obtained from Rochester University Medical Center; informed consent for biopsy samples was obtained from all individual participants.
Ethical approval: All procedures performed in this study were in accordance with the ethical standards of Akdeniz University Institutional Ethical Committee (Decision Number: 284, Date: 28.10.2015) and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

References

1. Hewitt, JE, Lyle, R, Clark, LN, Valleley, EM, Wright, TJ, Wijmenga, C, et al.. Analysis of the tandem repeat locus D4Z4 associated with facioscapulohumeral muscular dystrophy. Hum Mol Genet 1994;3:1287–95. https://doi.org/10.1093/hmg/3.8.1287.Search in Google Scholar

2. Lemmers, RJ, Wohlgemuth, M, Frants, RR, Padberg, GW, Morava, E, van der Maarel, SM. Contractions of D4Z4 on 4qB subtelomeres do not cause facioscapulohumeral muscular dystrophy. Am J Hum Genet 2004;75:1124–30. https://doi.org/10.1086/426035.Search in Google Scholar

3. Lemmers, RJ, Tawil, R, Petek, LM, Balog, J, Block, GJ, Santen, GW, 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–4. https://doi.org/10.1038/ng.2454.Search in Google Scholar

4. Lyle, R, Wright, TJ, Clark, LN, Hewitt, JE. The FSHD-associated repeat, D4Z4, is a member of a dispersed family of homeobox-containing repeats, subsets of which are clustered on the short arms of the acrocentric chromosomes. Genomics 1995;28:389–97. https://doi.org/10.1006/geno.1995.1166.Search in Google Scholar

5. Bosnakovski, D, Toso, EA, Hartweck, LM, Magli, A, Lee, HA, Thompson, ER, et al.. The DUX4 homeodomains mediate inhibition of myogenesis and are functionally exchangeable with the Pax7 homeodomain. J Cell Sci 2017;130:3685–97. https://doi.org/10.1242/jcs.205427.Search in Google Scholar

6. Banerji, CRS, Panamarova, M, Hebaishi, H, White, RB, Relaix, F, Severini, S, et al.. PAX7 target genes are globally repressed in facioscapulohumeral muscular dystrophy skeletal muscle. Nat Commun 2017;8:2152. https://doi.org/10.1038/s41467-017-01200-4.Search in Google Scholar

7. Banerji, CRS, Zammit, PS. 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–36. https://doi.org/10.1093/hmg/ddz043.Search in Google Scholar

8. Osborne, RJ, Welle, S, Venance, SL, Thornton, CA, Tawil, R. Expression profile of FSHD supports a link between retinal vasculopathy and muscular dystrophy. Neurology 2007;68:569–77. https://doi.org/10.1212/01.wnl.0000251269.31442.d9.Search in Google Scholar

9. Cheli, S, Francois, S, Bodega, B, Ferrari, F, Tenedini, E, Roncaglia, E, et al.. Expression profiling of FSHD-1 and FSHD-2 cells during myogenic differentiation evidences common and distinctive gene dysregulation patterns. PLoS One 2011;6:e20966. https://doi.org/10.1371/journal.pone.0020966.Search in Google Scholar

10. Banerji, CR, Knopp, P, Moyle, LA, Severini, S, Orrell, RW, Teschendorff, AE, et al.. Beta-catenin is central to DUX4-driven network rewiring in facioscapulohumeral muscular dystrophy. J R Soc Interface 2015;12:20140797. https://doi.org/10.1098/rsif.2014.0797.Search in Google Scholar

11. Block, GJ, Narayanan, D, Amell, AM, Petek, LM, Davidson, KC, Bird, TD, et al.. Wnt/beta-catenin signaling suppresses DUX4 expression and prevents apoptosis of FSHD muscle cells. Hum Mol Genet 2013;22:4661–72. https://doi.org/10.1093/hmg/ddt314.Search in Google Scholar

12. Giacoppo, S, Thangavelu, SR, Calabro, RS, Bramanti, P, Emanuela, M. Role of Wnt/Β-catenin signaling in regulating the balance between oxidative stress and apoptosis in peripheral blood mononuclear cells from patients with facioscapulohumeral muscular dystrophy. J Cell Mol Genet Med 2016;10. https://doi.org/10.4172/1747-0862.1000213.Search in Google Scholar

13. Tawil, R. Facioscapulohumeral muscular dystrophy. Neurotherapeutics 2008;5:601–6. https://doi.org/10.1016/j.nurt.2008.07.005.Search in Google Scholar

14. Zatz, M, Marie, SK, Cerqueira, A, Vainzof, M, Pavanello, RC, Passos-Bueno, MR. The facioscapulohumeral muscular dystrophy (FSHD1) gene affects males more severely and more frequently than females. Am J Med Genet 1998;77:155–61. https://doi.org/10.1002/(sici)1096-8628(19980501)77:2<155::aid-ajmg9>3.0.co;2-r.10.1002/(SICI)1096-8628(19980501)77:2<155::AID-AJMG9>3.0.CO;2-RSearch in Google Scholar

15. Ciafaloni, E, Pressman, EK, Loi, AM, Smirnow, AM, Guntrum, DJ, Dilek, N, et al.. Pregnancy and birth outcomes in women with facioscapulohumeral muscular dystrophy. Neurology 2006;67:1887–9. https://doi.org/10.1212/01.wnl.0000244471.05316.19.Search in Google Scholar

16. Awater, C, Zerres, K, Rudnik-Schoneborn, S. Pregnancy course and outcome in women with hereditary neuromuscular disorders: comparison of obstetric risks in 178 patients. Eur J Obstet Gynecol Reprod Biol 2012;162:153–9. https://doi.org/10.1016/j.ejogrb.2012.02.020.Search in Google Scholar

17. Mul, K, Horlings, CGC, Voermans, NC, Schreuder, THA, van Engelen, BGM. Lifetime endogenous estrogen exposure and disease severity in female patients with facioscapulohumeral muscular dystrophy. Neuromuscul Disord 2018;28:508–11. https://doi.org/10.1016/j.nmd.2018.02.012.Search in Google Scholar

18. Hangul, C, Bozkurt, S, Bilge, U, Ozdem, S, Altunbas, H, Uysal, H, et al.. The ratios of estradiol and progesterone to testosterone influence the severity of facioscapulohumeral muscular dystrophy. Neurol Sci Neurophysiol 2020;37:190–6.10.4103/NSN.NSN_37_20Search in Google Scholar

19. Puma, A, Garibaldi, M, Teveroni, E, Deidda, G, Moretti, F, Sacconi, S. Estrogens as a potential disease modifier in FSHD: a retrospective clinical study. Neuromuscul Disord 2017;27:S200. https://doi.org/10.1016/j.nmd.2017.06.386.Search in Google Scholar

20. Teveroni, E, Pellegrino, M, Sacconi, S, Calandra, P, Cascino, I, Farioli-Vecchioli, S, et al.. Estrogens enhance myoblast differentiation in facioscapulohumeral muscular dystrophy by antagonizing DUX4 activity. J Clin Invest 2017;127:1531–45. https://doi.org/10.1172/jci89401.Search in Google Scholar

21. Collins, BC, Arpke, RW, Larson, AA, Baumann, CW, Xie, N, Cabelka, CA, et al.. Estrogen regulates the satellite cell compartment in females. Cell Rep 2019;28:368–81. e366. https://doi.org/10.1016/j.celrep.2019.06.025.Search in Google Scholar

22. Ricci, E, Galluzzi, G, Deidda, G, Cacurri, S, Colantoni, L, Merico, B, et al.. Progress in the molecular diagnosis of facioscapulohumeral muscular dystrophy and correlation between the number of KpnI repeats at the 4q35 locus and clinical phenotype. Ann Neurol 1999;45:751–7. https://doi.org/10.1002/1531-8249(199906)45:6<751::aid-ana9>3.0.co;2-m.10.1002/1531-8249(199906)45:6<751::AID-ANA9>3.0.CO;2-MSearch in Google Scholar

23. van Overveld, PG, Enthoven, L, Ricci, E, Rossi, M, Felicetti, L, Jeanpierre, M, et al.. Variable hypomethylation of D4Z4 in facioscapulohumeral muscular dystrophy. Ann Neurol 2005;58:569–76. https://doi.org/10.1002/ana.20625.Search in Google Scholar

24. Li, Y, Wen, Y, Green, M, Cabral, EK, Wani, P, Zhang, F, et al.. Cell sex affects extracellular matrix protein expression and proliferation of smooth muscle progenitor cells derived from human pluripotent stem cells. Stem Cell Res Ther 2017;8:156. https://doi.org/10.1186/s13287-017-0606-2.Search in Google Scholar

25. La Colla, A, Vasconsuelo, A, Milanesi, L, Pronsato, L. 17beta-estradiol protects skeletal myoblasts from apoptosis through p53, Bcl-2, and FoxO families. J Cell Biochem 2017;118:104–15. https://doi.org/10.1002/jcb.25616.Search in Google Scholar

26. Berio, E, Divari, S, Starvaggi Cucuzza, L, Biolatti, B, Cannizzo, FT. 17beta-estradiol upregulates oxytocin and the oxytocin receptor in C2C12 myotubes. PeerJ 2017;5:e3124. https://doi.org/10.7717/peerj.3124.Search in Google Scholar

27. Bradford, MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem 1976;72:248–54. https://doi.org/10.1016/0003-2697(76)90527-3.Search in Google Scholar

28. Snider, L, Geng, LN, Lemmers, RJ, Kyba, M, Ware, CB, Nelson, AM, et al.. Facioscapulohumeral dystrophy: incomplete suppression of a retrotransposed gene. PLoS Genet 2010;6:e1001181. https://doi.org/10.1371/journal.pgen.1001181.Search in Google Scholar

29. Tassin, A, Laoudj-Chenivesse, D, Vanderplanck, C, Barro, M, Charron, S, Ansseau, E, et al.. DUX4 expression in FSHD muscle cells: how could such a rare protein cause a myopathy? J Cell Mol Med 2013;17:76–89. https://doi.org/10.1111/j.1582-4934.2012.01647.x.Search in Google Scholar

30. Statland, JM, Tawil, R. Facioscapulohumeral muscular dystrophy. Continuum 2016;22:1916–31. https://doi.org/10.1212/con.0000000000000399.Search in Google Scholar

31. Suzuki, A, Scruggs, A, Iwata, J. The temporal specific role of WNT/beta-catenin signaling during myogenesis. J Nat Sci 2015;1:e143.Search in Google Scholar

32. Brzoska, E, Przewozniak, M, Grabowska, I, Janczyk-Ilach, K, Moraczewski, J. Pax3 and Pax7 expression during myoblast differentiation in vitro and fast and slow muscle regeneration in vivo. Cell Biol Int 2009;33:483–92. https://doi.org/10.1016/j.cellbi.2008.11.015.Search in Google Scholar

33. Wigmore, PM, Evans, DJ. Molecular and cellular mechanisms involved in the generation of fiber diversity during myogenesis. Int Rev Cytol 2002;216:175–232. https://doi.org/10.1016/s0074-7696(02)16006-2.Search in Google Scholar

34. Klinge, CM. Estrogen receptor interaction with estrogen response elements. Nucleic Acids Res 2001;29:2905–19. https://doi.org/10.1093/nar/29.14.2905.Search in Google Scholar

35. Joe, I, Ramirez, VD. Binding of estrogen and progesterone-BSA conjugates to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the effects of the free steroids on GAPDH enzyme activity: physiological implications. Steroids 2001;66:529–38. https://doi.org/10.1016/s0039-128x(00)00220-8.Search in Google Scholar

36. Gaillard, MC, Roche, S, Dion, C, Tasmadjian, A, Bouget, G, Salort-Campana, E, et al.. Differential DNA methylation of the D4Z4 repeat in patients with FSHD and asymptomatic carriers. Neurology 2014;83:733–42. https://doi.org/10.1212/wnl.0000000000000708.Search in Google Scholar

37. Denny, AP, Heather, AK. Are antioxidants a potential therapy for FSHD? A review of the literature. Oxid Med Cell Longev 2017;2017:7020295. https://doi.org/10.1155/2017/7020295.Search in Google Scholar

38. Hangul, C, Ozsoy, U, Hizay, A, Bozkurt, S, Bilge, U, Ozdem, S, et al.. Facioscapulohumeral muscular dystrophy (FSHD) has a talk with endocrinologic parameters: estradiol, progesterone and testosterone. 26th annual FSHD international research congress, 2019. Marseille, France: FSHD Society; 2019:011 p. www.fshdsociety.org.Search in Google Scholar

39. Banerji, CRS, Panamarova, M, Pruller, J, Figeac, N, Hebaishi, H, Fidanis, E, et al.. Dynamic transcriptomic analysis reveals suppression of PGC1alpha/ERRalpha drives perturbed myogenesis in facioscapulohumeral muscular dystrophy. Hum Mol Genet 2019;28:1244–59. https://doi.org/10.1093/hmg/ddy405.Search in Google Scholar

Received: 2020-07-11

Accepted: 2020-11-06

Published Online: 2021-01-05

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

https://doi.org/10.1515/tjb-2020-0351

Keywords for this article

β-catenin; DUX4; estradiol; facioscapulohumeral muscular dystrophy; FSHD; PAX3; PAX7; skeletal muscle

Creative Commons

BY 4.0