Honglei Ren1,2†
Shanlin Xiang
Aiguo Liu
Zhenming Hu- 1Department of Orthopedic Surgery, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of Orthopedic Surgery, ChongQing Red Cross Hospital (People’s Hospital of JiangBei District), Chongqing, China
- 3Department of Orthopedics, Shanghai Pudong Hospital, Fudan University Pudong Medical Center, Shanghai, China
- 4Department of Orthopedic Surgery, The First Affiliated Hospital of Henan University, Kaifeng, China
- 5Department of Orthopedic Surgery, The University-Town Hospital of Chongqing Medical University, Chongqing, China
Purpose: Osteosarcoma (OS) is the most common type of primary malignant bone tumor. Transducing a functional TP53 gene can effectively inhibit OS cell activity. Poly lactic acid-glycolic acid (PLGA) nanobubbles (NBs) mediated by focused ultrasound (US) can introduce exogenous genes into target cells in animal models, but this technique relies on the passive free diffusion of agents across the body. The inclusion of superparamagnetic iron oxide (SPIO) in microbubbles allows for magnetic-based tissue localization. A low-intensity-focused ultrasound (LIFU) instrument was developed at our institute, and different intensities of LIFU can either disrupt the NBs (RLI-LIFU) or exert cytocidal effects on the target tissues (RHI-LIFU). Based on these data, we performed US-magnetic-mediated TP53-NB destruction and investigated its ability to inhibit OS growth when combined with LIFU both in vitro and in vivo.
Methods: Several SPIO/TP53/PLGA (STP) NB variants were prepared and characterized. For the in vitro experiments, HOS and MG63 cells were randomly assigned into five treatment groups. Cell proliferation and the expression of TP53 were detected by CCK8, qRT-PCR and Western blotting, respectively. In vivo, tumor-bearing nude mice were randomly assigned into seven treatment groups. The iron distribution of Perls’ Prussian blue-stained tissue sections was determined by optical microscopy. TUNEL-DAPI was performed to examine apoptosis. TP53 expression was detected by qRT-PCR and immunohistochemistry.
Results: SPIO/TP53/PLGA NBs with a particle size of approximately 200 nm were prepared successfully. For in vitro experiments, ultrasound-targeted transfection of TP53 overexpression in OS cells and efficient inhibition of OS proliferation have been demonstrated. Furthermore, in a tumor-bearing nude mouse model, RLI-LIFU-magnetic-mediated SPIO/TP53/PLGA NBs increased the transfection efficiency of the TP53 plasmid, resulting in apoptosis. Adding RHI-LIFU to the treatment regimen significantly increased the apoptosis of OS cells in vivo.
Conclusion: Combining LIFU and US-magnetic-mediated SPIO/TP53/PLGA NB destruction is potentially a novel noninvasive and targeted therapy for OS.
Introduction
Osteosarcoma (OS) is the most common type of primary malignant bone tumor and is characterized by local invasion and distant metastasis. The current clinical treatment modalities for OS include surgical resection of primary tumors and systemic chemotherapy. Research has shown that the application of neoadjuvant chemotherapy increases the tumor-free survival of OS from 1% to 50%–65%. However, resistance of OS cells to chemotherapeutic drugs remains a major limiting factor for current treatment modalities for this disease. Additionally, surgical treatment for OS may cause significant physical trauma. Therefore, a new, less traumatic therapy would be of great significance.
In recent years, gene therapy for OS has emerged as an exciting novel research field. TP53 was found to be the most frequently mutated gene among the tumor suppressor genes, with approximately 60% of human tumors presenting TP53 mutations. Current literature shows that OS exhibits a wide range of gene mutations and molecular alterations. However, apart from mutations in the TP53 and/or retinoblastoma (RB) genes, no novel gene mutations have been identified. Several studies have demonstrated that transfection of a functional copy of the TP53 gene could inhibit the growth of OS cells (Bekeredjian et al., 2007; Chen et al., 2012; Chen et al., 2010).
The introduction of exogenous genes relies on gene vectors, and the stability of the gene will be relatively poor without an optimal gene vector. In recent years, great progress has been made in vector research, and high transfection efficiencies have been shown for viral vectors. However, due to the high immunogenicity of viruses, their application is limited. Nonviral vectors such as liposomes are safe and less toxic compared with viral vectors, but instability and inefficiency remain constraints that cannot be ignored (Chertok et al., 2007; Du et al., 2011; Elzeny et al., 2017; Hida et al., 2011). Therefore, finding a safe, efficient, and targeted gene therapy transfection system has become the focus of researchers.
Poly (lactic-co-glycolic acid) (PLGA) nanobubbles (NBs) are nonviral vectors that offer significant advantages, including slow release, penetration, and ‘targeting’. These characteristics make PLGA NBs an ideal choice for drug or gene carriers (Huang et al., 2012; Hynynen et al., 2006; Jiang and Dalby, 2023). However, the “targeting” of this delivery system is still passive and depends on the free diffusion of the medium across the body. Superparamagnetic iron oxide (SPIO, diameter <10 nm) nanoparticles are a type of magnetic nanoparticle that are physically sensitive to external magnetic fields, and magnetic targeting (MT) can be applied to actively enhance the concentration in the target area (Lu et al., 2003; Sirsi et al., 2012). Hence, we hypothesized that the incorporation of SPIO and MT into TP53/PLGA NBs could serve as an effective strategy to augment gene transfection in the targeted area.
Recently, ultrasound (US)-mediated microbubble destruction (UTMD) has undergone rapid development. UTMD has been demonstrated to be a safe and effective method for gene transfection by numerous studies (Sun et al., 2023; Tang et al., 2012; Wells, 2010; Wu et al., 2020). UTMD has been successfully used for gene transfection of the kidney, testis, heart, pancreas, lung, skin, uterus, brain, retina diseases, and spinal cord (Sitta and Howard, 2021; Steiniger et al., 2004; Xenariou et al., 2007). Low-intensity-focused ultrasound (LIFU) can disrupt target nanobubbles at low intensity (referred to as relatively low-intensity mode, RLI mode), whereas at higher intensities (referred to as relatively high-intensity mode, RHI mode), LIFU induces range-controllable cytocidal effects on tissues (Yamaguchi et al., 2011; Yin et al., 2014; Zhong et al., 2012). To the best of our knowledge, few studies have investigated gene therapy for OS through US-magnetic-mediated NB destruction, and direct OS treatment using RHI-mode LIFU has not been reported. Thus, it remains unclear whether US-magnetic-mediated gene transfection in combination with LIFU represents an effective therapeutic strategy for OS.
In this study, we performed US-magnetic-mediated TP53-NB destruction and evaluated its potential to inhibit OS growth when combined with LIFU in both in vitro and in vivo experiments.
Material and methods
Preparations of SPIO/PLGA-NBs (SP-NBs), TP53/PLGA-NBs (TP-NBs) and SPIO/TP53/PLGA-NBs (STP-NBs)
NBs were freshly prepared prior to use. SPIO nanoparticles (mean diameter = 5 nm, concentration = 25 mg/mL) were purchased from Ocean Nanotech Co., Ltd. The TP53 plasmid was purchased from OriGene Technologies, Inc. and extracted using an endotoxin-free plasmid extract kit (enhanced) (DP120; TIANGEN BIOTECH (BEIJING) Co., Ltd.). PLGA (50:50, molecular weight [MW] = 12,000) was purchased from Polyscitech Co., Ltd. Polyvinyl alcohol (PVA, MW = 30,000-70,000) was purchased from Sigma Co., Ltd.
As reported previously, STP-NBs were prepared by the double-emulsion method. (Kiwada et al., 1986). First, 50 mg PLGA was completely dissolved in 2 mL dichloromethane, 600 µg TP53 plasmids and 20 µL SPIO were added to the solution. The mixture was emulsified in an ice bath using an ultrasonic processor for 5 min at a power of 175 W. Then, 5 mL cold PVA solution (4% W/V) was added to the initial emulsion, and an ultrasonic processor was used to further emulsify the product in the ice bath with a power of 125 W for 3 min. Then, 10 mL of 1% isopropanol solution was added to the second emulsion. The final solution was stirred with a magnetic agitator for 3 h. The volatilized solution was centrifuged at a low temperature centrifuge at 4°C for 10 min at 10000×g. The precipitate was washed with deionized water, and the supernatant was collected for further detection. The centrifugal washing step was repeated 3 times. Finally, the precipitate after washing was temporarily stored at 4°C for further use. SP-NBs and TP-NBs were similarly prepared without the addition of either TP53 or SPIO.
Characterizations of TP-NBs and STP-NBS
A laser particle size analyzer system (Zeta SIZER, Malvern) was used to analyze the mean diameter, zeta potential, and polydispersity index (PDI) of the samples. Transmission electron microscopy (TEM; Hitachi H-7600) was used to evaluate the shape of the NBs, whereas scanning electron microscopy (SEM; JEOL JSM-7800F) was used to assess the morphological characteristics of the nanoparticles. (Zhao et al., 2015). The iron content of the STP-NBs was determined by the o-phenanthroline method. The encapsulation efficiency of SPIO is calculated as follows: encapsulation efficiency = W1/W2×100%, where W1 is the amount of SPIO in STP-NBs and W2 is the total amount of SPIO used to prepare STP-NBs. (Kiwada et al., 1986; Ting et al., 2012). We compared SPIO encapsulation efficiencies at five SPIO total inputs (0.1, 0.3, 0.5, 0.7 and 0.9 mg). The SPIO encapsulation efficiency test was repeated three times. The gene encapsulation efficiency was determined by an ultraviolet spectrophotometer (Thermo NanoDrop 2000; Thermo Fisher Scientific). The encapsulation efficiency of TP53 is calculated as follows: encapsulation efficiency = (W2-W1/W2 ×100)%, where W1 is the number of genes in the supernatant after centrifugation and W2 is the total amount of TP53 used to prepare STP-NBs. (Zhang et al., 2014; Wang et al., 2019). We compared TP53 encapsulation efficiencies at five TP53 plasmid total inputs (0.1, 0.2, 0.4, 0.8, and 1 mg). The TP53 encapsulation efficiency test was repeated three times. To assess gene release behavior, in vitro sonication was performed using an RLI-mode LIFU. The frequency was 0.95 MHz, the power was modulated at a pulse modulation of 300 Hz, and the depth of focus of the probe was 15 mm. The ultrasonic power is 1.5 W and lasts 60 s. STP-NBs (TP-NBs) from 10 mL phosphate buffered saline (PBS) were transferred to dialysis bags (molecular weight cut-off: 10,000 Da) and placed in a 100 mL PBS reservoir after ultrasound. After the appropriate time interval, 1 mL dialysate was taken, and the gene concentration was determined by ultraviolet spectrophotometry. One milliliter of fresh PBS was added to the reservoir to keep the liquid volume constant. The cumulative release ratio of TP53 released was calculated. To investigate the magnetization capacity, STP-NBs were placed under a magnetic field, and the adsorption time was recorded.
In vitro studies on MG63 and HOS osteosarcoma cells
In vitro focused ultrasound (FUS) sonication setup
The MG63 and HOS cell lines used in this study were acquired from Professor Tong-Chuan He (Chongqing Medical University, Chongqing, China). To assess gene release behavior, in vitro sonication was performed by using RLI-mode LIFU. The sonication parameters of RLI-mode were set at 1.5 W for 5 min (Figure 1A).
Figure 1. (A) Schematic of the in vitro experimental setup, (B) Schematic of the in vivo experimental setup.
HOS and MG63 cells were randomly assigned into five treatment groups:
1. Control group: 50 µL PBS was added to 24-well plates, and 10 µL PBS was added to 96-well plates.
2. TP53 + SP-NBs group: 25 µL TP53 plasmid (1 µg) and 25 µL SP-NB solution were added to 24-well plates; 5 µL TP53 plasmid (0.2 µg) and 5 µL SP-NB solution were added to 96-well plates.
3. TP53 + US group: 25 µL TP53 plasmid (1 µg) and 25 µL PBS solution were added to 24-well plates, which were then subjected to RLI-mode LIFU irradiation; 5 µL TP53 plasmid (0.2 µg) and 5 µL PBS solution were added to 96-well plates, which were then subjected to RLI-mode LIFU irradiation (Figure 1A).
4. TP-NBs + US group: 50 µL TP-NB (TP53, 1 µg) solution was added to 24-well plates, which were then subjected to RLI-mode LIFU irradiation; 10 µL TP-NB (TP53, 0.2 µg) solution was added to 96-well plates, which were then subjected to RLI-mode LIFU irradiation (Figure 1A).
5. STP-NBs + US group: 50 µL STP-NB (TP53, 1 µg) solution was added to 24-well plates, which were then subjected to RLI-mode LIFU irradiation; 10 µL STP-NB (TP53, 0.2 µg) solution was added to 96-well plates, which were then subjected to RLI-mode LIFU irradiation (Figure 1A).
Cell proliferation assay
The proliferative capacity of cells was detected using a cell counting kit-8 assay (CCK8, MedChemExpress, United States). MG63 and HOS cells were seeded in 96-well plates at a density of 2000 cells/well and cultured for 36 h. Then, the cells were grouped, treated as specified above and cultured for 2 h. After that, the cells were washed three times with PBS and cultured for another 36 h in fresh cell culture medium. Finally, 10 µL CCK8 reagent was added to the wells. After the plates were incubated for 0.5–4 h, cell proliferation was detected by measuring the absorbance at 450 nm with a microplate reader. Cell viability (%) relative to that of the control cells was calculated according to the instructions of the CCK-8 assay. (Livak and Schmittgen, 2001).
Real-time PCR
MG63 and HOS cells were seeded in 24-well plates at a density of 1*105 cells/well and cultured for 36 h. Then, the cells were grouped, treated as specified above and cultured for 2 h. After that, the cells were washed three times with PBS and cultured for another 36 h in fresh cell culture medium. Finally, total RNA was isolated by using TRIzol reagent (TAKARA BIO INC) following the manufacturer’s protocols. The PrimeScript™ RT reagent kit was used to reverse transcribe RNA to complementary DNA (cDNA). SYBR Green Super Mixture (TAKARA BIO INC) was used to amplify the resulting cDNA. The P53 primer sequences were as follows: 5′- GCCATCTACAAGCAGTCACAGC-3′ (sense) and 5′- CATCCAAATACTCCACACGCAA-3′ (antisense). The GAPDH primer sequences were as follows: 5′-TCAAGAAGGTGGTGAAGCAGG-3′ (sense) and 5′-AGCGTCAAAGGTGGAGGAGTG-3′ (antisense). The CFX-Connect real-time PCR system (Bio-Rad) was used to conduct qRT-PCR in 96-well plates. The amplification conditions were as follows: 95°C for 3 min, followed by 45 cycles of 95°C for 10 s and 58°C for 30 s. Each qRT-PCR analysis was carried out in triplicate, and the 2−ΔΔCT method was used to analyze the data.
Western blotting analysis
MG63 and HOS cells were seeded in 24-well plates at a density of 1*105 cells/well and cultured for 36 h. Then, the cells were grouped, treated as specified above and cultured for another 2 h. After that, the cells were washed three times with PBS and cultured for an additional 36 h in fresh cell culture medium. Finally, proteins were homogenized in lysis buffer and phenylmethanesulfonyl fluoride (Beyotime, China). The bicinchoninic acid assay (Beyotime Institute of Biotechnology, China) was used to determine the protein concentrations. The lysates were centrifuged at a low temperature centrifuge at 4°C for 10 min at 10000×g, and the supernatants were collected and transferred into separate enzyme-free tubes. The same amounts of protein were separated by SDS-PAGE and transferred to a polyvinylidene fluoride (PVDF) membrane. The membranes were blocked in 5% skim milk for 1 h at room temperature and incubated with specific primary antibodies at 4°C overnight. The membranes were sealed in 5% skim milk at room temperature for 1 h and incubated with specific primary antibody (anti-P53, Anffinity, United States; 1:1,000; anti-GAPDH, Anffinity, United States; 1:10000) overnight at 4°C. After washing with TBST three times, membranes were incubated with the secondary antibodies at room temperature for 1 h. An Immobilon Western Chemiluminescent Kit (NCM Biotech, China) was used to visualize the labeled protein bands. The band density of P53 protein was normalized to GAPDH, and all western blots were quantified using ImageJ software.
In vivo studies on osteosarcoma
In vivo experimental setup
This study was conducted according to the guidelines for the care and use of experimental animals from the National Institutes of Health, and the experimental schemes were approved by the Animal Ethics Committee of Chongqing Medical University (Figure 1B).
HOS cells were digested and collected at 90% density and resuspended in serum-free 1,064 medium to 5 × 107 cells/mL. Each athymic nude mouse (4 weeks old, Beijing HFK Bioscience Corporation, China) was subcutaneously injected with 100 μL cell suspension (approximately 5 × 106 cells). (Lin et al., 2019). To assess gene release behavior, in vitro sonication was performed by using LIFU. The sonication parameters of the RLI-mode were 1.5 W for 5 min, and the parameters of the RHI-mode were set as follows: power, 12 W; sound intensity, 19.2 w/cm2 at focus; depth of focus of 1.5 cm. Two weeks after subcutaneous implantation, the tumor-bearing nude mice were randomly assigned into seven treatment groups:
1. Control group: 1 mL NS was injected via the tail vein.
2. TP53 + SP-NBs + MT group: 0.5 mL TP53 plasmid (50 µg) and 0.5 mL SP-NB solution were injected via the tail vein. For animals undergoing the MT procedure, a permanent magnet was placed tightly to the surface skin of the tumor nodules for 30 min.
3. TP53 + US group: 0.5 mL TP53 plasmid (50 µg) solution and 0.5 mL NS were injected via the tail vein and subjected to RLI-mode LIFU irradiation.
4. TP-NBs + MT + US group: 1 mL TP-NB (TP53, 50 µg) solution was injected via the tail vein, and a RuFe boron magnet was applied tightly to the surface skin of the tumor nodules for 30 min before the mice were subjected to RLI-mode LIFU irradiation (Figure 1B).
5. STP-NBs + US group: A 1 mL STP-NB (TP53, 50 µg) solution was injected via the tail vein, and then the mice were subjected to RLI-mode LIFU irradiation.
6. STP-NBs + US + MT group: A 1 mL STP-NB (TP53, 50 µg) solution was injected via the tail vein, and then a permanent magnet was placed tightly to the surface skin overlaying the tumor nodules for 30 min before the mice were subjected to RLI-mode LIFU irradiation (Figure 1B).
7. STP-NBs + US + MT + LIFU group: A 1 mL STP-NB (TP53, 50 µg) solution was injected via the tail vein, and then a RuFe boron magnet was applied tightly to the surface skin of the tumor nodules for 30 min. Next, the mice were subjected to RLI-mode LIFU irradiation followed by RHI-mode LIFU irradiation for 10 min (Figure 1B).
All mice received their respective drug injections every 24 h, which was continued for 5 days, and tumor nodules were collected 7 days after the last transfection. The retrieved tumor tissues were fixed in 4% paraformaldehyde for 24 h and embedded in paraffin. During the treatment, all nude mice were intraperitoneally injected with pentobarbital for deep anesthesia.
Perls’ Prussian blue assay
The iron distribution was determined by examining tissue sections stained with Perls’ Prussian blue by light microscopy. (Asano et al., 2006). Briefly, the tissue sections were dewaxed with xylene, rehydrated with reduced ethanol concentration, rinsed with distilled water, and incubated with the staining solution for 20 min. Then, the tissue sections were rinsed with distilled water for 5 min and counterstained with nuclear fast red for 5 min. Prussian blue-positive cells were stained blue. Sections were imaged under a microscope, and the number of positive cells was counted.
TUNEL assay
A terminal deoxynucleotidyl transferase-mediated dUTP-biotin nick end labeling (TUNEL) staining kit was used to analyze OS apoptosis. Briefly, tissue sections were dewaxed in xylene, rehydrated by reducing the concentration of ethanol, washed in distilled water, protease treated, and incubated with TUNEL reactive mixture. Nuclear reverse staining with 4′,6-diamidino-2-phenylindole (DAPI). The nuclei of TUNEL-positive cells were dense and showed green fluorescence, suggesting apoptosis. The sections were imaged by fluorescence microscopy.
Total TP53 expression in vivo after transfection
Quantification of gene expression was determined by qRT-PCR as previously described. TP53 protein expression was measured using immunohistochemistry. (Mungarndee et al., 2008). Briefly, tumor tissue sections were incubated at 4°C overnight with a primary antibody and then incubated at room temperature for 1 h with a secondary antibody. The expression of target proteins was observed by using chromogenic DAB substrate. Sections were imaged under a microscope.
Statistical analysis
SPSS 21 software (IBM, Armonk, NY, United States) was used for statistical analysis. All data are presented as the mean ± standard error of the mean of at least three independent samples. The differences among multiple groups were determined by using one-way ANOVA (p < 0.05). The differences between two groups at the same time point were analyzed by using Student’s t-test. A p-value of less than 0.05 indicated statistical significance.
Results
TP-NBs and STP-NBs were well prepared and with excellent functional characteristics
The TP-NB zeta potential was −6.75 ± 3.06 mV, the mean diameter was 174.7 ± 70.24 nm, and the PDI was 0.195 (Figures 2A, B). The STP-NB zeta potential was −5.92 ± 3.37 mV, the mean diameter was 243.6 ± 46.51 nm, and the PDI was 0.284 (Figures 2C, D). SEM and TEM analysis revealed that TP-NBs and STP-NBs were highly dispersed and had a well-defined spherical morphology, and SPIO was observed in the core of STP-NBs (Figures 2E–H). The o-phenanthroline method was used to measure the SPIO encapsulation efficiency of STP-NBs, and the SPIO encapsulation efficiency decreased with increasing total amount of SPIO. The results revealed that the SPIO content was highest when 0.5 mg of SPIO was used for encapsulation (Figure 2I). The ultraviolet spectrophotometry method was used to measure the TP53 plasmid encapsulation efficiency of TP-NBs and STP-NBs, and when 0.6 mg of TP53 plasmid was used, the TP53 plasmid encapsulation efficiency was highest (Figure 2J). The ultraviolet spectrophotometry method was used to measure the release efficiency of TP53 from TP-NBs and STP-NBs after ultrasonic irradiation. A total of 57.1% of TP53 was released from STP-NBs within 12 h after ultrasonic irradiation; 65.1% of TP53 was released at 48 h after ultrasonic irradiation. Regarding STP-NBs, 52.9% of the loaded TP53 was released from STP-NBs at 12 h after ultrasonic irradiation, and 62.4% was released at 48 h (Figure 2K). In addition, the suspended STP-NBs were attracted by permanent magnet when placed under a magnetic field for 4 min, indicating that STP-NBs have a high magnetization ability (Figure 2L).
Figure 2. Characterization of PP-NBs and SPP-NBS NBs. (A) The mean diameter of TP-NBs, (B) The zeta potential of TP-NBs, (C) The mean diameter of STP-NBs, (D) The zeta potential of STP-NBs, (E) SEM of TP-NBs, (F) TEM of TP-NBs, (G) SEM of STP-NBs, (H) TEM of STP-NBs, (I) SPIO encapsulation efficiency of STP-NBs, (J) TP53 plasmid encapsulation efficiency of TP-NBs and STP-NBs, (K) Cumulative gene release of TP-NBs and STP-NBs, (L) Suspended STP-NBs were attracted to the permanent magnet.
Ultrasound-targeted transfection of TP53 in two kinds of OS cell lines and the overexpression of P53, resulting in efficient inhibition of OS cells proliferation in vitro studies
The 2 -ΔΔCT method was used to analyze the qRT-PCR data. () TP53 mRNA expression in MG63 cells was significantly higher in the TP-NBs + US and STP-NBs + US groups than in the other groups at 36 h after transfection (p < 0.05) (Figure 3A). Furthermore, the expression patterns were not significantly different between the TP-NBs + US and STP-NBs + US groups. and were also not significantly different between the control group, TP53 + SP-NBs group and TP53 + US group (p > 0.05) (Figure 3A). Western blotting (WB) was used to detect TP53 protein expression. A double band was observed for each group. TP53 protein expression in MG63 cells 36 h after treatment was highest in the TP-NBs + US and STP-NBs + US groups (Figure 3C). There was no significant difference between the TP-NBs + US and STP-NBs + US groups (p > 0.05) and was also not significantly different among the control group, TP53 + SP-NBs group and TP53 + US group (p > 0.05) (Figure 3C). Similar results were shown for HOS cells (Figures 3B, D).
Figure 3. (A) qRT-PCR of MG63 cells, (B) qRT-PCR of HOS cells, (C) WB of MG63 cells, (D) WB of HOS cells, (E) Proliferation potential of MG63 cells, (F) proliferation potential of HOS cells. *, p < 0.05, compared with the Control group; NS = no significance; #, p < 0.05, compared with the TP53 + US group.
The proliferation potential of MG63 cells was detected 36 h after treatment. MG63 cell activity in the TP-NBs + US group and STP-NBs + US group decreased 36 h compared with the other groups (p < 0.05). There was no significant difference in proliferation rates between the TP-NBs + US group and the STP-NBs + US group (p > 0.05), and there was also no significant difference among the control group, TP53 + STP-NBs group and TP53 + US group (p > 0.05) (Figure 3E). Similar results were shown for HOS cells (Figure 3F). So the proliferation of OS cells was efficiently inhibited by the overexpression of TP53 in vitro.
Ultrasound-targeted transfection in magnetic field enhanced the transfection efficiency of the TP53 plasmid in local area, resulting in apoptosis on the osteosarcoma model in vivo
TP53 protein expression was detected by immunohistochemistry in vivo and showed increases in groups 4, 5, 6 and 7 compared with the other three groups (Figures 4A, B). Groups 6 and 7 showed the highest TP53 protein expression, but there was no significant difference (p > 0.05) between the two groups. In addition, there was also no significant difference in TP53 expression between groups 4 and 5(p > 0.05). Groups 1, 2 and 3 showed the least TP53 protein expression, the values of which were not significantly different among them (p > 0.05) (Figures 4A, B). The 2 -ΔΔCT method was used to analyze the qRT-PCR data. TP53 mRNA expression was increased in groups 4, 5, 6 and 7 compared with that in the other three groups (Figure 4C). Groups 6 and 7 showed the highest TP53 mRNA expression, but there was no significant difference (p > 0.05) between the two groups, and there was also no significant difference between groups 4 and 5 (p > 0.05). Groups 1, 2 and 3 showed the lowest TP53 mRNA expression, which was not significantly different among them (p > 0.05) (Figure 4C).
Figure 4. (A) Statistical chart of TUNEL-DAPI staining, (B) Statistical chart of immunohistochemistry, (C) Statistical chart of qRT-PCR. (D) Section staining of osteosarcoma tissue for Perls’ Prussian blue and the TUNEL assay, (E) Statistical chart of Perls’ Prussian blue assay, (F) Statistical chart of the TUNEL assay. *, p < 0.05, compared with the Control group; NS = no significance; #, p < 0.05, compared with the STP-NBs group; δ, p < 0.05, compared with the STP-NBs + US + MT group.
The iron distribution was detected using Perls’ Prussian blue assay 7 days after transfection, and the results showed that the iron content was significantly higher in groups 2, 6 and 7 than in the other groups (p < 0.05); furthermore, these three groups exhibited no significant differences (p > 0.05) (Figures 4D, E). In addition, group 5 showed the highest iron content among all the treatment groups. There was also no significant difference among groups 1, 3, and 4(p > 0.05) (Figures 4D, E).
The TUNEL assay was used to measure OS apoptosis. The TUNEL-DAPI assay showed significantly increased apoptosis in groups 4, 5, 6 and 7 and significantly higher apoptosis than the other three groups (p < 0.05) (Figures 4D, F). Group 7 showed the highest number of TUNEL-negative cells, followed by group 6 (p < 0.05). The degree of apoptosis in groups 4 and 5 was not significantly different (p > 0.05) (Figures 4D, F). A few TUNEL-negative cells were found in groups 1/2/3 and showed no significance (p > 0.05) (Figures 4D, F). So these results demonstrated that ultrasound-targeted transfection in magnetic field enhanced the transfection efficiency of the TP53 plasmid and the overexpression of TP53 in local area, resulting in apoptosis on the osteosarcoma model in vivo.
Discussion
Current clinical treatment modalities for OS include surgical resection of primary tumors and systemic chemotherapy. (Panez-Toro et al., 2023). In recent years, the therapeutic use of exogenous genes for treating OS has gained much interest. TP53 is one of the most well-studied tumor suppressor genes. Previous studies have shown that gene transduction of functional TP53 can effectively inhibit OS cell activity. The purpose of gene therapy is to correct abnormal protein expression at the gene level, and the introduction of exogenous TP53 to correct the reduction of TP53 level is a promising OS therapy method. (Tsuchiya et al., 2000; Ganjavi et al., 2006). The mechanism of TP53 is well understood in OS. (Song et al., 2015).
Magnetic nanoparticles could become magnetized and exhibit physical sensitivity to external magnetic fields. Thus, it has been considered to lead to an increase in the local drug concentration at the target site. It was initially used merely for enhancing magnetic resonance imaging in most of the current studies. (Kiwada et al., 1986; Liu et al., 2009; Ting et al., 2012). However, these characteristics signify that magnetic nanoparticles are one kind of most suitable material for gene therapies.
Ultrasound irradiation is an ideal physical gene transfection method to promote extracellular molecules to enter cells. (Feril, 2009). The LIFU device used in this study was developed and patented by our institute. Current studies have indicated that RLI-mode LIFU can destroy target nanobubbles, which not only does not destroy exogenous genes but also improves transfection efficiency; meanwhile, increasing the intensity (RHI-mode) of the LIFU device would exert range-controllable cytocidal effects on tissues. (Xu et al., 2009; Song et al., 2017). To date, limited research on the gene therapy of OS by ultrasound-magnetic-mediated NB destruction has been conducted, and OS treatment using RHI-mode LIFU directly has not been reported.
In an in vitro study, the results provide evidence of the successful preparation of STP-NBs. They had a mean diameter of 243.6 ± 46.51 nm, a PDI of 0.284, and a zeta potential of −5.92 ± 3.37 mV. The optimal total input of SPIO and TP53 plasmids was 0.5 mg and 0.6 mg, respectively. The STP-NBs showed good sustained release of genes. In addition, the physical sensitivity of STP-NBs to external magnetic fields has been proven in vitro.
In the cell experiments, we loaded STP-NBs with TP53 and administered RLI-mode LIFU to treat MG63 and HOS cells. TP53 protein and mRNA expression in cells were not significantly different between the TP-NBs + US and STP-NBs + US groups. The CCK8 results showed that cell activity was decreased in the TP-NBs + US and STP-NBs + US groups, but there was no significant difference between them. The TP53 + SP-NBs group was not significantly different from the control group for the TP53 + US group. Furthermore, it was also shown that STP-NBs have no cytotoxicity. These results suggest that STP-NBs were able to transfect MG63 cells and HOS cells. The TP53 plasmid was well expressed and inhibited cell proliferation.
In the animal experiments, we used both STP-NBs loaded with TP53 and RHL-mode LIFU to treat tumor-bearing nude mice. The Perls’ Prussian blue assay showed that the iron content was significantly higher in groups 2, 6 and 7 than in the other groups, but no significant difference was observed among these three groups. This suggests that STP-NBs can be concentrated in the targeting area in response to external magnetic fields in vivo.
The protein and mRNA expression levels of TP53 were increased in groups 4, 5, 6 and 7 compared with the other three groups. Groups 6 and 7 showed the highest TP53 expression, but there was no significant difference between these two groups; in addition, there was also no significant difference between groups 4 and 5. Groups 1, two and three showed the least TP53 expression, and there was no significant difference between them. These data indicate that STP-NBs were able to mediate TP53 transfection and expression, and the external magnetic field improved the transfection efficiency of STP-NBs in vivo.
The TUNEL assay showed significantly increased apoptosis in groups 4, 5, 6 and 7 compared with that of the other three groups. Group 7 showed the highest number of giant, disorganized, TUNEL-negative structures with large clusters of nuclei; group 6 had the next highest map. The degree of apoptosis in groups 4 and 5 was not significantly different. A few TUNEL-negative cells were found in groups 1, 2 and 3, and there was no observed statistical significance. The results demonstrated that following transfection and expression of the TP53 gene, the rate of apoptosis in OS tissue was increased in vivo. Meanwhile, this effect further progressed the addition of external magnetic fields and RHI-mode LIFU in vivo to our testing ro.
When interpreting these findings, several limitations must be considered. The best guidance method is an external focusing magnetic field located in the center of the tumor body. However, this study was limited by the available equipment, and all magnetic field guidance was generated by the RuFe boron magnet on the surface of the tumor. The strength of the magnetic field in different parts of the tumor was varied. A convergent magnetic field device that can accurately locate the focus will improve future studies of gene therapy. SPIO could enhance magnetic resonance imaging, which could potentially better demonstrate the concentration of nanoparticles in animals. (Fan et al., 2016). Due to equipment limitations, that study could not proceed. Furthermore, the LIFU waves used in this study may have released TP53 into other organs, causing additional biological effects. Previous studies only demonstrated the feasibility of the combination of LIFU and ultrasound-magnetic-mediated Spio-P53-nanobubble destruction at the cell and organism levels. (Elumalai et al., 2024). As a kind of high capacity of reflection biological tissues, less ultrasound energy is able to penetrate the inner portion of the bony tissue. Although high-intensity focused ultrasound (HIFU) experiments have proven that sound attenuation decreases and the sound beam easily passes through because of OS bone destruction, (Wu et al., 2001), the ossifying tumor of the human body needs further study. However, despite these limitations, this does not prevent us from trying some new research methods to treat OS more effectively or to exlpore new method to prepare for surgery.
Conclusion
We used LIFU combined with US-magnetic-mediated SPIO-TP53-NB destruction for OS treatment. Our findings indicate that the destruction of US-magnetic-mediated SPIO/P53/PLGA/NBs facilitates gene transfection and that LIFU combined with the destruction of these NBs may represent a novel, noninvasive and targeted therapy for OS.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
HR: Conceptualization, Writing–original draft. SX: Conceptualization, Writing–original draft. AL: Conceptualization, Writing–original draft. QW: Supervision, Writing–original draft. NZ: Supervision, Writing–review and editing. ZH: Supervision, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work has been financially supported by the National Natural Science Foundation of China (Nos. 81372003; 82102556). The Science and Technology Development Fund of Shanghai Pudong New Area (Grant No. PKJ2022-Y48), the Outstanding Leaders Training Program of Pudong Hospital affiliated to Fudan University (Grant No. LX202201), Program for the Outstanding Clinical Discipline Project of Shanghai Pudong (Grant No. PWYgy 2021-04), the Project of Key Medical Specialty and Treatment Center of Pudong Hospital of Fudan University (Zdzk 2020-02, Zdzk 2021-01).
Acknowledgments
The authors are grateful to American Journal Experts (AJE) for their help in language editing.
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.
References
Asano, Y., Meguro, R., Odagiri, S., Li, C., Iwatsuki, H., and Shoumura, K. (2006). Visualization of non-heme ferric and ferrous iron by highly sensitive non-heme iron histochemistry in the stress-induced acute gastric lesions in the rat. Histochem Cell Biol. 125, 515–525. doi:10.1007/s00418-005-0097-6
Bekeredjian, R., Kuecherer, H. F., Kroll, R. D., Katus, H. A., and Hardt, S. E. (2007). Ultrasound-targeted microbubble destruction augments protein delivery into testes. Urology 69, 386–389. doi:10.1016/j.urology.2006.12.004
Chen, S., Shimoda, M., Chen, J., Matsumoto, S., and Grayburn, P. A. (2012). Transient overexpression of cyclin D2/CDK4/GLP1 genes induces proliferation and differentiation of adult pancreatic progenitors and mediates islet regeneration. Cell Cycle 11, 695–705. doi:10.4161/cc.11.4.19120
Chen, Z. Y., Liang, K., and Qiu, R. X. (2010). Targeted gene delivery in tumor xenografts by the combination of ultrasound-targeted microbubble destruction and polyethylenimine to inhibit survivin gene expression and induce apoptosis. J. Exp. Clin. Cancer Res. 29, 152. doi:10.1186/1756-9966-29-152
Chertok, B., David, A. E., Huang, Y., and Yang, V. C. (2007). Glioma selectivity of magnetically targeted nanoparticles: a role of abnormal tumor hydrodynamics. J. Control Release 122, 315–323. doi:10.1016/j.jconrel.2007.05.030
Du, J., Shi, Q. S., Sun, Y., Liu, P. F., Zhu, M. J., Du, L. F., et al. (2011). Enhanced delivery of monomethoxypoly(ethylene glycol)-poly(lactic-co-glycolic acid)-poly l-lysine nanoparticles loading platelet-derived growth factor BB small interfering RNA by ultrasound and/or microbubbles to rat retinal pigment epithelium cells. J. Gene Med. 13, 312–323. doi:10.1002/jgm.1574
Elumalai, K., Srinivasan, S., and Shanmugam, A. (2024). Review of the efficacy of nanoparticle-based drug delivery systems for cancer treatment. Biomed. Technol. 5, 109–122. doi:10.1016/j.bmt.2023.09.001
Elzeny, H., Zhang, F., Ali, E. N., Fathi, H. A., Zhang, S., Li, R., et al. (2017). Polyphosphoester nanoparticles as biodegradable platform for delivery of multiple drugs and siRNA. Drug Des. Devel Ther. 22 (11), 483–496. doi:10.2147/DDDT.S128503
Fan, C. H., Cheng, Y. H., Ting, C. Y., Ho, Y. J., Hsu, P. H., Liu, H. L., et al. (2016). Ultrasound/magnetic targeting with SPIO-DOX-Microbubble complex for image-guided drug delivery in brain tumors. Theranostics 6, 1542–1556. doi:10.7150/thno.15297
Feril, L. B. (2009). Ultrasound-mediated gene transfection. Methods Mol. Biol. 542, 179–194. doi:10.1007/978-1-59745-561-9_10
Ganjavi, H., Gee, M., Narendran, A., Parkinson, N., Krishnamoorthy, M., Freedman, M. H., et al. (2006). Adenovirus-mediated p53 gene therapy in osteosarcoma cell lines: sensitization to cisplatin and doxorubicin. Cancer Gene Ther. 13, 415–419. doi:10.1038/sj.cgt.7700909
Hida, K., Lai, S. K., Suk, J. S., Won, S. Y., Boyle, M. P., and Hanes, J. (2011). Common gene therapy viral vectors do not efficiently penetrate sputum from cystic fibrosis patients. PLoS One 6, e19919. doi:10.1371/journal.pone.0019919
Huang, Q., Deng, J., Wang, F., Chen, S., Liu, Y., Wang, Z., et al. (2012). Targeted gene delivery to the mouse brain by MRI-guided focused ultrasound-induced blood-brain barrier disruption. Exp. Neurol. 233, 350–356. doi:10.1016/j.expneurol.2011.10.027
Hynynen, K., McDannold, N., Vykhodtseva, N., Raymond, S., Weissleder, R., Jolesz, F. A., et al. (2006). Focal disruption of the blood-brain barrier due to 260-kHz ultrasound bursts: a method for molecular imaging and targeted drug delivery. J. Neurosurg. 105, 445–454. doi:10.3171/jns.2006.105.3.445
Jiang, Z., and Dalby, P. A. (2023). Challenges in scaling up AAV-based gene therapy manufacturing. Trends Biotechnol. 41, 1268–1281. doi:10.1016/j.tibtech.2023.04.002
Kiwada, H., Sato, J., Yamada, S., and Kato, Y. (1986). Feasibility of magnetic liposomes as a targeting device for drugs. Chem. Pharm. Bull. (Tokyo) 34, 4253–4258. doi:10.1248/cpb.34.4253
Lin, Y., Lin, Y., Lin, X., Sun, X., and Luo, K. (2019). Combination of PET and CXCR4-targeted peptide molecule agents for noninvasive tumor monitoring. J. Cancer 10, 3420–3426. doi:10.7150/jca.31087
Liu, H. L., Hsu, P. H., Chu, P. C., Wai, Y. Y., Chen, J. C., Shen, C. R., et al. (2009). Magnetic resonance imaging enhanced by superparamagnetic iron oxide particles: usefulness for distinguishing between focused ultrasound-induced blood-brain barrier disruption and brain hemorrhage. J. Magn. Reson Imaging 29, 31–38. doi:10.1002/jmri.21599
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Lu, Q. L., Liang, H. D., Partridge, T., and Blomley, M. J. (2003). Microbubble ultrasound improves the efficiency of gene transduction in skeletal muscle in vivo with reduced tissue damage. Gene Ther. 10, 396–405. doi:10.1038/sj.gt.3301913
Mungarndee, S. S., Lundy, R. F. Jr, and Norgren, R. (2008). Expression of Fos during sham sucrose intake in rats with central gustatory lesions. Am. J. Physiol. Regul. Integr. Comp. Physiol. 295, R751–R763. doi:10.1152/ajpregu.90344.2008
Panez-Toro, I., Muñoz-García, J., Vargas-Franco, J. W., Renodon-Cornière, A., Heymann, M. F., Lézot, F., et al. (2023). Advances in osteosarcoma. Curr. Osteoporos. Rep. 21, 330–343. doi:10.1007/s11914-023-00803-9
Sirsi, S. R., Hernandez, S. L., Zielinski, L., Blomback, H., Koubaa, A., Synder, M., et al. (2012). Polyplex-microbubble hybrids for ultrasound-guided plasmid DNA delivery to solid tumors. J. Control Release 157, 224–234. doi:10.1016/j.jconrel.2011.09.071
Sitta, J., and Howard, C. M. (2021). Applications of ultrasound-mediated drug delivery and gene therapy. Int. J. Mol. Sci. 22, 11491. doi:10.3390/ijms222111491
Song, R., Tian, K., Wang, W., and Wang, L. (2015). P53 suppresses cell proliferation, metastasis, and angiogenesis of osteosarcoma through inhibition of the PI3K/AKT/mTOR pathway. Int. J. Surg. 20, 80–87. doi:10.1016/j.ijsu.2015.04.050
Song, W., Luo, Y., Zhao, Y., Liu, X., Zhao, J., Luo, J., et al. (2017). Magnetic nanobubbles with potential for targeted drug delivery and trimodal imaging in breast cancer: an in vitro study. Nanomedicine (Lond). 12, 991–1009. doi:10.2217/nnm-2017-0027
Steiniger, S. C., Kreuter, J., Khalansky, A. S., Skidan, I. N., Bobruskin, A. I., Smirnova, Z. S., et al. (2004). Chemotherapy of glioblastoma in rats using doxorubicin-loaded nanoparticles. Int. J. Cancer 109, 759–767. doi:10.1002/ijc.20048
Sun, Y., ZhangFang, C. Q., Zhang, W. Q., and Liu, W. (2023). Abnormal signal pathways and tumor heterogeneity in osteosarcoma. J. Transl. Med. 21, 99. doi:10.1186/s12967-023-03961-7
Tang, H. L., Wang, Z. G., Li, Q., Ran, H. T., Zheng, Y. Y., Ren, J. L., et al. (2012). Targeted delivery of bone mesenchymal stem cells by ultrasound destruction of microbubbles promotes kidney recovery in acute kidney injury. Ultrasound Med. Biol. 38 (4), 661–669. doi:10.1016/j.ultrasmedbio.2012.01.003
Ting, C. Y., Fan, C. H., Liu, H. L., Huang, C. Y., Hsieh, H. Y., Yen, T. C., et al. (2012). Concurrent blood-brain barrier opening and local drug delivery using drug-carrying microbubbles and focused ultrasound for brain glioma treatment. Biomaterials 33, 704–712. doi:10.1016/j.biomaterials.2011.09.096
Tsuchiya, H., Mori, Y., Ueda, Y., Okada, G., and Tomita, K. (2000). Sensitization and caffeine potentiation of cisplatin cytotoxicity resulting from introduction of wild-type p53 gene in human osteosarcoma. Anticancer Res. 20 (1A), 235–242.
Wang, Y., Yang, L., Chen, T., Liu, X., Guo, Y., Zhu, Q., et al. (2019). A novel lncRNA MCM3AP-AS1 promotes the growth of hepatocellular carcinoma by targeting miR-194-5p/FOXA1 axis. Mol. Cancer 18, 28. doi:10.1186/s12943-019-0957-7
Wells, D. J. (2010). Electroporation and ultrasound enhanced non-viral gene delivery in vitro and in vivo. Cell Biol. Toxicol. 26, 21–28. doi:10.1007/s10565-009-9144-8
Wu, F., Chen, W. Z., Bai, J., Zou, J. Z., Wang, Z. L., Zhu, H., et al. (2001). Pathological changes in human malignant carcinoma treated with high-intensity focused ultrasound. Ultrasound Med. Biol. 27, 1099–1106. doi:10.1016/s0301-5629(01)00389-1
Wu, Y., Sun, T., Tang, J., Liu, Y., and Li, F. (2020). Ultrasound-targeted microbubble destruction enhances the antitumor efficacy of doxorubicin in a mouse hepatocellular carcinoma model. Ultrasound Med. Biol. 46, 679–689. doi:10.1016/j.ultrasmedbio.2019.09.017
Xenariou, S., Griesenbach, U., Liang, H. D., Zhu, J., Farley, R., Somerton, L., et al. (2007). Use of ultrasound to enhance nonviral lung gene transfer in vivo. Gene Ther. 14, 768–774. doi:10.1038/sj.gt.3302922
Xu, H. X., Lu, M. D., Xie, X. H., Xie, X. Y., Xu, Z. F., Chen, L. D., et al. (2009). Three-dimensional contrast-enhanced ultrasound of the liver: experience of 92 cases. Ultrasonics 49, 377–385. doi:10.1016/j.ultras.2008.10.009
Yamaguchi, K., Feril, L. B., Tachibana, K., Takahashi, A., Matsuo, M., Endo, H., et al. (2011). Ultrasound-mediated interferon β gene transfection inhibits growth of malignant melanoma. Biochem. Biophys. Res. Commun. 411, 137–142. doi:10.1016/j.bbrc.2011.06.115
Yin, H., Kanasty, R. L., Eltoukhy, A. A., Vegas, A. J., Dorkin, J. R., and Anderson, D. G. (2014). Non-viral vectors for gene-based therapy. Nat. Rev. Genet. 15, 541–555. doi:10.1038/nrg3763
Zhang, X., Zheng, Y., Wang, Z., Huang, S., Chen, Y., Jiang, W., et al. (2014). Methotrexate-loaded PLGA nanobubbles for ultrasound imaging and Synergistic Targeted therapy of residual tumor during HIFU ablation. Biomaterials 35, 5148–5161. doi:10.1016/j.biomaterials.2014.02.036
Zhao, Y., Song, W., Wang, D., Ran, H., Wang, R., Yao, Y., et al. (2015). Phase-shifted PFH@PLGA/Fe3O4 nanocapsules for MRI/US imaging and photothermal therapy with near-infrared irradiation. ACS Appl. Mater Interfaces 7, 14231–14242. doi:10.1021/acsami.5b01873
Zhong, S., Shu, S., Wang, Z., Luo, J., Zhong, W., Ran, H., et al. (2012). Enhanced homing of mesenchymal stem cells to the ischemic myocardium by ultrasound-targeted microbubble destruction. Ultrasonics 52, 281–286. doi:10.1016/j.ultras.2011.08.013
Glossary
Keywords: osteosarcoma, low-intensity-focused ultrasound, nanobubbles, superparamagnetic iron oxide, PLGA
Citation: Ren H, Xiang S, Liu A, Wang Q, Zhou N and Hu Z (2024) A noval noninvasive targeted therapy for osteosarcoma: the combination of LIFU and ultrasound-magnetic-mediated SPIO/TP53/PLGA nanobubble. Front. Bioeng. Biotechnol. 12:1418903. doi: 10.3389/fbioe.2024.1418903
Received: 17 April 2024; Accepted: 10 June 2024;
Published: 28 June 2024.
Edited by:
Wei Nie, Wake Forest Baptist Medical Center, United StatesReviewed by:
Dongxu Ke, Wake Forest University, United StatesTao Ji, Shenzhen University, China
Xijian Liu, Shanghai University of Engineering Sciences, China
Copyright © 2024 Ren, Xiang, Liu, Wang, Zhou and Hu. 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: Qian Wang, Drwangqian23@163.com; Nian Zhou, zhounian_2008@126.com; Zhenming Hu, spinecenter@163.com
†These authors have contributed equally to this work