Gene-Level Regulation of Acupuncture Therapy in Spontaneously Hypertensive Rats: A Whole Transcriptome Analysis

Ma, Si-Ming; Yang, Jing-Wen; Tu, Jian-Feng; Yang, Na-Na; Du, Yu-Zheng; Wang, Xue-Rui; Wang, Lu; Huang, Jin; Liu, Cun-Zhi

doi:https://doi.org/10.1155/2019/9541079

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 9541079 | https://doi.org/10.1155/2019/9541079

Gene-Level Regulation of Acupuncture Therapy in Spontaneously Hypertensive Rats: A Whole Transcriptome Analysis

Si-Ming Ma,¹Jing-Wen Yang,¹Jian-Feng Tu,²Na-Na Yang,¹Yu-Zheng Du,³Xue-Rui Wang,²Lu Wang,¹Jin Huang,²and Cun-Zhi Liu¹

Academic Editor: Jana Sawynok

Received28 Sept 2018

Revised16 Dec 2018

Accepted09 Jan 2019

Published18 Feb 2019

Abstract

Hypertension is a global health problem. It has been reported that acupuncture at Taichong acupoints (LR3) decreases high blood pressure in spontaneously hypertensive rats. A transcriptome analysis can profile gene expression and its relationship with acupuncture. In this study, rats were treated with 2 weeks of acupuncture followed by regular recording of blood pressure (BP). The mRNA changes in the rostral ventrolateral medulla (RVLM) were evaluated to uncover the genetic mechanisms of acupuncture by using a whole transcript array (Affymetrix Rat Gene 1.0 ST array). BP measurements showed that acupuncture significantly decreased systolic blood pressure (SBP), mean arterial pressure (MAP), and heart rate (HR). In the bioinformatics results, 2371 differentially expressed genes (DEGs) were identified, where 83 DEGs were overlapped among Wistar-Kyoto rats (WKYs), spontaneously hypertensive rats (SHRs), and SHRs + acupuncture rats (SHRs+Acu). Gene ontology (GO) and pathway analysis revealed that 279 GO terms and 20 pathways with significant differences were related to oxidative stress, inflammation, and vascular endothelial function. In addition, coexpressed DEGs networks indicated that Cd4 and Il-33 might mediate the cascade of inflammation and oxidative stress responses, which could serve as a potential target of acupuncture treatment. In conclusion, our study demonstrated that acupuncture is a promising therapy for treating hypertension and could regulate multiple biological processes mainly involving oxidative stress, inflammation, and vascular endothelial function.

1. Introduction

Hypertension is a substantial public health problem, affecting over 1 billion people worldwide [1], and results in cardiovascular diseases and chronic kidney disease [2]. As a major independent risk factor for myocardial infarction, stroke, and end-stage renal disease, hypertension is a multifactor disorder that can be triggered by mechanisms in the central nervous systems (CNS) [3]. Neurogenic factors that induce hypertension might be a consequence of abnormal function of the autonomic nervous system, associated with stress, the renin–angiotensin system, insulin resistance, salt sensitivity, and genetic factors [4].

Among all brain regions, the rostral ventrolateral medulla (RVLM) is a critical node of blood pressure (BP) control [5]. Previous studies observed BP changes when RVLM neurons are bulk activated or inhibited with amino acids or GABA-mimetic compounds, while firing RVLM neurons with appropriate optogenetic or pharmacogenetic actuators could also produce the effects of BP fluctuation [5, 6]. Although other regions of the CNS contribute to the regulation of sympathetic tone, a series of experiments clearly showed that the RVLM is the main source of sympathetic activation that plays an irreplaceable role in inducing hypertension [7]. BP changes resulting primarily from variations in sympathetic nerve activity are presumably mediated by monosynaptic projections from the RVLM to preganglionic sympathetic neurons [8]. Abnormalities in the function and structure of the RVLM are closely related to the pathophysiological developing of hypertension [3].

Acupuncture was suggested to elicit integrative effects and modulate neural homeostasis in neurological disease, involving multiple neuroendocrine mechanisms [9, 10]. Previous clinical trials have indicated that acupuncture can effectively lower the systolic blood pressure (SBP) and diastolic blood pressure (DBP) in hypertensive patients [11, 12]. The mechanisms underlying antihypertensive effects of acupuncture involve apoptosis [13], endothelial dysfunction [14], apelinergic system [15], and the opioid system [16, 17]. Our team has recently demonstrated that high BP and oxidative stress in the RVLM of spontaneously hypertensive rats (SHRs) were decreased by acupuncture [18].

To validate the above theories, a bioinformatics study is needed with “big data” about genomics and its relationship with acupuncture, as well as downstream signaling pathways. Mapping and quantifying the transcriptome by microarray provide a digital measurement of the presence and prevalence of transcript [19]. In the RVLM, acupuncture’s effect on neurotransmitters has been profiled by a miRNA transcriptome study [20]. Given the inherent pleiotropic effects on gene expression of acupuncture, in the current study, a mRNA analysis was performed to uncover the potential genetic targets of acupuncture in antihypertensive effects.

2. Materials and Methods

2.1. Animals

All SHRs and Wistar-Kyoto rats (WKYs) were male adults purchased from Vital River Laboratory Animal Technology Co. Ltd (Beijing China). The animals were housed in 12 h light and 12 h dark cycles with a controlled temperature (24±0.5°C), with free access to food and tap water. All animal care and experiments conformed to the US National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee of the Capital Medical University.

2.2. Acupuncture Treatment

30 rats (12 weeks) were weighed and block-randomized with computer software (Graphpad Software, Inc., La Jolla, CA, USA) into 3 groups: Wistar-Kyoto rats (WKYs), spontaneously hypertensive rats (SHRs), and SHRs treated with acupuncture (SHRs+Acu). To eliminate the effect of constraint stress, an acupuncturist handled and fondled animals (1-2 hours/day, 2-3 days) before the acupuncture treatment. After the adaptive training, a 2-week period of acupuncture therapy was administrated to conscious rats at 3 p.m. to 4 p.m. each day. Before treatment, the bodies of rats were immobilized softly by an assistant and scrubbed with a 75% alcohol disinfectant at the acupoints. Then acupuncture was performed at bilateral LR3s with an intradermal needle by the acupuncturist. LR3 is located between the first and the second metatarsal bone of dorsal foot. The needles were inserted to a 5 mm depth and retained for 20 min before removal.

2.3. BP and Heart Rate Measurements

To investigate the effects of acupuncture on BP and sympathetic activity, a floating polyethylene catheter was inserted into the lower abdominal aorta, and the transmitter was implanted into the abdominal cavity. Seven days after the implantation, systolic blood pressure (SBP), mean arterial pressure (MAP), and heart rate (HR) were recorded for 2 h from 8 a.m. to 10 p.m. by a radio telemetry system (Data Sciences International, Saint Paul, Minnesota, USA) every 2 days. For BP and HR analysis, the system measured the mean values of these parameters for 2 hours in each measurement. The experimental staffs were blinded to group assignment.

2.4. Total RNA Extraction and mRNA Microarray

At the 14th day of acupuncture treatment, rats were euthanized with a sodium pentobarbital overdose (100mg/kg). Rat brains were removed and placed on dry ice. According to atlas, one 1 mm² region was obtained in the RVLM (bregma -12 to -12.96mm) from a 1 mm-thick slice at the level of brainstem via micropunch technique [21]. Total RNAs from RVLM were extracted using RNeasy Mini kit (Qiagen, Valencia, CA, USA). The amount and the integrity of total RNA were quantified and evaluated using a spectrophotometer (Beckman Coulter, Fullerton, CA, USA) and an Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA), respectively.

The total RNA was extracted using a TRIzol reagent (Invitrogen, Carlsbad, Canada). Hybridization and scanning of the chips were performed by Affymetrix Rat Gene 1.0 ST Array, which contained 59,302 gene-level probe sets.

2.5. Analysis of Differentially Expressed mRNA

To discern the genes that are differentially expressed, we chose p-value < 0.05 by ANOVA as well as the fold change <1.2 as the threshold screening among the 3 groups. The false discovery rate (FDR) was also computed in view that the smaller the FDR, the smaller the error in assessing the p-value.

2.6. Analysis of Series Test of Cluster (STC) of DEGs

Based on RVM (random variance model) corrective ANOVA, significant differentially expressed genes (DEGs) among the 3 groups were selected. According to the signal density tendency of genes in these differently treated groups, expression values were converted into ratio, and 16 unique model profiles were identified based on the actual number of genes assigned to each model. Fisher’s exact test and multiple comparison tests were performed to define the significant profiles.

2.7. qRT-PCR Analysis

In qRT-PCR analysis of validation, the total mRNA (2 μg) of samples was conducted to a reaction with AMV reverse transcriptase (Promega, Madison, WI) for first-strand cDNA synthesis. For each reaction, cDNA of sample was added to 2x SYBR master mix (Takara, Otsu, Shiga, Japan) and analyzed by using a BIORAD iCycler iQ5 (Bio-Rad, Hercules, CA). The sequences accessions of the primers are Cox5b (5′-GGAGGTGGTGTCCCTACTGA-3′ forward, 5′-GGAGGTGGTGTCCCTACTGA-3′ reverse), Sirt6 (5′-GCCGTCTGGTCATTGTCA-3′ forward; 5′-GCCGTCTGGTCATTGTCA-3′ reverse), Nf1 (5′-TTCGATACACTTGCGGAAAC-3′ forward; 5′-CACATTGGCAAGAGCCATAG-3′ reverse), Gabbr1 (5′-GGCTTTAGTCTGGGCTATG-3′ forward; 5′-GGCTTTAGTCTGGGCTATG-3′ reverse). Gene levels were normalized to that of β-actin. All samples were run in duplicate.

2.8. Gene Ontology (GO) Analysis and Pathway Analysis

Gene ontology (GO, www.geneontology.org/) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) pathway analysis were performed to identify the functions and associated enriched pathways of focused mRNAs. The z-score of threshold for mRNA expression was defined as ±2.0.

2.9. The Topological Analysis of Coexpressed DEGs

Gene coexpression networks were built to identify the interactions among 83 DEGs. We calculate the Pearson correlation and choose the significant correlation pairs to construct the networks. Degree centrality is defined as the link numbers that one gene has to the other. Besides, k-cores of networks were calculated to show the simplified topological characteristics.

2.10. Data and Statistical Analysis

Statistics in microarray and bioinformatic analyses were performed by 1-way analysis of variance (ANOVA) using Affymetrix® Expression Console™ TAC (Affymetrix® Expression Console™), followed by the least significant difference (LSD) test. The BP and HR data were analyzed by 2-way ANOVA with repeated measures, and followed by the Scheffe multiple-range test for post hoc assessment of individual means. Significant differences were considered at p < 0.05.

3. Results

3.1. The Effects of Acupuncture on BP and HR in Rats

To assess the antihypertensive effect of acupuncture in hypertension, radiotelemetry recording was performed before and during the experiment. As shown in Figure 1, acupuncture at LR3 in SHRs significantly decreased SBP and MAP from day 8 to the end of the treatment compared with SHRs (Figures 1(a) and 1(b)). HR also underwent a persistent decrease from day 4 in SHRs+Acu group (Figure 1(c)). But the normal BP and HR of WKYs were not altered by acupuncture.

(a)

(b)

(c)

3.2. DEGs Screening among 3 Groups

In the microarray analysis, genes expression that might be regulated by acupuncture treatment in the 3 groups was profiled. A threshold was defined as FD⩽0.001 and Ratio⩾1. As a result, a total of 2371 genes were acquired in all 3 groups (p < 0.05), the top 10 DEGs are Haghl, Fmo2, Mgst3, Dusp12, Cdkn1a, Glo1, Crot, Nmnat1, Eapp, Mettl7a (Table 1). To identify genes that were differentially regulated under the conditions of hypertension and acupuncture, a crossover comparison was performed in the following 3-paired groups: WKYs and SHRs, SHRs and SHRs+Acu, WKYs and SHRs+Acu. There was a number of 2169 (SHRs vs WKYs), 1221 (SHRs vs SHRs+Acu), and 2393 (WKYs vs SHRs+Acu) DEGs, as the Venn diagram demonstrates that they are coexpressed, respectively, in each pair of groups (Figure 2(b)), of which 83 DEGs were overlapped among 3 groups (Figure 2(a)). The number of up/downregulated DEGs was also calculated in each group (Figure 2(c)).

(a)

(b)

(c)

3.3. DEGs Clusters of Differentially Expression Patterns

To profile the DEGs and choose significant gene clusters with parallel expression patterns induced by acupuncture in SHR rats, 16 profiles with temporal expression patterns were defined in line with the trend of expression changes among the 3 groups. Temporal expression patterns and the significance of the genes were obtained by STEM software. The series test of cluster showed that 6 significant clusters were considered as potential profiles that could be affected by acupuncture (p < 0.05, Figure 3(a)). Notably, 226 genes were classified into the profile 2 and were gradually upregulated from WKYs to SHRs and SHRs+Acu. Meanwhile, profile 7 which included 197 genes was consistently decreased from the WKYs to SHRs and SHRs+Acu (Figure 3(b)).

(a)

(b)

(c)

3.4. Validation of Clusters Profiling Results by qRT-PCR

To verify the microarray profiling results, four significantly expressed genes that concern cardiovascular disease, Cox5b and Sirt6 in profile 2, Nf1 and Gabbr1 in profile 7, were evaluated by qRT-PCR after gene cluster analysis. The quantitative results of each expressed mRNA were consistent with the variation of its original clusters (Figure 3(c)).

3.5. GO and Pathway Analysis

To examine the characteristics of DEGs, functional classification of intersecting genes was performed using the gene ontology (GO) tool. A total of 1640 GO terms were identified, among which 279 GO terms were significant (p<0.05). In Figure 4(a), we show the 20 significantly overrepresented Go in terms of biological function. Top 10 GO terms are biological process, G-protein receptor signal transduction, Rho protein signal transduction, glutathione metabolic process, response to interleukin-1, regulation of glutamate secretion, positive regulation of vascular endothelial growth factor production, negative regulation of neuron apoptotic process, positive regulation of nuclear factor-kappaB (NF-kappaB) transcription factor, and toll-like receptor 4 (TLR4) signaling pathway. Accordingly, several representative genes (Fgfr2, Ptgsl, Ripk2, Vegfa, Nf1, Gabbr1, Gclm, Grm1) within specific GO terms are listed in Table 2.

(a)

(b)

To further understand the biological functions of the differentially expressed genes, we performed a KEGG pathway enrichment analysis on the basis of the intersection of genes (Figure 4(b)). A total of 177 pathways were identified, among which 20 pathways were significant (p<0.05). The top 10 significantly signaling pathways were, respectively, associated with metabolic pathways, a mitogen-activated protein kinases (MAPK) signaling pathway, insulin signaling pathway, focal adhesion, endocytosis, apoptosis, glutathione metabolism, Hippo signaling pathway, and inositol phosphate metabolism, and protein processing in endoplasmic reticulum was identified.

3.6. The Topological Analysis of Coexpressed DEGs

The dynamic real network was built to manifest the interaction of DEGs according to the data of a GO analysis (absolute value of interaction⩾0.800) (Figure 5(a)). We found that Cd4 and Il3 were connected with 39 genes (Degree =39) and occupied the leading position of modulation (k-core =12) in gene networks (Figure 5(b)). Since the expression of these 2 genes has been widely suggested to be of great importance to the formation of hypertension, we hypothesized that they may be involved in the antihypertensive effect of acupuncture.

(a)

(b)

4. Discussion

Acupoints selecting is essential for acquiring therapeutic effect in acupuncture therapy [22]. LR3 proved to be effective for treating hypertension. In stress-induced prehypertension model, needling at LR3 and Quchi (LI11) could reduce BP by targeting the gene expression of heart and hypothalamus [23, 24]. Consistent with previous studies, we selected LR3 and found that acupuncture significantly decreased the SBP, MBP, and HR in SHRs compared with WKYs, without any impact on normal BP. Moreover, some imaging studies showed that the effects of acupuncture on BP might involve altered cerebral or urinary metabolism [25, 26]. The molecular mechanism may be due to an increase in antioxidant enzyme in medulla according to a proteomic study [27]. Genetic information from in situ study was yielded by the microarray approach that was performed in the RVLM, and demonstrated that only a minority out of the 2371 DEGs were characterized as being associated with hypertension. Our analysis highlights once again that specific genes and biological processes may constitute the therapeutic effects of acupuncture on development of neurogenic hypertension. Of note, the bioinformatics results implicated that a large number of pathways and GO terms are associated with oxidative stress and inflammation. For instance, MAPK signaling pathway, JNK or ERK cascade, and glutathione metabolism are acknowledged as the fundamental processes mediating oxidative stress responses, while NF-kappa B and Toll-like receptor signaling pathways, and interleukin-1 could be the direct mediator of cellular inflammation or promotor of oxidative stress [28]. Notwithstanding the disperse evidence, these results provide a framework for further validation.

4.1. The BP Regulation of Acupuncture Might Be Involved in Multiple Gene Clusters

Cluster analysis was used to find out the transcriptional evidence of acupuncture’s effect on BP. Among 6 significant profiles, interestingly, the descending of profile 7 is mirrored by the ascending tendency of profile 2. We hypothesized that the phenomenon of their negative correlation might reflect the functional relation between 2 profiles in signaling pathway of hypertension.

To further answer the questions, 4 representative genes were selected (Cox5b, Sirt6, Nf1, and Gabbr1) to perform qRT-PCR due to their high correlation with oxidative stress or hypertension. Cytochrome c oxidase 5b (Cox5b) abnormal activity in hypertrophic hearts has been observed in SHRs [29]. Sirt6, predominantly expressed in nucleus, inhibits inflammation via reactive oxygen species (ROS) and Akt signaling pathway [30]. It is well known that Nf1 mutations lead to a common autosomal dominant disorder, neurofibromatosis type 1. The neurofibromatosis type 1 patients have an increased incidence of cardiovascular diseases, including obstructive vascular disorders and hypertension [31]. The gamma-aminobutyric acid B-type receptor unit 2 (Gabbr1) expressed in paraventricular nucleus (PVN) is important in controlling sympathetic activity, and its upregulation attenuated sympathoexcitation in chronic heart failure [32, 33]. Based on a series of elegant studies, such selections are relevant and sufficient for PCR validation. As expected, the expression scenario of these mRNA in 2 profiles is consistent with that in cluster analysis. This congruency suggests that the balance of some genetic subgroups, but not single gene, participates in maintaining blood pressure and could be regulated by acupuncture at LR 3.

GO analysis exhibited some redox-related genes that might be regulated by acupuncture. For example, Ripk2, also known as Rip2, is a kinase participating in inflammatory response signaling. Knocking down Ripk2 could reduce phosphorylation of p38 MAPK, ERK, and IkappaB Alpha and result in a decrease of IL-12 [34]. Moreover, suppression of endogenous Ripk2 significantly decreased apoptosis [35]. Fgfr2, a gene for encoding fibroblast growth factor receptor 2 (FGFR2), proved to regulate several growth-related signaling pathways in cancer [36]. The identification of specific mRNA sequence of Fgfr2 activates oncogenes through MAPK and PI3K/mTOR pathways in neuronal glial tumors [37]. Kobayashi T et al. demonstrated that Gclm knockout mice have long-term depletion of myocardial glutathione levels, which exacerbates myocardial oxidative stress [38] and dysfunction of the pressure-overloaded heart [39].

4.2. Oxidative Stress and Inflammation: The Key Combating in Antihypertensive Effect of Acupuncture

One main purpose of this study was to summarize the common and differential genetic factors that could be modulated by acupuncture through refining a huge dataset. Oxidative stress is an important contributor to hypertension, diabetes, and aging [40]. The ROS generated by NADPH oxidases (NOX) could mediate redox-sensitive signaling pathways through activating angiotensin II and cause endothelial dysfunction and vascular inflammation [41, 42]. This imbalance would be deleterious when inflammation is sustained. Although several signaling pathways affect the antihypertensive effects of treatment, the cross-talk between oxidative stress and inflammation might be a central mechanism of hypertension modulation for acupuncture.

Among 20 screened signaling pathways, MAPK signaling pathway is at the top of the KEGG list. As a critical inducer of oxidative stress, the activation of MAPK in RVLM leads to tissue and organ impairment in the form of phosphorylation [43]. Recent work of our group further revealed that acupuncture elicited antihypertensive effects and alleviated oxidative stress, especially via p38 MAPK and Erk1/2 signaling pathway in RVLM [18]. In addition, with the high response to oxidative stress, NF-KappaB signaling pathway is implicated in a series of metabolic syndromes. It has been well established that chronic inhibition of NF-KappaB in PVN attenuated sympathoexcitation by downregulating NOX and proinflammatory cytokine in hypertension [44]. This epigenetic changes induced by acupuncture now have been associated with MicroRNA-339 and Sirt2 (upstream of NFκB) [45].

Furthermore, 2 redox-sensitive pathways screened by microarray, including glutathione metabolism and the neurotrophin pathway, contribute to oxidative stress in hypertension. Glutathione is a major intracellular thiol-disulfide redox cofactor for many antioxidant enzymes and is considered to be a key node in combating increased blood pressure, not only for its capacity of free radical scavenging, but also by controlling nitric oxide (NO) bioavailability [46]. A proteomic study in hypertension indicates that glutathione S-transferase M5 could be upregulated by acupuncture [47]. Neurotrophin factors, such as brain-derived neurotrophin factor (BDNF), are acknowledged to be a sensitive marker of oxidative stress. They increase, prior to the BP rise, and enhance the survival of noradrenergic neurons hyperinnervated by sympathetic nerves [48, 49]. A miRNA study profiled that acupuncture could restore high BP through activating neurotrophin signaling pathway [20]. Therefore, our study reinforced the idea that acupuncture affects hypertension by regulating signaling pathways concerning oxidative stress.

4.3. Vascular Endothelial Dysfunction: Initiated by Oxidative Stress

Endothelial dysfunction induces the deterioration of endothelium-dependent vasodilation including the aberrant metabolism of nitric oxide and imbalance of vasomotor factors [50]. For example, as an angiogenic factor, vascular endothelial growth factor-A (VEGFA) regulates BP by promoting vasodilation. The increased value of VEGFA, accompanied by decreased NO, could be observed in hypertension. Ptgs1 is engaged in regulation of BP by encoding cyclooxygenase-1 (COX-1) [51, 52]. The COX-derived prostanoids further mediate vascular tone in programmed hypertension [19]. These emerging studies implicated the importance of endothelial dysfunction in hypertension, where the oxidative stress serves as the main factor of pathogenesis [28]. It has been evidenced that the overproduction of ROS is correlated with the dysfunctional endothelial NOS (eNOS) and results in the reduction of vascular relaxation [53]. During the course of pathological changes, acupuncture could attenuate BP by protecting endothelial function from oxidative stress in SHRs [14]. In accordance, GO annotation showed the beneficial effects of acupuncture in vascular endothelium and, additionally, implied that a tiny stimulation of acupuncture could cause multiple genes cascade due to the activation of oxidative stress in hypertension.

4.4. Unproved but Promising Pathways: All Roads Lead to Antioxidative Effect of Acupuncture

Our study also reported several redox-related signaling pathways that have not been mentioned in prior studies of the mechanisms of acupuncture for treating hypertension. First, increased BP is commonly accompanied by risk factors such as insulin resistance or dyslipidemia especially in SHRs and essential hypertension in human [54]. The characteristics of oxidative stress and the impaired PI3 pathway may favor the insulin resistance in rats and patients with hypertension [55, 56]. Second, the main function of the Hippo signaling pathway has been expanded to modulating cell proliferation, differentiation, and migration in many organs [57]. Its function of negative regulation of cardiomyocytes survival was proposed to be associated with inhibiting Yes-associated protein (YAP)-FoxO1 interaction and antioxidant gene expression [58]. Third, NOX4 upregulation would induce the unfolded protein response (UPR), causing endoplasmic reticulum (ER) stress in vascular smooth muscle cells, which was mentioned frequently in the last few years for its role in BP elevation and vascular injury in SHR [59]. Fourth, TGF-beta signaling pathway also activates NOX-induced oxidative stress and accelerates endothelial dysfunction in hypertension [60–62]. Thereby, these putative pathways could be some new targets for acupuncture treatment, and their function in acupuncture effectiveness and oxidative stress should be clarified by enough experiments.

4.5. Homeostasis: The Macroscopic Mechanism of Acupuncture Treatment

The topological analysis illustrated an unanticipated outcome; that is, Cd4 and Il33 take control of the gene network. Previous investigations indicated that the immune system mediates human essential hypertension via the activation of innate and adaptive immunity in periphery [63]. As a main T-regulatory lymphocytes (Treg), CD4⁺ lymphocytes that express choline acetyltransferase (ChAT) induce the vasodilation by relaying the signals from vagus nerve [64]. The gene expression of CD4⁺CDCD62 T helper cells would determine this distinct type of T-cell [65]. Furthermore, In vitro experiment identified a feedback loop, in which IL-33, in addition to its anti-inflammatory and protective effects in cardiovascular system, stimulates the expanding of tumorigenicity 2 (ST2)-expressing Treg during interactions with CD4⁺ [66].

In the CNS, however, the tight conjunction of blood brain barrier (BBB) makes it difficult for T-cells or IL-33 to penetrate into the brain parenchyma. A widely accepted hypothesis is that circumventricular organs (CVOs) such as the subfornical organ (SFO), which lack the normal BBB, are apt to be triggered by Ang II and induce the change of BBB filterability in other brain regions [67]. Indeed, T-lymphocyte infiltration was found in SFO tissue [68], and SFO-PVN-RVLM pathway induced by Ang II-ROS has been described in a few researches in neurogenic hypertension [67, 69]. This delineation of integrity that combines immunity, inflammation, and oxidative stress is consistent with the homeostasis theory in acupuncture treatment undoubtedly. We are cognizant that the expression of these genes might provide some small pieces of puzzle for mechanism of blood pressure maintaining and for mechanism of acupuncture therapeutic effects too.

5. Conclusion

A comprehensive bioinformatic analysis revealed that genes expression related to oxidative stress, inflammation, and vascular endothelial function and downstream signaling pathways it is involved in might be associated with the antihypertensive effects of acupuncture. Our present study enriched the understanding of the antihypertensive mechanisms of acupuncture. More studies are needed to address the relationship between these genes and hypertension, as well as the beneficial role of acupuncture in gene expression and pathway regulation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Si-Ming Ma analyzed the data and drafted most of the manuscript. Xue-Rui Wang, Jian-Feng Tu, and Yu-Zheng Du analyzed the results of the whole transcript array. Jing-Wen Yang, Lu Wang, Jin Huang and Na-Na Yang performed validation experiment. Cun-Zhi Liu conceived and designed the experiment and revised the paper with Si-Ming Ma.

Acknowledgments

This study was supported by Distinguished Young Scholars project of Beijing University of Chinese Medicine [Grant no. BUCM-2019-JCRC011].

Supplementary Materials

The supplementary information to support the results of transcriptome analysis is available in the online version of this paper. (Supplementary Materials)

References

K. T. Mills, J. D. Bundy, T. N. Kelly et al., “Global disparities of hypertension prevalence and control,” Circulation, vol. 134, no. 6, pp. 441–450, 2016.
View at: Publisher Site | Google Scholar
M. H. Forouzanfar, L. Alexander, H. R. Anderson et al., “Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks in 188 countries, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013,” Lancet (London, England), vol. 386, no. 10010, pp. 2287–2323, 2013.
View at: Google Scholar
A. F. Sved, S. Ito, and J. C. Sved, “Brainstem mechanisms of hypertension: role of the rostral ventrolateral medulla,” Current Hypertension Reports, vol. 5, no. 3, pp. 262–268, 2003.
View at: Publisher Site | Google Scholar
S. J. Mann, “Neurogenic essential hypertension revisited: The case for increased clinical and research attention,” American Journal of Hypertension, vol. 16, no. 10, pp. 881–888, 2003.
View at: Publisher Site | Google Scholar
N. Marina, A. P. L. Abdala, A. Korsak et al., “Control of sympathetic vasomotor tone by catecholaminergic C1 neurones of the rostral ventrolateral medulla oblongata,” Cardiovascular Research, vol. 91, no. 4, pp. 703–710, 2011.
View at: Publisher Site | Google Scholar
V. Geraldes, N. Goncalves-Rosa, B. Liu, J. F. R. Paton, and I. Rocha, “Essential role of RVL medullary neuronal activity in the long term maintenance of hypertension in conscious SHR,” Autonomic Neuroscience: Basic and Clinical, vol. 186, no. C, pp. 22–31, 2014.
View at: Publisher Site | Google Scholar
E. Colombari, M. A. Sato, S. L. Cravo, C. T. Bergamaschi, R. R. Campos Jr., and O. U. Lopes, “Role of the medulla oblongata in hypertension,” Hypertension, vol. 38, no. 3, Pt. 2, pp. 549–554, 2001.
View at: Publisher Site | Google Scholar
D. Les Brown and P. G. Guyenet, “Electrophysiological study of cardiovascular neurons in the rostral ventrolateral medulla in rats,” Circulation Research, vol. 56, no. 3, pp. 359–369, 1985.
View at: Publisher Site | Google Scholar
Z.-Q. Zhao, “Neural mechanism underlying acupuncture analgesia,” Progress in Neurobiology, vol. 85, no. 4, pp. 355–375, 2008.
View at: Publisher Site | Google Scholar
H. K. Shin, S.-W. Lee, and B. T. Choi, “Modulation of neurogenesis via neurotrophic factors in acupuncture treatments for neurological diseases,” Biochemical Pharmacology, vol. 141, pp. 132–142, 2017.
View at: Publisher Site | Google Scholar
L. A. Kalish, B. Buczynski, P. Connell et al., “Stop Hypertension with the Acupuncture Research Program (SHARP): clinical trial design and screening results,” Controlled Clinical Trials, vol. 25, no. 1, pp. 76–103, 2004.
View at: Publisher Site | Google Scholar
E. A. Macklin, P. M. Wayne, L. A. Kalish et al., “Stop Hypertension with the Acupuncture Research Program (SHARP): results of a randomized, controlled clinical trial,” Hypertension, vol. 48, no. 5, pp. 838–845, 2006.
View at: Publisher Site | Google Scholar
S.-N. Chang Lee, T.-J. Ho, M. A. Shibu et al., “Protective effects of electroacupuncture at LR3 on cardiac hypertrophy and apoptosis in hypertensive rats,” Acupuncture in Medicine: Journal of the British Medical Acupuncture Society, vol. 34, no. 3, pp. 201–208, 2016.
View at: Publisher Site | Google Scholar
S. B. Leung, H. Zhang, C. W. Lau, and Z.-X. Lin, “Attenuation of blood pressure in spontaneously hypertensive rats by acupuncture was associated with reduction oxidative stress and improvement from endothelial dysfunction,” Chinese Medicine, vol. 11, no. 1, article no. 38, 2016.
View at: Publisher Site | Google Scholar
C. Zhang, C. Xia, M. Jiang et al., “Repeated electroacupuncture attenuating of apelin expression and function in the rostral ventrolateral medulla in stress-induced hypertensive rats,” Brain Research Bulletin, vol. 97, pp. 53–62, 2013.
View at: Publisher Site | Google Scholar
M. Li, S. C. Tjen-A-Looi, Z.-L. Guo, and J. C. Longhurst, “Electroacupuncture modulation of reflex hypertension in rats: Role of cholecystokinin octapeptide,” American Journal of Physiology-Regulatory, Integrative and Comparative Physiology, vol. 305, no. 4, pp. R404–R413, 2013.
View at: Publisher Site | Google Scholar
M. Li, S. C. Tjen-A-Looi, Z.-L. Guo, and J. C. Longhurst, “Repetitive electroacupuncture attenuates cold-induced hypertension through enkephalin in the rostral ventral lateral medulla,” Scientific Reports, vol. 6, Article ID 35791, 2016.
View at: Publisher Site | Google Scholar
X. Wang, J. Yang, C. Ji et al., “Inhibition of NADPH oxidase–dependent oxidative stress in the rostral ventrolateral medulla mediates the antihypertensive effects of acupuncture in spontaneously hypertensive rats,” Hypertension, vol. 71, no. 2, pp. 356–365, 2017.
View at: Publisher Site | Google Scholar
W. B. Campbell and J. R. Falck, “Arachidonic acid metabolites as endothelium-derived hyperpolarizing factors,” Hypertension, vol. 49, no. 3, pp. 590–596, 2007.
View at: Publisher Site | Google Scholar
J.-Y. Wang, H. Li, C.-M. Ma I, J.-L. Wang, X.-S. Lai, and S.-F. Zhou, “MicroRNA profiling response to acupuncture therapy in spontaneously hypertensive rats,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 204367, 9 pages, 2015.
View at: Publisher Site | Google Scholar
T.-C. Yeh, C.-P. Liu, W.-H. Cheng et al., “Caffeine intake improves fructose-induced hypertension and insulin resistance by enhancing central insulin signaling,” Hypertension, vol. 63, no. 3, pp. 535–541, 2014.
View at: Publisher Site | Google Scholar
S. Lu, X. Cao, H. Ohara et al., “Common parameters of acupuncture for the treatment of hypertension used in animal models,” Journal of Traditional Chinese Medicine, vol. 35, no. 3, pp. 343–348, 2015.
View at: Google Scholar
Y. Guo, X. Xie, C. Guo, Z. Wang, and Q. Liu, “Effect of electro-acupuncture on gene expression in heart of rats with stress-induced pre-hypertension based on gene chip technology,” Journal of Traditional Chinese Medicine, vol. 35, no. 3, pp. 285–294, 2015.
View at: Google Scholar
X. Xie, Y. Guo, Q. Liu, Z. Wang, and C. Guo, “A Study on effect of electroacupuncture on gene expression in hypothalamus of rats with stress-induced prehypertension based on gene chip technology,” Evidence-Based Complementary and Alternative Medicine, vol. 2015, Article ID 621237, 12 pages, 2015.
View at: Publisher Site | Google Scholar
J. Li, Y. Wang, K. He et al., “Effect of acupuncture at LR3 on cerebral glucose metabolism in a rat model of hypertension: a 18F-FDG-PET study,” Evidence-Based Complementary and Alternative Medicine, vol. 2018, Article ID 5712857, 8 pages, 2018.
View at: Publisher Site | Google Scholar
Y. Wang, Y. Li, L. Zhou, and L. Guo, “Effects of acupuncture on the urinary metabolome of spontaneously hypertensive rats,” Acupuncture in Medicine, vol. 35, no. 5, pp. 374–382, 2017.
View at: Publisher Site | Google Scholar
X. Lai, J. Wang, N. R. Nabar et al., “Proteomic response to acupuncture treatment in spontaneously hypertensive rats,” PLoS ONE, vol. 7, no. 9, Article ID e44216, 2012.
View at: Publisher Site | Google Scholar
T. J. Kizhakekuttu and M. E. Widlansky, “Natural antioxidants and hypertension: promise and challenges,” Cardiovascular Therapeutics, vol. 28, no. 4, pp. e20–e32, 2010.
View at: Publisher Site | Google Scholar
W.-W. Kuo, C.-Y. Chu, C.-H. Wu et al., “The profile of cardiac cytochrome c oxidase (COX) expression in an accelerated cardiac-hypertrophy model,” Journal of Biomedical Science, vol. 12, no. 4, pp. 601–610, 2005.
View at: Publisher Site | Google Scholar
Y. He, Y. Xiao, X. Yang et al., “SIRT6 inhibits TNF-α-induced inflammation of vascular adventitial fibroblasts through ROS and Akt signaling pathway,” Experimental Cell Research, vol. 357, no. 1, pp. 88–97, 2017.
View at: Publisher Site | Google Scholar
K. Kobus, D. Hartl, C. E. Ott et al., “Double NF1 inactivation affects adrenocortical function in NF1Prx1 mice and a human patient,” PLoS ONE, vol. 10, no. 3, Article ID e0119030, 2015.
View at: Google Scholar
D. Li and H. Pan, “Role of gamma-aminobutyric acid (GABA)A and GABAB receptors in paraventricular nucleus in control of sympathetic vasomotor tone in hypertension,” The Journal of Pharmacology and Experimental Therapeutics, vol. 320, no. 2, pp. 615–626, 2006.
View at: Publisher Site | Google Scholar
R. Wang, Q. Huang, R. Zhou et al., “Sympathoexcitation in rats with chronic heart failure depends on homeobox D10 and MicroRNA-7b inhibiting GABBR1 translation in paraventricular nucleus,” Circulation: Heart Failure, vol. 9, no. 1, Article ID e002261, 2016.
View at: Google Scholar
N. Usluoglu, J. Pavlovic, K. Moelling, and G. Radziwill, “RIP2 mediates LPS-induced p38 and IκBα signaling including IL-12 p40 expression in human monocyte-derived dendritic cells,” European Journal of Immunology, vol. 37, no. 8, pp. 2317–2325, 2007.
View at: Publisher Site | Google Scholar
A. Krieg, G. Le Negrate, and J. C. Reed, “RIP2-β: A novel alternative mRNA splice variant of the receptor interacting protein kinase RIP2,” Molecular Immunology, vol. 46, no. 6, pp. 1163–1170, 2009.
View at: Publisher Site | Google Scholar
J. D. Holbrook, J. S. Parker, K. T. Gallagher et al., “Deep sequencing of gastric carcinoma reveals somatic mutations relevant to personalized medicine,” Journal of Translational Medicine, vol. 9, no. 1, p. 119, 2011.
View at: Google Scholar
P. Jain, L. F. Surrey, J. Straka et al., “Novel FGFR2-INA fusion identified in two low-grade mixed neuronal-glial tumors drives oncogenesis via MAPK and PI3K/mTOR pathway activation,” Acta Neuropathologica, vol. 136, no. 1, pp. 167–169, 2018.
View at: Google Scholar
T. Kobayashi, Y. Watanabe, Y. Saito et al., “Mice lacking the glutamate-cysteine ligase modifier subunit are susceptible to myocardial ischaemia-reperfusion injury,” Cardiovascular Research, vol. 85, no. 4, pp. 785–795, 2010.
View at: Publisher Site | Google Scholar
Y. Watanabe, K. Watanabe, T. Kobayashi et al., “Chronic depletion of glutathione exacerbates ventricular remodelling and dysfunction in the pressure-overloaded heart,” Cardiovascular Research, vol. 97, no. 2, pp. 282–292, 2013.
View at: Publisher Site | Google Scholar
G. R. Drummond, S. Selemidis, K. K. Griendling, and C. G. Sobey, “Combating oxidative stress in vascular disease: NADPH oxidases as therapeutic targets,” Nature Reviews Drug Discovery, vol. 10, no. 6, pp. 453–471, 2011.
View at: Publisher Site | Google Scholar
B. Lassègue, D. Sorescu, K. Szöcs et al., “Novel gp91^PHOX homologues in vascular smooth muscle cells: nox1 mediates angiotensin II-induced superoxide formation and redox-sensitive signaling pathways,” Circulation Research, vol. 88, no. 9, pp. 888–894, 2001.
View at: Publisher Site | Google Scholar
U. Landmesser, H. Cai, S. Dikalov et al., “Role of p47phox in vascular oxidative stress and hypertension caused by angiotensin II,” Hypertension, vol. 40, no. 4, pp. 511–515, 2002.
View at: Publisher Site | Google Scholar
S. H. H. Chan, K.-S. Hsu, C.-C. Huang, L.-L. Wang, C.-C. Ou, and J. Y. H. Chan, “NADPH oxidase-derived superoxide anion mediates angiotensin II-induced pressor effect via activation of p38 mitogen-activated protein kinase in the rostral ventrolateral medulla,” Circulation Research, vol. 97, no. 8, pp. 772–780, 2005.
View at: Publisher Site | Google Scholar
J. Qi, X. Yu, X. Shi et al., “NF-κB blockade in hypothalamic paraventricular nucleus inhibits high-salt-induced hypertension through NLRP3 and caspase-1,” Cardiovascular Toxicology, vol. 16, no. 4, pp. 345–354, 2016.
View at: Publisher Site | Google Scholar
J.-Y. Wang, H. Li, C.-M. Ma, J.-L. Wang, X.-S. Lai, and S.-F. Zhou, “Acupuncture may exert its therapeutic effect through microRNA-339/Sirt2/NFkappaB/FOXO1 axis,” BioMed Research International, vol. 2015, Article ID 249013, 9 pages, 2015.
View at: Publisher Site | Google Scholar
J. Robaczewska, K. Kedziora-Kornatowska, M. Kozakiewicz et al., “Role of glutathione metabolism and glutathione-related antioxidant defense systems in hypertension,” Journal of Physiology and Pharmacology, vol. 67, no. 3, pp. 331–337, 2016.
View at: Google Scholar
X. Zhang, B. Wu, K. Nie, Y. Jia, and J. Yu, “Effects of acupuncture on declined cerebral blood flow, impaired mitochondrial respiratory function and oxidative stress in multi-infarct dementia rats,” Neurochemistry International, vol. 65, pp. 23–29, 2014.
View at: Publisher Site | Google Scholar
D. Y. Atanasova and N. E. Lazarov, “Expression of neurotrophic factors and their receptors in the carotid body of spontaneously hypertensive rats,” Respiratory Physiology & Neurobiology, vol. 202, pp. 6–15, 2014.
View at: Publisher Site | Google Scholar
S.-H. Zhang and R. A. Rush, “Neurotrophin 3 is increased in the spontaneously hypertensive rat,” Journal of Hypertension, vol. 19, no. 12, pp. 2251–2256, 2001.
View at: Publisher Site | Google Scholar
D. Konukoglu and H. Uzun, “Endothelial dysfunction and hypertension,” Advances in Experimental Medicine and Biology, vol. 956, pp. 511–540, 2016.
View at: Publisher Site | Google Scholar
Y.-L. Tain, S. Leu, K. L. H. Wu, W.-C. Lee, and J. Y. H. Chan, “Melatonin prevents maternal fructose intake-induced programmed hypertension in the offspring: Roles of nitric oxide and arachidonic acid metabolites,” Journal of Pineal Research, vol. 57, no. 1, pp. 80–89, 2014.
View at: Publisher Site | Google Scholar
Y.-L. Tain, L.-T. Huang, J. Y. H. Chan, and C.-T. Lee, “Transcriptome analysis in rat kidneys: Importance of genes involved in programmed hypertension,” International Journal of Molecular Sciences, vol. 16, no. 3, pp. 4744–4758, 2015.
View at: Publisher Site | Google Scholar
E. Gerasimova, B. Audit, S. G. Roux et al., “Wavelet-based multifractal analysis of dynamic infrared thermograms to assist in early breast cancer diagnosis,” Frontiers in Physiology, vol. 5, p. 176, 2014.
View at: Google Scholar
T. J. Aitman, T. Gotoda, A. L. Evans et al., “Quantitative trait loci for cellular defects in glucose and fatty acid metabolism in hypertensive rats,” Nature Genetics, vol. 16, no. 2, pp. 197–201, 1997.
View at: Publisher Site | Google Scholar
J. Du, L. M. Fan, A. Mai, and J.-M. Li, “Crucial roles of Nox2-derived oxidative stress in deteriorating the function of insulin receptors and endothelium in dietary obesity of middle-aged mice,” British Journal of Pharmacology, vol. 170, no. 5, pp. 1064–1077, 2013.
View at: Publisher Site | Google Scholar
Y. Pan, Q. Qiao, L. Pan et al., “Losartan reduces insulin resistance by inhibiting oxidative stress and enhancing insulin signaling transduction,” Experimental and Clinical Endocrinology & Diabetes, vol. 123, no. 03, pp. 170–177, 2015.
View at: Publisher Site | Google Scholar
Z. Meng, T. Moroishi, and K.-L. Guan, “Mechanisms of Hippo pathway regulation,” Genes & Development, vol. 30, no. 1, pp. 1–17, 2016.
View at: Publisher Site | Google Scholar
D. Shao, P. Zhai, D. P. Del Re et al., “A functional interaction between Hippo-YAP signalling and FoxO1 mediates the oxidative stress response,” Nature Communications, vol. 5, article no. 3315, 2014.
View at: Publisher Site | Google Scholar
J. D. Malhotra and R. J. Kaufman, “Endoplasmic reticulum stress and oxidative stress: a vicious cycle or a double-edged sword?” Antioxidants & Redox Signaling, vol. 9, no. 12, pp. 2277–2293, 2007.
View at: Publisher Site | Google Scholar
A. Buday, P. Orsy, M. Godó et al., “Elevated systemic TGF-β impairs aortic vasomotor function through activation of NADPH oxidase-driven superoxide production and leads to hypertension, myocardial remodeling, and increased plaque formation in apoE -/- mice,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 299, no. 2, pp. H386–H395, 2010.
View at: Publisher Site | Google Scholar
J. L. Miguel-Carrasco, A. Baltanás, C. Cebrián et al., “Blockade of TGF-beta 1 signalling inhibits cardiac NADPH oxidase overactivity in hypertensive rats,” Oxidative Medicine and Cellular Longevity, vol. 2012, Article ID 726940, 8 pages, 2012.
View at: Publisher Site | Google Scholar
A. Alé, Y. Zhang, C. Han, and D. Cai, “Obesity-associated extracellular mtDNA activates central TGFβ pathway to cause blood pressure increase,” American Journal of Physiology-Endocrinology and Metabolism, vol. 312, no. 3, pp. E161–E174, 2017.
View at: Publisher Site | Google Scholar
E. L. Schiffrin, “Immune mechanisms in hypertension: How do T-regulatory lymphocytes fit in?” Journal of Hypertension, vol. 31, no. 10, pp. 1944-1945, 2013.
View at: Publisher Site | Google Scholar
M. Rosas-Ballina, P. S. Olofsson, M. Ochani et al., “Acetylcholine-synthesizing T cells relay neural signals in a vagus nerve circuit,” Science, vol. 334, no. 6052, pp. 98–101, 2011.
View at: Publisher Site | Google Scholar
P. S. Olofsson, B. E. Steinberg, R. Sobbi et al., “Blood pressure regulation by CD4+ lymphocytes expressing choline acetyltransferase,” Nature Biotechnology, vol. 34, no. 10, pp. 1066–1071, 2016.
View at: Publisher Site | Google Scholar
B. M. Matta, J. M. Lott, L. R. Mathews et al., “IL-33 is an unconventional alarmin that stimulates IL-2 secretion by dendritic cells to selectively expand IL-33R/ST2⁺ regulatory T cells,” The Journal of Immunology, vol. 193, no. 8, pp. 4010–4020, 2014.
View at: Publisher Site | Google Scholar
V. A. Braga, I. A. Medeiros, T. P. Ribeiro, M. S. França-Silva, M. S. Botelho-Ono, and D. D. Guimarães, “Angiotensin-II-induced reactive oxygen species along the SFO-PVN-RVLM pathway: implications in neurogenic hypertension,” Brazilian Journal of Medical and Biological Research, vol. 44, no. 9, pp. 871–876, 2011.
View at: Publisher Site | Google Scholar
D. P. Pollow, J. Uhrlaub, M. J. Romero-Aleshire et al., “Sex differences in T-lymphocyte tissue infiltration and development of Angiotensin II Hypertension,” Hypertension, vol. 64, no. 2, pp. 384–390, 2014.
View at: Publisher Site | Google Scholar
F. C. Nunes and V. A. Braga, “Chronic angiotensin II infusion modulates angiotensin II type i receptor expression in the subfornical organ and the rostral ventrolateral medulla in hypertensive rats,” Journal of the Renin-Angiotensin-Aldosterone System, vol. 12, no. 4, pp. 440–445, 2011.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Si-Ming Ma 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1799

Downloads

987

Citations