Introduction
Multiple myeloma (MM) is a hematological malignancy resulting from uncontrolled proliferation of plasma cells in the bone marrow. Myeloma bone disease (MBD) is the most common complication of MM, with approximately 80%–90% of patients showing bone lesions at diagnosis 1 . MBD mainly manifests as hypercalcemia, severe osteoporosis, severe bone pain, osteolytic lesions, pathological fractures, and vertebral compression fractures, which have negative effects on both the quality of life and survival of patients 2 . The pathogenesis of MBD is mainly a result of the formation of osteoclasts activated by MM cells or stromal cells, which leads to the enhancement of bone resorption function, accompanied by the inhibition of osteoblast function, which eventually leads to an imbalanced bone metabolism 1 , 3 . Extensive research on the function of osteoclasts has shown that the mechanism of osteoclast activation is complex, and plays a critical role in the pathogenesis of MBD 4 .
The complement system has an important relationship with bone development and homeostasis 5 . Complement proteins, including C3, factor B, C5, and C9, are found in the growth plate during bone development 6 . C1s are expressed by hypertrophic chondrocytes in the ossification center of the femur, but not by normal articular cartilage cells 7 , 8 . Ignatius et al. 9 found that the C3 and C5 complement components, the C3aR and C5aR anaphylatoxin receptors, the MCP (CD46), DAF (CD55) complement regulators, and MACIF (CD59) were expressed in differentiated and undifferentiated mesenchymal stem cells (MSCs) and osteoblasts. Moreover, osteoclasts express MACIF (CD59), C3aR, C5aR, and C3 9 , and the formation of osteoclasts is regulated by C3 (C3a), C5a, C3aR, and C5aR 10 – 12 . In addition, Teo et al. 13 found that C1q was produced by osteoclasts, and was responsible for enhancing osteoclast development.
In the present study, we analyzed correlations of serum levels of C3/C4/C4a in both newly diagnosed multiple myeloma (NDMM) patients and healthy controls, and bone disease stages, with the number of osteoblast precursors, osteoclast precursors (OCPs), and biological indicators of MBD. In a previous study, we reported the effects of C3a/C4a on the occurrence of MBD 14 , and found that serum levels of C3, C4, and C4a in NDMM patients were significantly positively correlated with the severity of bone disease, the number of OCPs, and the levels of the C-terminal cross-linking telopeptide biomarker of type I collagen/tartrate resistant acid phosphatase isoform-5b, which is related to bone destruction. We therefore proposed that serum levels of C3 and C4 may be sensitive biomarkers of myeloma bone disease, so we characterized the effects and potential mechanisms of C3a/C4a on osteoclasts in patients with MM.
Materials and methods
Patients and samples
The study cohort included 124 NDMM patients. All patients were in-patients at the Hematology Department of Tianjin Medical University General Hospital, from May 2017 to May 2019. Patients were diagnosed according to International Myeloma Working Group updated criteria for MM 15 , 16 . The clinical characteristics of the patients are shown in Table 1 . This study was approved by the Ethical Committee of Tianjin Medical University (Approval No. IRB2020-WZ-018). Samples (10 mL) of bone marrow were obtained from all NDMM patients after obtaining written informed consent from the participants.
Table 1
Baseline characteristics of patients
R-ISS: Revised International Staging System; Hb: hemoglobin.
| Characteristics | Newly diagnosed multiple myeloma patients ( n = 124) (%) |
| Gender | |
| Male | 78 (62.903) |
| Female | 46 (37.10) |
| Age (median) | 64 |
| Range | 34˜87 |
| R-ISS stage | |
| I | 15 (12.10) |
| II | 26 (20.97) |
| III | 83 (66.93) |
| M-component type | |
| IgG | 62 (50.00) |
| IgA | 25 (20.16) |
| IgM | 1 (0.8) |
| Light chain only | 26 (20.97) |
| Nonsecretory | 10 (8.07) |
| 2 -microglobulin | |
| < 5.5 mg/L | 41 (33.06) |
| ≥ 5.5 mg/L | 83 (66.94) |
| Ca | |
| > 2.75 mmol/L | 30 (24.19) |
| Serum creatinine | |
| < 177 μmol/L | 80 (64.52) |
| ≥ 177 μmol/L | 44 (35.48) |
| Hb | |
| < 100 g/L | 85 (68.55) |
| ≥ 100 g/L | 39 (31.45) |
Osteoclast cultures
Bone marrow mononuclear cells (BMMCs) were isolated from the bone marrow of patients using Ficoll-Paque Plus solution (Amersham Biosciences, Chalfont St. Giles, UK). Isolated BMMCs were incubated in α-MEM medium supplemented with 10% fetal bovine serum (FBS), 100 U/mL penicillin, and 100 mg/mL streptomycin in the presence of recombinant human receptor activator for nuclear factor-κ B ligand (RANKL) (150 ng/mL; Miltenyi Biotec, Sunnyvale, CA, USA) and recombinant human macrophage colony-stimulating factor (M-CSF) (50 ng/mL; Miltenyi Biotec, USA), followed by seeding in a 24-well plate at 1 × 10 6 cells/well, as described previously 17 . The cells of each patient were divided into 3 wells with different concentrations of recombinant human complement C3a (Emd Millipore, Burlington, MA, USA) at 10 μg/mL, 1 μg/mL, or 0 μg/mL DMSO), according to studies on the effect of human C3a on osteoblasts, osteoclasts, and MSCs 9 , 18 . The complement C4a experiment was similar, with a concentration of C4a (Fitzgerald, Crossville, TN, USA) of 10 μg/mL, 1 μg/mL, or 0 μg/mL; the concentration was set according to the serum level of C4a measured by an ELISA in NDMM patients. To study the effect of the SGK inhibitor (EMD638683) on osteoclasts, the concentration of EMD638683 (MedChemExpress, Monmouth Junction, NJ, USA) was set to 10 μM, 1 μM, or 0 μM, according to previous studies on EMD638683 19 , 20 . The 24-well plate was cultured for 14 days at 37 °C in a humidified incubator containing 5% CO 2 . Culture media with cytokines, recombinant human complement proteins, or EMD638683 were replaced every 3 days.
TRAP staining
Osteoclasts were identified using a TRAP staining kit (Sigma-Aldrich, St. Louis, MO, USA), according to manufacturer’s instructions. Briefly, 0.5 mL of Fast Garnet GBC base solution and 0.5 mL of sodium nitrite solution were added to a tube and mixed by gentle inversion for 30 s. Then, 100 μL of the resulting mixture was mixed with 4.5 mL of deionized water, 50 μL of naphthol AS-BI phosphate solution, 200 μL of acetate solution, and 100 μL of tartrate solution to a final volume of 10 mL in a tube. The 24-well plates with osteoclasts were fixed by immersing in fixative solution for 30 s. Finally, 500 μL of the mixed solution was added to each well and incubated 10 min in a water bath at 37 °C in the dark. After 10 min, the plates were rinsed thoroughly in deionized water and evaluated microscopically. Purple multinucleated cells (≥ 3 nuclei) were defined as osteoclasts.
Quantitative real-time PCR
Quantitative real-time PCR was performed as described previously 21 . RNA was extracted from osteoclasts using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and mRNA expression was quantified using the Bio-Rad iQ5 Real-time system (Bio-Rad, Hercules, CA, USA). The primer sequences of osteoclast-associated genes RANKL/OSCAR/TRAP/Cathepsin K and related genes of the C3a activating the osteoclasts pathway are shown in Supplementary Table S1 and include: SGK3/PIK3CA/POSTN/COL1A1/COL1A2/CREB1/MDM2/IKBKG/NTRK2/CDK2/TCL1A/PDK1/GAPDH.
SYBR Green (Invitrogen) was used as a double-strand DNA-specific dye. The amplification of RANKL/OSCAR/TRAP/Cathepsin K/SGK3/PIK3CA/POSTN/COL1A1/COL1A2/CREB1/MDM2/IKBKG/NTRK2/CDK2/TCL1A/PDK1 was conducted over 45 cycles at 95 °C for 30 s and 95 °C for 5 s, followed by extension at the temperatures shown in Supplementary Table S2 . Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as the housekeeping gene to standardize the targeted mRNA expression. The levels of RANKL/OSCAR/TRAP/Cathepsin K/SGK3/PIK3CA/POSTN/COL1A1/COL1A2/CREB1/MDM2/IKBKG/NTRK2/CDK2/TCL1A/PDK1 were calculated using the 2 -ΔΔCt method [(Ct, target gene Ct, GAPDH) sample - (Ct, target gene Ct, GAPDH) control] after normalizing the data to GAPDH mRNA expression.
Resorption pit assay
Isolated BMMCs were incubated in α-MEM medium supplemented with 10% FBS, 100 U/mL penicillin, and 100 mg/mL streptomycin in the presence of recombinant human RANKL (150 ng/mL), recombinant human M-CSF (50 ng/mL), recombinant human complement C3a/C4a, or EMD638683, followed by seeding on bovine cortical bone slices in a 24-well plate at 1 × 10 6 cells/well for 14 days at 37 °C in a humidified incubator containing 5% CO 2 . Culture media with cytokines, recombinant human complement proteins, or EMD638683 were replaced every 3 days. After culturing, the bovine cortical bone slices were fixed with 2.5% glutaraldehyde for 7 min, washed 3 times by sonication in 0.25 M aqua ammonia for 5 min, and subjected to ethanol gradient dehydration, drying, and spraying before observation using scanning electron microscopy (SEM).
RNA-seq analysis
Total RNA was extracted from the cell samples for use in RNA-seq analysis. The RNA samples were converted into individual cDNA libraries and generated using the BGISEQ-500 platform of the Beijing Genomics Institute (BGI) (Beijing, China) 22 , 23 . All of the generated raw sequencing reads were filtered to obtain clean reads before storing in the FASTQ format 24 . Clean reads were aligned to reference genes and genomes using Bowtie2 and HISAT, respectively 25 , 26 . The expression levels of the genes were normalized to FPKM using RSEM 27 . Genes were selected using the criteria of differentially expressed genes (DEGs) between the C3a group and control with a fold-change ≥ 1.5 and adjusted P- value ≤ 0.05. The relative enrichment and pathway annotation of genes for various functional associations were determined using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database.
Western blot
Western blot analysis was performed as described previously 28 to evaluate the PI3K/PDK1/SGK3/p-SGK3/AKT/p-AKT protein concentrations in the osteoclasts supplemented with recombinant human complement C3a or EMD638683. The cells were lysed on ice in a lysis buffer containing a phenylmethylsulfonyl fluoride protease inhibitor (Sigma-Aldrich), phosphatase inhibitor (Cell Signaling Technology, Danvers, MA, USA), and RIPA lysis buffer (Sigma-Aldrich) configured in a proportion of 1:1:100. The protein concentration was determined using a bicinchoninic acid (BCA) protein assay kit (Thermo Fisher Scientific, Waltham, MA, USA). The extracted proteins (30 μg/lane) were separated using SDS/PAGE (12% gels). The proteins were subsequently transferred to nitrocellulose membranes (Bio-Rad Laboratories, Hercules, CA, USA), which were blocked using 5% skim milk by incubating for 1 h at room temperature. The membranes were then incubated with primary antibodies against PI3K/PDK1/SGK3/p-SGK3/AKT/p-AKT/GAPDH [1:1000 dilution in 5% bovine serum albumin (BSA); Cell Signaling Technology] overnight at 4 °C. The membranes were washed with Tris-buffered saline containing 0.1% Tween-20 and incubated with horseradish peroxidase-conjugated anti-rabbit IgG sheep antibody (1:5000 dilution in 5% BSA; Abcam, Cambridge, MA, USA) for 1 h at room temperature. The protein bands were visualized using the hypersensitive ECL chemiluminescence kit (Cell Signaling Technology).
Statistical analysis
Data are expressed as numerical variables. The osteoclast area, absorption area of osteoclast resorption pit per view, and the relative gray scale of related protein expressions are expressed as the mean ± standard deviation (SD). The unpaired t -test and one-way analysis of variance were used to analyze the significance between different groups. The expressions of target genes are expressed as the median. Comparison between data were conducted using a nonparametric test. SPSS statistical software for Windows, version 21.0 (SPSS, Chicago, IL, USA) was used to conduct the statistical analysis. P < 0.05 was considered statistically significant.
Results
Complement C3a significantly promoted the formation and function of osteoclasts, while complement C4a did not
To evaluate the effect of C3a/C4a on osteoclasts in NDMM patients, we observed the formation and function of osteoclasts in different concentrations of C3a and C4a (1 μg/mL and 10 μg/mL).
In vitro , the osteoclasts area per view from 30 patients induced by 1 μg/mL (mean ± SD: 50.828 ± 12.984%) and 10 μg/mL (53.663 ± 12.685%) of C3a was significantly increased compared to the control group (0 μg/mL) (34.635 ± 8.916%) ( P < 0.001; P < 0.001) ( Figure 1 A and 1B ). There was no difference between 1 μg/mL (mean ± SD: 34.942 ± 9.920%) and 10 μg/mL (37.034 ± 8.964%) of the C4a and control groups (33.635 ± 6.639%) in 15 patients ( Figure 1 A and 1C ).

The relative mRNA expressions of the OSCAR/TRAP/RANKL/Cathepsin K genes from 30 patients were measured. The expressions of these genes induced by 1 μg/mL (median: 5.041, 3.726, 1.638, and 4.752, respectively) and 10 μg/mL (median: 5.140, 3.702, 2.250, and 5.172, respectively) in the C3a group was significantly increased compared to the control group (median: 3.137, 2.004, 0.573, and 2.257, respectively) (1 μg/mL: P = 0.001, P = 0.003, P < 0.001, and P = 0.008, respectively; 10 μg/mL: P < 0.001, P = 0.019, P < 0.001, and P = 0.002, respectively) ( Figure 1D ). There was no difference among the relative expressions of osteoclast-related genes (OSCAR/TRAP/RANKL/Cathepsin K, respectively) between 1 μg/mL (median: 2.672, 1.231, 2.056, and 1.115) and 10 μg/mL (median: 2.056, 1.084, 2.049, and 1.483) of the C4a vs. the control groups (median: 2.206, 1.341, 2.036, and 1.202) in 21 patients ( Figure 1E ).
The absorption area of the osteoclast resorption pit per view induced by 1 μg/mL (mean ± SD: 51.464 ± 11.983%) and 10 μg/mL (50.219 ± 12.067%) of C3a was also significantly increased (33.845 ± 8.331%) ( P < 0.001 and P < 0.001) in 30 NDMM patients compared to the control group ( Figure 1F and 1G ). There was no difference among the absorption areas of the osteoclast resorption pits between 1 μg/mL (mean ± SD: 35.950 ± 11.284%) and 10 μg/mL (36.279 ± 11.427%) of the C4a vs . the control group (33.594 ± 12.452%) in 15 patients ( Figure 1F and 1H ).
Together, these results suggested that C3a at concentrations of 1 μg/mL and 10 μg/mL significantly promoted the formation, differentiation, and functioning of osteoclasts in MM patients, while C4a did not.
Complement C3a promoted the formation and function of osteoclasts in NDMM patients by regulating the PI3K/PDK1/SGK3 pathway
To further identify related pathways, RNA-seq analysis was performed on the patient-derived osteoclasts of 4 NDMM patients treated with 1 μg/mL C3a and dimethyl sulfoxide. There were a total of 184 DEGs (fold-change ≥ 1.5) between the C3a and control groups, including 97 upregulated genes, and 87 downregulated genes ( Figure 2A and 2B ). Based on KEGG pathway analyses, the top 4 groups were genes involved in cytokine-cytokine receptor interaction (the main genes included CNTF/ CSF2/CSF3/IL-17A/CXCL10/CCL1/TNFSF12/TNFSF13/IFNA/IL-9 ), the phospholipase D signaling pathway (the main genes included PLA2G4B/IGFBP2/SHC1/SHC4/MS4A6A/GOLGA6L10/GOLGA8K/GPCR/RTK/IL-8 ), the phosphoinositide 3-kinase (PI3K) signaling pathway (including the PI3K-Akt signaling pathways and AKT-independent PI3K signaling pathways) (the main genes included PIK3CA/PDK1/SGK3/POSTN/COL1A1/COL1A2/CREB1/MDM2/IKBKG/NTRK2/CDK2/TCL1A/AKT3 ), and cell adhesion molecules (CAMs) (the main genes included HAPLN1/F13A1/ANGPTL2/CEACAM6/PTPRS/CLEC19A/ITGB1/ITGB2/SELE/NRCAM/CFPG2/SDC/PTPRF ) ( Figure 2C ). The PI3K signaling pathway mainly controls cell growth, transcription and translation, cell proliferation, cell movement, and glycogen metabolism, and is most likely to be related to osteoclast formation, differentiation, and maturation. DEGs on the PI3K signaling pathway that were upregulated in at least 3 patients were therefore selected and analyzed, as shown in the heat map ( Figure 2D ).

The upregulated DEGs were validated by qRT-PCR between the C3a and control group in 24 NDMM patients. There was no statistically significant difference with regards to the relative expressions of POSTN/COL1A1/COL1A2/CREB1/MDM2/IKBKG/NTRK2/CDK2/TCL1A/AKT3 between the C3a and control groups ( Figure 2E ). Fortunately, the relative expression levels of the PIK3CA/PDK1/SGK3 genes (median: 4.717, 2.078, and 4.428, respectively) in the C3a group (1 μg/mL) were significantly higher than in the control group (median: 2.135, 1.336, and 1.714, respectively) ( P = 0.022, P < 0.001, and P = 0.001, respectively) ( Figure 2F–2H ). We speculated that the promotional effect of C3a on osteoclasts in NDMM patients may have been related to upregulation of the PIK3CAA , PDK1 , and SGK3 genes.
The expression levels of PI3K p110α/Akt/p-Akt/PDK1/SGK3/p-SGK3 proteins in the PI3K signaling pathway were detected by Western blot. There was no significant difference in the relative grayscale of proteins Akt/p-Akt/SGK3 between the C3a (1 μg/mL) and control groups ( Figure 2I and 2J ). The relative grayscale levels of PI3K p110α/PDK1/p-SGK3 proteins (mean ± SD: 2.012 ± 0.286, 2.481 ± 0.228, and 1.956 ± 0.519, respectively) in the C3a group (1 μg/mL) were significantly higher than in the control group (mean ± SD: 0.967 ± 0.147, 1.341 ± 0.278, and 0.709 ± 0.494, respectively) ( P = 0.005, P = 0.005, and P = 0.039, respectively) ( Figure 2I and 2J ).
The expression levels of PI3K and PDK1 (both upstream of SGK3 and Akt), and p-SGK3 proteins increased, while the expressions of Akt and p-Akt proteins did not change, suggesting that C3a promoted the formation, differentiation, and functioning of osteoclasts in MM patients by regulating the PI3K/PDK1/SGK3 pathway.
Treatment with SGK inhibitor significantly inhibited the promoting effect of C3a on osteoclasts in NDMM patients
Ackermann et al. 19 and Basnet et al. 20 showed that the SGK inhibitor, EMD638683, could inhibit SGK3 up to 75% at a concentration of 1 μM. EMD638683 was therefore selected as the SGK3 inhibitor in this study.
The effects of EMD638683 on the expressions of SGK3 and p-SGK3 proteins in osteoclasts of NDMM patients were validated by Western blot. The results showed that the expressions of SGK3 and p-SGK3 proteins in the EMD638683 group were significantly lower than in the control group, confirming that EMD638683 inhibited SGK3 and p-SGK3 in the osteoclasts of NDMM patients ( Figure 3A ).

The SGK inhibitor (EMD638683) was then added to a culture system of osteoclasts activated by C3a (1 μg/mL) to evaluate whether the formation and function of osteoclasts were inhibited. Treatment with EMD638683 significantly reduced the osteoclasts areas per view in the 1 μM (mean ± SD: 39.244 ± 9.089%) and 10 μM (39.299 ± 9.587%) EMD638683-treated group, when compared to the control group (0 μM) (54.884 ± 12.837%) ( P < 0.001 and P < 0.001, respectively) ( Figure 3B and 3C ). The relative expressions of the osteoclast-related genes OSCAR/RANKL/TRAP /cathepsin K in the 1 μM (median: 0.869, 1.097, 0.902, and 1.328, respectively) and 10 μM (median: 0.703, 1.391, 0.843, and 1.418, respectively) EMD638683 treatment groups were significantly decreased compared to the control (median: 2.270, 3.024, 2.208, and 3.237, respectively) (1 μM: P = 0.015, P = 0.002, P = 0.003, and P = 0.015, respectively; 10 μM: P = 0.012, P = 0.006, P < 0.001, and P = 0.017, respectively) ( Figure 3D ). The absorption areas per view of osteoclast resorption pits in the 1 μM (mean ± SD: 35.383 ± 7.794%) and 10 μM (32.886 ± 8.993%) EMD638683-treated groups were significantly reduced compared to the control group (49.358 ± 11.856%) ( P < 0.001 and P < 0.001, respectively) ( Figure 3E and 3F ). Under the promotion of C3a, no difference was observed between the effects of 1 μM and 10 μM EMD638683 on the osteoclasts areas, relative expression of osteoclast-related genes, or absorption area of osteoclast resorption pits.
Discussion
MM is a hematological malignancy caused by the proliferation of plasma cells resulting in multiple organ function impairment. Osteoclast activation and osteoblast inhibition results in metabolism imbalance in bones. Over 80% of MM patients suffer from bone disease to different degrees at diagnosis, which seriously affects the quality of life and their prognoses 29 . The study of the pathogenesis of MBD is key for the identification of sensitive biomarkers for early diagnosis and the identification of new therapeutic targets for MBD. The complement system is an important part of the immune system and is associated with bone reconstruction and destruction 5 , 30 . MBD is associated with an imbalance between bone remodeling and destruction. Whether the complement system is related to this imbalance has not yet been reported. We found that serum levels of C3, C4, and C4a in NDMM patients were significantly correlated with the severity of bone disease, as described previously 14 . Many studies have focused on the function of osteoclasts; however, the mechanism of osteoclast activation is complex, and plays a critical role in the pathogenesis of MBD 3 , 4 . In this study, we characterized the effects and potential mechanisms of C3a/C4a on osteoclasts in MM patients.
We showed that C3a significantly promoted the formation and function of osteoclasts in MM patients, while the C4a complement did not. The C3a complement is a common focal point of three complement activation pathways, and has been extensively studied in bone metabolism. Several studies have indicated that C3a activates osteoclasts, consistent with our results on the effects of C3a on osteoclasts in MM patients. Sato et al. 11 found that the differentiation of osteoclasts in mouse bone marrow was inhibited after the addition of anti-C3 antibody. Using bone marrow cells from wild-type and C3-deficient mice (C3 −/− ), C3 −/− bone marrow cultures were found to generate fewer osteoclasts than wild-type marrow cells, with several features of osteoclast formation impaired in the bone marrow cultures of C3 −/− mice 10 . MacKay et al. 31 showed that C3 deficiency reduced bone loss at ovariectomy and may improve mechanical properties. Moreover, Matsuoka et al. 32 found that osteoclasts secreted C3a, which stimulated osteoblastogenesis and was involved in bidirectional communication between osteoblasts and osteoclasts.
No difference was observed between the 1 μg/mL and 10 μg/mL C3a groups regarding the effects on osteoclasts of MM patients. We hypothesized that C3a may function through the C3a receptor, which was already saturated at 1 μg/mL C3a. Schraufstatter et al. 18 found that chemoattractants of C3a for human MSCs showed increased concentrations of 0, 30 nM, and 100 nM, and peaked at concentrations of 100 nM, before gradually decreasing at concentrations of 300 nM and 1 μM. This suggested that the concentration dependence on C3a for the activation of osteoclasts in MM patients may occur between 0 and 1 μg/mL.
The effects of C4a on bone metabolism have not been previously reported. However, Zheng et al. 33 found that serum levels of C4 in patients with MM increased, suggesting that C4 was associated with the pathogenesis of MM. Dowling et al. 34 used a label-free mass spectrometry-based methodology to evaluate serum samples of monoclonal gammopathy of undetermined significance, smoldering multiple myeloma, MM patients with no bone disease, and MM patients with high bone disease. They found that the level of C4a increased in high bone disease MM patients, and identified C4a as a novel candidate biomarker associated with bone disease. However, these results were only obtained from correlation studies and not from functional experiments in vitro and in vivo . Our study is the first to report the effects of C4a on osteoclasts in MM patients. The results suggested that C4a at concentrations of 1 μg/mL and 10 μg/mL could not promote the formation, differentiation, or function of osteoclasts in MM patients. This may be due to C4a affecting osteoclasts by acting on other complement proteins, or may be simply due to C4a inhibiting osteoblasts, with little effect on osteoclasts. Barnum et al. 35 reported that activation of complements led to the generation of 3 anaphylatoxins: C3a, C4a, and C5a. C3a and C5a share a similar functional profile, which includes modulation of the innate and adaptive immune responses, cell homing, and tissue regeneration. In contrast, C4a has no similar function due to the absence of the C4a receptor and the inability of C4a to signal to the C3a and C5a receptors. This may also explain why C4a had no effect on osteoclasts.
Based on KEGG pathway analyses, the top 4 groups were genes involved in cytokine-cytokine receptor interaction, the phospholipase D signaling pathway, the PI3K signaling pathway, and CAMs. Cytokine-cytokine receptor interaction regulates cell growth and differentiation, and is involved in immunity, inflammation, and wound healing. The phospholipase D signaling pathway participates in various cellular signal transductions. CAMs are involved in the tight and adheren junctions of antigen-presenting cells, T cells, endothelial cells, leukocytes, neurons, Schwann cells, epithelial cells, ciliary bodies, and myoblasts. The PI3K signaling pathway mainly controls cell growth, transcription and translation, cell proliferation, cell movement, and glycogen metabolism, which is most likely related to osteoclast formation, differentiation, and maturation. We selected the PI3K signaling pathway and showed that C3a promoted the formation, differentiation, and function of osteoclasts in MM patients by regulating the PI3K/PDK1/SGK3 pathway. PI3K signaling plays a central role in cellular physiology, coordinating insulin signaling during organism growth and mediating critical cellular processes, such as glucose homeostasis, protein synthesis, cell proliferation, and survival 36 , 37 . The PI3K signaling pathway includes PI3K-Akt signaling pathways and AKT-independent PI3K signaling pathways (including the PDK1-mTORC2-SGK axis, Rac signaling, and TEC family kinase BTK) 38 . SGKs, which consist of three isoforms, SGK1, SGK2, and SGK3, are critical mediators of AKT-independent signaling downstream of PDK1 and mTORC2 in cancer. The phosphorylation of PDK1 is activated by mTORC2, which is a multi-protein complex consisting of mTOR, Rictor, mMM Sin1, mLST8, Protor1/2, and Deptor 39 . SGK3 is a unique member within the SGK family because it contains an N-terminal PX domain, which can regulate cell processes similar to AKT kinase signaling in survival, migration, and growth signaling 40 . The amplification and overexpression of SGK3 is more common than AKT in hepatocellular carcinoma, suggesting it may have a greater functional significance in the biology of this type of cancer 41 . The activation of SGK signaling has also recently emerged as an important mechanism of resistance to PI3K and AKT inhibitors. For example, the treatment of breast cancer cells with PI3K or AKT inhibitors results in increased expression and the activation of SGK3, which depends on hVps34 for activation by PDK1 and mTORC2 42 . In our study, PDK1 was found to be an important kinase for the phosphorylation of SGK3. Thus, it is reasonable that increased PDK1 expression led to increased phosphorylation of SGK3, indicating that C3a promoted the formation, differentiation, and function of osteoclasts in NDMM patients by the PI3K/PDK1/SGK3 signaling pathways ( Figure 4 ).

After treatment with a SGK inhibitor (EMD638683), up to 75% of SGK3 was inhibited at a concentration of 1 μM. In addition, the promotional effect of C3a on osteoclasts in NDMM patients was significantly inhibited. It has been reported that EMD638683 inhibits the growth and proliferation of colon carcinoma cells 43 . EMD638683 could therefore serve as a template for drugs counteracting hypertension in individuals with type 2 diabetes and metabolic syndrome 19 . In our study, treatment with EMD638683 induced a significant inhibitory effect on MM-derived osteoclasts, providing important evidence for the identification of new therapeutic targets and strategies for MBD patients. Further studies will be needed to validate its inhibitory effect on MBD in an animal model, such as mice.
Conclusions
C3a activated osteoclasts by regulating PI3K/PDK1/SGK3 pathways in MM patients, which was reduced by treatment with a SGK inhibitor. Overall, this study identified new potential therapeutic targets and strategies for MBD patients.
Supporting Information
Click here for additional data file.
Conflict of interest statement No potential conflicts of interest are disclosed.
Grant support
This work was supported by the National Natural Science Foundation of China (Grant Nos. 81770110, 81900131, and 82000219), the Anticancer Major Special Project of Tianjin (Grant No. 12ZCDZSY18000), the Tianjin Municipal Natural Science Foundation (Grant Nos. 18JCYBJC27200 and 18JCQNJC80400), the Tianjin Education Commission Research Project (Grant Nos. 2018KJ043 and 2018KJ045), the Tianjin Health and Family Planning Commission (Grant No. 15KG150), the Youth Incubation Fund of Tianjin Medical University General Hospital (Grant No. ZYYFY2019020), and the Tianjin Science and Technology Planning Project (Grant No. 20YFZCSY00060).
References
1. Roodman GD Mechanisms of bone metastasis N Engl J Med 2004 350 1655 64 15084698 10.1056/NEJMra030831
2. Terpos E Berenson J Cook RJ Lipton A Coleman RE Prognostic variables for survival and skeletal complications in patients with multiple myeloma osteolytic bone disease Leukemia 2010 24 1043 9 20376081 10.1038/leu.2010.62
3. Terpos E Ntanasis-Stathopoulos I Dimopoulos MA Myeloma bone disease: from biology findings to treatment approaches Blood 2019 133 1534 9 30760454 10.1182/blood-2018-11-852459
4. Oranger A Carbone C Izzo M Grano M Cellular mechanisms of multiple myeloma bone disease Clin Dev Immunol 2013 2013 289458 23818912 10.1155/2013/289458 PMC3681224
5. Schoengraf P Lambris JD Recknagel S Kreja L Liedert A Brenner RE Does complement play a role in bone development and regeneration? Immunobiology 2013 218 1 9 22464814 10.1016/j.imbio.2012.01.020
6. Andrades JA Nimni ME Becerra J Eisenstein R Davis M Sorgente N Complement proteins are present in developing endochondral bone and may mediate cartilage cell death and vascularization Exp Cell Res 1996 227 208 13 8831558 10.1006/excr.1996.0269
7. Sakiyama H Inaba N Toyoguchi T Okada Y Matsumoto M Moriya H Immunolocalization of complement C1s and matrix metalloproteinase 9 (92kDa gelatinase/type IV collagenase) in the primary ossification center of the human femur Cell Tissue Res 1994 277 239 45 8082118
8. Toyoguchi T Yamaguchi K Nakagawa K Fukusawa T Moriya H Sakiyama H Change of complement C1s synthesis during development of hamster cartilage Cell Tissue Res 1996 285 199 204 8766156 10.1007/s004410050637
9. Ignatius A Schoengraf P Kreja L Liedert A Recknagel S Kandert S Complement C3a and C5a modulate osteoclast formation and inflammatory response of osteoblasts in synergism with IL-1beta J Cell Biochem 2011 112 2594 605 21598302 10.1002/jcb.23186 PMC3158833
10. Tu Z Bu H Dennis JE Lin F Efficient osteoclast differentiation requires local complement activation Blood 2010 116 4456 63 20709903 10.1182/blood-2010-01-263590 PMC2996112
11. Sato T Abe E Jin CH Hong MH Katagiri T Kinoshita T The biological roles of the third component of complement in osteoclast formation Endocrinology 1993 133 397 404 8319587 10.1210/endo.133.1.8319587
12. Kovtun A Bergdolt S Hagele Y Matthes R Lambris JD Huber-Lang M Complement receptors C5aR1 and C5aR2 act differentially during the early immune response after bone fracture but are similarly involved in bone repair Sci Rep 2017 7 14061 10.1038/s41598-017-14444-3 PMC5656620 29070810
13. Teo BH Bobryshev YV Teh BK Wong SH Lu J Complement C1q production by osteoclasts and its regulation of osteoclast development Biochem J 2012 447 229 37 22812635 10.1042/BJ20120888
14. Jiang F Liu H Liu Z Song J Li L Ding K High serum levels of complements C3 and C4 as novel markers for myeloma bone disease Ann Hematol 2017 96 331 3 27924370 10.1007/s00277-016-2863-z
15. Rajkumar SV Dimopoulos MA Palumbo A Blade J Merlini G Mateos MV International Myeloma Working Group updated criteria for the diagnosis of multiple myeloma Lancet Oncol 2014 15 e538 48 25439696 10.1016/S1470-2045(14)70442-5
16. Palumbo A Avet-Loiseau H Oliva S Lokhorst HM Goldschmidt H Rosinol L Revised international staging system for multiple myeloma: a report from international myeloma working group J Clin Oncol 2015 33 2863 9 26240224 10.1200/JCO.2015.61.2267 PMC4846284
17. Jiang F Liu H Liu Z Yan S Chen J Shao Q Deficient invariant natural killer T cells had impaired regulation on osteoclastogenesis in myeloma bone disease J Cell Mol Med 2018 22 2706 16 29473714 10.1111/jcmm.13554 PMC5908096
18. Schraufstatter IU Discipio RG Zhao M Khaldoyanidi SK C3a and C5a are chemotactic factors for human mesenchymal stem cells, which cause prolonged ERK1/2 phosphorylation J Immunol 2009 182 3827 36 19265162 10.4049/jimmunol.0803055
19. Ackermann TF Boini KM Beier N Scholz W Fuchss T Lang F EMD638683, a novel SGK inhibitor with antihypertensive potency Cell Physiol Biochem 2011 28 137 46 21865856 10.1159/000331722
20. Basnet R Gong GQ Li C Wang MW Serum and glucocorticoid inducible protein kinases (SGKs): a potential target for cancer intervention Acta Pharm Sin B 2018 8 767 71 30245963 10.1016/j.apsb.2018.07.001 PMC6146383
21. Liu H Peng F Liu Z Jiang F Li L Gao S CYR61/CCN1 stimulates proliferation and differentiation of osteoblasts in vitro and contributes to bone remodeling in vivo in myeloma bone disease Int J Oncol 2017 50 631 9 28035364 10.3892/ijo.2016.3815
22. Huang J Liang X Xuan Y Geng C Li Y Lu H A reference human genome dataset of the BGISEQ-500 sequencer Gigascience 2017 6 1 9 10.1093/gigascience/gix024 PMC5467036 28379488
23. Huang J Liang X Xuan Y Geng C Li Y Lu H Erratum to: a reference human genome dataset of the BGISEQ-500 sequencer Gigascience 2018 7 giy144 10.1093/gigascience/giy144 PMC6266704 30500904
24. Cock PJ Fields CJ Goto N Heuer ML Rice PM The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants Nucleic Acids Res 2010 38 1767 71 20015970 10.1093/nar/gkp1137 PMC2847217
25. Langmead B Trapnell C Pop M Salzberg SL Ultrafast and memory-efficient alignment of short DNA sequences to the human genome Genome Biol 2009 10 R25 19261174 10.1186/gb-2009-10-3-r25 PMC2690996
26. Kim D Langmead B Salzberg SL HISAT: a fast spliced aligner with low memory requirements Nat Methods 2015 12 357 60 25751142 10.1038/nmeth.3317 PMC4655817
27. Li B Dewey CN RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome BMC Bioinformatics 2011 12 323 21816040 10.1186/1471-2105-12-323 PMC3163565
28. Liu Z Jing Q Wang Y Li Y Mi F Xiang C The short-term effect of histone deacetylase inhibitors, chidamide and valproic acid, on the NFkappaB pathway in multiple myeloma cells Int J Mol Med 2019 43 285 93 30387821 10.3892/ijmm.2018.3963 PMC6257846
29. Giuliani N Rizzoli V Roodman GD Multiple myeloma bone disease: pathophysiology of osteoblast inhibition Blood 2006 108 3992 6 16917004 10.1182/blood-2006-05-026112
30. Modinger Y Loffler B Huber-Lang M Ignatius A Complement involvement in bone homeostasis and bone disorders Semin Immunol 2018 37 53 65 29395681 10.1016/j.smim.2018.01.001
31. MacKay DL Kean TJ Bernardi KG Haeberle HS Ambrose CG Lin F Reduced bone loss in a murine model of postmenopausal osteoporosis lacking complement component 3 J Orthop Res 2018 36 118 28 28667799 10.1002/jor.23643 PMC5758433
32. Matsuoka K Park KA Ito M Ikeda K Takeshita S Osteoclast-derived complement component 3a stimulates osteoblast differentiation J Bone Miner Res 2014 29 1522 30 24470120 10.1002/jbmr.2187
33. Zheng MM Zhang Z Bemis K Belch AR Pilarski LM Shively JE The systemic cytokine environment is permanently altered in multiple myeloma PLoS One 2013 8 e58504 23544044 10.1371/journal.pone.0058504 PMC3609759
34. Dowling P Hayes C Ting KR Hameed A Meiller J Mitsiades C Identification of proteins found to be significantly altered when comparing the serum proteome from Multiple Myeloma patients with varying degrees of bone disease BMC Genomics 2014 15 904 25322877 10.1186/1471-2164-15-904 PMC4213504
35. Barnum SR C4a: an anaphylatoxin in name only J Innate Immun 2015 7 333 9 25659340 10.1159/000371423 PMC6738802
36. Engelman JA Targeting PI3K signalling in cancer: opportunities, challenges and limitations Nat Rev Cancer 2009 9 550 62 19629070 10.1038/nrc2664
37. Fruman DA Rommel C PI3K and cancer: lessons, challenges and opportunities Nat Rev Drug Discov 2014 13 140 56 24481312 10.1038/nrd4204 PMC3994981
38. Lien EC Dibble CC Toker A PI3K signaling in cancer: beyond AKT Curr Opin Cell Biol 2017 45 62 71 28343126 10.1016/j.ceb.2017.02.007 PMC5482768
39. Laplante M Sabatini DM mTOR signaling in growth control and disease Cell 2012 149 274 93 22500797 10.1016/j.cell.2012.03.017 PMC3331679
40. Bruhn MA Pearson RB Hannan RD Sheppard KE AKT-independent PI3-K signaling in cancer - emerging role for SGK3 Cancer Manag Res 2013 5 281 92 24009430 10.2147/CMAR.S35178 PMC3762672
41. Liu M Chen L Chan TH Wang J Li Y Li Y Serum and glucocorticoid kinase 3 at 8q13.1 promotes cell proliferation and survival in hepatocellular carcinoma Hepatology 2012 55 1754 65 22262416 10.1002/hep.25584
42. Bago R Sommer E Castel P Crafter C Bailey FP Shpiro N The hVps34-SGK3 pathway alleviates sustained PI3K/Akt inhibition by stimulating mTORC1 and tumour growth EMBO J 2016 35 1902 22 27481935 10.15252/embj.201693929 PMC5007552
43. Towhid ST Liu GL Ackermann TF Beier N Scholz W Fuchss T Inhibition of colonic tumor growth by the selective SGK inhibitor EMD638683 Cell Physiol Biochem 2013 32 838 48 24081014 10.1159/000354486