BioMed Research International
Volume 2016 (2016), Article ID 4682172, 6 pages
http://dx.doi.org/10.1155/2016/4682172

Research Article

MicroRNA Expression Signature in Degenerative Aortic Stenosis

Jing Shi, Hui Liu, Hui Wang, and Xiangqing Kong

Department of Cardiology, The First Affiliated Hospital of Nanjing Medical University, Nanjing 210029, China

Received 4 April 2016; Accepted 28 June 2016

Academic Editor: Dragos Cretoiu

Copyright © 2016 Jing Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Degenerative aortic stenosis, characterized by narrowing of the exit of the left ventricle of the heart, has become the most common valvular heart disease in the elderly. The aim of this study was to investigate the microRNA (miRNA) signature in degenerative AS. Through microarray analysis, we identified the miRNA expression signature in the tissue samples from healthy individuals () and patients with degenerative AS (). Six miRNAs (hsa-miR-193a-3p, hsa-miR-29b-1-5p, hsa-miR-505-5p, hsa-miR-194-5p, hsa-miR-99b-3p, and hsa-miR-200b-3p) were overexpressed and 14 (hsa-miR-3663-3p, hsa-miR-513a-5p, hsa-miR-146b-5p, hsa-miR-1972, hsa-miR-718, hsa-miR-3138, hsa-miR-21-5p, hsa-miR-630, hsa-miR-575, hsa-miR-301a-3p, hsa-miR-636, hsa-miR-34a-3p, hsa-miR-21-3p, and hsa-miR-516a-5p) were downregulated in aortic tissue from AS patients. GeneSpring 13.1 was used to identify potential human miRNA target genes by comparing a 3-way comparison of predictions from TargetScan, PITA, and microRNAorg databases. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to identify potential pathways and functional annotations associated with AS. Twenty miRNAs were significantly differentially expressed between patients with AS samples and normal controls and identified potential miRNA targets and molecular pathways associated with this morbidity. This study describes the miRNA expression signature in degenerative AS and provides an improved understanding of the molecular pathobiology of this disease.

1. Introduction

Degenerative aortic stenosis (AS), the most common valvular heart disease in the elderly, has become a leading cause for surgical valve replacement in industrialized countries. As a result of rising life expectancy and ageing populations, the global prevalence of AS is increasing and is expected to surpass 10 billion by 2100 [1–3].

Degenerative AS has long been considered a chronic process with gradual deposition of calcium phosphate in the valve occurring with age. However, emerging evidence has indicated that this condition is mediated by the interplay of complex biological processes that include the following: inflammation, cell apoptosis, lipids deposition, renin-angiotensin system activation, remodeling of the extracellular matrix, and bone formation [4–11]. However, the underlying mechanisms that regulate this process remain largely unknown [12, 13].

MicroRNAs (miRNAs) are small, 21–25-nucleotide, endogenous, single-stranded noncoding RNAs that regulate target gene expression by binding messenger RNAs (mRNAs) and inhibiting or reducing translation. A single miRNA can regulate numerous genes, while a single gene can be regulated by multiple miRNAs [14]. miRNAs play a critical role in many physiological processes, and there is a growing body of studies indicating that distinct patterns of altered miRNA expression are associated with specific disease processes [15–19]. In the present study, we explore the miRNA expression signature of degenerative AS to improve our understanding of the molecular alternations in this disease.

2. Methods

2.1. Tissue Samples Collection and RNA Isolation

Written informed consent was obtained from all participants of age and/or via their parents. All procedures in this study were approved by the Ethics Committees of the First Affiliated Hospital of Nanjing Medical University and conformed to the principles outlined in the Declaration of Helsinki. Tissue samples from four healthy controls were obtained from prospective multiorgan donors without cardiovascular pathology in cases in which technical reasons prevented transplantation. Aortic valves from patients with degenerative AS were obtained from four patients who underwent surgical valve replacement. General characteristics of these participants are displayed in Table 1. Gross and histological examination was performed to confirm the presence/absence of AS in each sample. Valve samples were immediately snap-frozen in liquid nitrogen upon collection. Total RNA was extracted using an RNeasy Mini Kit in accordance with the manufacturer’s protocol (Qiagen, Hilden, Germany).

Table 1: Characteristics of the participants.

Control
AS
value

Number
4
4
ns

Male
4
4
ns

Age
41.5 ± 10.4
58.3 ± 4.6
0.0429

LVEF < 50%
0
0
ns

Coronary disease
0
0
ns

Diabetes mellitus
0
0
ns

Medications
0
0
ns

Values are displayed as mean ± standard deviation.

LVEF: left ventricular ejection fraction.

2.2. RNA Labeling and Array Hybridization

miRNA expression was evaluated using Agilent miRNA arrays (V2), which included 723 human and 76 human viral miRNAs from the Sanger database v.10.1 (Agilent Technologies, Foster City, CA). Total RNA was dephosphorylated and ligated with pCp-Cy3 and subsequently hybridized to the arrays. After samples were washed and scanned, using Agilent Scan Control software, the Agilent Feature Extract software v9.5.3 was used to analyze the arrays. Each sample was evaluated in triplicate. miRNA expression data were normalized using a bead-based assay and the locally weighted smooth spline (LOWESS) method. After normalization, all expression values were transformed to a linear scale for statistical analysis.

2.3. Bioinformatic Analysis

GeneSpring 13.1 was used to identify potential human miRNA target genes, which compared TargetScan, PITA, and microRNAorg databases and created a Venn diagram to demonstrate relations among the databases. Gene Ontology (GO) analysis was performed to investigate the biological processes, cellular components, and specific molecular function of differentially expressed coding genes identified. Pathway analysis was used to determine the involvement of coexpressed genes in different biological pathways according to Kyoto Encyclopedia of Genes and Genomes (KEGG).

2.4. Statistical Analysis

Independent Student’s -test was used to determine whether there were any significant differences between the miRNA expression profiles between two groups. values less than 0.05 () were considered to be statistically significant. Significant data were further analyzed by cluster analysis, and the expression profiles were visualized with GeneSpring 10.0 (Agilent Technologies).

3. Results

3.1. Unsupervised Hierarchical Cluster Analysis of miRNA Microarray Data

miRNA microarray identified 20 miRNAs with significantly differential expression (>2.0-fold) in AS samples relative to normal controls. Six miRNAs were overexpressed (hsa-miR-193a-3p, hsa-miR-29b-1-5p, hsa-miR-505-5p, hsa-miR-194-5p, hsa-miR-99b-3p, and hsa-miR-200b-3p) and 14 (hsa-miR-3663-3p, hsa-miR-513a-5p, hsa-miR-146b-5p, hsa-miR-1972, hsa-miR-718, hsa-miR-3138, hsa-miR-21-5p, hsa-miR-630, hsa-miR-575, hsa-miR-301a-3p, hsa-miR-636, hsa-miR-34a-3p, hsa-miR-21-3p, and hsa-miR-516a-5p) were downregulated in aortic tissue from AS patients (Table 2). Unsupervised hierarchic clustering was performed based on the 20 differentially expressed miRNAs and displayed as heat map (Figure 1).

Table 2: Six overexpressed and 14 downregulated miRNAs in aortic tissue from AS patients compared to the control group.

Systematic name
value
FC (abs)
Regulation

hsa-miR-193a-3p
0.020637
2.067236
Up

hsa-miR-29b-1-5p
0.0222
11.56757
Up

hsa-miR-505-5p
0.030211
13.22606
Up

hsa-miR-194-5p
0.036807
14.95685
Up

hsa-miR-99b-3p
0.029772
15.98001
Up

hsa-miR-200b-3p
0.021638
20.1747
Up

hsa-miR-3663-3p
0.042986
2.063909
Down

hsa-miR-513a-5p
0.01703
2.238355
Down

hsa-miR-146b-5p
0.040094
2.244855
Down

hsa-miR-1972
0.033015
2.432029
Down

hsa-miR-718
0.033461
2.559791
Down

hsa-miR-3138
0.002648
2.69498
Down

hsa-miR-21-5p
0.004216
3.334316
Down

hsa-miR-630
0.020492
3.564175
Down

hsa-miR-575
0.029949
5.137385
Down

hsa-miR-301a-3p
0.011849
20.77046
Down

hsa-miR-636
0.010177
25.1902
Down

hsa-miR-34a-3p
0.013993
25.44595
Down

hsa-miR-21-3p
0.02304
32.11932
Down

hsa-miR-516a-5p
0.022169
41.76603
Down

Figure 1: Unsupervised hierarchical clustering identified two distinct groups (control versus AS) based on their miRNA expression profile. Sample names are listed at the top. The names of the significantly () differentially expressed miRNAs are shown on the right. Twenty miRNAs were expressed differently in target genes analysis.

microRNAorg, TargetScan, and PITA were used to predict the targets of differentially expressed miRNAs in CAVD samples using software GeneSpring 13.1. A Venn diagram was generated to highlight the relations among the 3 databases. There are 1010 overlapping genes identified by all 3 programs, which are most likely to be targets of miRNAs in patients with AS (Figure 2).

Figure 2: The red, green, and blue sets stand for target genes predicted by databases microRNAorg, TargetScan, and PITA, respectively.

3.2. GO and Pathway Analysis

Gene Ontology (http://geneontology.org/) was used to classify the function of up- and downregulated genes from 3 structured networks: biological processes, cellular components, and molecular function.

In this study, differentially expressed mRNAs were enriched in numerous biological processes including the following: cell adhesion, homophilic cell adhesion, positive regulation of transcription, negative regulation of JAK-STAT cascade, and positive regulation of G1/S transition of mitotic cell cycle associated with biological processes (Figure 3(a)). Similarly, the following cellular components were affected: nucleus, nucleoplasm, and cytoplasm linked with cellular components (Figure 3(b)), while affected molecular functions include the following: protein binding, zinc ion binding, and transcriptional activator activity involved in molecular functions (Figure 3(c)). Moreover, KEGG pathway analysis identified significantly () affected pathways including the following: cAMP signaling pathway, vascular smooth muscle contraction, regulation of actin cytoskeleton, neurotrophin signaling pathway, and cGMP-PKG signaling pathway (Figure 4).

(a)

(b)

(c)

Figure 3: GO analysis for differentially expressed mRNAs. (a)–(c) GO analysis according to biological process, cellular component, and molecular function, respectively, ranked by enrichment score ( ( value)).

Figure 4: Pathway analysis based on the KEGG database. Ranked by enrichment score ( ( value)).

4. Discussion

Recent studies have improved our understanding of the mechanisms underlying AS [12, 20]. However, there are not many reports that have investigated the function of miRNAs as they relate to the pathobiology of AS [21–26].

A microRNA expression signature provides a better understanding of the mechanisms of a disease [27]. In the present study, we explored miRNA expression signatures associated with degenerative AS using miRNA microarray analysis. Six overexpressed miRNAs (hsa-miR-193a-3p, hsa-miR-29b-1-5p, hsa-miR-505-5p, hsa-miR-194-5p, hsa-miR-99b-3p, and hsa-miR-200b-3p) and 14 downregulated miRNAs (hsa-miR-3663-3p, hsa-miR-513a-5p, hsa-miR-146b-5p, hsa-miR-1972, hsa-miR-718, hsa-miR-3138, hsa-miR-21-5p, hsa-miR-630, hsa-miR-575, hsa-miR-301a-3p, hsa-miR-636, hsa-miR-34a-3p, hsa-miR-21-3p, and hsa-miR-516a-5p) were identified in patients with AS, relative to normal controls, and their general characteristics and functional annotations were analyzed using bioinformatic tools.

There were other miRNA microarrays performed about AS before, and, in these studies, the main cause of AS is calcification of the aortic valves [20, 22]. However, as rheumatic fever remains to be the most important etiological factor in calcific aortic valve disease (CAVD) in China, we do not choose CAVD patients to reduce confounding factors when exploring the mechanisms that underlie degenerative AS.

As a limitation of our study, the specific pathways affected by miRNAs and which cause AS remain elusive. Further studies are required to functionally characterize the role of specific candidate miRNAs.

In conclusion, the present study improved our understanding of the role of miRNAs in degenerative AS. Our findings provide improved understanding of the molecular alterations in this disease and may provide potential targets for future clinical applications.

Competing Interests

The authors declare that there are no competing interests.

Authors’ Contributions

Jing Shi and Hui Liu contributed equally to this work.

Acknowledgments

Dr. Xiangqing Kong is a fellow at the Collaborative Innovation Center for Cardiovascular Disease Translational Medicine and this work was supported by the grants from the center.

References

J. J. Thaden, V. T. Nkomo, and M. Enriquez-Sarano, “The global burden of aortic stenosis,” Progress in Cardiovascular Diseases, vol. 56, no. 6, pp. 565–571, 2014.

N. Rezzoug, B. Vaes, A. Pasquet et al., “Prevalence and prognostic impact of valve area-gradient patterns in patients years with moderate-to-severe aortic stenosis (from the Prospective BELFRAIL Study),” American Journal of Cardiology, vol. 116, no. 6, Article ID 21248, pp. 925–932, 2015.

B. Bordoni, F. Saia, C. Ciuca et al., “Prevalence of degenerative aortic valve stenosis in the elderly: results of a large community-based epidemiological study,” Giornale Italiano di Cardiologia, vol. 14, no. 4, pp. 262–268, 2013.

K. Arishiro, M. Hoshiga, N. Negoro et al., “Angiotensin receptor-1 blocker inhibits atherosclerotic changes and endothelial disruption of the aortic valve in hypercholesterolemic rabbits,” Journal of the American College of Cardiology, vol. 49, no. 13, pp. 1482–1489, 2007.

F. C. Caira, S. R. Stock, T. G. Gleason et al., “Human degenerative valve disease is associated with up-regulation of low-density lipoprotein receptor-related protein 5 receptor-mediated bone formation,” Journal of the American College of Cardiology, vol. 47, no. 8, pp. 1707–1712, 2006.

A. Galeone, D. Paparella, S. Colucci, M. Grano, and G. Brunetti, “The role of TNF-α and TNF superfamily members in the pathogenesis of calcific aortic valvular disease,” The Scientific World Journal, vol. 2013, Article ID 875363, 10 pages, 2013.

N. Latif, P. Sarathchandra, A. H. Chester, and M. H. Yacoub, “Expression of smooth muscle cell markers and co-activators in calcified aortic valves,” European Heart Journal, vol. 36, no. 21, pp. 1335–1345, 2015.

L. Luna-Zurita, B. Prados, J. Grego-Bessa et al., “Integration of a Notch-dependent mesenchymal gene program and Bmp2-driven cell invasiveness regulates murine cardiac valve formation,” The Journal of Clinical Investigation, vol. 120, no. 10, pp. 3493–3507, 2010.

A. Alexopoulos, A. Kaoukis, H. Papadaki, and V. Pyrgakis, “Pathophysiologic mechanisms of calcific aortic stenosis,” Therapeutic Advances in Cardiovascular Disease, vol. 6, no. 2, pp. 71–80, 2012.

T. A. Pawade, D. E. Newby, and M. R. Dweck, “Calcification in aortic stenosis: the skeleton key,” Journal of the American College of Cardiology, vol. 66, no. 5, pp. 561–577, 2015.

D. A. Lerman, S. Prasad, and N. Alotti, “Calcific aortic valve disease: molecular mechanisms and therapeutic approaches,” European Cardiology Review, vol. 10, no. 2, pp. 108–112, 2015.

K. E. Yutzey, L. L. Demer, S. C. Body et al., “Calcific aortic valve disease: a consensus summary from the alliance of investigators on calcific aortic valve disease,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 34, no. 11, pp. 2387–2393, 2014.

I. Sathyamurthy and S. Alex, “Calcific aortic valve disease: is it another face of atherosclerosis?” Indian Heart Journal, vol. 67, no. 5, pp. 503–506, 2015.

L. A. Yates, C. J. Norbury, and R. J. C. Gilbert, “The long and short of microRNA,” Cell, vol. 153, no. 3, pp. 516–519, 2013.

T. Thum and G. Condorelli, “Long noncoding RNAs and microRNAs in cardiovascular pathophysiology,” Circulation Research, vol. 116, no. 4, pp. 751–762, 2015.

X. Liu, J. Xiao, H. Zhu et al., “miR-222 is necessary for exercise-induced cardiac growth and protects against pathological cardiac remodeling,” Cell Metabolism, vol. 21, no. 4, pp. 584–595, 2015.

Y. F. Melman, R. Shah, K. Danielson et al., “Circulating MicroRNA-30d is associated with response to cardiac resynchronization therapy in heart failure and regulates cardiomyocyte apoptosis a translational pilot study,” Circulation, vol. 131, no. 25, pp. 2202–2216, 2015.

H. Wang, Y. Bei, S. Shen et al., “miR-21-3p controls sepsis-associated cardiac dysfunction via regulating SORBS2,” Journal of Molecular and Cellular Cardiology, vol. 94, pp. 43–53, 2016.

D. Lv, J. Liu, C. Zhao et al., “Targeting microRNAs in pathological hypertrophy and cardiac failure,” Mini-Reviews in Medicinal Chemistry, vol. 15, no. 6, pp. 475–478, 2015.

N. M. Rajamannan, F. J. Evans, E. Aikawa et al., “Calcific aortic valve disease: not simply a degenerative process: a review and agenda for research from the National Heart and Lung and Blood Institute Aortic Stenosis Working Group. Executive summary: calcific aortic valve disease—2011 update,” Circulation, vol. 124, no. 16, pp. 1783–1791, 2011.

B. Yanagawa, F. Lovren, Y. Pan et al., “miRNA-141 is a novel regulator of BMP-2-mediated calcification in aortic stenosis,” Journal of Thoracic and Cardiovascular Surgery, vol. 144, no. 1, pp. 256–e2, 2012.

M. Zhang, X. Liu, X. Zhang et al., “MicroRNA-30b is a multifunctional regulator of aortic valve interstitial cells,” Journal of Thoracic and Cardiovascular Surgery, vol. 147, no. 3, pp. 1073–1080.e2, 2014.

P. Ohukainen, S. Syvaranta, J. Napankangas et al., “MicroRNA-125b and chemokine CCL4 expression are associated with calcific aortic valve disease,” Annals of Medicine, vol. 47, no. 5, pp. 423–429, 2015.

J. Rysä, “Gene expression profiling of human calcific aortic valve disease,” Genomics Data, vol. 7, pp. 107–108, 2016.

J. F. Esmerats, J. Heath, and H. Jo, “Shear-sensitive genes in aortic valve endothelium,” Antioxidants & Redox Signaling, 2016.

S. Coffey, M. J. Williams, L. V. Phillips, and G. T. Jones, “Circulating microRNA profiling needs further refinement before clinical use in patients with aortic stenosis,” Journal of the American Heart Association, vol. 4, no. 8, Article ID e002150, 2015.

J. Xiao, D. Liang, Y. Zhang et al., “MicroRNA expression signature in atrial fibrillation with mitral stenosis,” Physiological Genomics, vol. 43, no. 11, pp. 655–664, 2011.