- Open Access
The ACVR1 R206H mutation found in fibrodysplasia ossificans progressiva increases human induced pluripotent stem cell-derived endothelial cell formation and collagen production through BMP-mediated SMAD1/5/8 signaling
Stem Cell Research & Therapyvolume 7, Article number: 115 (2016)
The Activin A and bone morphogenetic protein (BMP) pathways are critical regulators of the immune system and of bone formation. Inappropriate activation of these pathways, as in conditions of congenital heterotopic ossification, are thought to activate an osteogenic program in endothelial cells. However, if and how this occurs in human endothelial cells remains unclear.
We used a new directed differentiation protocol to create human induced pluripotent stem cell (hiPSC)-derived endothelial cells (iECs) from patients with fibrodysplasia ossificans progressiva (FOP), a congenital disease of heterotopic ossification caused by an activating R206H mutation in the Activin A type I receptor (ACVR1). This strategy allowed the direct assay of the cell-autonomous effects of ACVR1 R206H in the endogenous locus without the use of transgenic expression. These cells were challenged with BMP or Activin A ligand, and tested for their ability to activate osteogenesis, extracellular matrix production, and differential downstream signaling in the BMP/Activin A pathways.
We found that FOP iECs could form in conditions with low or absent BMP4. These conditions are not normally permissive in control cells. FOP iECs cultured in mineralization media showed increased alkaline phosphatase staining, suggesting formation of immature osteoblasts, but failed to show mature osteoblastic features. However, FOP iECs expressed more fibroblastic genes and Collagen 1/2 compared to control iECs, suggesting a mechanism for the tissue fibrosis seen in early heterotopic lesions. Finally, FOP iECs showed increased SMAD1/5/8 signaling upon BMP4 stimulation. Contrary to FOP hiPSCs, FOP iECs did not show a significant increase in SMAD1/5/8 phosphorylation upon Activin A stimulation, suggesting that the ACVR1 R206H mutation has a cell type-specific effect. In addition, we found that the expression of ACVR1 and type II receptors were different in hiPSCs and iECs, which could explain the cell type-specific SMAD signaling.
Our results suggest that the ACVR1 R206H mutation may not directly increase the formation of mature chondrogenic or osteogenic cells by FOP iECs. Our results also show that BMP can induce endothelial cell dysfunction, increase expression of fibrogenic matrix proteins, and cause differential downstream signaling of the ACVR1 R206H mutation. This iPSC model provides new insight into how human endothelial cells may contribute to the pathogenesis of heterotopic ossification.
Diseases of heterotopic ossification, where bone forms at an abnormal site, provide valuable opportunities to examine the mechanisms that regulate osteogenesis. The bone growth can range from small incidental nodules to catastrophic paralysis. Heterotopic ossification is found in a wide variety of diseases, including traumatic injury, severe burns, brain injury, and invasive surgeries [1–3]. However, the diversity of triggers poses a significant challenge for dissecting the cellular and molecular mechanisms that cause heterotopic ossification. Thus, genetic conditions of abnormal bone formation provide a valuable model for identifying the key pathways and regulators of heterotopic ossification in soft tissues.
Fibrodysplasia ossificans progressiva (FOP) is a congenital disease characterized by a large amount of heterotopic ossification in postnatal muscles and tendons [4, 5] and has been used as a prototypical model for studying heterotopic ossification. Most patients have a conserved mutation in the Activin A type 1 (ACVR1/ALK2) gene [6–8], a bone morphogenetic protein (BMP) receptor. BMPs are major regulators of bone formation. They were initially identified based on their ability to induce bone formation in soft tissues such as muscle and tendon. The majority of ACVR1 mutations in FOP are localized to a single amino acid change (R206H) that is thought to increase ACVR1 signaling activity.
Histologically, the initial fibrocellular infiltrate at a bone formation site is fibrotic and contains inflammatory cells and endothelial cells (ECs) . ECs overexpressing ACVR1 R206H may contribute to heterotopic ossification , possibly via abnormal BMP signaling, by undergoing endothelial-to-mesenchymal transition (EndMT) . However, it is unclear if human FOP ECs can transdifferentiate into osteogenic cells and directly contribute to bone formation, if another mechanism leads to the increase in ECs in FOP lesions [7, 9], if ECs are crucial in the very early steps of heterotopic ossification, or if human FOP ECs also respond abnormally to BMP or Activin A signals acting on the ACVR1 R206H receptor [11, 12].
Studies on human skeletal development are hampered by a number of factors, including the rare access to fetal material of developing bones, the difficulty of obtaining primary cells for detailed laboratory analysis, the ubiquitous nature of many of the ligands, and the technical and ethical problems surrounding the potential human genetic experiments needed to test pathways and hypotheses. Given that primary tissues cannot be obtained from patients with FOP because of heterotopic ossification at the surgical site, patient-derived human induced pluripotent stem cells (hiPSCs)  provide a strategy to create FOP cells for study in vitro .
The recent development of hiPSCs has made genetic diseases easier to study. hiPSCs can be derived directly from patients with existing genetic mutations and can form any cell in the body. hiPSCs also allow us to study the effect of single mutations in the endogenous locus and in a cell-autonomous fashion without the inherent complications of transgenic or overexpression studies. Thus, iPSC-derived tissues can model early developmental events and provide diseased human tissues that cannot be obtained from a patient.
Here, we use hiPSCs from patients with FOP  as a model of cell-autonomous signaling induced by the ACVR1 R206H mutation to test if the mutation increases the formation of potential osteoprogenitors in the EC lineage, and if hiPSC-derived endothelial cells (iECs) expressing ACVR1 R206H show increased plasticity towards osteogenesis. We also investigate the cellular mechanisms by which ECs may contribute to heterotopic ossification.
Cell culture and differentiation
Pluripotent hiPSC lines derived from control (wild type, WT) and FOP fibroblasts previously described [14, 15] were maintained in mTeSR1 medium (StemCell Technologies) on irradiated SNL (mouse fibroblast STO cell line transformed with neomycin resistance and murine LIF genes) feeder cells . SNLs were removed by at least one passage in feeder-free conditions on growth-factor-reduced Matrigel (Corning)-coated plates (150–300 μg/ml, 40 min coating) before use in differentiation assays. ROCK inhibitor Y-27632 (10 μM, StemCell Technologies) dissolved in 100 % DMSO was added to mTeSR1 at the time of passaging and removed the following day.
hiPSC lines were differentiated into iECs according to protocols previously described . Embryoid bodies (EBs) were formed from hiPSCs and cultured in aggregation medium Stem-Pro-34 (Invitrogen) supplemented with 2 mM glutamine, 150 mg/ml transferrin, 1 mM ascorbic acid, and 0.4 μM monothioglycerol (Sigma-Aldrich). EBs were cultured in aggregation medium with 5 ng/ml human basic fibroblast growth factor (bFGF), 3.6 ng/ml Activin A, and 12 ng/ml BMP4 (Peprotech). On day 6 of differentiation, cells were sorted for endothelial progenitor markers PECAM and KDR and plated on fibronectin (Sigma-Aldrich) in endothelial cell medium (ECM, ScienCell). These iECs were able to be passaged up to three times. iEC cultures were also treated with the Activin A inhibitor follistatin at 1 μg/ml (Peprotech).
hiPSCs maintained in feeder-free conditions were plated in 20 % mTeSR1 mixed with 80 % OB (osteoblastic mineralization) medium [Dulbecco’s Modified Eagle’s Medium with 20 % fetal bovine serum, 2 mM GlutaMAXTM, 10 mM glycerol-2-phosphate, 1 nM dexamethasone, 0.1 mM 2-mercaptoethanol, and 50 μg/ml L-ascorbic acid 2-phosphate sesquimagnesium salt hydrate, 1 % non essential amino acids] and Y-27632 (10μM) at 400,000 cells per well in gelatin-coated 24-well plates . The medium was replaced on day 2 with 100 % OB medium. The medium was changed every other day until day 12. Samples for the alkaline phosphatase staining were fixed on day 15 in 95 % ethanol for 10 min to 1 h at room temperature and incubated with substrate BCIP/NBT (Sigma-Aldrich) for 10 min at 37 °C. Staining was quantified via Image J .
Reverse Transcriptase-PCR and quantitative expression analysis
Total RNA was prepared using TRI Reagent (Sigma-Aldrich) and treated with the Turbo DNA-Free kit (Ambion). One microgram of total RNA was reverse transcribed into cDNA with the iScript cDNA Synthesis kit (BioRad) according to manufacturer instructions. Real-time quantitative PCR was performed in triplicate with VeriQuest Probe qPCR Master Mix (Affymetrix) or ABI’s Sybr Green PCR Master Mix on a ViiA 7 thermocycler (Life Technologies). Some low-yield cDNA samples were pre-amplified and analyzed using the BioMark 48.48 Dynamic Array nanofluidic chip (Fluidigm Inc., USA) according to manufacturers’ instructions. Taqman primer and Sybr Green probe sets are listed in Additional file 1: Table S1 and Additional file 2: Table S2. GAPDH or β-actin was used for normalization as an endogenous control.
iECs were fixed with 4 % paraformaldehyde/phosphate-buffered saline for 10 min at room temperature, then blocked with 5 % bovine serum albumin. Cells were stained overnight with primary antibodies to PECAM (5 μg/ml, R&D Systems) and VE-Cadherin (2 μg/ml, R&D Systems). Secondary antibodies were from Life Technologies: Alexa488-conjugated goat anti-mouse IgG (1:500) and Alexa555-conjugated goat anti-rabbit IgG (1:500). Nuclei were stained with DAPI in the ProLong® Gold Antifade (Life Technologies) mounting media. Images were taken using a light microscope (Nikon Eclipse E800 or Leica DMI 4000B).
hiPSCs cultured in mineralization medium were dissociated into single cells with collagenase type I (Worthington) for 1 h at 37 °C and then isolated by a Ficoll gradient (Histopaque 1191, Sigma-Aldrich). Accutase was used to generate single-cell suspensions from EBs plated overnight on collagen IV-coated plates or from iECs grown on fibronectin-coated plates. Cells were stained with PECAM1-AF488, KDR-APC, and VE-Cadherin-PerCP-Cy5.5 antibodies for endothelial markers, with CD90-AF488, CD73-PE, and CD105-PerCP-Cy5.5 antibodies (all from BD Pharmingen) for mesenchymal stem cell (MSC) markers. ICAM-1-PE antibody was used in our TNFα induction assay. Fluorescence intensity was determined for 10,000 cells in total and percentages shown in figures are the percentage of living cells that fall within the gate shown.
Vascular tube formation assay
Cells were seeded at 2.5 × 105 per well on 24-well plates pre-coated with growth factor-reduced Matrigel (Corning) and incubated for 24 h at 37 °C. Images were taken using a light microscope (Nikon and Leica).
Activin A levels were measured in iEC culture supernatant using an immunoassay solid-phase ELISA (R&D Systems). Samples were assayed in biological triplicates.
iECs were plated on fibronectin following sorting at a density of 7.5 × 104 cells per well of a 6-well plate and grown for 3 days in ECM medium (ScienCell). iECs were serum-starved for 1 h before a 40 min treatment with either 50 ng/ml of BMP4 or Activin A (R&D Systems). Cells were harvested in RIPA buffer (Pierce, Thermo Scientific) supplemented with 1X protease and a phosphatase inhibitor cocktail (Roche). Whole-cell lysates were prepared in Laemmli buffer (BioRad) and resolved in 4–20 % tris-glycine gels (BioRad). Primary antibodies towards SMAD1/5/8 (Santa Cruz Biotechnology), phospho-SMAD1/5/8 (Cell Signaling), SMAD2/3 (Cell Signaling), and phospho-SMAD2/3 (Cell Signaling) were used at a dilution of 1:1000. Anti-GAPDH antibody (Thermo Scientific) was used at a dilution of 1:10,000. Binding was visualized with horseradish peroxidase-conjugated antibodies (Cell Signaling) and ECL (Enhanced ChemiLuminescence) substrate (Thermo Scientific). An ImageQuant LAS 4000 (GE Healthcare) was used to image the blots and quantifications were done using Image J software.
TNFα activation assay
Cells were seeded at 1 × 105 per well on 6-well plates and treated with 10 ng/ml TNFα overnight at 37 °C in 5 % CO2. Cells were then harvested and stained for ICAM-1 and PECAM for fluorescence-activated cell sorting (FACS) analysis as discussed below.
For the transwell assay, 1 × 105 cells were added onto each transwell inserts (8.0 μm pore, Sigma-Aldrich). The inserts were incubated for 24 h at 37 °C in 5 % CO2 with serum-free media, with or without 50 ng/ml VEGF in the lower chamber. Membranes from each insert were fixed for 20 min and stained with crystal violet for 1 h.
Because each hiPSC line was derived clonally, and thus may display different behaviors, we treated each cell line as an individual biological replicate and pooled our results into control or FOP categories. All statistical analyses were completed using GraphPad Prism software. P values were calculated using the Student’s t test. p values ≤ 0.05 were considered statistically significant and are shown in the figures. Non-significant p values are not indicated for figure clarity.
WT and FOP hiPSCs can form endothelial cells with equal efficiency in a directed differentiation protocol
Multiple signaling molecules, including BMPs , regulate EC formation. Also, human umbilical vein endothelial cells (HUVECs) overexpressing ACVR1 R206H can acquire MSC-like phenotypes via EndMT . We previously found that FOP hiPSCs show increased mineralization compared to control hiPSCs when cultured in mineralizing conditions . To determine if the ACVR1 R206H mutation could increase the formation of ECs during mineralization, we examined the expression of two EC markers, PECAM and VE-Cadherin (markers of mature ECs), at two time points during culture (Additional file 3: Figure S1A). PECAM and VE-Cadherin expression showed only a trend towards increased expression in the FOP cultures. Since this could occur from a dilution effect from non-ECs, we used FACS to quantitate cells with PECAM and KDR, which marks endothelial progenitors , and cells co-expressing MSC markers CD90, CD73, and CD105. We found increased PECAM+/KDR+ cells (ECs, Additional file 3: Figure S1B) but not increased CD90+/CDE73+/CD105+ cells (MSCs, Additional file 3: Figure S1C), suggesting that FOP hiPSCs cultured in mineralizing conditions increased EC lineages but not MSC lineages .
We next created iECs to test if FOP heterotopic ossification could arise from increased formation of skeletal precursors in the EC lineage without overexpression of ACVR1 R206H (Fig. 1a) . These conditions, which contain BMP4, produced a consistent yield of approximately 20 % phenotypic iECs (PECAM+/KDR+) from both WT and FOP hiPSC lines (Fig. 1b). These findings suggest that the R206H mutation does not impair or favor iEC progenitor formation from hiPSCs in our protocol.
WT and FOP iECs show no functional differences
To characterize our FOP iECs, we next purified the iECs and cultured them in vitro to assess their endothelial properties (Fig. 1a, c, d). After sorting, hiPSC-derived PECAM+/KDR+ iECs yielded more than 95 % PECAM+/VE-Cadherin+ proliferative cells (Fig. 1c), could be passaged up to three times (data not shown), and retained their endothelial phenotype (Fig. 1d). WT and FOP iECs formed vessel-like structures when cultured on Matrigel (Additional file 4: Figure S2A), showed increased ICAM-1 upon TNFα activation (Additional file 4: Figure S2B), and had similar migration properties to VEGF in transwells (Additional file 5: Figure S3). Cells showed equivalent mRNA levels of the endothelial markers PECAM, vWF, KDR, and TIE2 (Fig. 2a). VE-Cadherin expression was significantly lower in FOP iECs, suggesting a partial loss of mature endothelial marker expression. Conversely, CD34 mRNA, a marker for hematopoietic stem cell progenitors or vascular endothelial progenitors and also expressed on MSCs , was significantly increased in FOP iECs (Fig. 2b). These results indicate no significant differences in functionality between WT and FOP iECs but that FOP iECs show increased expression of an endothelial progenitor-like marker and decreased expression of one mature endothelial gene.
FOP iECs show increased expression of alkaline phosphatase and Collagen 1/2
Since overexpression of R206H ACVR1 can activate EndMT to create MSC-like cells from HUVECs [9, 20], and vascular-associated cells such as pericytes can form osteoblasts , we tested if FOP iECs could be directed towards osteogenesis. Chondrogenic and osteogenic gene expression profiling of purified iECs cultured in ECM showed elevated COL1A1 and COL2A1 (Fig. 2b), two gees that code for extracellular matrix proteins found in bone and cartilage. FOP iECs also expressed higher levels of the mesenchymal marker FSP-1 gene (Fig. 2b), consistent with a predisposition towards EndMT. Some cultures of iECs showed increased levels of alkaline phosphatase (ALPL) mRNA (Fig. 2b), which was unexpected as the iECs were maintained in endothelial culture conditions. Although ALPL can mark immature osteoblasts, expression of chondrogenic and osteogenic transcription factors ( SOX9, SP7, and RUNX2) were unchanged in FOP iECs cultured in ECM (Additional file 6: Figure S4).
We next cultured our iECs in mineralization conditions to test if FOP iECs could form mineral-depositing cells via EndMT. After 15 days of culture, ALPL staining was increased in FOP iECs (Fig. 3a, b), again suggesting presence of osteogenic precursors; however, neither WT nor FOP iECs showed significant mineral deposition by von Kossa staining (data not shown), suggesting that mineralizing osteoblastic cells were not present. Gene expression analysis showed no difference in osteogenic markers except for increased COL1A1 expression in FOP iECs (Fig. 3c). These results suggest that the ACVR1 R206H mutation may not directly increase osteogenic or chondrogenic potential in FOP iECs, but may increase extracellular matrix production and local tissue fibrosis.
FOP iECs show increased formation in conditions not permissive for WT cells
Because the ACVR1 R206H mutation is thought to up-regulate BMP pathway activity [14, 22, 23], we asked if FOP hiPSCs could form iECs in conditions not normally permissive for control cells. hiPSCs were cultured in varying concentrations of BMP4 during the initial EB formation (Fig. 1a). BMP4 removal significantly inhibited iEC formation from control hiPSCs, but the FOP hiPSCs still formed endothelial-like cells in those conditions (Fig. 4a). It is important to note that our differentiation protocol also contained Activin A. Thus, the formation of ECs at low BMP4 concentrations could result from increased sensitivity to BMP4 but also from aberrant activation of the SMAD1/5/8 pathway by Activin A in FOP cells [11, 12]. To test this possibility, we differentiated our WT and FOP hiPSCs in the absence of Activin A. We found no significant differences in the ability of WT or FOP hiPSCs to form iECs in the absence of or with the full dose (3.6 ng/ml) of Activin A (Fig. 4b). Together, these results suggest that FOP hiPSCs can form ECs even in conditions that are not normally permissive for WT cells and that the difference is likely due to BMP signaling rather than Activin A signaling.
FOP iECs show increased SMAD1/5/8 signaling
To test if the increased plasticity in FOP ECs was due to activation of the BMP pathway, we next examined if SMAD signaling was activated in FOP iECs. When we initially examined SMAD signaling in FOP iECs cultured in fully supplemented media, we found no significant differences in SMAD1/5/8 activation when WT or FOP iECs were stimulated with BMP4 (Additional file 7: Figure S5A). However, after serum starvation, FOP iECs showed significantly increased phosphorylation of SMAD1/5/8 compared to WT iECs. The addition of BMP4 significantly increased the phosphorylation of SMAD2 and SMAD1/5/8 in WT iECs but only SMAD1/5/8 in FOP iECs (Fig. 4c, d). These results are consistent with the fact that BMPs have been reported to induce expression of EndMT markers in ECs  and that several other studies have reported increased SMAD2 signaling after BMP exposure in some cancer or transformed cells [25, 26].
Expression of SMAD6, which inhibits the phosphorylation of SMAD1/5/8 , and SMAD7, an inhibitor of SMAD2 and SMAD1/5/8 phosphorylation , was not significantly different between WT and FOP hiPSCs. Interestingly, SMAD6 expression was significantly increased in FOP iECs compared to WT iECs. We also noticed that expression levels depended on cell type. Expression of SMAD7  was significantly lower in iECs compared to hiPSCs in WT cells but not in FOP cells (Fig. 5a). Conversely, expression of SMAD6  was significantly increased in iECs compared to hiPSCs. These findings suggest differential regulation of SMAD pathways dependent on cell type. In addition, SMAD2 and SMAD6 may have an unexpected autoregulatory effect on the SMAD1/5/8 pathway in FOP iECs.
Because recent studies [11, 12] showed that the ACVR1 R206H mutation may alter Activin A signaling in MSCs, we investigated if FOP iECs might produce or respond to Activin A. Although the expression of Activin A was significantly increased in iEC lines compared to hiPSC lines, we found no significant difference in Activin A mRNA or protein levels in culture media between WT and FOP iECs (Fig. 5b). FOP hiPSCs treated with Activin A showed increased SMAD1/5/8 phosphorylation compared to WT hiPSCs (Additional file 7: Figure S5B). In contrast, Activin A increased SMAD2 and SMAD1/5/8 phosphorylation in both WT and FOP iECs (Fig. 4c, d). FOP iECs treated with Activin A showed a non-significant trend towards increased SMAD1/5/8 phosphorylation as compared to WT iECs. We found no significant differences in SMAD2 phosphorylation between WT and FOP iECs treated with Activin A. This suggests that the BMP and Activin A ligands may differentially activate separate SMADs in FOP cells as compared to WT cells, thus leading to increased propensity for EndMT. To further elucidate if the up-regulation of mesenchymal genes found in FOP iECs could be caused by Activin A in the media, we treated our iEC cultures with the Activin A inhibitor follistatin. We found no significant differences between WT and FOP iECs treated with follistatin (Additional file 8: Figure S6). Thus, the up-regulation of mesenchymal gene expression in FOP iECs is not likely to be caused by Activin A.
TGFβ1 is known to play a major role in EC EndMT [28, 29]. Stimulating our iECs with TGFβ1 significantly increased SMAD2 phosphorylation in both WT and FOP iEC but did not change SMAD1/5/8 phosphorylation (Fig. 4d). This result suggests that FOP iECs respond in the same manner as WT iECs to TGFβ1 ligand.
ACVR1 and type II receptors are differentially expressed in hiPSCs and iECs
Our unexpected findings of mild activation of SMAD1/5/8 by Activin A and SMAD2 by BMP4 (Fig. 4c, d) raised the question of whether cell type-specific expression of the ACVR1 receptors might be occurring. Expression of WT ACVR1 was significantly increased in WT iECs and OB iECs compared to FOP iECs and FOP OB iECs (Fig. 5c). We also found a significant increase in WT ACVR1 in WT iECs and in OB iECs compared to WT hiPSCs. However, expression of R206H ACVR1 only showed a slight increase in iECs and the ratio of R206H to WT ACVR1 showed no significant differences between hiPSC, iECs, and OB iECs (Fig. 5c).
Finally, we investigated if the differences seen in SMAD activation among different cell types might be related to altered expression of type II receptors known to bind to ACVR1, including ACVR2A, ACVR2B, and BMPRII . Although WT and FOP cells expressed similar levels of ACVR2A, ACVR2B, and BMPRII (Fig. 5d), we again noted significant differences in the absolute expression levels between cells types. ACVR2A gene expression was significantly lower in WT iECs and OB iECs than in hiPSCs. ACVR2B was significantly lower only in OB iECs compared to WT and FOP hiPSCs and iECs. Interestingly, BMPRII expression levels were significantly lower in both WT and FOP hiPSCs in comparison to the other two cell types. Thus these differences in type II receptor expression between cell types, independent of the presence of the ACVR1 R206H allele, may account for the differences found in SMAD signaling in ECs compared to hiPSCs.
Aberrant BMP signaling can be found in many conditions of heterotopic ossification, including vascular calcification , atherosclerosis , and heterotopic ossification . However, further delineation of the mechanisms inducing the abnormal bone formation has been limited by the broad expression patterns of the BMP ligands, potential differences in cell function between rodent models and humans, and the diversity of triggers that can incite heterotopic ossification. Genetic models using hiPSCs allows us to create human cell lineages with endogenous levels of gene expression for detailed studies.
Here, we used iECs to determine how ECs may contribute to heterotopic bone formation as a result of ACVR1 R206H expression. We show that iECs derived from FOP hiPSCs, which express physiological levels of the R206H ACVR1 activating mutation, are useful tools for dissecting the cellular and molecular mechanisms that cause this debilitating disease. In addition, our study shows that the ACVR1 receptor may play a key role in regulating EC commitment, and that the ACVR1 R206H mutation may underlie the pathogenesis of abnormal bone formation by changing the propensity of ECs to undergo EndMT, possibly through variations in the cell’s responses to different ligands for ACVR1 (Fig. 6). We also show that ECs carrying the ACVR1 R206H receptor produce increased amounts of collagen proteins, likely contributing to the tissue fibrosis characteristic of heterotopic bone lesions.
Several studies suggest that multiple cell types, including ECs , mesenchymal progenitors , pericytes [20, 34], Tie2 cells [35, 36], and MSCs , may contribute to heterotopic ossification. Unlike a prior report , our detailed analysis revealed no evidence for impaired endothelial formation from FOP hiPSCs. In contrast, we found that FOP iECs could form in conditions not normally permissive for WT cells and that FOP iECs could express genes that may predispose them towards EndMT. Although our iEC model does not fully recapitulate the previously reported findings of increased mature osteogenesis by EC lineages in mice , our current in vitro conditions may not be optimized for creating mature mineral-depositing osteoblasts and thus limits our ability to understand how ACVR1 signaling affects the later steps of human osteogenesis. We also found an increase of ECs in FOP hiPSCs cultured in osteogenic media, suggesting that FOP ECs may not directly contribute to increased bone formation by transdifferentiation into osteoblasts but rather by acting as a supporting cell that helps establish a bone-forming niche for a bona fide osteogenic precursor. This possibility is supported by the significantly increased levels of collagen production found in FOP iECs. Furthermore, we were unable to create in vivo osteogenic implants made purely from the iECs (not shown), suggesting that a more complex cellular model may be needed. Addressing these important challenges is critical for future studies using in vitro directed differentiation approaches to study bone formation.
While our study cannot exclude roles for MSCs [22, 38] or pericytes  in FOP, PECAM+/KDR+ EC precursors could still be a potential cell source of osteoprogenitors. Our observed increase in some mesenchymal, chondrogenic, and osteogenic gene expression in FOP iECs suggests that ACVR1 R206H can induce changes associated with cell fate. Cell lineage identity is necessary for homeostasis of most adult tissue and even subtle shifts in cell identity can severely impact tissues as they undergo regeneration after stress or injury. This loss of lineage determination, combined with the increased ability to make EC precursors in conditions not normally permissive in WT cells, may be critical early steps for developing heterotopic ossification in FOP. Indeed, a recent study has shown that increased cell plasticity can increase fibrogenesis in skeletal muscle, thus impairing normal tissue regeneration . Our data that FOP iECs produce more extracellular matrix proteins provide supporting evidence that this may be happening in the early heterotopic ossification lesion, particularly given that fibrosis is a key early feature in FOP  and that changes in material properties can affect cell fate [41, 42]. We would speculate that the formation of fibrosis contributes to changes in cell fate [41, 42] and abnormal tissue regeneration in FOP via increased cell plasticity in skeletal muscle .
The recent identification of aberrant Activin A signaling by ACVR1 R206H [11, 38] using transgene expression models suggested that Activin A could be anti-angiogenic  and activated osteoblast differentiation and extracellular matrix mineralization . Our results do not exclude FOP ECs as responders to Activin A, but indicate that FOP ECs do not contribute significantly to Activin A production. However, we found that FOP iECs responded differently from WT iECs when stimulated with BMP4 but not with Activin A or TGFβ1. We also saw an unexpected increase in SMAD1/5/8 signaling when WT hiPSCs were stimulated with Activin A. This is likely due to the different mRNA levels of SMAD6 and SMAD7 seen in hiPSCs versus ECs, as well as increased expression of ACVR1 and BMPRII and decreased levels of ACVR2A and ACVR2B  in iECs versus hiPSCs. Surprisingly, we also found that BMP4 treatment of WT iECs increased SMAD2 phosphorylation. This could suggest that, despite BMP4’s low binding affinity to type II receptors , BMP4 may still be able to induce SMAD2 signaling in the absence of a high-affinity ligand such as Activin A. Indeed, a few studies have reported that BMPs could activate noncanonical SMAD2 [25, 26]. These findings suggest that individual cell types may respond differently to each ligand stimulus  and that specific clinical features of FOP may be caused by tissue-specific expression of either receptors or ligands.
Together, our results show that EC lineages may contribute to heterotopic bone formation in FOP by priming an injury site with increased EC formation and changes in extracellular matrix production and fibrosis. In addition, BMPs and Activin A have a more complex interaction than previously thought, and differences in signaling via the type II receptors needs to be studied further in abnormal bone formation, tissue fibrosis, or the abnormal injury response in the pathogenesis of heterotopic bone formation. The broad expression of BMP and Activin A, particularly after injury, suggests that these mechanisms may be more generally applicable to non-genetic forms of heterotopic ossification  and potentially identify a new role for endothelial lineage cells in both normal and pathogenic osteogenesis.
ACVR1, Activin A type 1 receptor; ALPL, alkaline phosphatase; bFGF, Basic fibroblast growth factor; BMP, Bone morphogenetic protein; EB, embryoid bodies; EC, endothelial cell; ECM, endothelial cell medium; EndMT, endothelial-to-mesenchymal transition; FOP, fibrodysplasia ossificans progressiva; hESC, human embryonic stem cell; hiPSC, human induced pluripotent stem cell; HUVEC, human umbilical vein endothelial cells; ICAM-1, Intercellular adhesion molecule; iEC, hiPSC-derived endothelial cell; KDR, Kinase insert domain receptor; MSC, mesenchymal stem cell; MTG, PBS, phosphate-buffered saline; PECAM, Platelet endothelial cell adhesion molecule; TGFβ1, Transforming growth factor β1; TNFα, Tumor necrosis factor α; VE-cadherin, Vascular endothelial cadherin; VEGF, Vascular endothelial cell growth factor; WT, wild type
Hendricks HT, Geurts AC, van Ginneken BC, et al. Brain injury severity and autonomic dysregulation accurately predict heterotopic ossification in patients with traumatic brain injury. Clin Rehabil. 2007;21:545–53.
Kocic M, Lazovic M, Mitkovic M, et al. Clinical significance of the heterotopic ossification after total hip arthroplasty. Orthopedics. 2010;33:16.
Leung C, Casey AT, Goffin J, et al. Clinical significance of heterotopic ossification in cervical disc replacement: a prospective multicenter clinical trial. Neurosurgery. 2005;57:759–63. discussion 759–63.
Connor JM, Evans DA. Fibrodysplasia ossificans progressiva. The clinical features and natural history of 34 patients. J Bone Joint Surg Br. 1982;64:76–83.
Kaplan FS, Le Merrer M, Glaser DL, et al. Fibrodysplasia ossificans progressiva. Best Pract Res Clin Rheumatol. 2008;22:191–205.
Shore EM, Xu M, Feldman GJ, et al. A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva. Nat Genet. 2006;38:525–7.
Pignolo RJ, Shore EM, Kaplan FS. Fibrodysplasia ossificans progressiva: clinical and genetic aspects. Orphanet J Rare Dis. 2011;6:80.
Kaplan FS, Xu M, Seemann P, et al. Classic and atypical fibrodysplasia ossificans progressiva (FOP) phenotypes are caused by mutations in the bone morphogenetic protein (BMP) type I receptor ACVR1. Hum Mutat. 2009;30:379–90.
Medici D, Shore EM, Lounev VY, et al. Conversion of vascular endothelial cells into multipotent stem-like cells. Nat Med. 2010;16:1400–6.
Le Bras A, Vijayaraj P, Oettgen P. Molecular mechanisms of endothelial differentiation. Vasc Med. 2010;15:321–31.
Hatsell SJ, Idone V, Wolken DM, et al. ACVR1R206H receptor mutation causes fibrodysplasia ossificans progressiva by imparting responsiveness to activin A. Sci Transl Med. 2015;7:303ra137.
Hino K, Ikeya M, Horigome K, et al. Neofunction of ACVR1 in fibrodysplasia ossificans progressiva. Proc Natl Acad Sci U S A. 2015;112:15438–43.
Takahashi K, Tanabe K, Ohnuki M, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.
Matsumoto Y, Hayashi Y, Schlieve CR, et al. Induced pluripotent stem cells from patients with human fibrodysplasia ossificans progressiva show increased mineralization and cartilage formation. Orphanet J Rare Dis. 2013;8:190.
Spencer CI, Baba S, Nakamura K, et al. Calcium transients closely reflect prolonged action potentials in iPSC models of inherited cardiac arrhythmia. Stem Cell Reports. 2014;3:269–81.
McMahon AP, Bradley A. The Wnt-1 (int-1) proto-oncogene is required for development of a large region of the mouse brain. Cell. 1990;62:1073–85.
White MP, Rufaihah AJ, Liu L, et al. Limited gene expression variation in human embryonic stem cell and induced pluripotent stem cell-derived endothelial cells. Stem Cells. 2013;31:92–103.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.
Sidney LE, Branch MJ, Dunphy SE, et al. Concise review: evidence for CD34 as a common marker for diverse progenitors. Stem Cells. 2014;32:1380–9.
Cai J, Orlova VV, Cai X, et al. Induced pluripotent stem cells to model human fibrodysplasia ossificans progressiva. Stem Cell Reports. 2015;5:963–70.
Collett GD, Canfield AE. Angiogenesis and pericytes in the initiation of ectopic calcification. Circ Res. 2005;96:930–8.
Culbert AL, Chakkalakal SA, Theosmy EG, et al. Alk2 regulates early chondrogenic fate in fibrodysplasia ossificans progressiva heterotopic endochondral ossification. Stem Cells. 2014;32:1289–300.
de la Pena LS, Billings PC, Fiori JL, et al. Fibrodysplasia ossificans progressiva (FOP), a disorder of ectopic osteogenesis, misregulates cell surface expression and trafficking of BMPRIA. J Bone Miner Res. 2005;20:1168–76.
Maddaluno L, Rudini N, Cuttano R, et al. EndMT contributes to the onset and progression of cerebral cavernous malformations. Nature. 2013;498:492–6.
Holtzhausen A, Golzio C, How T, et al. Novel bone morphogenetic protein signaling through Smad2 and Smad3 to regulate cancer progression and development. FASEB J. 2014;28:1248–67.
Zhang H, Tian S, Klausen C, et al. Differential activation of noncanonical SMAD2/SMAD3 signaling by bone morphogenetic proteins causes disproportionate induction of hyaluronan production in immortalized human granulosa cells. Mol Cell Endocrinol. 2016;428:17–27.
Derynck R, Zhang YE. Smad-dependent and Smad-independent pathways in TGF-beta family signalling. Nature. 2003;425:577–84.
Cooley BC, Nevado J, Mellad J, et al. TGF-beta signaling mediates endothelial-to-mesenchymal transition (EndMT) during vein graft remodeling. Sci Transl Med. 2014;6:227ra234.
Correia AC, Moonen JR, Brinker MG, et al. FGF2 inhibits endothelial-mesenchymal transition through microRNA-20a-mediated repression of canonical TGF-beta signaling. J Cell Sci. 2016;129:569–79.
Attisano L, Carcamo J, Ventura F, et al. Identification of human activin and TGF beta type I receptors that form heteromeric kinase complexes with type II receptors. Cell. 1993;75:671–80.
Thomas PS, Sridurongrit S, Ruiz-Lozano P, et al. Deficient signaling via Alk2 (Acvr1) leads to bicuspid aortic valve development. PLoS One. 2012;7, e35539.
Yao Y, Shao ES, Jumabay M, et al. High-density lipoproteins affect endothelial BMP-signaling by modulating expression of the activin-like kinase receptor 1 and 2. Arterioscler Thromb Vasc Biol. 2008;28:2266–74.
Yu PB, Deng DY, Lai CS, et al. BMP type I receptor inhibition reduces heterotopic [corrected] ossification. Nat Med. 2008;14:1363–9.
Crisan M, Deasy B, Gavina M, et al. Purification and long-term culture of multipotent progenitor cells affiliated with the walls of human blood vessels: myoendothelial cells and pericytes. Methods Cell Biol. 2008;86:295–309.
Lounev VY, Ramachandran R, Wosczyna MN, et al. Identification of progenitor cells that contribute to heterotopic skeletogenesis. J Bone Joint Surg Am. 2009;91:652–63.
Wosczyna MN, Biswas AA, Cogswell CA, et al. Multipotent progenitors resident in the skeletal muscle interstitium exhibit robust BMP-dependent osteogenic activity and mediate heterotopic ossification. J Bone Miner Res. 2012;27:1004–17.
Suda RK, Billings PC, Egan KP, et al. Circulating osteogenic precursor cells in heterotopic bone formation. Stem Cells. 2009;27:2209–19.
Matsumoto Y, Ikeya M, Hino K, et al. New protocol to optimize iPS cells for genome analysis of fibrodysplasia ossificans progressiva. Stem Cells. 2015;33:1730–42.
Pessina P, Kharraz Y, Jardi M, et al. Fibrogenic cell plasticity blunts tissue regeneration and aggravates muscular dystrophy. Stem Cell Reports. 2015;4:1046–60.
McKleroy W, Lee TH, Atabai K. Always cleave up your mess: targeting collagen degradation to treat tissue fibrosis. Am J Physiol Lung Cell Mol Physiol. 2013;304:L709–21.
Trappmann B, Gautrot JE, Connelly JT, et al. Extracellular-matrix tethering regulates stem-cell fate. Nat Mater. 2012;11:642–9.
Watt FM, Huck WT. Role of the extracellular matrix in regulating stem cell fate. Nat Rev Mol Cell Biol. 2013;14:467–73.
Kaneda H, Arao T, Matsumoto K, et al. Activin A inhibits vascular endothelial cell growth and suppresses tumour angiogenesis in gastric cancer. Br J Cancer. 2011;105:1210–7.
Eijken M, Swagemakers S, Koedam M, et al. The activin A-follistatin system: potent regulator of human extracellular matrix mineralization. FASEB J. 2007;21:2949–60.
Lowery JW, Intini G, Gamer L, et al. Loss of BMPR2 leads to high bone mass due to increased osteoblast activity. J Cell Sci. 2015;128:1308–15.
Sebald W, Nickel J, Zhang JL, et al. Molecular recognition in bone morphogenetic protein (BMP)/receptor interaction. Biol Chem. 2004;385:697–710.
Aykul S, Martinez-Hackert E. Transforming growth factor-beta family ligands can function as antagonists by competing for type ii receptor binding. J Biol Chem. 2016;291:10792–804.
Sullivan MP, Torres SJ, Mehta S, et al. Heterotopic ossification after central nervous system trauma: a current review. Bone Joint Res. 2013;2:51–7.
The authors would like to thank Eileen Shore, Anne Zovein, Koen Schepers, Mary Nakamura, and Debby Gawlitta for valuable comments and discussion; Susan Fisher for support from the UCSF Stem Cell Core; and Kate Jordan and Kristie Yu for technical assistance.
This work was supported in part by a National Institutes of Health (NIH) K08 Career Development Grant (AR056299 to ECH), Doris Duke Charitable Foundation Clinical Scientist Development Award (2014099 to ECH), March of Dimes (5-FY12-167 and 1-FY14-211 to ECH), UCSF Pilot Grant for Junior Investigators, and a FOP Developmental Grant from the Center for Research in FOP and Related Disorders (to ECH). AU and ECH received support from an NIH Diversity Supplement (R03AR060986-02). EB received support from the California Institute for Regenerative Medicine Fellowship Program to UCSF (TG2-01153). BMM received support through the CIRM Bridges Program to the City College of San Francisco (TB1-01188).
Availability of data and materials
The data supporting the conclusions of this article are included within the article and supporting files.
EB was involved in the design, collection, and assembly of data, data analysis and interpretation, and manuscript writing. BMM, WL, HK, AU assisted with the experiments and the differentiation of hiPSCs. SAW helped design and perform transwell experiments. MPW, CVT, and DS assisted with the development and analysis of endothelial differentiation from hiPSCs. ECH was responsible for the conception and design, financial support, and manuscript writing. All authors were involved in drafting the manuscript and revising it. All authors read and approved the final manuscript.
The authors declare that they have no competing interest. ECH receives research grant support from Clementia Pharmaceuticals for a clinical trial related to FOP but not related to this work.
Consent for publication
This manuscript does not contain any individual person’s data in any form.
Ethics approval and consent to participate
All human stem cell studies, procedures, and written consents were approved by the UCSF Committee on Human Research and the UCSF Gamete and Embryonic Stem Cells Research Committee.
Taqman Primers-gene expression. (DOCX 37 kb)
Sybr Green primers-gene expression. (DOC 27 kb)
Increased endothelial marker expression during osteoblast differentiation of hiPSCs and differentiation of hiPSCs into ECs. (A) Quantitation of PECAM and VE-Cadherin gene expression by qPCR during 12 days of culture in osteogenic medium. (B) FACS analysis of mean percentage of cells expressing both PECAM and KDR on day 6 and 12 of osteogenic differentiation of three WT and three FOP hiPSCs lines with at least three replicates for each cell line. Error bars represent mean ± one SD. * p < 0.05, ** p < 0.01 by Student’s t test. (C) Mean percentage of cells co-expressing MSC markers (CD90+/CD73+/CD105+) on day 12 of osteogenic differentiation of three WT and three FOP hiPSCs lines. Error bars represent mean ± one standard deviation of at least three independent replicates for each of the three WT and four FOP cell lines. MSCs were not significantly increased in FOP hiPSCs during osteogenic differentiation by Student’s t test. (TIF 116 kb)
Functional characterization of WT and FOP iECs. (A) Both WT and FOP iECs have the ability to form tube-like structures in an in vitro Matrigel assay. Scale bars, 300 μm. (B) WT and FOP iECs showed activation upon TNFα stimulation. WT and FOP iECs were stimulated with 10 ng/ml TNFα overnight. Cells co-expressing ICAM-1 and PECAM were quantified by FACS analysis. ** p < 0.01, *** p < 0.005, **** p < 0.0001 by Student’s t test. Error bars represent mean ± one standard deviation of at least three independent replicates for each HUVEC WT and FOP iEC line. (TIF 301 kb)
WT and FOP iECs show similar qualitative migration properties to VEGF in a transwell assay. Scale bars, 100 μm. (TIF 1177 kb)
Characterization of osteogenic and chondrogenic potential of FOP iECs. Gene expression analysis of SMA, SOX9, MMP10, SP7, and RUNX2 by qPCR of WT and FOP iECs. None of these gene expression levels were statistically different by Student’s t test. Error bars represent mean ± one standard deviation of at least three independent replicates for each of the three WT and four FOP iEC lines. (TIF 158 kb)
SMAD signaling activity of hiPSCs and iECs. (A) Quantitation of SMAD1/5/8 phosphorylation without serum starvation upon BMP4 stimulation in WT and FOP iECs. p-SMAD1/5/8 protein expression was normalized to GAPDH protein expression. At least three separate experiments were run for WT and FOP iECs for each group. Error bars represent mean ± one standard deviation. Mean values were not statistically significant by Student’s t test. (B) Representative western blot showing activation of SMAD2 and SMAD1/5/8 pathways upon BMP4 and Activin A stimulation in WT and FOP hiPSCs. Samples were derived from the same experiment and processed in parallel. (TIF 221 kb)
Inhibition of Activin A by follistatin does not down-regulate mesenchymal genes in FOP iECs. iECs were treated with follistatin. Error bars represent mean ± one standard deviation of at least two independent replicates for two WT and three FOP iEC lines; * p < 0.05, ** p < 0.01 by Student’s t test. (TIF 99 kb)