|Home | About | Journals | Submit | Contact Us | Français|
Recent studies in animal models and humans show that long non-coding RNAs (lncRNAs) are involved in the development of atherosclerosis, which contributes to the pathological foundation of coronary artery disease (CAD). LncRNAs in plasma and serum have been considered as promising novel biomarkers for diagnosis and prognosis of cardiovascular diseases, especially CAD. We here measured the circulating levels of 8 individual lncRNAs which are known to be relevant to atherosclerosis in the plasma samples from 300 patients with CAD and 180 control subjects by using quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) methods. We found that the plasma level of H19 and long intergenic non-coding RNA predicting cardiac remodeling (LIPCAR) were significantly increased in patients with CAD. The area under the receiver operating characteristic curve was 0.631 for H19 and 0.722 for LIPCAR. Multivariate logistic regression analyses indicated that plasma H19 and LIPCAR were independent predictors for CAD, even after adjustment for traditional cardiovascular risk factors. Further studies identified that plasma levels of H19 and LIPCAR were also increased in CAD patients with heart failure compared to those with normal cardiac function. Taken together, our results suggest that increased plasma levels of H19 and LIPCAR are associated with increased risk of CAD and may be considered as novel biomarkers for CAD.
Coronary artery disease (CAD) remains one of the major causes of mortality and morbidity in many countries, including China1. Numerous studies have identified several risk factors for CAD, including hypertension, dyslipidaemia, diabetes, obesity, smoking, dietary, gender, etc.2. Recently, genomics researches have revealed a series of new candidate biomarkers that may contribute to the pathogenesis of CAD3.
Long non-coding RNAs (lncRNAs), a novel class of non-coding RNAs, are defined as transcripts that are longer than 200 nt and lacking protein-encoding capacity4. LncRNAs play crucial roles in chromatin modification, imprinting, cell differentiation and proliferation, transcription, translation and other important biological processes5. Atherosclerosis, the main pathophysiological cause of CAD, is initiated by endothelial injury and activation, which leads to infiltration and proliferation of vascular smooth muscle cells (VSMC), leukocytes and other inflammatory cells in the arterial wall6. Recently, lncRNAs have emerged as important regulators in various pathological processes that contribute to the development of atherosclerosis7–9. For instance, the lincRNA metastasis associated lung adenocarcinoma transcript 1 (MALAT1) regulates blood vessel growth and MALAT1 inhibition prevents human endothelial cell proliferation and reduces vascular growth10. LncRNA H19 is highly expressed in the neo-intima after injuries and in human atherosclerotic lesions, but barely expressed in normal coronary arteries11, 12. LncRNA highly upregulated in liver cancer (HULC) and Apolipoprotein A1 antisense (APOA1-AS) can modulate multiple key lipometabolism-related genes and play important roles in lipid homeostasis13, 14. TNFα and hnRNPL related immunoregulatory LincRNA (THRIL) can promote the transcriptional process of TNF-α and induce inflammation15. By using microarray analysis in a rat model of ischemic heart disease, we also identified 331 pairs of differentially expressed lncRNAs and nearby coding genes, indicating that lncRNA might be also involved in the pathogenesis of CAD16.
To date, non-coding RNAs, including microRNAs and lncRNAs, in plasma and serum have been considered as promising novel biomarkers for diagnosis and prognosis of cardiovascular diseases17. We have previously demonstrated that lipometabolism-related microRNA-122 and microRNA-370 were associated with the risk and severity of CAD18. However, studies on circulating lncRNAs for the risk of CAD remain sparse. In the present study, we examined plasma levels of eight cardiac-related or atherosclerosis-related lncRNAs, including H19, long intergenic non-coding RNA predicting cardiac remodeling (LIPCAR), APOA1-AS, THRIL, HULC, SLC26A4-AS1, LincRNA-Cox2, LincRNA-p2111–15, 19–22, in patients with angiographically demonstrated CAD to investigate the possibility of these circulating lncRNAs as novel biomarkers for CAD.
From 2014 to 2016, 480 consecutive subjects (296 males and 184 females), aged 42–78 years, who underwent coronary angiography for suspected or known coronary atherosclerosis at the First Affiliated Hospital of Nanjing Medical University in China were enrolled in this study. Coronary artery disease (CAD) was defined as angiographic evidence of at least one segment of a major coronary artery, including the left anterior descending, left circumflex, or right coronary artery, with >50% organic stenosis. The severity of CAD was assessed by the Gensini score system based on the degree of luminal narrowing and its geographic importance23. Subjects with normal coronary arteries were considered as controls. Two cardiologists who were unaware of the patients included in this study assessed the angiograms. All subjects included in this study had no family history of CAD and no history of significant concomitant diseases, including hepatic failure, renal failure, hepatitis, cardiomyopathy, congenital heart disease, bleeding disorders, previous thoracic irradiation therapy, and malignant diseases. CAD patients were divided into three subgroups: stable angina pectoris (SAP), unstable angina pectoris (UAP), and acute myocardial infarction (AMI), which were defined as previously described24. The diagnosis of chronic heart failure (CHF) was made on the basis of typical symptoms and signs, and evidence of left ventricular enlargement and systolic functional impairment by echocardiography, according to the American College of Cardiology/American Heart Association guidelines25. Hypertension was defined as resting systolic blood pressure (SBP) above 140mmHg and/or diastolic blood pressure (DBP) above 90mmHg or in the presence of active treatment with antihypertensive agents. Diabetes mellitus was defined as fasting blood glucose (FBG) >7.0mmol/L or a diagnosis with diet adjustment or anti-diabetic drug therapy. Smoking was defined as >10 cigarettes per day. Written informed consent was obtained from each participant and this study was approved by the Ethics Committee of the First Affiliated Hospital of Nanjing Medical University. All of the experiments in the present study were performed in accordance with the relevant approved guidelines and regulations.
Fasting blood sample was collected from each subject and anticoagulated with ethylenediamine tetraacetic acid (EDTA) dipotassium salt in the early morning. The sample was separated immediately by centrifugation at 3000g for 15min at 4°C to retrieve plasma. The plasma was then stored at −80°C until assayed. Total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C) levels were measured enzymatically on a chemistry analyzer (Olympus Au2700, First Chemical Ltd., Tokyo, Japan). Glucose levels were measured by a glucose oxidase method (Reagent kit; Diagnostic Chemicals Ltd., Oxford, CT, UK).
Total RNA was isolated from 400μL of plasma using the mirVanaTM PARISTM Kit (Ambion, Austin, TX) according to the manufacturer’s instructions with modification. For normalization of sample-to-sample variation, 25 fmol of synthetic C.elegans miRNA cel-miR-39 (Qiagen, Germany) was added to each sample after addition of 2×Denaturing Solution (Ambion, Austin, TX) . RNA was dissolved in 100μL of RNase-free water, and then stored at −80°C until analysis. Total RNA was reverse transcribed using the DRR037A PrimeScript® RT Master mix (Takara, Dalian, P.R. China). The RT reaction was incubated at 37°C for 15min, at 85°C for 5s, and then held at 4°C. The cDNA product was stored at −20°C until analysis.
To detect plasma levels of lncRNA, 2μL of the cDNA product was used as template in 10μL reaction containing 5μL of TaqMan® Universal PCR Master Mix (Applied Biosystems, Foster, CA), 1μL of specific primer (Supplementary Table 1), 2μL of RNase-free water. qRT-PCR was performed with 7900HT real-time PCR system (Applied Biosystems, Foster, CA) at 95°C for 10min, followed by 40 cycles of 95°C for 15s and 60°C for 1min. Triplicate measurements were obtained for each sample on a 384-well plate. Data were analyzed with SDS Relative Quantification Software version 2.2.2 (Applied Biosystems, Foster, CA), with the automatic Ct setting for assigning baseline and threshold for Ct determination. The relative expression level of each individual lncRNA after normalization to cel-miR-39 was calculated using the 2−ΔΔCt method.
Normality of distribution was assessed using the Kolmogorov-Smirnov test. Comparison between two groups was performed with Student’s t tests or Mann–Whitney U tests. For comparison among more than 2 groups, one-way ANOVA or the Kruskal–Wallis test was used as appropriate. Pearson χ2 test was used to compare qualitative variables represented as frequencies. The correlations between plasma levels of lncRNAs and other variables were calculated using Spearman correlation coefficient. Univariate analysis and multi-variate logistic regression analysis were undertaken to determine the variables that independently contributed to the presence of CAD. Odds ratio (OR) and 95% confidence interval (CI) were also calculated. All tests were two-sided and P<0.05 was considered statistically significant. Receiver operating characteristic (ROC) curve and the area under ROC curve (AUC) were used to assess the sensitivity and specificity of lncRNA as a novel diagnostic tool for the detection of CAD.
Table 1 presents the characteristics of the study population. A total of 300 subjects with coronary artery disease (CAD) and 180 controls were enrolled in the study. Compared with the controls, patients with CAD had higher levels body mass index (BMI), TC, LDL-C, prevalence of hypertension, but lower HDL-C. No significant difference was found in age, gender, smoking, diabetes, creatinine, TG, and fasting blood glucose (FBG). The CAD group consists of 81 patients with stable angina pectoris (SAP), 189 patients with unstable angina pectoris (UAP), and 30 patients with acute myocardial infarction (AMI).
Compared with the control group, plasma levels of H19 and LIPCAR were higher in CAD patients (P <0.001, Figure 1A and B). However, no significant difference was observed in the plasma levels of THRIL, LincRNA-Cox2, LincRNA-p21, HULC, SLC26A4-AS1 and APOA1-AS between patients with CAD and controls (Fig. 1C–H). We further divided CAD patients into three subgroups, including stable angina pectoris (SAP), unstable angina pectoris (UAP), and acute myocardial infarction (AMI). As shown in Fig. 2, patients with AMI appeared to have the highest circulating levels of H19 and LIPCAR, although the difference did not reach statistical significance.
We further analyzed the correlations of plasma levels of H19 and LIPCAR with clinical characteristics in patients with CAD. As shown in Table 2, plasma levels of H19 was positively associated with BMI (R=0.121, P=0.022), LDL-C (R=0.134, P=0.012), and Gensini score (R=0.161, P=0.003), indicating that increased H19 level may correlated with the severity of CAD. In addition, we also found that the plasma level of LIPCAR was positively associated with age (R=0.201, P <0.001), but negatively associated with HDL-C (R=−0.203, P<0.001).
Univariate and multivariate logistic regression analysis revealed that plasma levels of H19 and LIPCAR were significantly associated with the presence of CAD, even after adjustment for age, gender, BMI, smoking, hypertension, diabetes, and blood lipid profiles, univariate and multivariate logistic regression analysis revealed that plasma levels of H19 and LIPCAR were significantly associated with the presence of CAD (Table 3).
Stratified analyses were conducted according to gender, age, diabetes and smoking status. As shown in Table 4, we found that the predictive effect of plasma H19 on the risk of CAD was more prominent in females (Adjusted OR=1.126, 95% CI=1.021–1.241, P=0.017), elderly (Adjusted OR=1.115, 95% CI=1.006–1.236, P=0.039) and non-diabetic subjects (Adjusted OR=1.092, 95%CI=1.013–1.177, P=0.021). While for LIPCAR, as shown in Table 5, the predictive effect on the risk of CAD was more prominent in younger subjects (Adjusted OR=1.306, 95% CI=1.061–1.607, P=0.012), non-diabetic subjects (Adjusted OR=1.227, 95%CI=1.090–1.382, P=0.001), and non-smoking subjects (Adjusted OR=1.682, 95% CI=1.198–2.361, P=0.003).
Having established that plasma H19 and LIPCAR are independent risk factors for CAD, we sought to determine the potential utility of plasma H19 and LIPCAR as diagnostic biomarkers of CAD. To this end, ROC analysis was performed to evaluate the predictive power of plasma H19 and LIPCAR for HF. Our results showed that the area under ROC curve (AUC) was 0.631 (95% CI=0.551–0.788) for H19 (Fig. 3A) and 0.722 (95% CI=0.669–0.782) for LIPCAR (Fig. 3B). The sensitivity and specificity at the optimal cut-off were 53.6% and 73.0% for H19, and 72.2% and 62.3% for LIPCAR, respectively.
Since both H19 and LIPCAR have been demonstrated to be involved in the pathological process of heart failure19, 26, we further investigated the differences of these two lncRNAs between CAD patienrs with and without CHF. As shown in Fig. 4, plasma levels of H19 and LIPCAR were both higher in patients with CHF (P=0.014 for H19; P=0.038 for LIPCAR, respectively).
In the present study, we investigated the plasma levels of a selected set of cardiac-related or atherosclerosis-related lncRNAs for their potential as novel biomarkers of coronary artery disease (CAD). We here showed that the plasma levels of two lncRNAs, H19 and LIPCAR were significantly increased in patients with CAD. Multivariate logistic regression analysis revealed that plasma levels of H19 and LIPCAR were independently associated with the risk of CAD, even after adjustment for traditional cardiovascular risk factors. Our results indicated that plasma H19 and LIPCAR might be served as promising candidate biomarkers for CAD.
The lncRNA H19 is a well-known imprinted gene which is abundantly expressed from the early stage of embryogenesis throughout fetal life, but is downregulated postnatally27. Accumulating data indicate that re-expression of H19 may play important roles in cardiovascular diseases28, 29. The expression levels of H19 were increased in human VSMCs treated with homocysteine and aortae of mice with hyperhomocysteinemia, which is a well-known independent risk factor for CAD30, 31. A recent study showed that over-expression of H19 could promote atherosclerosis by activating MAPK and NF-κB signaling pathway32. Our previous study also demonstrated that the common polymorphisms of H19 were associated with the risk and severity of CAD33. We here found that increased plasma level of H19 was independently correlated with the risk of CAD and positively associated with the severity of CAD evaluated by the Gensini score. Our stratified analysis revealed that the increased risk of CAD associated with the plasma levels of H19 was more prominent in subgroups of females, elderly, non-diabetic subjects. The most possible explanation for this sex and age difference is that as an imprinted gene, paternal allele of H19 is imprinted, while only the maternal allele is expressed34. Moreover, H19 is located in close proximity to the insulin-like growth factor 2 (IGF2) gene on human chromosome 11p15.5 and can downregulate the expression of IGF2 in cis and trans 35. IGF2 is one the key regulators in insulin signaling and is involved in the development of diabetes36. Thus, we speculated that abnormal insulin regulation in diabetes may also affect the predictive role of plasma H19 for CAD in our study. However, further studies are needed to confirm this hypothesis.
Another interesting finding in our study is that plasma H19 is increased in CAD patients with chronic heart failure (CHF). Our results are consistent with a recent study which showed that cardiac H19 level was increased in ischemic end-stage failing hearts37. Liu et al. showed that H19 could regulated the development of cardiac hypertrophy through miR-675/CaMKIIδ pathway26, indicating that dysregulation of H19 may also contribute to the pathogenesis of ischemic heart failure.
LIPCAR is considered as a mitochondria-derived lncRNA which can be readily detected in circulating19. Circulating level of LIPCAR is upregulated in patients with CHF independently of the pathogenesis, and that higher LIPCAR levels were associated with a higher risk of cardiovascular death19. We here also found that increased plasma LIPCAR level was independently correlated with increased risk of CAD. Further analyses showed that plasma LIPCAR levels were higher in CAD patients with CHF when compared to those without CHF. Our results are consistent with another study which demonstrated that circulating LIPCAR levels inversely correlated with echocardiographic E/A ratio, which is marker of LV dysfunction38. Our stratified analysis revealed that the increased risk of CAD associated with the plasma levels of LIPCAR was more prominent in subgroups of younger, non-diabetic, and non-smoking subjects. De Gonzalo-Calvo, D. et al. also showed that circulating LIPCAR is associated with left ventricular diastolic function and remodeling in patients with well-controlled type 2 diabetes, even after adjustment for possible confounding factors38. The mechanism underlying the correlation of LIPCAR with CAD remains unclear. It has been found that LIPCAR is strongly correlated with waist circumference, plasma fasting insulin, subcutaneous fat volume and HDL-C19. Our data also demonstrated that plasma level of LIPCAR was negatively associated with HDL-C. Thus, increased LIPCAR may induce metabolic dyshomeostasis, which in turn promotes atherosclerosis. Moreover, mitochondrial dysfunction is implicated in the etiology of cardiovascular diseases, especially CAD38. Thus, as a mitochondria-derived lncRNA, a potential function of LIPCAR in regulating mitochondrial pathways, such as oxidative phosphorylation and inflammasome activation, needs to be explored in future studies transcriptional.
Other limitations may also be addressed. Firstly, selection bias in the present study might have affected our results. Large-scale, multicenter studies are required to further elucidate the role of H19 and LIPCAR as a potential biomarker for the risk of CAD. Secondly, the underlining mechanisms of the association between up-regulated H19 and the severity of CAD need to be further studied. Thirdly, further experimental studies using animal or cellular model are needed to explore the mechanisms by which H19 and LIPCAR participate in the development of atherosclerosis.
In conclusion, our study showed that plasma levels of H19 and LIPCAR are associated with increased risk of CAD and may be considered as novel biomarkers for CAD.
This work was supported by grants from the National Natural Science Foundation of China (No. 81570363), the “333 Project” of Jiangsu Province (No. BRA2015326), a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the project of Nanjing Medical University (No. 2016NJMU060), and the Jiangsu Province Health Development Project with Science and Education (No. QNRC2016857).
Z.Z., G.W., and W.L.S. performed the data analyses and draft the manuscript. Z.Z., Z.J., and Y.Z.J. designed the experiments. L.Y.F., L.Q.Q. and L.D.C. collected the clinical data, which was supervised by W.L.S. G.W., L.Y.F. and Y.J.J. performed the experiments, which was supervised by Y.Z.J. All authors reviewed and approved the final manuscript.
The authors declare that they have no competing interests.
Zhen Zhang and Wei Gao contributed equally to this work.
Electronic supplementary material
Supplementary information accompanies this paper at doi:10.1038/s41598-017-07611-z
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.