Applying Statistical and Complex Network Methods to Explore the Key Signaling Molecules of Acupuncture Regulating Neuroendocrine-Immune Network

The mechanisms of acupuncture are still unclear. In order to reveal the regulatory effect of manual acupuncture (MA) on the neuroendocrine-immune (NEI) network and identify the key signaling molecules during MA modulating NEI network, we used a rat complete Freund's adjuvant (CFA) model to observe the analgesic and anti-inflammatory effect of MA, and, what is more, we used statistical and complex network methods to analyze the data about the expression of 55 common signaling molecules of NEI network in ST36 (Zusanli) acupoint, and serum and hind foot pad tissue. The results indicate that MA had significant analgesic, anti-inflammatory effects on CFA rats; the key signaling molecules may play a key role during MA regulating NEI network, but further research is needed.


Introduction
Acupuncture is a physical therapy of preventing or treating diseases by inserting needle into specific acupoints. As a kind of nonspecific physical stimulus, the effects of acupuncture are mediated by the regulatory systems in the body. This determines that the basic way of acupuncture effect is regulating the body's condition, with the characteristics of whole regulation.
Researches have shown that the body's inherent regulatory system is neuroendocrine-immune (NEI) network, including nervous system, endocrine system, and immune system, which is the biological basis to maintain the body's homeostasis. The three systems share common signaling molecules and their affiliated receptors, including some neuropeptides, neurotransmitters, cytokines, and hormones, and their receptors [1]. The cells in each system can secrete these signal molecules and at the same time the cells' surface has the moleculars' receptors. Hence, the common signaling molecules and their receptors constitute the molecular structural foundation of NEI network, being responsible for information communication and transmission between the three systems [1,2].
Some researches had demonstrated that acupuncture could regulate one system of the NEI network [3][4][5][6], in which most focus on the nervous system or immune system, but little studies pay attention to the regulatory effect of acupuncture on the whole NEI network [7]. Therefore, in order to reveal the regulatory effect of acupuncture on NEI network, in this study, we analyzed 55 common signaling molecules of NEI network in serum, supernatants from the ST36 acupoint, and hind foot pad tissue in rats with inflammatory pain 2.3. Inflammatory Pain Model. Rats per group were injected with either 0.1 mL CFA (Sigma, USA) or NS in the plantar surface of the right hind paw to induce intraplantar inflammation [10].

Measurement of Thermal
Hyperalgesia. Thermal hyperalgesia was assessed by hind paw withdrawal latency (PWL) to a noxious thermal stimulus using a plantar tester (BME-410C, Institute of Biological Medicine, Academy of Medical Science, China). Briefly, rats were placed in a clear plastic chamber (220 mm * 110 mm * 280 mm) and allowed to acclimatize for 30 min. A radiant heat stimulus was positioned under the glass floor directly beneath the right hind paw. When the rat withdrew its hind paw, we pressed the button to stop the heat stimulus, and the time was recorded as thermal PWL. Screening pain threshold before experiment: the rats with the PWL higher than 20 s or lower than 14 s were excluded from the experiment. A 30 s cut-off was used to prevent tissue injury. PWL was established by averaging the latency of 3 tests with a 5 min interval between each test. PWL was measured pre-CFA/normal saline injection and at D0 (after CFA injection), D1 (after MA), D7 (after MA), and D21 (after MA) at 14 o' clock to 16 o' clock.

Measurement of Hind Paw
Swelling. The swelling of rat's right hind paw was measured by volumetric method [11] with a self-made foot volume meter. The measured time was the same as PWL. The hind paw was immersed in a chamber containing PBS up to the boundary between hairy and nonhairy skin. The volume displacement represented the hind paw swelling and was determined by two observers. Paw volume was measured twice before CFA or NS injection (as basal paw volume) and at days 1, 7, and 21 after CFA or NS injection.
2.6. MA Treatment. Rats were immobilized in a holder and acupuncture needles (0.35 mm in diameter and 25 mm in length) were inserted to a depth of 5-7 mm at bilateral ST36 (Zusanli) acupoints. The needles were turned at a rate of 3 spins per second bidirectionally (1 spin consisted of clockwise rotation of 180 ∘ and a counterclockwise rotation of 180 ∘ ) for 2 min at Deqi, mild reinforcing and attenuating. The needles were manipulated every 5 min for a 30 min session. The manipulations were performed by the same person at 13 o' clock to 14 o' clock using metronome to keep the rhythm. MA treatment was given once a day for 7 consecutive days (day 1-day 7 after CFA injection) and then given every other day (day 8-day 21 after CFA injection), for a total of 14 sessions. In order to ensure the stability and repeatability of the manipulations, the operator practiced manipulation repeatedly at the ATP-II acupuncture manipulation parameter tester (which was manufactured by Shanghai University of Traditional Chinese Medicine Shang Xin Medical Technology Company) before and during the experiment. NS group and CFA group underwent grasping and fixation similar to CFA + MA group.

Sample Collection.
The hair on the right legs was removed from the skin using electronic hair clipper at 1 day before sample collection. After PWL measurement at D21 after CFA injection, rats were anesthetized with chloral hydrate (35 mg/kg, i.p.). Then the blood, local tissues in right ST36, and right hind footpad tissue were collected. The blood samples by abdominal aortic method was placed for 2 h under room temperature and then were centrifuged for 2000 rpm at 4 ∘ C for 10 min to get the serum. After blood collection, the tissues located in the right ST36 (1 cm in diameter and 0.5 cm thick, consisting of skin and subcutaneous and muscle tissues) and hind footpad tissue (skin and muscle) were collected immediately with an scapel. Next all tissues were triturated into homogenate with liquid nitrogen. 4 ∘ C and 2000 r/min for 10 min centrifugation was performed to get supernatants. The serum and supernatants were stored at −80 ∘ C before detecting the NEI common signaling molecules. The outline of the experimental protocol is summarized in Figure 1.  molecules completely, so we combined liquid chip, RIA, and ELISA to detect these signaling molecules.

Liquid
RIA and ELISA were conducted by Beijing Sinouk Institute of Biological Technology. Rat pituitary, rat stress hormone, rat thyroid, and rat neuropeptide liquid chip kit (Germany, Merck Millipore) were conducted by Beijing Institute of Hepatology. Bio-Plex Pro6 Rat Cytokine 24-plex Assay liquid chip kit (American, Bio-Rad) was conducted by Beijing Jian Yuan Wei Ye Technology Co., Ltd. Signaling molecules detection was carried out strictly according to the manufacturer's recommendations. The classification and detection methods of the signaling molecules were shown in Table 1.

Data
Analysis. The data about PWL and hind paw swelling were analyzed by statistical analysis method. The data about changes of common signaling molecules in NEI network of 3 kinds of samples induced by MA on day 21 were analyzed by 2 methods, including statistical analysis and complex network analysis.
2.9.1. Statistical Analysis. All statistical tests were conducted using SPSS 19.0 (SPSS Inc, Chicago, IL, USA) software. All statistical data were presented as the mean ± standard error. A value < 0.05 was considered to represent statistical significance.
PWL and hind paw volume data were analyzed using repeated measures analysis of variance (ANOVA), followed by Student-Newman-Keuls test which was used for post hoc analysis for differences between groups. If Mauchly's test of sphericity was not satisfied, One-Way ANOVA followed by LSD or Dunnett's T3 post hoc test were conducted. If data were not normally distributed or violated an assumption of a statistical test, they were transformed using commonly accepted methods or analyzed with a nonparametric test.
Signaling molecules were analyzed using One-Way ANOVA the same as above. The signaling molecules with statistical significance between groups (CFA group compared with NS group or CFA + EA group) were identified as the possible key signaling molecules in MA modulating the NEI network. All figures were generated using GraphPad Prism (GraphPad Software, La Jolla, CA).

Complex Network Analysis.
For more direct and visual analysis of the NEI changes induced by MA, all of the 55 common signaling molecules in the samples were analyzed by complex network methods. The complex network analysis method consists of 3 steps. In Step 1, based on the detection results, correlation coefficients between 55 NEI network signaling molecules of the serum, supernatants form the right ST36 point, and hind footpad tissue were calculated by Pearson correlation coefficient formula, and the correlation matrix was constructed by MATLAB software (Natick, America). In Step 2, the thresholds were used to filter the signaling molecules of the correlation coefficients ∈ [−1, −0.8], [0.8, 1]. In Step 3, screened signaling molecules were sorted by 3 complex network methods (i.e., node strength correlations, node degree, and node clustering coefficient) that could measure the importance of nodes in the network. Two or more methods in The first three nodes obtained by analysis with 2 or 3 methods mentioned above were considered as the key signaling molecules in the NEI network. The 3 complex network methods were as follows. A Node strength correlations: in complex network theory, signal molecules were viewed as nodes; the strength correlations were sum of absolute values of correlation coefficients of each node. Nodes were sorted from high to low by strength correlations, a node of higher strength correlations meaning it was more important in the network. B Node degree: the degree was the number of edges connected to each node. Nodes were sorted from large to small by degree, a node of larger degree meaning it was more important in the network. C Node clustering coefficient: it represented the possibility of connections between the other nodes that were connected to this one node. Nodes were sorted from low to high by the clustering coefficient, a node of lower strength correlations meaning the other nodes are with low possibility of connections except this one node, this also suggested that the node was important in the network. Figure 2, there were no differences in PWL under baseline condition (before CFA injection) among the three   resp.). However, the hind paw swelling significantly decreased ( < 0.05) in CFA + MA group (2.51 ± 0.17 ml) on day 21 compared with CFA group (2.84 ± 0.30 ml), when MA treatment persisted. These results suggested that MA could alleviate the CFA-induced hind paw swelling. resp.) were upregulated after MA. As shown in Figure 5, in the serum, the level of cytokine GRO/KC was upregulated and endocrine hormone PRL was downregulated after MA (515.95 ± 137.45 pg/ml, 6125.71 ± 4661.44 pg/ml, resp.). As shown in Figure 6, in the hind footpad tissue, neurotransmitter BDNF and cytokine, namely, RANTES and IL-2 (193.84 ± 45.32 pg/ml, 220.87±57.96 pg/ml, resp.), level were increased in the CFA group (43.77 ± 19.77 pg/ml). In addition, the level of endocrine hormone CRH was upregulated after MA (0.88 ± 0.06 pg/ml). These signaling molecules which were identified as statistically significant may be the possible key signaling molecules of NEI network.

The Changes of Common Signaling Molecules Based on
Complex Network Analysis. Correlation coefficients between the common signaling molecules of the NEI network in ST36 acupoint, serum, and lateral hind footpad tissue were calculated by Pearson correlation coefficient formula, and signaling moleculars association network (see Tables 2-7.). The first three nodes obtained by analysis with 2 or 3 methods mentioned above were considered as the key signaling molecules in the NEI network (see Table 8).

Discussion
In this study, the results showed that the PWL obviously decreased and hand paw swelling increased after CFA injections; MA could significantly increase the PWL and decrease hand paw swelling of the CFA rats; it indicated that MA had anti-inflammatory and antinociceptive effect on inflammatory pain in CFA rats. This is consistent with other studies [12][13][14].
The analysis of statistical results shows that, in CFA group, some common signaling molecules of NEI network in hind foot pad tissue were increased compared with NS Evidence-Based Complementary and Alternative Medicine group, including proinflammatory cytokines RANTES, IL-2, neuropeptide BDNF, indicating that hind footpad tissue after CFA injection was in inflammatory condition. Present studies indicated many similarities regarding the immunological changes and pathologic mechanisms existed between the CFA model and human rheumatoid arthritis (RA). So the CFA model is the widely used animal model for researching mechanisms and therapies of human RA [15]. Cytokines play an important role in pathogenesis of RA, for example, IL-2 could promote inflammatory response, activate macrophages and neutrophils, and inhibit Th2 lymphocyte proliferation in RA [16]. Other studies have found that the chemokine RANTES, secreted by monocytes/macrophages in the synovia of RA patients, was significantly increased [17] and could promote osteoclast formation [18], leading to increased inflammation. These were consistent with our results. The increased CRH in hind footpad tissue induced by MA should be further investigated.
Results about the common signaling molecules of NEI network in serum in CFA rats showed that GRO/KC increased and PRL decreased after MA treatment; that still needs to be further studied. Results about the common signaling molecules of NEI network in ST36 acupoint in CFA rats showed that MA induced the increasing expression of some hormones (TSH, melatonin, corticosterone, and FSH), proinflammatory cytokines (IL-1 , IL-6), and chemokine GRO/KC. As MA is a kind of noxious stimulation, so these changes of signaling molecules in ST36 acupoint maybe   transfer the acupuncture information [19]. Hi-Joon Par found that some neural and immune pathways, such as MAPK, B-cell receptor, T-cell receptor, and Toll-like receptor, in the local acupoint were involved in the anti-inflammatory and antinociceptive effect of acupuncture on inflammatory pain in CFA rats [20]. These findings support our results. The changes of these molecules may play a key role in the production and transmission of acupuncture information [21]. It still needs further research in the future. The analysis of complex network results shows that, in serum, IL-1 , IL-6, and IL-1 were the key signaling molecules in CFA rats, and T3, IL-13, and VEGF were the key signaling molecules in CFA + MA rats. Moreover, the key signaling molecules in hind footpad tissue of CFA rats were EPO, SP, and IL-1 and in CFA + MA rats were substance P, neurotensin, and oxytocin. These indicate that the key signaling molecules we acquired were consistent with the pathogenesis of RA and the analgesic effect and antiinflammatory mechanism of MA. Some researches reported that IL-2, IL-12, IL-1 , IL-1 , IL-6, and TNF-could induce acute phase reaction, stimulate the growth and differentiation of hematopoietic precursor cells, promote the proliferation of synovial fibroblasts, and cause joint damage [22] in the pathogenesis of RA. IL-13, produced by activated Th2 cells, had anti-inflammatory and immunomodulatory effects [23]. Therefore, these signaling molecules may play a key role in acupuncture regulating the NEI network but require further experimental confirmation. Other key signaling molecules identified, such as VEGF, EPO, SP, NT, and OT, need further experiments to explain the result. In the ST36 acupoint, the results of complex network showed that the key signaling molecules after MA were IFN-, MIP-3 , and TNF-, and the key signaling molecules identified by statistical methods included IL-1 , IL-6, and GRO/KC, which belong to proinflammatory cytokines or chemokines, so we hypothesize that acupuncture, as a physical stimulus, may cause an inflammatory reaction in the local acupoint, amplify the acupuncture information in cascade, and act on the NEI network, eventually producing acupuncture effects. The inflammatory response induced by acupuncture may be the starting point of acupuncture effect. In the future we will further research this part.
In this study, we detected 55 kinds of common signaling molecules of NEI network in three parts of CFA rats, including serum, ST36 acupoint, and hind footpad tissue. Since some signaling molecules in the three parts failed to be detected by the measurement techniques, the types and numbers of the signaling molecules detected in three parts were different. So it may affect the analysis results.

Conclusions
In conclusion, this study shows that MA has obvious analgesic and anti-inflammatory effects on CFA rats with inflammatory pain; the key signaling molecules of ST36 acupoint, serum, and hind footpad tissue acquired by statistical and complex network methods were all consistent with the pathogenesis of RA and the analgesic effect and anti-inflammatory mechanism of MA; these key signaling molecules in the three parts may play an important role in MA modulating NEI network; it still needs to be further studied in the future.