Evidence-Based Complementary and Alternative Medicine

Volume 2013 (2013), Article ID 890579, 7 pages

http://dx.doi.org/10.1155/2013/890579

Research Article

Gene Expression Profiles at Moxibustioned Site (ST36): A Microarray Analysis

Hai-Yan Yin,1,2,3 Yong Tang,1,2,3 Sheng-Feng Lu,⁴ Ling Luo,¹ Jia-Ping Wang,¹ Xu-Guang Liu,¹ and Shu-Guang Yu1,5

¹Acupuncture & Tuina School, Chengdu University of Traditional Chinese Medicine, Chengdu 610075, China

²Joint Laboratory of Biochip between Chengdu University of Traditional Chinese Medicine and CapitalBio Co. Ltd., 37 Shi-er Qiao Road, Chengdu 610075, China

³Key Laboratory for Acupuncture & Chronobiology of Sichuan Province, Chengdu 610075, China

⁴Acupuncture & Tuina School, Nanjing University of Traditional Chinese Medicine, Nanjing 210029, China

⁵Laboratory for Acupuncture & Systematic Biology, Chengdu University of Traditional Chinese Medicine, Chengdu 610075, China

Received 6 June 2013; Accepted 2 July 2013

Academic Editor: Yong-Qing Yang

Copyright © 2013 Hai-Yan Yin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

As a major alternative therapy in Traditional Chinese Medicine, it has been demonstrated that moxibustion could generate a series of molecular events in blood, spleen, and brain, and so forth. However, what would happen at the moxibustioned site remained unclear. To answer this question, we performed a microarray analysis with skin tissue taken from the moxibustioned site also Zusanli acupoint (ST36) where 15-minute moxibustion stimulation was administrated. The results exhibited 145 upregulated and 72 downregulated genes which responded immediately under physiological conditions, and 255 upregulated and 243 downregulated genes under pathological conditions. Interestingly, most of the pathways and biological processes of the differentially expressed genes (DEGs) under pathological conditions get involved in immunity, while those under physiological conditions are involved in metabolism.

1. Introduction

In acupuncture research, microarray analysis has been widely employed to uncover gene expression profiles at different tissues or organs [1–19]. Based on these gene expression profiles, researchers would be able to have the possibility to find out more potentially interesting targeted genes to conduct further experiment to explain the molecular events induced by acupuncture. Moxibustion, as one of the main therapies in acupuncture clinical practice, has been demonstrated to it could be useful for pain relief [20, 21] and generated a series of molecular events in blood [22, 23], spleen [24, 25], colonic mucosa [26], brain [27], and so forth, by utilizing moxa cone or stick to stimulate acupoint or some areas (also named moxibustioned site). However, none of gene expression profiles at moxibustioned site to date has been reported. Therefore, we proposed that moxibustion could, to a considerable extent, yield a great deal of differentially expressed genes (DEGs) at moxibustioned site, and we also anticipate to find out potential molecular targets to explain how moxibustion works at the stimulated site.

2. Material and Methods

2.1. Animals

Adult male Sprague-Dawley rats weighing 200–220 g obtained from Chengdu University of Traditional Chinese Medicine, Experimental Animal Centre, were utilized in this study. Maintained in animal room of automatically controlled day cycles (12 : 12 = light : dark cycle) at 24 ± 2°C, all rats were allowed to freely take food and water ad libitum and randomly assigned to the various experimental groups (, for each group). The experimental procedures were conducted in accordance with the National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals, and all experimental protocols were approved by the Animal Ethics Committee of Chengdu University of Traditional Chinese Medicine.

2.2. Experimental Design

In this study, we aimed to explore what would take place at the moxibustioned site in the view of potential molecular target under physiological and pathological conditions. Thus, we designed two different parts of microarray experiment. One is designed for uncovering the gene expression profile at physiological condition in which it consisted of 2 group: healthy control group (C) and healthy control with moxibustion stimulation group (CM). The other one is designed for revealing the gene expression profile at pathological condition in which it composed of model control group (M) and model control with moxibustion stimulation group (MM).

2.3. Intervention

2.3.1. Physiological Condition

Under physiological condition, the rats in CM group received moxibustion at the left acupoint Zusanli (ST36), at the depression below the knee from the anterior crest of the tibia [28] for 15 min. The moxibustion stimulation was manipulated with lighting moxa stick (length: 12 cm, diameter: 0.6 cm, Nanyang Hanyi Moxibustion Technology Development Co., Ltd., China) for 15 min (Figure 1). In case of skin burnt, the tip of moxa stick was kept about 2-3 cm from the skin.

Figure 1: Rat received moxibustion at ST36 (Zusanli).

2.3.2. Pathological Condition

Firstly, the pathological condition was established by injecting subcutaneously with 0.1 mL Freund’s Complete Adjuvant (FCA, Sigma, USA) into the plantar surface of the left hind paw of the rat [29]. The CFA injection immediately led to local inflammation, paw swelling and pain, which became apparent within 12 hours and persisted for at least 2 weeks after injection. In this experiment, the rats in MM group received moxibustion with the same procedure as mentioned above 1 week after injection.

2.4. RNA Extraction

Two hours after one time of moxibustion stimulation was completed, rats was euthanized by CO₂ inhalation. The cutaneous tissue (0.5 cm × 0.5 cm × 0.2 cm) located at moxibustioned site were immediately removed and preserved in RNAlater (Ambion, USA) to prevent RNA degradation. Total RNA was extracted using TRIzol Reagent (Invitrogen, USA) and purified with RNA clean-up Kit (MN, Germany) following the instructions of manufacturers, respectively. Total RNA was quantitated by spectrophotometry, and the integrity was assessed by formaldehyde denatured agarose gel electrophoresis.

2.5. Microarray Analysis

The microarray analysis service provided by CapitalBio Corporation (Beijing, China) was performed as described [30, 31]. Briefly, total RNA extracted from the samples was used to produce complementary RNA using in vitro transcription technique. Then cDNA was generated by reverse transcription and used as the template to synthesize the fluorescein-labeled cDNA by Klenow fragment polymerase. Universal rat reference RNA purchased from Stratagene was also labeled as common reference control. RNA from sample and common reference were fluorescently labeled by Cy5 or Cy3, respectively, and then were hybridized paired to 27K Rat Genome Array (CapitalBio, China). The array was comprised of 26,962 oligonucleotide probes covering 27,044 transcripts which represent about 22,012 genes. All arrays were scanned by LuxScan 10KA dual channel confocal laser scanner (CapitalBio, China). The obtained images were analyzed with LuxScan3.0 Image Analysis Software (CapitalBio, China), which employed the LOWESS normalization algorithm.

2.6. Data Analysis

2.6.1. Differentially Expressed Genes Selection

The detected signal intensities of all probes on the chip ≥400 were included for comparison analysis. We applied two-class unpaired algorithm in the Significant Analysis of Microarray software (SAM, Stanford) to identify significantly differentially expressed genes between CM and C groups, and MM and M groups. DEGs were determined with the threshold of false discovery rate, FDR ≤5% and fold change ≥2.0 or ≤0.5.

2.6.2. Pathway and Biological Processes Analysis of DEGs

We employed the online Molecule Annotation System (MAS) established by CapitalBio Corporation (http://bioinfo.capitalbio.com/mas3/) which integrated with KEGG and Gene Ontology (GO) database to perform pathway and GO Biological Process term enrichment analysis and calculate the statistical significance as described [32]. value <0.001 was considered statistically significant.

2.7. Real Time PCR Confirmation

To validate the expression patterns obtained from microarray data, we used quantitative real time polymerase chain reaction (qPCR) to detect the expression of four DEGs, Hspa1a, Mcpt8, Slpi, and Clqa, which were randomly selected from the 27K Rat Genome Array. Table 1 showed the primers designed for these genes and the housekeeping gene Gapdh. cDNA was prepared from DNase-treated total RNA using the First Strand SuperScript II Kit (Invitrogen, USA). qPCR was performed with DNA Master SYBR Green I Kit (Roche, Germany) and LightCycler machine (Roche, Germany) following the manufacturer’s protocols.

Table 1: The primer designed for validation.

Gene
Primer (5′-3′)
Temperature (°C)
Product size (bp)

Gapdh
FW: CCTTGTAAGGGCAAAACCAA

RV: ATGGCCTTCCGTGTTCCTAC
59
156

Hspa1a
FW: GGTGAACTACAAGGGCGAGA

RV: GCTGCGAGTCGTTGAAGTAG
58
152

Mcpt8
FW: CCAGGTCATCGCTGTTGTAAA

RV: CCCAGGTTTCACCCAGTCC
62
382

Slpi
FW: ACAGACAGGGGCTCTCTTGA

RV: CCTCCCAATAAGTGCCAGAA
60
216

Clqa
FW: AAGTGGGACCTTTGTCTGTCTATC

RV: CCCTGCTAACACCTGGAAGAG
59
108

3. Results

3.1. DEGs at Moxibustioned Site

Different numbers of DEGs at moxibustioned site were obtained from different condition. Under physiological condition, we obtained 145 upregulated and 72 downregulated DEGs (see Supplementary Table 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2013/890579). While under pathological condition, the results displayed 255 upregulated and 243 downregulated DEGs (Supplementary Table 2).

3.2. Enriched Pathways at Moxibustioned Site

Figure 2 showed us statistically significant pathways (value <0.001) at moxibustioned site. Under physiological condition (Figure 2(a)), it was found that 10 pathways were enriched based on all DEGs at moxibustioned Site. On the other hand, 21 enriched pathways were statistically significant under pathological condition (Figure 2(b)).

(a) Physiological condition

(b) Pathological condition

Figure 2: The statistically significant pathways ( value <0.001) involved in DEGs at moxibustioned site.

3.3. Enriched Biological Processes at Moxibustioned Site

From Figure 3, we would find out the biological processes with significantly statistical differences (value <0.001) at moxibustioned site. Under physiological condition (Figure 3(a)), it was found out that 9 biological processes were involved. Under pathological condition (Figure 3(b)), 29 biological processes were enriched.

(a) Physiological condition

(b) Pathological condition

Figure 3: The statistically significant biological processes ( value <0.001) involved in DEGs at moxibustioned site.

3.4. Validation of the 4 Selected Genes

To validate the results of the microarray, we selected 4 genes, Hspa1a, Mcpt8, Slpi, and Clqa, by qPCR. The results indicated that the expression levels of confirmed genes in microarray were similar with that in qPCR (Figure 4).

Figure 4: qPCR validation for 4 DEGs from microarray data.

4. Discussion

Here, we firstly reported the gene expression profiles at the moxibustioned site. The results imply that the moxibustioned site would also generate a setting of DEGs apart from those in other any tissue [16]. In other words, either under physiological condition or under pathological condition, we could find out a great number of DEGs responding to moxibustion stimulation at moxibustioned site. Based on these DEGs under different condition, we further seek out 22 coexpressed genes with similar expression tendency (Supplementary Table 3). To some extent, these genes (16 upnegulated and 6 downregulated genes) should be most likely to be considered as potential targets for continuing studies to determine which genes would be essential or critical in the role of moxibustion at moxibustioned site.

In this study, most (8/10) of the involved pathways at moxibustioned site were related to metabolism under physiological condition (Figure 2(a)). However, the significant pathways under pathological condition induced by FCA were most associated with immunity (Figure 2(b)), such as the pathway of antigen processing and presentation and natural killer cell-mediated cytotoxicity. It is suggested that different pathway at the moxibustioned site would get involved of the in different condition even with same moxibustion stimulation.

According to the biological processes analysis, we also could find out the difference from different state. Without FCA as the pathological stimulation, the biological progresses following moxibustion administration at moxibustioned site were most composed of oxidation reduction, potassium ion transport, and so forth. However, a variety of biological processes related to immunity, such as immune response, antigen processing and presentation, and antigen processing and presentation of peptide antigen via MHC class I, and were exhibited in this study.

Taken together, it seemed to be concluded that a series of molecular events would happen at moxibustioned site. Moreover, different pathways and biological progresses at moxibustioned site would be involved in different condition. However, this conclusion should be seriously taken for granted given the following limitations. Firstly, only one stimulation time point was used in this study. How about the time course or different time points, such as 5, 10, 15, 20, and 30 minutes, which were frequently practiced in moxibustion clinic, which needs to be answered in future study? Secondly, the pathological condition was induced by FCA injection in this study. To our knowledge, FCA injection will generate adjuvant arthritis through a series of immune actions [33, 34]. Moreover, previous studies also demonstrated that moxibustion would be useful to get better improved adjuvant arthritis [32, 35, 36]. Therefore, we cannot determine whether the DEGs at moxibustioned site under pathological condition would play important role in the moxibustion treatment of adjuvant arthritis. In this point of view, current data can only be used to explain what had happened at moxibustioned site under this condition. How about other pathological conditions? It is necessary to perform more researches to figure out the difference or similarity under different pathological conditions.

Additionally, in this study the tissue used for RNA extraction and microarray detection was taken from the skin of the moxibustioned site. So far, the skin has been also regarded as an important immune organ [37, 38] as well as a component of neuro-immuno-cutaneous system (NICS) [39–42]. In view of this aspect, we would be able to assure that the molecular event at moxibustioned site in this experiment will have the possibility to be applied to explain the initial mechanism of moxibustion activating the neuroimmune modulation which has been demonstrated in previous studies [43, 44]. This would be another potential mechanism of moxibustion apart from that it is currently considered as temperature-related and non-temperature-related work mechanisms [45].

5. Conclusions

The results suggested that a set of molecular events would have happened at moxibustioned site. Among those molecular events, different genes and different pathways and biological progresses at moxibustioned site would have got involved under different conditions.

Conflict of Interests

All authors manifest that there is no conflict of interests.

Authors’ Contribution

Hai-Yan Yin and Yong Tang contributed equally to this work.

Acknowledgments

This work was funded by the National Basic Research Program of China (973 Program, no. 2009CB522903), the National Natural Science Foundation of China (nos. 81173320, 81102667), Sichuan province Education Bureau (no. 12ZA037), and the Scientific Research Foundation of Sichuan Province of China for Returned Chinese Scholars.

References

S. K. Kim, J. Kim, E. Ko, et al., “Gene expression profile of the hypothalamus in DNP-KLH immunized mice following electroacupuncture stimulation,” Evidence-Based Complementary and Alternative Medicine, vol. 2011, Article ID 508689, 8 pages, 2011.

Y.-G. Choi, S. Yeo, Y.-M. Hong, and S. Lim, “Neuroprotective changes of striatal degeneration-related gene expression by acupuncture in an MPTP mouse model of Parkinsonism: microarray analysis,” Cellular and Molecular Neurobiology, vol. 31, no. 3, pp. 377–391, 2011.

Y.-G. Choi, S. Yeo, Y.-M. Hong, S.-H. Kim, and S. Lim, “Changes of gene expression profiles in the cervical spinal cord by acupuncture in an MPTP-intoxicated mouse model: microarray analysis,” Gene, vol. 481, no. 1, pp. 7–16, 2011.

Y.-G. Choi, S. Yeo, Y.-M. Hong, and S. Lim, “Neuroprotective changes of striatal degeneration-related gene expression by acupuncture in an MPTP mouse model of Parkinsonism: microarray analysis,” Cellular and Molecular Neurobiology, vol. 31, no. 3, pp. 377–391, 2011.

C. Tan, J. Wang, W. Feng, W. Ding, and M. Wang, “Preliminary correlation between warm needling treatment for knee osteoarthritis of deficiency-cold syndrome and metabolic functional genes and pathways,” Journal of Acupuncture and Meridian Studies, vol. 3, no. 3, pp. 173–180, 2010.

L.-H. Jiang and L.-L. Wang, “Gene chips-aided analysis on the profiles of hippocampal whole-genome expression in depression rats following electroacupuncture treatment,” Zhen Ci Yan Jiu, vol. 35, no. 2, pp. 83–89, 2010.

M. S. Hong, H.-K. Park, J.-S. Yang et al., “Gene expression profile of acupuncture treatment in 1-methyl-4-phenyl-1,2, 3,6- tetrahydropyridine-induced Parkinson's disease model,” Neurological Research, vol. 32, supplement 1, pp. S74–S78, 2010.

S.-H. Sohn, S. K. Kim, E. Ko et al., “The genome-wide expression profile of electroacupuncture in DNP-KLH immunized mice,” Cellular and Molecular Neurobiology, vol. 30, no. 4, pp. 631–640, 2010.

S. K. Kim, J. Y. Park, B. H. Koo et al., “Adenoviral gene transfer of acetylcholinesterase T subunit in the hypothalamus potentiates electroacupuncture analgesia in rats,” Genes, Brain and Behavior, vol. 8, no. 2, pp. 174–180, 2009.

H.-S. Shiue, Y.-S. Lee, C.-N. Tsai, Y.-M. Hsueh, J.-R. Sheu, and H.-H. Chang, “DNA microarray analysis of the effect on inflammation in patients treated with acupuncture for allergic rhinitis,” Journal of Alternative and Complementary Medicine, vol. 14, no. 6, pp. 689–698, 2008.

X.-Y. Wang, X.-L. Li, S.-Q. Hong, Y.-B. Xi-Yang, and T.-H. Wang, “Electroacupuncture induced spinal plasticity is linked to multiple gene expressions in dorsal root deafferented rats,” Journal of Molecular Neuroscience, vol. 37, no. 2, pp. 97–110, 2009.

Y.-Z. Gao, S.-Y. Guo, Q.-Z. Yin, T. Hisamitsu, and X.-H. Jiang, “An individual variation study of electroacupuncture analgesia in rats using microarray,” The American Journal of Chinese Medicine, vol. 35, no. 5, pp. 767–778, 2007.

L.-P. Yang, M.-C. Wang, W.-G. Liu, and M.-Q. Wang, “Effects of warming-needle therapy on gene expression pathways in the patient with knee osteoarthritis of deficiency-cold syndrome,” Zhongguo Zhen Jiu, vol. 27, no. 9, pp. 677–680, 2007.

M. Li and Y. Zhang, “Modulation of gene expression in cholesterol-lowering effect of electroacupuncture at Fenglong acupoint (ST40): a cDNA microarray study,” International Journal of Molecular Medicine, vol. 19, no. 4, pp. 617–629, 2007.

Y. Chae, H.-J. Park, D.-H. Hahm, S.-H. Yi, and H. Lee, “Individual differences of acupuncture analgesia in humans using cDNA microarray,” Journal of Physiological Sciences, vol. 56, no. 6, pp. 425–431, 2006.

H. G. . Wu, H. R. Liu, C. Zhao, et al., “Study on differentially expressed genes of ulcerative colitis in the rat treated by herbs-partitioned moxibustion,” Zhongguo Zhen Jiu, vol. 25, no. 5, pp. 359–365, 2005.

C. Keun Kim, S. C. Gi, D. O. Sang et al., “Electroacupuncture up-regulates natural killer cell activity: identification of genes altering their expressions in electroacupuncture induced up-regulation of natural killer cell activity,” Journal of Neuroimmunology, vol. 168, no. 1-2, pp. 144–153, 2005.

J.-C. Guo, H.-M. Gao, J. Chen et al., “Modulation of the gene expression in the protective effects of electroacupuncture against cerebral ischemia: a cDNA microarray study,” Acupuncture and Electro-Therapeutics Research, vol. 29, no. 3-4, pp. 173–186, 2004.

J. Ko, S. N. Doe, H. L. Young et al., “cDNA microarray analysis of the differential gene expression in the neuropathic pain and electroacupuncture treatment models,” Journal of Biochemistry and Molecular Biology, vol. 35, no. 4, pp. 420–427, 2002.

J.-H. Kim, H.-K. Kim, Y.-I. Park et al., “Moxibustion at ST36 alleviates pain in complete Freund's adjuvant-induced arthritic rats,” The American Journal of Chinese Medicine, vol. 34, no. 1, pp. 57–67, 2006.

L. Qi, H. R. Liu, T. Yi, et al., “Warming moxibustion relieves chronic visceral hyperalgesia in rats: relations to spinal dynorphin and orphanin-FQ system,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 920675, 10 pages, 2013.

H. Yamashita, Y. Ichiman, and Y. Tanno, “Changes in peripheral lymphocyte subpopulations after direct moxibustion,” The American Journal of Chinese Medicine, vol. 29, no. 2, pp. 227–235, 2001.

Y.-Y. Kung, F.-P. Chen, and S.-J. Hwang, “The different immunomodulation of indirect moxibustion on normal subjects and patients with systemic lupus erythematosus,” The American Journal of Chinese Medicine, vol. 34, no. 1, pp. 47–56, 2006.

G. S. Choi, J. B. Han, J. H. Park et al., “Effects of moxibustion to zusanli (ST36) on alteration of natural killer cell activity in rats,” The American Journal of Chinese Medicine, vol. 32, no. 2, pp. 303–312, 2004.

J.-B. Han, S.-D. Oh, K.-S. Lee et al., “The role of the sympathetic nervous system in moxibustion-induced immunomodulation in rats,” Journal of Neuroimmunology, vol. 140, no. 1-2, pp. 159–162, 2003.

E.-H. Zhou, H.-R. Liu, H.-G. Wu et al., “Down-regulation of protein and mRNA expression of IL-8 and ICAM-1 in colon tissue of ulcerative colitis patients by partition-herb moxibustion,” Digestive Diseases and Sciences, vol. 54, no. 10, pp. 2198–2206, 2009.

E.-H. Zhou, X.-M. Wang, G.-H. Ding et al., “Suspended moxibustion relieves chronic visceral hyperalgesia and decreases hypothalamic corticotropin-releasing hormone levels,” World Journal of Gastroenterology, vol. 17, no. 5, pp. 662–665, 2011.

S. Chen and S.-X. Ma, “Nitric oxide in the gracile nucleus mediates depressor response to acupuncture (ST36),” Journal of Neurophysiology, vol. 90, no. 2, pp. 780–785, 2003.

C. Stein, M. J. Millan, and A. Herz, “Unilateral inflammation of the hindpaw in rats as a model of prolonged noxious stimulation: alterations in behavior and nociceptive thresholds,” Pharmacology Biochemistry and Behavior, vol. 31, no. 2, pp. 445–451, 1988.

Y. Guo, H. Guo, L. Zhang et al., “Genomic analysis of anti-hepatitis B virus (HBV) activity by small interfering RNA and lamivudine in stable HBV-producing cells,” Journal of Virology, vol. 79, no. 22, pp. 14392–14403, 2005.

T. A. Patterson, E. K. Lobenhofer, S. B. Fulmer-Smentek et al., “Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project,” Nature Biotechnology, vol. 24, no. 9, pp. 1140–1150, 2006.

J.-Q. Fang, E. Aoki, A. Seto, Y. Yu, T. Kasahara, T. Hisamitsu, et al., “Influence of moxibustion on collagen-induced arthritis in mice,” In Vivo, vol. 12, no. 4, pp. 421–426, 1998.

C. M. Pearson and F. D. Wood, “Studies of arthritis and other lesions induced in rats by the injection of mycobacterial adjuvant. VII. Pathologic details of the arthritis and spondylitis,” The American Journal of Pathology, vol. 42, pp. 73–95, 1963.

C. M. Pearson, B. H. Waksman, and J. T. Sharp, “Studies of arthritis and other lesions induced in rats by injection of mycobacterial adjuvant. V. Changes affecting the skin and mucous membranes. Comparison of the experimental process with human disease,” The Journal of Experimental Medicine, vol. 113, no. 3, pp. 485–510, 1961.

Z. L. Tang, X. G. Song, F. Q. Zang, et al., “Study on the action mechanism of moxibustion in anti-inflammation and immunoregulation in rheumatoid arthritis rats,” Acupuncture Research, vol. 28, no. 4, pp. 292–298, 2003.

Y. . Shi and S. L. Zhou, “Effect of zusanli (ST36) acupoint on nerve-endocrine-immune network,” Journal of Jiangxi College of Traditional Chinese Medicine, vol. 15, no. 2, pp. 37–39, 2003.

J. C. Simon, “Skin as an immune organ,” Medizinische Monatsschrift für Pharmazeuten, vol. 18, no. 11, pp. 316–319, 1995.

J. K. Salmon, C. A. Armstrong, and J. C. Ansel, “The skin as an immune organ,” Western Journal of Medicine, vol. 160, no. 2, pp. 146–152, 1994.

L. Misery, “The neuro-immuno-cutaneous system (NICS),” Pathologie Biologie, vol. 44, no. 10, pp. 867–874, 1996.

J. D. Bos and M. L. Kapsenberg, “The skin immune system: progress in cutaneous biology,” Immunology Today, vol. 14, no. 2, pp. 75–78, 1993.

L. Misery, “The neuro-immuno-cutaneous system and ultraviolet radiation,” Photodermatology Photoimmunology and Photomedicine, vol. 16, no. 2, pp. 78–81, 2000.

F. O. Nestle, P. Di Meglio, J.-Z. Qin, and B. J. Nickoloff, “Skin immune sentinels in health and disease,” Nature Reviews Immunology, vol. 9, no. 10, pp. 679–691, 2009.

Y. Shi and S. L. Zhou, “Effect of Zusanli (ST36) Acupoint on nerve-endocrine-immune network,” Journal of Jiangxi College of Traditional Chinese Medicine, vol. 15, no. 2, pp. 37–39, 2003.

J. W. Li, J. M. Liu, Y. Y. Xiong, and S. Y. Xiang, “Advance in research on nerve-endocrine-immune network of rheumatoid arthritis mediated by acupuncture,” Traditional Chinese Medical Research, vol. 19, no. 3, pp. 57–60, 2006.

C. Jen-Hwey, “How does moxibustion possibly work?” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 198584, 8 pages, 2013.