Selenium Attenuates TBHP-Induced Apoptosis of Nucleus Pulposus Cells by Suppressing Mitochondrial Fission through Activating Nuclear Factor Erythroid 2-Related Factor 2

Intervertebral disc (IVD) degeneration (IDD), the leading cause of low back pain (LBP), remains intractable due to a lack of effective therapeutic strategies. Several lines of studies have documented that nucleus pulposus cell (NPC) death induced by excessive oxidative stress is a crucial contributor to IDD. However, the concrete role and regulation mechanisms have not been fully clarified. Selenium (Se), a vital prosthetic group of antioxidant enzymes, is indispensable for maintaining redox homeostasis and promoting cell survival. However, no light was shed on the role of Se on IDD progression, especially regulation on mitochondrial dynamics and homeostasis. To fill this research gap, the current study focuses on the effects of Se, including sodium selenite (SS) and selenomethionine (Se-Met), on IDD progression and the underlying mechanisms. In vitro, we found that both SS and Se-Met alleviated tert-butyl hydroperoxide- (TBHP-) induced oxidative stress, protected mitochondrial function, and inhibited apoptosis of NPCs. Further experiments indicated that Se suppressed TBHP-induced mitochondrial fission and rescued the imbalance of mitochondrial dynamics. Promoting mitochondrial fission by carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (FCCP) partially counteracted the cytoprotective effects of Se. Moreover, blocking nuclear factor erythroid 2-related factor 2 (Nrf2) with ML385 proved that the effect of Se on regulating mitochondrial dynamics was attributed to the activation of the Nrf2 pathway. In the puncture-induced rat IDD model, a supplement of Se-Met ameliorated degenerative manifestations. Taken together, our results demonstrated that Se suppressed TBHP-induced oxidative stress and mitochondrial fission by activating the Nrf2 pathway, thereby inhibiting the apoptosis of NPCs and ameliorating IDD. Regulation of mitochondrial dynamics by Se may have a potential application value in attenuating the pathological process of IDD.


Introduction
Low back pain (LBP), the leading cause of disability worldwide, imposes an enormous economic burden on both healthcare and social support systems [1,2]. Intervertebral disc (IVD) degeneration (IDD) has been reported to be the primary contributor to LBP [3,4]. Nucleus pulposus (NP), mainly constituted by NP cells (NPCs) and extracellular matrix (ECM), plays a pivotal role in maintaining the integrity and biomechanical equilibrium of IVD [5,6]. Previous studies have demonstrated that the decline of NPC quantity and quality is one of the major contributors to IDD [7][8][9]. Excessive death of NPCs resulted in the reduction of ECM synthesis with diminished water-binding ability, eventually triggering IDD. Nevertheless, the precise underlying mechanisms of IDD have not yet been fully elucidated, and effective interventions are still lacking.
Recently, ample compelling evidence indicated that oxidative stress is closely associated with the occurrence and development of IDD [10][11][12]. When the dynamic balance between reactive oxygen species (ROS) generation and antioxidant defense is disrupted in some pathological conditions, oxidative stress may occur and then lead to IDD [13][14][15]. Excessive oxidative stress triggers mitochondrial dysfunction, mitochondrial membrane potential collapse, and ATP production decline, finally contributing to apoptosis and senescence of NPCs [15][16][17]. Therefore, alleviating oxidative stress and maintaining mitochondrial homeostasis in NPCs are supposed to be promising therapeutic strategies for IDD.
Mitochondria are highly dynamic cellular organelles that undergo continuous fusion and fission, a process termed mitochondrial dynamics [18]. When cells suffer metabolic or environmental stress, mitochondrial fusion is activated to promote the complementation of damaged mitochondria, and fission contributes to quality control by removing the damaged mitochondria. However, inordinate cellular stress triggers excessive mitochondrial fission, which results in cell apoptosis [19]. Suppressing a high level of mitochondrial fission can rescue defective mitochondria and dramatically increase antioxidative defense [19]. Disrupted mitochondrial dynamics were proved to contribute to apoptosis and senescence of NPCs, while remodeling mitochondrial dynamics was beneficial to NPC survival [20]. Therefore, it is a novel and crucial insight to further explore IDD protection strategies based on regulating mitochondrial dynamics.
Selenium (Se) serves as an indispensable component of antioxidative enzymes, such as glutathione peroxidase (GPX) and thioredoxin reductase (TRXR) [29]. Se was reported to confer various protective effects in the body, such as anti-inflammatory, antioxidative, and antiapoptotic effects [30][31][32][33]. Zhang et al. revealed that Se shielded mouse kidneys from zearalenone-induced oxidative stress and apoptosis [34]. Moreover, Se has been proved to protect chicken hepatocytes against cadmium-induced autophagy and apoptosis through Nrf2-mediated protection [35]. Considering that oxidative stress and apoptosis of NPCs are major contributors to IDD, it is hypothesized that Se may serve as a promising therapeutic agent in alleviating IDD. To our best knowledge, studies regarding the association of Se and IDD are scarce.
To fill this research gap, the current study explored the antioxidative and antiapoptotic effects of Se against TBHPinduced damages in NPCs and whether the Nrf2 signaling pathway could be activated by Se to attenuate mitochondrial fission. Furthermore, we also evaluated whether Se-Met could attenuate puncture-induced IDD in vivo. Taken together, our study further discussed the mechanisms of regulation of oxidative stress and mitochondrial dynamics in NPCs and highlighted the therapeutic potential of Se for relieving IDD.

Rat Primary Nucleus Pulposus Cell Isolation and Culture.
Sprague-Dawley rats (weighing 200-250 g) were purchased from the Experimental Animal Center of Tongji Medical College, Huazhong University of Science and Technology (Wuhan, China). Primary NPCs were isolated from the tails of rats after being euthanized by intraperitoneal injection of an overdose of pentobarbital as described previously [36]. Briefly, NP tissues were mechanically cut into pieces and digested with 0.2% (m/v) type II collagenase (Sigma-Aldrich, St. Louis, MO, USA) at 37°C for 15 min. After centrifugation and resuspension, the NPCs were seeded into culture flasks in DMEM/F-12 (HyClone, Logan, UT, USA) supplemented with 10% fetal bovine serum (FBS; Gibco, Grand Island, NY, USA) and 1% penicillin/streptomycin and then cultured at 37°C in an incubator. The culture medium was changed every three days. NPCs were harvested by trypsinization after achieving 80-90% confluence for subculture. The second passage of cells was used for the following experiments.  Assay Kit (Beyotime) were applied to detect apoptosis levels, the production of cellular ROS and mitochondrial ROS, and mitochondrial membrane potential levels in NPCs, respectively. In brief, NPCs were harvested by trypsinization after treatment and then centrifuged and washed three times with PBS. Then, NPCs were labeled with the above-mentioned kits and detected by a flow cytometer, and the results were analyzed by the FlowJo software (Tree Star, Inc. San Carlos, CA, USA).

ATP Production
Assay. An Enhanced ATP Assay Kit (Beyotime) was used to measure the ATP content of NPCs after different treatments. NPCs were lysed with the ATP lysis buffer and then centrifuged at 12,000 g for 5 minutes.
After that, the supernatants were collected and reacted with the ATP detection working dilution. Meanwhile, cellular protein concentration was quantified by an Enhanced BCA Protein Assay Kit (Beyotime) as previously stated. Finally, luminescence spectrometry (PerkinElmer) was used to detect luminescence activity. The ATP level of NPCs was normalized to cellular protein concentration.
2.10. Malondialdehyde (MDA) Detection. The MDA concentration of NPCs was detected by the Lipid Peroxidation MDA Assay Kit (Beyotime). Briefly, after treatment, NPCs were lysed with lysis buffer and then centrifuged at 12,000 g for 10 minutes. Subsequently, the supernatants were collected and reacted with the MDA detection working dilution for 15 minutes at 100°C. After cooling to room temperature, the mixture of supernatants was centrifuged at 1000 g for 10 minutes. Finally, the supernatants were added to the 96-well 3 Oxidative Medicine and Cellular Longevity plate to measure the absorbance at 532 nm using a spectrophotometer (PerkinElmer). The MDA level of NPCs was normalized to cellular protein concentration.
2.11. Puncture-Induced Rat IDD Model. All animal operations were performed in strict compliance with the medical ethics committee of Tongji Medical College, Huazhong University of Science and Technology. 21 male Sprague-Dawley rats (weighing 200-250 g) were obtained from the Experimental Animal Center of Tongji Medical College, Huazhong University of Science and Technology (Wuhan, China). With the help of the digital table method, rats were randomly distributed into three groups: sham group (n = 7), IDD group (n = 7), and Se-Met group (n = 7). A previously reported puncture-induced rat IDD model was employed in the current study [37,38]. Rats in all three groups were anesthetized by intraperitoneal injection of pentobarbital (50 mg/kg), and then coccygeal 6/7 IVD of each rat tail in the IDD group and the Se-Met group were located by trial radiography. Subsequently, 27-gauge needles were used to perforate the whole layer of the AF, in parallel to the endplates. The needle was then rotated 360°and held in the disc for 1 minute. After puncturing, the rats in the Se-Met group were treated with Se-Met (2 mg/kg/day) by intragastric administration, and rats in the sham group and the IDD group were treated with an equal volume of deionized water every day for 4 weeks. After rats were euthanatized with an overdose of phenobarbital, the corresponding discs were harvested and fixed in 4% paraformaldehyde for 48 h, followed by decalcification with 10% ethylenediaminetetraacetic acid solution and embedded in paraffin.

Histological
Analysis. The paraffin-embedded discs were cut into 4 μm-thick sections. For histology, the sections were deparaffinized, rehydrated, and then stained with hematoxylin and eosin (HE) and Safranin O-Fast Green (SO). In addition, the histopathological grading scores were assessed according to the grading scale [39]. For IHC staining, the sections were deparaffinized, rehydrated, and incubated with 3% H 2 O 2 for 20 minutes at room temperature to eliminate endogenous peroxidase activity. Then, the antigen was retrieved by pressure cooking in 10 mmol/L citrate buffer, pH = 6. Subsequently, the sections were blocked with 10% goat serum albumin for 30 min at room temperature. The sections were then incubated with a primary antibody against Nrf2 (1 : 400, ProteinTech), collagen II (1 : 50, Affinity Biosciences), and MMP13 (1 : 300, ProteinTech) at 4°C overnight. Sections incubated with nonspecific IgG were used as negative controls. Next, the sections were incubated with horseradish peroxidase-conjugated goat anti-rabbit secondary antibody for 1 h at 37°C and counterstained with hematoxylin. For the gray-scale statistics of IHC staining, we used ImageJ to calculate the integrated optical density (IOD) of IHC images.

Statistical
Analysis. In the current study, all experiments were repeated at least three times independently. Data were analyzed using the GraphPad Prism V.8.0 software (GraphPad Software, San Diego, CA, USA), and the results were presented as mean ± standard deviation (SD). The student's t-test was used for comparisons between two independent groups, while for comparisons of more than two groups, one-way analysis of variance (ANOVA) was used and followed by Tukey's post hoc test. P values < 0.05 were taken to be statistically significant. The counteracting effects of both SS and Se-Met on TBHP (100 μM) cytotoxicity were evaluated in NPCs. As shown in Figure 1(d) 3 μM and 5 μM SS presented remarkable protective effect, while 75 μM and 100 μM Se-Met obviously increased the cell viability of NPCs ( Figure 1(e)). Therefore, the above-mentioned concentrations of Se were applied in the following experiments. In addition, live/dead cell staining showed that TBHP treatment decreased the ratio of alive cells (green fluorescence) and increased the ratio of dead cells (red fluorescence), whereas Se reversed the changes (Figure 1(f)).

Se Inhibited TBHP-Induced Oxidative Stress and
Restored Mitochondrial Function. Flow cytometry analysis indicated that the production of intracellular ROS in the TBHP group was significantly higher than the control group (P < 0:01), and Se greatly attenuated the elevation induced by TBHP (Figure 1(g) and Figure S1(a)). Consistently, Se markedly suppressed the TBHP-induced increase of MDA, an indicator of lipid peroxidation (Figure 1(h)). In addition, Se enhanced the expression of the antioxidant enzyme SOD2 in NPCs exposed to TBHP (Figures 1(i) and 1(j)).
Excessive oxidative stress has been widely known as the main cause of mitochondrial damage, which contributes to mitochondrial dysfunction. As shown in Figure 1(k) and Figure S1(b), Se reduced the excessive production of mitochondrial ROS caused by TBHP. Furthermore, the level of mitochondrial membrane potential was rescued from loss under Se pretreatment as indicated by the decreased JC-1 monomer ratio (Figure 1(l) and 1(m)). In addition, ATP production assay showed that Se markedly boosted ATP generation in TBHP-challenged NPCs (Figure 1(n)).

Se Ameliorated TBHP-Induced Apoptosis and ECM
Degeneration of NPCs. After treatment with TBHP (100 μM) for 6 h, the apoptosis rate of NPCs rose to over 30% (Figures 2(a) and 2(b)), which was significantly higher than the control group (P < 0:01). However, pretreatment with Se inhibited the increase of apoptosis rate obviously, as indicated by flow cytometry analysis (Figures 2(a) and 2(b)). 4 Oxidative Medicine and Cellular Longevity    Oxidative Medicine and Cellular Longevity fission is crucially involved in cell survival. IF staining of TOM-20 showed that TBHP treatment resulted in an increased level of mitochondrial fission, as indicated by shortened mitochondrial length compared to the control group. Nevertheless, the change of mitochondria was significantly ameliorated in cells pretreated with Se (Figure 3(a)).
In the TBHP group, the expression levels of proteins responsible for mitochondrial fission (DRP1, MFF, and Fis1) were upregulated, while the expression levels of OPA1, Mfn1, and Mfn2 were remarkably reduced compared to the control group. However, Se treatment relieved TBHPinduced alterations of the above proteins (Figures 3(b) and 3  (c)). Collectively, these results indicated that Se rectified the imbalance between mitochondrial fusion and fission caused by TBHP.  (Figure 4(a)) and western blot (Figures 4(b) and 4(c)). Besides, FCCP partly blunted the beneficial effects of Se on alleviating the production of intracellular ROS (Figure 4(d) and Figure S1(c)) and mitochondrial ROS (Figure 4(e) and Figure S1(d)) and improving mitochondrial membrane potential levels (Figures 4(f) and 4(g)) in TBHP-exposed NPCs.
Through blocking the protective effects of Se, FCCP facilitated apoptosis of NPCs regardless of Se intervention, as proved by flow cytometry analysis (Figures 4(h) and 4 (i)) and TUNEL staining ( Figure S2(a)). Moreover, FCCP dramatically upregulated the expression levels of cleavedcaspase 9, cleaved-caspase 3, and Bax and downregulated Bcl-2 compared to the Se group (Figures 4(j) and 4(k)).
Taken together, these results revealed that Se alleviated oxidative stress, maintained the balance of mitochondrial    3.6. Se Activated the Nrf2 Signaling Pathway in TBHP-Exposed NPCs. We further explored the involvement of the Nrf2 antioxidative signaling pathway in the protective effects of Se on NPCs. TBHP downregulated the protein expression of Nrf2 and its downstream target HO-1 compared with the control group, while Se increased the expression of them in NPCs significantly (Figures 5(a) and 5(b)). Subsequently, IF staining of Nrf2 indicated that Se treatment distinctly facilitated the translocation of Nrf2 into the nucleus, as well as the expression of Nrf2 compared with the TBHP group ( Figure 5(c)).
ML385 was employed to further confirm the effect of Se on activating the Nrf2 signaling pathway. ML385 greatly decreased the expression levels of Nrf2 and its downstream targets (HO-1 and SOD2) in comparison with the Se group (Figures 5(d) and 5(e)). Consistently, IF staining of Nrf2 showed that ML385 dramatically reduced the nucleus translocation of Nrf2, as well as the expression level of Nrf2 compared with the Se group ( Figure 5(f)).

Se Suppressed Mitochondrial Fission and Cell
Apoptosis through Activating the Nrf2 Pathway. Under the treatment of ML385, Se failed to rectify the imbalance of mitochondrial dynamics induced by TBHP (Figure 6(a)). Consistently, ML385 promoted the expression levels of fission-related proteins (DRP1, MFF, and Fis1), while repressed OPA1, Mfn1, and Mfn2 compared to the Se group (Figures 6(b) and 6(c)). Flow cytometry analysis further revealed that ML385 partly abolished the effects of Se on reducing the production of intracellular ROS (Figure 6(d) and Figure S1(e)) and mitochondrial ROS (Figure 6(e) and Figure S1(f)) and maintaining mitochondrial membrane potential levels (Figures 6(f) and 6(g)) in TBHP-exposed NPCs.
Through impairing the protective effects of Se, ML385 triggered the apoptosis of NPCs as proved by flow cytometry analysis (Figures 6(h) and 6(i)) and TUNEL staining ( Figure S2(b)). Additionally, ML385 partly abolished the antiapoptosis effect of Se on regulating the expressions of mitochondrial apoptosis-related proteins (cleaved-caspase 9, cleaved-caspase 3, Bax, and Bcl-2) (Figures 6(j) and 6(k)). Taken together, the Nrf2 pathway was involved in the cytoprotective effects of Se via regulating mitochondrial dynamics.

Se Ameliorated Puncture-Induced Rat IDD In Vivo.
To further verify the therapeutic effects of Se on IDD in vivo, we adopted a puncture-induced rat IDD model and administered Se-Met (2 mg/kg/day) intragastrically. As shown in Figure 7(a), HE and SO staining showed that obvious degenerative variations, including reduced NP size, collapsed disc height, significant loss of cells, and dense ECM, were observed in the IDD group compared to the sham group. Via assessing histological grading scores, we further proved that Se ameliorated puncture-induced IDD (Figure 7(b)). IHC staining showed that compared to the sham group, the levels of Nrf2 and collagen II decreased while the expression level of MMP13 increased in the IDD group. However, these changes were all weakened to a certain extent in the Se-Met group (Figures 7(a) and 7(c)).

Discussion
IDD is emerging as one of the most vital causes of LBP, which represents a major global health problem and leads to tremendous economic burdens. IVD is prone to degeneration due to the combination of environmental and genetic risk factors. Since NP serves as an important component of IVD, apoptosis and dysfunction of NPCs play crucial roles in the process of IDD. Despite extensive research, there remains an urgent need for exploring effective and specific therapeutic strategies to prevent or reverse the progression of IDD, owing to the incomplete elucidation of the underlying mechanisms. Hence, in the present study, we illustrated that Se ameliorated TBHP-induced oxidative stress and rectified the imbalance of mitochondrial dynamics, thereby protecting NPCs from mitochondrial apoptosis. Moreover, by the employment of ML385, we then clarified that the cytoprotective effects of Se on NPCs were related to the Se is an indispensable metalloid trace element that is essential for cell survival. Given the important role in maintaining the optimal redox environment in cells and tissues, a growing body of evidence has implicated that Se is involved in various physiological functions, such as antioxidation, immunity, metabolism, and reproduction, exhibiting beneficial antiapoptosis effects [31][32][33]. Sun et al. employed SS to treat hydrogen peroxide-exposed rat cardiomyoblast cells and found that SS supplementation alleviated hydrogen peroxide-induced oxidative stress [40]. In this study, both SS and Se-Met were proved to alleviate TBHP-induced overproduction of intracellular ROS and elevated MDA levels, as well as upregulate the expression level of SOD2. Moreover, existing research revealed that Se could defend LLC-PK (1) cells against cadmium-induced apoptosis [41]. Our results indicated that TBHP-induced excessive apoptosis of NPCs could be reversed by Se treatment. In addition to Se content, the speciation of Se is also important. Recently, in comparison to traditional inorganic Se (SS), organic Se (Se-Met) has attracted increasing attention on account of its superior properties, including higher bioavailability, more rapid and efficient incorporation, and greater antioxidative capacity [42][43][44]. In line with these findings, both SS and Se-Met exerted protective effects on NPCs in our study, but Se-Met showed a wider biosafety range and lower cytotoxicity. Hence, only Se-Met was then further used in vivo to verify its amelioration of disc degeneration. The HE and SO staining showed that Se-Met significantly attenuated punctureinduced disc degeneration.
It is well established that a healthy pool of mitochondria, with a dynamic equilibrium between consistent fission and fusion, is an important guarantee for cellular survival by supplying adequate energy. When encountering adverse stimuli, mitochondrial dynamics was disrupted and resulted in excessive mitochondrial fission [27]. In this study, TBHP treatment resulted in an increased level of mitochondrial fission. Excessive mitochondrial fission is closely related to mitochondrial dysfunction. Li et al. indicated that dynamin-like protein 1-dependent mitochondrial fragmentation contributed to mitochondrial dysfunction and then resulted in neuron apoptosis [45]. After hypoxia, mitochondrial DRP1 was found to contribute to the immoderate opening of the mitochondrial permeability transition pore, ultimately resulting in mitochondrial dysfunction [46]. Here, the addition of FCCP showed that imbalanced mitochondrial dynamics contributed to the overproduction of mitochondrial ROS and decreased mitochondrial membrane potential levels. Kang et al. found that compression resulted in excessive mitochondrial fission, which triggers the overproduction of cellular and mitochondrial ROS, mitochondrial dysfunction, and apoptosis in human NP cells [27]. Our results showed that FCCP induced excessive mitochondrial fission, thereby contributing to apoptosis of NPCs regardless of the cytoprotective effects of Se. Moreover, suppressing mitochondrial fission and promoting mitochondrial fusion were proven to be beneficial to the   survival of IVD cells [27,47]. Hence, agents targeting mitochondrial fission are promising strategies in the treatment of IDD [20,48].
Se was proved to provide antioxidative effects via activating the Nrf2 signaling pathway. A literature review reported that Se conferred protective effects on cardiovascular disease via several signaling pathways, among which the Nrf2 pathway plays a pivotal role in ameliorating oxidative stress [31]. In the in vitro study, Se notably upregulated the expression levels of Nrf2 and its downstream targets compared to the TBHP group. However, the employment of ML385 partially hindered the antioxidative effects of Se. Under stressed conditions, Nrf2 dissociates from the Keap1-Nrf2 complex and then translocates into the nucleus, eventually promoting the expression of cytoprotective genes that facilitate cell sur-vival [24]. Consistently, IF staining of Nrf2 showed that Se facilitated nucleus translocation of Nrf2 compared to the TBHP group, while ML385 inhibited this process. Notably, a large body of studies supported that activation of the Nrf2 pathway by various drugs protected IVD cells from oxidative stress-induced apoptosis in vitro and attenuated the progression of IDD in vivo [49,50]. Here, Se was proved to protect NPCs from TBHP-induced apoptosis. When ML385 was used to block the Nrf2 pathway, Se failed to rescue NPCs from apoptosis, indicating that Se ameliorated TBHP-induced oxidative stress and apoptosis via activating the Nrf2 pathway. Besides, Tang et al. revealed that the expression level of Nrf2 decreased in human degenerative NP tissue samples and that Nrf2 knockout also aggravated IDD [51]. Accordingly, the in vivo study showed that the (j, k) Representative western blot bands and quantitation of the expression levels of cleaved-caspase 9, cleaved-caspase 3, Bax, and Bcl-2 in NPCs treated with or without ML385. The data are represented as the mean ± SD from at least 3 independent experiments. Significant differences between groups are indicated as * * P < 0:01, * P < 0:05.

16
Oxidative Medicine and Cellular Longevity expression level of Nrf2 decreased in the IDD group compared with the sham group, while it significantly increased in the Se-Met group. The Nrf2 pathway is tightly associated with mitochondrial dynamics [52,53]. Several lines of evidence indicated that activation of the Nrf2 pathway suppressed mitochondrial fission and promoted mitochondrial fusion through the degradation of DRP1 [54][55][56]. Moreover, metformin was reported to ameliorate lead-induced mitochondrial fragmentation via activating the AMPK/Nrf2 pathway [57]. Immoderate mitochondrial fission leads to disturbance of mitochondrial homeostasis, thereby triggering cell death. Here, we confirmed that Se maintained the balance of mitochondrial dynamics via activating the Nrf2 pathway. However, blocking the Nrf2 signaling pathway by ML385 markedly enhanced mitochondrial fission and suppressed mitochondrial fusion compared to the Se group.
There are some limitations in the current study. Although NP can exchange substances with the blood circulation through the osmosis of annulus fibrosus and cartilage endplate, the distribution concentration of Se-Met in NP is not clear by intragastric administration. Moreover, X-ray