- Open Access
Targeted myocardial gene expression in failing hearts by RNA sequencing
Journal of Translational Medicine volume 14, Article number: 327 (2016)
Myocardial recovery with left ventricular assist device (LVAD) therapy is highly variable and difficult to predict. Next generation ribonucleic acid (RNA) sequencing is an innovative, rapid, and quantitative approach to gene expression profiling in small amounts of tissue. Our primary goal was to identify baseline transcriptional profiles in non-ischemic cardiomyopathies that predict myocardial recovery in response to LVAD therapy. We also sought to verify transcriptional differences between failing and non-failing human hearts.
RNA was isolated from failing (n = 16) and non-failing (n = 8) human hearts. RNA from each patient was reverse transcribed and quantitatively sequenced on the personal genome machine (PGM) sequencer (Ion torrent) for 95 heart failure candidate genes. Coverage analysis as well as mapping the reads and alignment was done using the Ion Torrent Browser Suite™. Differential expression analyses were conducted by empirical analysis of digital gene expression data in R (edgeR) to identify differential expressed genes between failing and non-failing groups, and between responder and non-responder groups respectively. Targeted cardiac gene messenger RNA (mRNA) expression was analyzed in proportion to the total number of reads. Gene expression profiles from the PGM sequencer were validated by performing RNA sequencing (RNAseq) with the Illumina Hiseq2500 sequencing system.
The failing sample population was 75% male with an average age of 50 and a left ventricular ejection fraction (LVEF) of 16%. Myosin light chain kinase (MYLK) and interleukin (IL)-6 genes expression were significantly higher in LVAD responders compared to non-responders. Thirty-six cardiac genes were expressed differentially between failing and non-failing hearts (23 decreased, 13 elevated). MYLK, Beta-1 adrenergic receptor (ADRB1) and myosin heavy chain (MYH)-6 expression were among those significantly decreased in failing hearts compared to non-failing hearts. Natriuretic peptide B (NPPB) and IL-6 were significantly elevated. Targeted gene expression profiles obtained from the Ion torrent PGM sequencer were consistent with those obtained from Illumina HiSeq2500 sequencing system.
Heart failure is associated with a network of transcriptional changes involving contractile proteins, metabolism, adrenergic receptors, protein phosphorylation, and signaling factors. Myocardial MYLK and IL-6 expression are positively correlated with ejection fraction (EF) response to LVAD placement. Targeted RNA sequencing of myocardial gene expression can be utilized to predict responders to LVAD therapy and to better characterize transcriptional changes in human heart failure.
Left ventricular assist devices are commonly being utilized to treat patients with advanced heart failure . Mortality due to infection, bleeding, stroke, right heart failure and device malfunction continues as a significant problem in this population but increasingly patients are receiving these devices as destination therapy . Cardiac transplantation is not an option for the majority of these patients due to limitations in donor supply . Myocardial recovery with LVAD therapy to the degree that allows device explant remains relatively rare for these patients [1, 3]. LVAD therapy may be associated with myocardial atrophy, changes in transcriptional profiles, and persistent abnormalities in protein phosphorylation [4, 5].
There is however increasing evidence that a subset of patients with LVAD therapy may remodel and recover sufficient myocyte function to allow device explant. Utilizing ventricular unloading, myocardial conditioning and guideline directed heart failure medical therapy it may feasible to explant a significant number of patients with LVADs . Myocardial recovery, while unusual in ischemic cardiomyopathies, may occur in 20–25% of non-ischemic cardiomyopathy patients .
Heart failure is characterized on a molecular level by transcriptional changes involving adrenergic signaling mechanisms, calcium handling, contractile proteins, inflammation, and metabolism . Next-generation sequencing allows for gene expression quantification, the study of epigenetic effects, identification of signal nucleotide polymorphisms and measurement of transcriptional variations. Understanding the transcriptome is essential to understanding disease development, progression, and response to therapeutics. Transcriptional profiling is currently being used to guide treatment in many forms of cancer and to diagnose cardiac rejection. It is unknown if gene expression profiles can predict response to LVAD therapy. Our primary goal was to identify transcriptional profiles in non-ischemic cardiomyopathy patients that predict myocardial recovery. We also sought to verify transcriptional differences between failing and non-failing human hearts utilizing targeted non-optical sequencing.
This study was approved by the IRB at the University of Nebraska Medical Center. Failing tissue was obtained from non-ischemic cardiomyopathy patients (n = 16) who met standard criteria for LVAD placement either as destination therapy or as bridge to transplant. Tissue was obtained at the time of LVAD placement and stored in RNA later prior to freezing at −80 °C. Non-failing tissue was obtained by endomyocardial biopsy in transplant patients with normal cardiac function during routine surveillance (n = 8). Cardiac function was evaluated by echocardiography.
Isolation of RNA
RNA was isolated from approximately 10–20 mg of cardiac tissue using the RNeasy® Mini Kit (Qiagen®). Quantitative measurement of the RNA was taken using Qubit® RNA assay kit by Qubit 2.0 fluorometer (Invitrogen, Life Technologies).
Constructing RNA Library
Ion AmpliSeq™ RNA Library Kit and Custom Panels were utilized in the construction of RNA library for this investigation, and all experiments were carried out in accordance with the instructions from the manufacturer (Ion Torrent, Life Technologies). 10 ng of RNA, isolated from the cardiac tissue, is reverse transcribed to synthesize complementary deoxyribonucleic acid (cDNA) using the Ion AmpliSeq™ RNA RT Module and Applied Biosystems thermal cycler. Next target sequences were amplified using Ion AmpliSeq™ RNA Custom Panels and Library Kit (life technologies). Our custom panel was designed to target 95 cardiac genes we previously have associated with reverse remodeling in heart failure utilizing the Ion AmpliSeq™ Designer . After amplification of target sequences, primer sequences were partially digested using FuPa reagent from Ion AmpliSeq™ RNA Library Kit. Ion AmpliSeq™ adaptors were then ligated to the targeted deoxyribonucleic acid (DNA) and purification of those DNA fragments followed. Target DNA sequences were purified in a two-round purification process with Dynabeads® Magnetic Beads where surplus primer sequences and high-molecular weight DNA are isolated and discarded from the solution. Following purification of the library, the library is amplified using Ion AmpliSeq™ RNA Library Kit and further purified using single-round purification with Dynabeads® Magnetic Beads. After amplification and purification of the library, Qubit® 2.0 fluorometer was used with Qubit® dsDNA HS Assay Kit to quantitatively measure the DNA library.
Preparation of targeted DNA template
Amplified stock library was diluted for appropriate Ion library preparation. Ion One Touch™ 2 System was utilized to amplify individual diluted libraries via emulsion polymerase chain reaction (PCR) on ion sphere particles (ISPs). The template-positive ISPs were then enriched using Ion One Touch™ Enrichment System.
Sequencing was performed for all patients on the Ion Torrent PGM utilizing the Ion 316 chip. Coverage analysis as well as mapping the reads and alignment was done using the Ion Torrent Browser Suite™. Gene expression profiles from the PGM sequencer were validated with a subsequent population by performing RNAseq with the Illumina Hiseq2500 sequencing system. RNAs from 7 non failing and 3 failing samples were sequenced using the Illumina HiSeq2500 sequencing system (University of Nebraska Medical Center, DNA core facility). All paired end sequencing adaptors were trimmed using fqtrim and any reads shorter than 36 bases were discarded. A quality filter using fqtrim −q 5 was then applied to the reads. Paired end reads that passing those filters were mapped to Homo sapiens (release 37) reference sequence (GRCh37/hg19); using BowTie2 and TopHat2. The raw read counts per gene were generated using the Rsubread package.
For the LVAD responder analysis the failing heart patients were divided into two groups: responders (∆LVEF ≥ 20%) and non-responders (∆LVEF < 20%). The expression of each gene on the cardiac panel was assessed in proportion to the total reads on the chip. The differential expression analysis was conducted using edgeR package in bioconductor developed by Robinson et al. [9, 10]. Differential gene expression between the non-failing heart patients (n = 8) and failing heart patients (n = 16) was analyzed in a similar fashion. The gene Calreticulin 3 (CALR3) was removed from the data due to 21 samples had zero counts for this gene. Spearman’s correlation coefficients were calculated to evaluate inter-platform consistency for raw counts and normalized counts per million (cpm) across two platforms (Ion Torrent and Illumina). The normalization factors for the data were estimated by the trimmed mean of the M-values normalization method in the edgeR package and used to adjust for varying sequencing depths and potentially other technical effects across samples.
To further evaluate the biological pathways, log fold change utilizing Ingenuity Pathway Analysis (IPA) was carried out on the targeted panel genes between failing and non-failing samples.
We analyzed a total 34 of human heart samples, 16 failing, and 8 non-failing. Patient clinical characteristics are shown on Table 1. The failing sample population was 75% male with an average age 50. The failing patients had a New York Heart Association (NYHA) functional classification of III-IV and diabetes mellitus was common. EFs were measured before and after 3 months LVAD implantation. Change in ejection fraction (ΔEF) greater than 20 signifies responder after LVAD placement. By these criteria, 10 patients were non responders and 6 patients were responders. The non- failing patients had EFs ranging from 60 to 65% whereas the EFs of LVAD patients were 15–25% at the time of implant. All patients in this study were recipients of the Heartmate II LVAD at the University of Nebraska Medical Center. Our center implanted 177 LVADs between June of 2006 and June of 2016. Most of these implants (55%) were for critical cardiogenic shock or a progressive declining clinical status despite inotropic support. The remainder were stable on inotropic therapy or with advanced resting symptoms. The overall clinical outcomes of LVAD therapy are shown in Fig. 1. Overall survival with LVAD as bridge to transplant or destination therapy at one year is approximately 80%. Very few patients historically have been explanted for recovery.
Responder vs non-responder RNA-Seq
Our study showed that MYLK (p = 0.005) and IL 6 (p < 0.0002) expression were significantly higher in LVAD responder patients than non-responder patients, Fig. 2. MYLK gene or myosin light chain kinase activates and regulates myosin light chain in the heart by phosphorylation. MYLK is necessary for normal cardiac function. Down regulation of MYLK expression increases the risk of heart failure [11–14]. Levels of IL-6 are found considerably higher in myocardium and serum of patients receiving LVAD and with advanced heart failure [15–17]. Further analysis of an additional 5 patients pre and post LVAD who did not recover cardiac function and required transplant shows that LVAD therapy decreased myocardial IL-6 (log fold change −2, p = 0.004) but did not change MYLK (p = 0.6).
Failing vs non-failing RNA-Seq results
Here we tried to understand the changes of cardiac genes expression during heart failure in association with contractile dysfunction. We evaluated mRNA expression of 95 cardiac gene using RNA-Sequence technology. It was found that 36 genes were differentially expressed when we compared failing versus non failing hearts (p < 0.05). Twenty-three genes were significantly decreased in failing hearts including MYLK, ADRB1 and MYH6, Fig. 3. In contrast 13 genes were up-regulated including NPPB and IL-6, Fig. 4. A complete list of genes studied, p values and false discovery rates  are listed in Appendix.
Functional and pathway analysis in IPA
To further analyze the biological pathways, log fold change utilizing IPA was carried out on the targeted panel genes between failing and non-failing hearts. The differentially expressed genes involved in cardiomyopathy and muscle formation were plotted. The top canonical pathways, molecular functions, upstream regulators and toxic effects are shown in Table 2. Figure 5 illustrates the network for genes that were differentially expressed, their molecular function and their relationship to different diseases including arrhythmia, hypertrophy of the heart, familial heart disease and necrosis of cardiac muscles. As shown in Fig. 5, eight genes of the targeted panel were involved in hypertrophy of the heart. Four genes clearly lead to activation of this pathway, up-regulation of IL-6 (an activator) and down regulation of SMAD4, PRKCB and SLC25A4 (inhibitors). These molecular changes collectively augment the hypertrophy signal of the heart. Six genes of the targeted panel were involved in the necrosis of cardiac muscle with a net z score of −2.115. Up-regulation of HK2 and IL6 (known inhibitors) and down regulation of ADRB1, PRKCB and TNF (known activators). The net effect of these changes are predicted to inhibit necrosis of the heart. Finally, nine molecules with a net of “not predicted effect” were involved in the familial cardiovascular disease and arrhythmia that could be furthered studied.
To validate the results obtained from Ion torrent PGM sequencer by RNAseq analysis, Illumina HiSeq sequencing was utilized. Ten patients were sequenced for selected 95 genes using RNAseq method on both sequencing platforms (Illumina HiSeq and Ion Torrent). Illumina data were processed to obtain raw read counts for genes by BowTie, TopHat, and Rsubread packages respectively [19–21]. To evaluate inter-platform consistency for gene expression measures, we computed Spearman’s correlation between raw counts and normalized cpm from two platforms for each individual, and the average spearman’s correlation across 10 individuals. The results demonstrated high inter-platform consistency for gene expression measures across two platforms. (Spearman R ≥ 0.83 for each patient sample, and average Spearman R = 0.84 using raw counts. Spearman R ≥ 0.84 for each patient sample and average Spearman R = 0.86 using normalized cpm).
The main findings of this study are that MYLK2 is diminished in heart failure, predicts response to LVAD therapy, and does not improve in advanced heart failure patients post LVAD who require transplantation. Our results further characterize transcriptional changes between failing and non-failing human hearts and identify gene expression profiles associated with myocardial recovery on LVAD therapy. Heart failure is a heterogeneous disease process. Etiology, disease duration, cardiac remodeling, medical therapy, and molecular signatures have all been associated with myocardial recovery in non-ischemic cardiomyopathies. Recent studies suggest that up to 20 percent of patients may recover their cardiac function within their first year post LVAD . Molecular phenotyping is increasingly being utilized as a clinical tool to risk stratify patients and guide clinical decision making. Transcriptional profiles are being utilized to screen heart transplant patient for rejection, identify patients with coronary disease, and guide treatment decisions for cancer patients. It may be feasible utilizing clinical characteristics and molecular markers to risk stratify patients for early transplant versus more prolonged aggressive attempts to recover myocardial function with medical and device therapies. Previous studies have identified myocardial microRNA (miRNA) signatures as potential biomarkers of recovery. Our results suggest that mRNA patterns may also be helpful in risk stratifying patients undergoing LVAD therapy.
Baseline MYLK and IL-6 levels were both abnormal in failing myocardium and associated with response to LVAD therapy. Myosin is comprised of two heavy chain and two pairs of light chains. Myosin light chain 2 (MLC-2) is a key regulatory component of myosin. Phosphorylation of MLC-2 is known to regulate both calcium sensitivity as well as sarcomere assembly [11, 23, 24]. MLC-2 can be phosphorylated by MYLK or protein kinase C and is dephosphorylated by light chain phosphatase. MLC-2 phosphorylation is known to be decreased in the end-stage failing human heart . Our results would suggest that down-regulation of MYLK expression maybe contributing to this process. The association of higher levels of MYLK with improved myocardial recovery on LVAD therapy is consistent with prior evidence supporting this gene as an adaptive mechanism which facilities sarcomere organization in response to hypertrophic stimuli. Persistent abnormalities in MYLK post LVAD therapy in patients who did not recover are consistent with ongoing stress and injury in this population.
IL-6 is known to be elevated systemically in heart failure and levels have been associated with severity of disease . Elevated expression of myocardial IL-6 in advanced heart failure could be contributing to this process. IL-6 can be produced by leukocytes, endothelial cells or vascular smooth muscle cells. Levels have been correlated with hemodynamic and neurohormonal variables. IL-6 is capable of producing myocardial dysfunction, vascular dilation, and muscle wasting. It is paradoxical that higher IL-6 levels at baseline are associated with enhanced recovery with LVAD therapy. IL-6 appears to be an important mediator in ventricular hypertrophy and does increase nitric oxide production . These mechanisms may play an important role in adapting to heart failure. IL-6 is extremely inducible in response to IL-1, TNFα, viral infection and angiotensin II peptide. It acts through two distinct mechanisms. The first one is a classic membrane receptor initiated pathway and the second one is trans-signaling pathway. It has been shown that IL-6 induces and activates signal transducer and activator of transcription 3 (STAT3) gene by engaging the suppressor of cytokine signaling 3 (SOCS3) gene. Multiple studies showed that the activation of STAT3 promotes cardiomyocyte survival and hypertrophy, as well as cardiac angiogenesis . Our results suggest LVAD therapy decreases myocardial IL-6 expression substantially even in patients who do not recover. Hence, it is unlikely to be contributing to the failure of myocardial recovery in patients post LVAD.
Heart failure is a complex molecular disease process involving more than IL-6 and MYLK. Numerous transcriptional changes have been described involving a complex gene network including contractile proteins, calcium handling, metabolism, signal transduction, apoptosis and inflammation [8, 29–31]. High throughput RNA sequencing now allows for the rapid quantification of gene expression in small amounts of tissue either globally or in a targeted fashion.
Most next generation sequencing technologies currently on the market use image capture, fluorescent detection or optical registration of some kind to capture sequencing data, which is then, converted to digital information. Our study utilized direct sequencing on a disposable semiconductor chip in a targeted fashion. This technology is relatively rapid, inexpensive and allows for quantification of mRNA over a broad range of gene expression levels. Between failing and non-failing human hearts we identify 36 genes that are differentially expressed. Other studies have previously identified these genes in association with heart failure and sequencing by optical and non-optical methodologies produced similar results. This suggests that RNA sequencing might provide a valuable tool to further characterize the molecular mechanisms of heart failure.
This study has several limitations. In regards to recovery, it is a single center retrospective analysis. There are several potential confounders that were not evaluated including genetic etiologies which could have impacted recovery independent of baseline gene expression. The study evaluated the transcriptome of only a limited number of genes. There are almost certainly other genes whose expression may be predictive of recovery or important mechanistically in heart failure. Finally, the transplant patients utilized as controls within this study may have different gene expression secondary to other factors such as immunosuppressive therapy independent of contractile dysfunction. For this reason we targeted candidate genes that previously have been identified in association with heart failure or myocardial remodeling.
While cardiac transplantation remains the gold standard for treating advanced heart failure, it is not an effective option for the majority of patients due to limitations in donor supply. The number of patients receiving left ventricular assist devices is growing exponentially. Less than 1% of patients are currently having their devices explanted for recovery. The rapid evolution in sequencing technologies may allow us to better classify LVAD patients who may recover. Prospective clinical trials to identify biomarkers of recovery are needed to help decision making in this population of patients. RNA signatures and clinical characteristics may ultimately help identify a subset of patients who can recover cardiac function.
MYLK2 is diminished in heart failure, predicts response to LVAD therapy, and does not improve in advanced heart failure patients post LVAD who require transplantation. Targeted next generation RNA sequencing allows rapid quantitative characterization of the transcriptional changes that occur in heart failure. Changes in this gene network involves contractile proteins, calcium handling, metabolism, signal transduction, apoptosis and inflammation. Pathway analysis suggests that these changes ultimately contribute to cardiac hypertrophy, cardiac necrosis, cell death, increased heart failure risk and cardiac fibrosis.
left ventricular assist device
personal genome machine
(empirical analysis of digital gene expression data in R)
left ventricular ejection fraction
myosin light chain kinase
Beta-1 adrenergic receptor
myosin heavy chain 6
natriuretic peptide B
polymerase chain reaction
ion sphere particles
change in left ventricular ejection fraction
counts per million
Ingenuity Pathway Analysis
New York Heart Association
change in ejection fraction
myosin light chain 2
signal transducer and activator of transcription 3
suppressor of cytokine signaling 3
Kirklin JK, Naftel DC, Pagani FD, Kormos RL, Stevenson LW, Blume ED, Miller MA, Timothy Baldwin J, Young JB. Sixth INTERMACS annual report: a 10,000-patient database. J Heart Lung Transplant. 2014;33:555–64.
Lund LH, Edwards LB, Kucheryavaya AY, Benden C, Christie JD, Dipchand AI, Dobbels F, Goldfarb SB, Levvey BJ, Meiser B, Yusen RD, Stehlik J, Transplantation ISoHaL. The registry of the International Society for Heart and Lung Transplantation: thirty-first official adult heart transplant report—2014; focus theme: retransplantation. J Heart Lung Transplant. 2014;33:996–1008.
Mancini DM, Beniaminovitz A, Levin H, Catanese K, Flannery M, DiTullio M, Savin S, Cordisco ME, Rose E, Oz M. Low incidence of myocardial recovery after left ventricular assist device implantation in patients with chronic heart failure. Circulation. 1998;98:2383–9.
Ambardekar AV, Walker JS, Walker LA, Cleveland JC, Lowes BD, Buttrick PM. Incomplete recovery of myocyte contractile function despite improvement of myocardial architecture with left ventricular assist device support. Circ Heart Fail. 2011;4:425–32.
Lowes BD, Zolty R, Shakar SF, Brieke A, Gray N, Reed M, Calalb M, Minobe W, Lindenfeld J, Wolfel EE, Geraci M, Bristow MR, Cleveland J. Assist devices fail to reverse patterns of fetal gene expression despite beta-blockers. J Heart Lung Transplant. 2007;26:1170–6.
Frazier OH, Baldwin AC, Demirozu ZT, Segura AM, Hernandez R, Taegtmeyer H, Mallidi H, Cohn WE. Ventricular reconditioning and pump explantation in patients supported by continuous-flow left ventricular assist devices. J Heart Lung Transplant. 2014;34:766–72.
Guglin M, Miller L. Myocardial recovery with left ventricular assist devices. Curr Treat Options Cardiovasc Med. 2012;14:370–83.
Kao D, Lowes B, Gilbert E, Minobe W, Epperson LE, Meyer L, Ferguson D, Volkman K, Zolty R, Borg D, Quaife R, Bristow M. Therapeutic molecular phenotype of β-blocker associated reverse-remodeling in nonischemic dilated cardiomyopathy. Circ Cardiovasc Genet. 2015;8:270–83.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40:4288–97.
Seguchi O, Takashima S, Yamazaki S, Asakura M, Asano Y, Shintani Y, Wakeno M, Minamino T, Kondo H, Furukawa H, Nakamaru K, Naito A, Takahashi T, Ohtsuka T, Kawakami K, Isomura T, Kitamura S, Tomoike H, Mochizuki N, Kitakaze M. A cardiac myosin light chain kinase regulates sarcomere assembly in the vertebrate heart. J Clin Invest. 2007;117:2812–24.
Warren SA, Briggs LE, Zeng H, Chuang J, Chang EI, Terada R, Li M, Swanson MS, Lecker SH, Willis MS, Spinale FG, Maupin-Furlowe J, McMullen JR, Moss RL, Kasahara H. Myosin light chain phosphorylation is critical for adaptation to cardiac stress. Circulation. 2012;126:2575–88.
Chan JY, Takeda M, Briggs LE, Graham ML, Lu JT, Horikoshi N, Weinberg EO, Aoki H, Sato N, Chien KR, Kasahara H. Identification of cardiac-specific myosin light chain kinase. Circ Res. 2008;102:571–80.
Toepfer C, Caorsi V, Kampourakis T, Sikkel MB, West TG, Leung MC, Al-Saud SA, MacLeod KT, Lyon AR, Marston SB, Sellers JR, Ferenczi MA. Myosin regulatory light chain (RLC) phosphorylation change as a modulator of cardiac muscle contraction in disease. J Biol Chem. 2013;288:13446–54.
Caruso R, Verde A, Campolo J, Milazzo F, Russo C, Boroni C, Parolini M, Trunfio S, Paino R, Martinelli L, Frigerio M, Parodi O. Severity of oxidative stress and inflammatory activation in end-stage heart failure patients are unaltered after 1 month of left ventricular mechanical assistance. Cytokine. 2012;59:138–44.
Caselli C, D’Amico A, Ragusa R, Caruso R, Prescimone T, Cabiati M, Nonini S, Marraccini P, Del Ry S, Trivella MG, Parodi O, Giannessi D. IL-33/ST2 pathway and classical cytokines in end-stage heart failure patients submitted to left ventricular assist device support: a paradoxic role for inflammatory mediators? Mediators Inflamm. 2013;2013:498703.
Birks EJ, Latif N, Owen V, Bowles C, Felkin LE, Mullen AJ, Khaghani A, Barton PJ, Polak JM, Pepper JR, Banner NR, Yacoub MH. Quantitative myocardial cytokine expression and activation of the apoptotic pathway in patients who require left ventricular assist devices. Circulation. 2001;104:I233–40.
Benjamini Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41:e108.
Drakos SG, Wever-Pinzon O, Selzman CH, Gilbert EM, Alharethi R, Reid BB, Saidi A, Diakos NA, Stoker S, Davis ES, Movsesian M, Li DY, Stehlik J, Kfoury AG, Investigators UUCRP. Magnitude and time course of changes induced by continuous-flow left ventricular assist device unloading in chronic heart failure: insights into cardiac recovery. J Am Coll Cardiol. 2013;61:1985–94.
Aoki H, Sadoshima J, Izumo S. Myosin light chain kinase mediates sarcomere organization during cardiac hypertrophy in vitro. Nat Med. 2000;6:183–8.
van der Velden J, Papp Z, Boontje NM, Zaremba R, de Jong JW, Janssen PM, Hasenfuss G, Stienen GJ. The effect of myosin light chain 2 dephosphorylation on Ca2+ - sensitivity of force is enhanced in failing human hearts. Cardiovasc Res. 2003;57:505–14.
van der Velden J, Papp Z, Zaremba R, Boontje N, de Jong J, Owen V, Burton P, Goldmann P, Jaquet K, Stienen G. Increased Ca2+ -sensitivity of the contractile apparatus in end-stage human heart failure results from altered phosphorylation of contractile proteins. Cardiovasc Res. 2003;57:37–47.
Tsutamoto T, Hisanaga T, Wada A, Maeda K, Ohnishi M, Fukai D, Mabuchi N, Sawaki M, Kinoshita M. Interleukin-6 spillover in the peripheral circulation increases with the severity of heart failure, and the high plasma level of interleukin-6 is an important prognostic predictor in patients with congestive heart failure. J Am Coll Cardiol. 1998;31:391–8.
Hirota H, Yoshida K, Kishimoto T, Taga T. Continuous activation of gp130, a signal-transducing receptor component for interleukin 6-related cytokines, causes myocardial hypertrophy in mice. Proc Natl Acad Sci USA. 1995;92:4862–6.
Hou T, Tieu BC, Ray S, Recinos Iii A, Cui R, Tilton RG, Brasier AR. Roles of IL-6-gp130 signaling in vascular inflammation. Curr Cardiol Rev. 2008;4:179–92.
Lowes BD, Minobe W, Abraham WT, Rizeq MN, Bohlmeyer TJ, Quaife RA, Roden RL, Dutcher DL, Robertson AD, Voelkel NF, Badesch DB, Groves BM, Gilbert EM, Bristow MR. Changes in gene expression in the intact human heart. Downregulation of alpha-myosin heavy chain in hypertrophied, failing ventricular myocardium. J Clin Invest. 1997;100:2315–24.
Lowes BD, Zolty R, Minobe WA, Robertson AD, Leach S, Hunter L, Bristow MR. Serial gene expression profiling in the intact human heart. J Heart Lung Transplant. 2006;25:579–88.
Lowes BD, Gilbert EM, Abraham WT, Minobe WA, Larrabee P, Ferguson D, Wolfel EE, Lindenfeld J, Tsvetkova T, Robertson AD, Quaife RA, Bristow MR. Myocardial gene expression in dilated cardiomyopathy treated with beta-blocking agents. N Engl J Med. 2002;346:1357–65.
BL, ER, JU, MM, DA and RZ were responsible for patient care and recruitment. KD, ER, JE and AM were responsible for laboratory analysis. JK was responsible for data analysis and FQ was responsible for statistical analysis. All authors read and approved the final manuscript.
This work was supported by the William Dodge Angle Endowment, NIH P20 HL101435, PI: Lowes, Brian, and the Free Radicals in Medicine Program at the University of Nebraska Medical Center. We also want to acknowledge Kathy Duncan for the clerical support she provided.
The authors declare that they have no competing interests.
Compliance with ethical guidelines
This study was approved by the IRB at the University of Nebraska Medical Center.
William Dodge Angle Endowment, P20 HL101435 (Brian Lowes PI).
See Table 3.
About this article
Cite this article
Dhar, K., Moulton, A.M., Rome, E. et al. Targeted myocardial gene expression in failing hearts by RNA sequencing. J Transl Med 14, 327 (2016) doi:10.1186/s12967-016-1083-6
- Left Ventricular Assist Device
- Myocardial Recovery
- Left Ventricular Assist Device Patient
- Advanced Heart Failure Patient
- Left Ventricular Assist Device Therapy