Dongqi Zhou
Changyan Zi
Gaofeng Gan1
Qiu Chen- 1Department of Traditional Chinese Medicine, Sichuan Taikang Hospital, Chengdu, Sichuan, China
- 2Department of Good Clinical Practice (GCP), Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 3Department of Endocrine, Hospital of Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
Background: There are complex interactions between osteoporosis and the immune system, and it has become possible to explore their causal relationship based on Mendelian randomization methods.
Methods: Utilizing openly accessible genetic data and employing Mendelian randomization analysis, we investigated the potential causal connection between 731 immune cell traits and the risk of developing osteoporosis.
Results: Ten immune cell phenotypes were osteoporosis protective factors and three immune cell phenotypes were osteoporosis risk factors. Specifically, the odds ratio (OR) of IgD+ CD24+ %B cell (B cell panel) risk on Osteoporosis was estimated to be 0.9986 (95% CI = 0.9978~0.9996, P<0.01). The OR of CD24+ CD27+ %B cell (B cell panel) risk on Osteoporosis was estimated to be 0.9991 (95% CI = 0.9984~0.9998, P = 0.021). The OR of CD33- HLA DR+AC (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9996 (95% CI = 0.9993~0.9999, P = 0.038). The OR of EM CD8br %CD8br (Maturation stages of T cell panel) risk on Osteoporosis was estimated to be 1.0004 (95% CI = 1.0000~1.0008, P = 0.045). The OR of CD25 on IgD+ (B cell panel) risk on Osteoporosis was estimated to be 0.9995 (95% CI = 0.9991~0.9999, P = 0.024). The OR of CD25 on CD39+ activated Treg+ (Treg panel) risk on Osteoporosis was estimated to be 1.001 (95% CI = 1.0001~1.0019, P = 0.038). The OR of CCR2 on CD62L+ myeloid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9992 (95% CI = 0.9984~0.9999, P = 0.048). The OR of CCR2 on CD62L+ plasmacytoid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9993 (95% CI = 0.9987~0.9999, P = 0.035). The OR of CD45 on CD33dim HLA DR+ CD11b- (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9988 (95% CI = 0.9977~0.9998, P = 0.031). The OR of CD45 on Mo MDSC (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9992 (95% CI = 0.9985~0.9998, P = 0.017). The OR of SSC-A on B cell (TBNK panel) risk on Osteoporosis was estimated to be 0.9986 (95% CI = 0.9972~0.9999, P = 0.042). The OR of CD11c on CD62L+ myeloid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9987 (95% CI = 0.9978~0.9996, P<0.01). The OR of HLA DR on DC (cDC panel) risk on Osteoporosis was estimated to be 1.0007 (95% CI = 1.0002~1.0011, P<0.01). No causal effect of osteoporosis on immune cells was observed.
Conclusions: Our study identified 13 unreported immune phenotypes that are causally related to osteoporosis, providing a theoretical basis for the bone immunology doctrine.
1 Introduction
Osteoporosis is a metabolic bone disease characterized by changes in the microstructure of bones, decreased bone density, and increased bone fragility (1). In recent years, the complex interaction between the immune system and osteoporosis has received attention. With the introduction of the concept of “immunoporosis” (2), the correlation between immune cells and osteoporosis has become a research hotspot.
A multitude of research indicates that the immune system, encompassing both its innate and adaptive branches, is critically involved in both the onset and persistence of the inflammatory state characteristic of osteoporosis (3–5). Cells of the immune system, along with an array of chemokines and cytokines, affect the metabolic dynamics of bones, influencing the growth, differentiation, and equilibrium of osteoclasts and osteoblasts (6, 7). Bone marrow macrophages (BMMs) can switch between M1 inflammatory state and M2 repair state depending on the surrounding environment (8). Research reports that M1 macrophages can disrupt osteoblast formation by secreting tumor necrosis factor-α (TNF-α) (9). Certain subsets of pro-inflammatory T cells may promote osteoblast differentiation during the early stages of fracture healing, but it is currently unclear whether T cells will exert the same effects in the persistent inflammatory environment of osteoporosis (10). B cells play a role in regulating the receptor activator of nuclear factor kappa-B ligand/osteoprotegerin (RANKL/OPG) axis, which affects osteoclast differentiation. They can inhibit osteoclast differentiation by secreting OPG (1). Breg cells, a subset of B cells, can inhibit the differentiation of Th17 cells by secreting interleukin-10 (IL-10) and transforming growth factor-beta (TGF-β). This inhibition leads to a decrease in the secretion of inflammatory factors by monocytes and dendritic cells (11). so, promoting B cell differentiation can have immunosuppressive effects and may also promote the maturation and growth of osteogenic precursor cells. It is hypothesized that there is a causal relationship between immune cells and osteoporosis. However, due to limitations in sample size, study design flaws, and confounding factors, it is currently not possible to clearly elucidate this relationship (12–15). Therefore, further research is necessary to clarify the causal relationship between immune cells, including B cells, and osteoporosis.
The method of Mendelian randomization (MR) combines instrumental variables (IVs) and Mendel’s laws of inheritance (16). It overcomes limitations of traditional randomized controlled trials by minimizing the impact of confounding factors (17). Compared to observational analyses, MR offers several advantages: (1) it maintains the proper sequence of exposure and outcome, preventing reverse causality; (2) it reduces the influence of confounding factors on the results; (3) it enables long-term studies of exposure and outcome (IVs remain valid throughout a lifetime); (4) it avoids regression dilution resulting from testing errors, thanks to the accurate genotype testing (18, 19). In this study, a comprehensive two-sample bidirectional MR analysis was conducted to establish the causal relationship between immune cell characteristics and osteoporosis.
2 Materials and methods
2.1 Study design
We employed a two-sample MR analysis to investigate the causality between 731 immune cell characteristics and osteoporosis. The IVs used in the analysis satisfied the following critical assumptions: genetic variation exhibits a direct association with exposure factors; genetic variation remains independent of possible confounders; and genetic variation does not affect outcome through pathways other than exposure.
2.2 Genome-wide association study data sources
The GWAS summary statistics for osteoporosis were obtained from the IEU openGWAS public database (Osteoporosis GWAS id: ebi-a-GCST90038656). The study included GWAS data from 484,598 individuals, with 7,751 cases and 476,847 controls. For each immune trait, GWAS summary statistics are available in the GWAS Catalog with accession numbers ranging from GCST0001391 to GCST0002121 (20, 21). The dataset includes a total of 731 immunophenotypes, which consist of various features such as absolute cell (AC) counts (118 traits), median fluorescence intensities (MFI) reflecting surface antigen levels (389 traits), morphological parameters (MP) (32 traits), and relative cell (RC) counts (192 traits). The MFI, AC, and RC features encompass various immune cell types, including B cells, CDCs, mature stages of T cells, monocytes, myeloid cells, TBNK (T cells, B cells, natural killer cells), and Treg panels. The MP feature includes CDC and TBNK panels. These immunophenotypes provide valuable information for investigating the relationship between immune cell characteristics and osteoporosis through two-sample MR analysis.
2.3 Screening mendelian randomized IVs
In the first step of the analysis, the IVs was set at 1×10-5, based on previous related studies (20, 22). The second step involved removing single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) with each other, with a parameter setting of r2 = 0.001 and kb=10,000. This step helps to avoid bias caused by the correlation between SNPs in close proximity to one another. Finally, if any of the selected SNPs were found to be associated with confounding factors that significantly correlated with the outcome (p value <5×10−8), they were removed from the list of selected SNPs. This step is important to ensure that the observed associations are not due to the influence of other factors that may affect both the genetic variant and the outcome of interest (23). The R package used in this study is available from a public database (URL: https://github.com/ZDQZBXZ/731-Immune-Cell-Code).
2.4 Statistical analysis
The analysis process was conducted using R 3.5.3 software (http://www.Rproject.org). A Mendelian Randomization (MR) analysis employing a two-sample approach was carried out on GWAS data to investigate the causal link between 731 immune cell types (as exposures) and osteoporosis (as the outcome), with the Wald ratio (WR) combined inverse variance weighting (IVW) strategy as the principal analytic method (24). Where MR results showed heterogeneity, the random effects version of IVW was applied; in its absence, a fixed effects model was adopted. To enhance the robustness and dependability of the findings, additional methods such as MR-Egger regression, Weighted Median Regression (WMR), and MR Pleiotropy RESidual Sum and Outlier (MR-PRESSO) were used alongside IVW (25). Specifically, the IVW method aims to assess the overall population by combining the effect estimates of each individual SNP in MR studies. It’s important to note that if a single SNP contributes more than 50% of the weight in all relevant SNP pairs, or if multiple SNPs collectively exceed half of the weight, it can lead to horizontal pleiotropy. Therefore, in this study, the WMR method was employed as an alternative strategy to address situations where IVW findings might lack precision (26). Indeed, the IVW method assumes its intercept must pass through zero (27). However, this assumption may not accommodate scenarios where the intercept deviates from zero. To address this limitation, the MR-Egger method was employed as an additional strategy (28). The use of the MR-PRESSO method aimed to detect and correct biases arising from pleiotropy (29). In this study, data adjusted for pleiotropy via MR-PRESSO were incorporated to ensure potential outliers were appropriately addressed. F-statistics were calculated as measures of IV strength, with values greater than 10 indicating minimal susceptibility to weak instrument bias (30). IVs failing to meet the criterion of F ≥ 10 were excluded, ensuring robust IV strength. Sensitivity checks including the Cochran’s Q test, MR-Egger intercept test, and MR-PRESSO global test were performed, ensuring the solidity of the results, evidenced by all MR-Egger intercept tests having p-values over 0.05, thus excluding the presence of a horizontal pleiotropy effect (31). In the event of inconsistency in MR findings, a reanalysis with stricter cutoffs for the p-value of instrumental variables and more rigorous selection criteria was executed. Scatter plots offered a graphical display of the reevaluated data, while funnel plots were employed to scrutinize any potential publication bias.
3 Results
In this study, we used a two-sample bidirectional MR method to explore the relationship between 731 immune cell phenotypes and osteoporosis. Regarding the causal analysis of immune cells for osteoporosis, we observed a total of 13 different immune cell phenotypes that have an effect on osteoporosis. Of these 3 immune cell phenotypes were osteoporosis risk factors and 10 immune cell phenotypes provided osteoporosis protection (Figure 1; Supplementary File 1). The results of IVW analysis were referenced as the primary assessment tool for the study. The IVW results for 13 immune cell phenotypes are shown in Table 1. Specifically, the odds ratio (OR) of IgD+ CD24+ %B cell (B cell panel) risk on Osteoporosis was estimated to be 0.9986 (95% CI = 0.9978~0.9996, P<0.01). The OR of CD24+ CD27+ %B cell (B cell panel) risk on Osteoporosis was estimated to be 0.9991 (95% CI = 0.9984~0.9998, P = 0.021). The OR of CD33- HLA DR+AC (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9996 (95% CI = 0.9993~0.9999, P = 0.038). The OR of EM CD8br %CD8br (Maturation stages of T cell panel) risk on Osteoporosis was estimated to be 1.0004 (95% CI = 1.0000~1.0008, P = 0.045). The OR of CD25 on IgD+ (B cell panel) risk on Osteoporosis was estimated to be 0.9995 (95% CI = 0.9991~0.9999, P = 0.024). The OR of CD25 on CD39+ activated Treg+ (Treg panel) risk on Osteoporosis was estimated to be 1.001 (95% CI = 1.0001~1.0019, P = 0.038). The OR of CCR2 on CD62L+ myeloid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9992 (95% CI = 0.9984~0.9999, P = 0.048). The OR of CCR2 on CD62L+ plasmacytoid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9993 (95% CI = 0.9987~0.9999, P = 0.035). The OR of CD45 on CD33dim HLA DR+ CD11b- (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9988 (95% CI = 0.9977~0.9998, P = 0.031). The OR of CD45 on Mo MDSC (Myeloid cell panel) risk on Osteoporosis was estimated to be 0.9992 (95% CI = 0.9985~0.9998, P = 0.017). The OR of SSC-A on B cell (TBNK panel) risk on Osteoporosis was estimated to be 0.9986 (95% CI = 0.9972~0.9999, P = 0.042). The OR of CD11c on CD62L+ myeloid DC (cDC panel) risk on Osteoporosis was estimated to be 0.9987 (95% CI = 0.9978~0.9996, P<0.01). The OR of HLA DR on DC (cDC panel) risk on Osteoporosis was estimated to be 1.0007 (95% CI = 1.0002~1.0011, P<0.01) (Figure 2; Supplementary File 2). Sensitivity analysis of 13 immune cells to osteoporosis was uploaded as a supplementary file (Supplementary File 3). There was no gene pleiotropy for 13 immune cell phenotypes, indicating that the results are reliable (Supplementary File 4). Publication bias and funnel plot results are also presented as supplementary documents (Supplementary File 5). No causal effect of osteoporosis on immune cells was observed (Supplementary File 6). The raw data of 731 immune cell phenotypes after excluding weak IVs, along with F-statistics for 13 immune cell phenotypes, are detailed in Supplementary File 7.
Figure 1 The causal relationship between 13 immune cell phenotypes and osteoporosis. An upward slope indicates exposure as a risk factor and a downward slope indicates exposure as a protective factor. The outcome is osteoporosis. (A) exposure factor is CCR2 on CD62L+ myeloid DC. (B) exposure factor is CCR2 on CD62L+ plasmacytoid DC. (C) exposure factor is CD11c on CD62L+ myeloid DC. (D) exposure factor is CD24+ CD27+ %B cell. (E) exposure factor is CD25 on CD39+ activated Treg. (F) exposure factor is CD25 on IgD+. (G) exposure factor is CD33- HLA DR+ AC. (H) exposure factor is CD45 on CD33dim HLA DR+ CD11b-. (I) exposure factor is CD45 on Mo MDSC. (J) exposure factor is EM CD8br %CD8br. (K) exposure factor is HLA DR on DC. (L) exposure factor is IgD+ CD24+ %B cell. (M) exposure factor is SSC-A on B cell.
Table 1 The IVW results for 13 immune cell traits.
Figure 2 The forest plot of the causal relationship between 13 immune cell phenotypes and osteoporosis. The inverse variance weighting method is used as the primary outcome p < 0.05 is considered causal relationship.
4 Discussion
Utilizing the genetic information from publicly accessible databases, our investigation delved into the potential causative links between 731 immune cell characteristics and osteoporosis via MR. As far as we are aware, Mendelian studies are the current focus of much scientific attention; however, the investigation into the connection between immune cells and diseases remains largely unexplored. Our MR analysis is pioneering in examining the causative links between a variety of immune phenotypes and osteoporosis, thereby laying some groundwork for the concept of immune-related osteoporosis. Our research identified 13 immune cell phenotypes spanning four categories of immune traits (MFI, RC, AC, and MP) that have a causal relationship with osteoporosis. Within these traits, 3 phenotypes were identified as risk factors, while 10 were established as protective factors, all with a P value of less than 0.05.
According to MR results, the traits IgD+ CD24+ %B cells, CD24+ CD27+ %B cells, and CD25 on IgD+ (B cell panel) have been identified as protective factors against osteoporosis. These three phenotypes all reflect the expression of surface-active substances on B cells. Typically, common markers found on the surface of B cells are utilized to isolate and characterize the various subgroups within the B cell population. To elaborate, IgD is a membrane-bound immunoglobulin on B cells that binds to antigens and participates in the B cell immune response (32, 33). CD24 is a cell surface marker involved in cell-to-cell interactions and interactions between cells and their external microenvironment, and it is associated with the adhesion and migration capabilities of B cells (33–35). The term IgD+ CD24+ %B cells indicates the percentage of B cell subgroups that co-express IgD and CD24 within the overall B cell population (the same principle applies to the other phenotypic characteristics mentioned in the text). Consequently, MR analyses suggest that a higher proportion of B cell subgroups expressing IgD and CD24 correlates with a reduced risk of osteoporosis. IgD plays a crucial role within the adaptive immune system, participating in the B-cell driven antibody-mediated immune reactions (36). Studies indicate that higher IgD levels and lower IgM levels on B cells correspond with reduced B cell activity in human peripheral blood (36–40). However, the connection between these observations and osteoporosis remains unclear. It’s hypothesized that the suppressive activity of B cells may interact with the development and expansion of bone-forming progenitor cells (11). Whether the control of B cell functionality through IgD surface expression is implicated in bone metabolism requires comprehensive investigation. CD24 participates in the osteogenic differentiation process of bone marrow mesenchymal stem cells (BMSCs) and acts as a signaling molecule within this context (41, 42). Study has demonstrated that CD24, located on the cell surface, displays a notable disparity in expression levels when BMSCs undergo differentiation into osteogenic or lipogenic lineages, with reduced CD24 expression being associated with diminished osteogenic differentiation (43). CD24 is currently proposed as a selective biomarker for a subpopulation of BMSCs with enhanced osteogenic potential (33). CD27, on the other hand, is a molecular marker on the cell membrane commonly used to differentiate various subgroups of memory B cells (44, 45). A previous MR study has reported a causal relationship between CD27 and bone density (46). CD25, also known as the IL-2 receptor alpha chain, is another cell surface marker that characterizes a cell surface receptor and plays a crucial role in regulating the activation, proliferation, and differentiation of T and B cells (47, 48). Evidence on CD25 expression on the surface of B cells and osteoporosis is limited, so this finding is revealing.
The term “SSC-A on B cell (TBNK panel)” refers to the quantification and assessment of the side scatter signal intensity (SSC-A) observed in B cells within the TBNK panel for cell analysis. SSC-A indicates the intensity of side scatter signals detected during flow cytometry, providing insights into cellular dimensions, shape, and intracellular complexity (21). Generally, a heightened SSC-A reading may indicate a greater cell population, though it does not directly quantify the actual number of cells. So, MR results from this research demonstrate a negative causal link between side scatter signal intensity on B cells (SSC-A) and osteoporosis. Moreover, the immunosuppressive characteristics of B cells, along with the presentation of certain proteins on their surface, play a significant role in providing defense against osteoporosis. However, further discussion is needed on the two main aspects of the potential relationship between B cells and osteoporosis: an osteoprotective function in physiological conditions and a bone-destructive function in pathological states. In the bone marrow microenvironment, B lymphocytes are unable to mature in the absence of secretory factors from mesenchymal and osteoblastic cells (49). Bruton’s tyrosine kinase (BTK) is an enzyme integral to signaling in pre-B cell receptors (pre-BCRs) (50). As part of their maturation process, B cells participate in the assembly of pre-BCRs via V(D)J recombination (49). Upon successful pre-BCR formation, BTK engages in subsequent signaling pathways (51). Interestingly, BTK signaling unexpectedly plays a crucial role in the differentiation of osteoclasts. When BTK signaling is impaired, osteoclasts become dysfunctional and bone resorption is diminished (52). OPG, belonging to the tumor necrosis factor (TNF) receptor superfamily, competitively interacts with RANKL, inhibiting the interaction of RANKL with its receptor RANK—a critical step in osteoclastogenesis and activation (53). There is previous evidence that B cells under physiological conditions can secrete OPG (54), and in the bone marrow, B cells contribute nearly half of the OPG (55), suggesting that B lymphocytes are crucial in maintaining bone integrity by mitigating excessive bone resorption. This means that B lymphocytes play a very important role in regulating bone health, particularly in preventing excessive bone resorption.
In instances of osteoporosis pathology, notably under conditions of chronic inflammation commonly associated with postmenopausal osteoporosis in women, B lymphocytes can contribute to bone destruction (56). Unlike senile osteoporosis, which is characterized by reduced bone resorption and even more significant decreases in bone formation (57), postmenopausal osteoporosis sees a marked increase in osteoclast numbers and a surge in bone resorption (58). More critically, the abrupt decline in estrogen triggers an inflammatory state and elevates inflammatory factors within the bone milieu, which in turn stimulates B lymphocytes to produce more granulocyte macrophage colony-stimulating factor (GM-CSF) (59). This factor has been shown to promote proliferation in osteoclast progenitor cells and an uptick in osteoclast numbers (60). Additionally, under inflammatory conditions, activated B cells secrete RANKL, which is instrumental in the activation of osteoclast formation (61). Besides osteoporosis in postmenopausal women, smoking and chronic diseases such as obesity, diabetes, and hypertension can also lead the body to a state of chronic inflammation. Compared to the physiological state where OPG secretion inhibits osteoclasts, B cells in pathological states promote osteoclasts. However, it is undeniable that B cells appear to play a significant role in the regulation of osteoclastogenesis through the RANKL/OPG signaling pathway. This study has identified only three B lymphocyte subpopulations and the SSC-A on B cell marker that are protective against osteoporosis. Due to MR analysis avoiding confounding bias and reverse causation, the results have some significance, but whether these cell subpopulations have potential roles in conjunction with the RANKL/OPG signaling system remains to be further investigated.
The term “EM CD8br %CD8br (Maturation stages of T cell panel)” denotes a specific subset of T cells within the T cell differentiation process that are characterized by an effective memory function, high CD8 expression, and comprise a certain percentage of the overall T cell population. This subset of T cells, known as memory CD8+ T cells, is found persistently in the peripheral blood of patients who suffer from delayed bone healing. Memory CD8+ T cells have been observed to hinder the activity and osteogenic differentiation of bone marrow stromal cells by excessively producing interferon-gamma (IFN-γ) (62). This was also suggested as an osteoporosis risk factor in the MR results. CD39 is prevalently located on the surface of Treg cells and is part of the Nucleoside Triphosphate Diphosphohydrolase (NTPDase) superfamily, and it plays a significant role in T cell functionality, storage, and proliferation (63–65). The term “activated Treg+” refers to a subset of Treg cells that express both CD39 and CD25 and are in an activated state, exhibiting regulatory properties in response to antigen stimulation. Evidence found that CD39 mediates osteoporosis by regulating the balance between osteoclasts and osteoblasts through the Wnt/β-collagen pathway (66). Additional evidence corroborates that soluble pro-inflammatory cytokines IL-17A and IL-17F, secreted by T cells (Th17 cells), are vital mediators for the maturation and differentiation of osteoblasts, and play a crucial role in enhancing bone formation as well as facilitating the healing of fractures (67, 68). Despite the known association of various cytokines, proteins, and T cells with osteoporosis, the intricate processes related to bone metabolism and remodeling involving CD25 and CD39, which are expressed on the T cells, remain a mystery (69).
In our study, CD11c on CD62L+ myeloid DC, CCR2 on CD62L+ myeloid DC, and CCR2 on CD62L+ plasmacytoid DC were all observed to be protective factors against osteoporosis. CD62L, also known as L-selectin, is an adhesion molecule found primarily on the surface of lymphocytes and indicated by certain dendritic cells (70, 71). Down-regulation of CD62L antigen has been suggested as a possible mechanism for neutrophilia during Immune response (72, 73). A previous study reported decreased CD62L expression in the spleen and increased expression in the bone marrow during inflammation-induced osteoporosis (74). CCR2 serves as a chemokine receptor on the cell membrane, playing a critical role in the migration of monocytes and the management of inflammation-associated activities (75). It has been associated with a decline in bone mass and the facilitation of bone degradation (76). CCR2 can also be found on preosteoblasts and osteoblasts, where its interaction with the chemokine ligand 2 (CCL2) is thought to contribute to the fusion and development of osteoblasts (77–79). Research indicates that the removal or blocking of CCR2 results in a reduction in both the quantity and the functionality of fully formed osteoclasts, as well as mitigated bone resorption (80). According to MR findings from the current study, CCR2 activity on both plasmacytoid and conventional dendritic cells acts as a defensive mechanism against osteoporosis, aligning with the prevailing body of evidence.
MR results show that CD45 on CD33dim HLA DR+ CD11b-, CD45 on Mo MDSC, and CD33- HLA DR+ AC were protective factors against osteoporosis, and all three were associated with Myeloid cells. Where CD33dim indicates a lower level of CD33 expression and CD11b- indicates no expression of CD11b.Mo-MDSC stands for monocyte-derived myeloid-derived suppressor cells, which play an important role in the regulation of immune responses, inflammation, and the tumor microenvironment (81). CD45 is a common leukocyte co-antigen that plays an important role in immune responses and cell signaling (82). For T cells, the extracellular structural domain of CD45 is expressed in several different isoforms, and the expression of specific isoforms depends on specific cellular subpopulations. The nature of the ligands for the different isoforms of CD45 is still uncertain, and the exact mechanism by which the potential ligands regulate CD45 function is not known (83). CD45 is inextricably linked to the immune system, and current studies have elucidated significant associations between CD45 in systemic lupus erythematosus, rheumatoid arthritis, and HIV (84). Dysfunction of CD45 has also been associated with hematologic malignancies and Alzheimer’s disease (85). Some evidence suggests that perturbations in CD45 activity may lead to the development of autoimmune diseases (86, 87). Despite the abundant expression of CD45, its role in myeloid cells is unclear, and there is no definitive proof indicating an association with osteoporosis (88). HLA-DR is an MHC class II cell surface receptor encoded by the human leukocyte antigen complex on chromosome 6, region 6P21. A statistically significant association was found between the number of HLA-DR mismatches and the diagnosis of osteoporosis in renal transplant patients (89). Increased immune response due to HLA-DR mismatch may be associated with osteoporosis and hip fracture development (90). CD33 is a myeloid-specific antigen that is highly expressed in myeloid cells. CD33 plays an important role in immune regulation and cell adhesion functions (91). CD11b is an adhesion molecule expressed in monocytes and myeloid cells, which facilitates cell adhesion and migration (92, 93). The low expression of CD11b and CD33 in myeloid cells may limit their cell adhesion and migration functions, suggesting a negative causal effect on osteoporosis, which is enlightening.
Before this study, research (Cao’s study) had already reported causal relationships between immune cell phenotypes and bone mineral density (BMD). Therefore, it is important to discuss the necessity of this study. Cao’s study reported a total of 53 causal relationships between different immune cell phenotypes and BMD (see Cao’s study Supplementary Table 2) (46). Similar to this study, neither found causal effects of osteoporosis or BMD on immune cells. CD45 on CD33dim HLA DR+ CD11b- was the same exposure factor found in both studies, one study defined its outcome as total body BMD, while this study defined its outcome as osteoporosis. Using BMD as an outcome may lean more towards assessing the overall impact of immune cells on bone density. Using osteoporosis as an outcome aims to better understand the specific role of immune cells in the development of osteoporosis. Therefore, the causal relationship between the other 12 immune cell phenotypes and osteoporosis, which was found in this study, is unreported and holds clinical significance. In this study, we identified and extensively discussed the complex causal role of B cells in osteoporosis, whereas Cao’s study focused on the correlation between the CD40/CD40L system and bone metabolism. In addition, based on existing literature evidence, this study discusses in detail the potential links between each immune phenotype traits and osteoporosis The results of these two studies complement each other, providing insights from different perspectives into the potential causal relationships between osteoporosis and immune cells.
Through our study, we found that four immune cell characteristics, including absolute cell counts, median fluorescence intensities reflecting surface antigen levels, morphological parameters, and relative cells all had different causal relationships with osteoporosis. CD24, CD27, CD25 from B cells are protective factors against osteoporosis, and similar findings have been discussed, but more in-depth mechanisms of action remain to be investigated. The signal intensity of side scatter (SSC-A) on B cells reflects the protective effect of osteoporosis, which is related to the comprehensive assessment of B cell size, number, complexity, and granular content. Therefore, several active factors on B cells and a certain degree of B cells may be beneficial for osteoporosis. It must be mentioned that there are two risk factors associated with T cells. The EM CD8br %CD8br (Maturation stages of T cell panel) primarily represents effector T cells that have undergone stimulation and differentiation in immune response and possess the ability for rapid response. These cells are characterized by high expression of CD8. Another factor is a subset of regulatory T cells (Tregs), which are characterized by the surface expression of CD25 and CD39. Additionally, these Treg cells display an activated state. However, the regulation of T cells in osteoporosis is complex and diverse, and further evidence is still worth exploring.
4.1 Limitations
Due to the self-reported nature of the diagnostic data for osteoporosis in the GWAS database with the identifier ebi-a-GCST90038656, there is a potential for a certain degree of bias in the results. Due to database limitations, we were unable to stratify osteoporosis by gender, age, race, and whether or not they smoked, and we will work to address this issue in our upcoming studies.
5 Conclusions
In conclusion, we demonstrated a causal relationship between 13 immune phenotypes and osteoporosis through comprehensive bidirectional MR analysis, highlighting the complex pattern of interactions between the immune system and osteoporosis. In addition, our study significantly reduced the effects of unavoidable confounders, reverse causality, and other factors. This may provide new avenues for researchers to explore the biological mechanisms of osteoporosis and help to explore early intervention and therapeutic approaches.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
DZ: Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. CZ: Writing – review & editing. GG: Writing – review & editing. DZ: Writing – review & editing. QC: Data curation, Methodology, Software, Supervision, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Hospital Natural Science Foundation of Sichuan Taikang Hospital 2023, Project number: SCTK2023ZR02.
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.2024.1341002/full#supplementary-material
Supplementary File 1 | Results of the causal effect of immune cells on osteoporosis.
Supplementary File 2 | Sensitivity analysis results of immune cells on osteoporosis.
Supplementary File 3 | Scatterplots of immune cells on osteoporosis.
Supplementary File 4 | Pleiotropy and Heterogeneity of immune cells on osteoporosis.
Supplementary File 5 | Supplementary figures.
Supplementary File 6 | Causal effects of osteoporosis on immune cells.
Supplementary File 7 | Removal of weak IVs and F-statistics.
References
1. Ferbebouh M, Vallières F, Benderdour M, Fernandes J. The pathophysiology of immunoporosis: innovative therapeutic targets. Inflammation Res. (2021) 70:859–75. doi: 10.1007/s00011-021-01484-9
2. Srivastava RK, Dar HY, Mishra PK. Immunoporosis: immunology of osteoporosis-role of T cells. Front Immunol. (2018) 9:657. doi: 10.3389/fimmu.2018.00657
3. Lorenzo J, Horowitz M, Choi Y. Osteoimmunology: interactions of the bone and immune system. Endocr Rev. (2008) 29:403–40. doi: 10.1210/er.2007-0038
4. Walsh MC, Takegahara N, Kim H, Choi Y. Updating osteoimmunology: regulation of bone cells by innate and adaptive immunity. Nat Rev Rheumatol. (2018) 14:146–56. doi: 10.1038/nrrheum.2017.213
5. Li JY, D’Amelio P, Robinson J, Walker LD, Vaccaro C, Luo T, et al. IL-17A is increased in humans with primary hyperparathyroidism and mediates PTH-induced bone loss in mice. Cell Metab. (2015) 22:799–810. doi: 10.1016/j.cmet.2015.09.012
6. Pietschmann P, Mechtcheriakova D, Meshcheryakova A, Föger-Samwald U, Ellinger I. Immunology of osteoporosis: A mini-review. Gerontology. (2016) 62:128–37. doi: 10.1159/000431091
7. Du D, Zhou Z, Zhu L, Hu X, Lu J, Shi C, et al. TNF-α suppresses osteogenic differentiation of MSCs by accelerating P2Y(2) receptor in estrogen-deficiency induced osteoporosis. Bone. (2018) 117:161–70. doi: 10.1016/j.bone.2018.09.012
8. Loi F, Córdova LA, Zhang R, Pajarinen J, Lin TH, Goodman SB, et al. The effects of immunomodulation by macrophage subsets on osteogenesis. vitro. Stem Cell Res Ther. (2016) 7:15. doi: 10.1186/s13287-016-0276-5
9. Li CJ, Xiao Y, Sun YC, He WZ, Liu L, Huang M, et al. Senescent immune cells release grancalcin to promote skeletal aging. Cell Metab. (2021) 33:1957–1973.e1956. doi: 10.1016/j.cmet.2021.08.009
10. Abe Y, Kawakami A, Nakashima T, Ejima E, Fujiyama K, Kiriyama T, et al. Etidronate inhibits human osteoblast apoptosis by inhibition of pro-apoptotic factor(s) produced by activated T cells. J Lab Clin Med. (2000) 136:344–54. doi: 10.1067/mlc.2000.109757
11. Sapra L, Bhardwaj A, Mishra PK, Garg B, Verma B, Mishra GC, et al. Regulatory B cells (Bregs) inhibit osteoclastogenesis and play a potential role in ameliorating ovariectomy-induced bone loss. Front Immunol. (2021) 12:691081. doi: 10.3389/fimmu.2021.691081
12. Sponholtz TR, Zhang X, Fontes JD, Meigs JB, Cupples LA, Kiel DP, et al. Association between inflammatory biomarkers and bone mineral density in a community-based cohort of men and women. Arthritis Care Res (Hoboken). (2014) 66:1233–40. doi: 10.1002/acr.22270
13. Berglundh S, Malmgren L, Luthman H, McGuigan F, Åkesson K. C-reactive protein, bone loss, fracture, and mortality in elderly women: a longitudinal study in the OPRA cohort. Osteoporos Int. (2015) 26:727–35. doi: 10.1007/s00198-014-2951-7
14. Kania DM, Binkley N, Checovich M, Havighurst T, Schilling M, Ershler WB. Elevated plasma levels of interleukin-6 in postmenopausal women do not correlate with bone density. J Am Geriatr Soc. (1995) 43:236–9. doi: 10.1111/j.1532-5415.1995.tb07328.x
15. Cauley JA, Barbour KE, Harrison SL, Cloonan YK, Danielson ME, Ensrud KE, et al. Inflammatory markers and the risk of hip and vertebral fractures in men: the osteoporotic fractures in men (MrOS). J Bone Miner Res. (2016) 31:2129–38. doi: 10.1002/jbmr.2905
16. Li C, Liu J, Lin J, Shang H. COVID-19 and risk of neurodegenerative disorders: A Mendelian randomization study. Transl Psychiatry. (2022) 12:283. doi: 10.1038/s41398-022-02052-3
17. Grover S, Sharma M. Sleep, pain, and neurodegeneration: A mendelian randomization study. Front Neurol. (2022) 13:765321. doi: 10.3389/fneur.2022.765321
18. Timpson NJ, Wade KH, Smith GD. Mendelian randomization: application to cardiovascular disease. Curr Hypertens Rep. (2012) 14:29–37. doi: 10.1007/s11906-011-0242-7
19. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. (2014) 23:R89–98. doi: 10.1093/hmg/ddu328
20. Orrù V, Steri M, Sidore C, Marongiu M, Serra V, Olla S, et al. Complex genetic signatures in immune cells underlie autoimmunity and inform therapy. Nat Genet. (2020) 52:1036–45. doi: 10.1038/s41588-020-0684-4
21. Wang C, Zhu D, Zhang D, Zuo X, Yao L, Liu T, et al. Causal role of immune cells in schizophrenia: Mendelian randomization (MR) study. BMC Psychiatry. (2023) 23:590. doi: 10.1186/s12888-023-05081-4
22. Yu XH, Yang YQ, Cao RR, Bo L, Lei SF. The causal role of gut microbiota in development of osteoarthritis. Osteoarthritis Cartilage. (2021) 29:1741–50. doi: 10.1016/j.joca.2021.08.003
23. Hautakangas H, Winsvold BS, Ruotsalainen SE, Bjornsdottir G, Harder AVE, Kogelman LJA, et al. Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles. Nat Genet. (2022) 54:152–60. doi: 10.1038/s41588-021-00990-0
24. Jayaraj RL, Azimullah S, Beiram R. Diabetes as a risk factor for Alzheimer’s disease in the Middle East and its shared pathological mediators. Saudi J Biol Sci. (2020) 27:736–50. doi: 10.1016/j.sjbs.2019.12.028
25. Chalitsios CV, Meena D, Manou M, Papagiannopoulos C, Markozannes G, Gill D, et al. Multiple long-term conditions in people with psoriasis: alatent class and bidirectional Mendelian randomization analysis. Br J Dermatol. (2023) 190(3):364–73. doi: 10.1093/bjd/ljad410
26. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. (2016) 40:304–14. doi: 10.1002/gepi.21965
27. Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum Mol Genet. (2018) 27:R195–r208. doi: 10.1093/hmg/ddy163
28. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. (2015) 44:512–25. doi: 10.1093/ije/dyv080
29. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. (2018) 50:693–8. doi: 10.1038/s41588-018-0099-7
30. Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
31. Nimptsch K, Aleksandrova K, Pham TT, Papadimitriou N, Janke J, Christakoudi S, et al. Prospective and Mendelian randomization analyses on the association of circulating fatty acid binding protein 4 (FABP-4) and risk of colorectal cancer. BMC Med. (2023) 21:391. doi: 10.1186/s12916-023-03104-1
32. Gutzeit C, Chen K, Cerutti A. The enigmatic function of IgD: some answers at last. Eur J Immunol. (2018) 48:1101–13. doi: 10.1002/eji.201646547
33. van de Peppel J, Schaaf GJ, Matos AA, Guo Y, Strini T, Verschoor W, et al. Cell surface glycoprotein CD24 marks bone marrow-derived human mesenchymal stem/stromal cells with reduced proliferative and differentiation capacity. In Vitro. Stem Cells Dev. (2021) 30:325–36. doi: 10.1089/scd.2021.0027
34. Liu Y, Zheng P. CD24-Siglec interactions in inflammatory diseases. Front Immunol. (2023) 14:1174789. doi: 10.3389/fimmu.2023.1174789
35. Altevogt P, Sammar M, Hüser L, Kristiansen G. Novel insights into the function of CD24: A driving force in cancer. Int J Cancer. (2021) 148:546–59. doi: 10.1002/ijc.33249
37. Zikherman J, Parameswaran R, Weiss A. Endogenous antigen tunes the responsiveness of naive B cells but not T cells. Nature. (2012) 489:160–4. doi: 10.1038/nature11311
38. Quách TD, Manjarrez-Orduño N, Adlowitz DG, Silver L, Yang H, Wei C, et al. Anergic responses characterize a large fraction of human autoreactive naive B cells expressing low levels of surface IgM. J Immunol. (2011) 186:4640–8. doi: 10.4049/jimmunol.1001946
39. Duty JA, Szodoray P, Zheng NY, Koelsch KA, Zhang Q, Swiatkowski M, et al. Functional anergy in a subpopulation of naive B cells from healthy humans that express autoreactive immunoglobulin receptors. J Exp Med. (2009) 206:139–51. doi: 10.1084/jem.20080611
40. Goodnow CC, Crosbie J, Jorgensen H, Brink RA, Basten A. Induction of self-tolerance in mature peripheral B lymphocytes. Nature. (1989) 342:385–91. doi: 10.1038/342385a0
41. Barkeer S, Chugh S, Batra SK, Ponnusamy MP. Glycosylation of cancer stem cells: function in stemness, tumorigenesis, and metastasis. Neoplasia. (2018) 20:813–25. doi: 10.1016/j.neo.2018.06.001
42. Ohl C, Albach C, Altevogt P, Schmitz B. N-glycosylation patterns of HSA/CD24 from different cell lines and brain homogenates: a comparison. Biochimie. (2003) 85:565–73. doi: 10.1016/S0300-9084(03)00107-X
43. Schäck LM, Buettner M, Wirth A, Neunaber C, Krettek C, Hoffmann A, et al. Expression of CD24 in human bone marrow-derived mesenchymal stromal cells is regulated by TGFβ3 and induces a myofibroblast-like genotype. Stem Cells Int. (2016) 2016:1319578. doi: 10.1155/2016/1319578
44. Borst J, Hendriks J, Xiao Y. CD27 and CD70 in T cell and B cell activation. Curr Opin Immunol. (2005) 17:275–81. doi: 10.1016/j.coi.2005.04.004
45. Beckers L, Somers V, Fraussen J. IgD(-)CD27(-) double negative (DN) B cells: Origins and functions in health and disease. Immunol Lett. (2023) 255:67–76. doi: 10.1016/j.imlet.2023.03.003
46. Cao RR, Yu XH, Xiong MF, Li XT, Deng FY, Lei SF. The immune factors have complex causal regulation effects on bone mineral density. Front Immunol. (2022) 13:959417. doi: 10.3389/fimmu.2022.959417
47. Luo H, Zhu Y, Guo B, Ruan Z, Liu Z, Fan Z, et al. Causal relationships between CD25 on immune cells and hip osteoarthritis. Front Immunol. (2023) 14:1247710. doi: 10.3389/fimmu.2023.1247710
48. Yu Q, Xu M, Yu F, Jin Y. CD4(+)CD25(+) regulatory T cells as a therapeutic target in rheumatoid arthritis. Cent Eur J Immunol. (2014) 39:100–3. doi: 10.5114/ceji.2014.42133
49. Frase D, Lee C, Nachiappan C, Gupta R, Akkouch A. The inflammatory contribution of B-lymphocytes and neutrophils in progression to osteoporosis. Cells. (2023) 12:1744. doi: 10.3390/cells12131744
50. Hendriks RW, Middendorp S. The pre-BCR checkpoint as a cell-autonomous proliferation switch. Trends Immunol. (2004) 25:249–56. doi: 10.1016/j.it.2004.02.011
51. Aoki Y, Isselbacher KJ, Pillai S. Bruton tyrosine kinase is tyrosine phosphorylated and activated in pre-B lymphocytes and receptor-ligated B cells. Proc Natl Acad Sci U.S.A. (1994) 91:10606–9. doi: 10.1073/pnas.91.22.10606
52. Shinohara M, Koga T, Okamoto K, Sakaguchi S, Arai K, Yasuda H, et al. Tyrosine kinases Btk and Tec regulate osteoclast differentiation by linking RANK and ITAM signals. Cell. (2008) 132:794–806. doi: 10.1016/j.cell.2007.12.037
53. Boyce BF, Xing L. Functions of RANKL/RANK/OPG in bone modeling and remodeling. Arch Biochem Biophys. (2008) 473:139–46. doi: 10.1016/j.abb.2008.03.018
54. Yun TJ, Chaudhary PM, Shu GL, Frazer JK, Ewings MK, Schwartz SM, et al. OPG/FDCR-1, a TNF receptor family member, is expressed in lymphoid cells and is up-regulated by ligating CD40. J Immunol. (1998) 161:6113–21. doi: 10.4049/jimmunol.161.11.6113
55. Walsh MC, Choi Y. Biology of the RANKL-RANK-OPG system in immunity, bone, and beyond. Front Immunol. (2014) 5:511. doi: 10.3389/fimmu.2014.00511
56. Fischer V, Haffner-Luntzer M. Interaction between bone and immune cells: Implications for postmenopausal osteoporosis. Semin Cell Dev Biol. (2022) 123:14–21. doi: 10.1016/j.semcdb.2021.05.014
58. Brown JP. Long-term treatment of postmenopausal osteoporosis. Endocrinol Metab (Seoul). (2021) 36:544–52. doi: 10.3803/EnM.2021.301
59. Toni R, Di Conza G, Barbaro F, Zini N, Consolini E, Dallatana D, et al. Microtopography of immune cells in osteoporosis and bone lesions by endocrine disruptors. Front Immunol. (2020) 11:1737. doi: 10.3389/fimmu.2020.01737
60. Zhang Z, Yuan W, Deng J, Wang D, Zhang T, Peng L, et al. Granulocyte colony stimulating factor (G-CSF) regulates neutrophils infiltration and periodontal tissue destruction in an experimental periodontitis. Mol Immunol. (2020) 117:110–21. doi: 10.1016/j.molimm.2019.11.003
61. Li Y, Terauchi M, Vikulina T, Roser-Page S, Weitzmann MN. B cell production of both OPG and RANKL is significantly increased in aged mice. Open Bone J. (2014) 6:8–17. doi: 10.2174/1876525401406010008
62. Reinke S, Geissler S, Taylor WR, Schmidt-Bleek K, Juelke K, Schwachmeyer V, et al. Terminally differentiated CD8+ T cells negatively affect bone regeneration in humans. Sci Transl Med. (2013) 5:177ra136. doi: 10.1126/scitranslmed.3004754
63. Timperi E, Barnaba V. CD39 regulation and functions in T cells. Int J Mol Sci. (2021) 22:8068. doi: 10.3390/ijms22158068
64. Takenaka MC, Robson S, Quintana FJ. Regulation of the T cell response by CD39. Trends Immunol. (2016) 37:427–39. doi: 10.1016/j.it.2016.04.009
65. Bono MR, Fernández D, Flores-Santibáñez F, Rosemblatt M, Sauma D. CD73 and CD39 ectonucleotidases in T cell differentiation: Beyond immunosuppression. FEBS Lett. (2015) 589:3454–60. doi: 10.1016/j.febslet.2015.07.027
66. Wu W, Xiao Z, Chen Y, Deng Y, Zeng D, Liu Y, et al. CD39 produced from human GMSCs regulates the balance of osteoclasts and osteoblasts through the Wnt/β-catenin pathway in osteoporosis. Mol Ther. (2020) 28:1518–32. doi: 10.1016/j.ymthe.2020.04.003
67. Nam D, Mau E, Wang Y, Wright D, Silkstone D, Whetstone H, et al. T-lymphocytes enable osteoblast maturation via IL-17F during the early phase of fracture repair. PloS One. (2012) 7:e40044. doi: 10.1371/journal.pone.0040044
68. Croes M, Öner FC, van Neerven D, Sabir E, Kruyt MC, Blokhuis TJ, et al. Proinflammatory T cells and IL-17 stimulate osteoblast differentiation. Bone. (2016) 84:262–70. doi: 10.1016/j.bone.2016.01.010
69. Dar HY, Azam Z, Anupam R, Mondal RK, Srivastava RK. Osteoimmunology: The Nexus between bone and immune system. Front Biosci. (2018) 23:464–92. doi: 10.2741/4600
71. Savage ND, Harris SH, Rossi AG, De Silva B, Howie SE, Layton GT, et al. Inhibition of TCR-mediated shedding of L-selectin (CD62L) on human and mouse CD4+ T cells by metalloproteinase inhibition: analysis of the regulation of Th1/Th2 function. Eur J Immunol. (2002) 32:2905–14. doi: 10.1002/1521-4141(2002010)32:10<2905::AID-IMMU2905>3.0.CO;2-6
72. Lobach AR, Uetrecht J. Clozapine promotes the proliferation of granulocyte progenitors in the bone marrow leading to increased granulopoiesis and neutrophilia in rats. Chem Res Toxicol. (2014) 27:1109–19. doi: 10.1021/tx500184c
73. Zhang D, Chen G, Manwani D, Mortha A, Xu C, Faith JJ, et al. Neutrophil ageing is regulated by the microbiome. Nature. (2015) 525:528–32. doi: 10.1038/nature15367
74. Li L, Wang Y, Zhang N, Zhang Y, Lin J, Qiu X, et al. Heterozygous deletion of LRP5 gene in mice alters profile of immune cells and modulates differentiation of osteoblasts. Biosci Trends. (2018) 12:266–74. doi: 10.5582/bst.2018.01013
75. Mader TL, Novotny SA, Lin AS, Guldberg RE, Lowe DA, Warren GL. CCR2 elimination in mice results in larger and stronger tibial bones but bone loss is not attenuated following ovariectomy or muscle denervation. Calcif Tissue Int. (2014) 95:457–66. doi: 10.1007/s00223-014-9914-z
76. Binder NB, Niederreiter B, Hoffmann O, Stange R, Pap T, Stulnig TM, et al. Estrogen-dependent and C-C chemokine receptor-2-dependent pathways determine osteoclast behavior in osteoporosis. Nat Med. (2009) 15:417–24. doi: 10.1038/nm.1945
77. Eraltan H, Cacina C, Kahraman OT, Kurt O, Aydogan HY, Uyar M, et al. MCP-1 and CCR2 gene variants and the risk for osteoporosis and osteopenia. Genet Test Mol Biomarkers. (2012) 16:229–33. doi: 10.1089/gtmb.2011.0216
78. Xing Z, Lu C, Hu D, Yu YY, Wang X, Colnot C, et al. Multiple roles for CCR2 during fracture healing. Dis Model Mech. (2010) 3:451–8. doi: 10.1242/dmm.003186
79. Sul OJ, Ke K, Kim WK, Kim SH, Lee SC, Kim HJ, et al. Absence of MCP-1 leads to elevated bone mass via impaired actin ring formation. J Cell Physiol. (2012) 227:1619–27. doi: 10.1002/jcp.22879
80. Barros SP, Arce RM, Galloway P, Lawter R, Offenbacher S. Therapeutic effect of a topical CCR2 antagonist on induced alveolar bone loss in mice. J Periodontal Res. (2011) 46:246–51. doi: 10.1111/jre.2011.46.issue-2
81. Ren X, Xiao J, Zhang W, Wang F, Yan Y, Wu X, et al. Inhibition of CCL7 derived from Mo-MDSCs prevents metastatic progression from latency in colorectal cancer. Cell Death Dis. (2021) 12:484. doi: 10.1038/s41419-021-03698-5
82. Rheinländer A, Schraven B, Bommhardt U. CD45 in human physiology and clinical medicine. Immunol Lett. (2018) 196:22–32. doi: 10.1016/j.imlet.2018.01.009
83. Altin JG, Sloan EK. The role of CD45 and CD45-associated molecules in T cell activation. Immunol Cell Biol. (1997) 75:430–45. doi: 10.1038/icb.1997.68
84. Penninger JM, Irie-Sasaki J, Sasaki T, Oliveira-dos-Santos AJ. CD45: new jobs for an old acquaintance. Nat Immunol. (2001) 2:389–96. doi: 10.1038/87687
85. Tan J, Town T, Mori T, Wu Y, Saxe M, Crawford F, et al. CD45 opposes beta-amyloid peptide-induced microglial activation via inhibition of p44/42 mitogen-activated protein kinase. J Neurosci. (2000) 20:7587–94. doi: 10.1523/JNEUROSCI.20-20-07587.2000
86. Takeuchi T, Pang M, Amano K, Koide J, Abe T. Reduced protein tyrosine phosphatase (PTPase) activity of CD45 on peripheral blood lymphocytes in patients with systemic lupus erythematosus (SLE). Clin Exp Immunol. (1997) 109:20–6. doi: 10.1046/j.1365-2249.1997.4371334.x
87. Blasini AM, Alonzo E, Chacón R, Riera R, Stekman IL, Rodriguez MA. Abnormal pattern of tyrosine phosphorylation in unstimulated peripheral blood T lymphocytes from patients with systemic lupus erythematosus. Lupus. (1998) 7:515–23. doi: 10.1191/096120398678920604
88. Hermiston ML, Xu Z, Weiss A. CD45: a critical regulator of signaling thresholds in immune cells. Annu Rev Immunol. (2003) 21:107–37. doi: 10.1146/annurev.immunol.21.120601.140946
89. Opelz G, Döhler B. Association of mismatches for HLA-DR with incidence of posttransplant hip fracture in kidney transplant recipients. Transplantation. (2011) 91:65–9. doi: 10.1097/TP.0b013e3181fa94d6
90. Nikkel LE, Hollenbeak CS, Fox EJ, Uemura T, Ghahramani N. Risk of fractures after renal transplantation in the United States. Transplantation. (2009) 87:1846–51. doi: 10.1097/TP.0b013e3181a6bbda
91. Maakaron JE, Rogosheske J, Long M, Bachanova V, Mims AS. CD33-targeted therapies: beating the disease or beaten to death? J Clin Pharmacol. (2021) 61:7–17. doi: 10.1002/jcph.1730
92. Schmid MC, Khan SQ, Kaneda MM, Pathria P, Shepard R, Louis TL, et al. Integrin CD11b activation drives anti-tumor innate immunity. Nat Commun. (2018) 9:5379. doi: 10.1038/s41467-018-07387-4
Keywords: osteoporosis, immunocyte phenotype, mendelian randomization, causal connection, B cells
Citation: Zhou D, Zi C, Gan G, Tang S and Chen Q (2024) An exploration of the causal relationship between 731 immunophenotypes and osteoporosis: a bidirectional Mendelian randomized study. Front. Endocrinol. 15:1341002. doi: 10.3389/fendo.2024.1341002
Received: 19 November 2023; Accepted: 03 July 2024;
Published: 17 July 2024.
Edited by:
Xingming Shi, Augusta University, United StatesReviewed by:
Rajeev Aurora, Saint Louis University, United StatesXingyu He, University of Cincinnati, United States
Copyright © 2024 Zhou, Zi, Gan, Tang and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiu Chen, chenqiu1005@cdutcm.edu.cn