- 1Department of Anatomy, School of Chinese Medicine, Beijing University of Chinese Medicine, Beijing, China
- 2Dongfang Hospital, Beijing University of Chinese Medicine, Beijing, China
- 3Center for Scientific Research, School of Chinese Medicine, Beijing University of Chinese Medicine, Beijing, China
As a typical traditional Chinese medicine, Bu-Yin-Qian-Zheng Formula (BYQZF) has been shown to have neuroprotective effects in patients with Parkinson’s disease (PD), particularly by ameliorating mitochondrial dysfunction and regulating expression of the parkin protein. However, the underlying mechanisms by which BYQZF affects mitochondrial function through parkin are unclear. Accordingly, in this study, we evaluated the mechanisms by which BYQZF ameliorates mitochondrial dysfunction through parkin in PD. We constructed a parkin-knockdown cell model and performed fluorescence microscopy to observe transfected SH-SY5Y cells. Quantitative real-time reverse transcription polymerase chain reaction and western blotting were conducted to detect the mRNA and protein expression levels of parkin. Additionally, we evaluated the cell survival rates, ATP levels, mitochondrial membrane potential (ΔΨm), mitochondrial morphology, parkin protein expression, PINK1 protein expression, and mitochondrial fusion and fission protein expression after treatment with MPP+ and BYQZF. Our results showed that cell survival rates, ATP levels, ΔΨm, mitochondrial morphology, parkin protein levels, PINK1 protein levels, and mitochondrial fusion protein levels were reduced after MPP+ treatment. In contrast, mitochondrial fission protein levels were increased after MPP+ treatment. Moreover, after transient transfection with a negative control plasmid, the above indices were significantly increased by BYQZF. However, there were no obvious differences in these indices after transient transfection with a parkin-knockdown plasmid. Our findings suggest that BYQZF has protective effects on mitochondrial function in MPP+-induced SH-SY5Y cells via parkin-dependent regulation of mitochondrial dynamics.
Introduction
Parkinson’s disease (PD) is a common progressive neurodegenerative disease manifested by motor and nonmotor symptoms (Kalia and Lang, 2015) with degeneration of dopaminergic neurons in the substantia nigra pars compacta (Song et al., 2010; Ham et al., 2020). PD mainly occurs in people over 60 years of age (Grünewald et al., 2019) and is a major burden on social and economic development, particularly as society ages. Therefore, the pathogenesis of PD has been widely studied (Golpich et al., 2017; Obeso et al., 2017; Chen et al., 2018). Mitochondrial dysfunction has been shown to play a central role in the pathogenesis of PD (Exner et al., 2012). Mitochondria are distributed in a highly dynamic tubular network that is constantly reshaped by fusion and fission. Changes in fusion or fission lead to mitochondrial fragmentation, limiting mitochondrial motility and decreasing ATP production, thereby resulting in mitochondrial dysfunction (Luo et al., 2015; Rao et al., 2020). Additionally, the gene encoding parkin, which is related to mitochondrial function and mitophagy, is a genetic susceptibility gene in patients with PD (Kitada et al., 1998). Mutation or deletion of the gene encoding parkin can trigger accumulation of dysfunctional mitochondria (Bankapalli et al., 2020) and cause the death of dopaminergic neurons in PD (Zhu et al., 2018). Because of the substantial role of mitochondrial dysfunction in PD, mitochondrion-targeted therapeutics may facilitate the development of novel drugs for treating PD (Camilleri and Vassallo, 2014; Park et al., 2018; Grünewald et al., 2019; Bai et al., 2020), and parkin may be an important mitochondrial target for expanding our understanding of the molecular mechanisms of mitochondrial dysfunction in PD therapy.
Currently, symptomatic treatment with drugs that stimulate dopamine receptors is used to treat PD (Kalia and Lang, 2015); however, these drugs have major side effects (Kordower, 2015). To reduce and counteract these effects, traditional Chinese medicine has been studied and applied (Chua et al., 2017; Liu et al., 2020). The traditional Chinese medicine Bu-Yin-Qian-Zheng formula (BYQZF) has been widely used to treat PD and its complications. BYQZF is composed of Da-Bu-Yin-Wan (Great Yin Tonic Pill) and Qian-Zheng-San (Symmetry Leading Powder). Our previous studies showed that Da-Bu-Yin-Wan ameliorates MPTP-induced behavioral impairment, increases the expression of tyrosine hydroxylase in the substantia nigra, and prevents loss of substantia nigra dopaminergic neurons in a PD mouse model (Gong et al., 2014). Moreover, Da-Bu-Yin-Wan and/or Qian-Zheng-San significantly increase ATP levels and complex I activity in the midbrain tissue of PD model mice (Zhang et al., 2013), and BYQZF reduced mitochondrial damage and improved the mitochondrial morphological structure in a cellular model of PD (Zhang et al., 2018; Gai et al., 2019). These studies suggest that BYQZF had neuroprotective effects on mitochondrial function and may be applicable as a mitochondrion-targeted therapeutic drug for treating PD. However, the mechanisms by which BYQZF ameliorates mitochondrial dysfunction via parkin are unclear.
Accordingly, in this study, we evaluated the potential mechanisms by which BYQZF ameliorated mitochondrial dysfunction through parkin in MPP+-induced SH-SY5Y cells. Our findings provide important insights into the roles of BYQZF and parkin in MPP+-induced SH-SY5Y cells.
Materials and Methods
Chemical Teagents and Antibodies
SH-SY5Y human neuroblastoma cells were acquired from Shanghai GeneChem Co., Ltd. 1-Methyl-4-phenylpyridinium (MPP+) was obtained from Sigma Corporation (St. Louis, MO, United States), and a Cell Counting Kit-8 (CCK-8) was obtained from Dojindo Molecular Technologies. Lipofectamine 3,000, MitoTracker® Red CMXRos (MTR), and tetramethylrhodamine were obtained from Life Technologies and the ATP detection reagent kit was obtained from Beyotime Biotechnology Co., Ltd. The BCA Protein Assay Kit and enhanced chemiluminescence (ECL) autoradiography kit were obtained from Applygen Technologies, Inc.
Anti-mitofusin (MFN) 1, anti-dynamin-related protein 1 (DRP1), and anti-optic atrophy 1 (OPA1) antibodies were obtained from Abcam (Cambridge, United Kingdom). Anti-parkin, anti-MFN2, anti-fission mitochondrial 1 (FIS1), and anti-β-actin polyclonal antibodies were obtained from Proteintech. All secondary antibodies (1:2000) were obtained from Zhong-Shan Golden Bridge Biotechnology.
Preparation and Analysis for Bu-Yin-Qian-Zheng Formula
BYQZF was composed of Rehmannia glutinosa (Gaertn.) DC. (Shu-Di-Huang), Phellodendron chinense C.K.Schneid. (Huang-Bai), Anemarrhena asphodeloides Bunge (Zhi-Mu), Plastrum Testudinis (Gui-Ban),Sauromatum giganteum (Engl.) Cusimano & Hett. (Bai-Fu-Zi), Bombyx Batryticatus (Jiang-Can), and Scorpio (Quan-Xie)(Chan et al., 2012). All traditional Chinese medicine herbs were purchased from Tong-Ren-Tang Drugstore (Beijing, China) and authenticated by experts in pharmacognosy, as shown in Table 1. To prepare the decoctions, exact amounts of herbs were weighed according to classic percentages and mixed well. The mixture was soaked in distilled water for 30 min and then boiled in eight volumes of water (v/w) for 1 h; this extraction procedure was carried out twice (Zhang et al., 2013). Finally, the decoction was concentrated to 0.5 g/ml, filtered through a 0.22-μm membrane, divided, and stored as a stock solution at –70°C.
Identification and quantification of the marker compounds in BYQZF were performed according to the method described in the updated Pharmacopeia of People’s Republic of China (Chinese Pharmacopoeia Commission, 2020). Briefly, the marker compounds were analyzed using the ultra-performance liquid chromatography tandem–mass spectrometry (UPLC-MS/MS). The separation was carried out on an Acquity BEH C18 (100 mm × 2.1 mm, 1.7 μm; Waters, Milford, MA, United States) analytical column. The column was maintained at 35°C. Water acidified with acetonitrile as Solvent A and 0.1% of formic acid (v/v) as Solvent B were used for the gradient elution at 0.3 ml/min. The gradient program was: 5–50% A (from 0 to 4.5 min), 50% A (from 4.5 to 6 min), 50–100% A (from 6 to 8 min), 100 to 5% A (from 8 to 8.1 min), 5% A (from 8.1 to 10 min) and injection volume of 2 µl. The experiment was performed on both the positive and negative ion modes. Ion spray voltage: positive ion modes, 5 kV; negative ion modes, −2.8 kV; lens voltage: positive ion modes, 50 V; negative ion modes, −50 V; capillary temperature, 320°C; auxiliary gas temperature, 350°C; sheath gas (N2) voltage, 40 arb; auxiliary gas (N2) voltage, 10 arb. Full scans were acquired in the range of 200–1,200 m/z with a resolution of 70,000 in Full MS and the resolution of 17,500 in dd-MS2. Fragmentation amplitude was set at 20, 40, 60 eV.
Cell Culture, Transfection, and Treatments
SH-SY5Y cells were cultured in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin in a humidified incubator (37°C, 5% CO2). For differentiation experiments, cells were replated into 6-, 12-, or 96-well plates and grown to 70–80% confluence after 24 h.
The human wild-type parkin knockdown plasmid was designed by GeneChem. SH-SY5Y cells were plated at a density of 1 × 105 cells/mL and transfected with parkin-knockdown plasmid using the optimum amount of Lipofectamine 3,000 (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s protocols. After transfection for 48 h, the expression of green fluorescent protein (GFP) was observed under a fluorescence microscope. Cells were cultured in medium containing MPP+ (1 mM) with or without BYQZF (5 μg/ml) for 48 h.
Experimental Groups
The experiment was divided into seven groups: Control group (C), Negative control group (NC), Negative model group (NM), Negative treatment group (NT), Parkin-knockdown control group (PKC), Parkin-knockdown model group (PKM) and Parkin-knockdown treatment group (PKT). Experimental groups and treatments are shown in Table 2.
Quantitative Real-Time Reverse Transcription Polymerase Chain Reaction
The mRNA expression of parkin was verified by qRT-PCR using a One-Step SYBR PrimeScript RT-PCR Kit II and quantified on an ABI StepOne Plus sequence detection system under the following conditions: 5 min at 42°C and 10 s at 95.0°C for initial denaturation, followed by 39 cycles of 5 s at 95°C and 20 s at 60°C. The unknown template was calculated using the 2−ΔΔCT method for quantitative analysis. Sequences of the qRT-PCR-specific primers are listed in Table 3.
Western Blotting
Total protein samples were obtained from SH-SY5Y cells using a BCA Protein Assay kit. First, equal amounts of protein were loaded into the wells and separated by electrophoresis. The proteins were then transferred to polyvinylidene difluoride membranes (Millipore, Billerica, MA, United States), and nonspecific binding was blocked for 1 h with skim milk. The membranes were then incubated with primary antibodies overnight at 4°C followed by incubation at room temperature with horseradish peroxidase-conjugated goat anti-rabbit/mouse IgG (H + L) secondary antibodies. Finally, protein bands were visualized using an ECL autoradiography kit.
Assessment of Cell Survival Rate
CCK-8 assays were performed to evaluate cell survival rates (Ishiyama et al., 1997). Briefly, 10 μl CCK-8 reagent and 90 μl Dulbecco’s modified Eagle’s medium were added to each well of a 96-well plate. After incubation for 1 h (37°C, 5% CO2), the absorbance at 450 nm was measured using a full-wavelength fluorescence microplate (Tecan Safire2; Tecan, Switzerland). The survival rate was calculated as follows:
where
Mitochondrial Membrane Potential (ΔΨm)
Mitochondrial staining was performed using the mitochondrial probe tetramethylrhodamine in living cells to detect the mitochondrial membrane potential followed by observation with an Olympus FV1000 confocal microscope (559 nm excitation/572 nm emission). ImageJ software (Schneider et al., 2012) was used to process the images and measure the mitochondrial fluorescence intensity (OD value) of each image.
Adenosine Triphosphate Levels
Cellular ATP content was detected using luciferase assays following the manufacturer’s instructions. Briefly, ATP concentrations in the samples were calculated based on a standard curve, and the protein concentration of the sample was detected using a BCA Protein Concentration Assay Kit. Finally, the ATP concentration was converted into nmol/mg protein.
Confocal Fluorescence Microscopy for Mitochondrial Morphology
Cells were seeded into a 20-mm confocal culture dish. Briefly, the cells were incubated with 100 nM MTR for 15 min at 37°C in the dark before visualization (Pendergrass et al., 2004). Fluorescence was detected (559 nm excitation/603 nm emission) using an Olympus FV1000 confocal microscope. Acquired images were quantified by ImageJ software. The form factor was calculated using the following equation: form factor = P2/4πA (P, perimeter; A, area), as previously described (Mortiboys et al., 2008). Mitochondrial fluorescence intensity (OD value) and the mean length and number of mitochondrial network branches were also measured using ImageJ software.
Statistical Analysis
The data are presented as the means ± standard deviations and were analyzed by Two-Two Comparisons in one-way analysis of variance (One-way ANOVA). LSD(L) and S-N-K(S) were performed to determine p values using SPSS20.0 software. Results with p values of less than 0.05 between two groups were considered as statistically significant.
Results
Analysis of Marker Compounds of Bu-Yin-Qian-Zheng Formula
The marker compounds in BYQZF were analyzed with UPLC-MS/MS. Chromatograms of the BYQZF analyzed with relative reference standards are shown in Figure 1. Table 4 shows the equations of the calibration curves and limits of quantification of those components determined. All calibration curves showed good linear regression (r2 > 0.99). Table 5 shows the results of precision and accuracy tests. Calculated by the calibration curves using an internal standard, it is clearly shown that the concentrations of the six compounds calculated by the calibration curves.
FIGURE 1. UPLC-MS/MS analysis of BYQZF. Six reference standards were used to identify the marker compounds for BYQZF. (A) Chromatograms of BYQZF. (B) Secondary debris ion of BYQZF.
TABLE 4. Linear range, r2 and limits of quantification of calibration curve used in the determination of the six analytes.
Expression of Green Fluorescent Protein in Transfected Cells
SH-SY5Y cells were transfected with the negative control plasmid and parkin-knockdown plasmid, both of which were labeled with GFP. The expression of GFP was observed by fluorescence microscopy at 48 h after transfection.
As shown in Figure 2, the number of cells was high, and the cells showed a plump morphology in both the control and negative control groups. In contrast, in the parkin-knockdown group, the number of cells was reduced, and the cells showed a generally normal morphology. Increased GFP expression was observed in the negative control group and parkin-knockdown control group (Figure 2), suggesting that the plasmid was successfully transfected into the cells.
FIGURE 2. Expression of GFP in SH-SY5Y cells. GFP expression was observed by fluorescence microscopy at 48 h after transfection in SH-SY5Y cells. PKC, parkin-knockdown control group; NC, negative control group.
Analysis of mRNA and Protein Expression Levels of Parkin
Next, we investigated parkin expression at 48 h after transfection by qRT-PCR and western blotting. As shown in Figures 3Ac, parkin mRNA expression was decreased in the parkin-knockdown control group (p < 0.05). There were no differences between the control and negative control groups (p > 0.05). Additionally, melting curve analysis revealed that no primer-dimers were formed during the 40 real-time PCR amplification cycles (Figures 3Aa,b).
FIGURE 3. mRNA and protein expressions of parkin in SH-SY5Y cells. Parkin expression was evaluated at 48 h after transfection in SH-SY5Y cells. (A) mRNA expression of parkin after transient transfection with a parkin-knockdown plasmid. (A) Amplification cycles. (B) Melt curve analysis. (C) Fold-change in parkin mRNA expression, as quantified using the 2−ΔΔCT method. (B) Protein expression of parkin after transient transfection with parkin-knockdown plasmid. (A) Western blotting bands. (B) Quantification of parkin expression; the bar graph shows the data in arbitrary units (AU) from three independent experiments after normalization to β-actin expression. Data are expressed as means ± standard deviations (n = 3). Compared with the control (C) group, ∗p < 0.05; compared with the negative control (NC) group, #p < 0.05.
As shown in Figure 3B, the expression of parkin protein was visibly decreased in the parkin-knockdown control group (p < 0.05). There were no differences between the control and negative control groups (p > 0.05). These findings indicate that the parkin-knockdown SH-SY5Y cell line was successfully constructed.
Effects of Bu-Yin-Qian-Zheng Formula on Cell Survival Rates, Adenosine Triphosphate Levels, and ΔΨm in Parkin-Knockdown Parkinson’s Disease Cells
Parkin, which is a common genetic susceptibility gene in patients with PD, plays key roles in cell death induced by mitochondrial dysfunction (Bai et al., 2020). Therefore, we evaluated the effects of BYQZF on mitochondrial function in PD cells after transfection with a parkin-knockdown plasmid. After treatment with MPP+, the cell survival rate (Figure 4A), ATP level (Figure 4B), and relative fluorescence intensity of mitochondria (Figure 4C) were obviously decreased in each model group (p < 0.01). Additionally, the cell survival rate, ATP level, and relative fluorescence intensity of mitochondria decreased significantly (p < 0.01) in parkin-knockdown SH-SY5Y cells. Notably, treatment with BYQZF markedly enhanced these parameters in the negative treatment group (p < 0.01; p < 0.05). However, the cell survival rate, ATP level, and relative fluorescence intensity of mitochondria did not increase after treatment with BYQZF in the parkin-knockdown model group (p > 0.05), suggesting that the effects of BYQZF on mitochondrial function were significantly lower after transfection with the parkin-knockdown plasmid.
FIGURE 4. Analysis of cell survival rates and mitochondrial function. (A) Effects of BYQZF on the survival rate of PD cells after transfection with a parkin-knockdown plasmid. (B) Effects of BYQZF on ATP levels in PD cells after transfection with a parkin-knockdown plasmid. (C) Effects of BYQZF on the ΔΨm of PD cells after transfection with a parkin-knockdown plasmid. (A) The mitochondria of living cells were labeled with tetramethylrhodamine. (B) Effects of BYQZF on the relative fluorescence intensity of mitochondria in PD cells after with transfection with a parkin-knockdown plasmid. Data are expressed as the means ± standard deviations (n = 3). Compared with the negative control (NC) group, **p < 0.01; compared with the negative model (NM) group, #p < 0.05, ##p < 0.01; compared with the negative treatment (NT) group, △p < 0.05, △△p < 0.01; compared with the parkin-knockdown control (PKC) group, ▲p < 0.05, ▲▲p < 0.01.
Effects of Bu-Yin-Qian-Zheng Formula on Mitochondrial Morphology in Parkin-Knockdown Parkinson’s Disease Cells
Mitochondrial morphology is highly variable in different biological environments. Mitochondria are fragmented and spherical and show a lack of dynamics and network formation when damaged. In this study, the probe MTR was used to detect mitochondrial morphology. As shown in Figure 5A, the cells were full of clear nuclei, and mitochondria were distributed around the nucleus with long protrusions in the negative control group. Compared with the negative control group, the number of cells was reduced and mitochondria were shorter and rounder in each model group and in the parkin-knockdown control group. BYQZF alleviated these effects in the negative treatment group. However, there were no significant changes in mitochondrial morphology between the parkin-knockdown treatment group and parkin-knockdown model group.
FIGURE 5. Effects of BYQZF on the mitochondrial morphology of PD cells after transfection with a parkin-knockdown plasmid. (A) Mitochondria were labeled with MTR and assessed by confocal fluorescence microscopy (100×). (B) Form factor of mitochondria. (C) Mean length of mitochondrial network branches. (D) Mean number of mitochondrial network branches. (E) Analysis of mitochondrial activity. Data are expressed as the means ± standard deviations (n = 3). Compared with the negative control (NC) group, *p < 0.05, **p < 0.01; compared with the negative model (NM) group, #p < 0.05, ##p < 0.01; compared with the negative treatment (NT) group, △p < 0.05, △△p < 0.01; compared with the parkin-knockdown control (PKC) group, ▲▲p < 0.01.
As shown in Figures 5B–E, the mitochondrial form factor, mean length and number of mitochondrial network branches, and mitochondrial activity were significantly reduced in each model group after treatment with MPP+ (p < 0.01). Moreover, the form factor and mitochondrial activity decreased significantly (p < 0.01) in parkin-knockdown SH-SY5Y cells; the mean length and number of mitochondrial network branches also decreased; however, these differences were not significant (p > 0.05). After treatment with BYQZF, the mitochondrial form factor, mean length and number of mitochondrial network branches, and mitochondrial activity were significantly increased in the negative treatment group (p < 0.01). Compared with the parkin-knockdown model group, there were no obvious differences in any of these parameters in the parkin-knockdown treatment group (p > 0.05).
Effects of Bu-Yin-Qian-Zheng Formula on Expression of PINK1 and Parkin Protein in Parkin-Knockdown Parkinson’s Disease Cells
As shown in Figure 6, the western blotting results showed that the expressions of parkin and PINK1 protein in each model group decreased noticeably after MPP+ treatment (p < 0.01). In addition, the expression of parkin protein was significantly reduced in parkin-knockdown SH-SY5Y cells (p < 0.01). BYQZF increased the expression of PINK1 and parkin protein in the negative treatment group (p < 0.01; Figures 6B,C). However, the expression of parkin protein did not increase after treatment with BYQZF in the parkin-knockdown model group (p > 0.05; Figure 6C).
FIGURE 6. Effects of BYQZF on the expression of PINK1 and parkin proteins in PD cells after transfection with a parkin-knockdown plasmid. (A) Western blotting to detect PINK1 and parkin proteins expression. (B) Quantification of PINK1 protein levels. (C) Quantification of parkin protein levels. Data are expressed as the means ± standard deviations (n = 3). Compared with the negative control (NC) group, **p < 0.01; compared with the negative model (NM) group, ##p < 0.01; compared with the negative treatment (NT) group, △△p < 0.01; compared with the parkin-knockdown control (PKC) group, ▲▲p < 0.01.
Effects of Bu-Yin-Qian-Zheng Formula on Expression of Mitochondrial Fusion and Fission Proteins in Parkin-Knockdown Parkinson’s Disease Cells
Mitochondrial dynamics are critical for regulating the mitochondrial morphology, number, and function (Meyer et al., 2017). Therefore, we next investigated the effects of BYQZF on mitochondrial dynamics in PD cells. Western blotting results showed that the expression levels of mitochondrial fusion proteins (MFN1, MFN2, and OPA1) were significantly reduced in each model group after MPP+ treatment (p < 0.01, p < 0.05; Figure 7A). The expression of MFN2 protein decreased in parkin-knockdown SH-SY5Y cells (p < 0.05). Additionally, BYQZF increased the expression of mitochondrial fusion proteins (MFN1, MFN2, and OPA1) in the negative treatment group (p < 0.01, p < 0.05) and expression of OPA1 protein in the parkin-knockdown model group (p < 0.05). However, the expression of mitochondrial fusion proteins (MFN1 and MFN2) did not increase after treatment with BYQZF in the parkin-knockdown model group (p > 0.05).
FIGURE 7. Analysis of expression of mitochondrial fusion and fission proteins. (A) Effects of BYQZF on the expression of mitochondrial fusion proteins after transfection with a parkin-knockdown plasmid. (a) Western blotting to detect MFN1, MFN2, and OPA1. Quantification of MFN1 (b), MFN2 (c), and OPA1 (d) levels. (B) Effects of BYQZF on expression of mitochondrial fission proteins after transfection with a parkin-knockdown plasmid. (a) Western blotting to detect DRP1 and FIS1. Quantification of DRP1 (b) and FIS1 (c) levels. Data are expressed as the means ± standard deviations (n = 3). Compared with the negative control (NC) group, *p < 0.05, **p < 0.01; compared with the negative model (NM) group, #p < 0.05, ##p < 0.01; compared with the negative treatment (NT) group, △p < 0.05, △△p < 0.01; compared with the parkin-knockdown control (PKC) group, ▲p < 0.05, ▲▲p < 0.01; compared with the parkin-knockdown model (PKM) group, ☆p < 0.05.
As shown in Figure 7B, MPP+ increased the expression of mitochondrial fission proteins (DRP1 and FIS1) in each model group (p < 0.01). After parkin knockdown, the expression of mitochondrial fission proteins (DRP1 and FIS1) increased significantly (p < 0.01, p < 0.05). Moreover, after BYQZF treatment, the expression of DRP1 and FIS1 was decreased in the negative treatment group (p < 0.01, p < 0.05). Compared with the parkin-knockdown model group, there were no significant changes in the expression levels of DRP1 and FIS1 in the parkin-knockdown treatment group (p > 0.05).
Discussion
The protective effects of BYQZF on mitochondrial function have been confirmed in in vitro and in vivo studies (Zhang et al., 2013; Zhang et al., 2018; Gai et al., 2019); parkin was found to play an important role in this process. However, the mechanisms by which BYQZF ameliorates mitochondrial dysfunction via parkin are unclear. In the current study, BYQZF showed reduced protective effects on mitochondrial function when parkin was knocked down. This is the first study to demonstrate that BYQZF ameliorated mitochondrial dysfunction through parkin.
In this study, the SH-SY5Y cell line has been used as a model for PD since it has the machinery to synthesize DA (Krishna et al., 2014; Xicoy et al., 2017). MPP+ is one of the most frequently used neurotoxins in creating cellular models and studying critical aspects of PD in SH-SY5Y cells (Yang et al., 2018; Chanthammachat and Dharmasaroja, 2019). Therefore, MPP+-induced SH-SY5Y neurotoxicity has been broadly utilized to create cellular models and study the mechanisms of PD.
Mitochondria play key roles in energy modulation (Nunnari and Suomalainen, 2012). The maintenance of mitochondrial functionality is essential for neuronal survival in PD (Nunnari and Suomalainen, 2012; Schapira, 2012; Bankapalli et al., 2020) because MPP+ inhibits mitochondrial complex I and thus causes mitochondrial dysfunction (Bender et al., 2006; Rao et al., 2020). The electrochemical gradient (ΔΨm) across the inner mitochondrial membrane is important for ATP production (Finkel et al., 2015). These factors reflect the state of mitochondrial function. Moreover, several studies have shown that parkin plays a direct role in mitochondrial biogenesis (Sliter et al., 2018; Noda et al., 2020), and mutation or deletion of parkin can lead to changes in mitochondrial function and trigger PD-like symptoms (Dodson and Guo, 2007; Zhu et al., 2018). For example, Stevens et al. (2015) found that loss of parkin may damage mitochondrial health by impairing mitochondrial biogenesis and eventually lead to cell death. Additionally, Lim et al. (2012) found that impaired mitochondria may contribute to DA neurodegeneration in parkin-related cases. In this study, a parkin-knockdown cell model was successfully constructed, and the ΔΨm and ATP levels were lower in the parkin-knockdown control group. Taken together, these data suggest that parkin maintains the normal state of the mitochondria and can be used as a mitochondrion-targeted protein for applications in PD therapy (Plotegher and Duchen, 2017).
BYQZF is composed of Amur cork tree bark (Huangbai), prepared Rehmannia root (Shudihuang), common anemarrhena (Zhimu), tortoise plastron (Guiban), Chinese scorpion (detoxicated; Quanxie), stiff silkworm (Jiangcan), and giant typhonium tuber (Baifuzi). Recent studies have shown that some active ingredients in BYQZF can protect the activity of substantia nigra DA neurons. For example, catalpol, a component of shú dì huáng, alleviates MPTP-induced degeneration of DA neurons in vivo (Wang et al., 2019). Additionally, mangiferin, a chemical constituent of zhî mu, significantly upregulates substantia nigra dopamine concentrations and decreases the effects of neurotoxins, thereby affecting mitochondrial function in PD (Kavitha et al., 2013; Feng et al., 2019). Extracts from guî bãn promote the differentiation of neural stem cells into dopaminergic neurons by regulating OTX2 expression (Zhou et al., 2017). Notably, we previously found that BYQZF ameliorated mitochondrial dysfunction and promoted the expression of parkin protein (Gai et al., 2019). However, further studies are required to determine whether parkin is a mitochondrial target for BYQZF for treating PD.
We found that BYQZF rescued MPP+-induced mitochondrial alterations, including ΔΨm and ATP levels in the negative treatment group, consistent with a previous study (Gai et al., 2019). However, BYQZF did not rescue mitochondrial dysfunction in MPP+-induced SH-SY5Y cells transfected with the parkin-knockdown plasmid. Accordingly, parkin may be a mitochondrial target for BYQZF and can protect mitochondrial function in the treatment of PD. However, the exact mechanisms through which BYQZF affects mitochondrial functions are poorly understood, and identification of these mechanisms may facilitate the development of novel treatments of PD. Thus, further studies are required.
Mitochondrial morphology, which can be regulated by the balance between fusion and fission (Li et al., 2016), is intimately linked with mitochondrial function (Nunnari and Suomalainen, 2012) and facilitates specific cellular activities (Lisowski et al., 2018). Notably, parkin can maintain mitochondrial morphology and function by regulating the state of fission and fusion. Additionally, most studies in mammalian cells have shown that parkin acts in the direction of fusion and that parkin overexpression leads to the growth of interconnected mitochondrial networks (Tang et al., 2016). Outer mitochondrial membrane fusion is mediated by MFNs (MFN1 and MFN2), whereas inner membrane fusion is mediated by OPA1 (Wallace, 2012; Pernas and Scorrano, 2016). Mitochondrial fission is mediated by DRP1 and FIS1 (Schapira, 2012; Pernas and Scorrano, 2016). In this study, we found that the mitochondrial form factor, mean length and number of mitochondrial network branches, mitochondrial activity, and expression of MFN2 protein were decreased, whereas DRP1 and FIS1 expression levels were increased in the parkin-knockdown control group. These results are consistent with those of previous studies demonstrating that parkin downregulation exacerbates mitochondrial fragment via regulation of mitochondrial fission and fusion (Pinto and Nissanka, 2018; Noda et al., 2020). Overall, these findings highlight the importance of parkin not only as a genetic marker of PD but also as an important factor in its pathogenesis, via a mechanism related to the regulation of mitochondrial dynamics. The precise mechanisms through which BYQZF affects mitochondrial functions via parkin may related to the regulation of mitochondrial dynamics.
Parkin is involved in the balance between fusion and fission, functioning to maintain mitochondrial morphology. MFN1, MFN2, DRP1, and FIS1 are the main substrates of parkin (Kandimalla and Reddy, 2016). In this study, we found that parkin protein expression was downregulated in MPP+-induced SH-SY5Y cells; however, treatment with BYQZF blocked this effect, promoting parkin protein expression after transfection of cells with the negative plasmid. In contrast, BYQZF did not enhance parkin expression following transfection with a parkin-knockdown plasmid. Moreover, changes in mitochondrial morphology were consistent with parkin protein expression. Additionally, we found that changes in the expression of MFN1 and MFN2 were consistent with those in the expression of parkin protein, whereas changes in DRP1 and FIS1 expression were opposite those of parkin protein. These findings suggest that knockdown of parkin suppressed the regulatory effects of BYQZF on mitochondrial fusion and fission proteins. Accordingly, parkin may mediate the effects of BYQZF on mitochondrial function and biogenesis by regulating fusion and fission processes. In this study, we evaluated the potential mechanisms by which BYQZF ameliorated mitochondrial dysfunction by regulating mitochondrial dynamics via parkin in PD. However, some hypotheses regarding the mechanism by which parkin mutations or loss of function lead to PD have suggested that parkin deficiency is associated with mitochondrial autophagy. We have not yet investigated the effect of BYQZF on mitochondrial autophagy. Whether the mechanism through which BYQZF affects mitochondrial functions via parkin is related to the regulation of mitochondrial autophagy? This is exactly what we need to study further.
Conclusion
In this study, we assessed the mechanisms through which parkin and BYQZF affect mitochondrial function. Our results showed that the effects of BYQZF on mitochondrial dysfunction were dependent on parkin expression. Thus, we identified parkin as a mitochondrial target for the treatment with BYQZF in PD. Additionally, we found that the mechanism through which BYQZF affected mitochondrial function involved regulation of mitochondrial dynamics via parkin. Taken together, our findings demonstrate that BYQZF ameliorated mitochondrial dysfunction by regulating mitochondrial dynamics via parkin, suggesting that this traditional Chinese medicine can be applied for treating PD.
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 authors.
Author Contributions
H-JM and CG contributed to the interpretation of the results and writing of the manuscript. H-MS contributed to the study design and interpretation of the results. Y-SG and Z-YG instructed experiments. YC, W-DF, C-CC, J-KZ, Y-XZ, L-PY, and Z-YG conducted the experiments. All authors reviewed the manuscript.
Funding
This work is supported by the National Natural Science Foundation of China (Grant Nos. 81774110 and 81573773).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.577017/full#supplementary-material.
References
Bai, H., Ding, Y., Li, X., Kong, D., Xin, C., Yang, X., et al. (2020). Polydatin protects SH-SY5Y in models of Parkinson’s disease by promoting Atg5-mediated but parkin-independent autophagy. Neurochem. Int. 134, 104671. doi:10.1016/j.neuint.2020.104671
Bankapalli, K., Vishwanathan, V., Susarla, G., Sunayana, N., Saladi, S., Peethambaram, D., et al. (2020). Redox-dependent regulation of mitochondrial dynamics by DJ-1 paralogs in Saccharomyces cerevisiae. Redox Biol. 32, 101451. doi:10.1016/j.redox.2020.101451
Bender, A., Krishnan, K. J., Morris, C. M., Taylor, G. A., Reeve, A. K., Perry, R. H., et al. (2006). High levels of mitochondrial DNA deletions in substantia nigra neurons in aging and Parkinson disease. Nat. Genet. 38, 515–517. doi:10.2174/1874609808666150122150850
Camilleri, A., and Vassallo, N. (2014). The centrality of mitochondria in the pathogenesis and treatment of Parkinson’s disease. CNS Neurosci. Ther. 20, 591–602. doi:10.1111/cns.12264
Chan, K., Shaw, D., Simmonds, M. S., Leon, C. J., Xu, Q., Lu, A., et al. (2012). Good practice in reviewing and publishing studies on herbal medicine, with special emphasis on traditional Chinese medicine and Chinese materia medica. J. Ethnopharmacol. 140, 469–475. doi:10.1016/j.jep.2012.01.038
Chanthammachat, P., and Dharmasaroja, P. (2019). Metformin restores the mitochondrial membrane potentials in association with a reduction in TIMM23 and NDUFS3 in MPP+-induced neurotoxicity in SH-SY5Y cells. EXCLI J. 18, 812–823. doi:10.17179/excli2019-1703
Chen, Z., Chen, S., and Liu, J. (2018). The role of T cells in the pathogenesis of Parkinson’s disease. Prog. Neurobiol. 169, 1–23. doi:10.1016/j.pneurobio.2018.08.002
Chua, K. K., Wong, A., Chan, K. W., Lau, Y. K., Bian, Z. X., Lu, J. H., et al. (2017). A randomized controlled trial of chinese medicine on nonmotor symptoms in Parkinson’s disease. Parkinsons Dis. 2017, 1902708. doi:10.1155/2017/1902708
Dodson, M. W., and Guo, M. (2007). Pink1, Parkin, DJ-1, and mitochondrial dysfunction in Parkinson’s disease. Curr. Opin. Neurobiol. 17, 331–337. doi:10.1016/j.conb.2007.04.010
Exner, N., Lutz, A. K., Haass, C., and Winklhofer, K. F. (2012). Mitochondrial dysfunction in Parkinson’s disease: molecular mechanisms and pathophysiological consequences. EMBO J. 31, 3038–3062. doi:10.1038/emboj.2012.170
Feng, S. T., Wang, Z. Z., Yuan, Y. H., Sun, H. M., Chen, N. H., and Zhang, Y. (2019). Mangiferin: a multipotent natural product preventing neurodegeneration in Alzheimer’s and Parkinson’s disease models. Pharmacol. Res. 146, 104336. doi:10.1016/j.phrs.2019.104336
Finkel, T., Menazza, S., Holmström, K. M., Parks, R. J., Liu, J., Sun, J., et al. (2015). The ins and outs of mitochondrial calcium. Circ. Res. 116, 1810–1819. doi:10.1161/CIRCRESAHA.116.305484
Gai, C., Feng, W. D., Qiang, T. Y., Ma, H. J., Chai, Y., Zhang, S. J., et al. (2019). Da-Bu-Yin-Wan and Qian-Zheng-San ameliorate mitochondrial dynamics in the Parkinson’s disease cell model induced by MPP. Front. Pharmacol. 10, 372. doi:10.3389/fphar.2019.00372
Golpich, M., Amini, E., Mohamed, Z., Azman Ali, R., Mohamed Ibrahim, N., and Ahmadiani, A. (2017). Mitochondrial dysfunction and biogenesis in neurodegenerative diseases: pathogenesis and treatment. CNS Neurosci. Ther. 23, 5–22. doi:10.1111/cns.12655
Gong, X. G., Sun, H. M., Zhang, Y., Zhang, S. J., Gao, Y. S., Feng, J., et al. (2014). Da-bu-yin-wan and qian-zheng-san to neuroprotect the mouse model of Parkinson’s disease. Evid. Based Complement. Alternat. Med. 2014, 729195. doi:10.1155/2014/729195
Grünewald, A., Kumar, K. R., and Sue, C. M. (2019). New insights into the complex role of mitochondria in Parkinson’s disease. Prog. Neurobiol. 177, 73–93. doi:10.1016/j.pneurobio.2018.09.003
Ham, S. J., Lee, D., Yoo, H., Jun, K., Shin, H., and Chung, J. (2020). Decision between mitophagy and apoptosis by parkin via VDAC1 ubiquitination. Proc. Natl. Acad. Sci. U.S.A. 117, 4281–4291. doi:10.1073/pnas.1909814117
Ishiyama, M., Miyazono, Y., Sasamoto, K., Ohkura, Y., and Ueno, K. (1997). A highly water-soluble disulfonated tetrazolium salt as a chromogenic indicator for NADH as well as cell viability. Talanta 44, 1299–1305. doi:10.1016/s0039-9140(97)00017-9
Kalia, L. V., and Lang, A. E. (2015). Parkinson’s disease. Lancet 386, 896–912. doi:10.1016/S0140-6736(14)61393-3
Kandimalla, R., and Reddy, P. H. (2016). Multiple faces of dynamin-related protein 1 and its role in Alzheimer’s disease pathogenesis. Biochim. Biophys. Acta 1862, 814–828. doi:10.1016/j.bbadis.2015.12.018
Kavitha, M., Nataraj, J., Essa, M. M., Memon, M. A., and Manivasagam, T. (2013). Mangiferin attenuates MPTP induced dopaminergic neurodegeneration and improves motor impairment, redox balance, and Bcl-2/Bax expression in experimental Parkinson’s disease mice. Chem. Biol. Interact. 206, 239–247. doi:10.1016/j.cbi.2013.09.016
Kitada, T., Asakawa, S., Hattori, N., Matsumine, H., Yamamura, Y., Minoshima, S., et al. (1998). Mutations in the parkin gene cause autosomal recessive juvenile parkinsonism. Nature 392, 605–608. doi:10.1038/33416
Kordower, J. H. (2015). Gene therapy for Parkinson’s disease: still a hot topic?. Neuropsychopharmacology 40, 255–256. doi:10.1038/npp.2014.235
Krishna, A., Biryukov, M., Trefois, C., Antony, P. M., Hussong, R., Lin, J., et al. (2014). Systems genomics evaluation of the SH-SY5Y neuroblastoma cell line as a model for Parkinson’s disease. BMC Genom. 15 (1), 1154. doi:10.1186/1471-2164-15-1154
Li, H., Ruan, Y., Zhang, K., Jian, F., Hu, C., Miao, L., et al. (2016). Mic60/Mitofilin determines MICOS assembly essential for mitochondrial dynamics and mtDNA nucleoid organization. Cell Death Differ. 23, 380–392. doi:10.1038/cdd.2015.102
Lim, K., Ng, X., Grace, L. G., and Yao, T. (2012). Mitochondrial dynamics and Parkinson’s disease: focus on parkin. Antioxid. Redox Signal. 16, 935–949. doi:10.1089/ars.2011.4105
Lisowski, P., Kannan, P., Mlody, B., and Prigione, A. (2018). Mitochondria and the dynamic control of stem cell homeostasis. EMBO Rep. 19. doi:10.15252/embr.201745432
Liu, X., Lu, J., Liu, S., Huang, D., Chen, M., Xiong, G., et al. (2020). Huangqi-Danshen decoction alleviates diabetic nephropathy in db/db mice by inhibiting PINK1/Parkin-mediated mitophagy. Am. J. Transl. Res. 12, 989–998.
Luo, Y., Hoffer, A., Hoffer, B., and Qi, X. (2015). Mitochondria: a therapeutic target for Parkinson’s disease?. Int. J. Mol. Sci. 16, 20704–20730. doi:10.3390/ijms160920704
Meyer, J. N., Leuthner, T. C., and Luz, A. L. (2017). Mitochondrial fusion, fission, and mitochondrial toxicity. Toxicology 391, 42–53. doi:10.1016/j.tox.2017.07.019
Mortiboys, H., Thomas, K. J., Koopman, W. J., Klaffke, S., Abou-Sleiman, P., Olpin, S., et al. (2008). Mitochondrial function and morphology are impaired in parkin-mutant fibroblasts. Ann. Neurol. 64, 555–565. doi:10.1002/ana.21492
Noda, S., Sato, S., Fukuda, T., Tada, N., Uchiyama, Y., Tanaka, K., et al. (2020). Loss of Parkin contributes to mitochondrial turnover and dopaminergic neuronal loss in aged mice. Neurobiol. Dis. 136, 104717. doi:10.1161/CIRCRESAHA.116.305484
Nunnari, J., and Suomalainen, A. (2012). Mitochondria: in sickness and in health. Cell 148, 1145–1159. doi:10.1016/j.cell.2012.02.035
Obeso, J., Stamelou, M., Goetz, C., Poewe, W., Lang, A., Weintraub, D., et al. (2017). Past, present, and future of Parkinson’s disease: a special essay on the 200th anniversary of the shaking palsy. Mov. Disord. 32, 1264–1310. doi:10.1002/mds.27115
Park, J. S., Davis, R. L., and Sue, C. M. (2018). Mitochondrial dysfunction in Parkinson’s disease: new mechanistic insights and therapeutic perspectives. Curr. Neurol. Neurosci. Rep. 18, 21. doi:10.1007/s11910-018-0829-3
Pendergrass, W., Wolf, N., and Poot, M. (2004). Efficacy of mitotracker Green™ and CMXrosamine to measure changes in mitochondrial membrane potentials in living cells and tissues. Cytometry 61A, 162–169. doi:10.1002/cyto.a.20033
Pernas, L., and Scorrano, L. (2016). Mito-morphosis: mitochondrial fusion, fission, and cristae remodeling as key mediators of cellular function. Annu. Rev. Physiol. 78, 505–531. doi:10.1146/annurev-physiol-021115-105011
Pinto, M., and Nissanka, N. (2018). Lack of parkin anticipates the phenotype and affects mitochondrial morphology and mtDNA levels in a mouse model of Parkinson’s disease. J. Neurosci. 38, 1042–1053. doi:10.1523/JNEUROSCI.1384-17.2017
Plotegher, N., and Duchen, M. R. (2017). Mitochondrial dysfunction and neurodegeneration in lysosomal storage disorders. Trends Mol. Med. 23, 116–134. doi:10.1016/j.molmed.2016.12.003
Rao, S. P., Sharma, N., and Kalivendi, S. V. (2020). Embelin averts MPTP-induced dysfunction in mitochondrial bioenergetics and biogenesis via activation of SIRT1. Biochim. Biophys. Acta Bioenerg. 1861, 148157. doi:10.1016/j.bbabio.2020.148157
Schapira, A. H. (2012). Mitochondrial diseases. Lancet 379, 1825–1834. doi:10.1016/S0140-6736(11)61305-6
Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. doi:10.1038/nmeth.2089
Sliter, D. A., Martinez, J., Hao, L., Chen, X., Sun, N., Fischer, T. D., et al. (2018). Parkin and PINK1 mitigate STING-induced inflammation. Nature 561, 258–262. doi:10.1038/s41586-018-0448-9
Song, J. X., Shaw, P. C., Sze, C. W., Tong, Y., Yao, X. S., Ng, T. B., et al. (2010). Chrysotoxine, a novel bibenzyl compound, inhibits 6-hydroxydopamine induced apoptosis in SH-SY5Y cells via mitochondria protection and NF-κB modulation. Neurochem. Int. 57, 676–689. doi:10.1016/j.neuint.2010.08.007
Stevens, D. A., Lee, Y., Kang, H. C., Lee, B. D., Lee, Y. I., Bower, A., et al. (2015). Parkin loss leads to PARIS-dependent declines in mitochondrial mass and respiration. Proc. Natl. Acad. Sci. U.S.A. 112 (37), 11696–11701. doi:10.1073/pnas.1500624112
Tang, J., Hu, Z., Tan, J., Yang, S., and Zeng, L. (2016). Parkin protects against oxygen-glucose deprivation/reperfusion insult by promoting Drp1 degradation. Oxid. Med. Cell. Longev. 2016, 8474303. doi:10.1155/2016/8474303
Wang, L. Y., Yu, X., Li, X. X., Zhao, Y. N., Wang, C. Y., Wang, Z. Y., et al. (2019). Catalpol exerts a neuroprotective effect in MPTP mouse model of Parkinson’s disease. Front. Aging Neurosci. 11, 316. doi:10.3389/fnagi.2019.00316
Xicoy, H., Wieringa, B., and Martens, G. J. (2017). The SH-SY5Y cell line in Parkinson’s disease research: a systematic review. Mol. Neurodegener. 12 (1), 10. doi:10.1186/s13024-017-0149-0
Yang, S. J., Yang, J. W., Na, J. M., Ha, J. S., Choi, S. Y., and Cho, S. W. (2018). 3-(Naphthalen-2-yl(propoxy)methyl) azetidine hydrochloride attenuates MPP+-induced cytotoxicity by regulating oxidative stress and mitochondrial dysfunction in SH-SY5Y cells. BMB Rep. 51 (11), 590–595. doi:10.5483/BMBRep.2018.51.11.123
Zhang, Y., Sun, H. M., He, X., Wang, Y. Y., Gao, Y. S., Wu, H. X., et al. (2013). Da-bu-yin-wan and qian-zheng-san, two traditional chinese herbal formulas, up-regulate the expression of mitochondrial subunit NADH dehydrogenase 1 synergistically in the mice model of Parkinson’s disease. J. Ethnopharmacol. 146, 363–371. doi:10.1016/j.jep.2013.01.005
Zhang, Y., Gong, X. G., Sun, H. M., Guo, Z. Y., Hu, J. H., Wang, Y. Y., et al. (2018). Da-Bu-Yin-Wan improves the ameliorative effect of DJ-1 on mitochondrial dysfunction through augmenting the Akt phosphorylation in a cellular model of Parkinson’s disease. Front. Pharmacol. 9, 1206. doi:10.3389/fphar.2018.01206
Zhou, X., Liu, Y., Chen, D., and Li, C. (2017). Study of Plastrum testudinis extracts in regulating orthodenticle homeobox 2 gene to promote neural stem cells differentiating into dopaminergic neurons. Tradit. Chin. Drug Res. Clin. Pharmacol. 28, 171–176.
Keywords: mitochondria, MPP, parkin, Bu-Yin-Qian-Zheng formula, Parkinson’s disease
Citation: Ma H-J, Gai C, Chai Y, Feng W-D, Cheng C-C, Zhang J-K, Zhang Y-X, Yang L-P, Guo Z-Y and Sun H-M (2020) Bu-Yin-Qian-Zheng Formula Ameliorates MPP+-Induced Mitochondrial Dysfunction in Parkinson’s Disease via Parkin. Front. Pharmacol. 11:577017. doi: 10.3389/fphar.2020.577017
Received: 28 June 2020; Accepted: 16 November 2020;
Published: 18 December 2020.
Edited by:
Yue Liu, Xiyuan Hospital, ChinaReviewed by:
Peijing Rong, China Academy of Chinese Medical Sciences, ChinaPing Li, Beijing Institute of Transfusion Medicine, China
Copyright © 2020 Ma, Gai, Chai, Feng, Cheng, Zhang, Zhang, Yang, Guo and Sun. 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: Hong-Mei Sun, hongmei_sun@yeah.net
†These authors have contributed equally to this work and share first authorship