key: cord-0003532-hp8psl4s authors: Chen, Xiantao; Zhang, Leilei; Liang, Dawei; Li, Jing; Liu, Fenzhi; Ma, Hongxia title: 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 date: 2018-12-03 journal: Front Physiol DOI: 10.3389/fphys.2018.01684 sha: a69fd3fb91451692a0d45c306bd5c8c7dbc9debd doc_id: 3532 cord_uid: hp8psl4s Aims: The purpose of this study was to assess the relationship between genetic variants and steroid-induced osteonecrosis of the femoral head (SONFH) in steroid use populations. Methods: We searched the public databases up to April 15, 2018. This study analyzed only the single-nucleotide polymorphisms (SNPs) that have appeared in more than three studies and assessed the level of evidence by classifying the outcomes according to the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) approach. Results: The ABCB1 rs1045642 C>T mutation had a protective effect against SONFH in the allelic model (I(2) = 50.2%; OR: 0.74; 95% CI: 0.55–1.00; p = 0.046). The rs2032582 mutation in the ABCB1 gene showed no relationship to SONFH (allelic model: I(2) = 63.4%; OR: 0.85; 95% CI: 0.58–1.23; p = 0.382). In ApoB rs693, four models showed that mutations can increase SONFH risk, but the allelic model did not. The ApoB rs1042031 mutation increased SONFH risk in the dominant model (I(2) = 50.3%; OR: 2.90; 95% CI: 1.49–5.66; p = 0.002). Conclusion: An allelic model of ABCB1 rs1045642 showed that mutations have a protective effect against SONFH at a very low level of evidence. The mutations in ApoB rs693 and rs1042031 increase the SONFH risk with moderate levels of evidence. Osteonecrosis of the femoral head (ONFH) is caused by the obstruction or interruption of local blood supply, resulting in tissue ischemia, necrosis, and finally bone collapse (Kim et al., 2011) . The core decompression method only relieved the pressure of the pulp cavity and delayed the progression of necrosis (Yu et al., 2018) . At present, long-term and/or large-dose use of steroids have become the most important risk factors for non-traumatic ONFH (Shigemura et al., 2011) . The pathogenesis of steroid-induced ONFH (SONFH) is not yet clear and may be related to an imbalance in lipid metabolism and abnormal microcirculation (Johnson et al., 2004) . The abnormal blood supply leads to apoptosis of osteocytes and osteoblasts, which causes bone loss and reduces bone mineral density (Luo et al., 2018) . In addition, disordered lipid metabolism is another crucial pathogenesis that leads to increases in circulating lipid levels, microvascular fat embolisms, and lipid accumulation in the pulp cavity and ultimately causes osteocyte death. Genetic polymorphisms have been found to be related to SONFH due to individual differences. Studies exploring the association between single-nucleotide polymorphisms (SNPs) and SONFH will help to identify the high SONFH risk population who use steroids. For these patients, it is necessary to avoid steroid application or to change the application strategy, increase the frequency of detection of SONFH, and enable early intervention. Furthermore, genetic variant-related studies can help to understand the pathogenesis of SONFH. Presently, several meta-analyses have been published that indicate that PAI-1 rs1799768 (Gong et al., 2013) , ABCB1 rs1045642 (Gong et al., 2013; Li et al., 2014; Zhou et al., 2015) , and CYP3A variants (Guo and Deng, 2017) are related to the occurrence of SONFH. However, there is still a dispute over ABCB1 rs2032582 (Gong et al., 2013; Li et al., 2014; Zhou et al., 2015) . This study will explicitly exclude studies in which the control group includes a healthy population or other types of ONFH; only comparison studies between ONFH and non-ONFH in the steroid use population are included. Because all studies are based on casecontrol and cohort designs with low evidence levels, we also used the Grading of Recommendations Assessment, Development, and Evaluation (GRADE) guidelines to assess the outcomes (Guyatt et al., 2008) . Overall, the relationship between genetic variants and steroid-induced osteonecrosis of the femoral head in steroid use populations will be assessed in this study. This meta-analysis was performed in accordance with the Preferred Reporting Items for Systematic Reviews (PRISMA) guidelines (Moher et al., 2009 ). Two authors independently performed a literature search using PubMed, Embase, the Cochrane Library, and the Chinese public databases, including the China National Knowledge Infrastructure (CNKI), the China Biology Medicine (CBM) Database, the China Science Periodical Database (CSPD, Wanfang Database), and the VIP Journal Integration Platform (VJIP). Searches were performed for studies published up to April 15, 2018. The following terms were used in the search strategy: hormone, glucocorticoid, steroid, corticosteroid, osteonecrosis, necrosis, femoral, femur, femoris, whirlbone, polymorphism, SNP, genetic, mutation, genotype, allele, allelic, and variation. We also conducted manual searches of the reference lists of relevant reviews to avoid omissions. The following studies were included in our meta-analysis if they fulfilled the inclusion criteria: (1), study has a case-control or cohort study design; (2), steroid-using patients are included; (3), study compares ONFH and non-ONFH patients after steroid use; (4), study assesses the association between SNPs and SONFH; and (5), study indicates the frequencies of specific alleles or the effect sizes of individual genotypes between cases and controls. The exclusion criteria included the following: (1), study compares patients with SONFH to healthy populations or populations with other types of ONFH; (2), study is about family heredity; (3), study is not SNP-related; and (4), study does not report data pertaining to allelic frequencies or calculable effect size. In addition, conference reports, editor comments, reviews, and academic dissertations were also excluded from the analysis. Two authors independently extracted the following information from each eligible study: the first author's name, publication year, research location, sample size, average age, primary disease, diagnostic mode, and genes of interest. The methodological quality of the included studies was evaluated using the Newcastle-Ottawa Scale (NOS), a validated tool for evaluating the quality of observation studies that includes the following 3 subscales: selection, comparability, and exposure (Stang, 2010) . We assessed all results for the level of evidence using the GRADE approach (Andrews et al., 2013) , the average NOS score, and the number of patients included and then presented the data in a threedimensional Manhattan plot. The Hardy-Weinberg equilibrium (HWE) of the included populations was calculated to assess the consistencies of allele frequencies between generations (Chen and Chatterjee, 2007) . Association analysis was performed using five genetic models (Areeshi et al., 2013) . The effect size was estimated by calculating the summary odds ratio (OR) and its 95% confidence intervals (CIs). The I 2 statistic was used to estimate the degree of heterogeneity among the studies. If the I 2 ≥ 50% (Q test, p < 0.1), a random-effect model was used; otherwise, a fixed-effect model was used. Subgroup analysis was also performed according to the ethnicity and location of the population. We assessed publication bias using the Begg and Egger tests. All tests were two-tailed, and a p-value of less than 0.05 was deemed statistically significant. We analyzed the data using the R program (Version 3.3.1) and STATA (Version 14.0). We used GRADE to assess the level of evidence with four levels graded from high (best) to very low (worst). All objective studies used an observational study design, which downgrades the quality of evidence. Inconsistency was evaluated using the I 2 and p-values of the Q test, and the outcome was imprecise if the standard error (SE) was larger than 0.2. Publication bias (reporting bias) also reduced the level of evidence. Finally, a large (OR > 2 or <0.5) or very large (OR > 5 or <0.2) effect upgraded the quality of evidence. Our research returned 278 English articles and 285 Chinese articles after removing duplications. After screening the titles and abstracts, 482 of these articles were excluded. The full texts of 81 articles were assessed, among which the control group did not receive steroid therapy (25); the article was a review (7); the study was non-SNP-related (5); the research was basic (3); the study did not report frequencies or effect size results (2); the study did not include steroid-related ONFH (2); the study was a case report (2); the research was microRNA related (2); the report was a duplicate (2); and the research was of heredity SONFH (1). Finally, we collected 30 trials assessing 7,553 patients for our meta-analysis (Table 1, Figure 1) . All studies used a case-control design except one that used a cohort design (Karol et al., 2015) . In addition, Weibao Fang's two studies (Wei, 2011a,b) , and Kuribayashi and Tamura's studies (Tamura et al., 2007; Kuribayashi et al., 2008) included the same patient population but did not assess the same SNPs. We analyzed different SNP results and excluded duplicated results. The NOS scores were six to nine points, and the overall quality was ideal in observational studies (Table 1) . Six SNPs in four genes were included in the meta-analysis. In the ABCB1 gene, rs1045642 is also known as C3435T, is located in the coding region, and is a synonymous mutation. In this research, the pooled results for rs1045642 showed that the C > T mutation protected against SONFH in the allelic model (I 2 = 50.2%; OR: 0.74; 95% CI: 0.55-1.00; p = 0.046) (Figure 2, Table 2 ). In the ApoB gene, rs693 is located in the coding region. The results showed that there was no significant relationship between this mutation and SONFH in the allelic model (I 2 = 58%; OR: 2.63; 95% CI: 0.92-7.54; p = 0.072). However, other models showed that this mutation could increase SONFH risk (heterozygous model: I 2 = 54.5%; OR: 2.46; 95% (Figure 2 , Table 2 ). The level of evidence of the above positive results was assessed according to GRADE guidelines. The allelic model of ABCB1 rs1045642 showed a protective effect against SONFH. In these results, 1,141 patients were included; the average NOS score was 7.88. Although the SE was relatively small, there was heterogeneity (I 2 = 50.2%, p = 0.02). Thus, the level of evidence was very low (Figure 3 ). The heterozygous model and dominant model of ApoB rs693 showed that the mutation increases SONFH risk. The SE of the two models was smaller than 0.2, and heterogeneity was not obvious (heterozygous model: I 2 = 54.5%, p = 0.111; dominant model: I 2 = 31.40%, p = 0.212). The level of evidence was moderate because of the large effect (OR>2) (Figure 3) . The homozygous and recessive models of ApoB rs693 also showed that the mutation increases disease risk without heterogeneity (homozygous model: I 2 = 24.4%, p = 0.25; recessive model: I 2 = 32.1%, p = 0.225). Although the SEs were larger than 0.2, the OR was greater than 5. Thus, the evidence levels were also moderate (Figure 3) . In addition, the average NOS scores of these four models were 8 to 8.25. The dominant model of ApoB rs1042031 indicated that the mutation increases SONFH risk without heterogeneity (I 2 = 50.3%, p = 0.11). The SE of the result was less than 0.2, and the OR was greater than 2. Thus, the level of evidence was moderate (Figure 3) . However, the genetic linkage disequilibrium that exists in the above positive results might affect stability. Our study analyzed the association between genetic polymorphisms and SONFH risk by comparing ONFH and non-ONFH in steroid use populations. In the SNP results of the included studies, more than three allelic models of ABCB1 rs1045642 showed that mutations protect against SONFH. The ApoB rs693 and rs1042031 mutations increase SONFH risk. According to the GRADE guidelines, the evidence levels for ApoB rs693 and rs1042031 were moderate, and ABCB1 rs1045642 was very low. In a previous meta-analysis, Gong et al analyzed 23 studies and 35 genes and indicated that PAI-1 rs1799768 and ABCB1 rs1045642, but not MTHFR rs1801133 or ABCB1 rs2032582, were related to the risk of osteonecrosis (Gong et al., 2013) . However, in this study, the control population included a healthy group and a population with other types of ONFH, which might affect the accuracy. In our study, there was no relationship between PAI-1 rs1799768 after adjusting for the control population. Guo et al assessed the correlation between SONFH and hepatic CYP3A activity and high found that CYP3A activity could reduce SONFH risk. However, the conclusion was based on animal model results, and it is unclear whether the risk is reduced in humans (Guo and Deng, 2017) . ABCB1 is a member of a superfamily of ATP-binding cassette transporters that transport variant molecules across cellular membranes. Zhou et al. researched the association between ABCB1 polymorphisms and SONFH and found that rs1045642 and rs2032582 are associated with SONFH (Zhou et al., 2015) . However, this study also included a healthy population in the control population. In contrast, a Chinese review showed that rs1045642 could reduce SONFH risk, while rs2032582 is not related to SONFH . For rs1045642, our results are the same as those of a previous study, but with a very low level of evidence. The mutation does not cause changes in the amino acid sequence; thus, the specific mechanism of SONFH risk change caused by this mutation is still unknown. However, this mutation has been shown to be related to the metabolism of some drugs and to the prognosis of some cancers. This relationship may be related to differences in the characteristics of nonexpressed proteins, such as codon bias, expression efficiency, and differences in mRNA characteristics. Therefore, this mutation is more suitable for clinical use as a predictor of SONFH risk for steroid users to guide individual drug application. People with the rs1045642 T allele have a low risk of SONFH with steroid application, while individuals with C allele have a high SONFH risk and should avoiding long-term high-dose steroid application, with increased duration of follow-up and frequency of diagnosis. For rs2032582, our results found that the mutation is not related to SONFH (allelic model: p = 0.382) with statistical heterogeneity (I 2 = 63.4%, p = 0.012). Additionally, we performed a post hoc analysis of Seth E. Karol's research that excluded SNP results with very low-level evidence (Karol et al., 2015) . Then, we matched potential SNPs with GO/KEGG annotations ( Table 3) . According to a gene set enrichment analysis based on overrepresentation enrichment analysis, only the molecular functions of ABCB1 and ApoB showed obvious relationships with lipid transporter activity (GO: 0005319) (Zhang et al., 2005) . This finding indicates that at present, SNPs involved in lipid transporter activity are more related to SONFH risk than other SNPs, which supports the lipid metabolism disorder theory. ApoB is an important factor in normal lipid metabolism, the mutation of which can also be used as a predictor of SONFH occurrence to assess the risk of SONFH and to guide steroid clinical application. Although rs693 is a synonymous variant, rs693 is related to the circulating concentration of LDL cholesterol. T allele carriers have high levels of TG, TC, and LDL-C and low HDL-C (Sandhu et al., 2008) . In addition, rs1042031 is a coding sequence variant that causes glutamate to lysine mutation (C > T), or the introduction of a stop codon in the amino acid sequence (C > A). The rs1042031 mutation is important for regulating the binding of apolipoprotein B to the LDL cholesterol receptor (Liu et al., 2013) . In this work, these two mutations were found to be related to SONFH. Therefore, we hypothesized that the increase of circulating LDL induced by mutation increased the risk of SONFH, supporting the hypothesis that lipid metabolism disorders are associated with SONFH. Abnormal lipid metabolism leads to bone marrow adipogenesis of the femoral head that is mainly due to hyperlipidemia and inhibition of mitochondrial dehydrogenation by steroid application. The abnormal lipid metabolism then leads to fat embolism, which blocks local blood perfusion and exacerbates local inflammation. Therefore, the question remains: does reducing circulating LDL effectively reduce the risk of SONFH with steroid application? In animal models, it has been suggested that reduced circulating lipid levels can reduce the risk of SONFH (Yang Z. et al., 2016) , but the clinical effect remains to be confirmed. Additionally, in the location subgroup analysis of rs1042031 results, the heterogeneity decreased significantly from 50 to 0%. The pool results of the three Chinese studies had low heterogeneity, and the population was concentrated in Guangxi and Shandong in China (Wei, 2011a; Zeng et al., 2014; Wei et al., 2015) . Another study analyzed the Japanese population (Hirata et al., 2007b) . Therefore, the centralization of the included population may be the main reason for the reduction in heterogeneity. The present study still has several limitations. First, our study was performed at the study level instead of the individual level. Second, our study did not consider the impact of primary disease or steroid therapy strategies on the results. Third, there is a large amount of heterogeneity in the results of this study. Although random effect models are used for heterogeneous results, these heterogeneities can still have a potential impact on the results. Fourth, our study only analyzed the SNPs that were examined in more than three studies, but with new research and evidence, more SNPs may be studied. Therefore, this study was limited by the research available at the time. An allelic model of ABCB1 rs1045642 showed that mutations had a protective effect against SONFH at a very low level of evidence. The mutations in ApoB rs693 and rs1042031 increased the SONFH risk with moderate levels of evidence. XC performed the conception and design of the work. XC and LZ drafted and revised the work critically. DL and JL analyzed data for work. FL and HM acquisited data for this work. All authors provide approval for publication of the content. GRADE guidelines: 14. Going from evidence to recommendations: the significance and presentation of recommendations CD14−159 C>T gene polymorphism with increased risk of tuberculosis: evidence from a meta-analysis ABCB1 C3435T and G2677T/A polymorphism decreased the risk for steroid-induced osteonecrosis of the femoral head after kidney transplantation Genetic analysis of steroid-induced osteonecrosis of the femoral head Relationship between postrenal transplant osteonecrosis of the femoral head and gene polymorphisms related to the coagulation and fibrinolytic systems in Japanese subjects Polymorphism in the PAI-1 (SERPINE1) gene and the risk of osteonecrosis in children with acute lymphoblastic leukemia Association of corticosteroids and factor V, prothrombin, and MTHFR gene mutations with avascular osteonecrosis in renal allograft recipients Exploiting Hardy-Weinberg equilibrium for efficient screening of single SNP associations from case-control studies Association of apolipoprotein A5 genetic polymorphisms with steroid-induced osteonecrosis of femoral head in a Chinese Han population Thrombophilia and avascular necrosis of femoral head in kidney allograft recipients Association of plasminogen activator inhibitor-1 genotype with avascular osteonecrosis in steroid-treated renal allograft recipients A PAI-1 (SERPINE1) polymorphism predicts osteonecrosis in children with acute lymphoblastic leukemia: a report from the Children's Oncology Group Genetic risk factors for glucocorticoid-induced osteonecrosis: a meta-analysis Correlation between steroid-induced osteonecrosis of the femoral head and hepatic CYP3A activity: a systematic review and meta-analysis GRADE: an emerging consensus on rating quality of evidence and strength of recommendations Incidence of genetic polymorphisms involved in lipid metabolism among Chinese patients with osteonecrosis of the femoral head Genetic susceptibility of corticosteroid-induced osteonecrosis of femoral head ApoB C7623T polymorphism predicts risk for steroid-induced osteonecrosis of the femoral head after renal transplantation Low molecular weight phenotype of Apo(a) is a risk factor of corticosteroid-induced osteonecrosis of the femoral head after renal transplant Vascular anatomy and microcirculation of skeletal zones vulnerable to osteonecrosis: vascularization of the femoral head Genetics of glucocorticoid-associated osteonecrosis in children with acute lymphoblastic leukemia Significant associations of PAI-1 genetic polymorphisms with osteonecrosis of the femoral head Combination analysis of three polymorphisms for predicting the risk for steroid-induced osteonecrosis of the femoral head Association between polymorphisms of steroid induced osteonecrosis of the femoral head and traditional Chinese medicine syndrome types ABCB1 gene polymorphisms and glucocorticoid-induced avascular necrosis of the femoral head susceptibility: a meta-analysis Associations of polymorphisms of rs693 and rs1042031 in apolipoprotein B gene with risk of breast cancer in Chinese The role of autophagy in steroid necrosis of the femoral head: a comprehensive research review Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement Influence of BCL2L11 polymorphism on osteonecrosis during treatment of childhood acute lymphoblastic leukemia Pharmacogenetic risk factors for osteonecrosis of the hip among children with leukemia LDL-cholesterol concentrations: a genome-wide association study Incidence of osteonecrosis associated with corticosteroid therapy among different underlying diseases: prospective MRI study Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses Genetic association of a polymorphism of the cAMPresponsive element binding protein-binding protein with steroid-induced osteonecrosis after kidney transplantation Roles of TNF-alpha gene polymorphisms in the occurrence and progress of SARS-Cov infection: a case-control study Association of a polymorphism in PON-1 gene with steroid-induced osteonecrosis of femoral head in Chinese Han population 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 The correlation between apolipoprotein gene polymorphisms and steroid-induced osteonecrosis of the femoral head MDR1 gene polymorphisms are associated with glucocorticoid-induced avascular necrosis of the femoral head in a Chinese population Correlation between ABCB1C3435T and G2677T/A polymorphism and steroid-induced osteonecrosis of femoral head:a meta-analysis MDR1(ABCB1) gene polymorphisms associated with steroid-induced osteonecrosis of femoral head in systemic lupus erythematosus The efficacy of statins in preventing glucocorticoid-related osteonecrosis in animal models: a meta-analysis Effectiveness of various hip preservation treatments for non-traumatic osteonecrosis of the femoral head: a network meta-analysis of randomized controlled trials Relationship between apolipoprotein B gene polymorphism and tendon-vessel stagnation syndrome of steroid induced osteonecrosis of the femoral head WebGestalt: an integrated system for exploring gene sets in various biological contexts Genetic association of the P-glycoprotein gene ABCB1 polymorphisms with the risk for steroid-induced osteonecrosis of the femoral head in Chinese population Genetic polymorphisms in plasminogen activator inhibitor-1 predict susceptibility to steroid-induced osteonecrosis of the femoral head in Chinese population Polymorphisms in the glucocorticoid receptor gene and associations with glucocorticoidinduced avascular osteonecrosis of the femoral head Association of ABCB1/MDR1 polymorphisms in patients with glucocorticoid-induced osteonecrosis of the femoral head: evidence for a meta-analysis