[Frontiers in Bioscience, Landmark, 22, 783-794, January 1, 2017]

Blockage of hemichannels alters gene expression in osteocytes in a high magneto-gravitational environment

Huiyun Xu 1 , Dandan Ning 1 , Dezhi Zhao 1 , Yunhe Chen 1 , Dongdong Zhao 1 , Sumin Gu 2 , Jean X. Jiang 2 , Peng Shang 1

1Key Laboratory for Space Biosciences and Biotechnology, School of Life Sciences, Northwestern Polytechnical University, Xi’an, Shaanxi Province, 710072, China, 2Department of Biochemistry, University of Texas Health Science Center at San Antonio, Texas, 78229, USA

TABLE OF CONTENTS

1. Abstract
2. Introduction
3. Materials and methods
    3.1. Exposure to HMGE
    3.2. Cell culture
    3.3. RNA extraction and labeling
    3.4. Array hybridization and scanning
    3.5. Microarray data analysis
    3.6. Real time PCR
4. Results
    4.1. HMGE and hemichannel inhibition affected gene expression profile of MLO-Y4 cells
    4.2. Osteocytic hemichannels regulate signaling pathways and biological processes
    4.3. Impairment of Cx43 hemichannels greatly affected MAPK signaling, apoptosis, cell cycle and ubiquitin-mediated proteolysis pathways
    4.4. Real-time PCR validated the genes responsive to Cx43 hemichannel inhibition
5. Discussion
6. Acknowledgement
7. References

1. ABSTRACT

Osteocytes, the most abundant cells in bone, are highly responsive to external environmental changes. We tested how Cx43 hemichannels which mediate the exchange of small molecules between cells and extracellular environment impact genome wide gene expression under conditions of abnormal gravity and magnetic field. To this end, we subjected osteocytic MLO-Y4 cells to a high magneto-gravitational environment and used microarray to examine global gene expression and a specific blocking antibody was used to assess the role of Cx43 hemichannels. While 3 hr exposure to abnormal gravity and magnetic field had relatively minor effects on global gene expression, blocking hemichannels significantly impacted the expression of a number of genes which are involved in cell viability, apoptosis, mineral absorption, protein absorption and digestion, and focal adhesion. Also, blocking of hemichannels enriched genes in multiple signaling pathways which are enaged by TGF-beta, Jak-STAT and VEGF. These results show the role of connexin hemichannels in bone cells in high magneto-gravitational environments.

2. INTRODUCTION

During space flight, astronauts experience a variety of physiological changes. Among them, the reduction of bone remodeling is commonly observed (1). Previous studies have been undertaken to investigate the causes of the bone loss induced by microgravity in space; however the underlying mechanism still remains largely elusive (2).

In more recent years, numerous evidence suggest that osteocyte is a primary responser to the changes of mechanical environment in bone tissue (3-5). Several cellular components such as cytoskeleton and membrane channels are shown to be involved in biological responses of osteocytes to mechanical loading (6-8). As the most abundant connexin subtype in bone tissue in forming gap junctions and hemichannels, connexin 43 (Cx43) plays a key role in sensing and transmitting the mechanical signal to chemical signals in bone cells (9,10). Gap junctions and hemichannels permit the passage of small molecules below 1.2. kDa between adjacent cells and the cell and the extracellular environment, respectively (11). Earlier studies have shown that Cx43-deficiency mice desensitize bone to the effects of mechanical disuse induced by hindlimb suspension (12,13). Also, our previous results have shown that parabolic flight decreases the expression of Cx43(14). However, there is limited knowledge regarding how weightlessness and hemichanels influence global gene expression pattern and regulatory networks in osteocytes. Due to the technical challenges associated with conducting experiments in the outer space, many studies have been conducted using simulated equipment on the ground. In our labortoary, high magneto-gravitational environment (HMGE) has been used to simulate the effects of weightlessness on mammalian cells, plant callus, silkworm eggs and microorganism (15). The data generated from our previous studies have demonstrated that the simulator is an invaluable tool which can be used to mimic weightless environment in the outer space (16-21).

In this study, microarray analysis was conducted to identify the alterations of global gene expression profile of osteocytic MLO-Y4 cells under HMGE. Furthermore, the effects of Cx43 hemichannels on gene expression were investigated by using a potent, Cx43 hemichannel-blocking antibody Cx43(E2) (22). The specificity ad potency of this antibody on inhibition of Cx43 hemichannels has been demonstrated by us and other groups (23-27).

3. MATERIALS AND METHODS

3.1. Exposure to HMGE

The superconducting magnet used in this study was manufactured by Japan Superconductor Technology, Inc (JASTEC, Seishin, Japan), and several accessories used for cell experiments were constructed in our laboratory, including temperature control system, object stage, gas control system, and observation system. Thediamagnetic materials were placed into thebore of the superconducting magnet. In the center ofthe bore the magnetic field strength is maximal (16 Tesla (T)). Here the magnetic field gradient is zero and the magnetic force is also zero. The diamagnetic materials are only affected by gravity and corresponding position called 1g (16T). At a position where the direction of magnet force and gravity are opposite and the magnitude is equal, the magnet force can balance out the gravity, leading to stable levitation. This position is called µg (hypo-gravity, diamagnetic levitation, 12T). For the same reason, there is a position called 2g (hyper-gravity, 12T) where the direction and magnitude of magnet force and gravity are the same. Figure 1. Diagram illustrates three positions of different magnetic force and gravity in the center of superconducting magnet. For external control (C), the cells were cultured under the identical culture condition in an incubator, but away from the magnet to avoid the influence of magnetic field. We obtained pure gravity effects by comparing the results at 2g and µg, because the magnetic fields of the two positions were identical. Magnetic effects were obtained by comparing the results at 1g and external control, by which the gravity was identical.

3.2. Cell culture

The MLO-Y4 cell line was provided by Dr.Lynda Bonewald (University of Missouri, Kansas City, MO, USA) (28). The cells were cultured in MEM (Gibco, Paisley, UK) with 2.5 percent fetal bovine serum (FBS) and 2.5 percent calf serum including 1 percent penicillin/streptomycin. Cells (5×105 cells/dish) in the logarithmic phase were seeded in 35 mm culture dish (Nunc, Inc. Roskilde, Danemark) for 24 hrs, and then treated by irrelevant rabbit IgG antibody (CWBIO, Beijing, China, 1:50 dilution) in an incubator at 37°C with 5 percent CO2 for 30 mins, then separated into 4 groups: control (C), 1g, 2g and μg groups (Table 1). The cells of C group were cultured in the same incubator for another 3 hrs, and at the meantime the cells of 1g, 2g and μg groups were placed on sample holder and cultured on the HMGE platform with temperature controlled at 37±0.5.°C and 5percent CO2 for 3 hrs (15). For antibody treated group, Cx43(E2) antibody (1:50 dilution, 28µg/ml) was added into the media in 35mm culture dish for 30 min, the identical culture condition with previous study (22), and then separated into 4 groups: E2-C, E2-1g, E2-2g and E2-μg groups (Table 1). Cells of E2-1g, E2-2g and E2-μg groups were cultured in the given positions on HMGE platform, while cells of E2-C group were cultured in the same incubator with C group.

3.3. RNA extraction and labeling

Total RNA was extracted using TRIzol Reagent (Life technologies, Carlsbad, CA, USA) following the manufacturer’s instructions, integrity of isolated RNA was examined by an Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, USA). Total RNA was further purified by RNeasy micro kit and RAase-Free DNase Set (QIAGEN, GmBH, Germany). Then total RNA was amplified, labeled and purified using Ambion WT Expression kit (Ambion, USA) and Gene chip WT Terminal labeling Kit (Affymetrix, Santa Clara, CA, USA) following the manufacturer’s instructions to obtain biotin-labeled cDNA.

3.4. Array hybridization and scanning

The whole genome transcriptional profile was determined by 2 repeats of independent hybridization assay using Affymetrix Mouse Gene 1.0 ST Microarrays. Array hybridization and wash were performed using GeneChip Hybridization, Wash and Stain Kit in Hybridization Oven 645 and Fluidics Station 450 based on the manufacturer’s instructions (Affymetrix). Slides were scanned by GeneChip Scanner 3000 (Affymetrix) operated by Command Console Software 3.1 (Affymetrix) with default settings. Raw data were normalized by Gene Spring Software 11.0 (Agilent Technologies).

3.5. Microarray data analysis

SBC Analysis System (SAS, Shanghai Biotechnology Corporation, shanghai, CN) was performed to filter differentially expressed genes (DEGs) under different conditions, and the fold changes of any two groups more than 1.5 or less than 0.65 were selected. The Gene Ontology (GO) project was used to provide controlled codes to describe gene and gene product attributes in any organism (http://www.geneontology.org). KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway (http://www.genome.jp/kegg/) analysis was performed to cluster the DEGs according to biological processes. The p-value (EASE-score, Fisher p value or Hypergeometric pvalue) was used to denot the significance of the pathway correlated to the conditions. The network of pathway (Path-Act-Net) was built according to the relationship among the changed pathways from KEGG database, and the network of gene (Gene-Act-Network) was built according to the relationship among the DEGs. The clustered heatmap was presented using gplots package of R programming tools.

3.6. Real time PCR

Primers were designed and synthesized (Sangon, Shanghai, China) based on the gene sequences available in GenBank. These primer sequences are shown in Table 2. The total RNA was extracted using Trizol (Invitrogen, Carlsbad, CA, USA). Real-time relative quantitative RT-PCR analysis was performed with SYBR® Premix Ex Taq™ II based on the manufacturer’s protocol (TaKaRa, Dalian, China). Each sample was measured in triplicate and normalized to GAPDH. The 2-ΔΔCt method was used to estimate the relative change in gene expression. The differences between μg (microgravity) vs (versus) control, μg vs 2g, 1g vs control, and antibody-treated groups vs control groups were statistically analyzed by one-way ANOVA of GraphPad Prism. Pvalue of less than 0.1 was regarded as being significant.

4. RESULTS

4.1. HMGE and hemichannel inhibition affected gene expression profile of MLO-Y4 cells

The genome-wide gene expression analysis using microarray (GEO databases online: GSE85968) was performed to obtain a list of DEGs after exposure to HMGE and hemichannel inhibition for 3 hrs in MLO-Y4 cells (Table 3). Among a total of 35556 detected genes, only 4 genes showed significant alteration in the 2g position, compared with µg position (Figure 2A). These genes were Slc30a1, Zfp942, AA987161 and an unknown gene. The results showed that gravity exposure for 3 hrs did not significantly affect overall gene expression profile. However, incubation with Cx43(E2) antibody (E2-2g vs E2-µg) altered 56 genes (46 upregulated and 10 downregulated). Interestingly, Slc30a1, Zfp942 and AA987161 and unknown genes altered under HMGE condition were not changed by Cx43(E2) antibody. These data indicate that the responsive genes to HMGE are unlikely regulated by Cx43 hemichannels. In 1g position, 49 genes were significantly changed (4 upregulated and 45 downregulated) compared with the external control position with normal gravity. After the treatment with Cx43(E2) antibody (E2-1g vs E2-c), 36 genes changed (3 upregulated and 33 downregulated). As described in detail, magnetic effect is more predominant at 1g position. Among the DEGs affected by magnetic field, 9 genes changed independent of the antibody treatment (Zfp942, Zfp943, Dzip3, 4932438A13Rik, Rock1, Dennd4a, Stag2, Setd2 and Myo9a), while the other 40 genes were not changed when compared E2-1g to E2-c, therefore the changes of the 40 genes were eliminated by antibody treatment. GO analysis showed that the 40 genes were mainly involved in various cellular functions including RNA and protein synthesis/degradation (Twistnb, Cul4b, Setd2, Thoc2, Eea1), bone metabolism (Arhgap5, Trpm7, Angpt1), mitogen-activated protein kinase (MAPK) signaling pathway (Taok1) and inositol phosphate metabolism (Pik3c2a). For the 49 DEGs affected by magnetic effect, GO analysis showed that there were 40 genes pertained to cellular structural components, and 30 and 22 genes were related to cellular and metabolic processes, respectively. Also, 28 and 13 genes involved in the categories of binding and localization, respectively (Figure 2B).

Interestingly, 801 genes were significantly changed, when comparing Cx43(E2) antibody-treated groups with the control groups (E2-c vs Con-c), of which 156 genes were up-regulated and 654 genes were down-regulated. GO analysis showed 629 DEGs enriched in cellular components and 428 genes in binding, 351 and 157 genes were related to metabolic and developmental processes, respectively. Also 273, 153 and 146 genes involved in catalytic activity, localization and biological regulation, respectively (Figure 2C).

The gene expression heatmap and cluster analysis on all of DEGs also confirmed that Cx43 hemichannel was a principal factor in regulating gene expression patterns, more DEGs were detected after blockage of hemichannel (comparing E2 groups with control groups), while gravity and magnetic play relatively minor roles here (Figure 3A).

4.2. Osteocytic hemichannels regulate signaling pathways and biological processes

KEGG pathway analysis was performed to cluster the DEGs to signaling pathways and biological processes (P value < 0.5). Because only 4 or 49 genes changed the expression after gravity or magnetic treatment, we focused more on the genes altered by Cx43(E2) antibody treatment. Table 4 showed 20 biological processes and pathways with the minimum P value influenced by blocking Cx43 hemichannels. Among them, protein/RNA synthesis and digestion were clustered, including ribosome biogenesis, aminoacyl-tRNA biosynthesis, ubiquitin mediated proteolysis, amino acid metabolism (glycine, serine and threonine metabolism and valine, leucine and isoleucine degradation), RNA transport and degradation, etc. Also, we found enriched KEGG pathways related to bone metabolism, such as mineral absorption, osteoclast differentiation and other related signal pathways including MAPK, transforming growth factor beta (TGF-beta), vascular endothelial growth factor (VEGF), Janus kinase-signal transducer and activator of transcription (Jak-STAT) and certain cellular processes including apoptosis, protein digestion and absorption and focal adhesion, etc.

The heatmap (Figure 3B) also showed the alterations concerning bone metabolism-related genes after blocking Cx43 hemichannels compared with the control groups including Tgfbr1, IL6, Jak2, Ctsk, Csf, Xiap and Mapkapk3, etc.

4.3. Impairment of Cx43 hemichannels greatly affected MAPK signaling, apoptosis, cell cycle and ubiquitin-mediated proteolysis pathways

The network of genes (Gene-Act-Network) and pathways (Path-Act-Net) was built according to the relationship among the DEGs and pathways from the KEGG database. There were direct (i.e. activation, binding or inhibition) or indirect interaction between changed genes related with bone metabolism in response to the inhibition of Cx43 hemichannels. The more one gene interacts with others, the more important the gene is in the Gene-Act network. That means this gene regulate or are regulated by more other genes. From Figure 4, the most central 15 genes were Pdgfra, Kras, Jak2, Figf, Arhgap5, Csf2rb, IL6, Xiap, Col3a1, Rock1, Csf3, Tgfbr1, IL8, Ccl3 and IL12rb1.

Path-Act-Net analysis showed that MAPK signaling pathway, apoptosis, cell cycle and ubiquitin-mediated proteolysis were the mostly affected cellular processes by Cx43 hemichannel inhibition (Figure 5).

4.4. Real-time PCR validated the genes responsive to Cx43 hemichannel inhibition

Real-time PCR was used to verify the results of the gene microarray. Figure 6 showed the relative expression of 15 central genes in Gene-Act-Net analysis, including Pdgfra, Kras, Jak2, Figf, Arhgap5, Csf2rb, IL6, Xiap, Col3a1, Rock1, Csf3, Tgfbr1, IL8, Ccl3 and IL12rb1. These genes showed significant difference after Cx43(E2) treatment, which were consistent with the changes of gene expression pattern obtained by gene microarray analysis.

5. DISCUSSION

Our previous results have shown that HMGE treatment for 24 hrs affects osteoblast gene expression profile and cytoskeleton-related gene expression (16, 20). 2595 genes were altered in 1g position vs normal gravity control (magnetic effect), while 1729 genes were altered in 2g position vs µg position (gravity effect) (29). Also, Hammer et al. (2009) have shown that 48 hrs of gravitational stress lead to up-or down-regulation of hundreds of genes in MC3T3-E1 osteoblast cells (30). In osteocytic MLO-Y4 cells, we have shown earlier that 48 hrs of HMGE exposure affect the expression of energy metabolism related genes, including enzyme, peptide hormone, G-protein coupled receptors and glucose metabolic process (31). To further understand the global gene expression profile of osteocytes under short duration of abnormal magnetic/gravity conditions, in this study, we subjected MLO-Y4 cells to HMGE for 3 hrs. Our results showed that only 4 genes expression was altered in 2g vs µg position (gravity effect). It is likely that 3 hrs of exposure to abnormal gravity is not long enough to cause major changes in gene expression of MLO-Y4 cells. Also, 3 hrs of exposure to high magnetic field (16T) only affect the expression of 49 genes (1g vs external 1g control) with down-reguation of 45 genes, which suggests a minor influence of magnetic field on cellular and metabolic processes. These results suggest that 3 hrs of HMGE treatment do not exert a significant effect on cellular function of MLO-Y4 cells.

The treatment of Cx43(E2) antibodies which specifically inhibits the opening of hemichannels (22,32) has a great impact on the gene expression of MLO-Y4 cells (E2 groups vs control groups, 801 DEGs). This suggests a quick response of the cell to the inhibition of hemichannels regardless of HGME treatment. Cx43 hemichannels in osteocytes have been shown to mediate the release of metabolic molecules in response to stress stimulation, including prostaglandins (PGs), ATP, etc (33-35). Blocking of hemichannels inhibits the exchange of molecules between cells and extracellular matrix, and the response of intracellular signal transduction, which may ultimately affect the gene expression and cellular process. Interestingly, except 2g group with similar up-and down-regulated genes, the other three groups (E2-C vs C, E2-1g vs 1g, E2-µg vs µg) consistently showed significantly more down-regulation genes than up-regulation. This data suggests that hemichannels primarily sever as a positive regulator by enhancing gene expression, and inhibition of the hemichannels has an adverse effect on gene expression in MLO-Y4 cells.

In this study KEGG pathways analysis showed that blocking of hemichannels led to the enrichment of DEGs on signaling pathways of TGF-beta, Jak-STAT and VEGF, and cellular processes of apoptosis, mineral absorption, protein digestion and absorption, and focal adhesion. These pathways and processes have been thought to involve in bone metabolism, especially osteocyte viability and osteoclast formation/differentiation.

Published studies have shown that osteocyte apoptosis is a prelude of osteoclast recruitment and followed by bone loss (36,37). Mechanical stimulation increases the production of factors required for cell survival, and thus inhibits osteocytes apoptosis (38). Our previous studies have demonstrated that Cx43 hemichannels play a predominant role in protecting osteocytes againist cell apoptosis and death. We show that the apoptotic osteocytes increased in Cx43(Δ130-136) transgenic mice with the impairment of both hemichannels and gap junctions in osteocytes, but not in Cx43(R76W) mice with the impairment of only gap junctions (39). Also, Cx43 hemichannels protect osteocytes against oxidative stress-induced cell death (25). The mechanism is likely a result of the release of small molecules and/or the activation of other signaling pathways (40). In this study, based on enriched KEGG pathways, we identified several DEGs involved in cellular apoptosis (Csf2rb2, Il1r1, Csf2rb, Casp6, Xiap and Atm) after Cx43 hemichannel blockage. Moreover, Gene-Act-Net analysis showed the down-regulation of some apoptosis related genes, such as Jak2, IL6, Csf2rb2, Csf2rb, Csf2 and Csf3, etc (41,42). Furthermore, Gene-Act-Net analysis showed up-regulation of genes related to collagen production (Col6a2, Col6a1 and Col3a1), which may help increase cell-matrix adhesion. And Path-Act-Net analysis also showed that cell viability related pathways were greatly altered after Cx43 hemichannel inhibition, such as MAPK signaling pathway, apoptosis, cell cycle and ubiquitin-mediated proteolysis. All these gene-related changes support the notion that functional Cx43 hemichannels are crucial to survival and apoptosis status of osteocytes. Furthermore, the gene profile data provided by this study point to potential, related downstream pathways for further investigation.

Taken together, this study showed that short exposure of osteocytes to abnormal gravity and high magnetic field had minor effects on gene expression profile of osteocytic MLO-Y4 cells. Also, the alterations of 40 out of 49 genes under high magnetic conditions were eliminated by Cx43 hemichannel inibition. This result suggests that Cx43 hemichannels mediate biological responses of osteocytes to short term exposure of high magnetic field. Importantly, inhibition of Cx43 hemichannels significantly changed the expression of over 800 genes in MLO-Y4 cells. Our data provide further evidence that Cx43 hemichannels play an important role in regulating global gene expression profile of osteocytes, in particular the genes involved in viability and apoptosis of osteocytes. These results provide important information and guidance for our further understanding the role of connexin hemichannels in bone cells in future.

6. ACKNOWLEDGEMENTS

We thank the grants from the National Natural Science Foundation of China (81472090 & 31328016 to HX, JXJ and PS), Natural Science Basic Research Plan in Shaanxi Province (2015JM8443 to HX), the Fundamental Research Funds for the Central Universities (3102014JKY15012 & 3102016ZY036 to HX), and National Institute of Health grant CA196214 and Welch Foundation grant AQ-1507 to JXJ.

7. REFERENCES

1. A. LeBlanc, V. Schneider, L. Shackelford, S. West, V. Oganov, A. Bakulin and L. Voronin: Bone mineral and lean tissue loss after long duration space flight. J Musculoskelet Neuronal Interact, 1(2), 157-60 (2000)

2. Y. Arfat, W.-Z. Xiao, S. Iftikhar, F. Zhao, D.-J. Li, Y.-L. Sun, G. Zhang, P. Shang and A.-R. Qian: Physiological Effects of Microgravity on Bone Cells. Calcified Tissue Int, 94(6), 569-579 (2014)
DOI: 10.1007/s00223-014-9851-x

3. L. F. Bonewald: Osteocytes as dynamic multifunctional cells. Ann N Y Acad Sci, 1116, 281-90 (2007)
DOI: 10.1196/annals.1402.018

4. L. F. Bonewald and M. L. Johnson: Osteocytes, mechanosensing and Wnt signaling. Bone, 42(4), 606-615 (2008)
DOI: 10.1016/j.bone.2007.12.224

5. L. I. Plotkin, I. Mathov, J. I. Aguirre, A. M. Parfitt, S. C. Manolagas and T. Bellido: Mechanical stimulation prevents osteocyte apoptosis: requirement of integrins, Src kinases, and ERKs. Am J Physiol Cell Physiol, 289(3), C633-643 (2005)
DOI: 10.1152/ajpcell.00278.2004

6. H. Xu, J. Zhang, J. Wu, Y. Guan, Y. Weng and P. Shang: Oscillatory fluid flow elicits changes in morphology, cytoskeleton and integrin-associated molecules in MLO-Y4 cells, but not in MC3T3-E1 cells. Biol Res, 45(2), 163-9 (2012)
DOI: 10.4067/S0716-97602012000200008

7. J. G. McGarry, J. Klein-Nulend and P. J. Prendergast: The effect of cytoskeletal disruption on pulsatile fluid flow-induced nitric oxide and prostaglandin E2 release in osteocytes and osteoblasts. Biochem Biophys Res Commun, 330(1), 341-8 (2005)
DOI: 10.1016/j.bbrc.2005.02.175

8. M. M. Thi, T. Kojima, S. C. Cowin, S. Weinbaum and D. C. Spray: Fluid shear stress remodels expression and function of junctional proteins in cultured bone cells. Am J Physiol Cell Physiol, 284(2), C389-403 (2003)
DOI: 10.1152/ajpcell.00052.2002

9. S. Burra and J. X. Jiang: Connexin 43 hemichannel opening associated with Prostaglandin E(2) release is adaptively regulated by mechanical stimulation. Commun Integr Biol, 2(3), 239-40 (2009)
DOI: 10.4161/cib.2.3.8154

10. A. E. Loiselle, J. X. Jiang and H. J. Donahue: Gap junction and hemichannel functions in osteocytes. Bone, 54(2), 205-212 (2013)
DOI: 10.1016/j.bone.2012.08.132

11. D. A. Goodenough and D. L. Paul: Beyond the gap: functions of unpaired connexon channels. Nat Rev Mol Cell Biol, 4(4), 285-94 (2003)
DOI: 10.1038/nrm1072

12. S. A. Lloyd, A. E. Loiselle, Y. Zhang and H. J. Donahue: Connexin 43 deficiency desensitizes bone to the effects of mechanical unloading through modulation of both arms of bone remodeling. Bone, 57(1), 76-83 (2013)
DOI: 10.1016/j.bone.2013.07.022

13. S. A. Lloyd, G. S. Lewis, Y. Zhang, E. M. Paul and H. J. Donahue: Connexin 43 deficiency attenuates loss of trabecular bone and prevents suppression of cortical bone formation during unloading. J Bone Miner Res, 2012(19) (1002)
DOI: 10.1002/jbmr.1687

14. S. M. Di, A. R. Qian, L. N. Qu, W. Zhang, Z. Wang, C. Ding, Y. H. Li, H. G. Ren and P. Shang: Graviresponses of osteocytes under altered gravity. Adv Space Res, 48(6), 1161-1166 (2011)
DOI: 10.1016/j.asr.2011.05.030

15. A. Qian, D. Yin, P. Yang, B. Jia, W. Zhang and P. Shang: Development of a ground-based simulated experimental platform for gravitational biology. IEEE Trans Appl Supercond, 19, 42-46 (2009)
DOI: 10.1109/TASC.2009.2012422

16. A. R. Qian, P. F. Yang, L. F. Hu, W. Zhang, S. M. Di, Z. Wang, J. Han, X. Gao and P. Shang: High magnetic gradient environment causes alterations of cytoskeleton and cytoskeleton-associated genes in human osteoblasts cultured in vitro. Adv Space Res, 46(6), 687-700 (2010)
DOI: 10.1016/j.asr.2010.04.012

17. A. R. Qian, X. Gao, W. Zhang, J. B. Li, Y. Wang, S. M. Di, L. F. Hu and P. Shang: Large gradient high magnetic fields affect osteoblast ultrastructure and function by disrupting collagen I or fibronectin/alphabeta1 integrin. PLoS One, 8(1), e51036 (2013)
DOI: 10.1371/journal.pone.0051036

18. A. R. Qian, D. C. Yin, P. F. Yang, Y. Lv, Z. C. Tian and P. Shang: Application of Diamagnetic Levitation Technology in Biological Sciences Research. IEEE Trans Appl Supercond, 23(1), 3600305-3600305 (2013)
DOI: 10.1109/TASC.2012.2232919

19. S. Di, Z. Tian, A. Qian, J. Li, J. Wu, Z. Wang, D. Zhang, D. Yin, M. L. Brandi and P. Shang: Large gradient high magnetic field affects FLG29.1. cells differentiation to form osteoclast-like cells. Int J Radiat Biol, 88(11), 806-13 (2012)
DOI: 10.3109/09553002.2012.698365

20. A. R. Qian, L. Wang, X. Gao, W. Zhang, L. F. Hu, J. Han, J. B. Li, S. M. Di and P. Shang: Diamagnetic levitation causes changes in the morphology, cytoskeleton, and focal adhesion proteins expression in osteocytes. IEEE Trans Biomed Eng, 59(1), 68-77 (2012)
DOI: 10.1109/TBME.2010.2103377

21. A. R. Qian, L. F. Hu, X. Gao, W. Zhang, S. M. Di, Z. C. Tian, P. F. Yang, D. C. Yin, Y. Y. Weng and P. Shang: Large gradient high magnetic field affects the association of MACF1 with actin and microtubule cytoskeleton. Bioelectromagnetics, 30(7), 545-55 (2009)
DOI: 10.1002/bem.20511

22. A. J. Siller-Jackson, S. Burra, S. Gu, X. Xia, L. F. Bonewald, E. Sprague and J. X. Jiang: Adaptation of connexin 43-hemichannel prostaglandin release to mechanical loading. J Biol Chem, 283(39), 26374-82 (2008)
DOI: 10.1074/jbc.M803136200

23. H. Xu, R. Liu, D. Ning, J. Zhang, R. Yang, M. A. Riquelme, J. Li, J. X. Jiang and P. Shang: Biological responses of osteocytic connexin 43 hemichannels to simulated microgravity. J Orthop Res (2016)
DOI: 10.1002/jor.23224

24. N. Batra, S. Burra, A. J. Siller-Jackson, S. Gu, X. Xia, G. F. Weber, D. DeSimone, L. F. Bonewald, E. M. Lafer, E. Sprague, M. A. Schwartz and J. X. Jiang: Mechanical stress-activated integrin alpha5beta1 induces opening of connexin 43 hemichannels. Proc Natl Acad Sci U S A, 109(9), 3359-64 (2012)
DOI: 10.1073/pnas.1115967109

25. R. Kar, M. A. Riquelme, S. Werner and J. X. Jiang: Connexin 43 channels protect osteocytes against oxidative stress-induced cell death. J Bone Miner Res, 28(7), 1611-21 (2013)
DOI: 10.1002/jbmr.1917

26. J. A. Orellana, K. F. Shoji, V. Abudara, P. Ezan, E. Amigou, P. J. Saez, J. X. Jiang, C. C. Naus, J. C. Saez and C. Giaume: Amyloid beta-induced death in neurons involves glial and neuronal hemichannels. J Neurosci, 31(13), 4962-77 (2011)
DOI: 10.1523/JNEUROSCI.6417-10.2011

27. X. Fang, T. Huang, Y. Zhu, Q. Yan, Y. Chi, J. X. Jiang, P. Wang, H. Matsue, M. Kitamura and J. Yao: Connexin43 hemichannels contribute to cadmium-induced oxidative stress and cell injury. Antioxid Redox Signal, 14(12), 2427-39 (2011)
DOI: 10.1089/ars.2010.3150

28. L. F. Bonewald: Establishment and characterization of an osteocyte-like cell line, MLO-Y4. J Bone Miner Metab, 17(1), 61-5 (1999)
DOI: 10.1007/s007740050066

29. A. Qian, S. Di, X. Gao, W. Zhang, Z. Tian, J. Li, L. Hu, P. Yang, D. Yin and P. Shang: cDNA microarray reveals the alterations of cytoskeleton-related genes in osteoblast under high magneto-gravitational environment. Acta Biochim Biophys Sin (Shanghai), 41(7), 561-77 (2009)
DOI: 10.1093/abbs/gmp041

30. B. E. Hammer, L. S. Kidder, P. C. Williams and W. W. Xu: Magnetic Levitation of MC3T3 Osteoblast Cells as a Ground-Based Simulation of Microgravity. Microgravity Sci Technol, 21(4), 311-318 (2009)
DOI: 10.1007/s12217-008-9092-6

31. Y. Wang, Z.-H. Chen, C. Yin, J.-H. Ma, D.-J. Li, F. Zhao, Y.-L. Sun, L.-F. Hu, P. Shang and A.-R. Qian: GeneChip expression profiling reveals the alterations of energy metabolism related genes in osteocytes under large gradient high magnetic fields. PloS One, 10(1), e0116359 (2015)
DOI: 10.1371/journal.pone.0116359

32. M. A. Riquelme, R. Kar, S. Gu and J. X. Jiang: Antibodies targeting extracellular domain of connexins for studies of hemichannels. Neuropharmacology, 75, 525-32 (2013)
DOI: 10.1016/j.neuropharm.2013.02.021

33. P. P. Cherian, A. J. Siller-Jackson, S. Gu, X. Wang, L. F. Bonewald, E. Sprague and J. X. Jiang: Mechanical strain opens connexin 43 hemichannels in osteocytes: a novel mechanism for the release of prostaglandin. Mol Biol Cell, 16(7), 3100-6 (2005)
DOI: 10.1091/mbc.E04-10-0912

34. D. C. Genetos, C. J. Kephart, Y. Zhang, C. E. Yellowley and H. J. Donahue: Oscillating fluid flow activation of gap junction hemichannels induces atp release from MLO-Y4 osteocytes. J Cell Physiol, 212(1), 207-214 (2007)
DOI: 10.1002/jcp.21021

35. A. L. Rath, L. F. Bonewald, J. Ling, J. X. Jiang, M. E. Van Dyke and D. P. Nicolella: Correlation of cell strain in single osteocytes with intracellular calcium, but not intracellular nitric oxide, in response to fluid flow. J Biomech, 43(8), 1560-4 (2010)
DOI: 10.1016/j.jbiomech.2010.01.030

36. J. I. Aguirre, L. I. Plotkin, S. A. Stewart, R. S. Weinstein, A. M. Parfitt, S. C. Manolagas and T. Bellido: Osteocyte Apoptosis Is Induced by Weightlessness in Mice and Precedes Osteoclast Recruitment and Bone Loss. J Bone Miner Res, 21(4), 605-615 (2006)
DOI: 10.1359/jbmr.060107

37. N. Basso and J. N. M. Heersche: Effects of hind limb unloading and reloading on nitric oxide synthase expression and apoptosis of osteocytes and chondrocytes. Bone, 39(4), 807-814 (2006)
DOI: 10.1016/j.bone.2006.04.014

38. S. D. Tan, A. D. Bakker, C. M. Semeins, A. M. Kuijpers-Jagtman and J. Klein-Nulend: Inhibition of osteocyte apoptosis by fluid flow is mediated by nitric oxide. Biochem Bioph Res Co, 369(4), 1150-1154 (2008)
DOI: 10.1016/j.bbrc.2008.03.007

39. H. Xu, S. Gu, M. A. Riquelme, S. Burra, D. Callaway, H. Cheng, T. Guda, J. Schmitz, R. J. Fajardo, S. L. Werner, H. Zhao, P. Shang, M. L. Johnson, L. F. Bonewald and J. X. Jiang: Connexin 43 Channels are Essential for Normal Bone Structure and Osteocyte Viability. J Bone Miner Res, 30(3), 550-62 (2015)
DOI: 10.1002/jbmr.2374

40. N. Batra, R. Kar and J. X. Jiang: Gap junctions and hemichannels in signal transmission, function and development of bone. Biochimica et Biophysica Acta, 1818(8), 1909-18 (2012)
DOI: 10.1016/j.bbamem.2011.09.018

41. W.-Y. Cheung, C. A. Simmons and L. You: Osteocyte apoptosis regulates osteoclast precursor adhesion via osteocytic IL-6 secretion and endothelial ICAM-1 expression. Bone, 50(1), 104-110 (2012)
DOI: 10.1016/j.bone.2011.09.052

42. Y. Fujikawa, A. Sabokbar, S. D. Neale, I. Itonaga, T. Torisu and N. A. Athanasou: The effect of macrophage-colony stimulating factor and other humoral factors (interleukin-1, -3, -6, and -11, tumor necrosis factor-α, and granulocyte macrophage-colony stimulating factor) on human osteoclast formation from circulating cells. Bone, 28(3), 261-267 (2001)
DOI: 10.1016/S8756-3282(00)00453-1

Abbreviations: HMGE: high magneto-gravitational environment; Cx43: connexin 43; T: Tesla; DEGs: differentially expressed genes; GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; MAPK: mitogen-activated protein kinase; TGF-beta: transforming growth factor beta; VEGF: vascular endothelial growth factor; Jak-STAT: Janus kinase-signal transducer and activator of transcription; PGs: prostaglandins

Key Words: Gene Expression Profile, HMGE, Connexin 43, Hemichannel, MLO-Y4 cells

Send correspondence to: Peng Shang,127 Youyi Xilu, Xi’an, Shaanxi Province, 710072, China, Tel: 86-29-88460391, Fax: 86-29-88461671, E-mail: celldon@126.com