- Open Access
Long noncoding RNA dysregulation in ischemic heart failure
Journal of Translational Medicinevolume 14, Article number: 183 (2016)
Long noncoding RNAs (lncRNAs) are non-protein coding transcripts regulating a variety of physiological and pathological functions. However, their implication in heart failure is still largely unknown. The aim of this study is to identify and characterize lncRNAs deregulated in patients affected by ischemic heart failure.
LncRNAs were profiled and validated in left ventricle biopsies of 18 patients affected by non end-stage dilated ischemic cardiomyopathy and 17 matched controls. Further validations were performed in left ventricle samples derived from explanted hearts of end-stage heart failure patients and in a mouse model of cardiac hypertrophy, obtained by transverse aortic constriction. Peripheral blood mononuclear cells of heart failure patients were also analyzed. LncRNA distribution in the heart was assessed by in situ hybridization. Function of the deregulated lncRNA was explored analyzing the expression of the neighbor mRNAs and by gene ontology analysis of the correlating coding transcripts.
Fourteen lncRNAs were significantly modulated in non end-stage heart failure patients, identifying a heart failure lncRNA signature. Nine of these lncRNAs (CDKN2B-AS1/ANRIL, EGOT, H19, HOTAIR, LOC285194/TUSC7, RMRP, RNY5, SOX2-OT and SRA1) were also confirmed in end-stage failing hearts. Intriguingly, among the conserved lncRNAs, h19, rmrp and hotair were also induced in a mouse model of heart hypertrophy. CDKN2B-AS1/ANRIL, HOTAIR and LOC285194/TUSC7 showed similar modulation in peripheral blood mononuclear cells and heart tissue, suggesting a potential role as disease biomarkers. Interestingly, RMRP displayed a ubiquitous nuclear distribution, while H19 RNA was more abundant in blood vessels and was both cytoplasmic and nuclear. Gene ontology analysis of the mRNAs displaying a significant correlation in expression with heart failure lncRNAs identified numerous pathways and functions involved in heart failure progression.
These data strongly suggest lncRNA implication in the molecular mechanisms underpinning HF.
It is estimated that heart failure (HF) is the cardiovascular disease with the worse rate of morbidity, mortality, accounting for one of the highest health care costs in the western world. Indeed, the 1-year-survival rate of patients with end-stage HF is about 50 % .
The adverse left ventricle (LV) remodeling process, leading to the clinical syndrome of HF, involves several deregulated proteins and is characterized in the adult heart by the reactivation of fetal cardiac gene expression . This scenario of transcriptional control is also complicated by the addition of epigenetic mechanisms. The encyclopedia of DNA elements (ENCODE) project indicates that at least 80 % of the genome is functional and is transcribed both into protein coding RNAs (about 1.1–1.5 %) and in a much larger quantity of non-coding, regulatory RNAs, arbitrarily divided into long (lncRNAs, >200 nt), and short (<200 nt) ncRNAs [3, 4].
While dysregulation and role of short ncRNAs, in particular of miRNAs, has been extensively explored [5–7], the involvement of lncRNAs in specific physiological and pathological processes [5, 8], as well as in cardiovascular diseases [9–16], is still in its early stages of study.
In this study, in order to investigate the molecular mechanisms underpinning HF during its progression, we profiled the expression of 83 lncRNAs, known to be implicated in human diseases, in LV biopsies of non end-stage HF patients.
Patient selection and tissue collection
The investigation was conducted in conformity with the principles outlined in the Helsinki Declaration and with the Italian laws and guidelines, and was authorized by local Ethics Committee (protocol #2438, 27/01/2009). All specimens were taken after informed consent disclosing future use for research. Left ventricle (LV) cardiac biopsies from patients affected by non end-stage dilated hypokinetic ischemic cardiomyopathy where collected during surgical ventricular restoration procedure performed as described previously [17, 18]. For each patient, two biopsies were collected from the non-ischemic and dysfunctional remote myocardium: one was immediately immersed in RNAlater (Qiagen GmbH) and stored at 4 °C for <24 h before RNA extraction, and the other one was fixed in 10 % neutral buffered formalin (NBF) and paraffin embedded for RNA in situ hybridization and immunohistochemistry assays. HF patients’ characteristics are described in Table 1. Peripheral blood mononuclear cells (PBMC) were isolated from the peripheral blood of HF patients the day before surgery or from healthy controls by Histopaque Ficoll (Sigma Diagnostics, St. Louis, MO, USA) gradient centrifugation as described previously .
We have also analyzed 11 post-ischemic end-stage heart failure patients [10 males and 1 female, aged 59 ± 2.4 years (median ± SE)]. Among these, 8 LV samples were collected by an expert pathologist from explanted hearts of patients that underwent cardiac transplantation; the remaining 3 end-stage samples were collected during the implantation of LV assist device (LVAD). Care was taken to avoid necrotic or fibrotic areas. Samples of ~2 mm3 were snap frozen in liquid nitrogen and kept at −80 °C until they were analyzed. The end-stage HF patients were classified as NYHA class 3 (82 %) and class 4 (18 %), showed a LV transverse diameter of 130 ± 5.3 mm (median ± SE), measured by the pathologist after explantation, and ejection fraction of 22 ± 2.9 % (median ± SE). Sample harvesting was conducted after approval of the Ethics Committee of Udine (2 August 2011, reference number 47,831), in accordance with the Declaration of Helsinki, once written informed consent was obtained from each enrolled patient.
The control group was composed by age- and sex-matched subjects, who had died for causes different from stroke, ischemia, or cachexia for chronic diseases (n = 17, females/males = 6/10; 58.3 ± 3.4 years old). Left ventricle samples were excised and processed with less than 30 min cold ischemic time and snap frozen (XpressBANK, Asterand Biosciences).
Mouse transverse aortic constriction (TAC)
All experimental procedures complied with the Guidelines of the Italian National Institutes of Health and with the Guide for the Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources, National Academy of Sciences, Bethesda, MD, USA) and were approved by the institutional Animal Care and Use Committee (IACUC n. 709). TAC was performed in 2 months old C57BL/6 J male mice with a 27-gauge needle . Before all surgical procedures, mice were anesthetized with an intraperitoneal injection of 10 mg/kg Xilazine (Intervet Farmaceutici, Italy) and 100 mg/kg Ketamine (Ketavet 100; Intervet Farmaceutici, Italy). After the surgery, mice were allowed to recover at 37 °C. Sham operated mice received the same procedure except for the ligation of the aorta. Only mice showing a pressure gradient >60 mmHg measured by Doppler echocardiography 7 days after TAC were included in the analysis. Sham-operated mice were used as controls .
Formalin-fixed and paraffin-embedded (FFPE) genomic DNA isolation and single nucleotide polymorphism (SNP) analysis
Genomic DNA was isolated from twenty 10 μm-thick paraffin-embedded tissue sections. Following deparaffinization twice for 5 min in xylene, DNA was extracted using QIAamp DNA FFPE Tissue kit (Qiagen, USA) following manufacturer’s instructions. Four CDKN2-AS SNPs (RS1333040, RS1333049, RS10757278 and RS2383207) were assessed by TaqMan SNP Genotyping Assay on ABI 7900 real time PCR platform (Life Technologies, Thermo Fisher Scientific Inc., MA, USA). 20 ng of DNA were amplified in a reaction volume of 25 μl, containing 12.5 μl of 2× TaqMan Master Mix, 1.25 μl of 20× Assay working stock solution and 6.25 μl of nuclease-free water. Real time PCR conditions were: one cycle of 2 min at 50 °C, one cycle of 10 min at 95 °C; 40 cycles of 15 s at 95 °C and 1 min at 60 °C.
RNA isolation, lncRNA profiling and RT-qPCR
Total RNA from tissues was extracted using TRIzol (Life Technologies, Thermo Fisher Scientific Inc., MA, USA) as described previously [17, 22]. Sizing, quantitation and quality of the extracted RNAs was checked by Nanodrop ND-1000 (Nanodrop, Thermo Fisher Scientific Inc., MA, USA) and Bioanalyzer 2100 (Agilent Technologies). Long noncoding RNA expression profiles were measured using the disease-related human lncRNA Profiler (SBI, System Biosciences) (Additional file 1: Table S1). One µg of total RNA was retro-transcribed using the SuperScript III Reverse Transcriptase kit (Life Technologies, Thermo Fisher Scientific Inc., MA, USA) according to the manufacturer’s instructions. cDNAs were analyzed using the SYBR-GREEN qPCR method (Life Technologies, Thermo Fisher Scientific Inc., USA) according to the manufacturer’s instructions. Data are deposited in Gene Expression Omnibus repository (GEO GSE77399 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE77399). After median Ct value normalization, relative RNA expression was calculated with the 2−ΔΔCt method . Significantly (p < 0.05) modulated lncRNAs with a Ct ≤ 33 in at least one group and displaying a 2−ΔΔCt fold change ≥|1| were used for further validation. Validation was carried out with newly designed specific primers by SYBR-GREEN qPCR (Additional file 1: Table S2), using UBC for normalization.
RNA in situ hybridization assay and immunohistochemistry
In situ hybridization was performed using the QuantiGene ViewRNA system (Affymetrix) according to the manufacturer’s instructions. Briefly, 3-μm-thick sections were derived from FFPE-embedded biopsies and mounted on Superfrost Plus Gold glass slides (Thermo Fisher Scientific). Next, deparaffinized sections were hybridized with oligonucleotide probes conjugated with alkaline phosphatase (LP-AP) type 1, followed by staining with fast blue substrate, counterstaining with Hoechst 3332 (Sigma-Aldrich Co.) and mounted with VectaMount AQ medium (Vector Lab., USA).
To identify endothelial cells, serial LV sections were deparaffinized, microwave-treated and incubated with anti-human CD31 mouse monoclonal antibody (M0823, Dako). After incubation with biotinylated secondary antibodies and with ABC complex (Vectastain), the reaction was revealed with diaminobenzidine (DAB, vector) and counterstained with Hematoxylin.
For both in situ hybridization and immunohistochemistry, at least two blinded readers carried out the analysis and random images were acquired using an Axio Imager M.1 microscope equipped with an Axiocam MRc5 camera (Zeiss) and AxioVision software (Zeiss).
To analyze the relationship between lncRNAs and mRNAs differentially expressed in the same HF samples, we used the dataset GEO GSE26887 . We identified the mRNA transcripts significantly correlated to the differentially expressed lncRNAs by using Co-LncRNA platform  and Spearman’s rank correlation test, considering p ≤ 0.05 as statistical significance threshold. The lists of the correlated mRNAs were used to predict the enriched pathways by using WebGestalt .
For neighboring gene analysis, the 300 kb of genomic sequence upstream and downstream the HF modulated lncRNAs were identified by using UCSC Genome Browser (GRCh38/hg38 assembly) . Significant differential expression of neighboring genes lying in these regions was determined according to GSE26887 dataset .
For the comparison of transcriptomic changes in LV and PBMCs of HF patients, the following GEO datasets were analyzed: for LV, GSE26887 , and for PBMCs, GSE9128  and GSE1869 . Enriched pathways were obtained from each list of significantly (p ≤ 0.01) differentially expressed mRNAs and common deregulated pathways were plotted as a Venn’s diagram using BioVenn web application .
Continuous variables are expressed as mean ± standard error. Gaussian distribution was tested by using the Kolmogorov–Smirnov test. For group-wise comparisons, Mann–Whitney test (two groups) or t test (two groups) were used as appropriate. All tests were performed 2-sided and a p ≤ 0.05 was considered as statistically significant. For statistical analysis, GraphPad Prism v.4.03 software (GraphPad Software Inc.) was used. For CDKN2B-AS SNPs analysis, observed allele frequencies were compared to reported allele frequencies in HapMap–CEU European using an exact multinomial test using R package EMT .
Identification of lncRNAs deregulated in HF patients
In order to identify lncRNAs deregulated in HF, LV biopsies of 13 HF patients and 12 age-and sex-matched controls were analyzed. Worth noting is that myocardial biopsies were harvested from the non-ischemic portion of the LV (remote zone) in non end-stage HF patients, allowing to investigate the molecular mechanisms underpinning the disease during its progression. As expected, histological analysis revealed cardiomyocyte hypertrophy (Additional file 2: Figure S1B), in keeping with decreased alpha-myosin heavy chain (MYH6) and increased Natriuretic Peptide A (NPPA) mRNA levels (Additional file 2: Figure S1C, D).
Total RNA was extracted from LV biopsies of HF patients and controls, checked for quality (Additional file 2: Figure S2) and the gene expression levels of 83 disease-related lncRNAs (Additional file 1: Table S1), were measured by RT-qPCR. LncRNAs profiling showed that 53 and 33 lncRNAs were detectable (Ct ≤ 33) in HF patients and controls, respectively. We found that 27 lncRNAs were significantly modulated (20 up- and 7 down) at least onefold in HF compared with CTR subjects (Additional file 2: Figure S3A). Since the profiling results relied on commercial primers with undisclosed sequences, after designing new couples of unique primers, interrogating the Reference Sequence  of each lncRNA (Additional file 1: Table S2), significantly modulated lncRNAs were validated by RT-qPCR assays in more patients and controls (18 and 17, respectively). We found that 13 lncRNAs (10 up- and 3 downregulated) were significantly modulated in HF patients compared with CTR subjects, identifying a HF lncRNA signature (Fig. 1 and Additional file 2: Figure S3B). LOC285194 (also known as TUSC7) was down-modulated, but exhibited only borderline significance (p < 0.055). One additional limitation applies to LOC285194, as well as to EGOT, since the primer couples interrogating the relevant reference sequences did not interrogate certain shorter isoforms (Additional file 1: Table S2).
Unfortunately, the number of patients analyzed precluded further stratification or a correlation analysis with clinical parameters.
CDKN2B-AS gene (also known as ANRIL) is located in a region with several SNPs that correlate to increased genetic susceptibility to coronary artery diseases and type 2 diabetes [32, 33]. In genomic DNA extracted from 26 HF patients FFPE sections, we analyzed the expression of four CDKN2B-AS SNPs (RS1333040, RS1333049, RS10757278 and RS2383207) already known to be associated with cardiovascular diseases [34, 35]. No statistically significant difference between observed and HapMap allele frequencies (Additional file 1: Table S3) or correlation with RNA expression levels (data not shown) were found, possibly due to the low number of subjects analyzed.
Validation in chronic end-stage HF patients
In order to validate the identified HF lncRNA signature, we analyzed 11 RNAs derived from left ventricle samples of patients with chronic ischemic end-stage HF. These patients exhibited more severe clinical conditions and LV dilation compared to the patients used for lncRNA profiling.
We found that 9 lncRNAs (CDKN2B-AS1, EGOT, H19, HOTAIR, LOC285194, RMRP, RNY5, SOX2-OT and SRA1) were significantly modulated in a concordant manner in both end- and non end-stage HF patients (Fig. 2; Additional file 1: Table S5). HOXA11-AS was similarly modulated as well, but failed to reach statistical significance. Conversely, the expression of LOC285194 displaying a borderline significant decrease in non end-stage patients was strongly and very significantly (p < 0.001) inhibited in end-stage patients. Finally, RMST expression level was increased in non end-stage HF patients and significantly decreased in end-stage HF patients, underlining the differences between the two patient groups.
HF lncRNAs modulation in a mouse model of cardiac hypertrophy
Cardiomyocyte hypertrophy is a hallmark of post-ischemic HF patients (Additional file 2: Figure S1). Thus, we used a mouse model of cardiac hypertrophy obtained by transverse aortic constriction (TAC) to assay the expression of the identified lncRNAs. Gene expression was analyzed at day 7 after TAC, in order to examine early responses to hypertrophy, avoiding other potentially confounding effects. As expected [36, 37], SERCA2 (atp2a2) transcript levels decreased, while natriuretic peptides A and B (nppa and nppb) as well as Actin α1 (acta) transcripts increased in TAC hearts compared to controls (Additional file 2: Figure S4). When the expression of mouse-conserved lncRNAs (H19, HOXA11-AS, RMST, RMRP, SOX2OT and HOTAIR, Additional file 2: Figure S5) was assayed, we found that h19, rmrp and hotair were significantly modulated in a manner similar to the human counterpart (Fig. 3).
HF lncRNAs expression in the peripheral blood
Peripheral blood mononuclear cells (PBMCs) are a particularly attractive biomarker source because of the accessibility of peripheral blood and its straightforward preparation. Furthermore, inflammation and the underlying cellular and molecular mechanisms seem to play a crucial pathological role in the progression toward HF [38, 39]. Accordingly, gene ontology analysis of the transcriptomic alterations observed in the LV and PBMCs of HF patients showed that the majority of the dysregulated pathways and functions were in common between the two tissues (Additional file 1: Table S8; Additional file 2: Figure S7).
Thus, we measured the expression of the HF deregulated lncRNAs in PBMCs from 25 HF patients and 18 healthy individuals by RT-qPCR. Interestingly, we found that CDKN2B-AS1, HOTAIR and LOC285194 showed similar modulation in PBMCs and heart tissue, suggesting a potential as disease biomarkers (Fig. 4).
Cell localization of HF lncRNAs
In situ hybridization assay can provide important information about the cellular and sub-cellular distribution of lncRNAs. However, this assay is technically challenging, especially for low expressed transcripts [40, 41]. We managed to detect RMRP, one of the most abundant among the HF lncRNAs, according to the Ct values in RT-qPCR (19.8 ± 0.3). Indeed, while absolute Ct values are not quantitative, low Ct levels generally correspond to high expression levels . Figure 5 shows that RMRP displayed a ubiquitous distribution and nuclear intracellular localization.
Global H19 levels in the heart did not seem particularly high (25.8 ± 0.3). However, we were able to detect H19 in left ventricle sections, possibly due to its localized accumulation, mainly interstitial; indeed, the staining of serial sections with H19 and anti-CD31 antibody indicated a likely vascular localization (Fig. 6). Accordingly, H19 was readily detectable by RT-qPCR in cultured endothelial cells . In keeping with results obtained in other tissues [44, 45], H19 RNA staining was both cytoplasmic and nuclear.
Interactions between coding and noncoding transcriptomic changes
In order to gain insight into the role played by the HF lncRNAs, we took advantage of a previous transcriptomic analysis performed on the same RNAs using microarray (GEO GSE26887) .
mRNAs and HF lncRNAs were used for correlation analysis, which identified 10,257 and 8852 transcripts displaying significant direct or inverse correlation, respectively, with at least one HF lncRNA (Additional file 1: Table S6; Additional file 2: Figure S5). Next, we used the list of mRNAs correlating to each HF lncRNA for gene ontology analysis, and then selected those enriched pathways found in at least three HF lncRNAs. Indeed, we reasoned that different HF lncRNAs might have an additive effect in modulating specific pathways. We identified several important pathways for cardiovascular disease, such as insulin, IGF1 and glucose signaling pathways, TGF-β pathway, as well as hypoxic and oxygen homeostasis regulation of HIF-1-α (Fig. 7, and for a complete list, Additional file 1: Table S7).
LncRNAs can act in cis to regulate neighboring genes . Thus we studied the expression of coding genes located 300 kb upstream or downstream the deregulated lncRNAs in non end-stage HF LV samples. To this aim, we analyzed the differentially expressed mRNAs in GSE26887 dataset  and identified 19 lncRNA/mRNA couples potentially involved in HF disease mechanisms (Additional file 1: Table S4).
Albeit correlative, these analyses strongly suggest lncRNA involvement in the molecular mechanisms underpinning HF.
Several lncRNAs have been shown to be involved in specific physiological and pathological processes [5, 8], as well as in cardiovascular diseases [9–16]. We identified 14 lncRNAs deregulated in non end-stage HF patients. The validity of these findings was confirmed by the fact that nine of these lncRNAs displayed a concordant deregulation in an independent group of end-stage HF patients. Interestingly, LOC285194 that was weakly modulated in non end-stage patients was strongly down modulated in explanted failing hearts, and RMST displayed an inverse modulation in end-stage and non end-stage patients. Little is known about LOC285194 and RMST functions, but it is possible to speculate that differential regulation in the two HF groups might be linked to differences in disease progression. LOC285194 was previously shown to suppress tumor cell growth  and is the antisense transcript of LSAMP that has been reported to be a tumor suppressor gene . Interestingly, LSAMP down-regulation in coronary artery diseases is associated to atherosclerosis burden . RMST physically interacts with the transcription factor SOX2 and together regulate a large pool of downstream genes implicated in neurogenesis .
It is well known that, while miRNAs seem to be highly conserved, longer transcripts are under diverse levels of evolutionary constraints in mammalian. As a matter of fact, some lncRNAs show a high level of nucleotide sequence conservation (>60 %) in mammalian , and others do not display extensive evolutionary conservation , precluding them from being studied using mouse models. Nevertheless, in order to investigate the potential involvement of the identified lncRNAs in one important aspect of HF, we used a mouse model of cardiac hypertrophy. In spite of the multiple differences, we found that RMRP and H19 levels increased while HOTAIR levels decreased in both human failing and mouse hypertrophic hearts, thus opening the way for further functional studies.
Also interesting is the deregulation of CDKN2B-AS1, HOTAIR and LOC285194 in PBMCs derived from HF patients. Given the importance of inflammation in HF progression , it is tempting to speculate that these lncRNAs might respond to inflammatory stimuli in both myocardium and PBMCs. Moreover, it is also possible that these concomitant regulations are, at least in part, due to lncRNA transfer from one cell type to another via exosomes or other vesicles [52, 53].
Since peripheral blood can be obtained with minimally invasive procedures, CDKN2B-AS1, HOTAIR and LOC285194 might have a potential as circulating biomarkers of HF. Further studies in larger patients groups are necessary to assess their validity as disease biomarkers.
Pathway enrichment analysis highlighted several important cardiovascular categories enriched in HF lncRNAs-correlated mRNAs. Although correlative , this analysis fits well with previously published functional studies, suggesting that the identified HF lncRNAs might have a role in HF progression and prompting validation with more direct approaches. Among the significantly downregulated lncRNAs, HOTAIR (HOX antisense intergenic RNA) is transcribed from the antisense strand of homeobox C gene locus, and, in coordination with chromatin modifying enzymes, regulates gene silencing [55, 56]. Down-regulation of HOTAIR transcript has been found in aortic valve cells exposed to cyclic stretch, a modulation mediated through WNT/β-CATENIN pathway . Accordingly, we have observed the decrease of HOTAIR in a mouse model of cardiac hypertrophy due to pressure overload, indicating a possible involvement of HOTAIR in the left ventricle remodeling associated to HF. Moreover, both WNT and cardiac hypertrophy pathways were enriched categories in the gene ontologies analysis of mRNAs correlated to HOTAIR. Equally interesting is another downregulated lncRNA, the cyclin kinase Cdk9 inhibitor 7SK, which has been recently demonstrated to be involved in cardiac hypertrophy . Accordingly, histological and biomarkers analysis showed cardiac hypertrophy in the non end-stage HF patients studied. Unfortunately, lack of conservation between humans and mice precluded the analysis of 7SK expression in TAC mice.
CDKN2B-AS1 is the antisense gene of CDKN2B (or p15ink4a), which is located in the 50 kb chromosomal locus 9p21 . CDKN2B-AS1 represses the transcription of the genes in the INK4 locus by direct binding to the INK4b transcript and by recruiting the polycomb repressor complex (PRC) [60, 61].
CDKN2B-AS1 is the strongest genetic marker of human atherosclerosis and is generally considered the ‘golden standard’ for any genome wide association study of atherosclerosis-related traits [32, 33]. We analyzed the expression of four well-known cardiovascular diseases-associated variants, but we were not able to find any frequency association with the HF status, probably due to the limited number of patients analyzed.
We found that CDKN2B-AS1 RNA levels increased in HF patients compared to controls, but we did not find any correlation with the analyzed SNPs. Indeed, the association of CDKN2B-AS1 RNA expression levels and 9p21.3 risk alleles is still controversial . Nonetheless, CDKN2B-AS1 RNA expression has been shown to stimulate cell proliferation, adhesion, and to reduce apoptosis, providing a potential atherosclerosis disease mechanism . Noteworthily, we found that CDKN2B-AS1 RNA was also upregulated in PBMCs of HF patients, further implicating CDKN2B-AS1 in HF pathogenic mechanisms.
The steroid receptor RNA activator 1 (SRA1) gene generates both steroid receptor RNA activator protein (SRAP) as well as several non-coding SRA transcripts, depending on alternative transcription start site usage and alternative splicing [64, 65]. Friedrichs et al.  identified a gene cluster including SRA1 on a 600-kb linkage disequilibrium block on chromosome 5q31. 2–3 associated with human dilated cardiomyopathy in three independent Caucasian populations. Moreover, a role in heart development has been observed in zebrafish . Accordingly, knock-down of SRA1 impairs cardiac function in zebrafish . SRA1 is also known to stimulate cell proliferation as well as apoptosis in vivo , suggesting that SRA1 may be involved in HF pathogenesis.
H19 is a developmentally regulated gene with putative tumor suppressor activity . Here we showed the increase of H19 levels in HF patients, both end-stage and non, and confirmed its upregulation in a mouse model of cardiac hypertrophy . H19 is expressed during development of rat aorta, decreases in adult, and, interestingly, increases after vascular injury both in vivo and in vitro , as well as upon hypoxic stimulus [43, 71, 72]. Moreover, hyperhomocysteinemia, an independent risk factor for coronary artery diseases (CAD), increases the expression of H19 in aorta and vascular smooth muscle cells [73, 74], indicating that upregulated H19 may participate in the progression of CAD, the most common cause of HF. Recently, it has also been shown that polymorphisms in H19 are correlated with CAD  and CAD risk factors, such as obesity, high birth weight, and hypertension [76, 77]. In this respect, the mainly vascular localization of H19 in the heart appears particularly significant.
Among the protein coding genes neighbor to the lncRNAs, it is interesting the concordant modulation of HOTAIR and SMUG1, which is related to free-radicals response . Additionally, 7SK is inversely modulated compared to GSTA4 and 5, that both are involved in gluthatione metabolism . Moreover, we observed that Cathepsin D was modulated concordantly to its neighbor lncRNA H19. Cathepsin D is an autophagy-related enzyme and it has been found up-regulated in a hamster model of dilated cardiomyopathy .
We identified 14 lncRNAs that are dysregulated in non end-stage HF patients. Validation in other patient groups and in animal experimental models corroborated the findings, and the deregulation of some of them in the peripheral blood suggests a potential as disease biomarkers. While future investigations of the identified HF lncRNA actions and dysfunctions are required, analysis of the correlated genes indicates their implication in the molecular mechanisms underpinning HF progression.
long noncoding RNA
peripheral blood mononuclear cells
mouse transverse aortic constriction
formalin-fixed and paraffin-embedded
single nucleotide polymorphism
coronary artery diseases
Hunt SA, Abraham WT, Marshall H, Feldman AM, Francis GS, Ganiats TG, Jessup M, Konstam MA, Mancini DM, Michl K, Oates JA, Peter S, Silver MA, Stevenson LW, Yancy CW, Antman EM, Smith SC, Adams CD, Anderson JL, David P, Fuster V, Halperin JL, Hiratzka LF, Hunt SA, Jacobs AK, Nishimura R, Ornato JP, Page RL, Riegel B. ACC/AHA practice guidelines ACC/AHA 2005 guideline update for the diagnosis and management of chronic heart failure in the adult. Circulation. 2005;112:154–235.
Creemers EE, Wilde AA, Pinto YM. Heart failure: advances through genomics. Nat Rev Genet. 2011;12:357–62.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.
Mattick JS, Makunin IV. Non-coding RNA. Hum Mol Genet. 2006;15:17–29.
Greco S, Gorospe M, Martelli F. Noncoding RNA in age-related cardiovascular diseases. J Mol Cell Cardiol. 2015;83:142–55.
van Rooij E, Olson EN. MicroRNA therapeutics for cardiovascular disease: opportunities and obstacles. Nat Rev Drug Discov. 2012;11:860–72.
Almeida MI, Reis RM, Calin GA. MicroRNA history: discovery, recent applications, and next frontiers. Mutat Res. 2011;717:1–8.
Taft RJ, Pang KC, Mercer TR, Dinger M, Mattick JS. Non-coding RNAs: regulators of disease. J Pathol. 2010;220:126–39.
Dimmeler S, Marian AJ, Olson E, Thum T, Condorelli G. Long noncoding RNAs and microRNAs in cardiovascular pathophysiology. Circ Res. 2015;116:751–63.
Papait R, Kunderfranco P, Stirparo GG. Long noncoding rna: a new player of heart failure? J Cardiovasc Transl Res. 2013;6:876–83.
Beltrami C, Angelini TG, Emanueli C. Noncoding RNAs in diabetes vascular complications. J Mol Cell Cardiol. 2014;89:1–9.
Di Salvo TG, Guo Y, Su YR, Clark T, Brittain E, Absi T, Maltais S, Hemnes A. Right ventricular long noncoding RNA expression in human heart failure. Pulm Circ. 2015;5:135–61.
Di Salvo TG, Yang KC, Brittain E, Absi T, Maltais S. Right ventricular myocardial biomarkers in human heart failure. J Card Fail. 2015;21:398–411.
Matkovich SJ, Edwards JR, Grossenheider TC, Strong CDG, Ii GWD. Epigenetic coordination of embryonic heart transcription by dynamically regulated long noncoding RNAs. Proc Natl Acad Sci. 2014;111:1–6.
Devaux Y, Zangrando J, Schroen B, Creemers EE, Pedrazzini T, Chang CP, Dorn GW 2nd, Thum T, Heymans S. Long noncoding RNAs in cardiac development and ageing. Nat Rev Cardiol. 2015;12:415–25.
Ounzain S, Burdet F, Ibberson M, Pedrazzini T. Discovery and functional characterization of cardiovascular long noncoding RNAs. J Mol Cell Cardiol. 2015;89:17–26.
Greco S, Fasanaro P, Castelvecchio S, D’Alessandra Y, Arcelli D, Di Donato M, Malavazos A, Capogrossi MC, Menicanti L, Martelli F. MicroRNA dysregulation in diabetic ischemic heart failure patients. Diabetes. 2012;61:1633–41.
Castelvecchio S, Menicanti L, Di Donato M. Surgical ventricular restoration to reverse left ventricular remodeling. Curr Cardiol Rev. 2010;6:15–23.
Voellenkle C, van Rooij J, Cappuzzello C, Greco S, Arcelli D, Di Vito L, Melillo G, Rigolini R, Costa E, Crea F, Capogrossi MC, Napolitano M, Martelli F. MicroRNA signatures in peripheral blood mononuclear cells of chronic heart failure patients. Physiol Genom. 2010;42:420–6.
Hu P, Zhang D, Swenson L, Chakrabarti G, Abel ED, Litwin SE. Minimally invasive aortic banding in mice: effects of altered cardiomyocyte insulin signaling during pressure overload. Am J Physiol Heart Circ Physiol. 2003;285:H1261–9.
Tanaka N, Dalton N, Mao L, Rockman HA, Peterson KL, Gottshall KR, Hunter JJ, Chien KR, Ross JJ. Transthoracic echocardiography in models of cardiac disease in the mouse. Circulation. 1996;94:1109–17.
Greco S, De Simone M, Colussi C, Zaccagnini G, Fasanaro P, Pescatori M, Cardani R, Perbellini R, Isaia E, Sale P, Meola G, Capogrossi MC, Gaetano C, Martelli F. Common micro-RNA signature in skeletal muscle damage and regeneration induced by Duchenne muscular dystrophy and acute ischemia. FASEB J. 2009;23:3335–46.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods. 2001;25:402–8.
UCSC Genome bioinformatics [https://genome.ucsc.edu/].
Cappuzzello C, Napolitano M, Arcelli D, Melillo G, Melchionna R, Di Vito L, Carlini D, Silvestri L, Brugaletta S, Liuzzo G, Crea F, Capogrossi MC. Gene expression profiles in peripheral blood mononuclear cells of chronic heart failure patients. Physiol Genom. 2009;38:233–40.
Kittleson MM, Minhas KM, Irizarry RA, Ye SQ, Breton E, Conte JV, Tomaselli G, Garcia JGN, Hare JM, Kittleson MM, Minhas KM, Irizarry RA, Ye SQ, Edness G, Breton E, Conte JV, Tomaselli G, Garcia JGN, Hare JM, Ye SQ, Edness G, Breton E, Conte JV, Garcia JGN, Hare JM. Gene expression analysis of ischemic and nonischemic cardiomyopathy: shared and distinct genes in the development of heart failure. Physiol Genom. 2007;21:299–307.
Menzel U. EMT: exact multinomial test: goodness-of-fit test for discrete multivariate data. R package version 1.1. 2013.
RefSeq. NCBI reference sequence database [http://www.ncbi.nlm.nih.gov/refseq/].
Cugino D, Gianfagna F, Santimone I, de Gaetano G, Donati MB, Iacoviello L, Di Castelnuovo A. Type 2 diabetes and polymorphisms on chromosome 9p21: a meta-analysis. Nutr Metab Cardiovasc Dis. 2012;22:619–25.
Holdt LM, Teupser D. From genotype to phenotype in human atherosclerosis–recent findings. Curr Opin Lipidol. 2013;24:410–8.
Lin H-F, Tsai P-C, Lin R-T, Khor G-T, Sheu S-H, Juo SHH. Sex differential genetic effect of chromosome 9p21 on subclinical atherosclerosis. PLoS ONE. 2010;5:e15124.
Scheffold T, Kullmann S, Huge A, Binner P, Ochs HR, Schols W, Thale J, Motz W, Hegge FJ, Stellbrink C, Dorsel T, Gulker H, Heuer H, Dinh W, Stoll M, Haltern G. Six sequence variants on chromosome 9p21.3 are associated with a positive family history of myocardial infarction: a multicenter registry. BMC Cardiovasc Disord. 2011;11:9.
Garcia-Menendez L, Karamanlidis G, Kolwicz S, Tian R. Substrain specific response to cardiac pressure overload in C57BL/6 mice. Am J Physiol Heart Circ Physiol. 2013;305:H397–402.
Liu Y, Shen Y, Zhu J, Liu M, Li X, Chen Y, Kong X, Song G, Qian L. Cardiac-specific PID1 overexpression enhances pressure overload-induced cardiac hypertrophy in mice. Cell Physiol Biochem. 2015;35:1975–85.
Di Stefano V, Zaccagnini G, Capogrossi MC, Martelli F. MicroRNAs as peripheral blood biomarkers of cardiovascular disease. Vasc Pharmacol. 2011;55:111–8.
Anker SD, von Haehling S. Inflammatory mediators in chronic heart failure: an overview. Heart. 2004;90:464–70.
Levsky JM, Singer RH. Fluorescence in situ hybridization: past, present and future. J Cell Sci. 2003;116(Pt 14):2833–8.
Itzkovitz S, van Oudenaarden A. Validating transcripts with probes and imaging technology. Nat Methods. 2011;8(4 Suppl):S12–9.
Bustin SA, Benes V, Nolan T, Pfaffl MW. Quantitative real-time RT-PCR–a perspective. J Mol Endocrinol. 2005;34:597–601.
Voellenkle C, Garcia-Manteiga JM, Pedrotti S, Perfetti A, De Toma I, Da Silva D, Maimone B, Greco S, Fasanaro P, Creo P, Zaccagnini G, Gaetano C, Martelli F. Implication of long noncoding RNAs in the endothelial cell response to hypoxia revealed by RNA-sequencing. Sci Rep. 2016;6:24141.
Berteaux N, Aptel N, Cathala G, Genton C, Coll J, Daccache A, Spruyt N, Hondermarck H, Dugimont T, Curgy J-J, Forne T, Adriaenssens E. A novel H19 antisense RNA overexpressed in breast cancer contributes to paternal IGF2 expression. Mol Cell Biol. 2008;28:6731–45.
van Heesch S, van Iterson M, Jacobi J, Boymans S, Essers PB, de Bruijn E, Hao W, MacInnes AW, Cuppen E, Simonis M. Extensive localization of long noncoding RNAs to the cytosol and mono- and polyribosomal complexes. Genome Biol. 2014;15:R6.
Greco S, Gorospe M, Martelli F. Noncoding RNA in age-related cardiovascular diseases. J Mol Cell Cardiol. 2015;83:142–55.
Liu Q, Huang J, Zhou N, Zhang Z, Zhang A, Lu Z, Wu F, Mo Y-Y. LncRNA loc285194 is a p53-regulated tumor suppressor. Nucleic Acids Res. 2013;41:4976–87.
Chen J, Lui W-O, Vos MD, Clark GJ, Takahashi M, Schoumans J, Khoo SK, Petillo D, Lavery T, Sugimura J, Astuti D, Zhang C, Kagawa S, Maher ER, Larsson C, Alberts AS, Kanayama H, Teh BT. The t(1;3) breakpoint-spanning genes LSAMP and NORE1 are involved in clear cell renal cell carcinomas. Cancer Cell. 2003;4:405–13.
Ng S, Bogu GK, Soh BS, Stanton LW. The long noncoding RNA RMST interacts with SOX2 to regulate neurogenesis. Mol Cell. 2013;51:349–59.
Louro R, Smirnova AS, Verjovski-almeida S. Genomics Long intronic noncoding RNA transcription: expression noise or expression choice? Genomics. 2009;93:291–8.
Ulitsky I, Bartel DP. LincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154:26–46.
Vlassov AV, Magdaleno S, Setterquist R, Conrad R. Biochimica et Biophysica Acta Exosomes: current knowledge of their composition, biological functions, and diagnostic and therapeutic potentials. BBA Gen Subj. 2012;1820:940–8.
Loyer X, Vion AC, Tedgui A, Boulanger CM. Microvesicles as cell-cell messengers in cardiovascular diseases. Circ Res. 2014;114:345–53.
Khatri P, Draghici S. Ontological analysis of gene expression data: current tools, limitations, and open problems. Bioinformatics. 2005;21:3587–95.
Marchese FP, Huarte M. Long non-coding RNAs and chromatin modifiers: their place in the epigenetic code. Epigenetics. 2014;9:21–6.
Bhan A, Hussain I, Ansari KI, Kasiri S, Bashyal A, Mandal SS. Antisense transcript long noncoding RNA (lncRNA) HOTAIR is transcriptionally induced by estradiol. J Mol Biol. 2013;425:3707–22.
Carrion K, Dyo J, Patel V, Sasik R, Mohamed SA, Hardiman G, Nigam V. The long non-coding HOTAIR is modulated by cyclic stretch and WNT/b-CATENIN in human aortic valve cells and is a novel repressor of calcification genes. PLoS ONE. 2014;9:1–7.
Sano M, Schneider MD. Cyclins that don’t cycle–cyclin T/cyclin-dependent kinase-9 determines cardiac muscle cell size. Cell Cycle. 2003;2:99–104.
Holdt LM, Beutner F, Scholz M, Gielen S, Gabel G, Bergert H, Schuler G, Thiery J, Teupser D. ANRIL expression is associated with atherosclerosis risk at chromosome 9p21. Arter Thromb Vasc Biol. 2010;30:620–7.
Kotake Y, Nakagawa T, Kitagawa K, Suzuki S, Liu N, Kitagawa M, Xiong Y. Long non-coding RNA ANRIL is required for the PRC2 recruitment to and silencing of p15(INK4B) tumor suppressor gene. Oncogene. 2011;30:1956–62.
Yap KL, Li S, Munoz-Cabello AM, Raguz S, Zeng L, Mujtaba S, Gil J, Walsh MJ, Zhou M-M. Molecular interplay of the noncoding RNA ANRIL and methylated histone H3 lysine 27 by polycomb CBX7 in transcriptional silencing of INK4a. Mol Cell. 2010;38:662–74.
Pilbrow AP, Folkersen L, Pearson JF, Brown CM, McNoe L, Wang NM, Sweet WE, Tang WHW, Black MA, Troughton RW, Richards AM, FrancoCereceda A, Gabrielsen A, Eriksson P, Moravec CS, Cameron VA. The chromosome 9p21.3 coronary heart disease risk allele is associated with altered gene expression in normal heart and vascular tissues. PLoS ONE. 2012;7:1–11.
Holdt LM, Hoffmann S, Sass K, Langenberger D, Scholz M, Krohn K, Finstermeier K, Stahringer A, Wilfert W, Beutner F, Gielen S, Schuler G, Gabel G, Bergert H, Bechmann I, Stadler PF, Thiery J, Teupser D. Alu elements in ANRIL non-coding RNA at chromosome 9p21 modulate atherogenic cell functions through trans-regulation of gene networks. PLoS Genet. 2013;9:e1003588.
Scheuermann JC, Boyer LA. Getting to the heart of the matter: long non-coding RNAs in cardiac development and disease. EMBO J. 2013;32:1805–16.
Cooper C, Vincett D, Yan Y, Hamedani MK, Myal Y, Leygue E. Steroid receptor RNA Activator bi-faceted genetic system: heads or tails? Biochimie. 2011;93:1973–80.
Friedrichs F, Zugck C, Rauch G, Friedrichs F, Zugck C, Ivandic B, Weichenhan D, Mu M, Meder B, Eddine N, Mokhtari E, Regitz-zagrosek V, Hetzer R, Scha A, Schreiber S, Chen J, Neuhaus I, Ji R, Siemers NO, Frey N, Rottbauer W, Katus HA, Stoll M. HBEGF, SRA1, and IK: three cosegregating genes as determinants of cardiomyopathy. Genome Res. 2009;19:395–403.
Lanz RB, Chua SS, Barron N, So BM, Demayo F, Malley BWO. Steroid receptor RNA activator stimulates proliferation as well as apoptosis in vivo. Mol Cell Biol. 2003;23:7163–76.
Brunkow ME, Tilghman SM. Ectopic expression of the H19 gene in mice causes prenatal lethality. Genes Dev. 1991;5:1092–101.
Lee JH, Gao C, Peng G, Greer C, Ren S, Wang Y, Xiao X. Analysis of transcriptome complexity through RNA sequencing in normal and failing murine hearts. Circ Res. 2011;109:1332–41.
Kim D, Zhang L, Dzau VJ, Pratt RE. H19, a developmentally regulated gene, is reexpressed in rat vascular smooth muscle cells after injury. J Clin Invest. 1994;93:355–60.
Ayesh S, Matouk I, Schneider T, Ohana P, Laster M, Al-Sharef W, De-Groot N, Hochberg A. Possible physiological role of H19 RNA. Mol Carcinog. 2002;35:63–74.
Matouk IJ, Mezan S, Mizrahi A, Ohana P, Abu-lail R, Fellig Y, Galun E, Hochberg A. Biochimica et Biophysica Acta The oncofetal H19 RNA connection: hypoxia, p53 and cancer. BBA - Mol Cell Res. 2010;1803:443–51.
Devlin AM, Bottiglieri T, Domann FE, Lentz SR. Tissue-specific changes in H19 methylation and expression in mice with hyperhomocysteinemia. J Biol Chem. 2005;280:25506–11.
Li L, Xie J, Zhang M, Wang S. Homocysteine harasses the imprinting expression of IGF2 and H19 by demethylation of differentially methylated region between IGF2/H19 genes. Acta Biochim Biophys Sin (Shanghai). 2009;41:464–71.
Gao W, Zhu M, Wang H, Zhao S, Zhao D, Yang Y, Wang Z-M, Wang F, Yang Z-J, Lu X, Wang L-S. Association of polymorphisms in long non-coding RNA H19 with coronary artery disease risk in a Chinese population. Mutat Res Mol Mech Mutagen. 2015;772:15–22.
Hernandez-Valero MA, Rother J, Gorlov I, Frazier M, Gorlova O. Interplay between polymorphisms and methylation in the H19/IGF2 gene region may contribute to obesity in Mexican-American children. J Dev Orig Health Dis. 2013;4:499–506.
Tragante V, Barnes MR, Ganesh SK, Lanktree MB, Guo W, Franceschini N, Smith EN, Johnson T, Holmes MV, Padmanabhan S, Karczewski KJ, Almoguera B, Barnard J, Baumert J, Chang YPC, Elbers CC, Farrall M, Fischer ME, Gaunt TR, Gho JMIH, Gieger C, Goel A, Gong Y, Isaacs A, Kleber ME, Leach IM, McDonough CW, Meijs MFL, Melander O, Nelson CP, et al. Gene-centric meta-analysis in 87,736 individuals of European ancestry identifies multiple blood-pressure-related loci. Am J Hum Genet. 2014;94:349–60.
Dizdaroglu M, Jaruga P, Birincioglu M, Rodriguez H. Free radical-induced damage to DNA: mechanisms and measurement. Free Radic Biol Med. 2002;32:1102–15.
Conklin DJ, Guo Y, Jagatheesan G, Kilfoil PJ, Haberzettl P, Hill BG, Baba SP, Guo L, Wetzelberger K, Obal D, Rokosh DG, Prough RA, Prabhu SD, Velayutham M, Zweier JL, Hoetker JD, Riggs DW, Srivastava S, Bolli R, Bhatnagar A. Genetic deficiency of glutathione S-transferase p increases myocardial sensitivity to ischemia-reperfusion injury. Circ Res. 2015;117:437–49.
Takemura G, Miyata S, Kawase Y, Okada H, Maruyama R, Fujiwara H. Autophagic degeneration and death of cardiomyocytes in heart failure. Autophagy. 2006;2:212–4.
SG and FM conceived the project. LM, SC, APB and NF performed the patient’s biopsies and data collection. GZ was responsible for TAC mouse model experiments and for immunohistochemistry analysis. RV performed the SNPs analysis. PF and AP performed qPCRs of hypertrophy markers and lncRNAs, respectively. CV performed the gene ontology analysis. SG performed all the other experiments. SG and FM wrote the paper. Critical discussion of the data and revision of the manuscript by all the authors. All authors read and approved the final manuscript.
We thank F. Ambrogi, Policlinico San Donato-IRCCS, for statistical analysis, P. Schiavone, Policlinico San Donato-IRCCS, for histological analysis, P. Carullo, Humanitas, Rozzano (MI) for his help to set up TAC mouse model, A. Esposito and L. Perani, Preclinical Imaging Facility, Experimental Imaging Center, San Raffaele Scientific Institute, Milan for TAC mouse ultrasound imaging. Editing revision by Silvia Marchetti, Republic of San Marino, is gratefully acknowledged.
Availability of data and materials
Data are deposited in Gene Expression Omnibus repository (GEO GSE77399 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE77399).
The authors declare that they have no competing interests.
Consent to publish
All human specimens were taken after informed consent from the participants disclosing future use for research and for publishing.
Ethics approval and consent to participate
All specimens were taken after informed consent disclosing future use for research. The investigation was conducted in conformity with the principles outlined in the Helsinki Declaration and with the Italian laws and guidelines. The study on non end-stage HF patients was authorized by local Ethics Committee (protocol #2438, 27/01/2009). The study on end-stage specimens was conducted after approval of the Ethics Committee of Udine (2 August 2011, reference number 47,831).
All mouse experimental procedures complied with the Guidelines of the Italian National Institutes of Health and with the Guide for the Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources, National Academy of Sciences, Bethesda, MD, USA) and were approved by the institutional Animal Care and Use Committee (IACUC n. 709).
The financial support of Ministero della Salute (Ricerca Corrente, RF-2011-02347907 and PE-2011-02348537), Telethon-Italy (Grant No. GGP14092), AFM-Telethon (Grant No. 18477) and Cariplo Foundation (Grant No. 2013-0887) is gratefully acknowledged.