Purpose:Acupuncture exerts cardioprotective effect on myocardial I/R injury. In order to elucidate the potential mechanisms, RNA-seq by next generation sequencing was used to identify the rat genome-wide alterations by EA at Neiguan acupoint pretreatment in this study Methods: Adult male Sprague Dawley rats (250-300g) were divided into four groups: Sham operation(SO), I/R, electro-acupuncture at Neiguan pretreatment (EA) and electro-acupuncture at non-acupoint (NA) groups and myocardium mRNA profiles of rats in each group were generated by deep sequencing,using Illumina Hiseq 2000. The sequence reads that passed quality filters were analyzed at the transcript isoform level with two methods: Burrows–Wheeler Aligner (BWA) followed by ANOVA (ANOVA) and TopHat followed by Cufflinks. Results: Using an optimized data analysis workflow, we mapped about 20~43 million sequence reads per sample to the mouse genome (UCSC RN4) and identified 26190 22895 25160 22852 transcripts in the SO, I/R, EA and NA rats with TopHat and Cufflinks workflow respectively. Approximately 32% and 22% of the transcripts showed differential expression between the SO and I/R, I/R and EA , I/R and NA, respectively with a |log2 fold change |≥1. Conclusions: Our results for the first time generated genome-wide gene expression profiles both in the rat I/R model and in acupuncture treatment by high throughput sequencing. The optimized data analysis workflows reported here should provide a framework for acupuncture investigations of expression profiles.

Hspa1b, Mir3074, Kbtbd5, Fam46b, Hist1h1d, Tnfrsf12a, Fos, Atf3, Has1, Camp, Mir3556b, Vdac3, Trh, Hmox1, Cyr61, Pgf, Slc7a5, Klk12, Sphk1, Cblc
Nr4a1, Numbl, Mt1a, Flnc, Lcn2, Gal, Alox15, Akr1b8, Mt2A, Cnot3, LOC367975, Mir145, Klre1, Ky, Clecsf6, RGD1306750, Cybb, Ptplad2, Cyp2e1, Sh2d1a
Tas2r120, Epsti1, C6, F13a1 Bex4, Tfec, Msr1, Ptprc, Clec4a2, Cd69, Agr2, Igsf6, Lilra5, Klri1, Ubd, Ahr, Ccr2, Pkhd1l1, Mpeg1, Dock11, Actc1
Myl3, Cyc1, Cox6a1, Myh6, Tnni3, Cacna1c, Cox5b, Kcnmb4, Adora2a, Rhoc, Pla2g5, Myl9, Atp5d, Ndufb7, Ndufb8, Ndufa7, Atp5g2, Atp5g1, Atp5g3, Ndufb2
Ndufs7, Ndufs8, Cos6a1, Atp5l, Cav3, Bcar1, Col5a3, Col5a1, Vegfb, Lama5, Col6a2, Col1a2, Col1a1, Myh7b, Nppb, Tnfsf10, Xiap, IL1rap, Pik3r1, Chuk
Casp8, Casp3, AIF, Endo G, Bcl2a1d, Gja1, Bcl-2, Adora1, GAPDH