Full-length transcriptome analysis reveals the mechanism of acupuncture at PC6 improves cardiac function in myocardial ischemia model
Chinese Medicine volume 16, Article number: 55 (2021)
The pathological process of myocardial ischemia (MI) is very complicated. Acupuncture at PC6 has been proved to be effective against MI injury, but the mechanism remains unclear. This study investigated the mechanism that underlies the effect of acupuncture on MI through full-length transcriptome.
Adult male C57/BL6 mice were randomly divided into control, MI, and PC6 groups. Mice in MI and PC6 group generated MI model by ligating the left anterior descending (LAD) coronary artery. The samples were collected 5 days after acupuncture treatment.
The results showed that treatment by acupuncture improved cardiac function, decreased myocardial infraction area, and reduced the levels of cTnT and cTnI. Based on full-length transcriptome sequencing, 5083 differential expression genes (DEGs) and 324 DEGs were identified in the MI group and PC6 group, respectively. These genes regulated by acupuncture were mainly enriched in the inflammatory response pathway. Alternative splicing (AS) is a post-transcriptional action that contributes to the diversity of protein. In all samples, 8237 AS events associated with 1994 genes were found. Some differential AS-involved genes were enriched in the pathway related to heart disease. We also identified 602 new genes, 4 of which may the novel targets of acupuncture in MI.
Our findings suggest that the effect of acupuncture on MI may be based on the multi-level regulation of the transcriptome.
Myocardial ischemia (MI) is a severe cardiovascular disease that is the leading cause of morbidity and mortality . Since the pathological process and mechanism of MI injury is complicated, which involved in cardiac fibrosis , inflammation , angiogenesis , apoptosis , the regulation of this disease is also multifaceted. As a complementary and alternative medical approach, acupuncture has been widely used to relieve the symptoms of angina and palpitations [24, 45]. Many studies have been proved the protective effect of acupuncture on MI [8, 17]. However, because of the extensive regulation effect of Acupuncture on various pathological processes of MI injury, the entire mechanism is also unclear.
Gene regulation is a fundamental process in developing and disease progression, and transcriptome research a necessary tool for understanding the gene expression profile. Therefore, we could like to research the mechanism of acupuncture's effect on MI from the transcriptome. The regulation of acupuncture on ischemia–reperfusion has been proved by some researches in transcriptome analysis . However, the effect of acupuncture on genome-wide gene expressions in myocardial ischemia models has not been reported. Moreover, most researches on the transcriptome of acupuncture effect on heart diseases are relatively limited and primarily focuses on differential genes. In this study, full-length transcriptome sequencing was used to obtain a more comprehensive transcriptome profile to analyze the pathogenesis of MI and the regulatory effect of acupuncture from a more perspective. In fact, the transcriptome is complex, such as alternative splicing, LncRNA, and polyadenylation, which can also affect gene function [33, 46]. As a primary mechanism for producing protein diversity, alternative splicing has been demonstrated to be involved in proliferation, apoptosis, angiogenesis, and so on . Alternative splicing also plays an essential role in heart disease [9, 25]. Research proved that alternative splicing in titin (TTN) could modify sarcomeric stiffness to influence cardiac function . Mutations in a splicing factor of TTN, RBM20, have been proved involved in dilated cardiomyopathy . However, the role of alternative splicing in myocardial ischemia has not been well understood. Based on full-length transcriptome, alternative splicing can be detected, which allows us to study the role of alternative splicing in MI. Moreover, some new genes which structure and function have not been reported can also be identified by full-length transcriptome, these new genes may be the critical targets for acupuncture.
To investigated the mechanism of acupuncture against MI injury from a more comprehensive perspective, we studied the transcriptome from the aspects of differential genes, alternative splicing, and new genes. This study aimed to investigate the mechanism of acupuncture against MI from the multi-level transcriptome. This research will help find new therapeutic targets for MI and a unique perspective for acupuncture research.
Experimental animals and grouping
Six to 8-week-old SPF (specific pathogen-free) C57BL/6 male mice (20–23 g) were purchased from Byrness Weil biotech Ltd (Hunan, China) and used in the experiments. Animals were fed adaptively for one week in the environment with a temperature of 23–25 °C and relative humidity of 45–55%. The study was performed according to the guidelines of the National Institutes of Health Animal Care and Use of Laboratory Animals (Bethesda, MD, USA). All the mice were divided into the Control group, MI group, and PC6 group randomly.
Induction of MI in mice
The MI models were produced by ligating the left anterior descending coronary artery (LAD). Mice were anesthetized by inhalation of 5% isoflurane and maintained with 2% isoflurane (RWD, Shenzhen, China). The skin was incised along the left side of the sternum, and the muscles were bluntly separated, then incised between the 3rd and 4th intercostal spaces and permanent LAD ligation. In the control group, the same procedure was performed except for the LAD ligation. The Lead II Electrocardiogram was monitored before and after the operation.
The PC6 group mice underwent acupuncture treatment after 24 h of MI operation bilateral PC6 acupoint once a day under awake, for a total of 5 days. Fold the acupuncture needles into an “L” shape, then inserted needles into PC6 acupoint about 1–2 mm. After the needles were fixed with adhesive tape, mice were put into the cages and moved freely. Then we pressed the needles every 5 min to ensure the needles would not fall off and strengthen the stimulation of the acupuncture intervention. The needles were removed after 20 min of treatment.
Electrocardiogram (ECG) and echocardiography
Before and after the MI operation, the electrocardiogram (lead II) of each group of mice was recorded by the Power Lab system (AD Instruments, Australia). Then, 5 min after LAD surgery, the ST-segment was elevated to the myocardial ischemic indication. 5-day after MI, the cardiac function was assessed by transthoracic echocardiography (Beijing Yeeran Technology Co., Ltd). Mice have maintained anesthesia with 1.5% isoflurane, and heart rate was maintained at about 450–550 bpm during the examination. The indices included left ventricular end-diastolic diameter (LVEDD) and left ventricular end-systolic diameter (LVESD). The left ventricular ejection fraction (LVEF) and left ventricular fractional shortening (LVFS) were calculated according to LVEDD and LVESD. The LVEF and LVFS were calculated as follows: EF = [(LV Vol: d-LV Vol: s)/LV Vol: d] × 100%; FS = [(LVEDDLVESD)/LVEDD] × 100%.
Triphenyitertrazolium chloride (TTC) staining
TTC staining was used to assess the infarct area in heart tissue. Mice were sacrificed five days after acupuncture treatment, and the heart tissues were frozen for 20 min at – 20 °C and cut into five pieces below the ligature. The slices were stained in 2% TTC solution (Sigma Chemical Co., USA), incubated at 37 °C for 20 min and fixed in 4% paraformaldehyde for 4 h. The infarct area was assessed by the use of Image-Pro Plus 6.0.
Collection of tissue and serum samples
After 5 days of acupuncture treatment, the mice were sacrificed under deep anesthesia with pentobarbital sodium (100 mg/kg) administered by intraperitoneal injection. Blood was collected from the orbital vein and then centrifuged for 10 min at 4000 rpm, 4 °C, to obtain the serum. The serums were stored at − 80 °C for myocardial enzymes and inflammatory cytokines determination.
Determination of serum myocardial enzymes and inflammatory cytokines
The levels of cTnT, cTnI, TNF-α, IL-1β, and IL-6 were quantified with the ELISA kits according to the manufacture’s protocol (E-EL-M1801c, E-EL-M1203c, E-EL-M0049c, E-EL-M0037c and E-EL-M2453c, Elabscience, China).
RNA extraction, library construction and sequencing
Total RNA was extracted from 9 heart samples and prepared for cDNA libraries using the protocol provided by Oxford Nanopore Technologies (ONT). Briefly, SuperScript IV First-Strand Synthesis System (Invitrogen) was used for full-length mRNA reverse transcription and following cDNA PCR for 14 circles with LongAmp Tag (NEB). The PCR products were then subjected to FFPE DNA repair and end-repair (NEB) steps and following adaptor ligation using T4 DNA ligase (NEB). Agencourt XP beads were used for DNA purification according to ONT protocol. The final cDNA libraries were added to FLO-MIN109 flowcells and were run on the PromethION platform at Biomarker Technology Company (Beijing, China).
Analysis for RNA-seq data
Raw reads were first filtered with minimum average read quality score = 7 and minimum read length = 500 bp. Ribosomal RNA was discarded after mapping to the rRNA database. Next, full-length, non-chimeric (FLNC) transcripts were determined by searching for primer at both ends of reads. Clusters of FLNC transcripts were obtained after mapping to the reference genome with mimimap2, and consensus isoforms were obtained after polishing within each cluster by pinfish. Transcripts were validated against known reference transcript annotations with gffcompare. AS events including IR, ES, AD, AA, and MEE were identified by the AStalavista tool. Differential expression analysis of two groups was performed using the DESeq R package (1.18.0). DESeq provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The resulting P values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with a P value < 0.05 and fold change ≥ 1.5 found by DESeq were assigned as differentially expressed. Gene Ontology (GO) enrichment analysis of the differentially expressed genes (DEGs) was implemented by the GOseq R packages based on Wallenius non-central hypergeometric distribution, which can adjust for gene length bias in DEGs. KEGG analysis was used KOBAS software to test the statistical enrichment of differential expression genes in KEGG pathways.
Hematoxylin–eosin (HE) staining
The heart tissues in the control, MI, and PC6 group were collected and fixed with buffered 4% formalin for two days. All the samples were sectioned at a thickness of 5 µm and stained with hematoxylin and eosin for histopathologic observation by light microscopy (Olympus Corporation, Japan).
Real-time fluorescence quantitative polymerase chain reaction (RT-qPCR)
Total mRNA was isolated from the heart tissues using MolPure® TRIeasy™ Plus Total RNA Kit (YEASEN, Shanghai, China), and cDNA was prepared using Hifair® III 1st Strand cDNA Synthesis SuperMix for qPCR (YEASEN, Shanghai, China) according to the manufacturer’s instructions. The mRNA levels were assessed on a bio-rad cfx maestro system Real-time PCR system (Bio-Rad, USA) by qPCR using Hieff® qPCR SYBR® Green Master Mix (YEASEN, Shanghai, China). The relative expression of mRNA was calculated by △△Ct according to standard methods. The primer sequences as follows:
Cxcl2, forward: 5′-CCAACCACCAGGCTACAGG-3′, reverse: 5′-GCGTCACACTCAAGCTCTG-3′; Saa3, forward: 5′-TGCCATCATTCTTTGCATCTTGA-3′, reverse: 5′-CCGTGAACTTCTGAACAGCCT-3′; Smad1, forward: 5′-GGGGGATCCGTAATGTGACAAGTTTATTTTC-3′, reverse: 5′-TTTGCGGCCGCTCAAGATACAGATGAAATAG-3′; Arg1, forward: 5′-AGCTCTGGGAATCTGCATGG-3′, reverse: 5′-ATGTACACGATGTCTTTGGCAGATA-3′; GAPDH, forward: 5′-TGTGTCCGTCGTGGATCTGA-3′, reverse: 5′-CCTGCTTCACCACCTTCTTGA-3′.
Statistical analysis was completed using SPSS 19.0 software (IBM, Chicago, USA). Values are presented as mean ± SEM. The two-group comparison was determined using the unpaired Student’s t-test, and one-way ANOVA made multiple-group comparisons. P < 0.05 was considered statistically significant.
Acupuncture treatment at the PC6 acupoint effectively protected the myocardium from MI injury
To investigate the effects of acupuncture at PC6 acupoint on myocardial injury, we first assessed the cardiac function and infarct area. Echocardiography results indicated that both EF and FS significantly decreased after MI compared with the control group. After treatment at PC6 acupoint, both EF and FS increased compared to the MI group (Fig. 1A–C). The infarct area was assessed at 5d after acupuncture treatment by TTC staining. The results showed that acupuncture treatment significantly decreases the myocardial injury size (Fig. 1D, E). The levels of myocardial-specific serum enzymes, including cardiac troponin T (cTnT) and cardiac troponin I (cTnI), which reflect acute myocardial injury, were measured by ELISA. The results indicated that the levels of cTnT and cTnI were increased after MI operation, acupuncture significantly decreased the levels of serum enzymes (Fig. 1F, G).
Acupuncture at the PC6 acupoint altered genes expressions genome-widely
To explore the mechanism of acupuncture treatment on MI, we completed gene expression profiling of different groups using Nanopore. RNA-seq results showed that compared to the control group, 5083 genes were differentially expressed in the MI group; of these, 2774 genes were up-regulated, and 2309 genes were down-regulated (Fig. 2A; Additional file 1: Table S1). We also found 324 differentially expressed genes (DEGs) in the PC6 group compared to the MI group; 126 differentially expressed genes were up-regulated, and 198 genes were down-regulated (Fig. 2B; Additional file 2: Table S2).
We further analyzed the gene expression patterns by clustering ONT data. The heatmaps were created with the up-regulated genes and down-regulated genes in the control group, MI group, and PC6 group. As Fig. 2C showed, the red color was higher, and the green color illustrated lower expression genes. Then, we analyze overlaps of the differential genes with Venn Diagrams (Fig. 2D). In the 2773 up-regulated DEGs of the MI group, 65 DEGs were down-regulated by acupuncture treatment. On the other side, in the 2308 down-regulated DEGs of the MI group, 19 DEGs were reversed by acupuncture treatment. We speculate that these 84 DEGs may be the key targets of the acupuncture effect against myocardial ischemia. Additionally, the DEGs co-regulated by MI and PC6 were subjected to protein–protein interaction (PPI) network analysis based on the STRING database (Fig. 2E). And it showed that proteins were clustered into four categories: inflammatory response, Regulation of actin cytoskeleton, Type I interferon-mediated signaling pathway, and Transcription-dependent tethering of RNA polymerase II gene DNA at the nuclear periphery. Our data set offered an insight into the interactions of proteins in the MI pathological processes and provided a guidance for the target genes of acupuncture.
The validations of inflammatory response
We found that most of the DEGs are primarily enriched in the inflammatory response signaling pathway so that we observed inflammation by HE staining, Elisa and qPCR. HE staining showed a massive accumulation of infiltrated cells in the MI group, and the inflammatory cell infiltration was significantly decreased by acupuncture (Fig. 3A). We also detected the levels of inflammatory cytokines in serum. As revealed in Fig. 3B–D, the levels of TNF-α (tumor necrosis factor-α), IL-1β (interleukin-1β), and IL-6 (interleukin-6) in serum were increased in the MI group compared with the control group, and the levels of these inflammatory cytokines were decreased by acupuncture. Then, we verified some of the DEGs, including Cxcl2, Saa3, Smad1, and Arg1, which involved inflammatory response by qPCR. The qPCR confirmed a similar tendency as shown in RNA-seq data in Additional files 1, 2: Tables S1, S2 (Fig. 3E–H). These results proved the RNA-seq data was reliable and confirmed the regulation effect of acupuncture on inflammation.
Integrating alternative splicing events in myocardial ischemia
Alternative splicing (AS) is one of the primary mechanisms for producing protein diversity. Previous researches have reported that AS is involved in the development of myocardial ischemia . We also identified a large number of AS in each group. There are five main types of AS, including Exon skipping (ES), Intron retention (IR), Alternative 3′ splice site (A3SS), Alternative 5′ splice site (A5SS), and Mutually exclusive exon (ME) (Fig. 4A). In all samples detected by RNA-seq, 8237 AS events associated with 1994 genes were found. We detected 4518 ES-types AS events involving 859 genes, 1394 A3SS-type AS events containing 402 genes, 1148 A5SS-type AS events involving 398 genes, 956 RI-types AS events involving 280 genes, 221 ME-types AS events involving 55 genes, as shown in Fig. 4B. The results indicated that one gene might have several types of mRNA splicing events; some genes may be expressed by 3 AS types. The proportion of the five types of AS were shown in Fig. 4C. We found that more than half of the AS events were ES events. Then we analyzed differentially expressed genes and AS events involved genes (Fig. 4D). The results showed many DEGs underwent AS events, indicating that many genes may play a role in the pathological process of MI and acupuncture treatment through alternative splicing. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of AS involved genes found that the pathways including mitochondrion organization, proteasomal protein catabolic process, RNA splicing, mRNA splicing via spliceosome, and metabolic process were the central pathways in which the AS involved genes enriched (Fig. 4E).
The analysis of differential alternative splicing events in myocardial ischemia
Differential AS analysis revealed that 347 differential AS events were identified in this research, of which 174 were differentially upregulated AS events, 173 were differentially downregulated AS events, the numbers of upregulated AS events were equal to downregulated AS events (Fig. 5A, B). There were 129 differentially upregulated AS events and 147 differentially downregulated AS events in the MI group compared to the control group. There were 45 differentially upregulated AS events and 26 differentially downregulated AS events in the PC6 group compared to the MI group (Fig. 5A, B). To further explore the genes involved in differential AS events, GO annotations and KEGG pathways were used. These genes were annotated through three independent ontologies in the GO database, including biological process (BP), molecular function (MF), and the cellular component (CC). The top 10 GO terms of BP, MF, and CC were showed in Fig. 5C. Through the KEGG enrichment analysis, we found that the top 10 pathways which these genes enriched including Ribosome, HIF-1 signaling pathway, Hypertrophic cardiomyopathy, Dilated cardiomyopathy (Fig. 5D). The results suggested that these genes involved in differential AS events were associated with heart disease.
The analysis of new genes detected in full-length transcriptome
Full-length sequencing can accurately identify the structure of transcripts. Gffcompare is used to align transcripts obtained by full-length sequencing to known transcripts of the genome, so novel genes and transcripts can be obtained to supplement genome annotation. In this study, we found 602 new genes. Most of these new genes were shared by all 9 samples, which suggested that these new genes can be stably expressed in all samples (Fig. 6A). We also noticed that some new genes were expressed only in the MI group and PC6 group, indicated the gene expression pattern was changed under the pathological condition (Fig. 6A). The heatmap was created with all the new genes of each group. As Fig. 6B showed, the new gene expression patterns of control were significantly different with the MI group and PC6 group; however, the MI group and PC6 group had similar new gene expression patterns. We analyzed overlaps of differential new genes with Venn Diagrams. And found that 3 up-regulated new genes in the MI group were down-regulated by acupuncture; only 1 down-regulated new gene in the MI group was up-regulated by acupuncture (Fig. 6C). The expression and sequence of these 4 new genes were showed in Fig. 5D and Additional file 3: Table S3.
Myocardial ischemia is the leading cause of disability and death worldwide . As a traditional Chinese treatment, acupuncture has been proven to have a protective effect against myocardial ischemic injury [38, 39, 41]. Previous researches have indicated that acupuncture at PC6 can attenuate ischemia injury by reducing apoptosis , correcting arrhythmia , improving energy metabolism . Our results found that acupuncture at PC6 could protect the heart from myocardial ischemia injury. In this study, MI models were performed by ligating the LAD, which was a classic mouse MI model . And then applied acupuncture intervention to these MI models. We evaluated myocardial injury by using ECG, echocardiography, TTC staining, and myocardial enzyme levels. Echocardiographic EF and FS mainly estimated systolic cardiac function ; these results showed that both EF and FS were significantly decreased in the MI group. TTC staining showed a significant myocardial infraction and high levels of cTnT and cTnI in MI models. Acupuncture at PC6 not only increased EF and FS, and decreased the injury area and the levels of the myocardial enzyme. These findings were consistent with the previous research [8, 17, 42].
MI dramatically modified the gene expression pattern of the tissue environment, and gene expression is fundamental and essential in disease progression and treatment. Therefore, transcriptome research is one of the necessary tools for understanding life processes. To investigated the mechanism of acupuncture against MI injury from a more comprehensive perspective, we studied the transcriptome from the aspects of differential genes, alternative splicing, and new genes by full-length transcriptome. In previous studies on MI, Second-generation sequencing technology (SGS) has been widely used [11, 23, 40, 48]. However, the short-read RNA-seq method of SGS is limited to identify and quantify complex transcript isoforms. Identification of all transcripts produced by a gene is essential to research the disease progression. Thus, the ability to produce longer reads would be a better choice to understand transcriptome complexity, including AS, polyadenylation, and fusion genes [1, 30, 43]. This study applied the third-generation sequencing technology (TGS), commercialized by Oxford Nanopore Technologies (ONT) MinION, to research the transcriptome complexity in MI and acupuncture treatment. Compared with SGS, the read length of TGS can be about 1.5 kbp, and the transcriptome information is more comprehensive.
In this study, we found that the gene expression patterns were significantly changed after MI, acupuncture also could regulate some of these genes. By analyzing differentially expressed genes in MI and PC6 group, we found 5083 DEGs (2774 up-regulated and 2309 down-regulated) in the MI group compared to the control group; 324 DEGs (126 up-regulated and 198 down-regulated) in the PC6 group compared to the MI group. The analysis of DEGs, which co-regulated in two groups, indicated that these genes were mainly enriched in Inflammatory response, Regulation of actin cytoskeleton, Type I interferon-mediated signaling pathway, and Transcription-dependent tethering of RNA polymerase II gene DNA at the nuclear periphery. The effects of acupuncture on some of these signaling pathways have been reported in previous research [36, 42], suggesting that these signaling may be the key targets of acupuncture. Since most of these co-regulated genes are involved in the inflammatory response, and the inflammatory response plays a crucial role in scar formation and myocardial remodeling in MI  so that we verified the inflammatory reaction in MI and the regulatory effect of acupuncture on it from the morphological, gene and protein levels. We firstly verified the DEGs involved in the inflammatory response by qPCR. These genes, including Cxcl2, Arg1, and Saa3 play a key role in inflammatory response [13, 18, 21]. The role of Cxcl2 and Smad1 in myocardial ischemia also have been reported [12, 26, 32]. We found that qPCR results were consistent with RNA-seq, indicating that our RNA-seq data are reliable. HE staining showed that the inflammatory infiltration of the ischemic myocardium was reduced by acupuncture. TNF-α, IL-1β, and IL-6 are the major pro-inflammatory cytokines that aggravate the inflammatory response after ischemia injury. Our results also confirmed that acupuncture could decrease the level of these inflammatory cytokines in serum, further proving the anti-inflammatory effect of acupuncture.
In recent years, alternative splicing, another factor that causes protein diversity, has also been widely studied in many fields [4, 27, 33]. The pre-mRNA generated by gene transcription has many splicing methods. Different exons are selected to produce different mature mRNAs, then they are translated into different proteins, and these proteins can bring diversity to biological traits [5, 47]. The process of pre-mRNA treatment is called alternative splicing. There are five main types of AS, including ES, RI, ME, A3SS, A5SS, and ES, the most common AS in mammals. In recent years, there are many studies on AS in the cardiovascular field. In the heart, alternative splicing of sarcomeric genes, ion channels, and cell signaling proteins can cause cardiomyopathies and arrhythmias [19, 20]. Research on cardiovascular disease found that abnormalities in alternative splicing of apoptotic genes, including the Bcl-2 family, Caspases, Binp3, and Nix, are associated with many heart diseases . Alternative splicing disorder of CaMKIIδ can cause LDH release, ATP reduction, and ROS accumulation to aggravate hypoxia-reoxygenation injury of cardiomyocytes . Alternative splicing factor-regulated alternative splicing of CaMKIIδ also plays a vital role in pathological cardiac remodeling in heart failure . Although alternative splicing is essential in heart disease research, its role in myocardial ischemia is still not well understood, and the effect of acupuncture on alternative splicing has not been reported. This study first analyzed AS events involved genes and found that these genes were mainly enriched in pathways associated with RNA splicing and energy metabolism. Moreover, a considerable number of AS events involved genes showed no difference in expression through the analysis of overlap. These genes enhanced the diversity of protein by AS events but do not affect gene expression. Analysis of differential alternative splicing found that DAS involved genes were enriched in signaling pathways associated with heart disease, including the hypoxia-inducible transcription factor-1 (HIF-1) signaling pathway, Hypertrophic cardiomyopathy, Dilated cardiomyopathy. These results indicated that the regulation of AS events might also be a possible target of acupuncture. In addition, HIF-1 is also a global regulator of macrophage and neutrophil inflammatory [49. It also have been reported to play a role in inflammatory response after MI by interacting with Pyruvate kinase isozyme type M2 (PKM2) to influence macrophage M2 polarization . This finding is consistent with the regulation of inflammation by DEGs. Since DEGs and alternative splicing both have a relationship with inflammation after MI, the inflammatory response may be a target to research the interaction between DEGs and alternative splicing.
Full-length sequencing can identify new genes. In this study, we found a total of 602 new genes detected in 9 samples sequenced. More than half of these new genes were expressed in all samples, suggest that these genes were relatively stable. Some genes were only expressed in the MI and/ or PC6 groups, indicating that these genes may be induced under the pathological condition. Finally, we detected 89 new genes were up-regulated in the MI group; 3 of these genes were down-regulated in the PC6 group; 83 new genes were down-regulated in the MI group, 1 of these genes was up-regulated in the PC6 group. The structure and function of these genes have not been reported in previous research, but the expression of these genes changes in MI condition and also can be regulated by acupuncture, we speculate that these 4 new genes may be a novel target of acupuncture.
This study provided a more comprehensive profile of the transcriptome regulated by acupuncture in MI models; however, in this study, only the DEGs were validated. Further research on critical genes regulated by alternative splicing and new genes is need.
We demonstrate that acupuncture at PC6 could improve cardiac function and reduce infarct injury. By using Full-length transcriptome, we found that acupuncture can regulate differential expression genes and alternative splicing which may affect the structure and function of the protein. We also detected 4 new genes, which may be the target of acupuncture in the MI situation. Based on these results, we concluded that acupuncture might against myocardial ischemia through the regulation of multifaceted of the transcriptome.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Differential expression genes
Cardiac troponin I
Left ventricular end diastolic diameter
Left ventricular end systolic diameter
Alternative 3′ splice site
Mutually exclusive exon
Serum amyloid antigen 3
Mothers against decapentaplegic homolog 1
Tumor necrosis factor-α
Left anterior descending
Cardiac troponin T
Left ventricular ejection fraction
Left ventricular fractional shortening
Alternative 5′ splice site
CXC motif chemokine ligand 2
Abdel-Ghany SE, Hamilton M, Jacobi JL, Ngam P, Devitt N, Schilkey F, Ben-Hur A, Reddy ASN. A survey of the sorghum transcriptome using single-molecule long reads. Nat Commun. 2016;7:1–11. https://doi.org/10.1038/ncomms11706.
Chen H, Jing X-Y, Shen Y-J, Wang T-L, Chen Ou, Sheng-Feng Lu, Cai Y, et al. Stat5-dependent cardioprotection in late remote ischaemia preconditioning. Cardiovasc Res. 2018;114(5):679–89. https://doi.org/10.1093/cvr/cvy014.
Cheng Y, Feng Y, Xia Z, Li X, Rong J. ω-Alkynyl arachidonic acid promotes anti-inflammatory macrophage M2 polarization against acute myocardial infarction via regulating the cross-talk between PKM2, HIF-1α and INOS. Biochim Biophys Acta Mol Cell Biol Lipids. 2017;1862(12):1595–605. https://doi.org/10.1016/j.bbalip.2017.09.009.
Dlamini Z, Tshidino SC, Hull R. Abnormalities in alternative splicing of apoptotic genes and cardiovascular diseases. Int J Mol Sci. 2015;16(11):27171–90. https://doi.org/10.3390/ijms161126017.
Dvinge H, Jamie G, Porter PL, Bradley RK. RNA components of the spliceosome regulate tissue- and cancer-specific alternative splicing. BioRxiv. 2018;29:1591–604. https://doi.org/10.1101/326983.
Frangogiannis NG. Regulation of the inflammatory response in cardiac repair. Circ Res. 2012;110(1):159–73. https://doi.org/10.1161/CIRCRESAHA.111.243162.
Frangogiannis NG. Pathophysiology of myocardial infarction. Compr Physiol. 2015;5(4):1841–75. https://doi.org/10.1002/cphy.c150006.
Fu SP, He SY, Bin X, Hu CJ, Lu S-F, Shen WX, Huang Y, et al. Acupuncture promotes angiogenesis after myocardial ischemia through H3K9 acetylation regulation at VEGF gene. PLoS ONE. 2014;9(4):e94604. https://doi.org/10.1371/journal.pone.0094604.
Gang H, Hai Y, Dhingra R, Gordon JW, Yurkova N, Aviv Y, Li H, et al. A novel hypoxia-inducible spliced variant of mitochondrial death gene Bnip3 promotes survival of ventricular myocytes. Circ Res. 2011;108(9):1084–92. https://doi.org/10.1161/CIRCRESAHA.110.238709.
Gao J, Zhang L, Wang Y, Bo Lu, Cui H, Weixing Fu, Wang H, Youhua Yu, Xiaochun Yu. Antiarrhythmic effect of acupuncture pretreatment in rats subjected to simulative global ischemia and reperfusion—involvement of adenylate cyclase, protein kinase A, and L-type Ca2+ channel. J Physiol Sci. 2008;58(6):389–96. https://doi.org/10.2170/physiolsci.RP007108.
Gast M, Rauch BH, Haghikia A, Nakagawa S, Haas J, Stroux A, Schmidt D, et al. Long noncoding RNA NEAT1 modulates immune cell functions and is suppressed in early onset myocardial infarction patients. Cardiovasc Res. 2019;115(13):1886–906. https://doi.org/10.1093/cvr/cvz085.
Guo LY, Yang F, Peng LJ, Li YB, Wang AP. CXCL2, a new critical factor and therapeutic target for cardiovascular diseases. Clin Exp Hypertens. 2020;42(5):428–37. https://doi.org/10.1080/10641963.2019.1693585.
Haskill S, Peace A, Morris J, Sporn SA, Anisowicz A, Lee SW, Smith T, Martin G, Ralph P, Sager R. Identification of three related human GRO genes encoding cytokine functions. Proc Natl Acad Sci USA. 1990;87(19):7732–6. https://doi.org/10.1073/pnas.87.19.7732.
He J, Yao J, Sheng H, Zhu J. Involvement of the dual-specificity tyrosine phosphorylation-regulated kinase 1A-alternative splicing factor-calcium/calmodulin-dependent protein kinase IIδ signaling pathway in myocardial infarction-induced heart failure of rats. J Cardiac Fail. 2015;21(9):751–60. https://doi.org/10.1016/j.cardfail.2015.05.015.
Hoogenhof MMG, Den V, Beqqali A, Amin AS, Van Der Made I, Aufiero S, Khan MAF, Schumacher CA, et al. RBM20 mutations induce an arrhythmogenic dilated cardiomyopathy related to disturbed calcium handling. Circulation. 2018;138(13):1330–42. https://doi.org/10.1161/CIRCULATIONAHA.117.031947.
Hu C, Tkebuchava T, Dayi Hu. Managing acutemyocardial infarction in China. Eur Heart J. 2019;40(15):1179–81. https://doi.org/10.1093/eurheartj/ehz182.
Huang Y, Sheng Feng Lu, Chen Jun Hu, Shu Ping Fu, Shen WX, Liu WX, Li Q, et al. Electro-acupuncture at neiguan pretreatment alters genome-wide gene expressions and protects rat myocardium against ischemia-reperfusion. Molecules. 2014;19(10):16158–78. https://doi.org/10.3390/molecules191016158.
Jung M, Yonggang M, Rugmani PI, DeLeon-Pennell KY, Andriy Y, Garrett MR, Lindsey ML. IL-10 improves cardiac remodeling after myocardial infarction by stimulating M2 macrophage polarization and fibroblast activation. Basic Res Cardiol. 2017;112(3):33. https://doi.org/10.1007/s00395-017-0622-5.
Kong SW, Hu YW, Ho JWK, Sadakatsu I, Sean P, Ranjit J, Hall JL, et al. Heart failure-associated changes in RNA splicing of sarcomere genes. Circ Cardiovasc Genet. 2010;3(2):138–46. https://doi.org/10.1161/CIRCGENETICS.109.904698.
Lara-Pezzi E, Gómez-Salinero J, Gatto A, García-Pavía P. The alternative heart: impact of alternative splicing in heart disease. J Cardiovasc Transl Res. 2013;6(6):945–55. https://doi.org/10.1007/s12265-013-9482-z.
Lee JY, Hall JA, Kroehling L, Wu L, Najar T, Nguyen HH, Lin WY, Yeung ST, Silva HM, Li D, Hine A, Loke P, Hudesman D, Martin JC, Kenigsberg E, Merad M, Khanna KM, Littman DR. Serum amyloid a proteins induce pathogenic Th17 cells and promote inflammatory disease. Cell. 2020;183(7):2036–9. https://doi.org/10.1016/j.cell.2020.12.008
Lindsey ML, Bolli R, Canty JM, Xiao Jun Du, Frangogiannis NG, Frantz S, Gourdie RG, et al. Guidelines for experimental models of myocardial ischemia and infarction. Am J Physiol Heart Circ Physiol. 2018;314(4):H812–38. https://doi.org/10.1152/ajpheart.00335.2017.
Liu Y, Morley M, Brandimarto J, Sridhar Hannenhalli YuHu, Ashley EA, Wilson Tang WH, et al. RNA-Seq identifies novel myocardial gene expression signatures of heart failure. Genomics. 2015;105(2):83–9. https://doi.org/10.1016/j.ygeno.2014.12.002.
Liu Y, Meng HY, Khurwolah MR, Liu JB, Tang H, Aa N, Yang ZJ. Acupuncture therapy for the treatment of stable angina pectoris: an updated meta-analysis of randomized controlled trials. Complement Ther Clin Pract. 2019;34:247–53. https://doi.org/10.1016/j.ctcp.2018.12.012.
Luo H, Song S, Chen Y, Xu M, Sun L, Meng G, Zhang W. Inhibitor 1 of protein phosphatase 1 regulates Ca2+/calmodulin-dependent protein kinase II to alleviate oxidative stress in hypoxia-reoxygenation injury of cardiomyocytes. Oxid Med Cell Longev. 2019. https://doi.org/10.1155/2019/2193019.
Masaki M, Izumi M, Oshima Y, Nakaoka Y, Kuroda T, Kimura R, Sugiyama S, et al. Smad1 protects cardiomyocytes from ischemia-reperfusion injury. Circulation. 2005;111(21):2752–9. https://doi.org/10.1161/CIRCULATIONAHA.104.490946.
Müller T, Boileau E, Talyan S, Kehr D, Varadi K, Busch M, Most P, Krijgsveld J, Dieterich C. Updated and enhanced pig cardiac transcriptome based on long-read RNA sequencing and proteomics. J Mol Cell Cardiol. 2021;150:23–31. https://doi.org/10.1016/j.yjmcc.2020.10.005.
Pulkkinen HH, Kiema M, Lappalainen JP, Toropainen A, Beter M, Tirronen A, Holappa L, Niskanen H, Kaikkonen MU, Laakkonen JP. BMP6/TAZ-Hippo signaling modulates angiogenesis and endothelial cell response to VEGF. Angiogenesis. 2020. https://doi.org/10.1007/s10456-020-09748-4.
Ruiz-villalba A, Romero JP, Hernandez SC. Single-cell RNA-Seq analysis reveals a crucial role for collagen triple helix repeat containing 1 (CTHRC1) cardiac fibroblasts after myocardial infarction. Circulation. 2020;2020(142):1831–47.
Suzuki A, Suzuki M, Mizushima-Sugano J, Frith MC, Makałowski W, Kohno T, Sugano S, Tsuchihara K, Suzuki Y. Sequencing and phasing cancer mutations in lung cancers using a long-read portable sequencer. DNA Res. 2017;24(6):585–96. https://doi.org/10.1093/dnares/dsx027.
Tang HM, Haisong Ju, Zhao S, LaDuke C, Hahn S, Glick J, Carey C, Friedrichs GS. Translational assessment of cardiac contractility by echocardiography in the telemetered rat. J Pharmacol Toxicol Methods. 2016;77:24–32. https://doi.org/10.1016/j.vascn.2015.09.005.
Tarzami ST, Cheng R, Miao W, Kitsis RN, Berman JW. Chemokine expression in myocardial ischemia: MIP-2 dependent MCP-1 expression protects cardiomyocytes from cell death. J Mol Cell Cardiol. 2002;34(2):209–21. https://doi.org/10.1006/jmcc.2001.1503.
Urbanski LM, Leclair N, Anczuków O. Alternative-splicing defects in cancer: splicing regulators and their downstream targets, guiding the way to novel cancer therapeutics. Wiley Interdiscip Rev RNA. 2018;9(4):1–36. https://doi.org/10.1002/wrna.1476.
Watanabe T, Kimura A, Kuroyanagi H. Alternative splicing regulator RBM20 and cardiomyopathy. Front Mol Biosci. 2018;5:1–11. https://doi.org/10.3389/fmolb.2018.00105.
Williams AL, Khadka V, Tang M, Avelar A, Schunke KJ, Menor M, Shohet RV. HIF1 mediates a switch in pyruvate kinase isoforms after myocardial infarction. Physiol Genomics. 2018;50(7):479–94. https://doi.org/10.1152/physiolgenomics.00130.2017.
Wu X, Kun W, Shuai C, Shengbing W, Guoqi Z, Meiqi Z. Electroacupuncture ameliorates acute myocardial ischemia: a potential role of the locus coeruleus. Evid Based Complement Alternat Med. 2020. https://doi.org/10.1155/2020/4298657.
Xiao Y, Chen W, Zhong Z, Ding L, Bai H, Chen H, Zhang H, Yihuang Gu, Shengfeng Lu. Electroacupuncture preconditioning attenuates myocardial ischemia-reperfusion injury by inhibiting mitophagy mediated by the MTORC1-ULK1-FUNDC1 Pathway. Biomed Pharmacother. 2020;127: 110148. https://doi.org/10.1016/j.biopha.2020.110148.
Yang Q, Mao H, Chen X, Zhang Y, Zhang X, Liu Z, Jiang G, Huang W. Neiguan (PC6)-based acupuncture pretreatment for myocardial ischemia reperfusion injury: a protocol for preclinical systematic review and meta-analysis. Medicine. 2020;99(28): e20792. https://doi.org/10.1097/MD.0000000000020792.
Zeng Q, He H, Xian BW, Yu QZ, Hong XL, Zhi PT, Shang FH, Guo ZH. Electroacupuncture preconditioning improves myocardial infarction injury via enhancing AMPK-dependent autophagy in rats. BioMed Res Int. 2018. https://doi.org/10.1155/2018/1238175.
Zhai TY, Bao HC, Yang Z, Xin YX, Lei Z, Xin L, Xiao SZ, et al. Exosomes released from CaSR-stimulated PMNs reduce ischaemia/reperfusion injury. Oxid Med Cell Longev. 2021. https://doi.org/10.1155/2021/3010548.
Zhang HR, Tao JL, Bai H, Yang EM, Zhong ZH, Liu XT, Yi Huang Gu, Sheng Feng Lu. Changes in the serum metabolome of acute myocardial ischemia rat pretreatment with electroacupuncture. Am J Chin Med. 2019;47(5):1025–41. https://doi.org/10.1142/S0192415X19500526.
Zhang T, Yang W-X, Wang Y-L, Yuan J, Qian Yi, Sun Q-M, Mei-ling Yu, Shu-ping Fu, Bin Xu, Sheng-feng Lu. Electroacupuncture preconditioning attenuates acute myocardial ischemia injury through inhibiting NLRP3 inflammasome activation in mice. Life Sci. 2020;248: 117451. https://doi.org/10.1016/j.lfs.2020.117451.
Xiujuan Z, Li G, Jiang H, Li L, Ma J, Li H, Chen J. Full-length transcriptome analysis of litopenaeus vannamei reveals transcript variants involved in the innate immune system. Fish Shellfish Immunol. 2019;87:346–59. https://doi.org/10.1016/j.fsi.2019.01.023.
Zhao J, Chang L, Xianen Gu, Liu J, Sun B, Wei Xi. Systematic profiling of alternative splicing signature reveals prognostic predictor for prostate cancer. Cancer Sci. 2020;111(8):3020–31. https://doi.org/10.1111/cas.14525.
Zhao L, Li D, Zheng H, Chang X, Cui J, Wang R, Shi J, et al. Acupuncture as adjunctive therapy for chronic stable angina: a randomized clinical trial. JAMA Intern Med. 2019;179(10):1388–97. https://doi.org/10.1001/jamainternmed.2019.2407.
Zhao S. Alternative splicing, RNA-Seq and drug discovery. Drug Discov Today. 2019;24(6):1258–67. https://doi.org/10.1016/j.drudis.2019.03.030.
Zhu LY, Zhu YR, Dai DJ, Wang X, Jin HC. Epigenetic regulation of alternative splicing. Am J Cancer Res. 2018;8(12):2346–58.
Zhuo LA, Wen YT, Wang Y, Liang ZF, Gang Wu, Nong MD, Miao L. LncRNA SNHG8 is identified as a key regulator of acute myocardial infarction by RNA-Seq analysis. Lipids Health Dis. 2019;18(1):1–17. https://doi.org/10.1186/s12944-019-1142-0.
Zinkernagel AS, Johnson RS, Nizet V. Hypoxia inducible factor (HIF) function in innate immunity and infection. J Mol Med. 2007;85(12):1339–46. https://doi.org/10.1007/s00109-007-0282-2.
This study was supported by the National Key R&D Program of China (No. 2019YFC1709001), the National Natural Science Foundation of China (No. 81774434), the National Natural Science Foundation of China (No. 81704187), the Department of Science and Technology of Sichuan Province (No. 2019YJ0587), the Department of Science and Technology of Sichuan Province (No. 2018JY0482).
Ethics approval and consent to participate
This study was approved by the Committee on Ethical Use of Animals of Chengdu University of Traditional Chinese Medicine.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yuan, J., Wang, JM., Li, ZW. et al. Full-length transcriptome analysis reveals the mechanism of acupuncture at PC6 improves cardiac function in myocardial ischemia model. Chin Med 16, 55 (2021). https://doi.org/10.1186/s13020-021-00465-8