Next Article in Journal
Endocrine Disruptors and Endometrial Cancer: Molecular Mechanisms of Action and Clinical Implications, a Systematic Review
Previous Article in Journal
The Impact of Sub-Millisecond Damage Fixation Kinetics on the In Vitro Sparing Effect at Ultra-High Dose Rate in UNIVERSE
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis

1
Musculoskeletal Research Group, Institute of Cellular Medicine, Newcastle University, Framlington Place, Newcastle upon Tyne NE2 4HH, UK
2
Blood Sciences, South Tees Hospitals NHS Foundation Trust, Middlesbrough TS4 3BW, UK
3
Center for Psychopharmacology, Diakonhjemmet Hospital, 0319 Oslo, Norway
4
Department of Engineering, Faculty of Science, Durham University, Durham DH1 3 LE, UK
5
Department of Medical Biochemistry, Oslo University Hospital, 0424 Oslo, Norway
6
Unger-Vetlesen Institute, Lovisenberg Diaconal Hospital, 0440 Oslo, Norway
7
School of Biomedical, Nutritional and Sport Sciences, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne NE2 4HH, UK
8
Blood Sciences, The Newcastle upon Tyne Hospitals NHS Foundation Trust, Royal Victoria Infirmary, Newcastle upon Tyne NE2 4HH, UK
9
Department of Plastic and Reconstructive Surgery, Oslo University Hospital, 0424 Oslo, Norway
*
Author to whom correspondence should be addressed.
Current address: Division of Trauma & Orthopaedic Surgery, University of Cambridge, Cambridge CB2 2QQ, UK.
Int. J. Mol. Sci. 2022, 23(6), 2957; https://doi.org/10.3390/ijms23062957
Submission received: 4 January 2022 / Revised: 2 March 2022 / Accepted: 5 March 2022 / Published: 9 March 2022
(This article belongs to the Section Molecular Immunology)

Abstract

:
Mechanical loading exerts a profound influence on bone density and architecture, but the exact mechanism is unknown. Our study shows that expression of the neurological transcriptional factor zinc finger of the cerebellum 1 (ZIC1) is markedly increased in trabecular bone biopsies in the lumbar spine compared with the iliac crest, skeletal sites of high and low mechanical stress, respectively. Human trabecular bone transcriptome analyses revealed a strong association between ZIC1 mRNA levels and gene transcripts characteristically associated with osteoblasts, osteocytes and osteoclasts. This supposition is supported by higher ZIC1 expression in iliac bone biopsies from postmenopausal women with osteoporosis compared with age-matched control subjects, as well as strongly significant inverse correlation between ZIC1 mRNA levels and BMI-adjusted bone mineral density (BMD) (Z-score). ZIC1 promoter methylation was decreased in mechanically loaded vertebral bone compared to unloaded normal iliac bone, and its mRNA levels correlated inversely with ZIC1 promoter methylation, thus linking mechanical stress to epigenetic control of gene expression. The findings were corroborated in cultures of rat osteoblast progenitors and osteoblast-like cells. This study demonstrates for the first time how skeletal epigenetic changes that are affected by mechanical forces give rise to marked alteration in bone cell transcriptional activity and translate to human bone pathophysiology.

1. Introduction

Bone is a metabolically active and dynamic organ that remodels throughout life by coordinated activities of bone-resorbing osteoclasts and bone-forming osteoblasts. Mechanical loading plays an important role in the maintenance of bone density, but knowledge of the cellular and molecular biology of mechanosensing remains limited. Mature bone is removed by remodeling that occurs as an adaptive response to different factors, including mechanical loading and replacement of damaged bone following micro- and macrofractures. Studies [1,2] have uncovered the initial molecular mechanisms of ZIC1 osteocyte actions and implicated a role in human bone pathophysiology. Previously, ZIC1 had been related to craniofacial developmental defects [3]. Its function was further assessed in a large-scale GWAS meta-analysis focusing specifically on BMD variation in the skull, finding 4 novel loci among 59 with genome-wide significance, including ZIC1 [4]. The gene was further screened in zebrafish crispants that displayed low skull BMD and abnormal suture patterning [4].
Both osteogenic cell types, osteoblasts and osteocytes, act as mechanosensors and transduce mechanical loading forces into inter- and intracellular biochemical signals [5]. Osteocytes, the most abundant cell type, have many unique characteristics that make them ideal mechanosensors [6,7,8]. In addition to their large number and special distribution in bone, the cells show special morphological features and cell–cell connectivity with each other and several other types of bone cells [7]. The molecular details of how external mechanical load is transmitted to individual cells and specifically to osteocytes have not yet been established. Furthermore, how mechanical loading contributes to the etiology of postmenopausal osteoporosis, a condition of increased bone resorption and fragile bones, is not yet fully elucidated.
The relative significance and contribution of the force, arising from fluid and tissue tethering, in transmitting macroscopic mechanical loading to the microscopic cellular level is uncertain. In addition, mechanical forces acting upon bone generate electrical potentials [9], but their possible physiological importance in influencing the activity of bone cells is unknown.
Skeletal cells, osteoblasts and especially osteocytes respond to mechanical loading by producing a variety of cascading biochemical signals, a process whereby mechanical stress is translated to a biological response [5,10,11]. Mechanically activated osteocytes and osteoblasts modulate several signaling systems; among these, the Wnt/β-catenin pathway clearly plays a powerful role in the regulation of bone mass and new bone formation [12,13]. Previously, it was reported that the zinc finger protein of the cerebellum (ZIC1) may mediate fluid-flow-shear-stress-induced modulation by activating the TCF/LEF transcription factor in cell culture [1]. ZIC proteins have been shown to be potential modulators of the hedgehog-mediated signaling pathway and can interact with all GLI proteins, including the repressive form. GLI proteins function downstream of the hedgehog signaling pathway and act as both transcriptional activators and repressors [14,15]. How mechanical loading influences the epigenetic and genetic mechanisms and its downstream effect on the skeletal cell transcriptome and protein translation is, however, not known.
This study investigated whether modulation of the ZIC1 transcription factor at skeletal sites of high stress in postmenopausal healthy and osteoporotic women was accompanied by DNA methylation changes measured in trabecular bone biopsies. ZIC1 promoter methylation was significantly lower in bone taken from the lumbar spine compared to the relatively less mechanically stressed os ilium. The extent of DNA methylation of ZIC1 was linearly negatively correlated with ZIC1 expression. ZIC1 mRNA levels in biopsies from the iliac of postmenopausal women showed significant inverse correlation with BMD and were significantly higher in women with osteoporosis than in healthy controls. ZIC1 promoter methylation was also correlated with ZIC1 mRNA levels in bone cell cultures, which corroborates its intracellular role in mechanosensing and bone turnover.

2. Results

2.1. Ex Vivo Studies

2.1.1. ZIC1 mRNA Expression and Promoter Methylation Differences in Men

Previous initial transcriptome analysis compared gene expression in male human bone biopsy samples taken from lumbar spine (LS) and iliac bone (IB) sites that experience high and low levels of mechanical stress, respectively [2]. These studies showed that the ZIC1 transcript was profoundly upregulated in the trabecular biopsies from the lamina of the L5 LS vertebrae (77-fold; p < 109) in comparison with IB. The present study shows that the degree of DNA methylation also differs markedly between mechanically stressed LS and unstressed IB. The question as to whether ZIC1 promoter methylation correlated with ZIC1 expression at the two skeletal sites was addressed using male bone biopsies from LS and IB (Supplementary Table S1). Even if the degree of methylation varied between the denoted positions, CpG1–CpG20, there was nevertheless a significantly higher degree of methylation in the trabecular biopsies from the IB than from LS (Table 1A). Male trabecular bone ZIC1 mRNA levels and the degree of promoter methylation at positions CpG1 to CpG5, showed a strong negative correlation (Table 1A, Figure 1A). The most intensely methylated promoter region (CpG1) had a 55% higher degree of methylation in IB compared to LS vertebrae (49.2 vs. 31.8; p < 0.001) (Figure 1C). All but one (CpG 19) of the 20 CpG methylation sites in the ZIC1 promoter region showed strong significant inverse correlation between the degree of methylation and ZIC1 expression (Pearson r = −0.33 to −0.74, p = 9 × 10−4 to 3 × 10−3) (Table 1A). Furthermore, ZIC1 mRNA levels and degree of DNA methylation in each of 20 promoter CpG sites correlate strongly (r = 0.72, p = 0.0009 for the CpG site in position 1 (CpG1)) (Table 1A).
Interestingly, of the three Affymetrix HU133plus2.0 microarray probe sets for ZIC1 array expression profile, (236896_at, 234716_at and 206373_at), the probe set 206373_at showed the most variability between IB and LS (Table 1B). The same probe set showed differential expression between osteoporotic and control samples and revealed an inverse association between ZIC1 expression and total femoral Z-score adjusted for BMI (Figure 1B).

2.1.2. Relationship between ZIC1 Bone DNA Methylation, mRNA Levels and Postmenopausal Osteoporosis

ZIC1 expression determined by Affymetrix microarray and PCR analysis in the iliac bone of postmenopausal women with a range of BMD showed a markedly significant inverse correlation with the total BMD adjusted for BMI. A negative association between ZIC1 expression was seen with total body BMD Z-score that had been adjusted for body mass index (BMI) (r = −0.398; p = 0.003) and of BMI-adjusted total femoral BMD (r = −0.383; p = 0.003) and femoral neck BMD (r = −0.331; p = 0.010) (Table 2, Figure 1B). The demographic and other relevant data for the cohort are given in Supplementary Table S2. There was also a significant inverse correlation with the adjusted lumbar spine L1–L4 BMD (r = −0.337; p = 0.010). However, no significant correlation was seen with lean body or fat mass (Table 2). Further refined analysis compared ZIC1 expression in postmenopausal osteoporotics (n = 23), defined as T-score ≤ −2.5 at the lumbar spine or hip region with at least one fragility fracture, and age balanced female controls (n = 30) with T-score ≥ −0.3), and showed significantly lower ZIC1 mRNA levels in these patients (Mean ± SD 3.88 ± 0.14 vs. 4.03 ± 0.26, p = 0.008).
In view of the increased mRNA expression of ZIC1 in stressed and high-turnover bone, a possible difference in ZIC1 DNA methylation, as well as an association between ZIC1 expression, BMD and osteoporosis was sought in a previously characterised cohort comprised of Norwegian postmenopausal osteoporotic women (n = 40) and healthy control subjects (n = 29) (Supplementary Table S3) [16,17]. Interestingly, significant differences were observed in CpG methylation levels between osteoporotic (OP) and healthy subjects, as well as after adjustment for multiple testing (Table 3).

2.1.3. Correlation of ZIC1 Expression with Key Bone-Related Transcripts

The correlation between the key bone-related transcripts and ZIC1 expression in loaded trabecular bone (LS) compared to unloaded (IB) bone from male subjects was investigated (Supplementary Table S1), selecting bone-cell-characteristic or widely employed markers of osteoblasts, osteoclasts and osteocytes. The selected transcripts reflecting osteoblast matrix-forming activities were COL1A1, IBSP, SPARC, BGLAP and CDH11 (Figure 2A and Supplementary Figure S1A) and, as makers for osteoblast differentiation, BMP2 and RUNX2 (Supplementary Figure S1). These transcripts showed strong positive correlation with ZIC1 in LS but weaker or no correlation in IB (Figure 2A and Supplementary Figure S1A). The PTHR1 transcript, used to indicate osteoblast number, also showed a strong positive correlation with ZIC1 expression (Supplementary Figure S1A). In contrast, transcripts related to osteoclast numbers, such as CALCR and OSCAR, had significant negative correlation with ZIC1 expression in the LS but weaker or no correlation, respectively, in the IB, notably due to fewer samples (Supplementary Figure S1B). On the other hand, the transcripts reflecting osteoclast activity, ACP5 and CTSK, showed strong positive correlation with ZIC1 expression (Figure 2B and Supplementary Figure S1B). These observations are consistent with increased ZIC1 expression being positively associated with both enhanced osteoblast number and increased cellular activity, whereas its relation to osteoclast activity and number is inverse and appears more complex. Although osteoclast activity was elevated, as reflected by positive correlation with ACP5 and CTSK transcripts, the actual number of osteoclasts, as reflected by CALCR and OSCAR transcripts, was suppressed (Supplementary Figure S1B). The osteoclastic marker CALCR also showed a significant inverse correlation to ZIC1 in LS (r = −0.54, p = 0.008) but, again, with no association in IB. Expression of PDPN, MEPE and SOST, osteocyte characteristic transcripts, showed strong correlation with ZIC1 mRNA at both high (LS) and low (IB) bone turnover sites (Figure 2C and Supplementary Figure S1C).

2.2. In Vitro Studies

2.2.1. ZIC1 Induction by Fluid Shear Stress in Human SaOS2 Cells and Human Primary Skin Fibroblasts (HSF)

Further in vitro studies were conducted to determine possible association between the degree of ZIC1 promoter methylation and ZIC1 mRNA expression. For this purpose, a rapidly proliferating osteoblast-like human osteosarcoma cell line (SAOS2) was subjected to 24 h oscillatory fluid flow, resulting in fluid shear stress (FSS). Basal ZIC1 mRNA expression levels in SAO2 cells, as determined by qRT-PCR, was found to be 11.4-fold higher than in HSFs (ΔCT 6.98 vs. 10.39). The differences in ZIC1 expression closely mirrored the degree of ZIC promoter methylation in the two cell types, which was considerably higher in SaOS2 cells than in fibroblasts (Table 4A,B), thereby indicating an important regulatory role of promoter methylation in ZIC1 transcription. When SAOS2 were subjected to 24 h of FSS, there was a significant further rise in ZIC1 expression, which could be demonstrated despite relatively high pre-existing basal expression (Table 4B). The FSS-induced induction of ZIC1 in SAOS2 cells was accompanied by an increased expression of ciliary protein KIF3A, which is a subunit of the kinesin II intraflagellar transport (IFT) protein belonging to the kinesin-2 family, and of SP7 transcription factor (Osterix) (Table 4B). In contrast to SAOS2 cells, in HSF cells, exposure to fluid shear stress failed to induce a significant increase in ZIC1 expression despite the low basal expression (Table 4C). However, neither cell type responded with statistically significant changes in the degree of ZIC1 promotor methylation.

2.2.2. Induction of ZIC1 mRNA by FSS in Rat Osteoprogenitor Cells

Further in vitro investigations were designed to elucidate the role of ZIC1 in the differentiation of osteoblasts from rat osteoprogenitor cells (ROPs) and to determine whether mechanical stresses influence this process. ROPs were cultured in osteogenic medium and subjected to FSS for three days, followed by qRT-PCR analyses of osteogenic markers and Zic1. The results showed that FSS led to induction of the key osteogenic markers Runx2, Alpl and Col1a1 when compared to controls (Figure 3A,B). However, the mRNA level of Taz1, an important marker of osteogenic potential, was unaltered (Figure 3A,B). When the rat osteoprogenitor cells were cultured in normal medium but subjected to similar magnitude of FSS, no induction of Zic1 or Col1a1 was seen over three days (Figure 3C). However, a further exposure of these cells to FSS in normal medium led to induction of both Col1a1 and Zic1 (Figure 3C).

2.2.3. Immunocytochemical Localization of Rat ZIC1 Protein

Immunocytochemical localization of ZIC1 protein in rat osteoprogenitor cells cultured in osteogenic medium showed greater overall concentration of ZIC1 protein in rat osteoprogenitor cells subjected to FSS than static controls (Figure 4 and Figure 5), mirroring the increased expression of Zic1 mRNA (Figure 3). In static control cells, RUNX2 showed considerable perinuclear staining, with a proportion of cells lacking RUNX2 in the nucleus (Figure 4 and Figure 5). The fluorescence intensity of both RUNX2 and ZIC1 in the cytoplasm and nuclei of cells subjected to FSS compared to static controls showed increased overall protein expression, as well as evident translocation from cytoplasm to nuclei (Figure 4A–F). A plot of the intensity of ZIC1 and RUNX2 fluorescence in nuclei of cells subjected to FSS and static cells shows clearly that the nuclei of cells subjected to FSS reveal greater nuclear translocation of and RUNX2 (Figure 5A,B). The graph shows a linear relationship, suggesting a direct relationship between translocation of ZIC1 and RUNX2 in the nucleus after application of FSS (Figure 5C).

3. Discussion

It is generally accepted that external loading influences bone cells in two main ways: as a shear force and through cellular deformation [18,19]. The bone matrix, a heterogeneous structure, may act as a stress conductor and amplify the strain experienced by the osteocytes by as much as 15-fold compared to that experienced by the macroscopic bone [18,19,20].
We therefore studied bone biopsies taken from high or low mechanical loading sites representing elevated and suppressed bone turnover, respectively, and found a close negative correlation between ZIC1 promoter methylation and gene expression. At the site of high loading and increased bone turnover, ZIC1 mRNA expression was elevated, and ZIC1 promoter methylation was reduced. In contrast, in the iliac bone, a site of low mechanical loading, ZIC1 mRNA expression was substantially lower. A strong correlation between ZIC1 expression and transcripts characteristically expressed in osteoblasts, osteocytes and osteoclasts was documented. We then investigated a possible association between the transcription factor and osteoporosis and BMD in transiliac bone biopsies from postmenopausal Norwegian women and found a negative correlation between total, spinal and hip BMD Z-scores, respectively, and ZIC1 expression. The correlation was even stronger when BMD was adjusted for BMI variation. Therefore, this is the first study to link skeletal epigenetic changes with mechanical stresses and associated increase in bone turnover. The study also demonstrates an association between ZIC1 expression, the degree of ZIC1 promoter methylation, bone turnover and BMD.
To corroborate ex vivo observations in human bone biopsies, we studied the effects of mechanical stresses (FSS) on the modulation of key osteogenic genes, ZIC1 mRNA expression and its promoter methylation in cultured cells. We asked whether mechanical stress could upregulate ZIC1 mRNA in isolated rat osteoprogenitor cells, whcih possess osteogenic potential and affect the osteogenic response. The transcriptional response of the cells to FSS in a two-armed study (mechanically strained vs. static cells) was carried out to rule out factors other than FSS affecting the osteogenic response of the cells. Following the application of mechanical strain, the expression of genes associated with osteogenic differentiation (Alp, Runx2 and Col1a1) and genes of interest that might play a role in mediating the effect of FSS (Taz, Gli3 and Zic1) were determined (Figure 3A–C). Immunocytochemical localization of ZIC1 and RUNX2 in rat osteoprogenitor cells subjected to FSS was further compared with controls and showed increased ZIC1 and RUNX2 levels. Intense staining for RUNX2 was found in the nuclei, and we also demonstrated a significant ZIC1 co-localisation (Figure 4A–F). The fluorescent intensity of ZIC1 and RUNX2 in the cells was found to be significantly positively correlated (Figure 5C). A comparative analysis of the expression of osteogenic marker genes and Zic1 in rat osteoprogenitor cells subjected to FSS with static controls suggested that FSS was essential for the induction, promotion and maintenance of an osteogenic phenotype in ROPs (rat osteoprogenitors). Furthermore, FSS leads to ZIC mRNA induction.
It was notable that ZIC1 promoter methylation was markedly distinct at 20 CpG sites between loaded LS and unloaded IB bone and presented a strong negative correlation with ZIC1 mRNA levels (Table 1A,B) (Figure 1A). The expression of ZIC1 was also strongly inversely correlated with age and BMI-adjusted BMD in postmenopausal women. This may be explained by ZIC1 as an important sensor molecule in coupling bone stress to biological responses regulating bone remodeling. ZIC1 may respond to increased load with enhanced gene transcription. Osteocytes/osteoblasts in osteoporotic subjects (OP) with thinner and fewer trabeculae endure enhanced load and possibly microfractures and respond to this by increased ZIC1 transcriptional activity, possibly as a result of hypo methylation. We attempted to corroborate this hypothesis with our experiments using osteogenic cell cultures.
The importance and relevance of DNA methylation of genes involved in bone metabolism are underscored by several papers demonstrating associations. For example, we have previously shown that DNA methylation of SOST in bone is associated with fracture risk PMID: 25155887. Additionally, as reviewed PMID: 33569383, PMID: 33921902 several osteogenic differentiation markers and genes involved in bone metabolism have been shown to be regulated by DNA methylation, including RUNX2, SP7(OSX), ALP, RANKL, and OPG.
It has been previously shown that a shear-stress-induced increase in T-cell factor/lymphoid enhancer factor transcriptional activity in osteogenic cells was abolished by ZIC1 silencing, implying that ZIC1 has an important role in shear flow mechanotransduction [1]. In order to establish whether observed changes in ZIC1 epigenetic promoter and accompanying increased ZIC1 mRNA expression were causally linked to mechanical loading, we subjected rat osteogenic cells and human osteoblast-like cells (SAOS2) to FSS.
SAOS2 and HSF cells were used to study ZIC1 promotor methylation before and after 12 h of fluid shear stress, conditions proven to activate the osteogenic response, but induced no statistical differences in either cell type. For the reasons unknown, the degree of basal ZIC1 promoter methylation differed markedly between the two cell types, as did ZIC1 expression. It is possible that FSS, over hours, may not represent a suitable in vitro model for DNA methylation caused by inductive stress. However, in vivo bone anabolic response can never be fully mimicked in in vitro cell culture experiments. This is particularly the case when it involves regulation of complex signal–transductions reactions as initiated by ZIC1 activation of osteocytes located in bone lacunae.
Despite considerable progress in identifying mechanosensing and mechanotransduction-signaling molecules and pathways, there is a lack of knowledge as to how the transmitted information from the membrane is converted into transcriptional changes. An earlier comparison [2] of global gene expression using microarrays (human Affymetrix gene chip) in human bone samples from sites that experience high and low levels of mechanical stress, lumbar spine and iliac crest, respectively, revealed that ZIC1 was significantly upregulated in the lumbar spine compared to the iliac crest. More recently, we showed that that the zinc finger protein of the cerebellum (ZIC1) plays an important role in shear flow mechanotransduction in osteocytes and osteoblasts, as well as activation of the primary cilia to mechanotransduction in bone [1]. Therefore, we investigated whether stress-related increased expression of ZIC1 was accompanied by epigenetic changes by comparing ZIC1 promoter DNA methylation in male iliac crest and lumbar spine. Trabecular bone was used, as it has an extensive surface representing the main site of bone metabolic activity. ZIC1 promoter methylation was significantly lower in the lumbar spine than in the iliac crest, and the extent of DNA demethylation was linearly negatively correlated with ZIC1 gene expression (Table 1A, Figure 1A).
The ZIC1 expression in human trabecular bone was positively correlated with osteoinductive bone morphogenetic protein 2 (BMP2) (Supplementary Figure S1A). The possible role of ZIC1 in promoting osteogenesis was underlined by its significant positive correlation with transcripts characteristically associated with osteogenic cells, osteocytes (SOST, DKK1 and PDPN) and osteoblasts (COL1A1, BGLAP, SPP1, CDH11 and PTHR1) and common to both cell types (RUNX2, IBSP and MEPE) (Figure 2 and Supplementary Figure S1). A few transcripts associated with osteoclast activity showed positive correlation (CTSK and ARC5) with ZIC1 expression, whereas OSCAR showed inverse correlation (Supplementary Figure S1B). The observations reported in this paper combined with those from earlier work establish that ZIC1 upregulation at sites of high stress and high bone turnover has an important role in mechanical-stress-induced increases in osteogenic differentiation and osteoblast and osteocyte activity. The present data suggest that this increase in ZIC1 expression in human trabecular bone was the result of reduced methylation at its promoter region.
Finally, the results also suggest that increased mechanical loading and stress present in the fewer and thinner trabeculae in osteoporotic bone with ensuing reduced ZIC1 promoter methylation is likely a compensatory mechanism aiming to increase ZIC1 transcription and bone formation. Additionally, the expression of SOST/sclerostin, an important negative regulator of the Wnt receptor anabolic pathway, is compensatory reduced at the mRNA and protein level in osteoporotics, accompanied by hypermethylation of its promotor [21]. Since ZIC1 expression is increased, whereas SOST expression is reduced at low BMD, a positive correlation between SOST and ZIC1 mRNA levels may appear paradoxical. However, the results make sense if both transcripts change their expression as compensatory mechanisms attempting to counteract the bone loss.

4. Material and Methods

4.1. Male Participants

All participants were Caucasian men (n = 13) from the northeast of England, U.K. The vertebral biopsies from these subjects were collected whilst they were undergoing spinal fusions or spinal decompression laminectomies of the lumbar spine (LS). The study was conducted as per ethical principles for medical research involving human subjects expressed in the World Medical Association Declaration of Helsinki. The study was commenced after ethical approval was obtained from the Institutional Review Board of Northumberland Local Research Ethics Committee, Blyth, Northumberland NE24 2AG (REC 09/05/2005; Reference: 04/Q0902/29). Informed written consent was obtained from all the study participants. All possible secondary causes of bone loss were excluded by relevant detailed medical history. Medical examination and a range of laboratory investigation are detailed in Supplementary Materials. Patients with a history of oral glucocorticoids, bone preservation agents, anticoagulants or anticonvulsants were excluded [2]. The mean ± SD (range) of anthropometric indices and BMD, given as areal density, is presented in Supplementary Table S1.

4.2. Female Participants

The female participants were a randomly selected subset (n = 57) from the ‘Osteogene Study’ cohort comprised of 84 Norwegian women (50–86 years) with a varying range of BMDs who were free of metabolic bone disease. All these subjects were recruited at the outpatient clinic of Lovisenberg Diaconal Hospital, Oslo (Supplementary Table S4).
Transiliac bone biopsies were obtained from ilia at the same location as described previously [16]. The power of the cohort has been previously reported [16,17].

4.3. RNA Purification and Gene Expression Analysis

Microarray analysis of total RNA on HG-U133 plus 2.0 chips (Affymetrix Santa Clara, CA, USA) was performed as previously described [2]. TaqMan gene expression analysis (real-time RT-PCR validation) of selected transcripts (n = 21) (as per the manufacturer’s instructions), together with the ribosomal protein L41 (RPL41) as internal standard [17], was used to control the reproducibility of the analysis on the chips. The rationale for selecting ribosomal protein L41 (RPL41) as an internal standard was that Affymetrix analyses had shown a similar magnitude of signal in the female bones (10,895 ± SD 315). The relevant female data were submitted to the EMBL-EBI (European Bioinformatics Institute) ArrayExpress repository (accession number: E-MEXP-1618); the male data can be accessed at the same site (accession number: E-MEXP-2219).

4.4. Bone Biopsies

To avoid possible anatomical-site-related variability, all bone biopsies were trabecular and were taken from a specific iliac site, namely two centimetres distal to the iliac crest and two centimetres posterior from the anterior superior iliac spine [17]. To avoid interoperator variability, all biopsies were collected under the supervision of one named senior surgeon, who was instructed to avoid attached tissue and muscle.
Likewise, male biopsies from the lumbar spine were also trabecular. These were primarily from the lamina process of lumbar vertebrae (L2, L3 or L4), whereas 19 biopsies from the vertebral lamina of the LS and 5 biopsies were from the iliac crest (IB).
Affymetrix ZIC1 signal values were correlated with signal values for the transcripts reflecting osteocyte (SOST, PDPN and MEPE), osteoblast (COL1A1, IBSP, SPARC, BGLAP and CDH11), osteoblast differentiation (BMP2 and RUNX2), and osteoclast activity (CALCR and OSCAR).

4.5. Global Methylation Analyses

Analysis of the DNA obtained from iliac bone biopsies was performed using an Infinium HumanMethylation450 BeadChip (Illumina) as per the manufacturer’s instructions. This allowed for detection of the methylation status of 485,000 individual CpGs. On average, there were 17 CpG sites per gene region, which were found to be distributed across the first exon, gene body, promoter, 5′UTR and 3′UTR. The quantitative measurement of the methylation for each CpG (β value) was obtained from the fluorescence data in BeadStudio (Illumina). The data were preprocessed using minfi and normalized using BMIQ [22,23].

4.6. ZIC1-Specific Methylation Analysis

The methylation pattern of 20 consecutive CpGs within the promoter region of ZIC1 was measured using pyrosequencing. DNA from female bone biopsies (~100 ng) and/or human osteoblast-like cells (SAOS2) (~1 µg) was bisulphite-treated and purified using an EpiTect bisulphite kit (Qiagen, Germantown, MD, USA). A 207 bp amplicon (chr3:147,126,936–147,127,142; GRCh37/hg19 assembly) was amplified from ~10 ng of bisulfite-treated DNA in 25 µL reactions using 1× PyroMark PCR Master Mix, 1× CoralLoad Concentrate (both from Qiagen) and 0.2 µM primers. Primers were designed using PyroMark Assay design 2.0 (Qiagen); forward primer: 5′-TGGTTTGTTAAAAGGGGATGTT-3′ and biotinylated reverse primer 5′-B-ACACCCTCCCCCCCTTAATAA-3′. Cycling conditions were an initial denaturation step at 95 °C for 15 min, followed by 50 cycles of 94 °C for 30 s, 56 °C for 30 s, 72 °C for 30 s and a final extension period at 72 °C for 10 min.
A total of 20 µL of PCR products was added to 40 µL binding buffer (Qiagen), 2 µL streptavidin–Sepharose high-performance beads (GE Healthcare, Chicago, IL, USA) and 18 µL water and stirred for 5–10 min at 14,000 rpm. Single-stranded biotinylated templates were isolated using PyroMark Vacuum Prep WorkStation (Qiagen) and dispensed onto a PyroMark Q24 plate containing 25 µL of 0.3 µM sequencing primer (5′-AAGAGTTTTATAATATTTGGGATTG-3′) and annealing buffer (Qiagen). The plates were incubated for 2 min at 80 °C and subsequently cooled at room temperature for at least 5 min. Sequencing was performed in a PyroMark Q24 instrument with PyroGold reagents (Qiagen). Results were analysed using PyroMark Q24 2.0.6 software, Qiagen, Germantown, MD, USA.

4.7. Data Analysis

Data analysis was performed by importing CEL files (raw gene expression data) into Partek Genomics Suite software (Partek, Inc., St. Louis, MO, USA). Generation of signal values and normalization was performed by the application of the robust multichip analysis (RMA) algorithm. A one-way analysis of variance (ANOVA) model was used for the expression comparisons of different groups. The results are given as fold change.
In accordance with the manufacturer’s recommendation, raw gene expression data were normalized and processed using the Affymetrix GCOS software module and MAS 5.0 in ArrayAssist (Stratagene, La Jolla, CA, USA), which generated a list of genes ‘present’ in a given sample. Following the compilation of a raw data set identification of differentially expressed genes was performed by ArrayAssist software. Statistical significance of microarrays was assessed by unpaired t-test, which identified differential gene expression. Normalized data were filtered on expression, with a cutoff value of two-fold change in expression (p ≤ 0.05). Gene ontology/enrichment analysis was then performed on the differentially expressed clusters.

4.8. Cell Lines and Culture Conditions

Immortalized human SAOS2 osteoblast-like cells (ATCC, Manassas, Virginia, USA) were cultured in 25 cm2 flasks and incubated at 37 °C in 5% CO2. Briefly, in each experiment, SAOS2 cells were treated with trypsin-EDTA solution (Sigma-Aldrich, Merck KGaA, Darmstadt, Germany), removed from the flask surface and grown in McCoy’s 5a medium (Sigma-Aldrich) supplemented with 10% foetal bovine serum (FBS). For primary human skin fibroblast (HSF) cells, Dulbecco’s modified Eagle’s medium (DMEM, Sigma-Aldrich) supplemented with 5% Calf Serum (CS), 100 μg/mL Streptomycin and 100 U/mL penicillin was used. Both cell lines were being subcultured when they reached 70–80% confluence.

4.9. mRNA Extraction and Quantification

Cells were grown in original 25 cm2 culture flasks until they reached 80% confluence. Normal medium was removed, and the flasks were washed with PBS to ensure that no residual media or growth factors were present. Cells were then lysed by directly adding 1 mL of Trizol reagent (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) in the flask, passed through the pipette several times to homogenize the lysates, transferred to RNAse-free Eppendorf tubes and left for 5 min to completely dissociate the nucleoprotein complexes. An amount of 200 μL of chloroform was then added to the tubes, shaken vigorously for 10–15 s and left for another 3 min. Tubes were then centrifuged at 13,000 rpm for 15 min, and the colourless upper aqueous phase containing RNA was removed from all tubes, transferred to new tubes and precipitated by mixing with an equal volume of isopropyl alcohol. Samples were then left for 10 min, centrifuged at 13,000 rpm for 15 min and supernatant was removed. Remaining pellets were washed once with 1 mL of 70% ethanol, mixed by vortexing and centrifuged at 13,000 rpm for 3 min. Ethanol was removed, and RNA pellets were left to dry completely for 5 min and finally redissolved by adding 50 μL of RNase-free water, passed through the pipette tip several times and incubated for 10 min at 60 °C. The absorbance of each sample was then measured at 260 nm by a spectrophotometer (NanoDrop, Thermo Fisher Scientific, Waltham, MA, USA) to calculate the concentration of mRNA in ng/μL. In addition, calculating the 260/230 nm and 260/280 nm absorbance ratios determined the integrity and purity of extracted RNA. Prior to each measurement, the spectrophotometer was calibrated with a blank (RNase-free water) at 260 nm to provide a more accurate reading.
The first-strand cDNA was synthesized from extracted mRNA by reverse transcription (RT). A total of 4 μg of total RNA was used as starting template for each reaction and included oligo dT, dNTP (10 mM). The volume of water added for each reaction was adjusted depending on the concentration and volume of mRNA added to give a final volume of 14 μL.

4.10. Immunofluorescence Detection of ZIC1

Primary rat osteoprogenitor cells were isolated from 3–4-day-old rat pup calvariae. The pups were sprayed and cleaned with 70% alcohol before the calvariae were dissected and washed with balanced Hanks solution, cut into small 2–3 mm wide pieces and digested with 1 mg/mL collagenase (type II) on a rocking bed. The supernatant was discarded, and the collagenase reaction was repeated for a further 30 min, following which the supernatant was isolated and kept as fraction 1. Calvariae were then washed with PBS that was added to fraction 1. In the next step, fraction 2 was collected by washing calvariae with PBS containing 4 mM EDTA for 15 min at 37 °C on a rocking bed. The calvariae were washed once more with Hanks solution and added to fraction 2. The collagenase reaction was repeated for 30 min, and the supernatant was kept as fraction 3. The calvariae were once again washed with Hanks solution, and this was added to fraction 3. In the final step, all the fractions were pooled and spun at 800× g (1000 rpm) for 3 min, and the resulting cell populations were cultured in T75 flasks for further experiments. The rat calvarial osteoblasts used were a mixed primary culture that contained cells at a number of points in the continuum from progenitors to mature cells. Previous characterization of the cells showed their response to osteoinduction, with increased levels of alkaline phosphatase and deposition of mineralised matrix, as well as appropriate changes in expression of alkaline phosphatase, collagen I, osteopontin and osteocalcin [24].
Cells were harvested and grown on sterile 13 mm glass cover slips (in a 24-well plate) overnight at 37 °C to let them settle down and adhere on the surface. Then, cells were briefly washed with PBS, and freshly prepared (<3 months) preheated (37 °C) 4% paraformaldehyde in PBS was added for 10 min to fix cells on the surface. Cover slips were thoroughly washed with PBS/0.1% Tween three times, and specimens were then incubated with 3% goat serum dissolved in PBS/0.1% Tween for approximately 30 min. Goat serum was used as a blocking serum to prevent non-specific binding of IgG, as secondary antibody was originally raised in goat. After a brief wash with PBS/0.1% Tween, 150 μL of rabbit polyclonal anti-ZIC1 antibody (diluted at 2.5 μg/mL in PBS, Abcam, Cambridge, UK) was added and incubated at room temperature for 2 h.
After the 2 h incubation time, cover slips were extensively washed with PBS/0.1% Tween three times (each wash for 15 min), and 2.5μL of anti-rabbit IgG (whole molecule, FITC conjugate, Sigma-Aldrich) was added per cover slip (in a well containing 500 μL PBS) and incubated at room temperature for 1 h in a dark chamber. After three extensive washes with PBS/0.1% Tween, cover slips were dried and mounted with Vectashield with DAPI (Vector) on microscope slides and examined on a microscope (DMLB, Leica, Wetzlar, Germany) with a fluorescent lamp (ebq 100 isolated, Leistungselektronik, Jena, Germany).

4.11. In Vitro Fluid Shear Stress Studies

To address the question of whether mechanical stress can upregulate ZIC1 in isolated rat osteoprogenitor cells and if so, whether ZIC1 was associated with enhanced osteogenesis (rat progenitors), we performed studies on freshly isolated primary rat osteoprogenitor cells. Cells were seeded on sterile 13 mm glass coverslips (in a 24-well plate) or in culture flasks incubated at 37 °C overnight and were then serum-starved for another 12 h. Serum-free medium of the culture flaks was then replaced with 5 ml of normal medium supplemented with 5.958 g/L HEPES and transferred to an orbital shaker incubator (Gallenkamp; Fisher Scientific, Loughbrough, UK) to apply shear stress and incubated for variable required time at 37 °C and 50 rpm.
The transcriptional response of the cells to fluid shear stress in a two-armed study (mechanically strained vs. static cells) was carried out to rule out factors other than mechanical stress affecting the osteogenic response of the cells. The magnitude of mechanical strain applied was restricted to 10–12 dynes/cm2 and calculated using the formula τ = a(ηρ(2πf)3)1/2, where a is the orbital radius of rotation of the shaker, ρ is the density of the culture medium (0.9973 g/mL), η is the viscosity of the medium (0.0101) and f is the frequency of rotation (rotation/s). Following the application of mechanical strain, the expression of genes associated with osteogenic differentiation (ALP, RUNX2 and COL1A1) and genes of interest that might play a role in mediating the effect of FSS (TAZ, GLI3 and ZIC1) were determined.
The immunostaining intensity (i.e., ZIC1 and RUNX2) was measured for a total of 20 cells per treatment using ImageJ 1.42q software (U.S. National Institutes of Health, Bethesda, MD, USA), for a total of 20 cells per treatment. A straight line of a of distance of 100 pixels was drawn along the long axis of each cell such that it passed through the middle of the nucleus and passed through the cytoplasm at the start and at the end. For each pixel for a given cell, the colour intensity was given as mean and SD, and it was plotted on a graph against the x-axis as the pixel distance.

4.12. Quantitative RT-PCR

Validation of the Affymetrix gene expression data was carried out for selected transcripts using TaqMan gene expression assays and the Applied Biosystems Prism 7900 HT sequence detection system. From each donor, total RNA (500 ng) was reverse-transcribed (Omniscript, Qiagen Ltd.), and resulting cDNA, which represented 2.5 ng total RNA, was the used in the PCR reaction. The analyses were performed in duplicate. The gene expression results are given as RQ (relative quantitation), and the relative changes of a given gene transcript were calculated using the 2(ΔΔC(T)) method. The ribosomal protein L41 (RPL41) served as an endogenous control; this was selected as internal standard because Affymetrix analyses had shown similar signal values in all samples (10,895 ± SD 315).

5. Conclusions

This study unravels the molecular sensory machinery and downstream effector molecules in osteoblasts and osteocytes, identifying ZIC1 as an important anabolic factor to induce bone formation during mechanical stress and microfractures. The present results suggest the clinical relevance of ZIC1 as a mechanosensory regulator to modulate osteogenic anabolic activity, which, after further study, may provide a potential therapeutic target.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23062957/s1.

Author Contributions

H.K.D. and K.M.G. conceived, designed and supervised the overall project. M.A.B., H.K.D. and J.J.W. designed and supervised the in vitro studies. H.K.D. and M.A.B. carried out analysis of the in vitro data. M.K.K., G.S. and O.K.O. performed the experiments. S.R., S.J.C., O.K.O. and H.K.D. carried out the analysis of the male and female transcriptome data. K.M.G., V.T.G., S.R., M.P., S.P.T. and H.K.D. carried out the registration of postmenopausal and male bone donors, obtained transiliac biopsies, and purified bone DNA and RNA. H.K.D. wrote the paper with main input from K.M.G., S.R. and S.P.T., along with all co-authors. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the National Osteoporosis Society (now known as the Royal Society for Osteoporosis) (# BH121971), EPSRC (EP/K036939/1). Knowledge Transfer Account from Lovisenberg Diaconal Hospital Oslo, South-Eastern Norway Regional Health Authority (Oslo University Hospital) (Project 29750104) and the European Union FP6 Project Osteogene (FP6—502491).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Northumberland Local Research Ethics Committee, Blyth, Northumberland NE24 2AG (REC 09/05/2005; Reference: 04/Q0902/29).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Data Availability Statement

The data for postmenopausal biopsies have been submitted to the European Bioinformatics Institute (EMBL-EBI; ArrayExpress repository, ID: E-MEXP-1618). The data on male biopsies are accessible through accession number E-MEXP-2219. The scripts for analysis are available at https://github.com/Bioinformatics-Support-Unit/python-scripts/tree/master/zic1 (accessed on 26 April 2016).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kalogeropoulos, M.; Varanasi, S.S.; Olstad, O.K.; Sanderson, P.; Gautvik, V.T.; Reppe, S.; Francis, R.M.; Gautvik, K.M.; Birch, M.A.; Datta, H.K. Zic1 transcription factor in bone: Neural developmental protein regulates mechanotransduction in osteocytes. FASEB J. 2010, 24, 2893–2903. [Google Scholar] [CrossRef] [PubMed]
  2. Varanasi, S.S.; Olstad, O.K.; Swan, D.C.; Sanderson, P.; Gautvik, V.T.; Reppe, S.; Francis, R.M.; Gautvik, K.M.; Datta, H.K. Skeletal site-related variation in human trabecular bone transcriptome and signaling. PLoS ONE 2010, 5, e10692. [Google Scholar] [CrossRef] [PubMed]
  3. Mansour, T.A.; Lucot, K.; Konopelski, S.E.; Dickinson, P.J.; Sturges, B.K.; Vernau, K.L.; Choi, S.; Stern, J.A.; Thomasy, S.M.; Döring, S. Whole genome variant association across 100 dogs identifies a frame shift mutation in DISHEVELLED 2 which contributes to Robinow-like syndrome in Bulldogs and related screw tail dog breeds. PLoS Genet. 2018, 14, e1007850. [Google Scholar] [CrossRef] [PubMed]
  4. Medina-Gomez, C.; Mullin, B.H.; Chesi, A.; Prijatelj, V.; Kemp, J.P.; Shochat-Carvalho, C.; Trajanoska, K.; Wang, C.; Joro, R.; Evans, T.E. Genome Wide Association Metanalysis of Skull Bone Mineral Density Identifies Loci Relevant for Osteoporosis and Craniosynostosis. medRxiv 2021. [Google Scholar] [CrossRef]
  5. Bonewald, L.F. The amazing osteocyte. J. Bone Miner. Res. 2011, 26, 229–238. [Google Scholar] [CrossRef]
  6. Bonewald, L. Osteocytes: A proposed multifunctional bone cell. J. Musculoskelet. Neuronal Interact. 2002, 2, 239–241. [Google Scholar]
  7. Dallas, S.L.; Prideaux, M.; Bonewald, L.F. The osteocyte: An endocrine cell… and more. Endocr. Rev. 2013, 34, 658–690. [Google Scholar] [CrossRef] [Green Version]
  8. Guo, D.; Bonewald, L.F. Advancing our understanding of osteocyte cell biology. Ther. Adv. Musculoskelet. Dis. 2009, 1, 87–96. [Google Scholar]
  9. Trumbore, D.; Heideger, W.; Beach, K. Electrical potential difference across bone membrane. Calcif. Tissue Int. 1980, 32, 159–168. [Google Scholar] [CrossRef]
  10. Bonewald, L.F.; Johnson, M.L. Osteocytes, mechanosensing and Wnt signaling. Bone 2008, 42, 606–615. [Google Scholar] [CrossRef] [Green Version]
  11. Burra, S.; Nicolella, D.P.; Francis, W.L.; Freitas, C.J.; Mueschke, N.J.; Poole, K.; Jiang, J.X. Dendritic processes of osteocytes are mechanotransducers that induce the opening of hemichannels. Proc. Natl. Acad. Sci. USA 2010, 107, 13648–13653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Datta, H.; Ng, W.; Walker, J.; Tuck, S.; Varanasi, S. The cell biology of bone metabolism. J. Clin. Pathol. 2008, 61, 577–587. [Google Scholar] [CrossRef] [PubMed]
  13. Teti, A. Bone development: Overview of bone cells and signaling. Curr. Osteoporos. Rep. 2011, 9, 264. [Google Scholar] [CrossRef] [PubMed]
  14. Zhong, J.; Chen, S.; Xue, M.; Du, Q.; Cai, J.; Jin, H.; Si, J.; Wang, L. ZIC1 modulates cell-cycle distributions and cell migration through regulation of sonic hedgehog, PI 3 K and MAPK signaling pathways in gastric cancer. BMC Cancer 2012, 12, 290. [Google Scholar] [CrossRef] [Green Version]
  15. Merzdorf, C.S.; Sive, H.L. The zic1 gene is an activator of Wnt signaling. Int. J. Dev. Biol. 2004, 50, 611–617. [Google Scholar] [CrossRef]
  16. Jemtland, R.; Holden, M.; Reppe, S.; Olstad, O.K.; Reinholt, F.P.; Gautvik, V.T.; Refvem, H.; Frigessi, A.; Houston, B.; Gautvik, K.M. Molecular disease map of bone characterizing the postmenopausal osteoporosis phenotype. J. Bone Miner. Res. 2011, 26, 1793–1801. [Google Scholar] [CrossRef]
  17. Reppe, S.; Refvem, H.; Gautvik, V.T.; Olstad, O.K.; Høvring, P.I.; Reinholt, F.P.; Holden, M.; Frigessi, A.; Jemtland, R.; Gautvik, K.M. Eight genes are highly associated with BMD variation in postmenopausal Caucasian women. Bone 2010, 46, 604–612. [Google Scholar] [CrossRef]
  18. Jacobs, C.R.; Temiyasathit, S.; Castillo, A.B. Osteocyte mechanobiology and pericellular mechanics. Annu. Rev. Biomed. Eng. 2010, 12, 369–400. [Google Scholar] [CrossRef]
  19. Hoffman, B.D.; Grashoff, C.; Schwartz, M.A. Dynamic molecular processes mediate cellular mechanotransduction. Nature 2011, 475, 316. [Google Scholar] [CrossRef]
  20. Liu, C.; Zhao, Y.; Cheung, W.-Y.; Gandhi, R.; Wang, L.; You, L. Effects of cyclic hydraulic pressure on osteocytes. Bone 2010, 46, 1449–1456. [Google Scholar] [CrossRef] [Green Version]
  21. Reppe, S.; Lien, T.G.; Hsu, Y.-H.; Gautvik, V.T.; Olstad, O.K.; Yu, R.; Bakke, H.G.; Lyle, R.; Kringen, M.K.; Glad, I.K. Distinct DNA methylation profiles in bone and blood of osteoporotic and healthy postmenopausal women. Epigenetics 2017, 12, 674–687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Teschendorff, A.E.; Marabita, F.; Lechner, M.; Bartlett, T.; Tegner, J.; Gomez-Cabrero, D.; Beck, S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 2012, 29, 189–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kirmizidis, G.; Birch, M.A. Microfabricated grooved substrates influence cell–cell communication and osteoblast differentiation in vitro. Tissue Eng. Part A 2008, 15, 1427–1436. [Google Scholar] [CrossRef]
Figure 1. Correlation observed between ZIC1 promotor methylation and ZIC1 expression. (A) Correlation between ZIC1 promotor methylation in CpG1 and ZIC1 expression in the male lumbar spine trabecular bone (Pearson r = −0.79, p = 2.4 × 10−6) (Table 1). (B) ZIC1 expression in the iliac bone of postmenopausal women showed a significant negative correlation with BMD. An inverse association between ZIC1 expression and total femoral BMD Z-score that had been adjusted for variation in body mass index (BMI) (r = −0.383; p = 0.0033) was observed. (C) A difference was observed in the methylation levels of ZIC1 promoter at CpG1 from iliac (IB) and lumbar regions (L3–L5) (expressed as %), skeletal sites of low and high mechanical stress, respectively.
Figure 1. Correlation observed between ZIC1 promotor methylation and ZIC1 expression. (A) Correlation between ZIC1 promotor methylation in CpG1 and ZIC1 expression in the male lumbar spine trabecular bone (Pearson r = −0.79, p = 2.4 × 10−6) (Table 1). (B) ZIC1 expression in the iliac bone of postmenopausal women showed a significant negative correlation with BMD. An inverse association between ZIC1 expression and total femoral BMD Z-score that had been adjusted for variation in body mass index (BMI) (r = −0.383; p = 0.0033) was observed. (C) A difference was observed in the methylation levels of ZIC1 promoter at CpG1 from iliac (IB) and lumbar regions (L3–L5) (expressed as %), skeletal sites of low and high mechanical stress, respectively.
Ijms 23 02957 g001
Figure 2. Affymetrix ZIC1 signal values were correlated with signal values for the transcripts reflecting osteoblast matrix-forming activity. Transcripts for (A) COL1A1, IBSP*, SPARC*, BGLAP*, CDH11* and osteoblast differentiation (BMP2* and RUNX2*) showed strong positive correlation with ZIC1 in the lumbar spine but had weaker or no correlation in IB. The PTHR1* transcript used to indicate osteoblast number also showed a strong positive correlation with ZIC1 expression. (B) ACP5 and CTSK** transcripts that reflect osteoclast activity showed strong positive correlation with ZIC1 expression. Osteoclast-associated transcripts CALCR** and OSCAR** showed inverse correlation with ZIC1 expression in the lumbar spine but weaker and no correlation, respectively, in the iliacus, probably due to fewer samples. The expression of osteocyte-associated transcripts (C) SOST, PDPN*** and MEPE*** at both high (LS) and low (IB) bone turnover sites showed strong correlation with ZIC1 mRNA. (Please see Supplementary Figure S1. A*, B** and C*** for correlations).
Figure 2. Affymetrix ZIC1 signal values were correlated with signal values for the transcripts reflecting osteoblast matrix-forming activity. Transcripts for (A) COL1A1, IBSP*, SPARC*, BGLAP*, CDH11* and osteoblast differentiation (BMP2* and RUNX2*) showed strong positive correlation with ZIC1 in the lumbar spine but had weaker or no correlation in IB. The PTHR1* transcript used to indicate osteoblast number also showed a strong positive correlation with ZIC1 expression. (B) ACP5 and CTSK** transcripts that reflect osteoclast activity showed strong positive correlation with ZIC1 expression. Osteoclast-associated transcripts CALCR** and OSCAR** showed inverse correlation with ZIC1 expression in the lumbar spine but weaker and no correlation, respectively, in the iliacus, probably due to fewer samples. The expression of osteocyte-associated transcripts (C) SOST, PDPN*** and MEPE*** at both high (LS) and low (IB) bone turnover sites showed strong correlation with ZIC1 mRNA. (Please see Supplementary Figure S1. A*, B** and C*** for correlations).
Ijms 23 02957 g002aIjms 23 02957 g002b
Figure 3. The effect of fluid shear stress (FSS) on osteogenic genes. (A) Rat osteoprogenitor cells were cultured in osteogenic medium and subjected to FSS (F). Control cells are labelled S (stationary). The cells subjected to FSS for three days (lane 3F) had induction of the key osteogenic marker genes Runx2, ALP and Col1a1 when compared to controls (lane C). (B) The induction was confirmed by qRT-PCR in the FSS and control stationary cells. (C) Cells cultured in normal non-osteogenic medium, subjected to similar magnitude of FSS, showed no evidence of Zic1 or Col1a1 induction over three days (lane 3F vs. 3S). However, after prolonged shear stress, there was an induction of Col1a1, as well as ZIC1, in these cells (lanes 5F, 5S, 7F and 7S).
Figure 3. The effect of fluid shear stress (FSS) on osteogenic genes. (A) Rat osteoprogenitor cells were cultured in osteogenic medium and subjected to FSS (F). Control cells are labelled S (stationary). The cells subjected to FSS for three days (lane 3F) had induction of the key osteogenic marker genes Runx2, ALP and Col1a1 when compared to controls (lane C). (B) The induction was confirmed by qRT-PCR in the FSS and control stationary cells. (C) Cells cultured in normal non-osteogenic medium, subjected to similar magnitude of FSS, showed no evidence of Zic1 or Col1a1 induction over three days (lane 3F vs. 3S). However, after prolonged shear stress, there was an induction of Col1a1, as well as ZIC1, in these cells (lanes 5F, 5S, 7F and 7S).
Ijms 23 02957 g003
Figure 4. Rat osteoprogenitors (ROPs) cultured in osteogenic medium and maintained as static controls were fixed on 25 mm coverslips and immunostained at the end of day 3 with (A) anti-RUNX2 (green fluorescence, plate 1) and (B) anti-ZIC1 (red fluorescence, plate 2) antibodies. (C) Plates 1 and 2 were merged together to identify co-localisation of RUNX2 and ZIC1 antibodies (Plate 3). In the representative field of vision, many nuclei showed nuclear translocation of RUNX2, suggesting relevant osteogenic activity in these cells. A few nuclei showed co-localisation of RUNX2 and ZIC1, and it appears that expression of RUNX2 (yellow) was accompanied by some ZIC1 (red/green). The effect of FSS on the intracellular distribution of RUNX2 and ZIC1 in ROPs is seen in (D) plate 4 and (E) plate 5, respectively. (F) Plates 4 and 5 were merged together to identify co-localisation of RUNX2 and ZIC1 antibodies (Plate 6).
Figure 4. Rat osteoprogenitors (ROPs) cultured in osteogenic medium and maintained as static controls were fixed on 25 mm coverslips and immunostained at the end of day 3 with (A) anti-RUNX2 (green fluorescence, plate 1) and (B) anti-ZIC1 (red fluorescence, plate 2) antibodies. (C) Plates 1 and 2 were merged together to identify co-localisation of RUNX2 and ZIC1 antibodies (Plate 3). In the representative field of vision, many nuclei showed nuclear translocation of RUNX2, suggesting relevant osteogenic activity in these cells. A few nuclei showed co-localisation of RUNX2 and ZIC1, and it appears that expression of RUNX2 (yellow) was accompanied by some ZIC1 (red/green). The effect of FSS on the intracellular distribution of RUNX2 and ZIC1 in ROPs is seen in (D) plate 4 and (E) plate 5, respectively. (F) Plates 4 and 5 were merged together to identify co-localisation of RUNX2 and ZIC1 antibodies (Plate 6).
Ijms 23 02957 g004
Figure 5. (A,B) The fluorescence intensity of RUNX2 and ZIC1 nuclei (expressed as mean, grey values) from ROPs cultured in osteogenic medium, maintained either as static control or subjected to FSS at the end of day 3. (C) The plot shows a direct relationship between translocation of RUNX2 and ZIC1 to the nucleus when subjected to FSS.
Figure 5. (A,B) The fluorescence intensity of RUNX2 and ZIC1 nuclei (expressed as mean, grey values) from ROPs cultured in osteogenic medium, maintained either as static control or subjected to FSS at the end of day 3. (C) The plot shows a direct relationship between translocation of RUNX2 and ZIC1 to the nucleus when subjected to FSS.
Ijms 23 02957 g005aIjms 23 02957 g005b
Table 1. (A) Comparison of ZIC1 promoter methylation in trabecular bone biopsies from male iliac crest (IB) and the lumbar spine (L5), as well as the correlation between the degree of ZIC1 promoter methylation at CpG positions and ZIC1 expression. (B) Differential expression of ZIC1. The mean relative expression levels of ZIC1 mRNA in the lumbar spine vertebrae and iliac crest trabecular bone from men.
Table 1. (A) Comparison of ZIC1 promoter methylation in trabecular bone biopsies from male iliac crest (IB) and the lumbar spine (L5), as well as the correlation between the degree of ZIC1 promoter methylation at CpG positions and ZIC1 expression. (B) Differential expression of ZIC1. The mean relative expression levels of ZIC1 mRNA in the lumbar spine vertebrae and iliac crest trabecular bone from men.
(A)
ZIC1
Methylation Site
Pos. 1Pos. 2Pos. 3Pos. 4Pos. 5Pos. 6Pos. 7Pos. 8Pos. 9Pos. 10Pos. 11Pos. 12Pos. 13Pos. 14Pos. 15Pos. 17Pos. 18Pos. 19Pos. 20Mean
Mean IC
(n = 7)
49.223.72516.223.518.619.911.711.324.512.59.413.816.19.59.98.358.816.7
Mean L5
(n = 7)
31.815.418.511.715.312.612.77.57.516.88.66.610.9116.36.964.15.311.3
p-value
(L5 vs. IB)
9 × 10−40.0030.0060.0090.0050.0090.0060.0110.010.010.0150.0110.0420.0070.0040.0040.0190.0550.0320.004
Pearson r (ZIC1 vs. methylation)−0.719−0.71−0.56−0.47−0.71−0.64−0.74−0.56−0.59−0.72−0.62−0.455−0.409−0.677−0.682−0.316−0.346−0.327−0.343−0.67
(B)
ZIC1Probe set IDIliac CrestL3L4L5ANOVA p-Value
206373_at188652858565.56 × 10−10
234716_at373939520.853
236896_at81814230.0003
Table 2. Correlation between ZIC1 expression level in iliac bone of postmenopausal women with BMD, lean body and fat mass (n = 57).
Table 2. Correlation between ZIC1 expression level in iliac bone of postmenopausal women with BMD, lean body and fat mass (n = 57).
Pearson rp-Value
Total hip Z-score, BMI adj−0.3830.003
Total hip Z-score−0.3690.005
Femoral neck Z-score, BMI adj−0.3310.012
Femoral neck Z-score−0.3230.011
L1–L4 Z-score, BMI adj−0.3370.01
L1–L4 Z-score *−0.3370.01
Total body Z-score, BMI adj *−0.3980.003
Total body Z-score−0.3950.004
Lean mass (g)0.118NS
Fat mass (g)0.108NS
Probe set for ZIC1 expression level: 206373_at * n = 51.
Table 3. Differential ZIC1 DNA methylation levels observed between osteoporotic (OP) and healthy subjects at iliac bone site.
Table 3. Differential ZIC1 DNA methylation levels observed between osteoporotic (OP) and healthy subjects at iliac bone site.
I.D.Mean Healthy
(n = 40)
±SDMean OP
(n = 29)
±SDp-Value
cg172803460.240.040.300.070.0004 *
cg012275370.140.030.170.060.0043
cg166366710.110.030.140.040.0116
cg234496960.340.040.370.040.0138
cg057166710.100.020.110.030.0211
cg150022940.070.010.080.030.0224
cg053715780.190.030.210.050.0257
cg025197510.190.040.230.080.0294
cg147509480.190.040.220.060.0344
cg162096640.180.030.190.040.0413
cg144566830.260.030.280.050.0480
cg190291810.250.060.280.060.0499
This table shows CpGs with DNA methylation levels differing between OP and healthy (nominally significant). p-values were calculated using Student’s t-test. * Shows significance after adjustment for multiple testing (p = 0.00041* 20 ZIC1 CpGs = 0.008).
Table 4. (A) Comparison between ZIC1 promoter DNA methylation at five positions in primary human skin fibroblasts (HSFs) and osteoblast-like human osteosarcoma cell line (SAOS2 cells). (B) Effect of FSS on modulation of gene transcription/expression in osteoblast-like human osteosarcoma cell line (SaOS2 cells). (C) Effect of FSS on transcription in primary human skin fibroblasts (HSFs).
Table 4. (A) Comparison between ZIC1 promoter DNA methylation at five positions in primary human skin fibroblasts (HSFs) and osteoblast-like human osteosarcoma cell line (SAOS2 cells). (B) Effect of FSS on modulation of gene transcription/expression in osteoblast-like human osteosarcoma cell line (SaOS2 cells). (C) Effect of FSS on transcription in primary human skin fibroblasts (HSFs).
(A)
Cell TypePos. 1Pos. 2Pos. 3Pos. 4Pos. 5
HSFs (Mean ± SD, n = 3)43.99 ± 3.628.89 ± 5.0326.72 ± 12.4211.73 ± 8.796.57 ± 3.54
SAOS (Mean ± SD, n = 3)3.34 ± 0.621.88 ± 0.162.55 ± 1.121.15 ± 0.1470.96 ± 0.38
(B)
GeneStatic
(Delta Ct Mean ± SD)
(n = 3)
FSS
(Delta Ct Mean ± SD)
(n = 3)
Fold Change (ln)
FSS vs. Static
p-Value
RUNX23.787 ± 0.0813.622 ± 0.0111.120.0049
ZIC16.983 ± 0.0346.694 ± 0.0961.22NS
SP76.088 ± 0.0895.517 ± 0.0581.490.0003
KIF3A6.877 ± 0.0376.608 ± 0.0851.200.0045
GLI36.000 ± 0.0095.979 ± 0.1111.01NS
GLI27.240 ± 0.0087.166 ± 0.0541.05NS
GLI19.617 ± 0.1079.205 ± 0.1691.330.0176
PTCH16.878 ± 0.0176.890 ± 0.170−1.01NS
COL1A1−1.645 ± 0.056−1.571 ± 0.174−1.05NS
ALPL0.530 ± 0.0720.546 ± 0.209−1.01NS
SMO6.295 ± 0.0176.369 ± 0.238−1.05NS
(C)
GeneStatic
(Delta Ct Mean ± SD)
(n = 3)
FSS
(Delta Ct Mean ± SD)
(n = 3)
Fold Change (ln)
FSS vs. Static
p-Value
RUNX27.650 ± 0.2667.341 ± 0.2271.14NS
ZIC110.391 ± 0.22710.204 ± 0.1921.22NS
SP712.455 ± 0.45212.331 ± 0.9431.09NS
KIF3A7.252 ± 0.1066.750 ± 0.0261.420.0006
GLI36.7434 ± 0.1556.410 ± 0.0451.260.0174
GLI26.966 ± 0.2456.624 ± 0.0221.27NS
GLI112.1314 ± 0.22911.181 ± 0.3231.930.0098
PTCH19.602 ± 0.1449.597 ± 0.0941.0NS
COL1A1−3.255 ± 0.189−3.0435 ± 0.192−1.16NS
ALPL13.475 ± 0.58511.475 ± 0135−4.000.0025
SMO8.121 ± 0.1657.666 ± 0.119−1.370.0128
Numbers indicate methylation levels expressed as percentage (%).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Datta, H.K.; Kringen, M.K.; Tuck, S.P.; Salpingidou, G.; Olstad, O.K.; Gautvik, K.M.; Cockell, S.J.; Gautvik, V.T.; Prediger, M.; Wu, J.J.; et al. Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis. Int. J. Mol. Sci. 2022, 23, 2957. https://doi.org/10.3390/ijms23062957

AMA Style

Datta HK, Kringen MK, Tuck SP, Salpingidou G, Olstad OK, Gautvik KM, Cockell SJ, Gautvik VT, Prediger M, Wu JJ, et al. Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis. International Journal of Molecular Sciences. 2022; 23(6):2957. https://doi.org/10.3390/ijms23062957

Chicago/Turabian Style

Datta, Harish K., Marianne K. Kringen, Stephen P. Tuck, Georgia Salpingidou, Ole K. Olstad, Kaare M. Gautvik, Simon J. Cockell, Vigdis T. Gautvik, Michael Prediger, Jun Jie Wu, and et al. 2022. "Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis" International Journal of Molecular Sciences 23, no. 6: 2957. https://doi.org/10.3390/ijms23062957

APA Style

Datta, H. K., Kringen, M. K., Tuck, S. P., Salpingidou, G., Olstad, O. K., Gautvik, K. M., Cockell, S. J., Gautvik, V. T., Prediger, M., Wu, J. J., Birch, M. A., & Reppe, S. (2022). Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis. International Journal of Molecular Sciences, 23(6), 2957. https://doi.org/10.3390/ijms23062957

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