Skip to main content

Use of serial analysis of gene expression to reveal the specific regulation of gene expression profile in asthmatic rats treated by acupuncture

Abstract

Background

Asthma has become an important public health issue and approximately 300 million people have suffered from the disease worldwide. Nowadays, the use of acupuncture in asthma is increasing. This study intended to systematically analyze and compare the gene expression profiles between the asthmatic and acupuncture-treated asthmatic rat lung, and tried to gain insight into the molecular mechanism underlying the early airway response (EAR) phase of asthma treated by acupuncture.

Methods

Four tag libraries of serial analysis of gene expression (SAGE) were established from lung tissues of control rats (CK), asthmatic rats (AS), asthmatic rats treated by acupuncture (ASAC), and control rats treated by acupuncture (CKAC). Bioinformatic analyses were carried out by using the methods including unsupervised hierarchical clustering, functional annotation tool of the database for annotation, visualization, and integrated discovery (DAVID), gene ontology (GO) tree machine, and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis.

Results

There were totally 186 differentially expressed tags (P < 0.05, PCK/AS) between the libraries of CK and AS, 130 differentially expressed tags between libraries of AS/ASAC (P < 0.05, PAS/ASAC), and 144 differentially expressed tags between libraries of CK/CKAC (P < 0.05, PCK/CKAC). The gene expression profiles of AS and ASAC were more similar than other libraries via unsupervised SAGE clustering. By comparison of PCK/AS and PAS/ASAC, the DAVID genes functional classification was found to be changed from "immune response" to "response to steroid hormone stimulus", and the GO term "antigen processing and presentation of peptide antigen" disappeared in PAS/ASAC. Totally 3 same KEGG pathways were found among the three groups. Moreover, 21 specific tags of the acupuncture in treating asthma were detected using Venn diagrams.

Conclusion

Our SAGE research indicates that the gene expression profile of the EAR phase of asthma could be effectively and specifically regulated by acupuncture, which suggests that the gene expression of immune response and steroid hormone may play an important role in the treatment.

Background

Asthma is a complex syndrome involving potentially permanent airway obstruction, airway hyperresponsiveness, and multicellular inflammation. It is estimated that approximately 300 million people have suffered from asthma worldwide and the burden of this disease in countries as well as families is increasing [1]. Although inhaled steroids can significantly improve the symptoms, curative therapies are not yet available [2]. Moreover, there are also significant concerns regarding the potential side effects from the long term use of conventional drugs, such as corticosteroids. Thus, an effective, low-risk, and non-drug strategy would provide a valuable and adjunctive treatment in asthma management [3]. Complementary and alternative medicine (CAM), such as traditional Chinese herbal remedies, homeopathy as well as acupuncture, is widely applied in the asthma management. It is reported that CAM has been used among 59% of patients with asthma or rhinosinusitis in the United Kingdom, 41% in the United States, 26.5% in Germany, and 27.2% in Singapore [4].

Acupuncture literally means to puncture with a needle, which is an important therapy in traditional Chinese medicine (TCM) for at least 2,500 years [5]. The theory of acupuncture holds that there are different types of energy flow (qi) in the human body, and the disruptions of qi are believed to be responsible for diseases. Acupuncture practitioners may use thin, solid, metallic needles to correct the imbalances of the flow via the stimulation of special points in the body, which are manipulated manually or by electrical stimulation. The therapy is developed as a relatively global system of medicine and is utilized to treat many diseases. The World Health Organization listed asthma and other 42 indications for acupuncture in 1979 [6] and classified the diseases treated by acupuncture into four categories, 107 illnesses in 2002 [7]. The National Institutes of Health (NIH) has accepted the validity of acupuncture treatment [8] and recommended it as an adjunctive treatment in comprehensive management programs of addiction, stroke rehabilitation, and asthma, etc [5, 9]. Acupuncture has traditionally been used to treat asthma in China and has been shown to be beneficial in acute asthma in short term [10, 11]. However, well-designed scientific researches in this field are needed and encouraged, which are not only important for elucidating the mechanism, but also useful for exploring new pathways in a systematic manner.

Acupuncture is a complex intervention on diseases and many studies have demonstrated that acupuncture can cause multiple biological responses and regulate many cellular and physiological processes, which could lead to changes of gene expression [5]. These processes could occur at either the proximal or distal ends of acupuncture application, which may be mediated mainly by neural and humoral mechanisms [5, 12]. High-throughput technologies, such as microarray and serial analysis of gene expression (SAGE), may help to reveal and clarify the possible mechanism of acupuncture. In the nerve system, acupuncture was reported to regulate brain aging related genes [13], suppress several genes in the nociceptive pathways [14], and up-regulate genes in the spinal cord injury [15] by using microarray. In the immune and endocrine system, results of microarray have revealed that acupuncture could regulate genes to increase activities of natural killer cell [16] and superoxide dismutase [17], to keep the cytokines balance between Th1 and Th2 [18], and to induce hypocholesterolemic effects [19]. However, these researches focused mainly on the nerve and immune systems, and the analyses of the generated data were limited. SAGE, as a powerful expression profiling method and much more accessible for dissecting the complex system, was applied to qualitatively and quantitatively evaluate the transcription of the genes via particular length without the prerequisite of a hybridization probe for each transcript [20]. However, there has been no SAGE report concerning the mechanism of acupuncture in treating asthma so far.

In this study, four SAGE tag libraries were established respectively from lung tissues of control rats, asthmatic rats, asthmatic rats treated by acupuncture, and control rats treated by acupuncture. The study aimed to systematically analyze and compare the gene expression profiles in the lung of the four different groups, and attempted to examine the molecular mechanism of early airway response (EAR) phase of asthma treated by acupuncture.

Methods

Animal and experimental asthma model

Pathogen-free, male Sprague-Dawley (SD) rats (4 weeks old, 110–130 g, SLAC Laboratory Animal Co. Ltd., Shanghai, China), raised in a pathogen-free rodent facility and provided with food and water ad libitum, were randomly divided into four groups (each group contained 8 rats): control rats (CK), asthmatic rats (AS), asthmatic rats treated by acupuncture (ASAC), control rats treated by acupuncture (CKAC). The protocol of SD rat model of asthma was described as previously [21] and rats of CK and CKAC were sensitized and boosted to normal saline instead of OVA. Rats were kept in animal facilities approved by the Shanghai Committee for Accreditation of Laboratory Animal and the animal experiment conformed to the regulations of the State Science and Technology Commission.

Acupuncture treatment

Manual acupuncture was performed for two weeks from the first day after the sensitization, once every other day (Figure 1). The acupuncture points, Dazhui (GV 14, between C7 and T1 vertebrae), bilateral Fengmen (BL12, foveola laterally between T2 and T3 vertebrae) and bilateral Feishu (BL13, foveola laterally between T3 and T4 vertebrae), were selected based on the theory of TCM in treating asthma [22]. In consideration of the specific and non-specific effects of acupuncture [23], the same points in normal rats were selected and acted as control of ASAC, which could distinguish the specificity of acupuncture in treating asthma. The stainless needles (0.30 × 13 mm, Suzhou Medical Appliance Factory, Suzhou, China) were inserted about 5 mm deep into the skin. The needles were twisted approximately 360° evenly at the rate of 60 times/min for 20 s, manipulated every 5 min and withdrawn after 20 min. For the convenient manipulation of the acupuncture points on the back, rat was placed on the suspended shelf (50 × 45 mm, about 50 cm high from the ground, see additional file 1), which could easily make it calm and stand still without anesthesia. Acupuncture was performed by the same experienced practitioner and the animals were handled while awake, with special care to minimize stress. Both rats in ASAC and CKAC received the same acupuncture treatment.

Figure 1
figure 1

Protocol and measurements for the rat model of allergic asthma. Male Sprague-Dawley rats were sensitized intraperitoneally at day 0 with 1 mg of ovalbumin (OVA) precipitated with 10 mg of aluminium hydroxide gel in 1 ml normal saline. Manual acupuncture was performed regularly from the first day after the sensitization and given every other day over a period of two weeks. On day 14, the sensitized animals were challenged with 1 ml/kg of 5% OVA in normal saline (5 mg/kg) by injection into the external jugular vein over 10 s. The respiratory parameters of pulmonary resistance, dynamic compliance and respiratory rate were immediately recorded for 10 min. Rats were killed immediately after the measurement. Rats of CK and CKAC were sensitized and challenged with normal saline instead of OVA. Rats of CKAC received the same acupuncture treatment as those of ASAC.

Measurement of pulmonary function

The measurements of pulmonary resistance (RL), dynamic compliance (Cdyn), and respiratory rate (RR), were modified from Glaab T et al [24, 25]. Briefly, a rat was placed in supine position on a wood plate warmed by an incandescent lamp after anaesthesia. At the upper part of the trachea, a T-shape cutting was made and a T-shape cannula, which was directly attached to a heater controlled pneumotachograph (Series 3850A, Hans Rudolph, USA), was gently inserted into the trachea. Tidal flow was determined by the pneumotachograph connected to a differential pressure transducer (600D-011, AutoTran, USA). To measure transpulmonary pressure, a water-filled PE-90 tubing was inserted into the esophagus to the level of the midthorax (lower one-third of the esophagus) and coupled to a pressure transducer (PT14MX, Jialong Teaching Equipment, Shanghai). The pneumotachograph tidal flow signal was integrated with time to obtain tidal volume. RL and Cdyn were calculated over a complete respiratory cycle using an integration method over flows, volumes and pressures, and were continuously recorded with software (Shanghai Medical College, Fudan University) for physiology experiments. Respiratory parameters were averaged in 60 s segments and maximum RL, minimum Cdyn and change of RR values were taken and calculated as differential values subtracted from the corresponding baseline values (Figure 1). Then the rats were sacrificed. The lungs were excised right away, rinsed in ice-cold normal saline, dissected free from surrounding tissues in an ice-bath and frozen immediately in liquid nitrogen.

Construction, annotation, and confirmation of the SAGE libraries

Construction and annotation of the SAGE libraries were described as previously [21]. The confirmation of the four SAGE libraries was performed by Quantitative Real-Time PCR (qRT-PCR) on an Applied Biosystems 7300 Real-Time PCR System using TOYOBO Realtime PCR Master Mix (Toyobo, Osaka, Japan). The threshold cycle number was determined using SDS v1.4 Software and the reactions were performed in triplicate. Total RNA (5 μg) was reversely transcribed into cDNA by using the RevertAid First Strand cDNA synthesis kit (Cat. No. K1622; Fermentas, EU). For qRT-PCR of the cDNA, primer pairs were designed to generate intron-spanning products of 101–150 bp (Primer sequences were listed in additional file 2). The generation of specific PCR products was confirmed by the melting curve and gel analysis. The expressional ratio was calculated according to the formula 2(Rt-Et)/2(Rn-En) as described previously [26]. Transcripts with a twofold increase in expression were considered to be up-regulated and those with a 0.5-fold decrease in expression were considered to be down-regulated.

Bioinformatic analysis of SAGE tags

To identify genes preferentially regulated in the four groups, the two-way unsupervised clustering method based on tag copies was applied. This unguided approach allowed pattern discovery for subsequent supervised functional analysis. To reduce the magnitude effects of the extreme data, the genes with total tag counts less than 20 were filtered. The dendrogram of differentially expressed tags was created by using the TIGR MultiExperiment Viewer 4.0 http://www.tm4.org/mev.html[27], mainly with the average clustering and Euclidean distance.

Differentially expressed genes (P < 0.05) between SAGE libraries were functionally annotated and classified by using the functional annotation tool of database for annotation, visualization, and integrated discovery (DAVID) http://david.abcc.ncifcrf.gov/[28], which provided integrated solutions for the annotation and analysis of genome-scale datasets derived from high-throughput technologies.

Key regulatory processes in asthma were analyzed by Gene Ontology (GO) Tree Machine http://bioinfo.vanderbilt.edu/gotm/[29]. GO Tree Machine generated a directed acyclic graph (DAG) for input gene sets, which was made to identify the most important GO categories and to suggest their potential biological importance.

The Kyoto encyclopedia of genes and genomes (KEGG) pathway is a collection of manually drawn pathway maps of the molecular interaction and reaction networks. The KEGG pathways of the differentially expressed genes between SAGE libraries were matched by using the DAVID Functional Annotation Tool.

Statistical analysis

One-way ANOVA (analysis of variance) followed by the least significant difference (LSD) test for post hoc analysis was used to analyze the significance of RL, Cdyn and RR among the four groups. Statistical analysis for the significance of each of the four SAGE libraries was made using Monte Carlo analysis. The enrichments of GO Tree Machine were statistically significant as determined by the hypergeometric test [29].

Results

Measurements of pulmonary functions after OVA challenge

The allergen-specific early airway response to ovalbumin (OVA) in sensitized rats showed significant increases in RL and significant decreases in simultaneously measuring Cdyn and RR compared with those of the controls, thus indicating an allergen-specific EAR to OVA. By One-way ANOVA among the four groups, the significant differences (P < 0.05) were at 1–7 min for RL, 2 and 5 min for Cdyn, 2–4 min for RR after challenge. The LSD test for post hoc analysis indicated that RL in ASAC group was significantly decreased at 1–4 min after challenge in comparison with AS group (P < 0.05, Table 1), and Cdyn and RR in ASAC group were significantly increased at 2, 5 min and 2, 3 min respectively after challenge when compared with those of AS group (P < 0.05, see additional file 3). The results indicated that the immediate effects (1–3 min) of acupuncture were the most significant in the EAR phase of asthma.

Table 1 The comparisons of pulmonary resistance (kPa/ml/s) of the four groups

General analysis of SAGE libraries

The four SAGE libraries of rat lung were deposited in the SAGEmap database at National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/geo, Accession numbers are GSM45195, GSM119459, GSM279944, and GSM279945). The information of the matched genes and ESTs of the 4 libraries was listed in the Table 2.

Table 2 Summary of serial analysis of gene expression data for the four libraries

By Comparing with the SAGE data between the libraries of CK and AS, there were totally 186 differentially expressed tags (P < 0.05, PCK/AS). Similarly, there were 130 differentially expressed tags between libraries of AS and ASAC (P < 0.05, PAS/ASAC), and 144 differentially expressed tags between libraries of CK and CKAC (P < 0.05, PCK/CKAC). See additional file 4 for the lists of differentially expressed tags of PCK/AS, PAS/ASAC, PCK/CKAC.

Confirmation of SAGE Results by qRT-PCR

To confirm the expression profiles among the four SAGE libraries, three differentially expressed genes of interest were chosen and their expression levels were evaluated by qRT-PCR. Dusp1 encodes a protein that catalyzes the dephosphorylation and inactivation of MAP kinase. S100A9 is a calcium binding protein that may be associated with acute inflammatory processes. MT-2 has a high content of cysteine residues that bind various heavy metals and could be transcriptionally regulated by both heavy metals and glucocorticoids. The results indicate that the expression profiles of representative genes in qRT-PCR analysis correspond to the SAGE data (Figure 2), which supports the validation of our quantitative data for further bioinformatic analysis.

Figure 2
figure 2

Real-time PCR confirmation of differentially regulated genes of interest predicted by SAGE. * meant the expression levels were up-regulated more than 2 folds compared with CK, # meant the expression levels were down-regulated more than 0.5 fold compared with AS. The expression levels of the three genes of CKAS were up-regulated less than 2 folds compared with CK.

Unsupervised SAGE clustering based on tag expression

Based on the global expression matrix, the dendrogram clearly identified CK from the other 3 libraries, while AS and ASAC displayed more similar patterns than CKAC on the gene expression profiles (Figure 3). The genes with over fivefold change or less than 0.2-fold between AS and ASAC were listed at the right side of the hierarchical clustering dendrogram. Seven genes were ranked with an expression over 5-fold increase: Fxyd3 (Rn.3896), Col6a2 (Rn.11889), Uba52 (Rn.4300), Tm9sf2 (Rn.11839), CD24 (Rn.6007), Mgp (Rn.2379), and ribosomal protein L31 (Rpl31, Rn.1101). Four genes were listed with an expression over 5-fold decrease: Ap1m1 (Rn.139185), Mgst2 (Rn.7854), Atp1b1 (Rn.8925), and Ywhaq (Rn.2502).

Figure 3
figure 3

The dendrogram of the hierarchical cluster analysis. The genes with over fivefold change or less than 0.2-fold between AS and ASAC were listed at the right side of the graph. Abbreviations of the seven up-regulated genes: Fxyd3: FXYD domain-containing ion transport regulator 3; Col6a2: procollagen, type VI, alpha 2; Uba52: ubiquitin A-52 residue ribosomal protein fusion product 1; Tm9sf2: transmembrane 9 superfamily member 2; CD24: CD24 antigen; Mgp: matrix Gla protein; Rpl31: ribosomal protein L31. Abbreviations of the four down-regulated genes: Ap1m1: Adaptor-related protein complex AP-1, mu subunit; Mgst2: Microsomal glutathione S-transferase 2; Atp1b1: ATPase, Na+/K+ transporting, beta 1 polypeptide; Ywhaq: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, theta polypeptide.

The comparisons of the DAVID gene functional classification of PCK/AS, PAS/ASAC and PCK/CKAC

Six DAVID gene functional classifications of PAS/ASAC with the enrichment score equal to or higher than 1.0 were found. They were cellular biosynthetic process, homeostatic process, response to steroid hormone stimulus, cell migration, cellular process, and cellular lipid metabolic process. The related gene functional classifications among PCK/AS, PAS/ASAC and PCK/CKAC were compared (Figure 4). The classification of homeostasis process was only found in PAS/ASAC. The DAVID gene functional classification was found to be changed from "immune response" in PCK/AS to "response to steroid hormone stimulus" in PAS/ASAC. (See additional file 5 for the gene lists of DAVID gene functional classification of PCK/AS, PAS/ASAC, PCK/CKAC)

Figure 4
figure 4

The comparison of the selected DAVID gene functional classification among P CK/AS , P AS/ASAC and P CK/CKAC . Six DAVID gene functional classifications with enrichment score equal to or higher than 1 were selected. There is no classification homeostatic process in PCK/AS and PCK/CKAC, so the enrichment score is zero.

DAG analysis of gene categories of PCK/AS, PAS/ASAC and PCK/CKAC

There were 25, 5 and 10 enriched GO categories of biological process in PCK/AS, PAS/ASAC and PCK/CKAC respectively (Figure 5, 6, 7). Biosynthesis was only one overlapping GO term among them. The GO terms, such as antigen processing and presentation, cell migration, acetylcholine receptor signaling, muscarinic pathway disappeared in PAS/ASAC, and new GO terms, like regulation of biosynthesis, regulation of liquid surface tension, epidermal growth factor receptor signaling pathway appeared in PAS/ASAC when compared with PCK/AS. Regulation of liquid surface tension was the only shared GO term between PAS/ASAC and PCK/CKAC.

Figure 5
figure 5

The DAG graphs of the P CK/AS . After uploading 186 differentially expressed tags between the libraries of CK and AS to the GO Tree Machine, the directed acyclic graph was generated automatically. Totally 25 enriched GO categories of biological process were found in PCK/AS.

Figure 6
figure 6

The DAG graphs of the P AS/ASAC . 130 differentially expressed tags between libraries of AS and ASAC were uploaded to the GO Tree Machine and totally 5 enriched GO categories of biological process were found in PAS/ASAC.

Figure 7
figure 7

The DAG graphs of the P CK/CKAC . 144 differentially expressed tags between libraries of CK and CKAC have matched 10 enriched GO categories of biological process by using the GO Tree Machine.

Finding KEGG pathways of PCK/AS, PAS/ASAC and PCK/CKAC

For understanding functional roles of the differentially expressed genes, KEGG pathway analysis was assigned by applying the DAVID annotation tool. There were 15, 8, and 2 matched KEGG pathways (counts ≥ 3) in PCK/AS, PAS/ASAC and PCK/CKAC respectively (see Table 3), in which tight junction pathway kept the same. Three pathways, MAPK signaling pathway, T cell receptor signaling pathway, focal adhesion, were identical between PCK/AS and PAS/ASAC. Although there are the same matching pathways among three groups, the involved genes were different. (See additional file 6)

Table 3 Kyoto encyclopedia of genes and genomes (KEGG) pathways of differentially expressed tags

The Venn diagram among PCK/AS, PAS/ASAC, and PCK/CKAC

By using the Venn diagram, 21 tags were attributed to the same group between PCK/AS and PAS/ASAC, 7 tags were found involved in the three groups, 6 tags were identified between PAS/ASAC and PCK/CKAC, and 25 tags were commonly clustered between PCK/AS and PCK/CKAC(Figure 8). (See additional file 7 for the gene lists of each group)

Figure 8
figure 8

The Venn diagram among P CK/AS , P AS/ASAC , and P CK/CKAC . The diagram shows 21 common tags between PCK/AS and PAS/ASAC, 7 common tags in the three groups, 25 common tags between PCK/AS and PCK/CKAC, and 6 common tags between PAS/ASAC and PCK/CKAC.

Discussion

The gene expression profiles of acupuncture in treating asthma were not yet well studied before. Our SAGE study, in which acupuncture served as a kind of biological perturbation of systems biology, has provided a molecular base of EAR phase of asthma treated by acupuncture and may provide clues for the further research.

The cluster analysis of gene expression profiles of the four libraries and genes with marked changes

Based on the global expression matrix, the result of the hierarchical cluster analysis indicated that the gene expression profiles of different libraries changed a lot when receiving different perturbation, such as OVA sensitization and acupuncture. It suggested that acupuncture could perturb the biological condition by regulating a number of genes. Moreover, the gene expression profiles of ASAC and CKAC were different, and the profiles of AS and ASAC were more similar than those in other groups, which indicated the specific effect of acupuncture in the EAR phase of asthma.

The genes with marked changes in expression levels probably have potential importance in treating asthma. CD24, referred to as heat-stable antigen, is an important co-stimulatory molecule in immunity [30]. Previous studies have suggested that CD24 was involved in cell adhesion and signal transduction via phosphorylation of intracellular proteins, intracellular calcium mobilization, and activation of transcription factors [31]. In our libraries, the expression of CD24 was found to be down-regulated 12 times in asthma while up-regulated 7 times by acupuncture. The result suggested that the gene may correlate with the immune-modulating effects of acupuncture in EAR phase of asthma.

The gene expressions of immune response and steroid hormone were regulated by acupuncture in the EAR phase of asthma

Several immune-related GO terms, such as antigen processing and presentation, antigen processing and presentation of peptide antigen, antigen processing and presentation of exogenous peptide antigen, disappeared in PAS/ASAC via GO Tree Machine. Moreover, the DAVID genes functional classification was found to shift from "immune response" in PCK/AS to "response to steroid hormone stimulus" in PAS/ASAC. The changes suggested that the immune response of EAR phase of asthma was attenuated by acupuncture, which may be realized through the release of the endogenous steroid hormone.

There were 14 genes in the DAVID genes functional classification "response to steroid hormone stimulus" in PAS/ASAC, in which the genes of MGP (Rn.2379) and Srd5a2 (Rn.9938) were interesting. MGP is a vitamin K-dependent protein that serves as a substrate for the enzyme γ-carboxylase [32]. It plays a role in lung growth and development [33] and its expression was incited by dexamethasone [34]. In our libraries, the expression of MGP was up-regulated by acupuncture, which may result from the release of endogenous steroid hormones in consideration of the DAVID gene classification. Srd5a2, an isozyme of the 5-alpha-reductase family, is present in a large number of cells. It plays a key role in the conversion of testosterone to dihydrotestosterone and in the removal of excess of potentially neurotoxic steroids [35]. In our libraries, the expression of Srd5a2 was up-regulated in asthma but down-regulated by acupuncture. The changed expression suggested that Srd5a2 served as an adjustable gene target of acupuncture, which may relate with the genesis of steroid by acupuncture in the EAR phase of asthma.

Pathways regulated by acupuncture in the EAR phase of asthma

Pathway-level visualization of omics data provides an essential means for systems biology to capture the systematic properties of the inner activities of cells. Eight pathways were involved in the PAS/ASAC, in which tight junction pathway was commonly shared among PCK/AS, PAS/ASAC, and PCK/CKAC. Three pathways were the same between PCK/AS and PAS/ASAC, and 4 pathways were unique in PAS/ASAC.

The tight junction pathway was considered to be the background of OVA sensitization and acupuncture regulation. MAPK cascade is one of the three same pathways between PCK/AS and PAS/ASAC. It contributes to the amplification and specificity of the transmitted signals and plays discrete yet complementary roles in accentuating allergic airway inflammation [36]. It is also reported that the MAPK pathway could regulate steroidogenesis [37]. GnRH signaling pathway is one of the four unique pathways in PAS/ASAC, which is the central regulator of the reproductive hormonal cascade [38] and could enhance the basal steroidogenesis [39]. Several signaling pathways are activated by GnRH, including MAPK, and protein kinase C [40], which could lead to the increase of corticotropin-releasing hormone-binding protein at mRNA level [41]. In our study, the data suggested that GnRH signaling pathway may interact with other pathways and participate in the genesis of steroid by acupuncture in the EAR phase of asthma.

The specific and non-specific tags of acupuncture

From the Venn diagrams, the 21 tags were found to represent the specific genes of acupuncture in treating asthma. By using the DAVID classification analysis tool, 18 of the 21 tags were divided into 7 groups. The classification with highest score was "response to steroid hormone stimulus" (enrichment score = 1.71), which was in accordance with the classification of the PAS/ASAC. The genes involved in the classification were: MGP, Abca2, and Sult1a1, which demonstrated that the transportation and production of steroid hormone may be regulated by acupuncture in the EAR phase of asthma at the level of transcription.

Abca2 expressed in a broad range of tissues and previous studies indicated that the gene participated in the transport of steroids [42]. Besides, the elevated expression of Abca2 was considered to be a conserved mechanism of cell survival [43]. In our study, Abca2 was up-regulated in asthma, which may correlate with the cell survival in the EAR phase of asthma. Sult1a1 is a member of sulfotransferase families, primarily localized in the trans-Golgi apparatus and associated with the sulfation of steroids and various hormones [44, 45]. Our SAGE data demonstrated that Sult1a1 was up-regulated in asthma, which indicated the gene may catalyze the sulfation process of steroidgensis. However, the two genes both belong to the DAVID classification "homeostatic process" and the gene expressions of Abca2 and Sult1a1 after acupuncture were almost equivalent to those of the controls. It suggested that acupuncture could alleviate the OVA challenge by maintaining the internal equilibrium.

Seven common tags were found among the different groups when receiving the chemical and acupuncture stimulation both in normal and asthmatic conditions. Twenty-five tags were found to be the same when receiving the OVA sensitization or acupuncture treatment. Furthermore, 6 tags were identified as the non-specific genes of acupuncture in PAS/ASAC and PCK/CKAC. The above-mentioned genes belong mainly to cellular process, and cell communication, which suggested that they may serve as the background genes in EAR phase of asthma and acupuncture. For example, Syntaxin 5 is a Golgi-localized SNARE protein required for endoplasmic reticulum-Golgi traffic in yeast and Golgi reassembly following cell division in mammalian cells [46].

Conclusion

We have presented evidences of gene expression profiles in the lung of asthmatic rats and those treated by acupuncture via SAGE. The study indicates that the gene expression profile of the EAR phase of asthma could be effectively and specifically regulated by acupuncture at the transcriptional level, which suggests that the gene expression of immune response and steroid hormone may play an important role in the treatment.

References

  1. Masoli M, Fabian D, Holt S, Beasley R: The global burden of asthma: executive summary of the GINA Dissemination Committee report. Allergy. 2004, 59: 469-478. 10.1111/j.1398-9995.2004.00526.x.

    Article  PubMed  Google Scholar 

  2. Chu EK, Drazen JM: Asthma: one hundred years of treatment and onward. Am J Respir Crit Care Med. 2005, 171: 1202-1208. 10.1164/rccm.200502-257OE.

    Article  PubMed  Google Scholar 

  3. McCarney RW, Lasserson TJ, Linde K, Brinkhaus B: An overview of two Cochrane systematic reviews of complementary treatments for chronic asthma: acupuncture and homeopathy. Respir Med. 2004, 98: 687-696. 10.1016/j.rmed.2004.05.005.

    Article  CAS  PubMed  Google Scholar 

  4. Györik SA, Brutsche MH: Complementary and alternative medicine for bronchial asthma: is there new evidence?. Curr Opin Pulm Med. 2004, 10: 37-43. 10.1097/00063198-200401000-00007.

    Article  PubMed  Google Scholar 

  5. NIN Consensus Development Panel on Acupuncture: Acupuncture. JAMA. 1998, 280: 1518-1524. 10.1001/jama.280.17.1518.

    Article  Google Scholar 

  6. Use of acupuncture in modern health care. WHO Chron. 1980, 34: 294-301.

  7. World Health Organization: Acupuncture:review and analysis of reports on controlled clinical trials. 2002

    Google Scholar 

  8. Culliton BJ: NIH says "yes" to acupuncture. Nat Med. 1997, 3: 1307-

    CAS  PubMed  Google Scholar 

  9. Acupuncture: NIH Consens Statement 1997 Nov 35. 1997, 15: 1-34.

    Google Scholar 

  10. Chu KA, Wu YC, Lin MH, Wang HC: Acupuncture resulting in immediate bronchodilating response in asthma patients. J Chin Med Assoc. 2005, 68: 591-594.

    Article  PubMed  Google Scholar 

  11. Takishima T, Mue S, Tamura G, Ishihara T, Watanabe K: The bronchodilating effect of acupuncture in patients with acute asthma. Ann Allergy. 1982, 48: 44-49.

    CAS  PubMed  Google Scholar 

  12. Han JS: Physiology of acupuncture: review of thirty years of research. J Altern Complement Med. 1997, 3: S101-108. 10.1089/acm.1997.3.101.

    Google Scholar 

  13. Ding X, Yu J, Yu T, Fu Y, Han J: Acupuncture regulates the aging-related changes in gene profile expression of the hippocampus in senescence-accelerated mouse (SAMP10). Neurosci Lett. 2006, 399: 11-16. 10.1016/j.neulet.2006.01.067.

    Article  CAS  PubMed  Google Scholar 

  14. Ko J, Na DS, Lee YH, Shin SY, Kim JH, Hwang BG, Min BI, Park DS: cDNA microarray analysis of the differential gene expression in the neuropathic pain and electroacupuncture treatment models. J Biochem Mol Biol. 2002, 35: 420-427.

    Article  CAS  PubMed  Google Scholar 

  15. Wang XY, Li XL, Hong SQ, Xi-Yang YB, Wang TH: Electroacupuncture Induced Spinal Plasticity is Linked to Multiple Gene Expressions in Dorsal Root Deafferented Rats. J Mol Neurosci. 2009, 37: 97-110. 10.1007/s12031-008-9095-1.

    Article  CAS  PubMed  Google Scholar 

  16. Kim CK, Choi GS, Oh SD, Han JB, Kim SK, Ahn HJ, Bae H, Min BI: Electroacupuncture up-regulates natural killer cell activity Identification of genes altering their expressions in electroacupuncture induced up-regulation of natural killer cell activity. J Neuroimmunol. 2005, 168: 144-153. 10.1016/j.jneuroim.2005.07.005.

    Article  CAS  PubMed  Google Scholar 

  17. Rho SW, Choi GS, Ko EJ, Kim SK, Lee YS, Lee HJ, Hong MC, Shin MK, Min BI, Kee HJ: Molecular changes in remote tissues induced by electro-acupuncture stimulation at acupoint ST36. Mol Cells. 2008, 25: 178-183.

    CAS  PubMed  Google Scholar 

  18. Shiue HS, Lee YS, Tsai CN, Hsueh YM, Sheu JR, Chang HH: DNA Microarray Analysis of the Effect on Inflammation in Patients Treated with Acupuncture for Allergic Rhinitis. J Altern Complement Med. 2008, 14: 689-98. 10.1089/acm.2007.0669.

    Article  PubMed  Google Scholar 

  19. Li M, Zhang Y: Modulation of gene expression in cholesterol-lowering effect of electroacupuncture at Fenglong acupoint (ST40): a cDNA microarray study. Int J Mol Med. 2007, 19: 617-629.

    CAS  PubMed  Google Scholar 

  20. Velculescu VE, Zhang L, Zhou W, Vogelstein J, Basrai MA, Bassett DE, Hieter P, Vogelstein B, Kinzler KW: Characterization of the yeast transcriptome. Cell. 1997, 88: 243-251. 10.1016/S0092-8674(00)81845-0.

    Article  CAS  PubMed  Google Scholar 

  21. Yin LM, Jiang GH, Wang Y, Wang Y, Liu YY, Jin WR, Zhang Z, Xu YD, Yang YQ: Serial analysis of gene expression in a rat lung model of asthma. Respirology. 2008, 13: 972-982.

    PubMed  Google Scholar 

  22. Shao JM, Ding YD: Clinical observation on 111 cases of asthma treated by acupuncture and moxibustion. J Tradit Chin Med. 1985, 5: 23-25.

    CAS  PubMed  Google Scholar 

  23. Endres HG: Acupuncture: specific and non-specific effects. Forsch Komplement Med. 2008, 15: 6-8. 10.1159/000113416.

    Article  Google Scholar 

  24. Glaab T, Ziegert M, Baelder R, Korolewitz R, Braun A, Hohlfeld JM, Mitzner W, Krug N, Hoymann HG: Invasive versus noninvasive measurement of allergic and cholinergic airway responsiveness in mice. Respir Res. 2005, 6: 139-10.1186/1465-9921-6-139.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Glaab T, Mitzner W, Braun A, Ernst H, Korolewitz R, Hohlfeld JM, Krug N, Hoymann HG: Repetitive measurements of pulmonary mechanics to inhaled cholinergic challenge in spontaneously breathing mice. J Appl Physiol. 2004, 97: 1104-1111. 10.1152/japplphysiol.01182.2003.

    Article  PubMed  Google Scholar 

  26. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  27. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.

    CAS  PubMed  Google Scholar 

  28. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.

    Article  PubMed  Google Scholar 

  29. Zhang B, Schmoyer D, Kirov S, Snoddy J: GOTree Machine (GOTM): a web-based platform for interpreting sets of interesting genes using Gene Ontology hierarchies. BMC Bioinformatics. 2004, 5: 16-23. 10.1186/1471-2105-5-16.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Li O, Zheng P, Liu Y: CD24 expression on T cells is required for optimal T cell proliferation in lymphopenic host. J Exp Med. 2004, 200: 1083-1089. 10.1084/jem.20040779.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Magnaldo T, Barrandon Y: CD24 (heat stable antigen, nectadrin), a novel keratinocyte differentiation marker, is preferentially expressed in areas of the hair follicle containing the colony-forming cells. J Cell Sci. 1996, 109 (Pt 13): 3035-3045.

    CAS  PubMed  Google Scholar 

  32. Kirfel J, Kelter M, Cancela LM, Price PA, Schule R: Identification of a novel negative retinoic acid responsive element in the promoter of the human matrix Gla protein gene. Proc Natl Acad Sci USA. 1997, 94: 2227-2232. 10.1073/pnas.94.6.2227.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Gilbert KA, Rannels SR: Matrix GLA protein modulates branching morphogenesis in fetal rat lung. Am J Physiol Lung Cell Mol Physiol. 2004, 286: L1179-1187. 10.1152/ajplung.00188.2003.

    Article  CAS  PubMed  Google Scholar 

  34. Gilbert KA, Rannels SR: Glucocorticoid effects on vitamin K-dependent carboxylase activity and matrix Gla protein expression in rat lung. Am J Physiol Lung Cell Mol Physiol. 2003, 285: L569-577.

    Article  CAS  PubMed  Google Scholar 

  35. Compagnone NA, Mellon SH: Neurosteroids: biosynthesis and function of these novel neuromodulators. Front Neuroendocrinol. 2000, 21: 1-56. 10.1006/frne.1999.0188.

    Article  CAS  PubMed  Google Scholar 

  36. Duan W, Wong WS: Targeting mitogen-activated protein kinases for asthma. Curr Drug Targets. 2006, 7: 691-698. 10.2174/138945006777435353.

    Article  CAS  PubMed  Google Scholar 

  37. Miyoshi T, Otsuka F, Inagaki K, Otani H, Takeda M, Suzuki J, Goto J, Ogura T, Makino H: Differential regulation of steroidogenesis by bone morphogenetic proteins in granulosa cells: involvement of extracellularly regulated kinase signaling and oocyte actions in follicle-stimulating hormone-induced estrogen production. Endocrinology. 2007, 148: 337-345. 10.1210/en.2006-0966.

    Article  CAS  PubMed  Google Scholar 

  38. Millar RP, Lu ZL, Pawson AJ, Flanagan CA, Morgan K, Maudsley SR: Gonadotropin-releasing hormone receptors. Endocr Rev. 2004, 25: 235-275. 10.1210/er.2003-0002.

    Article  CAS  PubMed  Google Scholar 

  39. Leung PC, Steele GL: Intracellular signaling in the gonads. Endocr Rev. 1992, 13: 476-498.

    CAS  PubMed  Google Scholar 

  40. Ferris HA, Shupnik MA: Mechanisms for pulsatile regulation of the gonadotropin subunit genes by GNRH1. Biol Reprod. 2006, 74: 993-998. 10.1095/biolreprod.105.049049.

    Article  CAS  PubMed  Google Scholar 

  41. Westphal NJ, Seasholtz AF: Gonadotropin-releasing hormone (GnRH) positively regulates corticotropin-releasing hormone-binding protein expression via multiple intracellular signaling pathways and a multipartite GnRH response element in alphaT3-1 cells. Mol Endocrinol. 2005, 19: 2780-2797. 10.1210/me.2004-0519.

    Article  CAS  PubMed  Google Scholar 

  42. Mack JT, Beljanski V, Soulika AM, Townsend DM, Brown CB, Davis W, Tew KD: "Skittish" Abca2 knockout mice display tremor, hyperactivity, and abnormal myelin ultrastructure in the central nervous system. Mol Cell Biol. 2007, 27: 44-53. 10.1128/MCB.01824-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Mack JT, Brown CB, Tew KD: ABCA2 as a therapeutic target in cancer and nervous system disorders. Expert Opin Ther Targets. 2008, 12: 491-504. 10.1517/14728222.12.4.491.

    Article  CAS  PubMed  Google Scholar 

  44. Falany CN: Sulfation and sulfotransferases. Introduction: changing view of sulfation and the cytosolic sulfotransferases. FASEB J. 1997, 11: 1-2.

    CAS  PubMed  Google Scholar 

  45. Gamage N, Barnett A, Hempel N, Duggleby RG, Windmill KF, Martin JL, McManus ME: Human sulfotransferases and their role in chemical metabolism. Toxicol Sci. 2006, 90: 5-22. 10.1093/toxsci/kfj061.

    Article  CAS  PubMed  Google Scholar 

  46. Xu H, Brill JA, Hsien J, McBride R, Boulianne GL, Trimble WS: Syntaxin 5 is required for cytokinesis and spermatid differentiation in Drosophila. Dev Biol. 2002, 251: 294-306. 10.1006/dbio.2002.0830.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (No. 30873299, 90409014, 30701123), and the Shanghai Leading Academic Discipline Project (S30304). We are grateful to Ping Wan, Lei-Jun Zhuang for their excellent technical assistance and professor Zhi-Yong Ming, Nian-Qing Ding and Ka Bian for helpful advice in the preparation of this manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yong-Qing Yang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YL was responsible for the qRT-PCR validation, bioinformatic analyses and draft of the manuscript. JG, WY, LY, & JW duplicated the rat asthmatic model and established four SAGE tag libraries, WY, ZQ, XY, & YY contributed to the design of the study and draft of the manuscript. All authors read and approved of the final manuscript.

Electronic supplementary material

12929_2008_46_MOESM1_ESM.eps

Additional file 1: Rat on the suspended shelf. For the convenient manipulation of the acupuncture points on the back, rat was placed on the suspended shelf, which could easily make it calm and stand still without anesthesia. (EPS 304 KB)

12929_2008_46_MOESM2_ESM.doc

Additional file 2: The intron-spanning primer pairs of genes of quantitative real-time PCR confirmation. The primer pairs of 3 differentially expressed genes of interest (Dusp1, S100A9, and MT-2) and the reference gene (GAPDH). (DOC 29 KB)

12929_2008_46_MOESM3_ESM.doc

Additional file 3: The comparisons of Cdyn and RR of the four groups. The table included the 10 min measurement of Cdyn and RR of the four groups. (DOC 43 KB)

12929_2008_46_MOESM4_ESM.xls

Additional file 4: The lists of differentially expressed tags of P CK/AS , P AS/ASAC , P CK/CKAC . The details of the 186, 130 and 144 differentially expressed tags (P < 0.05) of the PCK/AS, PAS/ASAC, and PCK/CKAC. This table included the tag sequence, tag copies of different libraries, UniGene number and the gene description. (XLS 84 KB)

12929_2008_46_MOESM5_ESM.xls

Additional file 5: The gene lists of DAVID gene functional classification of P CK/AS , P AS/ASAC , P CK/CKAC . The complete gene lists of 6 functional groups with enrichment scores equal to or higher than 1 were shown for the Figure 4. (XLS 63 KB)

12929_2008_46_MOESM6_ESM.xls

Additional file 6: The gene lists of the KEGG pathways. The complete gene lists of the matched KEGG pathways were demonstrated for the Table 3. (XLS 27 KB)

12929_2008_46_MOESM7_ESM.xls

Additional file 7: The gene lists of Venn diagram. The details of the different gene groups of Venn diagram, which included the tag sequence, tag copies of different libraries, UniGene number and the gene description. (XLS 28 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Yin, LM., Jiang, GH., Wang, Y. et al. Use of serial analysis of gene expression to reveal the specific regulation of gene expression profile in asthmatic rats treated by acupuncture. J Biomed Sci 16, 46 (2009). https://doi.org/10.1186/1423-0127-16-46

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1423-0127-16-46

Keywords