Abstract

Danqi pill (DQP) is a widely prescribed traditional Chinese medicine (TCM) in the treatment of cardiovascular diseases. The objective of this study is to systematically characterize altered gene expression pattern induced by myocardial ischemia (MI) in a rat model and to investigate the effects of DQP on global gene expression. Global mRNA expression was measured. Differentially expressed genes among the sham group, model group, and DQP group were analyzed. The gene ontology enrichment analysis and pathway analysis of differentially expressed genes were carried out. We quantified 10,813 genes. Compared with the sham group, expressions of 339 genes were upregulated and 177 genes were downregulated in the model group. The upregulated genes were enriched in extracellular matrix organization, response to wounding, and defense response pathways. Downregulated genes were enriched in fatty acid metabolism, pyruvate metabolism, PPAR signaling pathways, and so forth. This indicated that energy metabolic disorders occurred in rats with MI. In the DQP group, expressions of genes in the altered pathways were regulated back towards normal levels. DQP reversed expression of 313 of the 516 differentially expressed genes in the model group. This study provides insight into the multitarget mechanism of TCM in the treatment of complex diseases.

1. Introduction

Acute myocardial ischemia (MI) occurs when blood flow stops to a part of the heart causing damage to the myocardial tissue, and MI is one of the leading causes of death worldwide [13]. How to reduce MI-caused mortality is a major challenge to the entire medical community. Conventional managements of MI include intravenous thrombolysis, percutaneous coronary intervention, optimization of oxygenation, and pain control [4]. Danqi pill (DQP), composed of Salvia miltiorrhiza Bunge and Panax notoginseng, is a widely prescribed traditional Chinese medicine (TCM) in the treatment of a variety of cardiovascular diseases [5, 6]. DQP has been used as an alternative or complementary medicine in the prevention and treatment of MI in China. There are multiple components of DQP that have potential regulative effects on multiple targets in the treatment of cardiovascular diseases [710]. However, the overall regulative effect of DQP has not been explored yet. Investigations into the effects of DQP on global gene expressions will further our understanding on the mechanisms of TCM in treating complex diseases.

In recent years, transcriptome sequencing technologies have been developing very fast. Digital Gene Expression (DGE) technique is one of the RNA-sequencing methods for the analysis of differentially expressed genes in different samples [11]. DGE has more accuracy in quantification of gene expression levels compared with microarray technologies. It is able to provide quantitative readout of mRNA expression levels in samples [12]. Illumina’s sequencing platform was applied in our present study.

Acute MI model was induced by ligation of left anterior descending (LAD) coronary artery in rats. The rats were divided into three groups: sham-operated group, model group, and DQP-treated group. Gene expressions in the three groups were measured by RNA-sequencing technology, and differentially expressed genes and signaling pathways among different groups were analyzed. This study will expand our knowledge on the multitarget mechanism of DQP in the treatment of ischemic heart disease.

2. Materials and Methods

2.1. Grouping of Animals

Sprague-Dawley (SD) male rats, weighing 220 ± 10 g, were randomly divided into three groups: sham-operated group, acute MI model group, and DQP-treated group. Each group contained 10 rats. The animals were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd. This study complied with the China Physiological Society’s “Guiding Principles in the Care and Use of Animal” and got the approval of Animal Care Committee of Beijing University of Chinese Medicine.

MI model was induced in model and DQP group by direct LAD artery ligation as previously described [13, 14]. Briefly, SD rats were anaesthetized intraperitoneally with pentobarbital sodium (1%, 50 mg/kg). The left anterior descending coronary artery was ligated proximal to its main branching point with a 5–0 polypropylene suture. Sham-operated rats went through identical thoracotomy procedure but their coronary arteries were not ligated. After the operations, all of the rats were fed with a standard diet and were maintained on a 12 h light and dark cycle for 28 days. Animals in the DQP group were treated with concentrated DQP (purchased from Tongren Tang, Beijing, China) dissolved in pure water, with the daily dosage of 1.5 g/kg, for 28 consecutive days beginning from the day after the operation. At the end of the study, all rats were anaesthetized using pentobarbital sodium following an overnight fast. Rats were sacrificed and the left ventricle was carefully dissected to keep only the viable myocardium in the marginal zone of the infarct region in the DQP group and MI model group. The same region in the sham group was also dissected. Dissected heart tissues were frozen and stored in a freezer at −80°C before RNA extraction.

2.2. Echocardiographic Assessment of Heart Function

Echocardiography (Vevo 2100, Visual Sonics, Canada) was applied to assess the cardiac function-related parameters, including left ventricular ejection fraction (LVEF), left ventricular internal diameter at end-diastole (LVID;d) and at end-systole (LVID;s), and left ventricular fractional shortening (LVFS). LVEF and LVFS were calculated automatically by the software.

2.3. RNA Preparation

Heart tissues of the 3 rats in each group were homogenized in liquid nitrogen for RNA sequencing. Total RNA of the heart tissues were extracted using TRIzol Reagent® (Invitrogen, Carlsbad, CA), following the instruction of the manufacturer. Extracted RNA was treated with DNase to remove potential genomic DNA contamination. The quality of RNA was examined by Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA).

2.4. cDNA Library Preparation and RNA Sequencing

Sequence tag was prepared using Illumina’s Digital Gene Expression Tag Profiling Kit according to the manufacturer’s protocol. Briefly, mRNA was purified from total RNA by binding mRNA to a magnetic oligo(dT) bead. Oligo(dT) was then used as a primer to synthesize cDNA. The bead-bound cDNA was subsequently digested with restriction enzyme NlaIII, which recognizes and cuts off the CATG sites. The fragments apart from the 3 cDNA fragments connected to oligo(dT) beads are washed away and the Illumina adaptor 1 is ligated to the sticky 5 end of the digested bead-bound cDNA fragments. Mme I, a type of endonuclease with separated recognition sites and digestion sites, was used to cut at 17 bp downstream of the CATG site, producing tags with adaptor 1. After removing 3 fragments with magnetic beads precipitation, Illumina adaptor 2 is ligated to the 3 ends of tags. A tag library with different adaptors of both ends was thus formed.

After PCR amplification of 12 cycles, fragments were purified by a 6% Novex TBE PAGE Gel electrophoresis. In the end, the purified cDNA tags were sequenced using Illumina HiSeq™ 2000 at BGI-Shenzhen.

2.5. Analysis of Sequencing Data and Tag Mapping

During the quality control steps, Agilent 2100 Bioanalyzer was used in quantification and qualification of the sample library. Low-quality tags and tags of copy number less than two were filtered to produce clean tags. The clean tags were classified according to their copy number in the library. Clean tags were aligned to the rat reference sequences and annotated. Clean tag numbers corresponding to each gene were counted in each sample.

2.6. Analysis of Differentially Expressed Genes

The number of clean tags in each sample was normalized to transcripts per million (TPM) and TPM values were used to calculate fold change and false discovery rate across different groups. To avoid possible noise signal from high-throughput sequencing, the genes with TPM less than 3 in 2 or more samples were excluded. Analysis of differentially expressed genes among different groups was performed in R (version 3.0.2) with edgeR Bioconductor package [15]. EdgeR uses empirical Bayes estimation and exact tests based on the negative binomial distribution. The resulting values for all genes were corrected for multiple tests using a FDR adjustment. In this study, the fold change larger than 2 and FDR less than 0.01 were used to define the differentially expressed gene.

2.7. Gene Functional and Enrichment Analysis

The gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of all differentially expressed genes were carried out using the DAVID Functional Annotation Tool [16]. According to the transcription factor database (TFdb) (http://genome.gsc.riken.jp/TFdb/), transcription factors that had altered expressions induced by ischemia were selected. The target genes of these transcription factors were searched through databases, including Transcriptional Regulatory Element Database (http://rulai.cshl.edu/TRED), Human Transcriptional Regulation interaction database (http://www.lbbc.ibb.unesp.br/htri/index.jsp), and Transcription factor checkpoint (http://www.tfcheckpoint.org/index.php/search) [17]. Regulation networks were constructed using Cytoscape software [18].

2.8. Messenger RNA Expression by Quantitative Real-Time PCR

Real-time PCR was applied to validate differentially expressed genes in six samples in each group. First-strand cDNA was synthesized from total RNA with a RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific, USA, lot number: K1622) according to the manufacturer’s instruction. Quantitative real-time PCR assays were performed using C1000 Thermal Cycler PCR machine (Bio-Rad, USA). The reaction volume was 20 μl including 1 μl forward and reverse primer pairs, 2 μl cDNA, 10 μl FastStart Universal SYBR Green Master (Roche, Germany, lot number: 04913914001), and 7 μl RNase free water. The PCR procedures were as follows: 15 s at 95°C for denaturation and 1 min at 60°C for annealing and extension. Ct values were obtained after 40 cycles of reactions. Primer sequences of each gene were listed in Table 1. Ct values of targeted mRNA were normalized to the Ct values of GAPDH. Relative expressions of these genes were calculated by the 2−ΔΔCT method.

3. Results

3.1. Effects of DQP on Cardiac Function

Twenty-eight days after surgery, echocardiography showed that LVEF and LVFS of rats that underwent ligation in model group were downregulated significantly by 63.7% and 72.7%, compared with those of the sham group (Figure 1, ), indicating that cardiac function of rats in the model group was impaired and a MI model was established. LVID;s and LVID;d increased by 204.9% and 44.3% in the model group compared with those in the sham group. After treatment with DQP, the LVEF and LVFS were upregulated by 56.3% and 72.1%, compared with those in the model group (). LVID;s and LVID;d also decreased by 23.3% and 2.4% after treatment with DQP, suggesting that DQP could improve cardiac functions in the MI model (Figure 1).

3.2. High-Throughput Sequencing Data and Differentially Expressed Genes

In total, we obtained four million raw tags and over 3.7 high-quality clean tags from each sample. 51.8%, 47.8%, and 52.5% tags were mapped to annotated rat genomes in the DQP-treated, sham-operated, and model rats, respectively. About 5% of the tags were unknown ones. Altogether, 10,813 genes were detected, and 9537, 8907, and 9344 genes were detected and quantified in the DQP-treated, sham-operated, and model rats, respectively (Supplement Table 1).

Redundancy and heterogeneity are characteristics of mRNA expression. The majority of mRNAs have low expression level, whereas a minority of mRNAs has high abundance of expression. In this study, among the distinct clean tags in the nine libraries, the majority of distinct clean tags (60.21%–63.31%) had 2–5 copies. Only 3.54%–4.47% clean tags had more than 100 copies. The distribution of distinct clean tags showed a similar tendency, indicating that gene expression patterns among the libraries were similar.

Differentially expressed genes were analyzed using the following criteria: fold change > 2 and FDR < 0.01. The analysis showed that compared with the sham-operated group, 339 genes were upregulated and 177 genes were downregulated significantly in the model group, demonstrating that gene expression pattern was altered under ischemic condition. Compared with the model group, 886 genes were upregulated and 1082 genes were downregulated in the DQP-treated group.

Genes were expressed at different abundances, and expressions of some genes were altered remarkably in ischemic heart tissues. The top 20 upregulated and 20 downregulated genes in regard to abundance in the model compared with the sham-operated group and their expressions in DQP group were shown in Figure 2. NPPA, which encodes atrial natriuretic factor precursor, was upregulated greatly in ischemic heart. Insulin-like growth factor binding protein 7 (IGFBP7) was also upregulated. Other genes that were upregulated in great abundance encode collagen (COL3A1, COL1A1), light chain of myosin (MYL7), metallopeptidase (MMP2), and so forth. The top 20 downregulated genes in the model group are involved in fatty acid metabolism (ACAA2, ACSL1, ACADM, and DCI) and glycolysis (ALDOA, LDHB, MDH2, and PDHA1). Most of the products of the downregulated genes are located in the mitochondria. Among the top 20 genes that were upregulated or downregulated in the model group, DQP was able to regulate the expression of 19 and 14 genes back towards normal levels, respectively (Figure 2).

3.3. Biological Pathways and Processes Affected in Ischemia Heart

To investigate the global gene expression changes in ischemic heart in the model group, we studied the biological pathways and processes affected in ischemic heart by KEGG pathway enrichment and GO enrichment analysis. The upregulated KEGG pathways and biological processes mainly include complement and coagulation cascades, ECM-receptor interaction, dilated cardiomyopathy, response to wounding, coagulation, immune response, extracellular matrix organization, focal adhesion, skeletal system development, and cell growth (Tables 2 and 3). The downregulated biological processes and pathways mainly involve in nutrient metabolism and energy supply, such as fatty acid metabolism, pyruvate metabolism, amino acid metabolism, carbohydrate catabolism, citrate cycle, and oxidation reduction (Tables 2 and 3). This analysis illustrated that the energy supply of the heart was seriously compromised and the myocardial remodeling took place to compensate for the short supply of energy in the ischemic process.

3.4. Genes and Signaling Pathways Regulated by DQP

Expressions of 516 genes were significantly affected by MI. In the DQP-treated group, expressions of 313 out of these 516 genes were reversed by DQP. These overlapped genes are considered as DQP-responsive genes. Of the 339 genes that were upregulated in the model group, 212 were significantly downregulated in the DQP treatment group. These genes were enriched in infection and immunity, cell growth, extracellular matrix deposition, and so forth (Table 4). 46% of proteins encoded by the downregulated genes were located in extracellular space. Among the 177 downregulated genes in the model group, DQP significantly reversed expressions of 101 genes. The enriched pathways and biological processes were involved in energy embolisms, such as carbon metabolism, fatty acid degradation, and pyruvate metabolism (Table 5). 25% of the proteins encoded by the downregulated genes were located in the mitochondrion. These results demonstrated that DQP can regulate multiple pathways altered under ischemic stimulus.

3.5. Transcriptional Regulations by DQP

Transcription factors play critical roles in regulating expressions of target genes. According to the transcription database, we searched transcription factors that had altered expressions induced by ischemic stimuli. Among the dysregulated 516 genes in the MI model group, 12 were transcription factors, including CEBPD, CREB3, BCL3, SIX4, NFATC4, MAFF, ERF, NFIB, MITF, AEBP1, EGR2, and PRRX2. The target genes of these transcription factors were also searched out through databases, and the regulative networks were constructed (Figure 3). Expressions of SIX4, NFIB, and MITF were downregulated in the model group compared with the sham-operated group and the expressions of the nine genes were downregulated. The deregulated transcription factors play roles in inflammation, lipid metabolism, mitochondrial metabolism, and cardiac fibrosis. DQP treatment regulated expressions of the 12 transcription factors towards normal levels.

3.6. Real-Time PCR Validation Analyses

To validate the RNA-sequencing results, we measured the expression patterns of 10 genes by quantitative real-time PCR analyses. Five genes involved in fatty acid metabolism pathways were chosen for validations. Real-time PCR results showed that expressions of all of these five genes were downregulated in the model group and upregulated by DQP treatment, consistent with RNA-sequencing results (Figure 4(a)). Expressions of another five genes involved in ventricular remodeling, cellular apoptosis, and inflammation were also validated. Real-time PCP results showed that expressions of these five genes were upregulated in the model group and downregulated by DQP treatment, which were also consistent with RNA-sequencing results (Figure 4(b)).

4. Discussion

MI is one of the cardiovascular conditions that threaten people’s health, and TCM has been shown to be effective in attenuating symptoms of MI. The aim of this study was to systematically analyze the effects of DQP on altered global gene expression patterns induced by MI. An MI rat model was induced by coronary artery ligation, and Illumina’s RNA-seq platform was applied in this study. The results showed that heart functions were impaired in the MI model rat and DQP treatment protected heart functions at 28 days after the operation. 10,813 genes were detected and quantified in the heart tissues of three groups of rats. Compared with the sham-operated group, 516 genes were differentially expressed in the MI model group. In the DQP-treated group, expressions of 313 of the altered 516 genes in MI model were regulated back towards normal levels. DQP treatment could regulate the altered signaling pathways in ischemic heart tissues and exert an overall regulative effect on multiple targets in signaling pathways.

DQP is composed of Salvia miltiorrhiza Bunge and Panax notoginseng. Our results showed that DQP had a remarkable cardioprotective effect, as demonstrated by improved LVEF and LVFS. Previous studies have demonstrated that DQP has cardioprotective effects, and the mechanism is related to its regulation on energy metabolism pathways [5, 19, 20]. DQP has also been shown to have anti-inflammatory effect in heart failure model of animals [6]. Salvia miltiorrhiza Bunge and Panax notoginseng are the major components of DQP, and their effects on cardiovascular diseases have been investigated by numerous studies [2126]. The major potential effective components of Salvia miltiorrhiza Bunge and Panax notoginseng include salvianolic acids and Panax notoginseng saponins, and they may exert cardioprotective effects through multiple targets in a synergistic way [27, 28]. The effects of DQP on global gene expression patterns in ischemic heart tissues were investigated in this study.

Myocardial ischemic stimulus induced global gene expression changes in the border zone of the heart tissues. In response to a short supply of oxygen and nutrient, several signaling pathways were downregulated, including fatty acid metabolism, amino acid degradation, pyruvate metabolism, PPAR signaling pathway, and citrate cycle. The proteins encoded by these genes were enriched in mitochondria. The top genes with the greatest abundance reductions in the MI model group were enriched in fatty acid metabolism and glycolysis (Figure 2). For example, long-chain acyl-coenzyme A synthetase 1 (ACSL1) could interact with fatty acid transport proteins (FATP) and contribute to the efficient cellular uptake of long-chain fatty acids through vectorial acylation. Inhibition of ACSL1 activity in ischemia heart impairs fatty acid uptake [29]. Acyl-coenzyme A dehydrogenase for medium-chain fatty acids (ACADM) catalyzes the initial step of the mitochondrial fatty acid beta-oxidation pathway [30]. Reduced expressions of these genes suggested that energy metabolism is disrupted in ischemic heart tissues. Among the 177 downregulated genes, DQP treatment upregulated expressions of 101 of these genes, demonstrating that DQP could modulate energy metabolism under ischemic conditions. Our previous study showed that the major components of DQP, salvianolic acids, and Panax notoginseng saponins, had effects on energy metabolism [19].

In the MI model group, expressions of 339 genes were upregulated compared with the sham group. According to gene ontology terms, these genes were enriched in response to wounding, extracellular matrix organization, skeletal system development, bone development, regulation of cellular component size, and so forth. Expressions of gene encoding collagens (COL3A1, COL1A1) were greatly increased. Gene encoding natriuretic peptide A (NPPA) implicated in the control of extracellular fluid volume and electrolyte homeostasis was also greatly upregulated. Levels of metalloproteinases (MMP2, MMP23) were upregulated in the model group [31]. These results demonstrated that ventricular remodeling had occurred in response to MI. DQP treatment reversed expressions of these genes towards normal levels. Oxidative stress could produce reactive oxygen species (ROS). ROS could damage endothelial cells and heart tissues, contributing to myocardial fibrosis [32, 33]. Superoxide dismutase (SOD) is a major defense mechanism against ROS [32]. Mitochondrial SOD2 was repressed in ischemic rat models and was upregulated by DQP, suggesting that DQP has a protective antioxidative effect on the heart. Components of DQP, such as salvianolic acid A, salvianolic acid B, ginsenoside Rg1, and so on, have been shown to have antioxidative effects and attenuate hypertrophy [3437]. Transcription factors involved in inflammation, lipid metabolism, mitochondrial metabolism, and cardiac fibrosis were also regulated back towards normal levels by DQP.

In conclusion, this study demonstrated that DQP could exert cardioprotective effects in the MI model by regulating global gene expression pattern. One limitation of this study is that the sample size was small. The implications of this study warrant further studies with larger sample size to validate the synergistic effects of Chinese medicine. Furthermore, this kind of study will provide insight into the effective components of Chinese medicine and their respective targets in treating complex diseases.

Disclosure

The manuscript was presented as an abstract in the 26th Great Wall International Congress of Cardiology.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this manuscript.

Authors’ Contributions

Qiyan Wang, Hui Meng, and Qian Zhang contribute equally to this paper.

Acknowledgments

This work was financially supported, in part, by the grants from the National Natural Science Foundation of China (nos. 81503379, 81673712, 81530100, and 81473456), Fok Ying Tung Education Foundation (no. 151044), Beijing Nova program (Z171100001117028), and the excellent young scientist foundation of BUCM (nos. 2016-JYB-XJ003, 2015-JYB-QNJSZX001, and 2015-JYB-XYQ001).

Supplementary Materials

Supplemental Table 1: basic information of RNA-sequencing data. (Supplementary Materials)