The polymorphisms of MIR31HG gene is correlated with alcohol-induced osteonecrosis of the femoral head in Chinese Han male population

Introduction

Osteonecrosis of the femoral head (ONFH) is a common hip illness characterized by impaired microvascular circulation leading to the death of bone cells. It eventually causes structural changes, collapse of the femoral head, and joint dysfunction. It is reported that the age of patients with ONFH is mostly 30~50 years old, and the incidence rate is gradually increasing with the increase of age (1). There are 8.12 million ONFH patients in China, with approximately 150,000 to 200,000 newly diagnosed ONFH patients yearly (2). Some research showed that one of the common causes of ONFH was excessive drinking (3, 4). Besides, Cui et al. reported that alcoholic ONFH accounted for 30.7% of all femoral head necrosis in China (5). Nevertheless, it has been found in clinical practice that only some people who drink excessively suffer from femoral head necrosis, which suggested that genetic susceptibility may contribute to the occurrence of alcoholic ONFH. Meanwhile, a large number of reports confirmed that genetic variants were related to alcoholic ONFH predisposition (68).

Long non-coding RNAs (lncRNAs) refers to non-coding RNA with a length of over 200 nucleotides, which can regulate the physiological and pathological activities of organisms by participating in gene transcription and post-transcriptional regulation, epigenetic modification, and translation (9, 10). Recently, lncRNAs have aroused researcher’s concern in bone development, differentiation, and osteonecrosis. For example, Tang et al. found that lncRNA-OG could promote the differentiation of osteoblasts (11). Liu et al. indicated that lncRNA AK077216 facilitated bone resorption and osteoclastogenesis (12). Besides, a study performed by Xiang et al, which reported that lncRNA RP11-154D6 was reduced in ONFH patients and involved in the progression of ONFH (13).

MIR31HG (also named as lncHIFCAR, LOC554202) is described to participant in proliferation, differentiation, invasion of cancer cells (14, 15). In addition, researchers have discovered that MIR31HG silence promoted osteogenic differentiation and relieved the inflammation-induced inhibition of osteogenesis (16). However, the role of MIR31HG in alcoholic ONFH development remains elusive. However, the role of MIR31HG in alcoholic ONFH development remains elusive.

Here, we explored the association of MIR31HG genetic variations with alcoholic ONFH susceptibility among the Chinese Han male population. This will help to provide new understandings for MIR31HG into the pathogenesis of alcoholic ONFH.

MethodsStudy subjects

A total of 733 Chinese Han males (308 alcoholic ONFH patients and 425 healthy controls) were enrolled. All patients were randomly selected from Hainan Affiliated Hospital of Hainan Medical University. The patients met the criteria as following: 1) exceeded 400 mL/week (17) alcohol intake for more than 6 months; 2) diagnosed ONFH within one year after drinking alcohol; 3) no rheumatoid arthritis, hyperlipidemia, spinal cord cavitation, decompression sickness, osteoporosis, cardiovascular disease, and virus infection; 4) no history of steroid use; 5) The diagnosis of alcoholic ONFH by X-ray, nuclear magnetic resonance imaging (MRI), and computed tomography (CT). The stage of alcoholic ONFH patients were identified by the Ficat Classification system (18). Healthy people were selected based on 1) drinking habits or greater than 400 mL alcohol intake per week for more than 6 months; 2) no history of the traumatic disease (ONFH, rheumatoid arthritis, hyperlipidemia, spinal cord cavitation, decompression sickness, osteoporosis, cardiovascular disease, etc…), and no steroid use. The study protocol was approved by the Ethics Committee of the Hainan Affiliated Hospital of Hainan Medical University, in compliance with the declaration of Helsinki. And the signature consent of participants was received.

SNP genotyping

Seven SNPs in MIR31HG were chosen with a minor allele frequency (MAF) >5% in 1000 Genomes Chinese Han Beijing population. Total DNA was isolated from peripheral blood cells by GoldMag blood DNA Kit (GoldMag Co. Ltd, Xi’an, China), and the concentration were evaluated via NanoDrop 2000 (Thermo Scientific, USA). SNP genotyping was completed through the Agena MassARRAY platform and the data were analyzed using the Agena Typer 4.0 software.

Data analysis

Student t-test was applied to evaluate age difference between the two groups. Hardy-Weinberg equilibrium (HWE) in controls was calculated using χ2 test. The linkage between MIR31HG polymorphisms and alcoholic ONFH susceptibility was examined by odds ratio (OR) and 95% confidence interval (CI) through logistic regression. Linkage disequilibrium (LD) was analyzed by Haploview software. Multifactor dimensionality reduction (MDR) was conducted to evaluate the SNP-SNP interactions in the risk of ONFH. False positive report probability (FPRP) values and statistical power were calculated (19). P < 0.05 was identified a significant difference.

ResultsParticipants and selected SNPs in MIR31HG

As listed in Table 1, this study included 308 alcoholic ONFH patients and 425 healthy controls, with an average age of 43.37 ± 11.34 years and 42.73 ± 12.88 years, respectively. The age distribution between the two groups was well-matched (p=0.478).

www.frontiersin.org

Table 1 Characteristics of ONFH patients and controls in this study.

We chose and genotyped seven SNPs (rs1332184, rs72703442, rs2025327, rs55683539, rs2181559, rs10965059, rs10965064) in the intron region of MIR31HG, and all SNPs in line with HWE (p>0.05). And we found that rs10965059 could reduce the susceptibility of alcoholic ONFH in the allele model (p < 0.001, OR = 0.48, 95%CI=0.35-0.66, Table 2).

www.frontiersin.org

Table 2 Basic information of the selected SNPs in MIR31HG.

Alcoholic ONFH risk assessment

The association of seven SNPs in MIR31HG and alcoholic ONFH risk was evaluated (Table 3). The results revealed that MIR31HG- rs10965059 was related to a lower risk of alcoholic ONFH under codominant (p<0.001, OR=0.43, 95%CI=0.30-0.62), dominant (p<0.001, OR=0.42, 95%CI=0.30-0.61) and additive models (p<0.001, OR=0.45, 95%CI=0.32-0.62).

www.frontiersin.org

Table 3 Association of MIR31HG polymorphisms with alcohol-induced ONFH.

Age stratification showed that rs10965059 decreased the susceptibility to alcoholic ONFH individuals older than 40 years in the allele (p<0.001, OR=0.36), codominant (p<0.001, OR=0.45), dominant (p<0.001, OR=0.31), and additive models (p<0.001, OR=0.33, 9 Table 4). Rs10965064 only reduced the alcoholic ONFH susceptibility under the dominant model (p=0.049, OR=0.67).

www.frontiersin.org

Table 4 Relationships between MIR31HG SNPs and alcohol-induced ONFH susceptibility based on stratification by age.

Besides, the necrotic sites stratification results (Table 5) indicated that rs10965059 played a protective role in alcoholic bilateral ONFH in the allele (OR=0.45, p<0.001), codominant (OR=0.40, p<0.001), dominant (OR=0.39, p<0.001), and recessive model (OR=0.42, p<0.001).

www.frontiersin.org

Table 5 Association between MIR31HG polymorphisms and ONFH risk stratified by necrotic sites.

Haplotype analysis and MDR analysis

We analyzed the haplotype of MIR31HG gene. Table 6 represented that there was no linkage between haplotypes and alcoholic ONFH susceptibility (p>0.05). Besides, we found a LD block formed by rs72703442, rs2025327, and rs55683539 (Figure 1).

www.frontiersin.org

Table 6 Haplotype analysis of MIR31HG SNPs with alcohol-induced ONFH.

www.frontiersin.org

Figure 1 Linkage disequilibrium (LD) plots containing three polymorphisms from MIR31HG. Block 1 includes rs72703442, rs2025327, and rs55683539. The numbers inside the diamonds indicate the D’ for pairwise analyses.

The Fruchterman-Rheingold of SNP-SNP interaction was presented in Figure 2. The results of Supplemental Table 1 revealed that rs10965059 was the best single locus model to prediction the ONFH susceptibility (testing accuracy=0.581, CVC=10/10, p < 0.0001).

www.frontiersin.org

Figure 2 The Fruchterman-Reingold of SNP-SNP interactions. Each SNP is reported in per cent the value of Information Gain (IG), while numbers in the connections indicate the entropy-based IG for the SNP pairs. Orange bar indicate the high-level synergies on the phenotype, while the brown indicate a medium-level interaction, green and blue connections with negative IG values indicate redundancy or lack of synergistic interactions between the markers.

FPRP analysis

As shown in Supplemental Table 2, all significant results of rs10965059 remained noteworthy, at the prior probability of 0.001 and FPRP threshold of 0.2.

Discussion

We illustrated the relationship of MIR31HG polymorphisms with alcoholic ONFH susceptibility in this study. Results of our research indicated that MIR31HG-rs10965059 decreased the susceptibility of alcoholic ONFH overall. In the age stratification, we also observed that rs10965059 and rs10965064 had protective effect on alcoholic ONFH occurrence older than 40 years old. Besides, rs10965059 could reduce the alcoholic bilateral ONFH risk. These data underline the importance of MIR31HG in alcoholic ONFH occurrence and may serve as a new biomarker for the early prevention and treatment of alcoholic ONFH.

MIR31HG is located at the chromosomal locus 9q21.3 in humans. It has been described to participant in cell proliferation, differentiation in many diseases (2022). In recent years, the functional role of MIR31HG in bone-related diseases has been studied. For example, Sun et al. found that MIR31HG was higher in tissues of osteosarcoma patients compared with healthy controls and it could regulate osteosarcoma cell growth and migration via miR-361 (23). Ma et al. have shown that MIR31HG was elevated in chordoma patients and MIR31HG silence repressed the migration, growth, and invasion of chordoma cells (24). Besides, suppression of MIR31HG could facilitate the differentiation of osteoblast in human adipose-derived stem cell (16). Moreover, a genome-wide association studies uncovered that MIR31HG polymorphism was related to radius bone density and content in children (25). These evidences led us to hypothesize that MIR31HG may have pathogenic significance in alcoholic ONFH. Here, we firstly observed the contribution of MIR31HG polymorphisms to alcoholic ONFH susceptibility. Rs10965059 exerted a protective role in alcoholic ONFH occurrence in both overall and stratified analysis.

Rs10965059 polymorphism is in the intron region of MIR31HG. Zhao et al. confirmed that the common intronic WDFY4 rs877819 affects the expression of WDFY4 gene by affecting YY1 binding (26). Choi et al. also found that intronic SNP (rs2280964) significantly correlated with reduced the expression of CXCR3 gene, which resulted in changes of immune cell responses to chemokine-cytokine signaling in ex vivo and vitro (27). Based on the above studies, we speculatedthat rs10965059 may affect the susceptibility of alcoholic ONFH by altering the expression of MIR31HG gene. In follow-up studies, we will explore the functional consequence of intronic SNP rs10965059 in vitro and ex vivo to support our hypothesis.

Inevitably, there are some limitations to this study. First, we did not conduct a functional analysis, which is essential to understand the role of MIR31HG in alcoholic ONFH. Second, the subjects were all Chinese Han population, and there may be a certain selection bias. Therefore, we needed animal or cell experiments and more ethnic groups to verify our findings.

Conclusions

To sum up, we confirmed that MIR31HG variants have a significant correlation with the occurrence of alcoholic ONFH in a Chinese Han male population. This may provide new ideas for the prevention and diagnosis of alcoholic ONFH.

Data availability statement

The data presented in the study are deposited in the Zenodo repository, accession number DOI: 10.5281/zenodo.7349781.

Ethics statement

The studies involving human participants were reviewed and approved by Hainan Affiliated Hospital of Hainan Medical University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

JX designed this study protocol; WL drafted the manuscript; XW performed the DNA extraction and genotyping; JC performed the data analysis; FZ performed the sample collection and information recording. JX conceived and supervised the study. All authors read and approved the final manuscript.

Funding

This study was supported by Hainan Provincial Natural Science Foundation of China (No. 822RC808) and Project supported by the Education Department of Hainan Province (No. Hnky2022ZD-13).

Conflict of interest

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

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.976165/full#supplementary-material

Abbreviations

ONFH, Osteonecrosis of the femoral head; OR, odds ratio; CI, confidence intervals; MAF, minor allele frequency; HWE, Hardy-Weinberg Equilibrium; CT, computed tomography; MRI, nuclear magnetic resonance imaging; LD, linkage disequilibrium.

References

1. Zalavras CG, Lieberman JR. Osteonecrosis of the femoral head: evaluation and treatment. J Am Acad Orthopaedic Surgeons (2014) 22(7):455–64. doi: 10.5435/JAAOS-22-07-455

CrossRef Full Text | Google Scholar

2. Song Y, Du Z, Ren M, Yang Q, Wang Q, Chen G, et al. Association of gene variants of transcription factors PPARγ, RUNX2, osterix genes and COL2A1, IGFBP3 genes with the development of osteonecrosis of the femoral head in Chinese population. Bone (2017) 101:104–12. doi: 10.1016/j.bone.2017.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Hong G, Han X, He W, Xu J, Sun P, Shen Y, et al. Analysis of circulating microRNAs aberrantly expressed in alcohol-induced osteonecrosis of femoral head. Sci Rep (2019) 9(1):18926. doi: 10.1038/s41598-019-55188-6

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Yoon BH, Kim TY, Shin IS, Lee HY, Lee YJ, Koo KH. Alcohol intake and the risk of osteonecrosis of the femoral head in Japanese populations: a dose-response meta-analysis of case-control studies. Clin Rheumatol (2017) 36(11):2517–24. doi: 10.1007/s10067-017-3740-4

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cui L, Zhuang Q, Lin J, Jin J, Zhang K, Cao L, et al. Multicentric epidemiologic study on six thousand three hundred and ninety five cases of femoral head osteonecrosis in China. Int Orthopaedics (2016) 40(2):267–76. doi: 10.1007/s00264-015-3061-7

CrossRef Full Text | Google Scholar

6. Guo Y, Cao Y, Gong S, Zhang S, Hou F, Zhang X, et al. Correlation analysis between CARMEN variants and alcohol-induced osteonecrosis of the femoral head in the Chinese population. BMC Musculoskeletal Disord (2020) 21(1):547. doi: 10.1186/s12891-020-03553-2

CrossRef Full Text | Google Scholar

7. Wang J, Shi X, Yang H, Du J, Ouyang Y, Wang H, et al. Association between alcohol-induced osteonecrosis of femoral head and risk variants of MMPS in han population based on a case-control study. Oncotarget (2017) 8(38):64490–8. doi: 10.18632/oncotarget.16380

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Liu C, An F, Cao Y, Wang J, Tian Y, Wu H, et al. Significant association between RETN genetic polymorphisms and alcohol-induced osteonecrosis of femoral head. Mol Genet Genomic Med (2019) 7(8):e822. doi: 10.1002/mgg3.822

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Yang L, Lin C, Liu W, Zhang J, Ohgi KA, Grinstein JD, et al. ncRNA- and Pc2 methylation-dependent gene relocation between nuclear structures mediates gene activation programs. Cell (2011) 147(4):773–88. doi: 10.1016/j.cell.2011.08.054

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Zhou J, Yang L, Zhong T, Mueller M, Men Y, Zhang N, et al. H19 lncRNA alters DNA methylation genome wide by regulating s-adenosylhomocysteine hydrolase. Nat Commun (2015) 6:10221. doi: 10.1038/ncomms10221

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tang S, Xie Z, Wang P, Li J, Wang S, Liu W, et al. LncRNA-OG promotes the osteogenic differentiation of bone marrow-derived mesenchymal stem cells under the regulation of hnRNPK. Stem Cells (Dayton Ohio). (2019) 37(2):270–83. doi: 10.1002/stem.2937

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Liu C, Cao Z, Bai Y, Dou C, Gong X, Liang M, et al. LncRNA AK077216 promotes RANKL-induced osteoclastogenesis and bone resorption via NFATc1 by inhibition of NIP45. J Cell Physiol (2019) 234(2):1606–17. doi: 10.1002/jcp.27031

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Xiang S, Li Z, Weng X. The role of lncRNA RP11-154D6 in steroid-induced osteonecrosis of the femoral head through BMSC regulation. J Cell Biochem (2019) 120(10):18435–45. doi: 10.1002/jcb.29161

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wang R, Ma Z, Feng L, Yang Y, Tan C, Shi Q, et al. LncRNA MIR31HG targets HIF1A and P21 to facilitate head and neck cancer cell proliferation and tumorigenesis by promoting cell-cycle progression. Mol cancer. (2018) 17(1):162. doi: 10.1186/s12943-018-0916-8

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Nie FQ, Ma S, Xie M, Liu YW, De W, Liu XH. Decreased long noncoding RNA MIR31HG is correlated with poor prognosis and contributes to cell proliferation in gastric cancer. Tumour Biol (2016) 37(6):7693–701. doi: 10.1007/s13277-015-4644-z

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Jin C, Jia L, Huang Y, Zheng Y, Du N, Liu Y, et al. Inhibition of lncRNA MIR31HG promotes osteogenic differentiation of human adipose-derived stem cells. Stem Cells (Dayton Ohio) (2016) 34(11):2707–20. doi: 10.1002/stem.2439

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Matsuo K, Hirohata T, Sugioka Y, Ikeda M, Fukuda A. Influence of alcohol intake, cigarette smoking, and occupational status on idiopathic osteonecrosis of the femoral head. Clin orthopaedics related Res (1988) 234):115–23. doi: 10.1097/00003086-198809000-00021

CrossRef Full Text | Google Scholar

18. Ficat RP. Idiopathic bone necrosis of the femoral head. early diagnosis and treatment. J Bone Joint Surg Br volume. (1985) 67(1):3–9. doi: 10.1302/0301-620X.67B1.3155745

CrossRef Full Text | Google Scholar

19. Wacholder S, Chanock S, Garcia-Closas M, El Ghormli L, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst (2004) 96(6):434–42. doi: 10.1093/jnci/djh075

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Gao J, Chen F, Hua M, Guo J, Nong Y, Tang Q, et al. Knockdown of lncRNA MIR31HG inhibits cell proliferation in human HaCaT keratinocytes. Biol Res (2018) 51(1):30. doi: 10.1186/s40659-018-0181-8

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Yan S, Tang Z, Chen K, Liu Y, Yu G, Chen Q, et al. Long noncoding RNA MIR31HG inhibits hepatocellular carcinoma proliferation and metastasis by sponging microRNA-575 to modulate ST7L expression. J Exp Clin Cancer Res: CR (2018) 37(1):214. doi: 10.1186/s13046-018-0853-9

CrossRef Full Text | Google Scholar

22. He A, Chen Z, Mei H, Liu Y. Decreased expression of LncRNA MIR31HG in human bladder cancer. Cancer Biomarkers: Section Dis Markers (2016) 17(2):231–6. doi: 10.3233/CBM-160635

CrossRef Full Text | Google Scholar

23. Sun Y, Jia X, Wang M, Deng Y. Long noncoding RNA MIR31HG abrogates the availability of tumor suppressor microRNA-361 for the growth of osteosarcoma. Cancer Manage Res (2019) 11:8055–64. doi: 10.2147/CMAR.S214569

CrossRef Full Text | Google Scholar

24. Ma X, Qi S, Duan Z, Liao H, Yang B, Wang W, et al. Long non-coding RNA LOC554202 modulates chordoma cell proliferation and invasion by recruiting EZH2 and regulating miR-31 expression. Cell Proliferation (2017) 50(6):e12388. doi: 10.1111/cpr.12388

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Chesi A, Mitchell JA, Kalkwarf HJ, Bradfield JP, Lappe JM, McCormack SE, et al. A trans-ethnic genome-wide association study identifies gender-specific loci influencing pediatric aBMD and BMC at the distal radius. Hum Mol Genet (2015) 24(17):5053–9. doi: 10.1093/hmg/ddv210

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhao H, Yang W, Qiu R, Li J, Xin Q, Wang X, et al. An intronic variant associated with systemic lupus erythematosus changes the binding affinity of Yinyang1 to downregulate WDFY4. Genes Immunity (2012) 13(7):536–42. doi: 10.1038/gene.2012.33

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Choi JW, Park CS, Hwang M, Nam HY, Chang HS, Park SG, et al. A common intronic variant of CXCR3 is functionally associated with gene expression levels and the polymorphic immune cell responses to stimuli. J Allergy Clin Immunol (2008) 122(6):1119–26.e7. doi: 10.1016/j.jaci.2008.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

留言 (0)

沒有登入
gif