Next Article in Journal / Special Issue
MyD88 Is Not Required for Muscle Injury-Induced Endochondral Heterotopic Ossification in a Mouse Model of Fibrodysplasia Ossificans Progressiva
Previous Article in Journal
Triazole-Modified Nucleic Acids for the Application in Bioorganic and Medicinal Chemistry
Previous Article in Special Issue
Challenges and Opportunities for Drug Repositioning in Fibrodysplasia Ossificans Progressiva
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 9
Issue 6
10.3390/biomedicines9060629
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients

by
Ton Schoenmaker
1,*,
Michal Mokry
2,
Dimitra Micha
3,
Coen Netelenbos
4,
Nathalie Bravenboer
5,
Marjolijn Gilijamse
6,
E. Marelise W. Eekhoff
4 and
Teun J. de Vries
1
1
Department of Periodontology, Academic Centre for Dentistry Amsterdam (ACTA), University of Amsterdam and Vrije Universiteit, 1012 Amsterdam, The Netherlands
2
Division Laboratory, Pharmacy and Biomedical Genetics, UMC Utrecht, 3508 Utrecht, The Netherlands
3
Department of Clinical Genetics, Amsterdam Movement Sciences, Amsterdam UMC, Vrije Universiteit, 1012 Amsterdam, The Netherlands
4
Department of Internal Medicine Section Endocrinology, Amsterdam Movement Sciences, Amsterdam UMC, Vrije Universiteit Amsterdam, 1012 Amsterdam, The Netherlands
5
Department of Clinical Chemistry, Amsterdam Movement Sciences, Vrije Universiteit Amsterdam, 1012 Amsterdam, The Netherlands
6
Department of Oral and Maxillofacial Surgery and Oral Pathology, Amsterdam UMC, Vrije Universiteit, 1012 Amsterdam, The Netherlands
*
Author to whom correspondence should be addressed.
Biomedicines 2021, 9(6), 629; https://doi.org/10.3390/biomedicines9060629
Submission received: 16 March 2021 / Revised: 17 May 2021 / Accepted: 28 May 2021 / Published: 1 June 2021
(This article belongs to the Special Issue Fibrodysplasia Ossificans Progressiva: Studies on Disease Mechanism towards Novel Therapeutic Approaches)
Browse Figures
Versions Notes

Abstract

:
Fibrodysplasia Ossificans Progressiva (FOP) is a rare genetic disease characterized by heterotopic ossification (HO). It is caused by mutations in the Activin receptor type 1 (ACVR1) gene, resulting in enhanced responsiveness to ligands, specifically to Activin-A. Though it has been shown that capturing Activin-A protects against heterotopic ossification in animal models, the exact underlying mechanisms at the gene expression level causing ACVR1 R206H-mediated ossifications and progression are thus far unknown. We investigated the early transcriptomic changes induced by Activin-A of healthy control and patient-derived periodontal ligament fibroblasts (PLF) isolated from extracted teeth by RNA sequencing analysis. To study early differences in response to Activin-A, periodontal ligament fibroblasts from six control teeth and from six FOP patient teeth were cultured for 24 h without and with 50 ng/mL Activin-A and analyzed with RNA sequencing. Pathway analysis on genes upregulated by Activin-A in FOP cells showed an association with pathways involved in, among others, Activin, TGFβ, and BMP signaling. Differential gene expression induced by Activin-A was exclusively seen in the FOP cells. Median centered supervised gene expression analysis showed distinct clusters of up- and downregulated genes in the FOP cultures after stimulation with Activin-A. The upregulated genes with high fold changes like SHOC2, TTC1, PAPSS2, DOCK7, and LOX are all associated with bone metabolism. Our open-ended approach to investigating the early effect of Activin-A on gene expression in control and FOP PLF shows that the molecule exclusively induces differential gene expression in FOP cells and not in control cells.

1. Introduction

Fibrodysplasia ossificans progressiva (FOP) is a rare, one-in-two-million-occurring, autosomal dominant genetic disease characterized by progressive heterotopic bone formation (HO) where especially muscles, tendons, and ligaments are converted into bone [1,2,3]. The clinical manifestation of the heterotopic ossification is rather diverse. HO can occur after a flare up, during inflammation, after injury, or even spontaneously. In some patients, HO is more progressive than in others, and episodes exist with complete absence of HO formation [4]. Because of these differences in clinical manifestations, cell biological approaches that can shed light on the biochemical events that precede heterotopic ossification are of great importance as a prelude to therapy. Over the past decade, our understanding of the pathogenesis of the disease has improved considerably. The causative mutation in the gene encoding the Activin Receptor Type I receptor (ACVR1) was identified by the group of Shore et al. [5]. The most frequent R206H mutation changes an arginine to a histidine in the glycine–serine rich cytoplasmic domain of this bone morphogenetic protein (BMP) type I receptor [2,5]. This amino acid change makes the receptor more sensitive to BMP signaling and simultaneously results in decreased binding of the receptor’s natural inhibitor, FKBP12, likely resulting in leaky signaling of the receptor through SMADs, the main signal transducers of the transforming growth factor-beta (TGF-β) family [6,7,8,9,10]. More recently, the discovery by two independent groups that Activin-A can also function as an activator of the mutated ACVR1 [11,12] has given great impulse to the FOP research community, since now it is conceivable to specifically prevent heterotopic bone formation by inhibiting Activin-A. Under normal circumstances, Activin-A signals mainly via ACVR1B receptor complexes through SMAD 2/3 phosphorylation [13,14,15]. Upon binding to ACVR1, it normally inhibits SMAD 1/5/9 signaling and subsequent osteogenic differentiation [16], but, as mentioned before, when Activin-A binds to the mutated ACVR1, it activates SMAD1/5/9 phosphorylation and subsequent osteogenic differentiation [11,12]. This was demonstrated by Hatsell et al. in their FOP mouse model, but also in the currently conducted LUMINA-I trial investigating the Activin-A blocking antibody Garetosmab, which has revealed promising results in the prevention of newly formed HO [17,18]. Even though it has by now been clearly shown that capturing Activin-A protects against heterotopic ossification, this knowledge is mainly limited to mouse models for FOP. Moreover, the exact underlying mechanism at the gene expression level causing ACVR1 R206H-mediated ossification and progression of the disease is unknown. Thus, knowledge of the downstream transcriptomic changes after Activin-A-ACVR1 R206H binding using open-ended approaches such as RNA sequencing is warranted.
One of the challenges that FOP research faces is the lack of primary cells from patients, since surgical removal of the heterotopic bone, as a potential source of bone-forming osteoblasts, often results in the appearance of new heterotopic bone formation. Several alternative sources for bone-forming cells are currently being investigated as a tool for the unraveling of the FOP HO mechanism. Dermal fibroblasts have been demonstrated by Micha et al. to be a source of cells with osteoblastic properties. They showed that the FOP-derived fibroblasts have a high osteogenic potential [19]. Another source of human primary cells are cells isolated from extracted teeth. Many FOP patients develop a locked jaw, making oral care extremely difficult. In some cases, tooth extraction is the only dental treatment possible. Additionally, teeth are sometimes extracted to provide extra space in the mouth. There have been no reports of HO in the sockets at the extraction sites. Some groups have made use of the cells from the discarded primary teeth (SHED cells) [7,20], showing higher osteogenic differentiation in the FOP cells compared to control cells. Our group has recently shown that periodontal ligament cells obtained from FOP patients after tooth extraction show both osteogenic- as well as osteoclastogenic-inducing capacities, as do PLF from healthy individuals [21,22,23,24]. Here, we used these cells to investigate the Activin-A-induced transcriptome differences between controls and FOP patients. This study explored for the first time RNA sequencing as an open-ended approach for non-biased identification of potential new differences in early gene expression between control and FOP fibroblasts, directly following Activin-A binding. Given the unique signaling by Activin-A via the mutated ACVR-1, as described in mouse models [11] and iPS cell approaches [12], we hypothesize that Activin-A will induce early transcriptomic differences specifically in primary cells from FOP patients.

2. Materials and Methods

2.1. Periodontal Ligament Fibroblasts

Periodontal ligament cells were retrieved from 6 extracted teeth from 5 FOP patients (3 female patients age 26, 28, and 46 years and 2 male patients, age 21 and 69 years) and 6 healthy control donors (4 females age 19–26 years, 2 males, age 23–27 years). All FOP patients carried the classical R206H mutation. Written informed consent was obtained from each participant, and researchers were not able to trace the origin of the material to any individual, as required by Dutch law. Extractions were performed at the Department of Oral and Maxillofacial Surgery and Oral Pathology, Amsterdam UMC, Location VUmc. No differences in bone healing between the FOP patients and the control group following tooth removal were reported. Periodontal ligament cells were isolated, as previously described [21,25,26]. Shortly, the periodontal ligament was scraped off the middle one third of the root, and cells were allowed to grow out of the tissue by culturing them in culture medium, which consisted of Dulbecco’s modified eagle medium (DMEM, Gibco BRL, Paisley, Scotland) supplemented with 10% FCS (HyClone, Logan, UT, USA) and 1% antibiotics: 100 U/mL penicillin, 100 µg/mL streptomycin, and 250 ng/mL amphotericin B (Sigma, St. Louis, MO, USA). Cells were propagated, and 3rd passage cells were frozen and stored in liquid nitrogen. All experiments were performed with 4th passage cells.

2.2. Cell Culturing and RNA Isolation

Cells were first allowed to attach by overnight culturing in 6-well plates at a density of 3.2 × 104 cells/cm2 (3 × 105 cells/well) in DMEM + 10% FCS + 1% antibiotics. The following day media were replaced with media without or with 50 ng/mL Activin-A (Sigma, St Louis, MO, USA). After 24 h, RNA was isolated using the RNeasy Mini kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. RNA yield and quality were measured with the Synergy multi-mode reader (Biotek, Winooski, VT, USA).

2.3. RNA Sequencing

Sample and library preparation was performed using the NEXTflex Rapid Directional mRNA-Seq kit (Bioo Scientific, Austin, TX, USA). Single-end 75 bp RNA sequencing was performed on the Illumina NextSeq500 (Illumina, San Diego, CA, USA) at the Wilhelmina Children’s Hospital in Utrecht. Reads were aligned to the human reference genome GRCh37 using STAR version 2.4.2a. Picard’s AddOrReplaceReadGroups (v1.98) was used to add read groups to the BAM files, which were sorted with Sambamba v0.4.5, and transcript abundances were quantified with HTSeq-count version 0.6.1p1 using the union mode. Subsequently, reads per kilobase million reads sequenced (RPKMs) were calculated with edgeR’s RPKM function. Differentially expressed genes were identified using the DESeq2 package with standard settings. In the first analysis, the genes with a non-false discovery rate corrected p value of <0.01 were used for the pathway and GO analysis. To finally assess on the single gene level which genes were statistically significant differentially expressed, we performed a false discovery correction of 10%, generating p adjusted (p(adj)) values. For these p(adj) values, we used a cutoff of 0.1 (p(adj) < 0.1) to determine significance.

2.4. Pathway and Gene Enrichment Analysis

Pathway analysis and functional enrichment analysis were performed on the differentially expressed genes with p < 0.01, using the GeneMania application version 2017-03-14 at www.genemania.org [27,28,29] using standard settings. The program builds a network based on the interactions between the input genes (e.g., genes upregulated by Avtivin-A in the FOP cells) and adds genes that are relevant to that network. The pathway maps shown in Figure 2 are an adjustment to the enrichment map consolidated pathways described by Merico et al. [30]. Boxed genes are genes from the input list; the other genes are relevant genes added by GeneMania. Size of the nodes is dependent on the number of interactions that gene has with the network. The gene enrichment analysis was performed with the ToppGene suite application ToppFun version September 2020 [31,32] using standard settings.

3. Results

3.1. Different Donors Cluster Together

To compare the transcription profiles of the different donors, we first performed non-supervised Principal Component Analysis (PCA) of all samples that passed the QC filters. Two control libraries and one FOP (Control 2 + Activin-A, Control 5–Activin-A and FOP 5–Activin-A) library did not pass the Quality Control filter, and therefore these samples were excluded from further statistical analysis. Figure 1A depicts the inter-donor variability between the samples. Figure 1B depicts the variability induced by the experimental conditions (e.g., without or with Activin-A). This analysis shows that the two culture conditions of one donor tend to cluster together, indicating that the variability between donors is higher than the variability induced by the Activin-A treatment.

3.2. Activin-A Activates Distinct Pathways in FOP Cells

We next investigated which genes were influenced by Activin-A in the PLF cells. Using the non-FDR corrected gene expression data with p < 0.01, we found that 131 genes were differentially expressed in FOP cells (69 upregulated and 62 downregulated genes) as compared to 46 genes in the control cells (18 upregulated and 28 downregulated genes), see supplementary Tables S1 and S2. Network analysis using the GeneMania application showed that the genes upregulated by Activin-A in the FOP cells show an association with different pathways involved in, among others, TGF-beta signaling (MSigDB M2642) (Figure 2A), BMP receptor signaling (MsigDB M181) (Figure 2B), and signaling by Activin (Mysid M26965) (Figure 2C) showing Activin-A indeed induces its known signaling pathways also in this primary cell system. Additionally, functional enrichment analysis performed in GeneMania showed gene ontology terms enriched among the upregulated genes. These GO terms are involved in Activin, TGF-beta, and BMP signaling; the top 10 of the associated GO terms are listed in Table 1. None of these pathways or GO terms were found to be associated to Activin-A stimulation in the control cells (data not shown), indicating that Activin-A signals distinctly differently in the FOP cells as compared to in the control cells. Gene enrichment analysis on the FOP cells using ToppGene Suite showed that the differentially expressed genes in biological domain molecular functio, are mainly associated with gene ontology (GO) terms linked to cytoskeletal and actin binding (Table 2). In the biological domain biological process, the GO terms are linked to cell–cell and cell–substrate junctions and responses to injury and endogenous stimuli (Table 3). In the cellular component domain, the genes are associated with focal adhesions and cell–cell and cell–substrate junctions (Table 4). The observed association with the collagen gene family is possibly also related to bone metabolism (Table 5). None of these associations were found in the control cells.

3.3. Activin-A Induces Differential Gene Expression in FOP Cells

Subsequent analysis on individual gene expression level was performed only on samples with both libraries, without and with Activin-A, per donor and corrected for donor variances in the final model. A false discovery correction of 10% was performed generating p(adj) values. Under these criteria, the final analyses were performed on 4 control samples and 5 FOP samples. To graphically display the differential gene expression induced by Activin-A, MA plots were generated where the Log2 fold change (M) was plotted against the average expression across the samples (A) using a p-adjusted cut-off value of 0.1. No differentially expressed genes were observed in control samples (Figure 3A). Interestingly, using these same settings, differential gene expression was observed in the FOP samples, shown in the MA plot as red dots (Figure 3B).

3.4. Differential Gene Expression Is Induced by Activin-A in FOP Cells

Having demonstrated that Activin-A alters gene expression only in FOP cells, we next investigated which specific genes were differentially expressed under the influence of Activin-A in these cells. Median-centered supervised gene expression analysis showed distinct clusters of up- and downregulated genes in the untreated (−ActA) and Activin-A treated (+ActA) groups (Figure 3C). From this analysis, using a false detection rate (FDR) of 10%, a list of the highest up- and downregulated genes was constructed as shown in Table 6. The genes that are more than two-fold upregulated by Activin-A in the FOP cells are somehow involved in bone metabolism. For instance, SHOC2 is a positive regulator of the MAP/ERK pathway, possibly involved in pathologies with bone and skeletal defects [33,34]. TTC1 activates RAS signaling, which regulates osteoprogenitor cell proliferation [35,36]. PAPSS2 is needed for sulfation of extracellular matrix molecules during bone development; deficiencies result in osteochondrodysplasias [37,38]. DOCK7 is a guanidine nucleotide exchange factor. Misty mice, with a loss of function mutation in this gene, show uncoupled bone remodeling and reduced bone formation probably linked to reduced brown fat production [39,40]. LOX is involved in the crosslinking of collagen and elastin and knock-out mice show decreased osteoblast differentiation [41,42]. A link to more direct Activin-A signaling can be found in RAB27B, which is involved in pituitary hormone secretion like FSH, the secretion of which is stimulated by Activin-A [43,44].
Except for COL6A3, the 1.3- to 3.9-fold downregulated genes seem to be less clearly associated with muscle and bone metabolism. COL6A3, which encodes for the alpha3 chain of collagen typeVI, is associated with the extracellular matrix of skeletal muscle, skin, and cartilage [45,46]. PLEC is an important molecule in muscle fibers [47]. CD9 is involved in cell fusion processes, as seen in osteoclast formation and muscle cell fusion [48,49]. TGFBRAP1, finally, is a Smad4 chaperone that associates with inactive TGFβ and Activin receptor complexes [50].

4. Discussion

The discovery of the causative mutation in the ACVR1 gene resulting in the altered responsiveness of this BMP receptor [5] has paved the way to investigating the exact mechanism that underlies the heterotopic ossification seen in FOP. Especially during the last decade, great progress has been made, and the discovery of Activin-A as an activator of the mutated ACVR1 receptor [11,12] has stimulated the scientific FOP community to explore new therapeutic avenues. Here, we investigated the effect on early transcriptome differences induced by Activin-A in our recently established human primary FOP cell culturing system using periodontal ligament fibroblasts [21].
Our data show that indeed the mutated ACVR1 receptor has an altered responsiveness to Activin-A as compared to the non-mutated receptor. Initial gene expression analysis without false discovery rate correction showed that, when using a p-value cutoff of p < 0.01, 131 genes are differentially expressed in FOP cells as compared to 46 genes in the control cells. Pathway analysis of these genes using GeneMania shows an association with pathways involved in Activin, BMP, and TGFβ signaling only in the FOP cells. FOP is an autosomal dominant genetic disease, meaning that both the non-mutated as well as the mutated gene is being expressed, which has been demonstrated in FOP derived SHED cells and monocytes [22,51]. The association with the Activin and TGFβ signaling pathways indicates signaling via the non-mutated receptor, whereas the BMP signaling via SMAD1, 5, and 9 shows that the altered responsiveness of the mutated receptor to Activin-A as described by other groups [20,52] is also present in the periodontal ligament fibroblasts from FOP patients. Additionally, the gene enrichment analysis on these genes show a stronger reaction of the mutated receptor when exposed to Activin-A. There is an association with different GO terms, whereas no such association could be found in the control cells. HO in FOP often occurs after inflammation or injury [3], a process that is possibly mediated by an increase in Activin-A [3,11]. The altered gene expression profile seen here in the FOP cells after short exposure to Activin-A shows an association with GO terms GO 0009611 (response to wounding), GO:0009719 (response to endogenous stimulus), and GO:0071495 (cellular response to endogenous stimulus). This probably reflects the altered reaction described in FOP after injury or infection [53,54,55]. Under normal circumstances when Activin-A binds to the non-mutated ACVR1 receptor, it acts as an inhibitor of the canonical BMP-mediated signaling via ACVR1 and subsequent osteogenic differentiation. In FOP, however, the mutation results in activation of the mutated ACVR1 upon Activin-A binding, resulting in enhanced osteogenic differentiation. The association with the Activin-A-induced upregulated genes in the FOP cells with the collagen gene family observed in this study probably reflects the osteogenesis process taking place as a response to this altered signaling of the mutated ACVR1.
The inter-donor variability has a higher impact on the difference in overall expression compared to the experimental condition (e.g., without or with Activin-A), as described in Figure 1, suggesting that differences on the single-gene-expression level induced by Activin-A treatment could be difficult to detect above the “noise” induced by the above-mentioned inter-donor variability. Surprisingly, despite this fact, supervised gene expression analysis using a false discovery rate of 10% showed that differential expression was seen, but only in the FOP cells. Some of the genes that showed a higher than two-fold upregulation in the FOP cells under the influence of Activin-A can be related to bone metabolism or heterotopic ossification and can thus be regarded as potential targets to inhibit HO in FOP. These genes will be discussed next in the context of the process of endochondral bone formation.
Heterotopic ossification, as seen in FOP, is believed to be formed via endochondral bone formation where chondrogenic differentiation of stem cells results in cartilage formation that is subsequently converted to bone. We find an increase in DOCK7 expression in the FOP cells, possibly resulting in increased chondrogenesis and osteogenesis, ultimately resulting in enhanced HO. The mutated ACVR1 has been shown to play an important role in enhanced chondrogenic differentiation in a gain-of-function mouse model [56]. In mouse models of HO, one of the early steps in the formation of HO is the biogenesis of brown adipose tissue (BAT) [57]. This BAT creates a hypoxic environment that is suggested to induce chondrogenic differentiation of mesenchymal stem cells [58] and subsequent endochondral bone formation. Misty mice lack BAT and show decreased osteogenesis and increased osteoclast activity [40]. This uncoupling of bone formation and bone resorption is the result of a loss of function mutation in the DOCK7 gene [39,59]. Additionally, the 2.3-fold upregulation of PAPPS2 observed in the FOP cells here may relate to the increased early stage endochondral ossification. The cartilage matrix that is produced during the early stages of endochondral ossification contains high concentrations of sulfate groups bound to several extracellular matrix proteins. PAPPS2 is a phosphosulphate phosphatase that catalyzes the sulfation of these extracellular matrix proteins [38]. Expression of PAPPS2 is partly regulated by TGFβ [60]. Mutations in PAPPS2 result in different types of osteochondrodysplasias, examples of which are described by Wang et al. [37]. The upregulation of TTC1 expression suggests that the RAS signaling pathway, which belongs to the osteogenic differentiation program, is activated. Osteogenic differentiation may take place via different pathways. One of the signaling pathways involved in osteogenic differentiation is the RAS signaling [61,62]. TTC1 is an adaptor protein that binds to Gα proteins, resulting in an interaction with small GTPase RAS that activates the RAS signaling pathway [35,36]. Similarly, the higher expression of SHOC2 seen here is an implication of increased RAS-MAPK signaling. SHOC2 is a leucine-rich repeat scaffolding protein that is involved in the RAS-MAPK signaling [33,34]. Mutations in this gene are linked to RASopathies like Mazzaniti syndrome and Noonan syndrome [63,64], both of which result in reduced postnatal growth. At a later stage of endochondral bone formation, the formation of collagen cross links is an important step that determines the integrity and strength of the bone. One of the enzymes involved in the formation of these crosslinks is Lysil Oxidase (LOX), a copper-dependent enzyme that catalyzes the deamination of lysine and hydroxylisine residues, thereby initiating the crosslinking of collagen molecules [41,42,65,66]. LOX-deficient mice show lower osteoblast differentiation and activity [42], and lower numbers of crosslinks are related to bone fragility in, for instance, aging and osteoporosis [65]. Interestingly, expression of LOX is partly regulated by the hypoxia-induced transcription factor HIF1 [67,68]. HIF1 is, together with Activin-A, upregulated during injury [20,53], providing potential links between hypoxia, FOP, and LOX expression.
Taken together, the upregulated genes described above play a role in either the early, chondrogenic, stages of endochondral bone formation (DOCK7 and PAPPS2) or the later, more osteogenic, stages of this process (TTC1, SHOC2. and LOX). This indicates that the altered responsiveness of the mutated ACVR1 gene to Activin-A could induce heterotopic bone formation. Further studies to elucidate the exact role of athe bove-described genes in HO in FOP cells need to be performed. In these studies, blocking experiments with, for instance, follistatin will be included. We and other groups have previously shown that follistatin indeed inhibits the effect of Activin-A on other FOP derived primary cells [20,22]. We did not include such blocking experiments in this RNA sequencing experiment, because we first wanted to investigate which genes were influenced in the first place. Additionally, the subsequent data and statistical analysis where we would compare four different experimental conditions would become rather complex.
As discussed above, the up- and downregulated genes were identified only in FOP patient-derived cells. Nevertheless, and strikingly, the fold-change was relatively small, between 1.3- and 2.5-fold changes in the upregulated genes and 1.3- and 3.9-fold changes in the downregulated group. This suggests that the advanced heterotopic ossification in FOP could be the result of additive gene expression effects, rather than being caused by one dominantly over- or under expressed gene. QPCR experiments did not confirm the statistical differences, as seen in the RNAsequencing data analysis. This could be due to the fact that, as stated in Figure 1, the inter-donor variability is even higher than the differences induced by the experimental condition (e.g., addition of Activin-A), but probably more importantly, the fact that the average expression differences induced by Activin-A were below two-fold. In our experience, it is hard to reach statistical differences in QPCR experiments on human material with a low n as used here. In this paper, we use the periodontal ligament fibroblast as a cell biological tool to study early effects of Activin-A on osteogenic differentiation. Different groups have previously shown the osteoblast-like properties of these cells, making them a relatively accessible source of cells to study osteogenic differentiation in, for instance, FOP research [69,70,71,72]. However, as previously stated, following tooth extraction, no HO formation in the socket has been reported. This may imply that the periodontal ligament fibroblasts exhibit HO formation inhibiting properties, which may be the underlying reason for the relatively small fold changes witnessed, though an alternative possibility is the short incubation time with Activin-A. Longer exposure to the molecule might reveal more and possibly stronger effects on the transcriptome, resulting, when using osteogenic medium, in Activin-A-enhanced osteogenesis, as recently shown with tooth-associated cells by Wang et al. [20]. Gene enrichment and pathway analysis did not identify many or any new pathways associated to Activin-A signaling through the mutated receptor. This may be due to the relatively small number of differentially expressed genes, which does not allow for extensive pathway analysis. Nevertheless, when applying supervised gene expression analysis with a false discovery rate of 10%, our data show that Activin-A only induces differential gene expression in the FOP cells bearing the R206H mutation. The highest upregulated genes have all been linked to bone metabolism, and some even to heterotopic ossification in general, but none of these genes have, to our knowledge, previously been described to be linked to HO in FOP. The differential expression of genes somehow involved in bone formation by Activin-A exclusively in cells carrying the FOP mutation is in strong support of the therapeutic treatment rationale of inhibiting Activin-A. This study adds evidence to the notion that when disarming Activin-A’s malevolent effect, heterotopic bone formation can be tempered in the lives of people who suffer the daily anxiety of progressive ectopic bone formation with all its disabling consequences.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biomedicines9060629/s1, Table S1A: Non-FDR corrected Upregulated genes by Activin-A in Control cells; Table S1B: Non-FDR corrected Downregulated genes by Activin-A in Control cells; Table S2A: Non-FDR corrected Upregulated genes by Activin-A in FOP cells; Table S2B: Non-FDR corrected Downregulated genes by Activin-A in FOP cells.

Author Contributions

T.S. designed and conducted the study, performed the data analysis, and wrote the manuscript. M.M. performed the RNA sequencing and initial data analysis. D.M., C.N. and N.B. contributed to the design of the study and writing the manuscript. M.G. performed the extraction of the teeth. E.M.W.E. and T.J.d.V. supervised the study and contributed to the study design and writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the IFOPA FOP Competitive Research Grant.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Free University Amsterdam (protocol code 2012.467).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study, and researchers were not able to trace the origin of the material to any individual, as required by Dutch law.

Data Availability Statement

The RNAseq data presented in this study are available upon reasonable request. The data are not publicly available, because EU regulations do not allow us to put genetic or private data on the public domain (like GEO).

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Bravenboer, N.; Micha, D.; Triffit, J.T.; Bullock, A.N.; Ravazollo, R.; Bocciardi, R.; Di Rocco, M.; Netelenbos, J.C.; Dijke, P.T.; Sánchez-Duffhues, G.; et al. Clinical Utility Gene Card for: Fibrodysplasia ossificans progressiva. Eur. J. Hum. Genet. 2015, 23, 1431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Pignolo, R.J.; Shore, E.M.; Kaplan, F.S. Fibrodysplasia Ossificans Progressiva: Clinical and Genetic Aspects. Orphanet J. Rare Dis. 2011, 6, 80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kaplan, F.S.; Le Merrer, M.; Glaser, D.L.; Pignolo, R.J.; Goldsby, R.E.; Kitterman, J.A.; Groppe, J.; Shore, E.M. Fibrodysplasia ossificans progressiva. Best Pract. Res. Clin. Rheumatol. 2008, 22, 191–205. [Google Scholar] [CrossRef] [PubMed]
  4. Botman, E.; Teunissen, B.P.; Raijmakers, P.; De Graaf, P.; Yaqub, M.; Treurniet, S.; Schoenmaker, T.; Bravenboer, N.; Micha, D.; Pals, G.; et al. Diagnostic Value of Magnetic Resonance Imaging in Fibrodysplasia Ossificans Progressiva. JBMR Plus 2020, 4, e10363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Shore, E.M.; Xu, M.; Feldman, G.J.; A Fenstermacher, D.; Cho, T.-J.; Choi, I.H.; Connor, J.M.; Delai, P.; Glaser, D.L.; LeMerrer, M.; et al. A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva. Nat. Genet. 2006, 38, 525–527. [Google Scholar] [CrossRef] [PubMed]
  6. Song, G.-A.; Kim, H.-J.; Woo, K.M.; Baek, J.-H.; Kim, G.-S.; Choi, J.-Y.; Ryoo, H.-M. Molecular Consequences of the ACVR1R206H Mutation of Fibrodysplasia Ossificans Progressiva. J. Biol. Chem. 2010, 285, 22542–22553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Billings, P.C.; Fiori, J.L.; Bentwood, J.L.; O’Connell, M.P.; Jiao, X.; Nussbaum, B.; Caron, R.J.; Shore, E.M.; Kaplan, F.S. Dysregulated BMP Signaling and Enhanced Osteogenic Differentiation of Connective Tissue Progenitor Cells from Patients with Fibrodysplasia Ossificans Progressiva (FOP). J. Bone Miner. Res. 2007, 23, 305–313. [Google Scholar] [CrossRef]
  8. Kaplan, F.S.; Fiori, J.; De La Peña, L.S.; Ahn, J.; Billings, P.C.; Shore, E.M. Dysregulation of the BMP-4 Signaling Pathway in Fibrodysplasia Ossificans Progressiva. Ann. N. Y. Acad. Sci. 2006, 1068, 54–65. [Google Scholar] [CrossRef]
  9. Fiori, J.L.; Billings, P.C.; De La Peña, L.S.; Kaplan, F.S.; Shore, E.M. Dysregulation of the BMP-p38 MAPK Signaling Pathway in Cells from Patients with Fibrodysplasia Ossificans Progressiva (FOP). J. Bone Miner. Res. 2006, 21, 902–909. [Google Scholar] [CrossRef]
  10. De La Peña, L.S.; Billings, P.C.; Fiori, J.L.; Ahn, J.; Kaplan, F.S.; Shore, E.M. Fibrodysplasia Ossificans Progressiva (FOP), a Disorder of Ectopic Osteogenesis, Misregulates Cell Surface Expression and Trafficking of BMPRIA. J. Bone Miner. Res. 2005, 20, 1168–1176. [Google Scholar] [CrossRef]
  11. Hatsell, S.J.; Idone, V.; Wolken, D.M.A.; Huang, L.; Kim, H.J.; Wang, L.; Wen, X.; Nannuru, K.C.; Jimenez, J.; Xie, L.; et al. ACVR1R206H receptor mutation causes fibrodysplasia ossificans progressiva by imparting responsiveness to activin A. Sci. Transl. Med. 2015, 7, 303ra137. [Google Scholar] [CrossRef]
  12. Hino, K.; Ikeya, M.; Horigome, K.; Matsumoto, Y.; Ebise, H.; Nishio, M.; Sekiguchi, K.; Shibata, M.; Nagata, S.; Matsuda, S.; et al. Neofunction of ACVR1 in fibrodysplasia ossificans progressiva. Proc. Natl. Acad. Sci. USA 2015, 112, 15438–15443. [Google Scholar] [CrossRef] [Green Version]
  13. Ries, A.; Schelch, K.; Falch, D.; Pany, L.; Hoda, M.A.; Grusch, M. Activin A: An emerging target for improving cancer treatment? Expert Opin. Ther. Targets 2020, 24, 985–996. [Google Scholar] [CrossRef]
  14. Namwanje, M.; Brown, C.W. Activins and Inhibins: Roles in Development, Physiology, and Disease. Cold Spring Harb. Perspect. Biol. 2016, 8, a021881. [Google Scholar] [CrossRef]
  15. Pearsall, R.S.; Canalis, E.; Cornwall-Brady, M.; Underwood, K.W.; Haigis, B.; Ucran, J.; Kumar, R.; Pobre, E.; Grinberg, A.; Werner, E.D.; et al. A soluble activin Type IIA receptor induces bone formation and improves skeletal integrity. Proc. Natl. Acad. Sci. USA 2008, 105, 7082–7087. [Google Scholar] [CrossRef] [Green Version]
  16. Olsen, O.E.; Wader, K.F.; Hella, H.; Mylin, A.K.; Turesson, I.; Nesthus, I.; Waage, A.; Sundan, A.; Holien, T. Activin A inhibits BMP-signaling by binding ACVR2A and ACVR2B. Cell Commun. Signal. 2015, 13, 27. [Google Scholar] [CrossRef] [Green Version]
  17. Vanhoutte, F.; Liang, S.; Ruddy, M.; Zhao, A.; Drewery, T.; Wang, Y.; Ms, R.D.; Forleo-Neto, E.; Rajadhyaksha, M.; Herman, G.; et al. Pharmacokinetics and Pharmacodynamics of Garetosmab (Anti-Activin A): Results from a First-in-Human Phase 1 Study. J. Clin. Pharmacol. 2020, 60, 1424–1431. [Google Scholar] [CrossRef]
  18. Eekhoff, E.M.W. Garetosmab, an Inhibitor of Activin-A, Reduces Formation of Heterotopic Bone and Soft Tissue Flare-Ups in Patients with Fibrodysplasia Ossificans Progressiva. J. Bone Miner. Res. 2020, 35 (Suppl. S1), 17. [Google Scholar]
  19. Micha, D.; Voermans, E.; Eekhoff, M.E.W.; van Essen, H.W.; Zandieh-Doulabi, B.; Netelenbos, C.; Rustemeyer, T.; Sistermans, E.A.; Pals, G.; Bravenboer, N. Inhibition of TGFbeta signaling decreases osteogenic differentiation of fibrodysplasia ossificans progressiva fibroblasts in a novel in vitro model of the disease. Bone 2016, 84, 169–180. [Google Scholar] [CrossRef]
  20. Wang, H.; Shore, E.M.; Pignolo, R.J.; Kaplan, F.S. Activin A amplifies dysregulated BMP signaling and induces chondro-osseous differentiation of primary connective tissue progenitor cells in patients with fibrodysplasia ossificans progressiva (FOP). Bone 2018, 109, 218–224. [Google Scholar] [CrossRef]
  21. De Vries, T.J.; Schoenmaker, T.; Micha, D.; Hogervorst, J.; Bouskla, S.; Forouzanfar, T.; Pals, G.; Netelenbos, C.; Eekhoff, E.M.W.; Bravenboer, N. Periodontal ligament fibroblasts as a cell model to study osteogenesis and osteoclastogenesis in fibrodysplasia ossificans progressiva. Bone 2018, 109, 168–177. [Google Scholar] [CrossRef] [PubMed]
  22. Schoenmaker, T.; Botman, E.; Sariyildiz, M.; Micha, D.; Netelenbos, C.; Bravenboer, N.; Kelder, A.; Eekhoff, E.M.W.; De Vries, T.J. Activin-A Induces Fewer, but Larger Osteoclasts from Monocytes in Both Healthy Controls and Fibrodysplasia Ossificans Progressiva Patients. Front. Endocrinol. 2020, 11, 501. [Google Scholar] [CrossRef] [PubMed]
  23. Schoenmaker, T.; Deng, D.; De Vries, T.J. Tailored Teaching for Specialized (Para-)medical Students—Experience from Incorporating a Relevant Genetic Disease Throughout a Course of Molecular Cell Biology. Front. Public Health 2020, 8, 224. [Google Scholar] [CrossRef] [PubMed]
  24. Schoenmaker, T.; Wouters, F.; Micha, D.; Forouzanfar, T.; Netelenbos, C.; Eekhoff, E.M.W.; Bravenboer, N.; de Vries, T.J. The effect of Activin-A on periodontal ligament fibroblasts-mediated osteoclast formation in healthy donors and in patients with fibrodysplasia ossificans progressiva. J. Cell. Physiol. 2019, 234, 10238–10247. [Google Scholar] [CrossRef] [Green Version]
  25. Ruppeka-Rupeika, E.; Hogervorst, J.; Wouters, F.; Schoenmaker, T.; Forouzanfar, T.; de Vries, T.J. Osteogenic and osteoclastogenic potential of jaw bone-derived cells—A case study. J. Cell. Biochem. 2018, 119, 5391–5401. [Google Scholar] [CrossRef]
  26. De Vries, T.J.; Schoenmaker, T.; Wattanaroonwong, N.; Hoonaard, M.V.D.; Nieuwenhuijse, A.; Beertsen, W.; Everts, V. Gingival fibroblasts are better at inhibiting osteoclast formation than periodontal ligament fibroblasts. J. Cell. Biochem. 2006, 98, 370–382. [Google Scholar] [CrossRef]
  27. Franz, M.; Rodriguez, H.; Lopes, C.; Zuberi, K.; Montojo, J.; Bader, G.D.; Morris, Q. GeneMANIA update 2018. Nucleic Acids Res. 2018, 46, W60–W64. [Google Scholar] [CrossRef] [Green Version]
  28. Mostafavi, S.; Ray, D.; Warde-Farley, D.; Grouios, C.; Morris, Q. GeneMANIA: A real-time multiple association network integration algorithm for predicting gene function. Genome Biol. 2008, 9 (Suppl. S1), S4. [Google Scholar] [CrossRef] [Green Version]
  29. Warde-Farley, D.; Donaldson, S.L.; Comes, O.; Zuberi, K.; Badrawi, R.; Chao, P.; Franz, M.; Grouios, C.; Kazi, F.; Lopes, C.T.; et al. The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010, 38, W214–W220. [Google Scholar] [CrossRef]
  30. Merico, D.; Isserlin, R.; Stueker, O.; Emili, A.; Bader, G.D. Enrichment map: A network-based method for gene-set en-richment visualization and interpretation. PLoS ONE 2010, 5, e13984. [Google Scholar] [CrossRef]
  31. Chen, J.; Xu, H.; Aronow, B.J.; Jegga, A.G. Improved human disease candidate gene prioritization using mouse phenotype. BMC Bioinform. 2007, 8, 392. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, J.; Bardes, E.E.; Aronow, B.J.; Jegga, A.G. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009, 37, W305–W311. [Google Scholar] [CrossRef] [PubMed]
  33. Motta, M.; Chillemi, G.; Fodale, V.; Cecchetti, S.; Coppola, S.; Stipo, S.; Cordeddu, V.; Macioce, P.; Gelb, B.D.; Tartaglia, M. SHOC2 subcellular shuttling requires the KEKE motif-rich region and N-terminal leucine-rich repeat domain and impacts on ERK signalling. Hum. Mol. Genet. 2016, 25, 3824–3835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Matsunaga-Udagawa, R.; Fujita, Y.; Yoshiki, S.; Terai, K.; Kamioka, Y.; Kiyokawa, E.; Yugi, K.; Aoki, K.; Matsuda, M. The Scaffold Protein Shoc2/SUR-8 Accelerates the Interaction of Ras and Raf. J. Biol. Chem. 2010, 285, 7818–7826. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, A.M.; Lo, R.K.; Lee, M.M.; Wang, Y.; Yeung, W.W.; Ho, M.K.; Su, Y.; Ye, R.D.; Wong, Y.H. Galpha16 activates Ras by forming a complex with tetratricopeptide repeat 1 (TPR1) and Son of Sevenless (SOS). Cell. Signal. 2010, 22, 1448–1458. [Google Scholar] [CrossRef]
  36. Marty, C.; Browning, D.D.; Ye, R.D. Identification of Tetratricopeptide Repeat 1 as an Adaptor Protein That Interacts with Heterotrimeric G Proteins and the Small GTPase Ras. Mol. Cell. Biol. 2003, 23, 3847–3858. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, W.; Li, F.; Wang, K.; Cheng, B.; Guo, X. PAPSS2 Promotes Alkaline Phosphates Activity and Mineralization of Osteoblastic MC3T3-E1 Cells by Crosstalk and Smads Signal Pathways. PLoS ONE 2012. 7, e43475. [CrossRef] [Green Version]
  38. Kv, V.; Venkatachalam, K. Human 3′-phosphoadenosine 5′-phosphosulfate (PAPS) Synthase: Biochemistry, Molecular Biology and Genetic Deficiency. IUBMB Life 2003, 55, 1–11. [Google Scholar] [CrossRef]
  39. Le, P.T.; Bishop, K.A.; Maridas, D.E.; Motyl, K.J.; Brooks, D.J.; Nagano, K.; Baron, R.; Bouxsein, M.L.; Rosen, C.J. Sponta-neous mutation of Dock7 results in lower trabecular bone mass and impaired periosteal expansion in aged female Misty mice. Bone 2017, 105, 103–114. [Google Scholar] [CrossRef]
  40. Motyl, K.J.; Bishop, K.A.; DeMambro, V.E.; Bornstein, S.A.; Le, P.; Kawai, M.; Lotinun, S.; Horowitz, M.C.; Baron, R.; Bouxsein, M.L.; et al. Altered thermogenesis and impaired bone remodeling in Misty mice. J. Bone Miner. Res. Off. J. Am. Soc. Bone Miner. Res. 2013, 28, 1885–1897. [Google Scholar] [CrossRef] [Green Version]
  41. Cai, L.; Xiong, X.; Kong, X.; Xie, J. The Role of the Lysyl Oxidases in Tissue Repair and Remodeling: A Concise Review. Tissue Eng. Regen. Med. 2017, 14, 15–30. [Google Scholar] [CrossRef] [PubMed]
  42. Pischon, N.; Mäki, J.M.; Weisshaupt, P.; Heng, N.; Palamakumbura, A.H.; N’Guessan, P.; Ding, A.; Radlanski, R.; Renz, H.; Bronckers, T.A.L.J.J. Lysyl oxidase (lox) gene de-ficiency affects osteoblastic phenotype. Calcif. Tissue Int. 2009, 85, 119–126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, H.; Larson, M.; Jablonka-Shariff, A.; Pearl, C.A.; Miller, W.L.; Conn, P.M.; Boime, I.; Kumar, T.R. Redirecting in-tracellular trafficking and the secretion pattern of FSH dramatically enhances ovarian function in mice. Proc. Natl. Acad. Sci. USA 2014, 111, 5735–5740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Zhao, S.; Torii, S.; Yokota-Hashimoto, H.; Takeuchi, T.; Izumi, T. Involvement of Rab27b in the regulated secretion of pituitary hormones. Endocrinology 2002. 143, 1817–1824. [CrossRef]
  45. Bönnemann, C.G. The collagen VI-related myopathies: Muscle meets its matrix. Nat. Rev. Neurol. 2011, 7, 379–390. [Google Scholar] [CrossRef] [Green Version]
  46. Heumüller, S.E.; Talantikite, M.; Napoli, M.; Armengaud, J.; Mörgelin, M.; Hartmann, U.; Sengle, G.; Paulsson, M.; Moali, C.; Wagener, R. C-terminal proteolysis of the collagen VI α3 chain by BMP-1 and proprotein convertase(s) releases endotrophin in fragments of different sizes. J. Biol. Chem. 2019, 294, 13769–13780. [Google Scholar] [CrossRef]
  47. Smith, F.; Eady, R.; Leigh, I.; McMillan, J.; Rugg, E.; Kelsell, D.; Bryant, S.; Spurr, N.; Geddes, J.; Kirtschig, G.; et al. Plectin deficiency results in muscular dystrophy with epidermolysis bullosa. Nat. Genet. 1996, 13, 450–457. [Google Scholar] [CrossRef]
  48. Ishii, M.; Iwai, K.; Koike, M.; Ohshima, S.; Kudo-Tanaka, E.; Ishii, T.; Mima, T.; Katada, Y.; Miyatake, K.; Uchiyama, Y.; et al. RANKL-Induced Expression of Tetraspanin CD9 in Lipid Raft Membrane Microdomain Is Essential for Cell Fusion During Osteoclastogenesis. J. Bone Miner. Res. 2006, 21, 965–976. [Google Scholar] [CrossRef]
  49. Tachibana, I.; HemLer, M.E. Role of Transmembrane 4 Superfamily (Tm4sf) Proteins Cd9 and Cd81 in Muscle Cell Fusion and Myotube Maintenance. J. Cell Biol. 1999, 146, 893–904. [Google Scholar] [CrossRef]
  50. Wurthner, J.U.; Frank, D.B.; Felici, A.; Green, H.M.; Cao, Z.; Schneider, M.D.; McNally, J.G.; Lechleider, R.J.; Roberts, A.B. Transforming growth factor-beta receptor-associated protein 1 is a Smad4 chaperone. J. Biol. Chem. 2001, 276, 19495–19502. [Google Scholar] [CrossRef] [Green Version]
  51. Kaplan, J.; Kaplan, F.S.; Shore, E.M. Restoration of normal BMP signaling levels and osteogenic differentiation in FOP mesenchymal progenitor cells by mutant allele-specific targeting. Gene Ther. 2011, 19, 786–790. [Google Scholar] [CrossRef] [Green Version]
  52. Hino, K.; Horigome, K.; Nishio, M.; Komura, S.; Nagata, S.; Zhao, C.; Jin, Y.; Kawakami, K.; Yamada, Y.; Ohta, A.; et al. Activin-A enhances mTOR signaling to promote aberrant chondrogenesis in fibrodysplasia ossificans progressiva. J. Clin. Investig. 2017, 127, 3339–3352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hwang, C.; Pagani, C.A.; Das, N.; Marini, S.; Huber, A.K.; Xie, L.; Jimenez, J.; Brydges, S.; Lim, W.K.; Nannuru, K.C.; et al. Activin A does not drive post-traumatic heterotopic ossification. Bone 2020, 138, 115473. [Google Scholar] [CrossRef]
  54. Lees-Shepard, J.B.; Yamamoto, M.; Biswas, A.A.; Stoessel, S.J.; Nicholas, S.-A.E.; Cogswell, C.A.; Devarakonda, P.M.; Schneider, M.J., Jr.; Cummins, S.M.; Legendre, N.P.; et al. Activin-dependent signaling in fibro/adipogenic progenitors causes fibrodysplasia ossificans progressiva. Nat. Commun. 2018, 9, 471. [Google Scholar] [CrossRef] [Green Version]
  55. Dey, D.; Bagarova, J.; Hatsell, S.J.; Armstrong, K.A.; Huang, L.; Ermann, J.; Vonner, A.J.; Shen, Y.; Mohedas, A.H.; Lee, A.; et al. Two tissue-resident progenitor lineages drive distinct phenotypes of heterotopic ossification. Sci. Transl. Med. 2016, 8, 366ra163. [Google Scholar] [CrossRef]
  56. Culbert, A.L.; Chakkalakal, S.A.; Theosmy, E.G.; Brennan, T.A.; Kaplan, F.S.; Shore, E.M. Alk2 regulates early chon-drogenic fate in fibrodysplasia ossificans progressiva heterotopic endochondral ossification. Stem Cells 2014, 32, 1289–1300. [Google Scholar] [CrossRef] [Green Version]
  57. Olmsted-Davis, E.; Gannon, F.H.; Ozen, M.; Ittmann, M.M.; Gugala, Z.; Hipp, J.A.; Moran, K.M.; Fouletier-Dilling, C.M.; Schumara-Martin, S.; Lindsey, R.W. Hypoxic adipocytes pattern early heterotopic bone formation. Am. J. Pathol. 2007, 170, 620–632. [Google Scholar] [CrossRef] [Green Version]
  58. Schipani, E.; Ryan, H.E.; Didrickson, S.; Kobayashi, T.; Knight, M.; Johnson, R.S. Hypoxia in cartilage: HIF-1alpha is es-sential for chondrocyte growth arrest and survival. Genes Dev. 2001, 15, 2865–2876. [Google Scholar]
  59. Sviderskaya, E.V.; Novak, E.K.; Swank, R.T.; Bennett, D.C. The murine misty mutation: Phenotypic effects on mela-nocytes, platelets and brown fat. Genetics 1998, 148, 381–390. [Google Scholar] [CrossRef]
  60. Ramaswamy, G.; Sohn, P.; Eberhardt, A.; Serra, R. Altered responsiveness to TGF-β results in reduced Papss2 expression and alterations in the biomechanical properties of mouse articular cartilage. Arthritis Res. Ther. 2012. 14, R49. [CrossRef] [Green Version]
  61. Papaioannou, G.; Mirzamohammadi, F.; Kobayashi, T. Ras signaling regulates osteoprogenitor cell proliferation and bone formation. Cell Death Dis. 2016, 7, e2405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Schindeler, A.; Little, D.G. Ras-MAPK Signaling in Osteogenic Differentiation: Friend or Foe? J. Bone Miner. Res. 2006, 21, 1331–1338. [Google Scholar] [CrossRef] [PubMed]
  63. Cordeddu, V.; Di Schiavi, E.; A Pennacchio, L.; Ma’Ayan, A.; Sarkozy, A.; Fodale, V.; Cecchetti, S.; Cardinale, A.; Martin, J.; Schackwitz, W.; et al. Mutation of SHOC2 promotes aberrant protein N-myristoylation and causes Noonan-like syndrome with loose anagen hair. Nat. Genet. 2009, 41, 1022–1026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Tartaglia, M.; Zampino, G.; Gelb, B. Noonan Syndrome: Clinical Aspects and Molecular Pathogenesis. Mol. Syndr. 2010, 1, 2–26. [Google Scholar] [CrossRef] [Green Version]
  65. Saito, M.; Marumo, K. Collagen cross-links as a determinant of bone quality: A possible explanation for bone fragility in aging, osteoporosis, and diabetes mellitus. Osteoporos. Int. 2010, 21, 195–214. [Google Scholar] [CrossRef]
  66. Siegel, R.C.; Pinnell, S.R.; Martin, G.R. Cross-linking of collagen and elastin. Properties of lysyl oxidase. Biochemistry 1970, 9, 4486–4492. [Google Scholar] [CrossRef]
  67. Kumari, S.; Panda, T.K.; Pradhan, T. Lysyl Oxidase: Its Diversity in Health and Diseases. Indian J. Clin. Biochem. 2016, 32, 134–141. [Google Scholar] [CrossRef]
  68. Halberg, N.; Khan, T.; Trujillo, M.E.; Wernstedt-Asterholm, I.; Attie, A.D.; Sherwani, S.; Wang, Z.V.; Landskroner-Eiger, S.; Dineen, S.; Magalang, U.J.; et al. Hypoxia-inducible factor 1alpha induces fibrosis and insulin resistance in white adipose tissue. Mol. Cell. Biol. 2009, 29, 4467–4483. [Google Scholar] [CrossRef] [Green Version]
  69. Nohutcu, R.M.; McCauley, L.K.; Koh, A.J.; Somerman, M.J. Expression of extracellular matrix proteins in human periodontal ligament cells during mineralization in vitro. J. Periodontol. 1997, 68, 320–327. [Google Scholar] [CrossRef] [Green Version]
  70. Chien, H.H.; Lin, W.L.; Cho, M.I. Interleukin-1beta-induced release of matrix proteins into culture media causes in-hibition of mineralization of nodules formed by periodontal ligament cells in vitro. Calcif. Tissue Int. 1999, 64, 402–413. [Google Scholar] [CrossRef]
  71. Berendsen, A.D.; Smit, T.H.; Schoenmaker, T.; Walboomers, X.F.; Harris, S.E.; Everts, V.; Bronckers, A.L. Inorganic Phosphate Stimulates DMP1 Expression in Human Periodontal Ligament Fibroblasts Embedded in Three-Dimensional Collagen Gels. Cells Tissues Organs 2010, 192, 116–124. [Google Scholar] [CrossRef]
  72. Berendsen, A.D.; Smit, T.H.; Hoeben, K.A.; Walboomers, X.F.; Bronckers, A.L.; Everts, V. Alkaline phosphatase-induced mineral deposition to anchor collagen fibrils to a solid surface. Biomaterials 2007, 28, 3530–3536. [Google Scholar] [CrossRef]
Biomedicines 09 00629 g001 550
Figure 1. Donors group together. (A) PCA plot depicting inter-donor variability; (B) PCA plot depicting inter-experimental variability. Non-supervised principal component analysis shows that donors tend to group together, indicating that the variation between the control and FOP samples is larger than the comparison between control and experimental conditions. Inter-donor variability is shown in Figure 1A; each control and each FOP sample is represented by two same-colored dots, one indicating the control condition and one indicating the experimental condition (+Activin-A). Inter-experimental variability is shown in Figure 1B, each experimental condition is depicted by dots of the same color. Act:C = Control samples with Activin-A, ACT:F = FOP samples with Activin-A, Con: C = Control samples without Activin-A, Con:F = FOP samples without Activin-A. Samples Control 2 with Activin-A, Control 5 without Activin-A, and FOP 5 without Activin-A did not pass the QC filters and were not used in these and subsequent analysis.
Figure 1. Donors group together. (A) PCA plot depicting inter-donor variability; (B) PCA plot depicting inter-experimental variability. Non-supervised principal component analysis shows that donors tend to group together, indicating that the variation between the control and FOP samples is larger than the comparison between control and experimental conditions. Inter-donor variability is shown in Figure 1A; each control and each FOP sample is represented by two same-colored dots, one indicating the control condition and one indicating the experimental condition (+Activin-A). Inter-experimental variability is shown in Figure 1B, each experimental condition is depicted by dots of the same color. Act:C = Control samples with Activin-A, ACT:F = FOP samples with Activin-A, Con: C = Control samples without Activin-A, Con:F = FOP samples without Activin-A. Samples Control 2 with Activin-A, Control 5 without Activin-A, and FOP 5 without Activin-A did not pass the QC filters and were not used in these and subsequent analysis.
Biomedicines 09 00629 g001
Biomedicines 09 00629 g002 550
Figure 2. TGF-beta, BMP, and Activin signaling pathways are associated with upregulated genes by Activin-A. (A) Association with TGF-beta signaling pathway; (B) Association with BMP receptor signaling pathway. (C) Association with Activin signaling pathway. Pathway analysis on non-FDR-corrected differentially expressed genes with p < 0.01 performed with GeneMania (genemania.org, accessed on 21 September 2020). Figures are adjustments of the enrichment map of consolidated pathways TGF-beta signaling (MSigDB M2642) (A), BMP receptor signaling (MSigDB M181) (B), and signaling by Activin (MSigDB M26965) (C), Boxed genes are genes from the input list; the other genes are relevant genes added by GeneMania. Size of the nodes is dependent on the number of interactions that gene has with the network.
Figure 2. TGF-beta, BMP, and Activin signaling pathways are associated with upregulated genes by Activin-A. (A) Association with TGF-beta signaling pathway; (B) Association with BMP receptor signaling pathway. (C) Association with Activin signaling pathway. Pathway analysis on non-FDR-corrected differentially expressed genes with p < 0.01 performed with GeneMania (genemania.org, accessed on 21 September 2020). Figures are adjustments of the enrichment map of consolidated pathways TGF-beta signaling (MSigDB M2642) (A), BMP receptor signaling (MSigDB M181) (B), and signaling by Activin (MSigDB M26965) (C), Boxed genes are genes from the input list; the other genes are relevant genes added by GeneMania. Size of the nodes is dependent on the number of interactions that gene has with the network.
Biomedicines 09 00629 g002
Biomedicines 09 00629 g003 550
Figure 3. Activin-A induces differential gene expression only in FOP cells. Graphic display of differential gene expression (MA plot) in which the Log2 Fold Change induced by Activin-A treatment is plotted against the mean of normalized counts. Differentially expressed genes (p(adj) < 0.1) depicted as red dots. (A) MA plot of control samples shows no differentially expressed genes. (B) MA plot of FOP shows some differential gene expression induced by Activin-A, n = 4 for control and n = 5 for FOP samples. (C) Gene expression heatmap from median centered supervised gene expression analysis of the differentially expressed genes in the FOP cells. (p(adj) < 0.1, Z score is Log2 scale), n = 5.
Figure 3. Activin-A induces differential gene expression only in FOP cells. Graphic display of differential gene expression (MA plot) in which the Log2 Fold Change induced by Activin-A treatment is plotted against the mean of normalized counts. Differentially expressed genes (p(adj) < 0.1) depicted as red dots. (A) MA plot of control samples shows no differentially expressed genes. (B) MA plot of FOP shows some differential gene expression induced by Activin-A, n = 4 for control and n = 5 for FOP samples. (C) Gene expression heatmap from median centered supervised gene expression analysis of the differentially expressed genes in the FOP cells. (p(adj) < 0.1, Z score is Log2 scale), n = 5.
Biomedicines 09 00629 g003
Table
Table 1. GO terms associated with upregulated genes by Activin-A in FOP cells.
Table 1. GO terms associated with upregulated genes by Activin-A in FOP cells.
IDFunctionFDR B&HCoverage
GO:0032924Activin receptor signaling pathway8.03 × 10−1712/31
GO:0007178transmembrane receptor protein/threonine kinase signaling pathway4.06 × 10−1418/218
GO:0005072 *transforming growth factor beta receptor, cytoplasmic mediator activity2.25 × 10−117/10
GO:0090092regulation of transmembrane receptor protein/threonine kinase signaling pathway1.25 × 10−1013/132
GO:0030509BMP signaling pathway2.19 × 10−1011/79
GO:0032925regulation of activin receptor signaling pathway2.90 × 10−97/18
GO:0090100positive regulation of transmembrane receptor protein/threonine kinase signaling pathway2.53 × 10−78/53
GO:0060393regulation of pathway-restricted SMAD protein phosphorylation4.32 × 10−77/35
GO:0048185Activin binding7.92 × 10−75/10
GO:0060389pathway-restricted SMAD protein phosphorylation7.92 × 10−77/39
Top 10 GO terms enriched among the upregulated genes. Non-supervised functional enrichment analysis of GO terms associated with differentially upregulated genes by Activin-A in the FOP cells with p < 0.01 using GeneMania with default Benjamini and Hochberg’s False Discovery Rate (FDR B&H) settings. * GO:0005072 is replaced by GO:0000981.
Table
Table 2. Molecular functions regulated by Activin-A in FOP cells.
Table 2. Molecular functions regulated by Activin-A in FOP cells.
IDNameFDR B&HGenes from InputGenes in Annotation
GO:0044877protein-containing complex binding3.60 × 10−3231356
GO:0008092cytoskeletal protein binding4.73 × 10−2171061
actin filament binding
GO:0051015cell adhesion molecule binding4.73 × 10−27215
actin binding
GO:0050839protein-containing complex binding4.73 × 10−211525
GO:0003779cytoskeletal protein binding4.73 × 10−210451
actin filament binding
Non-supervised gene enrichment analysis of molecular function associated differentially expressed genes in the FOP cells with p < 0.01 using ToppGene suite using a Benjamini and Hochberg’s False Discovery Rate (FDR B&H) of 5%.
Table
Table 3. Biological processes regulated by Activin-A in FOP cells.
Table 3. Biological processes regulated by Activin-A in FOP cells.
IDNameFDR B&HGenes from InputGenes in Annotation
GO:0034330cell junction organization2.48 × 10−210309
GO:0034329cell junction assembly2.48 × 10−29245
GO:0007044cell-substrate junction assembly2.48 × 10−26103
GO:0009611response to wounding2.48 × 10−215756
GO:0009719response to endogenous stimulus2.80 × 10−2251834
GO:0071495cellular response to endogenous stimulus3.54 × 10−2221541
GO:0045196establishment or maintenance of neuroblast polarity3.54 × 10−223
GO:0045200establishment of neuroblast polarity3.54 × 10−223
Non-supervised gene enrichment analysis of biological processes associated differentially expressed genes in the FOP cells with p < 0.01 using ToppGene suite using a Benjamini and Hochberg’s False Discovery Rate (FDR B&H) of 5%.
Table
Table 4. Cellular components regulated by Activin-A in FOP cells.
Table 4. Cellular components regulated by Activin-A in FOP cells.
IDNameFDR B&HGenes from InputGenes in Annotation
GO:0005925focal adhesion6.75 × 10−817411
GO:0030055cell–substrate junction6.75 × 10−817421
GO:0005912adherens junction4.90 × 10−718560
GO:0070161anchoring junction5.55 × 10−718575
GO:0030054cell junction5.46 × 10−7241352
GO:0030424axon1.88 × 10−714817
GO:0030426growth cone2.16 × 10−27226
GO:0030427site of polarized growth2.16 × 10−27231
GO:0030496midbody2.82 × 10−26177
GO:0005911cell–cell junction3.08 × 10−210506
GO:0150034distal axon3.62 × 10−29432
Non-supervised gene enrichment analysis of differentially expressed genes in the FOP cells with p < 0.01 using ToppGene suite using a Benjamini and Hochberg’s False Discovery Rate (FDR B&H) of 5%.
Table
Table 5. Gene families regulated by Activin-A in FOP cells.
Table 5. Gene families regulated by Activin-A in FOP cells.
IDNameFDR B&HGenes from InputGenes in Annotation
939Plakins2.73 × 10−228
490Collagens2.73 × 10−2346
634Low density lipoprotein receptors2.73 × 10−2213
1149NADH:ubiquinone oxidoreductase core subunits2.73 × 10−2214
596Armadillo repeat containing|Importins3.63 × 10−2218
Non-supervised gene enrichment analysis of differentially expressed genes in the FOP cells with p < 0.01 using ToppGene suite using a Benjamini and Hochberg’s False Discovery Rate (FDR B&H) of 5%.
Table
Table 6. (A) Upregulated genes by Activin-A in FOP cells. (B) Downregulated genes by Activin-A in FOP cells.
Table 6. (A) Upregulated genes by Activin-A in FOP cells. (B) Downregulated genes by Activin-A in FOP cells.
A
ENSEMBL IDGene NameFold Changep(adj)-Value
ENSG00000041353RAB27B2.500.032
ENSG00000253626EIF5AL12.500.098
ENSG00000108061SHOC22.490.052
ENSG00000113312TTC12.400.099
ENSG00000198682PAPSS22.280.032
ENSG00000116641DOCK72.230.098
ENSG00000134108ARL8B1.880.098
ENSG00000089693MLF21.660.050
ENSG00000113083LOX1.470.060
ENSG00000198899MT-ATP61.310.032
B
ENSEMBL IDGene NameFold Changep(adj)-Value
ENSG00000163359COL6A3−1.250.074
ENSG00000178209PLEC−1.500.009
ENSG00000010278CD9−1.800.100
ENSG00000140264SERF2−1.850.089
ENSG00000160007ARHGAP35−1.900.032
ENSG00000117713ARID1A−2.420.098
ENSG00000117408IPO13−3.190.060
ENSG00000135966TGFBRAP1−3.900.021
Top 18 up and downregulated genes using supervised gene expression analysis. p(adj)-value are the p-values adjusted for an FDR of 10%.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schoenmaker, T.; Mokry, M.; Micha, D.; Netelenbos, C.; Bravenboer, N.; Gilijamse, M.; Eekhoff, E.M.W.; de Vries, T.J. Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients. Biomedicines 2021, 9, 629. https://doi.org/10.3390/biomedicines9060629

AMA Style

Schoenmaker T, Mokry M, Micha D, Netelenbos C, Bravenboer N, Gilijamse M, Eekhoff EMW, de Vries TJ. Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients. Biomedicines. 2021; 9(6):629. https://doi.org/10.3390/biomedicines9060629

Chicago/Turabian Style

Schoenmaker, Ton, Michal Mokry, Dimitra Micha, Coen Netelenbos, Nathalie Bravenboer, Marjolijn Gilijamse, E. Marelise W. Eekhoff, and Teun J. de Vries. 2021. "Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients" Biomedicines 9, no. 6: 629. https://doi.org/10.3390/biomedicines9060629

APA Style

Schoenmaker, T., Mokry, M., Micha, D., Netelenbos, C., Bravenboer, N., Gilijamse, M., Eekhoff, E. M. W., & de Vries, T. J. (2021). Activin-A Induces Early Differential Gene Expression Exclusively in Periodontal Ligament Fibroblasts from Fibrodysplasia Ossificans Progressiva Patients. Biomedicines, 9(6), 629. https://doi.org/10.3390/biomedicines9060629

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop