key: cord-263322-y4htkvux
authors: Yang, Jun; Jing, Ming; Yang, Xiaoge
title: Association between genetic polymorphisms and osteonecrosis in steroid treatment populations: a detailed stratified and dose-response meta-analysis
date: 2019-05-14
journal: Biosci Rep
DOI: 10.1042/bsr20190024
sha: 
doc_id: 263322
cord_uid: y4htkvux

Steroid treatment has become recognized as an important risk factor for avascular osteonecrosis of the femoral head. However, not all patients who receive long-term, high-dose steroids develop osteonecrosis, indicating that there are individual differences in occurrence. We explored the relationship between polymorphisms and steroid-induced osteonecrosis of the femoral head (SONFH) incidence with variables. We used a multilevel mixed-effects logistic regression model, which is an expansion of logistic regression, for each type of steroid, primary disease, drug dose, applied duration, and single-nucleotide polymorphism (SNP). We also conducted a dose-response meta-analysis to analyze the cumulative dosage and SONFH risk in mutation carriers. There were significant correlations between the ABCB1 rs1045642 mutant and SONFH in the prednisone-use and methylprednisolone/prednisone-use populations. The ABCB1 rs2032582 mutant homozygote had a protective effect in the methylprednisolone/prednisolone renal transplant population. For ApoB rs693, mutation increased the incidence of SONFH in prednisone-use and methylprednisolone/prednisolone-use populations and renal transplant patients. For ApoB rs1042031, mutation increased the risk of SONFH in the prednisone-use population. The PAI-1 rs1799768 mutation had a protective effect on the SONFH risk prednisone-use and renal transplant populations. ABCB1 rs1045642 mutations have a protective effect against SONFH, and ApoB rs693 and rs1042031 increase the SONFH risk. Cumulative dosage and treatment duration had little effect on the results. In addition, there was a dose-effect correlation in ABCB1 rs1045642 and rs2032582 mutation carriers.

Steroid-induced osteonecrosis of the femoral head (SONFH), which leads to collapse of the femoral head and articular dysfunction, has an incidence of 9-40% amongst patients receiving steroid treatment [1] . The exact pathology of SONFH is still unclear and might be related to lipid metabolism disorders, abnormal microcirculation, insufficient blood supply, inflammation, and bone marrow mesenchymal stem cell osteogenesis differentiation dysfunction. The abnormal blood supply leads to the apoptosis of osteocytes and osteoblasts, followed by bone loss and reduced bone mineral density [2] . Lipid metabolism dysfunction is also an important pathology, and steroid application may lead to an increase in blood lipid levels and then to the development of intravascular lipid embolism in the microvasculature [3] . Embolism accumulation might affect microcirculation and increase the pressure of the intramedullary cavity, eventually leading to the death of bone cells. Long-term or mass steroid application is the critical pathogenesis of SONFH [4] .

However, not all patients who receive long-term, high-dose steroids will develop SONFH, indicating that there are individual differences in the occurrence of SONFH.

To date, several meta-analyses have been published and shown that PAI-1 4G/5G (rs1799768) [5] , ABCB1 C3435T (rs1045642) [5] [6] [7] and CYP3A activity [8] are associated with SONFH incidence. However, the results of the ABCB1 G2677T/A polymorphism are still disputed [5] [6] [7] . Our previous research reported only the single-nucleotide polymorphisms (SNPs) that appeared in more than three studies and indicated that ABCB1 rs1045642 has a protective effect on SONFH in an allelic model and that the ApoB rs693 and rs1042031 mutations promote the pathogenesis of SONFH. ABCB1 rs2032582, MTHFR rs1801133, and PAI-1 rs1799768 were not correlated with SONFH incidence. However, heterogeneity still exists in the previous results, and further analysis of the characteristics of the included studies is needed. Therefore, this research further analyzes the effects of primary disease, type of steroids, cumulative steroid dosage, and treatment duration on the results.

We used the meta-analysis of observational studies in epidemiology guidelines in the present study [9] .

Two authors independently performed a literature search of the PubMed, Embase, Cochrane Library, and Chinese public databases, including the China National Knowledge Infrastructure, the China Biology Medicine Database, the China Science Periodical Database (Wanfang Database), and the VIP Journal Integration Platform. The search included studies published through 29 July 2018. The following terms were used in the search strategy: hormone, glucocorticoid, steroid, corticosteroid, osteonecrosis, femoral, femur, femoris, whirlbone, polymorphism, SNP, genetic, mutation, genotype, allele, allelic, and variation. The search strategy is shown in Supplementary Table S1. Divergence in the search results was resolved by discussion.

The studies were included in our meta-analysis if they met the following criteria: (1) case-control or cohort studies comparing a population that suffered SONFH with a population that did not suffer after steroid treatment, (2) studies assessing the associations between genetic polymorphisms and SONFH, and (3) studies reporting the frequencies of specific alleles or the effect sizes of individual genotypes between cases and controls. Studies were excluded from the analysis for the following reasons: (1) noncase-control or noncohort studies, (2) the case group included SONFH patients with other etiologies or SONFH patients who were not reported separately, (3) the control group included an ONFH population without steroid application or a healthy population, (4) non-SNP-related studies, (5) studies about family heredity, and (6) studies that did not report data pertaining to allelic frequencies or effect size. In addition, conference reports, editor comments, reviews, case reports, and academic dissertations were excluded from the analysis.

Two authors independently extracted the following data from each eligible study: first author's name, publication year, research location, sample size, average subject age, primary disease, steroid type, cumulative steroid dose, treatment duration, and genes of interest. If the average or median value of the cumulative steroid dosage in each group was not reported, it was estimated. The cumulative dose was generally calculated by multiplying the daily average dose and the treatment duration. For studies with a fixed treatment strategy, the cumulative dosage was calculated by multiplying the dosage of each treatment cycle by the total number of cycles. In the present study, we defined the average body surface area as 1.6 m 2 and the average body weight as 60 kg. For closed-interval categories, the value assigned to each dosage or duration, such as the average or median, was the midpoint. For open-interval categories, the value was estimated as the limit value multiplied by 1.2 [10] . Because of the different types of steroids used in the included studies, it was necessary to convert the doses of different types of steroids to the equivalent dose of prednisone. The conversion standard is 5 mg prednisone = 5 mg prednisolone = 4 mg methylprednisolone = 0.75 mg dexamethasone [11] . The Newcastle-Ottawa Scale (NOS), a validated tool for evaluating the quality of observation studies, was used to evaluate the methodological quality of the included studies; the scale includes the following three subscales: selection, comparability, and exposure [12] . 

Association analysis was performed using five genetic models: allelic (W vs M), dominant (WW+WM vs MM), recessive (WW vs WM+MM), heterozygous (WM vs MM), and homozygous (WW vs MM) models [13] . W is the wild allele and M is the mutation allele. The odds ratios (ORs) and their 95% CIs were used to assess the strength of the associations between polymorphisms and SONFH incidence. The I 2 statistic was used to estimate the degree of heterogeneity amongst the studies. Cochran's Q statistics were also taken as measures of between-study heterogeneity. Generally, a significance level of P<0.1 suggests the existence of significant heterogeneity to expand a test's sensitivity [13] . If the I 2 ≥ 50% (Q test, P<0.1), the random-effect model was used; otherwise, the fixed-effect model was used. Then, we attempted to use a multilevel mixed-effects logistic regression model, which is an expansion of logistic regression, for each type of steroid, primary disease, drug dose, applied duration, and SNP. The factors of different genotypes and other variables were included as fixed effects, and different studies were considered random effects.

Finally, we conducted a dose-response meta-analysis across the cumulative steroid dosage [14] . To derive the dose-response curve, we modeled the dose using restricted cubic splines with three knots at fixed percentiles of 10, 50, and 90% of the distribution [15] . All tests were two-tailed, and a p-value less than 0.05 was considered statistically significant except in Cochran's Q test.

Our research returned 287 English articles and 305 Chinese articles after removing duplicates. After screening the titles and abstracts, 509 of these articles were excluded. The full texts of 83 articles were assessed, amongst which studies were excluded for the following reasons: the control group did not receive steroid therapy (26); reviews (8);

non-SNP-related studies (5); basic studies (3); studies without frequencies or effect size results (2); non-steroid-related ONFH studies (2); case reports (2); microRNA-related studies (2); duplicated reports (2); and heredity SONFH study (1) . Ultimately, we included 30 studies that assessed a total of 7553 patients in our meta-analysis (Figure 1 , Table 1 ).

All studies were case-control designs, except one study, which was a cohort design [18] . The primary diseases included organ transplantation, hematological disease, systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), nephrotic syndrome, and ophthalmic disease. In addition, two studies by Wei He [30, 31] , two studies by Weibao Fang [27, 28] , and studies by Masaaki and Kyoko [32, 35] included the same patient groups but did not assess the same SNPs. The NOS scores ranged from six to nine points, and the overall quality of the observational studies was ideal ( Table  1 ). The details of NOS scores for each study are shown in Supplementary Table S2 .

Amongst the preliminary results of SNPs that were included in more than three studies, an allelic model showed that the ABCB1 rs1045642 mutation had a protective effect against SONFH (OR: 0.74; 95% CI: 0. Figure 2 ). Other SNPs were not included in the analysis because they did not appear in more than three studies. A further analysis was performed in the present study using a multilevel mixed-effects logistic regression model. The details of the included studies for each group analysis are listed in Supplementary Table S3 .

In the results of the ABCB1 rs1045642 multilevel mixed-effects regression model, the mutant homozygote and the mutant carrier could significantly protect against SONFH after steroid application (TT vs CC genotype: OR: 0.51; 95% CI: 0. 34 

primary disease showed no significant results. It should be noted that the studies that did not report primary disease type or that included several types of primary disease were not included in the categorical analysis of primary disease ( Figure 3) .

For ABCB1 rs2032582, the negative results did not change after considering cumulative dose and treatment duration. In the categorical analysis, the protective effect of the mutant was only found in the methylprednisolone/prednisolone-use renal transplant population (OR: 0.27; 95% CI: 0.07-0.97; P=0.046). However, this result was based on a single study and had a large standard error (SE = 0.66); thus, this result needs to be confirmed ( Figure 4) .

For ApoB rs693, mutation increased the incidence of SONFH (CT Figure  5 ). Due to a large number of zero events, the regression results did not converge after considering the cumulative dose. The results were not changed after considering treatment duration. In addition, in prednisone-use and methylprednisolone/prednisolone-use populations and renal transplant patients, mutation increased the risk of SONFH.

For ApoB rs1042031, mutation increased the risk of SONFH (GA vs GG: OR: 2.68; 95% CI: 1.43-5.01; P=0.002) ( Figure 6 ). The results were not obviously changed after considering cumulative dose and treatment duration. In the categorical analysis, mutation increased the risk of SONFH in the prednisone-use population (OR: 5 (Figure 8 ). These results differed from previous results because in previous research a random-effect model was used due to large heterogeneity, and a mixed-effect model was adopted in this research. The results did not change after considering cumulative dose and treatment duration. In the categorical analysis, the protective effect existed in the prednisone-use and renal transplant population. However, there was no significant correlation in the methylprednisolone/prednisolone-use population (4G5G vs 5G5G: OR: 0.44; 95% CI: 0. 18 Furthermore, we sought to analyze the cumulative dosage and SONFH risk in mutation carriers by a dose-response meta-analysis. SONFN/ONFH incidence in mutation carriers cannot be obtained when steroids are not used. Therefore, the minimum dose must be used as a reference, and whether increasing the dose affects SONFH incidence must be assessed. Two dose assessment methods were used: an absolute dose and multiple doses. For ABCB1 rs1045642, there was a correlation between the cumulative dose and the occurrence of SONFH in the mutation carriers. As the dose increased, the risk of SONFH decreased (absolute dose: P=0.008; multiple dose: P=0.0064) ( Figure 9A,B) . This negative correction also persists in homozygous mutation (TT genotype, absolute dose: P=0.0368; multiple dose: P=0.0435) and heterozygous mutation patients (CT genotype, absolute dose: P<0.001; multiple dose: P<0.001) (Supplementary Figure S1A-D) . In wild-type populations, there was no correlation between the dose and the risk of SONFH (absolute dose: P=0.098; multiple dose: P=0.098). Although there was no significant correlation between the ABCB1 rs2032582 mutation and the risk of SONFH, the cumulative dose in the mutant population was correlated with the occurrence of SONFH. As the dose increased, the risk of SONFH decreased (absolute dose: P=0.0231; multiple dose: P=0.0258) ( Figure 9C These findings may suggest that mutation carriers still have a risk of SONFH when using steroids, similar to the wild-type population. If SONFH does not occur after small-dose steroid application, increasing the steroid dose in mutation carriers will not increase but rather decrease the risk of SONFH.

For ApoB rs693, although SONFH was correlated with the dose in the mutation carriers (absolute dose: P=0.0064; multiple dose: P=0.0064), a correlation was also found in the wild-type population (absolute dose: P=0.0177; multiple dose: P=0.0013). For ApoB rs1042031, there was no correlation between the cumulative dose and SONFH risk in the mutation carriers (absolute dose: P=0.662; multiple dose: P=0.7078) or in the wild-type population (absolute dose: P=0.3106; multiple dose: P=0.0719). ApoB rs693 and rs1042031 have too few homozygous mutation patients to analyze separately. We did not perform a dose-response analysis for MTHFR rs1801133 and PAI-1 rs1799768 due to the small number of included studies.

ABCB1 rs1045642 mutations have a protective effect against SONFH, and mutations in ApoB rs693 and rs1042031 increase the SONFH risk. The present study further analyzed the effects of each specific characteristic, such as the primary disease of the included patients, the types of steroids, the cumulative dosage, and the treatment duration, on outcomes. However, with the increase in variable factors, the results will be more dependent on the results of single studies, which will reduce the stability of the results. However, the analyses performed here can provide additional clues and research directions for further basic and clinical research studies. In previous research, we analyzed the associations between SNPs and SONFH by traditional meta-analysis, and more than three SNPs were included. Ultimately, ABCB1 rs1045642, rs2032582, ApoB rs693, rs1042031, MTHFR rs1801133, and PAI-1 rs1799768 were included.

In previous meta-analyses, one study analyzed the association between ABCB1 gene polymorphism and SONFH and the 3435T (rs1045642) and 2677T/A (rs2032582) alleles of ABCB1 were identified, which may reduce the risk of GC-induced ONFH [6] . However, our study included more reports and found that rs2032582 had no correlation with SONFH. Our result is the same as that from another meta-analysis that analyzed the correlation between gene polymorphism and SONFH occurrence with the GRADE method to assess the level of evidence. The results showed that ABCB1 rs1045642 could reduce the occurrence of SONFH but had a very low level of evidence. ApoB rs693 and rs1042031 mutations increased the risk of SONFH and had a moderate level of evidence [46] . The results of our study are essentially the same as those of the above two studies, except for rs2032582. However, our work further analyzed the associations between SNP and SONFH based on steroid type, primary disease, drug dose, and applied duration using a mixed-effects logistic regression model. We also conducted a dose-response meta-analysis to analyze the cumulative dosage and SONFH risk in mutation carriers. Our work expanded the research results and found that cumulative steroid dosage and treatment duration had little effect on the results. There was a dose-effect correlation in ABCB1 rs1045642 and rs2032582 mutation carriers. Another study analyzed the associations between ONFH and genetic polymorphisms of vascular endothelial growth factor rs2010963, endothelial nitric oxide synthase rs2070744, and ABCB1 rs1045642. However, the included studies contained SONFH-and ONFH-related studies. The results showed that the mutations of rs2010963 and rs1045642 increased the risk of ONFH, and the mutation of rs2070744 increased the risk of ONFH in an allelic model. There were errors in the present study, as the authors proposed that the mutation of ABCB1 rs1045642 increased the risk of ONFH, but the data showed that it reduced the risk [47] . Our study analyzed SONFH risk, but not ONFH, which is different from the above study.

In the present study, there were significant correlations between the ABCB1 rs1045642 mutant and SONFH in the prednisone-use and methylprednisolone/prednisone-use populations. The ABCB1 rs2032582 mutant homozygote had a protective effect in the methylprednisolone/prednisolone renal transplant population. For ApoB rs693, mutation increased the incidence of SONFH in the prednisone-use and methylprednisolone/prednisolone-use populations and renal transplant patients. For ApoB rs1042031, mutation increased the risk of SONFH in the prednisone-use population. The PAI-1 rs1799768 mutation had a protective effect in the SONFH risk prednisone-use and renal transplant population.

There were different outcomes for the different types of steroid-use populations, possibly because the results were based on even fewer studies. These differences may also indicate that different types of steroids have subtle differences in affecting the risk of SONFH. It is generally believed that prednisone must be transformed to an active form through the liver, while prednisolone does not require this change, which may be the reason for the different results of the two drugs in ABCB1-related mutation states. In addition, more evidence is needed to determine whether the different results for the types of steroids in ApoB and PAI-1 gene polymorphisms were due to different effects on adipocytes.

Currently, a registered study is analyzing the pharmacokinetics of different steroids in children, but the results have not yet been reported (NCT02252237). An earlier study also reported that the cumulative dosage of methylprednisolone, but not prednisone, is correlated with avascular necrosis incidence [48] . Therefore, defined steroid types should be evaluated in further genetic polymorphism-and SONFH-related studies.

Although there are reports about the effect of primary disease on SONFH risk, most have focussed on the influence of steroids used in disease treatment. It was also reported that in renal transplant patients, cyclosporine versus tacrolimus as well as gender factors may be independent risk factors for avascular necrosis [49] . A meta-analysis of the avascular necrosis risk in SLE patients also showed, that in addition to hormonotherapy, arthritic, cushingoid, hypertension, cytotoxic drugs, etc., also affect the occurrence of SONFH [50] . Therefore, the primary disease and the treatment strategy are also possible risk factors for SONFH incidence.

It is generally believed that the occurrence of SONFH is due to long-term and/or high-dosage steroid application. A dose-response meta-analysis showed that the incidence of osteonecrosis is correlated with both the cumulative dose and treatment duration [51] . However, there is still no related research on specific mutant carriers. First, the present study found that the ABCB1 rs1045642 mutation has a protective effect on SONFH, and the risk will be further reduced with increasing cumulative steroid dosage. It can also be concluded that some mutant carriers will have SONFH when the cumulative dosage is low; however, if SONFH does not occur, the risk of SONFH will not increase when the cumulative dose of steroids is increased. Therefore, for mutant carriers, the critical time period for prevention and diagnosis of SONFH is during the early stage of steroid application, not the late stage. A negative correlation was also found in ABCB1 rs2032582 mutation carriers, even though the mutations were not significantly associated with SONFH in general.

The present study still has several limitations. First, our study was performed at the study level instead of the individual level. Second, in stratified analyses, with the increase in variable factors, the results become more dependent on the results of single studies, which will reduce the stability of the results. Third, in dose-response analyses, SONFN/ONFH incidence in mutation carriers without steroid application cannot be calculated; therefore, the minimum dose must be used as a reference. Fourth, our study only analyzed SNPs that were examined in more than three studies, but with emerging research and evidence, more SNPs may be studied in the future. Therefore, the present study was limited by the research available at the time.

The role of autophagy in steroid necrosis of the femoral head: a comprehensive research review

Vascular anatomy and microcirculation of skeletal zones vulnerable to osteonecrosis: vascularization of the femoral head

Incidence of genetic polymorphisms involved in lipid metabolism amongst Chinese patients with osteonecrosis of the femoral head

The pathogenesis of steroid-induced osteonecrosis of the femoral head: a systematic review of the literature

Genetic risk factors for glucocorticoid-induced osteonecrosis: a meta-analysis

Association of ABCB1/MDR1 polymorphisms in patients with glucocorticoid-induced osteonecrosis of the femoral head: evidence for a meta-analysis

ABCB1 gene polymorphisms and glucocorticoid-induced avascular necrosis of the femoral head susceptibility: a meta-analysis

Correlation between steroid-induced osteonecrosis of the femoral head and hepatic CYP3A activity: a systematic review and meta-analysis

Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group

Meta-analysis of epidemiologic dose-response data

Efficacy of glucocorticoids in rodents of severe acute pancreatitis: a meta-analysis

Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses

The SNP rs931794 in 15q25.1 is associated with lung cancer risk: a hospital-based case-control study and meta-analysis

Sleep duration and the risk of osteoporosis amongst middle-aged and elderly adults: a dose-response meta-analysis

Generalized least squares for trend estimation of summarized dose-response data

Polymorphisms in the glucocorticoid receptor gene and associations with glucocorticoid-induced avascular osteonecrosis of the femoral head

Influence of BCL2L11 polymorphism on osteonecrosis during treatment of childhood acute lymphoblastic leukemia

Genetics of glucocorticoid-associated osteonecrosis in children with acute lymphoblastic leukemia

The correlation between apolipoprotein gene polymorphisms and steroid-induced osteonecrosis of the femoral head

Genetic association of the P-glycoprotein gene ABCB1 polymorphisms with the risk for steroid-induced osteonecrosis of the femoral head in Chinese population

MDR1 gene polymorphisms are associated with glucocorticoid-induced avascular necrosis of the femoral head in a Chinese population

Association of apolipoprotein A5 genetic polymorphisms with steroid-induced osteonecrosis of femoral head in a Chinese Han population

Relationship between apolipoprotein B gene polymorphism and tendon-vessel stagnation syndrome of steroid induced osteonecrosis of the femoral head

Genetic polymorphisms in plasminogen activator inhibitor-1 predict susceptibility to steroid-induced osteonecrosis of the femoral head in Chinese population

Association of a polymorphism in PON-1 gene with steroid-induced osteonecrosis of femoral head in Chinese Han population

Association between polymorphisms of steroid induced osteonecrosis of the femoral head and traditional Chinese medicine syndrome types

Relationship between gene polymorphisms and syndrome types in steroid-induced femoral head osteonecrosis

Relationship between gene polymorphisms and hereditary susceptibility of steroid-induced femoral head osteonecrosis in Linyi City of Shandong Province

Polymorphism in the PAI-1 (SERPINE1) gene and the risk of osteonecrosis in children with acute lymphoblastic leukemia

Incidence of genetic polymorphisms involved in lipid metabolism amongst Chinese patients with osteonecrosis of the femoral head

Genetic susceptibility of corticosteroid-induced osteonecrosis of femoral head

Combination analysis of three polymorphisms for predicting the risk for steroid-induced osteonecrosis of the femoral head

A PAI-1 (SERPINE1) polymorphism predicts osteonecrosis in children with acute lymphoblastic leukemia: a report from the Children's Oncology Group

Roles of TNF-alpha gene polymorphisms in the occurrence and progress of SARS-Cov infection: a case-control study

Genetic association of a polymorphism of the cAMP-responsive element binding protein-binding protein with steroid-induced osteonecrosis after kidney transplantation

MDR1(ABCB1) gene polymorphisms associated with steroid-induced osteonecrosis of femoral head in systemic lupus erythematosus

Low molecular weight phenotype of Apo(a) is a risk factor of corticosteroid-induced osteonecrosis of the femoral head after renal transplant

ApoB C7623T polymorphism predicts risk for steroid-induced osteonecrosis of the femoral head after renal transplantation

Thrombophilia and avascular necrosis of femoral head in kidney allograft recipients

Association of corticosteroids and factor V, prothrombin, and MTHFR gene mutations with avascular osteonecrosis in renal allograft recipients

Pharmacogenetic risk factors for osteonecrosis of the hip amongst children with leukemia

Relationship between postrenal transplant osteonecrosis of the femoral head and gene polymorphisms related to the coagulation and fibrinolytic systems in Japanese subjects

Genetic analysis of steroid-induced osteonecrosis of the femoral head

ABCB1 C3435T and G2677T/A polymorphism decreased the risk for steroid-induced osteonecrosis of the femoral head after kidney transplantation

Association of plasminogen activator inhibitor-1 genotype with avascular osteonecrosis in steroid-treated renal allograft recipients

Lipid transporter activity-related genetic polymorphisms are associated with steroid-induced osteonecrosis of the femoral head: an updated meta-analysis based on the GRADE guidelines

VEGF, eNOS, and ABCB1 genetic polymorphisms may increase the risk of osteonecrosis of the femoral head

Avascular necrosis of bone after cardiac transplantation. Prevalence and relationship to administration and dosage of steroids

Cyclosporine use and male gender are independent determinants of avascular necrosis after kidney transplantation: a cohort study

The risk factors of avascular necrosis in patients with systemic lupus erythematosus: a meta-analysis

The authors declare that there are no sources of funding to be acknowledged.

NOS, Newcastle-Ottawa Scale; OR, odds ratio; RA, rheumatoid arthritis; SLE, systemic lupus erythematosus; SNP, single-nucleotide polymorphism; SONFH, steroid-induced osteonecrosis of the femoral head.

Y.J. designed the study; performed statistical analysis; edited and reviewed the manuscript; guarantor of integrity of the entire study. Y.X.G. performed literature research and data acquisition and analysis. J.M. performed literature research and data acquisition and analysis; prepared and reviewed the manuscript. All authors approved the final manuscript.

The authors declare that there are no competing interests associated with the manuscript.