Manual therapy regulates oxidative stress in aging rat lumbar intervertebral discs through the SIRT1/FOXO1 pathway

Chongjie Yao; Guangxin Guo; Ruixin Huang; Cheng Tang; Qingguang Zhu; Yanbin Cheng; Lingjun Kong; Jun Ren; Min Fang

doi:10.18632/aging.203949

Research Paper Volume 14, Issue 5 pp 2400—2417

Manual therapy regulates oxidative stress in aging rat lumbar intervertebral discs through the SIRT1/FOXO1 pathway

Chongjie Yao^{1,2,
,} , Guangxin Guo^{1,4,
,} , Ruixin Huang^2, , Cheng Tang^2, , Qingguang Zhu^2,3, , Yanbin Cheng^2,3, , Lingjun Kong^2,3, , Jun Ren^2, , Min Fang^1,2,3, ,

¹ School of Acupuncture-Moxibustion and Tuina, Shanghai University of Traditional Chinese Medicine, Shanghai 201203, P.R. China
² Yueyang Hospital of Integrated Traditional Chinese and Western Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai 200437, P.R. China
³ Research Institute of Tuina, Shanghai Academy of Traditional Chinese Medicine, Shanghai 200437, P.R. China
⁴ Shanghai Municipal Hospital of Traditional Chinese Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai 200071, P.R. China

* Equal contribution

Received: October 28, 2021 Accepted: March 1, 2022 Published: March 15, 2022

https://doi.org/10.18632/aging.203949
How to Cite

Copyright: © 2022 Yao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

With the increasing burden of a globally aging population, low back pain has become one of the most common musculoskeletal disorders, caused mainly by intervertebral disc (IVD) degeneration. There are currently several clinical methods to alleviate back pain, but there is scarce attention paid as to whether they can improve age-related IVD degeneration. It is therefore difficult to conduct an in-depth evaluation of these methods. A large number of clinical studies have shown that manual therapy (MT), a widely used comprehensive alternative method, has effects on pain, the mechanisms of which require further study. In this study, MT was performed on aging rats for 6 months, and their behaviors were compared with those of a non-intervention group of aging and young rats. After the intervention, all rats were examined by X-ray to observe lumbar spine degeneration, and the IVD tissues were dissected for detection, including pathological staining, immunofluorescence, Western bolt, etc. This study demonstrated the possibility that MT intervention delay the lumbar IVD degeneration in aging rats, specifically improving the motor function and regulating senescence-associated β-galactosidase, p53, p21, p16, and telomerase activity to retard the senescence of cells in IVDs. Moreover, MT intervention can modify oxidative stress, increase the expression of SIRT1 and FOXO1 in IVDs and decrease ac-FOXO1 expression, suggesting that MT can reduce oxidative stress through the SIRT1/FOXO1 pathway, thereby playing a role in delaying the aging of IVDs. This study shows that drug-free, non-invasive mechanical interventions could be of major significance in improving the physical function of the elderly.

Introduction

Low back pain (LBP) is a common symptom of intervertebral disc (IVD) degeneration, which is one of the most common musculoskeletal disorders and the leading cause of disability [1]. An IVD is composed of fibrocartilaginous tissue and consists of the nucleus pulposus (NP), annulus fibrosus (AF), and cartilage endplate (CEP) [2]. Unlike other tissues in the body, the IVD’s structural integrity begins to show a degenerative state in early stages of human life as the result of physiological aging, a state aggravated with age [3, 4]. Delaying senescence has therefore, been proposed as a feasible therapeutic target for degenerative disc disease [5]. Despite the reports of certain treatments for relieving pain and symptoms [6–8], these methods’ role in delaying IVD aging is rarely explained on the molecular level, which creates challenges to their clinical application.

Given that the accumulation of senescent disc cells is considered an important factor leading to IVD degeneration, the number of studies discussing disc cell senescence has been steadily increasing [9–11]. Various causes for disc cell senescence have been identified, with oxidative stress one of the major contributors [12]. According to the free radical theory of senescence, the decline of tissue and organ function is closely related to the oxidative stress induced by reactive oxide species (ROS), in particular mitochondrial-derived ROS [13]. With the aging process, ROS levels in IVD increase significantly, including those of hydrogen peroxide, hydroxyl radicals, superoxide anions, and nitric oxide, which are byproducts of cellular oxidative metabolism [14]. As an important mediator, ROS can regulate extracellular matrix metabolism, pro-inflammatory factor phenotype, autophagy, and apoptosis, indicating that oxidative stress in the degenerative disc microenvironment can play a crucial role in the pathological development of disc cell senescence [15, 16].

When cells are exposed to oxidative stress over time, they become senescent and lose the ability to proliferate, showing stagnation of the cell cycle in the G0/G1 phase and the activation of senescence-related p53-p21-phosphorylated retinoblastoma (p-Rb) and p16-pRb pathways [17, 18]. Silent mating type information regulation 2 homolog-1 (SIRT1) is a nicotinamide dinucleotide (NAD+)-dependent deacetylase, which protects against cell senescence by regulating p53, p21, and p16, as well as molecules involved in DNA damage and repair [19]. Previous studies have investigated the role of SIRT1 in delaying aging-related diseases and have explored the relationship between SIRT1 and oxidative stress. On one hand, SIRT1 has been shown to regulate cellular oxidative stress through the deacetylase function [20, 21]; on the other hand, oxidative stress can affect SIRT1 activity by regulating gene expression at the transcriptional level and through post-translational modifications [22, 23]. Recent studies have implied that SIRT1 might relieve oxidative stress induced senescence by regulating the SIRT1/forkhead box O (FOXO)1 pathway. Alvarez et al. [24] showed that decreasing FOXO1 expression led to cell senescence and tissue degeneration in IVDs, indicating that FOXO1 is an important regulator for maintaining IVD homeostasis during aging. As the upstream of FOXO1, SIRT1 can reinforce the resistance to oxidative damage of tissue and delay aging by inhibiting FOXO1 acetylation [25]. However, the functions of the SIRT1/FOXO1 pathway in the aging of NP cells (NPC) has not been well studied.

As a widely used complementary and alternative therapy, manual therapy (MT) is performed by the hands on external areas of the body, which can affect the internal physiological state [26]. In China, Tuina is a traditional style of MT that involves pressing, pinching, and kneading the body [27] and has demonstrated positive effects for treating LBP in clinical practice [28–30]. However, there have been few studies on the effects of MT on tissue degeneration beyond pain relief, although an association has been shown between LBP and age-related IVD degeneration [31]. Van et al. [32] reported that MT can synthesize free myofibroblast protein in muscle and can therefore be used as a potential therapeutic option for aged muscle impairment, indicating that MT might alleviate the tissue degeneration caused by aging. In addition, a number of studies [33, 34] have suggested that MT can reduce oxidative stress, which might be one of the mechanisms of the analgesic effect.

However, whether MT can inhibit the onset of age-related IVD degeneration via antioxidation remains unknown. The objective of this study is therefore to explore whether MT can be a useful intervention on the aging of lumbar IVDs by regulating oxidative stress and to determine the role that the SIRT1/FOXO1 pathway plays in the intervention.

Results

Effects of MT on the behavior of aging rats

Due to rats’ inability to accurately convey their feelings, a single behavioral assay might not precisely capture the presence of LBP in these animals. It is therefore, reasonable to use a combination of stimulation and non-stimulation to evaluate pain. The rats’ paw withdrawal threshold (PWT) can be determined through stimulation to assess their mechanical hyperalgesia. As shown in Figure 1A, the PWT of the young rats increased slowly as they aged from 0 to 6 months, an increase that might be related to the enhancement of muscle strength and adaptation to stimulation. This upward trend remained even in the aging rats; the PWT of the 12-month-old rats continued to increase until 15 months of age. From the age of 16 months, the PWT of the rats in the aging group (AG) was significantly lower than that for the previous month (P < 0.01), until it decreased to 65.00 ± 4.78 at the age of 18 months. The PWT of the manual therapy group (MTG) was significantly higher than that of the AG at the same time point after the first month of intervention (P < 0.05, P < 0.01), and the statistically significant difference remained until the age of 18 months, when the PWT of the MTG remained at 76.17 ± 5.07.

Effects of MT on the behaviors of aging rats. (A) The PWT of rats were obtained by stimulating to assess their mechanical hyperalgesia. Compared with the previous time point of this group, #P ##P *P **P B) The open field test was used to measure mobility and voluntary standing to indirectly reflect whether the rats had LBP. (C–F) Further quantitative analysis was used for the distance traveled, mean and max speed, and the standing time/ resting time of the rats. Values are expressed as the mean ±SD (N = 8). Compared with the YCG, ##P *P **P

Figure 1. Effects of MT on the behaviors of aging rats. (A) The PWT of rats were obtained by stimulating to assess their mechanical hyperalgesia. Compared with the previous time point of this group, ^#P < 0.05, ^##P < 0.01; compared with the AG at the same time point, ^*P < 0.05, ^**P < 0.01. Rats in the YCG were 6 months old at BL; rats in the AG and MTG were 12 months old at BL. (B) The open field test was used to measure mobility and voluntary standing to indirectly reflect whether the rats had LBP. (C–F) Further quantitative analysis was used for the distance traveled, mean and max speed, and the standing time/ resting time of the rats. Values are expressed as the mean ±SD (N = 8). Compared with the YCG, ^##P < 0.01; compared with the AG, ^*P < 0.05, ^**P < 0.01. PWT: paw withdrawal threshold, BL: baseline, YCG: young control group, AG: aging group, MTG: manual therapy group.

The open field test was used to measure mobility and voluntary standing, which is a non-stimulating method to indirectly reflect whether rats have LBP. As shown in Figure 1B, the 6-month-old rats preferred to explore the surrounding environment and run continuously around the field; while the 18-month-old rats seemed to prefer moving slowly at the edge of the field, especially in the corners. Further quantitative analysis (Figure 1C–1F) showed that during the 5-minute test, the distance traveled, the mean and max speed, and the standing time/resting time of the AG were significantly lower than those of the young group (P < 0.01). After 6 months of continuous MT intervention, the distance traveled, the mean speed, and the standing time/resting time of the MTG were significantly higher than those of non-intervened rats of the same age (P < 0.05, P < 0.01); however, the MT intervention did not seem to improve the maximum speed of the aging rats (P > 0.05). These results revealed that the aging rats spent less time standing and exploring, which might be partially improved by the MT intervention.

MT improved degenerative changes in the lumbar IVDs with aging

To intuitively evaluate the age-related changes in the lumbar IVD, we took radiographs of the rats’ lumbar spine in the coronal and sagittal planes (Figure 2A). The X-ray images indicated that the IVDs of the AG were characterized by height loss and blurring of the CEP boundary, which can be quantified by calculating the disc height index (DHI) (Figure 2B). Compared with the 6-month-old rat lumbar IVDs, the 18-month-old lumbar IVDs showed a significant reduction in DHI (P < 0.01), indicating that MT can significantly increase the DHI of aging rats (P < 0.01) (Figure 2C).

Figure 2. MT improved degenerative changes in the lumbar IVDs with aging. (A) Radiographs showed the sagittal and coronal view of representative L5/L6 vertebrae. (B) DHI was measured on sagittal radiographs with the equation below. (C) Quantification analysis was used for DHI. (D) Histological scores of disc degeneration were analyzed. (E) The tissue sections were stained with HE and Safranin O/Fast Green (bar = 500 and 100 μm). The white arrow indicates severe degenerative phenotypes, the black arrow indicates inner layers buckling inward, and the green arrow shows partial necrosis in the disc. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ^##P < 0.01; compared with the AG, ^*P < 0.05, ^**P < 0.01. DHI: disc height index, YCG: young control group, AG: aging group, MTG: manual therapy group.

To observe the degenerative disc phenotypes, the tissue sections were stained with hematoxylin and eosin and Safranin O/Fast Green. The results showed that the young control group (YCG) had healthy discs (Figure 2E), characterized by normal cellularity with a mix of large, vacuolated (notochordal) cells in the NP, and clear demarcation between NP and AF tissues with an intact and tightly arranged structure. On the other hand, the discs of the AG revealed severe degenerative phenotypes, including cell arrangement disorders, pronounced cell loss, and the presence of non-vacuolated cells (Figure 2E, white arrow). The clear distinction between the NP and AF tissue compartments was lost, and some of the inner AF layers showed inward buckling (black arrow). Tissue demarcation between the NP and CEP was indistinct and showed partial necrosis (green arrow). Similarly, the junction between the AF and NP was blurred in the discs of the MTG due to aging. Compared with the AG, however, the rats that underwent MT had improved degenerative disc phenotypes, including cell loss, altered cell distribution and morphology, thus maintaining IVD function. The histological scores revealed (Figure 2D) that the IVD scores of the AG were significantly higher than those of the YCG (P < 0.01), while the scores for the MTG were significantly lower than those of the AG (P < 0.05).

Taken together, these data suggest that aging rats develop degenerative IVDs, resulting in structural destruction and cellular disorder, and that MT can alleviate the age-related changes in IVDs.

MT altered the redox balance of lumbar IVDs by promoting the antioxidant capacity

To assess the effects of MT on the oxidative balance of IVDs, we measured ROS levels, superoxide dismutase (SOD) activity, and the malondialdehyde (MDA) and glutathione (GSH) content. The ROS expression levels in the NP tissues of the aging rats were significantly higher than those of the young rats (P < 0.01), suggesting that aging can cause more oxidative stress damage to tissues (Figure 3). As the product of lipid peroxidation, the MDA content of the aging rats was also significantly higher than that of the young rats (P < 0.01). In contrast, SOD and GSH are important antioxidant enzymes for countering the overproduction of active free radicals. However, the SOD activity of the NP in the AG was significantly lower than (almost half) that in the YCG (58.06 ± 17.86 vs. 115.71 ± 18.20; P < 0.01). Although the GSH content did not decrease substantially (11.38 ± 0.86 vs. 13.99 ± 1.07), the difference was statistically significant (P < 0.01). The results indicated that aging decreased IVD antioxidant capacity, disrupting the original redox balance.

MT modified the redox balance of lumbar IVDs by promoting antioxidant capacity. (A–D) The level of ROS, the activity of superoxide dismutase SOD, and the content of MDA and GSH were detected to assess the effects of MT on oxidative balance. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ##P **P

Figure 3. MT modified the redox balance of lumbar IVDs by promoting antioxidant capacity. (A–D) The level of ROS, the activity of superoxide dismutase SOD, and the content of MDA and GSH were detected to assess the effects of MT on oxidative balance. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ^##P < 0.01; compared with the AG, ^**P < 0.01. ROS: reactive oxide species, SOD: superoxide dismutase, MDA: malondialdehyde, GSH: glutathione, YCG: young control group, AG: aging group, MTG: manual therapy group.

The lack of antioxidant capacity and redox imbalance in the lumbar IVDs might be modified by the MT intervention. Compared with the non-intervened aging rats, the rats that underwent MT had enhanced SOD activity and GSH content and decreased ROS and MDA content, differences that were statistically significant (P < 0.01).

MT delayed NPC senescence in lumbar IVDs

To objectively evaluate the effects of MT on aging rat lumbar IVD tissue, we measured the cellular senescence-related indexes. Senescence-associated β-galactosidase (SA-β-Gal) staining showed that the blue positive cells (Figure 4B, black arrow) were densely arranged in the IVD tissues of the aging rats, compared with few positive cells in the young rats. Further statistical analysis (Figure 4B) revealed that the SA-β-Gal-positive cell rate in the IVD tissues of the young rats was significantly lower than that of the aging rats (P < 0.01). Moreover, NPC telomerase activity was indirectly reflected by the relative mRNA expression of telomerase reverse transcriptase detected by polymerase chain reaction (PCR). The results (Figure 4C) showed that the NPC telomerase activity in the AG was significantly lower than that in the YCG (P < 0.01). These two indicators can be used as relatively specific cell sensitivity markers. Using these markers, we found that MT can significantly reduce the SA-β-Gal-positive cell rate and improve telomerase activity (P < 0.01), which suggests that MT can play an important role in delaying cell aging.

MT retarded the senescence of NPCs in lumbar IVDs. (A) SA-β-Gal staining showed the cellular senescence in discs (bar = 100μm). The blue cells indicated by the black arrow are senescence-related positive cells. (B) The percentage of SA-β-Gal positive cells was calculated. (C) Telomerase activity of the NPCs was indirectly reflected by the relative mRNA expression of telomerase reverse transcriptase detected by PCR. (D–G) Western blot analysis showed the relative expressions of senescence-related p53, p21 and p16 proteins in NPCs. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ##P **P

Figure 4. MT retarded the senescence of NPCs in lumbar IVDs. (A) SA-β-Gal staining showed the cellular senescence in discs (bar = 100μm). The blue cells indicated by the black arrow are senescence-related positive cells. (B) The percentage of SA-β-Gal positive cells was calculated. (C) Telomerase activity of the NPCs was indirectly reflected by the relative mRNA expression of telomerase reverse transcriptase detected by PCR. (D–G) Western blot analysis showed the relative expressions of senescence-related p53, p21 and p16 proteins in NPCs. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ^##P < 0.01; compared with the AG, ^**P < 0.01. NPC: nucleus pulposus cell, SA-β-Gal: senescence-associated β-galactosidase, YCG: young control group, AG: aging group, MTG: manual therapy group.

The western blot analysis (Figure 4D–4G) showed that the expression of senescence-related p53, p21, and p16 proteins in the NPCs of the AG were significantly higher than those of the YCG (P < 0.01) and that MT can significantly lower the expression of these proteins when compared with the non-intervened aging rats (P<0.01).

Taken together, these data suggest that aging induces a decrease in NPC telomerase activity in the rat lumbar IVDs, which can activate p53, p21, and p16 proteins regulating cell senescence, thereby accelerating the aging process. MT intervention can inhibit the reduction of telomerase activity, reducing the expression of senescence-related proteins, thereby delaying NPC senescence in lumbar IVDs.

MT exerted its effects by regulating the SIRT1/FOXO1 pathway

To investigate the possible mechanism by which MT reduces lumbar IVD degeneration, we analyzed the changes in the SIRT1/FOXO1 pathway in the NPCs. Immunofluorescence staining showed that SIRT1 (Figure 5A, red fluorescent) and ac-FOXO1 (Figure 5A, green fluorescent) were co-expressed in the NPCs of the MTG (Figure 5A, white arrow). However, ac-FOXO1 was difficult to observe in the NPCs of the YCG, and there was a large amount of SIRT1. In contrast, there was only a small amount of SIRT1 in the AG, and ac-FOXO1 was mainly expressed in the NPCs. These results suggest that SIRT1, which is abundantly expressed in the NPCs of the young rats, can inhibit FOXO1 acetylation to maintain NPC homeostasis; the lack of SIRT1 in the aging rats can increase the degree of FOXO1 acetylation and disrupt the original homeostasis. The MT intervention can significantly increase SIRT1 expression in the NPCs and reduce ac-FOXO1 expression (P < 0.01), thereby reducing the effects of NPC aging (Figure 5B, 5C).

MT exerted its effects through regulating the SIRT1/FOXO1 pathway. (A) Immunofluorescence staining showed the expressions of SIRT1 (red fluorescent) and ac-FOXO1 (green fluorescent) (bar = 50 μm). The white arrow indicates that the two proteins were co-expressed in the NPCs. (B,C) The fluorescence intensity of SIRT1 and ac-FOXO1 was analyzed. (D–G) A Western blot analysis showed the relative expressions of SIRT1, FOXO1 and ac-FOXO1 proteins in NPCs. (H,I) PCR results showed the relative mRNA expressions of SIRT1 and ac-FOXO1. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ##P **P

Figure 5. MT exerted its effects through regulating the SIRT1/FOXO1 pathway. (A) Immunofluorescence staining showed the expressions of SIRT1 (red fluorescent) and ac-FOXO1 (green fluorescent) (bar = 50 μm). The white arrow indicates that the two proteins were co-expressed in the NPCs. (B,C) The fluorescence intensity of SIRT1 and ac-FOXO1 was analyzed. (D–G) A Western blot analysis showed the relative expressions of SIRT1, FOXO1 and ac-FOXO1 proteins in NPCs. (H,I) PCR results showed the relative mRNA expressions of SIRT1 and ac-FOXO1. Values are expressed as the mean ± SD (N = 8). Compared with the YCG, ^##P < 0.01; compared with the AG, ^**P < 0.01. SIRT1: silent mating type information regulation 2 homolog-1, FOXO1: forkhead box O1, YCG: young control group, AG: aging group, MTG: manual therapy group.

The western blot analysis (Figure 5E) also showed changes in the SIRT1, FOXO, and ac-FOXO1 proteins in the pathway. The results for SIRT1 and ac-FOXO1 were consistent with those of the immunofluorescence (Figure 5D, 5G). Compared with the YCG, the AG had significantly decreased FOXO1 expression in the NPCs (P < 0.01), which might have been caused by acetylation of a large amount of FOXO1, and MT might indirectly increase FOXO1 expression by inhibiting acetylation (P < 0.01). In addition, the PCR results for SIRT1 and FOXO1 further confirmed the above results (Figure 5H, 5I). Thus, MT might exert its effects by regulating the SIRT1/FOXO1 pathway.

Discussion

IVD degeneration is a common risk factor for LBP in the elderly, which is a characteristic of aging [31]. Human lumbar IVDs undergo age-related pathological changes, including collapse, herniation, and tears, which can be significant causes of related symptoms [35]. In particular, patients with LBP have shorter lumbar IVDs than those without LBP, which is particularly common in the lumbosacral level L5/S1, followed by L4/L5, which suggests that the vulnerability of these levels in humans might be partly due to the vertical pressure exerted by their upright posture [36]. From this point of view, finding a suitable animal model appears difficult. Most studies [37–39] on disc degeneration have therefore used surgically created animal models to cause similar changes. However, IVD degeneration also reveals that age-related changes in IVD structure might be more directly reflected in quadrupeds. We therefore, chose aging rats as the research subjects in this experiment. Despite the rats’ quadrupedal stance, they often need to stand upright to eat and explore their external environment, which shows that rat lumbar IVDs might have similar mechanical properties to those of humans.

With the growing problem of the aging global population, MT has become an effective means to improve the health of the elderly and even act as a developmental strategic plan in a number of countries [40]. However, there are few studies on the effects of MT on age-related diseases. Previous studies have suggested that MT can improve the age-related symptoms of postmenopause [41], Parkinson’s disease [42], and dementia [43], and the effects are more obvious with long-term interventions. These effects by MT are thought to be achieved through the stimulation of baroreceptors on the body surface, transmitting signals to the vagus nerve, so as to inhibit the function of the hypothalamic-pituitary-adrenocortical axis [44], which is also a vital method for delaying aging [45]. However, we still have no direct evidence of the effectiveness of MT, which makes its promotion difficult due to the lack of evidence.

In this study, we observed the changes in the mechanical stimulation threshold and spontaneous behavior of rats of different ages. The rats exhibited behavioral changes related to IVD degeneration with increasing age, results that were partially similar to those of a previous study [31]. The rats’ sensitivity to mechanical stimulation was evaluated by measuring their PWT on a monthly basis. The PWT of the non-intervened rats showed an upward trend from 0 to 6 months and from 12 to 15 months, and decreased significantly compared with the previous month from the age of 15 months, indicating that their PWT began to decrease. We also compared the changes in spontaneous behavior between the 18-month-old and 6-month-old rats and found that the young rats had good motor ability and often stood on their hind limbs to free their front paws for exploring. The aged rats kept moving and exploring novel environments but opted to remain on all four limbs on the ground. According to a previous report [31], the aged rats’ reluctance to stand is unlikely to be due to muscle atrophy or loss of strength because muscle mass increases with age and remains proportional to body weight. Instead, the reluctance might be related to discomfort due to LBP or other musculoskeletal problems when standing. Taking into consideration the exact efficacy of MT in alleviating LBP [7, 46], we observed the effects of MT on the aging of rat IVDs by pressing and kneading the Shenshu (BL23) points. We selected BL23 points as the MT site because they are located on both sides of the spine, so mechanical loading can directly relax the lumbar muscles, and they have been reported to delay the effects of aging [47–49]. Through behavioral testing, we found that the pain tolerance and behavior style of the aged rats undergoing the MT intervention were significantly improved. However, we were unable to determine whether MT alone alleviated LPB or simultaneously relieved lumbar IVD degeneration; we therefore conducted the following study.

A longitudinal study of the elderly revealed the association between LBP and loss of IVD height [50]. In our study, radiographs showed the sagittal and coronal views of representative L5/L6 vertebrae, indicating loss of IVD height, increased IVD width, and reduced DHI, which emerge during old age in rats, in a similar manner as in humans [51]. We also observed histological changes at 18 months in the rats’ lumbar IVDs. As the IVD ages, severe degenerative structural defects such as cell arrangement disorders, pronounced cell loss, and the presence of non-vacuolated cells become more common, indicating that progressive degenerative changes occur spontaneously with advanced age in rats despite the absence of injury and lifestyle-related risk factors that cause IVD degeneration in humans. Our results support the effects of MT on the structural improvement of rat lumbar IVDs, indicating that MT can not only alleviate LBP but also delay IVD degeneration.

Senescence in IVD cells, especially NPCs, has recently emerged as a major driver of disc aging, leading to changes in IVD structure and function [52]. SA-β-gal is a relatively specific cell senescence marker that reflects increased expression of lysosomal β-galactosidase protein in senescent cells [53]. Compared with undegenerated IVDs, degenerated discs show an increase in SA-β-gal-positive cells [54]. Moreover, telomerase activity is another marker for evaluating cellular senescence [55], because decreased telomerase activity can activate the expression of tumor suppressor gene p53 and p16 by mimicking DNA damage signals, which in turn reduces the degree of pRB protein through the p53-p21-pRB and p16-pRB pathways, eventually inducing cellular senescence. In our study, the increase in SA-β-gal-positive cells and senescence-related p53, p21, and p16 protein expression and the decrease in telomerase activity in the NPs of the aged rats suggest cellular senescence in degenerative IVDs. We found that the degree of IVD aging in the aged rats that underwent 6 months of MT was lower than that of the rats that did not undergo the intervention, and the expression of the senescence-related proteins was also significantly reduced, suggesting that MT delayed NPC senescence in the lumbar IVDs. Oxidative stress induced by various risk factors has been considered closely associated with NPC senescence, and excessive oxidative stress results in accelerated aging and IVD degeneration [56]. Our study verified the results of previous studies [2, 57] and showed the antioxidant capacity and redox imbalance in aged discs, which was specifically manifested in the large increase in ROS and decrease in antioxidant enzymes, resulting in the accumulation of peroxide products. Our study confirmed that MT can modify the redox balance of lumbar IVDs by promoting the antioxidant capacity, which might be one of the important mechanisms of MT in delaying IVD degeneration.

A previous study [24] observed a reduction in FOXO1 expression in mouse lumbar IVD in early aging, which compromised the ability of IVD cells to neutralize ROS, suggesting that FOXO1 expression and activity might have an impact on disc cell homeostasis by decreasing oxidative damage. However, the mechanism by which FOXO1 expression is reduced in aging IVD has not been studied. Taking into consideration the SIRT1 function in mitigating senescence-related dysfunction in NPCs [58], we assumed that FOXO1 might be regulated through the SIRT1/FOXO1 pathway. To confirm our hypothesis, we studied the changes in SIRT1, FOXO1, and ac-FOXO1 in the pathway. Immunofluorescence staining showed that SIRT1 was abundantly expressed in young NPCs, while ac-FOXO1 was difficult to find, which was in contrast to the expression trend in aged NPCs. The results of the western blot and PCR were consistent with the above trends. In the NPCs of the aged rats with low SIRT1 levels, FOXO1 expression was also decreased, while ac-FOXO1 expression was increased. Through the MT intervention, SIRT1 levels in the NPCs of the aged rats increased significantly, suggesting that the decrease in ac-FOXO1 might be due to SIRT1 deacetylation, which indirectly increased FOXO1 expression by promoting the nuclear translocation and activation of FOXO1.

There are a number of limitations that should be noted in this study. First, due to the short project time, the overall changes in the rats from 0 to 18 months were not be observed; instead, we compared the aged and young rats, which might have a certain effect on the results (especially the behavior tests). Second, the quantitative standard of the MT intervention referred to our previous study on LBP, and further study is needed to determine whether there was a more suitable standard for this research. While the present findings demonstrated that MT might alleviate oxidative stress through the SIRT1/FOXO1 pathway and delay IVD degeneration caused by aging, a causal relationship needed to be established. Thus, we need to improve the experimental design in future studies.

Nevertheless, our study confirmed the view that MT might be beneficial for patients with age-associated diseases, for which immunosuppressants are often ineffective; and physical therapy should be recommended as an alternative [59]. With the problem of the global aging population becoming increasingly serious, it is of major importance to develop external mechanical interventions to improve the physiological function of the elderly in a drug-free, non-invasive manner.

In conclusion, this study demonstrated the possibility that an MT intervention can delay lumbar IVD degeneration in aged rats, specifically improving the motor function and delaying NPC senescence. Moreover, the MT intervention can regulate oxidative stress, increase SIRT1 and FOXO1 expression in NPCs and decrease ac-FOXO1 expression, which suggests that MT can reduce oxidative stress through the SIRT1/FOXO1 pathway, thereby playing a role in delaying IVD aging.

Materials and Methods

Animals

Eight young male Sprague–Dawley rats (2 weeks) and 16 adult males (12 months, weight 700 ± 50 g) were obtained from Shanghai Jihui Laboratory Animal Co., Ltd. [certificate number: SCXK (Shanghai)2017-0012] and housed in the Laboratory Animal Center of Yueyang Hospital of Integrated Traditional Chinese and Western Medicine. All the rats were fed in standard laboratory conditions (22~27° C, 50%–70% indoor humidity) under a 12-h light-dark cycle with rat chow and water ad libitum for one week before the experiment. To alleviate pain and avoid injury, the animal feeding, caring, and intervention were conducted in accordance with procedures approved by the Experimental Animal Ethics Committee of Yueyang Hospital of Integrated Traditional Chinese and Western Medicine (ethics number: YYLAC-2020-095).

Grouping and intervention

The 16 adult male Sprague-Dawley rats were randomly divided into the aged group (AG) and the manual therapy group (MTG) (N = 8). The 8 immature male rats comprised the young control group (YCG). The MTG underwent manual therapy on the acupoint beside the spine from 12 months of age. The specific operation was as follows: 1) before the intervention, the rats were placed in the experimenters’ hands 30 min prior to the operation to adapt; (2) the right thumb was used to press and knead the Shenshu (BL23) points on the back of the rat, which are located on both sides of the second lumbar vertebrae and on the erector spinalis muscle (Figure 6); (3) to maintain the amount of stimulation during the manipulation, the right thumb was operated with a wireless tactile measurement finger sleeve (FingerTPS, Pressure Profile System, USA), with a pressure of 5 N and a frequency of 2 Hz [27]; 4) the rats underwent the intervention alternately on the left or right side of the spinal column, one acupoint each time for 10 min, once a day, 6 days a week for 6 months; and 5) after the MT intervention, the rats were pacified before being put into the cage. The AG and YCG underwent no intervention but were grabbed and fixed in the same way as the MTG.

Figure 6. Location of Shenshu (BL23) point and IVD levels used for each outcome. IVD: intervertebral disc.

Behavioral testing

Behavioral tests were conducted for 6 months during the light cycle between 9 am and 5 pm. The rats were habituated to the experimental environment for 2 h daily 3 days before testing. The mechanical hyperalgesia test was performed on the rats every month, once a day for 3 consecutive days, and the mean results of the 3 tests were taken as the final behavioral test results of the month. The open field test was performed after the 6-month intervention, repeated tests and data processing method was the same as above. The tests were implemented, recorded and supervised by 3 experimenters, respectively, and rotated after each test to reduce the experimental error as much as possible.

An electric von Frey aesthesiometer (IITC/Life Science, USA) was applied to test the mechanical hyperalgesia of the rats, as previously reported [60]. The tip head was aimed at the central part of the left hind foot of the rat and was raised at a constant speed for stimulation. When the aesthesiometer recorded the maximum amount of stimulation making the rat withdraw, and the value of this force was assigned to the paw withdrawal threshold (PWT). To prevent sensory sensitization caused by frequent stimulation, the right foot should be stimulated at least 20 s after the measurement of the left foot. The left and right feet of each rat were repeatedly measured 3 times. The maximum and minimum of the 6 values were deleted, and the mean of the remaining 4 values was deemed as the rat’s PWT in one test [61].

The open field test consisted of a black square field (100×100 cm) surrounded by 30-cm high walls made of translucent plastic, which was placed in a quiet and well-lit room. The test was initiated by placing a single rat in the center of the field and allowing it move freely for 5 min while being photographed by a camera. The field was cleaned with low alcohol after every test, and the next test was carried out after the alcohol evaporated completely. The behaviors of the rats were measured by the experimenter and by analysis software (SMART v3.0); the main parameters were as follows: (1) distance traveled: the total travel distance of the rats in the field; (2) mean speed and max speed: the average speed of the rats in the field and the maximum speed at a certain time point; (3) resting time: the time during which the rats did not move horizontally (including standing time); (4) standing time: the time of rats standing on their hind legs (including wall climbing).

Radiography and relative height of disc measurement

After the final behavioral testing, images of the lumbar spine were captured with X-ray imaging and irradiation systems (Samsung XGEO, South Korea). Digital images were obtained at 55 kVp under identical imaging conditions using identical acquisition parameters. The rats were anesthetized with 2% pentobarbital sodium and placed in the X-ray chamber on their ventral side with all limbs extended away from the midline in a natural position to obtain coronal images. The rats were then placed on their side in a natural position with the tail extended to obtain sagittal images. Disc height indices (DHI) were measured using ImagePro plus 6.0, by averaging measurements from the anterior, middle, and posterior portions of the IVD (L5/L6) from a sagittal plane and dividing it by the average measurements of the anterior, middle, and posterior portions of the adjacent vertebral bones [62].

Histological staining

After 6 months of intervention, all rats were euthanized. Intact lumbar spinal segments and IVD tissues from L5/L6 were collected and fixed in 4% formaldehyde for 48h, decalcified in 10% ethylenediaminetetraacetic acid (EDTA) for 1 month. Tissues were then dehydrated in gradient alcohol and embedded in paraffin. 4 μm sections were cut and stained with either hematoxylin and eosin (HE) or Safranin O/Fast Green according to standard procedures. Images were acquired using a microscope (Leica, Germany) and were scored by 2 independent observers using the modified Thompson grading scale [63] to evaluate the degeneration of IVDs.

Detection of redox balance in IVDs

Fresh NP tissues dissected from L2-L5 were ground and homogenized then lysed thoroughly in RIPA Lysis Buffer (Beyotime, China) at a concentration of 10ul/mg. The cell lysate was centrifuged at 12,000×g for 5 min, and the concentration of the supernatant was calculated by using a BCA protein concentration assay kit (Beyotime, China). Based on the concentration, the tissue homogenate was diluted to 1 mg/ml with PBS.

Test kits (Nanjing Jiancheng, China) were used to detect the level of ROS, the activity of SOD, and the content of MDA and GSH according to the instructions. Additionally, 2, 7- dichlorofluorescein diacetate (DCFH-DA) was used as a probe for ROS, and the fluorescence intensity of ROS was measured at 488nm/525nm with a microplate reader (Thermo Scientific, USA) to determine the level of ROS.

Senescence-associated β-galactosidase activity assay

Fresh IVDs were dissected from L6/S1, directly embedded in optimal cutting temperature compound, and stored at −80° C to detect the activity of SA-β-Gal in IVDs. Next, 5 μm sections were cut and stained with SA-β-Gal kit (Beyotime, China) according to the manufacturer’s protocol. Before mounting them with resin, sections were stained again with eosin for 10 sec. The number of cells positive (blue) for localization of SA-β-Gal and total cells (blue and red) in the field were counted using the ImagePro plus 6.0 under 400× magnification. The ratio was recorded as the percentage of positive cells. Three fields were randomly selected from each section to calculate the average percentage of SA-β-Gal positive senescent cells.

Immunofluorescence

After deparaffinization and rehydration, the prepared paraffin sections were treated with citric acid buffer (0.01 mol/L, pH6.0) at 98° C for 15 min for antigen retrieval and were cooled to room temperature for 1 h, followed by incubation with 1% hydrogen peroxide solution for 10 min to block endogenous peroxidase activity. The sections were then incubated in 5% bovine serum albumin (BSA) (Beyotime, China) at room temperature for 30 min, followed by primary antibody incubation (SIRT1, 1:50, CST, 8469s; ac-FOXO1, 1:50, Affinity Biosciences, AF2305) overnight at 4° C. After 12 h, the sections were rewarmed in an oven at 37° C for 45 min and incubated with fluorescein-conjugated secondary antibodies (1:100, Sanggon Biotech, China) at room temperature for 30 min in the dark. Nuclei were stained with 4,6-diamidino-2-phenylindole (DAPI) (Beyotime, China) for 2 min. Three fields were randomly selected from each section to calculate the average fluorescence intensity of the target protein under 400× magnification.

Western blot

The extraction protocol for total protein quantification of fresh NP tissues was described above. After adjusting the protein concentration, 20 μg tissue lysates were separated by 12% SDS-PAGE gel using electrophoresis and transferred to PVDF membranes. The membranes were incubated with primary antibodies (p53, 1:1000, Abcam, ab26; p21, 1:1000, Abcam, ab109520; p16, 1:1000, Abcam, ab51243; p53, 1:1000, Abcam, ab26; SIRT1, 1:1000, CST, 8469s; FOXO1, 1:1000, CST, 2880s; ac-FOXO1, 1:1000, Affinity Biosciences, AF2305; GAPDH, 1:5000, CST, 5174s) overnight at 4 after blocked by 5% BSA at room temperature for 2 h. The membranes were rinsed 3 times with TBST and incubated with horseradish peroxidase (HRP)-labeled secondary antibody (1:5000, Sanggon Biotech, China) at room temperature for 2 h. The final results were acquired using the ChemiDoc MP system (Bio-Rad, USA) with enhanced chemiluminescence (ECL) kit (Millipore, USA) and were quantified using the ImagePro plus 6.0.

Quantitative real-time PCR

Reverse transcription was performed with a PrimeScript RT reagent Kit with gDNA Eraser (TaKaRa Bio, Japan). Real-time polymerase chain reaction (RT-PCR) was performed on a PCR System (Bio-Rad, USA) with SYBR Premix Ex Taq (TaKaRa, Japan). The reaction conditions were a single cycle of 95° C for 5 min, and 40 cycles of 95° C for 15 sec, 60° C for 30 sec. GAPDH was used as a reference, each sample was tested in triplicate. The relative expression of the target gene was calculated using the 2^-ΔΔCt method. The primers used in the study are shown in Table 1.

Table 1. Primers used for RT-PCR analysis.

Gene	Forward (5′-3′)	Reverse (5′-3′)
GAPDH	AGGTGACCGCATCTTCTTGT	TACGGCCAAATCCGTTCACA
FOXO1	CCCAATCTCGGAGCGACAC	GCAGGCTCAGGTTGCTCATA
SIRT1	ATCTCCCAGATCCTCAAGCCA	CTTCCACTGCACAGGCACAT
TERT	TTCCTTCCACCAGGTGTCATC	AGCCAGCACATTCCTCTCAC

Statistical analysis

SPSS 21.0 software was used for the statistical analysis. Normality and homogeneity of variance tests were performed on the data. The data are presented as mean ± standard deviation (SD) if they were normally distributed; data that did not conform to the normal distribution are presented as the median (P25 and P75). If the data conformed to the normal distribution, the comparison of differences among groups was performed using one-way ANOVA. If the variances were homogeneous, a pairwise comparison was performed using the least significant difference test; if the variances were not homogenous, Dunnett’s T3 test was used for pairwise comparison. If the data did not conform to the normal distribution, a nonparametric Kruskal-Wallis H test were performed and using the Nemenyi method for pairwise comparison. Repeated measurement analysis of variance was used to analyze the data on the same observation index at various time points between groups. If the results were statistically different, a multivariate analysis of variance was used for pairwise comparison. The significance level of the statistical examination was α = 0.05. Therefore, P <0.05 indicated a significant difference.

Author Contributions

In this study, CJY and GXG conceived of the project and designed the experiments. CJY, RXH, CT, and JR performed the experiments and acquired the data. GXG, QGZ, YBC, LJK, and MF supervised the project and analyzed the data. CJY, GXG, and MF wrote the manuscript and revised the manuscript in response to reviewers' comments. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This work was supported by the National Natural Science Foundation of China (no. 82030121, 82004493), Innovation Team and Talents Cultivation Program of National Administration of Traditional Chinese Medicine (no. ZYYCXTD-C-202008), Qihuang Scholar of the Talent Project of Traditional Chinese Medicine (2018), Shanghai Post-doctoral Excellence Program (no. 2020371), and the Graduate Student Innovation Ability Project of Shanghai University of Traditional Chinese Medicine (no. Y2020017).

References

1. Zheng L, Cao Y, Ni S, Qi H, Ling Z, Xu X, Zou X, Wu T, Deng R, Hu B, Gao B, Chen H, Li Y, et al. Ciliary parathyroid hormone signaling activates transforming growth factor-β to maintain intervertebral disc homeostasis during aging. Bone Res. 2018; 6:21. https://doi.org/10.1038/s41413-018-0022-y [PubMed]
2. Jin LY, Lv ZD, Wang K, Qian L, Song XX, Li XF, Shen HX. Estradiol Alleviates Intervertebral Disc Degeneration through Modulating the Antioxidant Enzymes and Inhibiting Autophagy in the Model of Menopause Rats. Oxid Med Cell Longev. 2018; 2018:7890291. https://doi.org/10.1155/2018/7890291 [PubMed]
3. Oh CH, Yoon SH. Whole Spine Disc Degeneration Survey according to the Ages and Sex Using Pfirrmann Disc Degeneration Grades. Korean J Spine. 2017; 14:148–54. https://doi.org/10.14245/kjs.2017.14.4.148 [PubMed]
4. Rider SM, Mizuno S, Kang JD. Molecular Mechanisms of Intervertebral Disc Degeneration. Spine Surg Relat Res. 2018; 3:1–11. https://doi.org/10.22603/ssrr.2017-0095 [PubMed]
5. Feng C, Liu H, Yang M, Zhang Y, Huang B, Zhou Y. Disc cell senescence in intervertebral disc degeneration: Causes and molecular pathways. Cell Cycle. 2016; 15:1674–84. https://doi.org/10.1080/15384101.2016.1152433 [PubMed]
6. Froud R, Patterson S, Eldridge S, Seale C, Pincus T, Rajendran D, Fossum C, Underwood M. A systematic review and meta-synthesis of the impact of low back pain on people’s lives. BMC Musculoskelet Disord. 2014; 15:50. https://doi.org/10.1186/1471-2474-15-50 [PubMed]
7. Furlan AD, Giraldo M, Baskwill A, Irvin E, Imamura M. Massage for low-back pain. Cochrane Database Syst Rev. 2015; 2015:CD001929. https://doi.org/10.1002/14651858.CD001929.pub3 [PubMed]
8. Chou R, Deyo R, Friedly J, Skelly A, Hashimoto R, Weimer M, Fu R, Dana T, Kraegel P, Griffin J, Grusing S, Brodt ED. Nonpharmacologic Therapies for Low Back Pain: A Systematic Review for an American College of Physicians Clinical Practice Guideline. Ann Intern Med. 2017; 166:493–505. https://doi.org/10.7326/M16-2459 [PubMed]
9. Jiao S, Li J, Liu B, Yang M, Xiu J, Qu D. Nucleus pulposus cell apoptosis is attenuated by CDMP-2 through regulating oxidative damage under the hyperosmotic environment. Biosci Rep. 2018; 38:BSR20181176. https://doi.org/10.1042/BSR20181176 [PubMed]
10. Xu J, Li H, Yang K, Guo S, Wang J, Feng C, Chen H. Hyper-osmolarity environment-induced oxidative stress injury promotes nucleus pulposus cell senescence in vitro. Biosci Rep. 2019; 39:BSR20191711. https://doi.org/10.1042/BSR20191711 [PubMed]
11. Li X, Lin F, Wu Y, Liu N, Wang J, Chen R, Lu Z. Resveratrol attenuates inflammation environment-induced nucleus pulposus cell senescence in vitro. Biosci Rep. 2019; 39:BSR20190126. https://doi.org/10.1042/BSR20190126 [PubMed]
12. Qin R, Dai S, Zhang X, Liu H, Zhou B, Zhou P, Hu C. Danshen Attenuates Intervertebral Disc Degeneration via Antioxidation in SD Rats. Oxid Med Cell Longev. 2020; 2020:6660429. https://doi.org/10.1155/2020/6660429 [PubMed]
13. Kauppila TE, Kauppila JH, Larsson NG. Mammalian Mitochondria and Aging: An Update. Cell Metab. 2017; 25:57–71. https://doi.org/10.1016/j.cmet.2016.09.017 [PubMed]
14. Nasto LA, Robinson AR, Ngo K, Clauson CL, Dong Q, St Croix C, Sowa G, Pola E, Robbins PD, Kang J, Niedernhofer LJ, Wipf P, Vo NV. Mitochondrial-derived reactive oxygen species (ROS) play a causal role in aging-related intervertebral disc degeneration. J Orthop Res. 2013; 31:1150–7. https://doi.org/10.1002/jor.22320 [PubMed]
15. Schröder K. NADPH oxidases in bone homeostasis and osteoporosis. Free Radic Biol Med. 2019; 132:67–72. https://doi.org/10.1016/j.freeradbiomed.2018.08.036 [PubMed]
16. Yang RZ, Xu WN, Zheng HL, Zheng XF, Li B, Jiang LS, Jiang SD. Involvement of oxidative stress-induced annulus fibrosus cell and nucleus pulposus cell ferroptosis in intervertebral disc degeneration pathogenesis. J Cell Physiol. 2021; 236:2725–39. https://doi.org/10.1002/jcp.30039 [PubMed]
17. Dimozi A, Mavrogonatou E, Sklirou A, Kletsas D. Oxidative stress inhibits the proliferation, induces premature senescence and promotes a catabolic phenotype in human nucleus pulposus intervertebral disc cells. Eur Cell Mater. 2015; 30:89–102. https://doi.org/10.22203/ecm.v030a07 [PubMed]
18. Che H, Li J, Li Y, Ma C, Liu H, Qin J, Dong J, Zhang Z, Xian CJ, Miao D, Wang L, Ren Y. p16 deficiency attenuates intervertebral disc degeneration by adjusting oxidative stress and nucleus pulposus cell cycle. Elife. 2020; 9:e52570. https://doi.org/10.7554/eLife.52570 [PubMed]
19. Lee SH, Lee JH, Lee HY, Min KJ. Sirtuin signaling in cellular senescence and aging. BMB Rep. 2019; 52:24–34. https://doi.org/10.5483/BMBRep.2019.52.1.290 [PubMed]
20. Li Y, Xu W, McBurney MW, Longo VD. SirT1 inhibition reduces IGF-I/IRS-2/Ras/ERK1/2 signaling and protects neurons. Cell Metab. 2008; 8:38–48. https://doi.org/10.1016/j.cmet.2008.05.004 [PubMed]
21. Khan M, Ullah R, Rehman SU, Shah SA, Saeed K, Muhammad T, Park HY, Jo MH, Choe K, Rutten BP, Kim MO. 17β-Estradiol Modulates SIRT1 and Halts Oxidative Stress-Mediated Cognitive Impairment in a Male Aging Mouse Model. Cells. 2019; 8:928. https://doi.org/10.3390/cells8080928 [PubMed]
22. Caito S, Rajendrasozhan S, Cook S, Chung S, Yao H, Friedman AE, Brookes PS, Rahman I. SIRT1 is a redox-sensitive deacetylase that is post-translationally modified by oxidants and carbonyl stress. FASEB J. 2010; 24:3145–59. https://doi.org/10.1096/fj.09-151308 [PubMed]
23. Laemmle A, Lechleiter A, Roh V, Schwarz C, Portmann S, Furer C, Keogh A, Tschan MP, Candinas D, Vorburger SA, Stroka D. Inhibition of SIRT1 impairs the accumulation and transcriptional activity of HIF-1α protein under hypoxic conditions. PLoS One. 2012; 7:e33433. https://doi.org/10.1371/journal.pone.0033433 [PubMed]
24. Alvarez-Garcia O, Matsuzaki T, Olmer M, Masuda K, Lotz MK. Age-related reduction in the expression of FOXO transcription factors and correlations with intervertebral disc degeneration. J Orthop Res. 2017; 35:2682–91. https://doi.org/10.1002/jor.23583 [PubMed]
25. Zhang J, Chen Y, Liu C, Li L, Li P. N1-Methylnicotinamide Improves Hepatic Insulin Sensitivity via Activation of SIRT1 and Inhibition of FOXO1 Acetylation. J Diabetes Res. 2020; 2020:1080152. https://doi.org/10.1155/2020/1080152 [PubMed]
26. Yao C, Cheng Y, Zhu Q, Lv Z, Kong L, Fang M. Clinical Evidence for the Effects of Manual Therapy on Cancer Pain: A Systematic Review and Meta-Analysis. Evid Based Complement Alternat Med. 2021; 2021:6678184. https://doi.org/10.1155/2021/6678184 [PubMed]
27. Jiang S, Zhang H, Fang M, Zhang Y, Lu N, Zhu Q, Cheng Y, Ai J, Zhou N, Li J, Fang L, Yao F. Analgesic effects of Chinese Tuina massage in a rat model of pain. Exp Ther Med. 2016; 11:1367–74. https://doi.org/10.3892/etm.2016.3055 [PubMed]
28. Zheng Z, Wang J, Gao Q, Hou J, Ma L, Jiang C, Chen G. Therapeutic evaluation of lumbar tender point deep massage for chronic non-specific low back pain. J Tradit Chin Med. 2012; 32:534–7. https://doi.org/10.1016/s0254-6272(13)60066-7 [PubMed]
29. Li T, Li X, Huang F, Tian Q, Fan ZY, Wu S. Clinical Efficacy and Safety of Acupressure on Low Back Pain: A Systematic Review and Meta-Analysis. Evid Based Complement Alternat Med. 2021; 2021:8862399. https://doi.org/10.1155/2021/8862399 [PubMed]
30. Zhang Y, Tang S, Chen G, Liu Y. Chinese massage combined with core stability exercises for nonspecific low back pain: a randomized controlled trial. Complement Ther Med. 2015; 23:1–6. https://doi.org/10.1016/j.ctim.2014.12.005 [PubMed]
31. Vincent K, Mohanty S, Pinelli R, Bonavita R, Pricop P, Albert TJ, Dahia CL. Aging of mouse intervertebral disc and association with back pain. Bone. 2019; 123:246–59. https://doi.org/10.1016/j.bone.2019.03.037 [PubMed]
32. Van Pelt DW, Confides AL, Abshire SM, Hunt ER, Dupont-Versteegden EE, Butterfield TA. Age-related responses to a bout of mechanotherapy in skeletal muscle of rats. J Appl Physiol (1985). 2019; 127:1782–91. https://doi.org/10.1152/japplphysiol.00641.2019 [PubMed]
33. Salgado AS, Stramosk J, Ludtke DD, Kuci AC, Salm DC, Ceci LA, Petronilho F, Florentino D, Danielski LG, Gassenferth A, Souza LR, Rezin GT, Santos AR, et al. Manual Therapy Reduces Pain Behavior and Oxidative Stress in a Murine Model of Complex Regional Pain Syndrome Type I. Brain Sci. 2019; 9:197. https://doi.org/10.3390/brainsci9080197 [PubMed]
34. Kolberg C, Horst A, Moraes MS, Duarte FC, Riffel AP, Scheid T, Kolberg A, Partata WA. Peripheral oxidative stress blood markers in patients with chronic back or neck pain treated with high-velocity, low-amplitude manipulation. J Manipulative Physiol Ther. 2015; 38:119–29. https://doi.org/10.1016/j.jmpt.2014.11.003 [PubMed]
35. Hartvigsen J, Hancock MJ, Kongsted A, Louw Q, Ferreira ML, Genevay S, Hoy D, Karppinen J, Pransky G, Sieper J, Smeets RJ, Underwood M, and Lancet Low Back Pain Series Working Group. What low back pain is and why we need to pay attention. Lancet. 2018; 391:2356–67. https://doi.org/10.1016/S0140-6736(18)30480-X [PubMed]
36. de Schepper EI, Damen J, van Meurs JB, Ginai AZ, Popham M, Hofman A, Koes BW, Bierma-Zeinstra SM. The association between lumbar disc degeneration and low back pain: the influence of age, gender, and individual radiographic features. Spine (Phila Pa 1976). 2010; 35:531–6. https://doi.org/10.1097/BRS.0b013e3181aa5b33 [PubMed]
37. Zhang Q, Li J, Li Y, Che H, Chen Y, Dong J, Xian CJ, Miao D, Wang L, Ren Y. Bmi deficiency causes oxidative stress and intervertebral disc degeneration which can be alleviated by antioxidant treatment. J Cell Mol Med. 2020; 24:8950–61. https://doi.org/10.1111/jcmm.15528 [PubMed]
38. Yokozeki Y, Uchida K, Kawakubo A, Nakawaki M, Okubo T, Miyagi M, Inoue G, Itakura M, Sekiguchi H, Takaso M. TGF-β regulates nerve growth factor expression in a mouse intervertebral disc injury model. BMC Musculoskelet Disord. 2021; 22:634. https://doi.org/10.1186/s12891-021-04509-w [PubMed]
39. Nakawaki M, Uchida K, Miyagi M, Inoue G, Kawakubo A, Satoh M, Takaso M. Changes in Nerve Growth Factor Expression and Macrophage Phenotype Following Intervertebral Disc Injury in Mice. J Orthop Res. 2019; 37:1798–804. https://doi.org/10.1002/jor.24308 [PubMed]
40. Dryden T, Sumpton B, Shipwright S, Kahn J, Reece BF. Massage therapy and canadians’ health care needs 2020: proceedings of a national research priority setting summit. Int J Ther Massage Bodywork. 2014; 7:3–10. https://doi.org/10.3822/ijtmb.v7i1.228 [PubMed]
41. Taavoni S, Darsareh F, Joolaee S, Haghani H. The effect of aromatherapy massage on the psychological symptoms of postmenopausal Iranian women. Complement Ther Med. 2013; 21:158–63. https://doi.org/10.1016/j.ctim.2013.03.007 [PubMed]
42. Donoyama N, Suoh S, Ohkoshi N. Effectiveness of Anma massage therapy in alleviating physical symptoms in outpatients with Parkinson’s disease: a before-after study. Complement Ther Clin Pract. 2014; 20:251–61. https://doi.org/10.1016/j.ctcp.2014.07.010 [PubMed]
43. Rodríguez-Mansilla J, González-López-Arza MV, Varela-Donoso E, Montanero-Fernández J, Jiménez-Palomares M, Garrido-Ardila EM. Ear therapy and massage therapy in the elderly with dementia: a pilot study. J Tradit Chin Med. 2013; 33:461–7. https://doi.org/10.1016/s0254-6272(13)60149-1 [PubMed]
44. Field T. Massage therapy research review. Complement Ther Clin Pract. 2016; 24:19–31. https://doi.org/10.1016/j.ctcp.2016.04.005 [PubMed]
45. Herman JP, McKlveen JM, Ghosal S, Kopp B, Wulsin A, Makinson R, Scheimann J, Myers B. Regulation of the Hypothalamic-Pituitary-Adrenocortical Stress Response. Compr Physiol. 2016; 6:603–21. https://doi.org/10.1002/cphy.c150015 [PubMed]
46. Shipton EA. Physical Therapy Approaches in the Treatment of Low Back Pain. Pain Ther. 2018; 7:127–37. https://doi.org/10.1007/s40122-018-0105-x [PubMed]
47. Yao C, Zhao C, Zhang S, Liu S. Effect of Moxibustion on Testosterone Secretion and Apoptosis of Spermatogenic Cells in Aging Rats. Evid Based Complement Alternat Med. 2019; 2019:5186408. https://doi.org/10.1155/2019/5186408 [PubMed]
48. Yu CC, Wang J, Ye SS, Gao S, Li J, Wang L, Jiang T, Wang XS, Li BC, Shu Q, Lu W, Du YJ, Kong LH. Preventive Electroacupuncture Ameliorates D-Galactose-Induced Alzheimer's Disease-Like Pathology and Memory Deficits Probably via Inhibition of GSK3β/mTOR Signaling Pathway. Evid Based Complement Alternat Med. 2020; 2020:1428752. https://doi.org/10.1155/2020/1428752 [PubMed]
49. Wen T, Fan X, Li M, Han J, Shi X, Xing L. Changes of metallothionein 1 and 3 mRNA levels with age in brain of senescence-accelerated mice and the effects of acupuncture. Am J Chin Med. 2006; 34:435–47. https://doi.org/10.1142/S0192415X06003977 [PubMed]
50. Akeda K, Yamada T, Inoue N, Nishimura A, Sudo A. Risk factors for lumbar intervertebral disc height narrowing: a population-based longitudinal study in the elderly. BMC Musculoskelet Disord. 2015; 16:344. https://doi.org/10.1186/s12891-015-0798-5 [PubMed]
51. Karabulut Ö, Tuncer MC, Karabulut Z, Açlkgöz A, Hatipoglu ES, Akkus Z. Relationship between radiographic features and bone mineral density in elderly men. Folia Morphol (Warsz). 2010; 69:170–6. [PubMed]
52. Ngo K, Patil P, McGowan SJ, Niedernhofer LJ, Robbins PD, Kang J, Sowa G, Vo N. Senescent intervertebral disc cells exhibit perturbed matrix homeostasis phenotype. Mech Ageing Dev. 2017; 166:16–23. https://doi.org/10.1016/j.mad.2017.08.007 [PubMed]
53. Oh JH, Kim EY, Nam TJ. Phycoerythrin-Derived Tryptic Peptide of a Red Alga Pyropia yezoensis Attenuates Glutamate-Induced ER Stress and Neuronal Senescence in Primary Rat Hippocampal Neurons. Mol Nutr Food Res. 2018; 62:e1700469. https://doi.org/10.1002/mnfr.201700469 [PubMed]
54. Gruber HE, Ingram JA, Norton HJ, Hanley EN Jr. Senescence in cells of the aging and degenerating intervertebral disc: immunolocalization of senescence-associated beta-galactosidase in human and sand rat discs. Spine (Phila Pa 1976). 2007; 32:321–7. https://doi.org/10.1097/01.brs.0000253960.57051.de [PubMed]
55. Tang N, Dong Y, Chen C, Zhao H. Anisodamine Maintains the Stability of Intervertebral Disc Tissue by Inhibiting the Senescence of Nucleus Pulposus Cells and Degradation of Extracellular Matrix via Interleukin-6/Janus Kinases/Signal Transducer and Activator of Transcription 3 Pathway. Front Pharmacol. 2020; 11:519172. https://doi.org/10.3389/fphar.2020.519172 [PubMed]
56. Cheng X, Ni B, Zhang F, Hu Y, Zhao J. High Glucose-Induced Oxidative Stress Mediates Apoptosis and Extracellular Matrix Metabolic Imbalances Possibly via p38 MAPK Activation in Rat Nucleus Pulposus Cells. J Diabetes Res. 2016; 2016:3765173. https://doi.org/10.1155/2016/3765173 [PubMed]
57. Yin H, Wang K, Das A, Li G, Song Y, Luo R, Cheung JP, Zhang T, Li S, Yang C. The REDD1/TXNIP Complex Accelerates Oxidative Stress-Induced Apoptosis of Nucleus Pulposus Cells through the Mitochondrial Pathway. Oxid Med Cell Longev. 2021; 2021:7397516. https://doi.org/10.1155/2021/7397516 [PubMed]
58. Sun Y, Li X, Yang X, Chen B, Zhang W. Small Extracellular Vesicles Derived from Adipocytes Attenuate Intervertebral Disc Degeneration in Rats by Rejuvenating Senescent Nucleus Pulposus Cells and Endplate Cells by Delivering Exogenous NAMPT. Oxid Med Cell Longev. 2021; 2021:9955448. https://doi.org/10.1155/2021/9955448 [PubMed]
59. Seo BR, Payne CJ, McNamara SL, Freedman BR, Kwee BJ, Nam S, de Lázaro I, Darnell M, Alvarez JT, Dellacherie MO, Vandenburgh HH, Walsh CJ, Mooney DJ. Skeletal muscle regeneration with robotic actuation-mediated clearance of neutrophils. Sci Transl Med. 2021; 13:eabe8868. https://doi.org/10.1126/scitranslmed.abe8868 [PubMed]
60. Austin PJ, Wu A, Moalem-Taylor G. Chronic constriction of the sciatic nerve and pain hypersensitivity testing in rats. J Vis Exp. 2012; 3393. https://doi.org/10.3791/3393 [PubMed]
61. Li Y, Wang L, Zhang G, Qiao X, Zhang M. SIRT1 Mediates Neuropathic Pain Induced by Sciatic Nerve Chronic Constrictive Injury in the VTA-NAc Pathway. Pain Res Manag. 2020; 2020:4245968. https://doi.org/10.1155/2020/4245968 [PubMed]
62. Vo N, Seo HY, Robinson A, Sowa G, Bentley D, Taylor L, Studer R, Usas A, Huard J, Alber S, Watkins SC, Lee J, Coehlo P, et al. Accelerated aging of intervertebral discs in a mouse model of progeria. J Orthop Res. 2010; 28:1600–7. https://doi.org/10.1002/jor.21153 [PubMed]
63. Bedore J, Sha W, McCann MR, Liu S, Leask A, Séguin CA. Impaired intervertebral disc development and premature disc degeneration in mice with notochord-specific deletion of CCN2. Arthritis Rheum. 2013; 65:2634–44. https://doi.org/10.1002/art.38075 [PubMed]

Research Paper Volume 14, Issue 5 pp 2400—2417