Skip to main content

Specific circulating microRNA signature of bicuspid aortic valve disease

Abstract

Background

We aimed to determine the circulating miRNA expression profile associated with BAV and aortic dilation to provide diagnostic and prognostic biomarkers for BAV and/or aortic dilation.

Methods and results

We applied a miRNome-wide microarray approach using plasma samples (n = 24) from healthy tricuspid aortic valve individuals, BAV patients and BAV patients with aortic dilation to compare and identify the specific miRNAs associated with BAV and aortic dilation. In a second stage, the expression patterns of the miRNA candidates were validated by RT-qPCR in an independent cohort (n = 43). The miRNA microarray data and RT-qPCR analyses revealed that the expression levels of circulating miR-122, miR-130a and miR-486 are significantly influenced by the morphology of the aortic valve (bicuspid/tricuspid) and could be functionally involved in the regulation of TGF-β1 signalling. Furthermore, the expression pattern of miR-718 in the plasma was strongly influenced by dilation of the ascending aorta. miR-718 expression was inversely correlated with the aortic diameter (R = −0.63, p = 3.1 × 10−5) and was an independent predictor of aortic dilation (β = −0.41, p = 0.022). The genes targeted by miR-718 are involved in the regulation of vascular remodelling.

Conclusions

We propose that miR-122, miR-130a, miR-486 and miR-718 are new molecular features associated with BAV and aortic dilation principally by the activation of TGF-β1 pathway and vascular remodelling mediated by VEGF signalling pathways.

Background

Bicuspid aortic valve (BAV) represents the most common congenital cardiac malformation and is commonly associated with the development of aortic valve dysfunction and the progressive dilation of the ascending aorta. In patients with BAV, the dilation of the aorta is associated with aortic regurgitation [1] and the risk of aortic dissection and thus often requires prophylactic aortic surgery [26]. The causes of aorta dilation have been debated for several years, and changes in the flow characteristics within the ascending aorta have been suggested as the main cause [7, 8]. However, the mechanisms involved have not been fully elucidated.

MicroRNAs (miRNAs) are endogenously expressed 19–23nt-long noncoding RNAs (ncRNAs) that silence gene expression at the post-transcriptional level, largely via imperfect base-pairing interactions that occur preferentially within the 3′ untranslated regions (UTRs) of target mRNAs. Due to the ability of miRNAs to target hundreds of mRNAs, they are considered as potent post-transcriptional regulators, and in fact, it is suggested that a significant portion of the human genome is regulated by miRNAs. Aside from their intracellular function, miRNAs can be exported or released by cells into the circulating blood in highly stable forms [9]. miRNA signatures have been proposed as potential biomarkers with the potential to improve disease diagnosis and prognosis in clinical practice and have been identified as useful biomarkers for a wide range of cardiovascular diseases, including cardiovascular calcification [10], vulnerable coronary artery disease [11], heart failure [12], myocardial infarction [13], abdominal aortic aneurysm [14] and coronary artery disease [15]. Despite the large number of miRNAs identified thus far, the biological roles of most miRNAs and the molecular mechanisms underlying their function remain largely unknown.

In this study, we applied a miRNome-wide microarray approach using the plasma of healthy tricuspid aortic valve (TAV) individuals and BAV patients [with or without aortic dilation (BAVdil)] to compare and identify specific miRNAs associated with BAV and the dilation of the ascending aorta. In a second stage, the expression levels of the miRNA candidates were validated in an independent cohort, providing new insights into the complex mechanisms that underlie BAV and aortic dilation.

Methods

Study population

The participants included in this study belong to a cohort of prospective patients with BAV that were followed-up in our facilities. Upon enrolment, the participants were prospectively entered into a specific database, underwent a blood draw and provided informed written consent. The samples were stored until needed in our biological sample bank (Biobanc IISPV-HUSJR). BAV diagnosis was made when two aortic leaflets were clearly visualized, with or without a raphe, either on the parasternal short-axis view of a transthoracic echocardiogram [16], on a transesophageal echocardiogram [16] or on a cardiac magnetic resonance [17]. Explorations were performed or supervised by the same observer (JMA). This study was approved by the Institutional Review Board (the Clinical Ethics Committee) of our institution.

The study design was composed of two evaluations to determine the possible miRNAs associated with BAV and aortic dilation. First, we studied the association of BAV with the modulation of the circulating miRNOme profile by determining the miRNA expression profile in the plasma of patients diagnosed with BAV (BAVnon-dil; n = 9) or BAV and aortic dilation (BAVdil; n = 9) using microarrays, and then, we compared these with the expression profiles of healthy TAV control subjects (TAVnon-dil; n = 6). For this stage and to improve the study power, we selected three groups composed of patients with characteristics that were extremely homogeneous (Table 1) to exclude possible confounding factors, such as age, sex or BMI. In addition, patients with aortic stenosis, hypertension, diabetes mellitus or pharmacologic treatment (including statins, ACE/ARII and/or β-blockers) were excluded. An aortic diameter lower than 19 mm/m2 or higher than 21 mm/m2 was used as the criterion to differentiate between the non-dilated and dilated patients, respectively.

Table 1 Clinical and echocardiographic characteristics of the patients included in the microarray

In the second evaluation, we validated the miRNA candidates by RT-qPCR in a new cohort of individuals comprising TAV (n = 19) and BAV individuals (n = 24) (Table 2). The BAV group comprised BAVnon-dil (n = 12) and BAVdil (n = 12) patients. The TAV group comprised healthy controls TAVnon-dil (n = 12), and we also included a new group of patients with TAV and aortic dilation (TAVdil; n = 7) in this evaluation to determine the specific effects associated with valve morphology or aortic dilation per se. Patients with Marfan syndrome were excluded.

Table 2 Clinical and echocardiographic characteristics of the independent cohort used in the validation study

Blood sampling

Blood samples were collected in overnight fasting conditions in the early morning and were processed within 90 min after collection. The samples were centrifuged at 1500g for 15 min to obtain plasma. The plasma samples were stored at −80 °C in our biological sample bank until needed (Biobanc IISPV-HUSJR).

RNA extraction and miRNA microarray preparation

Total RNA was extracted from 250 μL of plasma using the TRIzol reagent according to the manufacturer’s instructions (Invitrogen) and purified using an RNeasy minikit (Qiagen). To increase RNA recovery, 1 μg of MS2 carrier RNA was added to each plasma sample. The quality of the total RNA isolated was determined using an Agilent 2100 Bioanalyser.

The plasma miRNA expression levels for 24 individuals were assessed using SurePrint G3 human 8 × 60 k miRNA microarrays (Agilent Technologies), covering 1205 human miRNAs (Sanger miRBase release 16). The miRNAs were dephosphorylated and labelled with cyanine 3-cytidine biphosphate, including a labelling spike-in solution (Agilent Technologies) to assess the labelling efficiency. The samples were hybridized on the arrays, including a hybridization spike-in solution (Agilent Technologies) to monitor the hybridization efficiency. The arrays were scanned with a G2565CA Microarray Scanner System with SureScan High Resolution Technology (Agilent Technologies) using Scan Control software. The Feature Extraction 11.5.11 (Agilent Technologies) and GeneSpring 12.6.1 software packages were used for data processing, and all steps were performed according to the manufacturer’s protocol.

Microarray analysis of the miRNA expression data

The data from microarrays were normalized using the robust multi-array average (RMA) method [18], implemented in the AgiMicroRna Bioconductor package, and the fold changes of the circulating miRNAs were determined using the linear model implemented in the limma Bioconductor package [19]. The Benjamini and Hochberg method was used to adjust the p values for multiple testing and to control the false discovery rate [20].

MicroRNA quantification by real-time qRT-PCR

The microarray expression data were validated in an independent patient cohort (n = 43), including healthy TAV controls, BAV patients, BAVdil patients, and an additional group of TAV individuals with aortic dilation. Reverse transcription was performed using the TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems) and the miRNA-specific reverse-transcription primers provided with the TaqMan MicroRNA Assay (Applied Biosystems) to analyse the expression of each miRNA. The final total RNA concentration used was 2.5 ng/µL. The reaction was performed at 16 °C for 30 min, 42 °C for 30 min and 85 °C for 5 min. We used 1.33 µL of the obtained cDNAs in a subsequent quantitative qRT-PCR amplification using the TaqMan Universal PCR master mix (Applied Biosystems) and the associated specific probe provided in the TaqMan MicroRNA Assay Kit (Applied Biosystems). Specific Taqman probes were used for each gene as follows: microRNA-122-5p (miR-122-5p; hsa-mir-122-5p), 5′-UGGAGUGUGACAAUGGUGUUUG-3′; microRNA-130a-3p (miR-130a-3p; hsa-mir-130a-3p), 5′-CAGUGCAAUGUUAAAAGGGCAU-3′; microRNA-486-5p (miR-486-5p; hsa-mir-486-5p), 5′-UCCUGUACUGAGCUGCCCCGAG-3′ and microRNA-718 (miR-718; hsa-mir-718), 5′-CUUCCGCCCCGCCGGGCGUCG-3′. The results were normalized to the expression of the U6 small nuclear RNA (U6 snRNA), which was used as an endogenous control. Amplification was performed using the ABI Prism 7300 SDS Real-Time PCR system (Applied Biosystems, Madrid, Spain) at 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. The fold change in the miRNA level was calculated by the log 2 scale according to the equation 2−ΔΔCt, where ΔCt = Ct miRNA-Ct U6 and ΔΔCt = ΔCt treated samples − ΔCt untreated controls.

Functional and pathway enrichment analysis of the miRNA target genes

The target genes of miR-122, miR-130a and miR-486 were predicted using DIANA-miRPath (v3; microrna.gr/tarbase) [21], mapped on the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways based on an enrichment analysis, and prioritized for the intracellular signal transduction pathways that could be regulated by the three-miRNAs.

Because of the limited target genes identified for miR-718, for this miRNA, we combined the targets genes predicted by four databases, namely, the TargetScan (v7.0; targetscan.org) [22], miRTarBase (2016 version; mirtarbase.mbc.nctu.edu.tw) [23], DIANA-microT-CDS (v3; microrna.gr/microT-CDS) [21] and miRDB (2016 version; mirdb.org) [24] databases. In this case, the functional classification of the target genes was carried out with Gene Ontology (GO) and pathway analysis using DAVID (david.ncifcrf.gov) [25] and WebGestalt (webgestalt.org) [26], respectively.

Cytoscape version 3.0.1 [27] was used to visualize the predicted miRNA-mRNA interaction. In the network, the intracellular signal transduction pathways are represented by hexagons-shaped nodes, the miRNAs are represented by rectangle-shaped nodes and the target genes are represented by round-shaped nodes.

Statistical analysis

The means of two groups were compared using Student’s t test. A Chi squared test or a Fisher’s exact test, when appropriate, was used to compare the frequencies of the categorical variables. The effects of the valve morphology and aortic dilation were assessed using analysis of variance (ANOVA). Tukey’s test was used for the pairwise comparisons. A Pearson’s correlation was used to identify linear relationships between the miRNA expression levels and the diameter of the ascending aorta. Backward linear regression models were used to identify the independent predictors of the expression of the circulating miRNAs. p values <0.05 were considered significant. The statistical analyses were performed using SPSS software version 21.0 (IBM, Chicago, IL, USA).

Results

Identification of the circulating miRNome profiles

Global miRNA expression profiling was performed on the plasma samples of the TAV controls and the BAVnon-dil and BAVdil patients to identify the miRNA candidates, and the expression of these miRNAs was altered in BAV because of the valve morphology and/or because of the aortic dilation. The clinical characteristics of the patients involved in this stage are provided in Table 1. After processing the microarray data corresponding to the miRNA expression profiles, 260 miRNAs were identified as expressed in at least 10% of the samples. We performed a paired comparison between the groups to identify miRNA candidates for not only BAV but also aortic dilation. In this manner, we first analysed the effect of the valve morphology (comparing patients with TAVnon-dil and BAVnon-dil) on the miRNome, and as shown in the volcano plot (Fig. 1a), the expression levels of miR-122 and miR-486 differed significantly from the remainder, indicating that the expression levels of these two miRNAs in plasma might be associated with the morphology of the aortic valve.

Fig. 1
figure 1

Volcano plot representing differences in the expression levels of the circulating miRNAs depending on a the bicuspid or tricuspid morphology of the aortic valve and b the morphology of the valve together with the dilation of the aorta. The log-fold change is plotted against the −log10 (p value)

In a second analysis, we compared the two most divergent groups (TAVnon-dil and BAVdil) to integrate the complexity associated with the valve morphology and the dilation of the ascending aorta. This analysis indicated that the expression of two other miRNAs, miR-130a and miR-718 (Fig. 1b), was significantly different between the plasma of the TAVnon-dil and BAVdil patients.

RT-qPCR validation of the miRNA candidates in an independent cohort

The miRNA candidates selected using the microarray approach were validated by RT-qPCR in a new set (n = 43), which comprised TAVnon-dil healthy controls and BAVnon-dil and BAVdil patients, and we also included a new group of patients with TAV and aortic dilation (TAVdil) to determine the specific pattern associated with the valve morphology or aortic dilation per se.

The expression levels of miR-122, miR-130a and miR-486 are associated with the morphology of the aortic valve

We compared the plasma expression of the miRNA candidates, miR-122, miR-130a, miR-486 and miR-718, which depended on the morphology of the valve and was independent of the presence or absence of dilation in the ascending aorta (Fig. 2a).

Fig. 2
figure 2

The expression levels of miR-122, miR-130a and miR-486 are influenced by the morphology of the aortic valve. a The validation of the miRNA candidates for aortic morphology was performed by RT-qPCR. Two asterisks Significant values (p < 0.001; unpaired t test). b The pathways of intracellular signal transduction were significantly enriched by the predicted target genes of miR-122, miR-130a and miR-486. c Interactive miRNA-gene network regarding the TGF-beta signalling pathway. TAV tricuspid aortic valve, BAV bicuspid aortic valve, TGF transforming growth factor

This validation study corroborated that the expression of three of the miRNA candidates, miR-122, miR-130a, and miR-486, differed significantly between the plasma of TAV individuals and that of BAV patients (BAV vs. TAV), and therefore, their expression levels in the circulation were influenced by the morphology of the valve. The bicuspid valve morphology was associated with decreased circulating levels of miR-122 and miR-486 instead of the increased levels reported for miR-130a.

We performed a functional categorization of the miRNA target genes to gain further insight into the biological role of miR-122, miR-130a and miR-486 in BAV (Fig. 2b). The targets of each of the selected miRNAs were predicted, mapped on KEGG pathways based on an enrichment analysis, and used to prioritize the intracellular signal transduction pathways that could be regulated by the three miRNAs. Using this approach, we identified the transforming growth factor-beta (TGF-β1) signalling pathway as the most significantly enriched pathway (p = 7.73 × 10−7), with a total of 32 genes targeted by the four selected miRNAs (Fig. 2c).

Effect of aortic dilation on the miR-130a and miR-718 expression levels in plasma

Next, we validated whether the modulation of the miRNA expression profile of miR-130a and miR-718 in the plasma was related to the dilation of the aorta. We first analysed the effect of the aortic dilation on miR-130a and miR-718, comparing the expression levels between the dilated and non-dilated patients independent of the aortic valve morphology (n = 43; Fig. 3A). The up-regulation of miR-130a expression in the plasma of the BAV patients was only influenced by the morphology of the aortic valve, and the dilation of the aorta had no modulatory effect. In contrast, the plasma expression of miR-718 was significantly affected by the dilation of the ascending aorta, and this dilation-promoted effect was masking the differential expression consequence of the valve morphology not observed in the comparisons used in the previous validation.

Fig. 3
figure 3

The expression of circulating miR-718 is affected by the dilation of the ascending aorta. A RT-qPCR validation of the effect of the ascending aortic dilation on the expression levels of miR-130a and miR-718 in the plasma. The specific effects of the valve morphology and aortic dilation were analysed for miR-130a (B) and miR-718 (C). Two asterisks Significant values (p < 0.001; unpaired t test). The bars with different letters are significantly different (one-way ANOVA, Tukey test, p < 0.05)

We explored these morphology- and aortic-dilation-mediated effects by comparing the expression of miR-130a (Fig. 3B) and miR-718 (Fig. 3C) between the four sub-groups depending on the valve morphology and the dilation of the ascending aorta (TAVnon-dil, TAVdil, BAVnon-dil and BAVdil). This analysis corroborated that up-regulated miR-130a expression was associated only with the bicuspid morphology of the aortic valve, and no effect of the aortic dilation was observed. Furthermore, the down-regulated expression of miR-718 observed in BAV patients was significantly more prevalent in the presence of aortic dilation compared with the expression in the non-aortic dilated patients and was independent of the tricuspid or bicuspid aortic valve morphology.

The expression of miR-718 in plasma is a predictor of aortic dilation

The role of miR-718 as a biomarker of aortic dilation was supported by the inverse correlation between miR-718 plasma expression and the diameter of the ascending aorta (R = −0.629, p = 3.1 × 10−5; Fig. 4a) in all of the patients in the validation cohort; this correlation was independent of valve morphology (TAV and BAV) or the presence or absence of aortic dilation (dilated and non-dilated). Furthermore, using this validation cohort, we built a multiple logistic and linear models using the miR-718 expression profiles as the predictive variables and the dilation of the aorta (represented as a categorical variable or as the diameter of the aorta) as the outcome variable in addition to age and hypertension as possible confounding factors (Table 3). This analysis confirmed that the expression of miR-718 in the plasma is an independent predictor of the dilation of the aorta.

Fig. 4
figure 4

Functional implications of miR-718 expression in ascending aortic dilation. a The expression of circulating miR-718 is inversely correlated with the diameter of the ascending aorta. b miR-718 target gene network resulting from the functional enrichment analysis

Table 3 Logistic and linear regression analyses of the miR-718 expression levels as a predictor of aortic dilation

To provide a better understanding of the possible functional relevance of miR-718 down-regulation, we characterized the putative biological functions of this miRNA based on the identification of the miR-718 target genes. The predicted miRNA-gene interaction network (Fig. 4b) was composed of miR-718 and the 167 target genes. We also performed a functional enrichment analysis to predict the potential biological function of the miR-718 target genes. DAVID was used to determine the GOs associated with biological processes that might be overrepresented within the miR-718 target genes. Four of the target genes were significantly enriched for the GO biological process “blood vessel remodelling” (p = 1.9 × 10−2). Using WebGestalt, we performed a functional analysis based on the pathways that might be associated with the miR-718 targets; 3 of the miR-718 target genes were significantly enriched for the “focal adhesion pathway” (p = 9.0 × 10−4).

Discussion

This study found that the miRNA expression pattern in plasma is influenced by the morphology of the aortic valve and the dilation of the ascending aorta. Using comprehensive miRNA expression profiling and further RT-qPCR validation, we identified a circulating miRNA signature, including miR-122, miR-130a and miR-486, whose regulation is dependent on the tricuspid or bicuspid morphology of the aortic valve. Furthermore, the plasma expression of miR-718 may be a biomarker of aortic dilation, that is, an increased ascending aorta diameter is associated with the down-regulation of miR-718. We explored the biological implications of the differential expression of these miRNAs based on target gene prediction and the putative biological functions and processes regulated by the target genes. In this manner, the differential modulation of the bicuspid valve morphology-associated miRNAs might affect the TGF-β1 signalling pathway, and the dysregulation of miR-718 might be associated with the progression of aortic dilation by modulating the focal adhesion and blood vessel remodelling processes.

The molecular features and mechanisms that underlie BAV and the dilation of the aorta are poorly understood, and to our knowledge, the miRNA expression profile and putative gene regulation in BAV disease and aortic dilation identified in this study by screening in plasma of BAV patients has not been previously reported. Previous studies have analysed the expression of miRNAs in aortic tissue segments of patients with BAV and animal models. Nigam et al. [28] determined that miR-26a, miR-30b and miR-195 were decreased in the aortic valve leaflets of patients requiring aortic valve replacement due to aortic stenosis compared to those requiring replacement due to aortic insufficiency, and both pathologies are associated with BAV. Using porcine valvular interstitial cells and human valve leaflets, Yanagawa et al. [29] determined that the down-regulation of miR-141 in stenotic bicuspid leaflets regulates BMP-2-mediated valve calcification. Wu et al. [30] performed a paired comparison of the miRNA expression between severely dilated and normal-appearing, less-dilated aortic samples, associating the differential regulation of the miR-17 gene cluster with the predisposition to dilation through the dysregulation of the matrix metalloproteinases-tissue inhibitors of matrix metalloproteinases (TIMP-MMP) pathway. However, although the determination of the miRNA expression profile in tissue samples provides valuable information regarding the pathophysiological mechanisms underlying diseases, this expression signature often cannot be extrapolated to the bloodstream. Thus, unlike these tissue-specific miRNAs, the determination of a circulating miRNA expression profile allows for the integration of the complexity associated with multiple tissues, as is in the case of complex diseases, in addition to serving as biomarkers of these diseases.

In clinical practice, the diagnosis of BAV and ascending aorta dilation is currently performed exclusively based on imaging techniques (i.e., echocardiography, TAC, and magnetic resonance) and the subsequent interpretation by an expert. Moreover, both BAV and the dilation of the aorta are normally presented as asymptomatic clinical conditions, and the molecular causes underlying both of these pathologies are unknown. Therefore, determining plasma molecular biomarkers of BAV and the progression of aortic dilation would be extremely useful in expanding the basic knowledge of the aortic dysfunction, as well as for the identification of new therapeutic targets that might significantly improve the diagnosis, prognosis and treatment of BAV and aortic dilation. Moreover, the determination of biomarkers in the plasma has an added value because in the plasma, the miRNA levels are reproducible and consistent among individuals, and the miRNAs appear to be protected from endogenous ribonuclease-induced degradation [31]. In addition, in contrast to other tissues, the plasma is more easily accessible by non-invasive means, which might improve the feasibility of these molecules as useful biomarkers in routine clinical practice.

The bicuspid morphology of the aortic valve is associated with an abnormal blood flow and shear stress at the ascending aorta that leads to damage and apoptosis of endothelial cells at the cellular level, involving the activation of multiple biological processes, such as endothelial dysfunction, inflammation, oxidative stress, angiogenesis and extracellular matrix remodelling [3234]. Considering this BAV-associated environment, we hypothesized that circulating miRNAs might be involved in the paracrine communication between cells and contribute to the regulation of disease progression, and thus, the circulating levels of miRNAs might reflect the hemodynamic alteration of the blood flow resulting from the bicuspid morphology of the aortic valve.

In addition to their role as biomarkers for various diseases, miRNAs can be considered potential dynamic master regulators of signalling pathways due to their capacity to specifically control the expression of key components of signal transduction at multiple levels. Therefore, based on a functional analysis, our data supported the crucial role of the dysregulation of genes involved in the TGF-β1 signalling pathway in BAV disease [35]. In addition, we propose that miR-122, miR-130a and miR-486 might serve as molecular effectors in the BAV associated dysregulation of TGF-β1 receptor. TGF-β1 is a crucial factor mediating tissue fibrosis and extracellular matrix (ECM) remodelling in ascending aorta. TGF-β1 is a pleiotropic cytokine secreted and stored in the ECM as a biologically inactive form that is complexed with latent TGF-β1 binding protein (LTBP-1) in the ECM. Activated TGF-β1 binds to a heterodimeric receptor complex comprising two subunits (TGFβR1 and TGFβR2) whose genes are targeted by miR-122 and miR-130a, respectively, activating Smad proteins, which mediates the intracellular signalling of TGF-β1 [36]. The phosphorylation of SMAD2 (the gene targeted by miR-122) and SMAD3 results in the formation of heterooligomeric complexes with SMAD4 (the gene targeted by miR-122, miR-130a and miR-486). TGF-β1 signalling has been reported to play divergent roles in aortic dilation depending on its location (i.e. abdominal or thoracic aortic aneurysms) or aortic valve morphology [37]. As limitation of this study, the role of these miRNAs in the regulation TGF-β1 signalling pathway was determined based on computational predictions and it was not experimentally confirmed.

We can speculate that the opposite regulation of mir-122 (decreased expression levels) and miR-130a (increased expression levels) that we observed could be associated with the imbalance of TGFβR1 and TGFβR2 subunits. This speculation is supported by previous findings where the postnatal disruption of the Tgfbr2 gene expression was associated with decreased canonical TGF-β1 signalling and accelerated aneurysm growth in a murine model of Marfan syndrome [38]. Moreover, Forte et al. [39], reported that ascending aorta aneurysms in BAV patients was associated with an imbalance between the TGFβR1 and TGFβR2 subunits, although, in that case the authors associated this imbalance by an intrinsic defect of TGFβR1 expression. This imbalance between both subunits of the TGF-β1 receptor could be indicating a redirection towards the activation of the non-canonical TGF-β1 signalling [37], which is related to the induction of the ECM degradation and MMP activity, both biological processes associated with BAV and aortic dilation.

In addition to the biological and cardiology-relevant implication of each miRNA individually, miR-122 accounts for 70% of the total miRNAs in the adult liver [40] and is considered a central player in liver biology and in the regulation of cholesterol and fatty acid metabolism [41]. Interestingly, Beaumont et al. [42] demonstrated that the down-regulation of miR-122 was associated with the up-regulation of TGF-β1 expression, increasing the severity of myocardial fibrosis in patients with aortic stenosis, which is the most frequent valvular complication in BAV patients [43]. miR-130a has important roles in the cardiovascular system due to its capacity to promote the proliferation of endothelial cells and vascular smooth muscular cells. In fact, miR-130a has been described as a regulator of the angiogenic phenotype of vascular endothelial cells through its ability to modulate the expression of the homebox genes GAX and HOXA5 [44]. Furthermore, the expression of miR-130a in the aorta is correlated with vascular remodelling in spontaneously hypertensive rats (SHRs) [45]. Regarding miR-486, Holliday et al. [46] demonstrated that the expression of miR-486 in human aortic valvular endothelial cells was shear dependent, and Navon et al. [47] suggested that a reduction of miR-486 supports tissue remodelling characterized by increased proliferation. Therefore, our findings regarding the miRNA signature associated with the morphology of the aortic valve support the notion that miR-122, miR-130a and miR-486 play a key role in BAV disease and add new knowledge regarding the molecular effectors of the TGF-β 1 pathway dysregulation associated with BAV.

We propose miR-718 as a new non-invasive biomarker of dilation of the ascending aorta. We demonstrated its potential as molecular feature associated with an aortic dilation in patients with a dilated ascending aorta that was independent of the tricuspid or bicuspid morphology of the aortic valve. We found only 4 references [4851] in PubMed when searching in the literature for information on the biological function of miR-718; to date, no one has referred to the biological function of miR-718 expression in the cardiovascular field, including its function in relation to aortic valve morphology or the dilation of the aorta.

Using a bioinformatics approach, we predicted the putative target genes of miR-718 and the biological function and processes that it might regulate. The miR-718 target genes were significantly enriched for the “blood vessel remodelling” process and “focal adhesion pathway”. Ascending aorta dilation is defined as a degenerative vascular disorder due to the destructive remodelling of the aortic wall and the degradation of the ECM proteins, leading to the recruitment and infiltration of immune system cells mediated by the secretion of adhesion molecules. Interestingly, vascular endothelial growth factor A (VEGFA) is one of the miR-718 target genes that was significantly enriched for the “focal adhesion pathway”. VEGFA acts on endothelial cells by increasing vascular permeability and angiogenesis, inducing vasculogenesis and endothelial cell growth, promoting cell migration, and inhibiting apoptosis [52]. These biological implications of the miR-718 target genes might explain the relation between miR-718 and ascending aorta dilation.

In BAV, as in other disease-associated complications, the emergence of the ascending aorta dilation cannot be considered as a discrete alteration in time; instead, it is a continuous process associated with the pathophysiological implications of BAV, and such a complex process involves multiple cell signalling pathways that might be integrated to reflect the complexity of the disease. In this manner, our results are consistent with previous findings that postulated the crosstalk between TGF-β1 and VEGF together with the Notch signalling pathways, to act coordinately in space and time in the regulation of vascular morphogenesis [53]. With this consideration, our circulating miRNA signature, including miR-122, miR-130a, miR-486 and miR-718 expression, integrates the complexity associated with not only the bicuspid morphology of the aortic valve but also the progressive dilation of the ascending aorta. Furthermore, whether these miRNAs are directly involved in the aortic wall pathogenesis of BAV disease or if they are a consequence of the increased aortic shear stress generated by the anomalous aortic flow caused by BAV remains unknown.

Conclusions

We proposed miR-122, miR-130a, miR-486 and miR-718 as new molecular features associated with BAV and aortic dilation by the regulation of TGF-β1 and vascular remodelling mediated by the VEGF signalling pathway.

Abbreviations

TAV:

tricuspid aortic valve

TAVnondil:

tricuspid aortic valve with non-dilated aorta

BAV:

bicuspid aortic valve

BAVnondil:

bicuspid aortic valve with non-dilated aorta

miRNAs:

microRNAs

LV:

left ventricular

WSS:

wall systolic stress

MMP:

matrix metalloproteinases

TIMP:

inhibitors of matrix metalloproteinases

ECM:

extracellular matrix

TGF-β1 :

transforming growth factor-beta

VEGFA:

vascular endothelial growth factor A

References

  1. Alegret JM, Vernis JM, Palazón O, Borràs X, Duran I. Factors related to aortic regurgitation in the presence of a dilated aortic root. Am J Cardiol. 2005;95:417–20.

    Article  PubMed  Google Scholar 

  2. Michelena HI, Desjardins VA, Avierinos J-F, Russo A, Nkomo VT, Sundt TM, et al. Natural history of asymptomatic patients with normally functioning or minimally dysfunctional bicuspid aortic valve in the community. Circulation. 2008;117:2776–84.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Tzemos N, Therrien J, Yip J, Thanassoulis G, Tremblay S, Jamorski MT, et al. Outcomes in adults with bicuspid aortic valves. JAMA. 2008;300:1317–25.

    Article  CAS  PubMed  Google Scholar 

  4. Tzemos N, Lyseggen E, Silversides C, Jamorski M, Tong JH, Harvey P, et al. Endothelial function, carotid-femoral stiffness, and plasma matrix metalloproteinase-2 in men with bicuspid aortic valve and dilated aorta. J Am Coll Cardiol. 2010;55:660–8.

    Article  CAS  PubMed  Google Scholar 

  5. Michelena HI, Khanna AD, Mahoney D, Margaryan E, Topilsky Y, Suri RM, et al. Incidence of aortic complications in patients with bicuspid aortic valves. JAMA. 2011;306:1104–12.

    Article  CAS  PubMed  Google Scholar 

  6. Alegret JM, Ligero C, Vernis JM, Beltrán-Debón R, Aragonés G, Duran I, et al. Factors related to the need for surgery after the diagnosis of bicuspid aortic valve: one center’s experience under a conservative approach. Int J Med Sci. 2013;10:176–82.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Elefteriades JA. Thoracic aortic aneurysm: reading the enemy’s playbook. Yale J Biol Med. 2008;81:175–86.

    PubMed  PubMed Central  Google Scholar 

  8. Bissell MM, Hess AT, Biasiolli L, Glaze SJ, Loudon M, Pitcher A, et al. Aortic dilation in bicuspid aortic valve disease: flow pattern is a major contributor and differs with valve fusion type. Circ Cardiovasc Imaging. 2013;6:499–507.

    Article  PubMed  Google Scholar 

  9. Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. 2008;105:10513–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Goettsch C, Hutcheson JD, Aikawa E. MicroRNA in cardiovascular calcification: focus on targets and extracellular vesicle delivery mechanisms. Circ Res. 2013;112:1073–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ren J, Zhang J, Xu N, Han G, Geng Q, Song J, et al. Signature of circulating microRNAs as potential biomarkers in vulnerable coronary artery disease. PLoS ONE. 2013;8:e80738.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Naga Prasad SV, Duan Z-H, Gupta MK, Surampudi VSK, Volinia S, Calin GA, et al. Unique microRNA profile in end-stage heart failure indicates alterations in specific cardiovascular signaling networks. J Biol Chem. 2009;284:27487–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. D’Alessandra Y, Devanna P, Limana F, Straino S, Di Carlo A, Brambilla PG, et al. Circulating microRNAs are new and sensitive biomarkers of myocardial infarction. Eur Heart J. 2010;31:2765–73.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kin K, Miyagawa S, Fukushima S, Shirakawa Y, Torikai K, Shimamura K, et al. Tissue- and plasma-specific MicroRNA signatures for atherosclerotic abdominal aortic aneurysm. J Am Heart Assoc. 2012;1:e000745.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Fichtlscherer S, De Rosa S, Fox H, Schwietz T, Fischer A, Liebetrau C, et al. Circulating microRNAs in patients with coronary artery disease. Circ Res. 2010;107:677–84.

    Article  CAS  PubMed  Google Scholar 

  16. Alegret JM, Palazón O, Duran I, Vernis JM. Aortic valve morphology definition with transthoracic combined with transesophageal echocardiography in a population with high prevalence of bicuspid aortic valve. Int J Cardiovasc Imaging. 2005;21:213–7.

    Article  PubMed  Google Scholar 

  17. Gleeson TG, Mwangi I, Horgan SJ, Cradock A, Fitzpatrick P, Murray JG. Steady-state free-precession (SSFP) cine MRI in distinguishing normal and bicuspid aortic valves. J Magn Reson Imaging. 2008;28:873–8.

    Article  PubMed  Google Scholar 

  18. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  PubMed  Google Scholar 

  19. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.

    Google Scholar 

  21. Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, et al. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43:W460–6.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.

    Article  PubMed Central  Google Scholar 

  23. Chou C-H, Chang N-W, Shrestha S, Hsu S-D, Lin Y-L, Lee W-H, et al. miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database. Nucleic Acids Res. 2015;44:D239–47.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43:D146–52.

    Article  PubMed  Google Scholar 

  25. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  Google Scholar 

  26. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013;41:W77–83.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Nigam V, Sievers HH, Jensen BC, Sier HA, Simpson PC, Srivastava D, et al. Altered microRNAs in bicuspid aortic valve: a comparison between stenotic and insufficient valves. J Heart Valve Dis. 2010;19:459–65.

    PubMed  PubMed Central  Google Scholar 

  29. Yanagawa B, Lovren F, Pan Y, Garg V, Quan A, Tang G, et al. miRNA-141 is a novel regulator of BMP-2-mediated calcification in aortic stenosis. J Thorac Cardiovasc Surg. 2012;144:256–62.

    Article  CAS  PubMed  Google Scholar 

  30. Wu J, Song H-F, Li S-H, Guo J, Tsang K, Tumiati L, et al. Progressive aortic dilation is regulated by miR-17-associated miRNAs. J Am Coll Cardiol. 2016;67:2965–77.

    Article  CAS  PubMed  Google Scholar 

  31. Oury C, Servais L, Bouznad N, Hego A, Nchimi A, Lancellotti P. MicroRNAs in valvular heart diseases: potential role as markers and actors of valvular and cardiac remodeling. Int J Mol Sci. 2016. doi:10.3390/ijms17071120.

    PubMed  PubMed Central  Google Scholar 

  32. Siu SC, Silversides CK. Bicuspid aortic valve disease. J Am Coll Cardiol. 2010;55:2789–800.

    Article  PubMed  Google Scholar 

  33. Abdulkareem N, Smelt J, Jahangiri M. Bicuspid aortic valve aortopathy: genetics, pathophysiology and medical therapy. Interact Cardiovasc Thorac Surg. 2013;17:554–9.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Mathieu P, Bossé Y, Huggins GS, Della Corte A, Pibarot P, Michelena HI, et al. The pathology and pathobiology of bicuspid aortic valve: state of the art and novel research perspectives. J Pathol Clin Res. 2015;1:195–206.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Rocchiccioli S, Cecchettini A, Panesi P, Farneti PA, Mariani M, Ucciferri N, et al. Hypothesis-free secretome analysis of thoracic aortic aneurysm reinforces the central role of TGF-β cascade in patients with bicuspid aortic valve. J Cardiol. 2016. doi:10.1016/j.jjcc.2016.05.007.

    PubMed  Google Scholar 

  36. Heldin CH, Miyazono K, ten Dijke P. TGF-beta signalling from cell membrane to nucleus through SMAD proteins. Nature. 1997;390:465–71.

    Article  CAS  PubMed  Google Scholar 

  37. Forte A, Galderisi U, Cipollaro M, De Feo M, Corte AD. Epigenetic regulation of TGF-1 signalling in dilative aortopathy of the thoracic ascending aorta. Clin Sci. 2016;130:1389–405.

    Article  PubMed  Google Scholar 

  38. Li W, Li Q, Jiao Y, Qin L, Ali R, Zhou J, et al. Tgfbr2 disruption in postnatal smooth muscle impairs aortic wall homeostasis. J Clin Invest. 2014;124:755–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Forte A, Della Corte A, Grossi M, Bancone C, Provenzano R, Finicelli M, et al. Early cell changes and TGFβ pathway alterations in the aortopathy associated with bicuspid aortic valve stenosis. Clin Sci (Lond). 2013;124:97–108.

    Article  CAS  Google Scholar 

  40. Jopling C. Liver-specific microRNA-122: biogenesis and function. RNA Biol. 2012;9:137–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, et al. miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006;3:87–98.

    Article  CAS  PubMed  Google Scholar 

  42. Beaumont J, López B, Hermida N, Schroen B, San José G, Heymans S, et al. microRNA-122 down-regulation may play a role in severe myocardial fibrosis in human aortic stenosis through TGF-β1 up-regulation. Clin Sci (Lond). 2014;126:497–506.

    Article  CAS  Google Scholar 

  43. Fedak PWM, Verma S, David TE, Leask RL, Weisel RD, Butany J. Clinical and pathophysiological implications of a bicuspid aortic valve. Circulation. 2002;106:900.

    Article  PubMed  Google Scholar 

  44. Chen Y, Gorski DH. Regulation of angiogenesis through a microRNA (miR-130a) that down-regulates antiangiogenic homeobox genes GAX and HOXA5. Blood. 2008;111:1217–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wu W-H, Hu C-P, Chen X-P, Zhang W-F, Li X-W, Xiong X-M, et al. MicroRNA-130a mediates proliferation of vascular smooth muscle cells in hypertension. Am J Hypertens. 2011;24:1087–93.

    Article  CAS  PubMed  Google Scholar 

  46. Holliday CJ, Ankeny RF, Jo H, Nerem RM. Discovery of shear- and side-specific mRNAs and miRNAs in human aortic valvular endothelial cells. Am J Physiol Heart Circ Physiol. 2011;301:H856–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Navon R, Wang H, Steinfeld I, Tsalenko A, Ben-Dor A, Yakhini Z. Novel rank-based statistical methods reveal microRNAs with differential expression in multiple cancer types. PLoS ONE. 2009;4:e8003.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Xue M, Yao S, Hu M, Li W, Hao T, Zhou F, et al. HIV-1 Nef and KSHV oncogene K1 synergistically promote angiogenesis by inducing cellular miR-718 to regulate the PTEN/AKT/mTOR signaling pathway. Nucleic Acids Res. 2014;42:9862–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Leng R, Zha L, Tang L. MiR-718 represses VEGF and inhibits ovarian cancer cell progression. FEBS Lett. 2014;588:2078–86.

    Article  CAS  PubMed  Google Scholar 

  50. Sugimachi K, Matsumura T, Hirata H, Uchi R, Ueda M, Ueo H, et al. Identification of a bona fide microRNA biomarker in serum exosomes that predicts hepatocellular carcinoma recurrence after liver transplantation. Br J Cancer. 2015;112:532–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Shi J, Liu H, Wang H, Kong X. MicroRNA expression signature in degenerative aortic stenosis. Biomed Res Int. 2016;2016:4682172.

    PubMed  PubMed Central  Google Scholar 

  52. Cheuk BLY, Cheng SWK. Identification and characterization of microRNAs in vascular smooth muscle cells from patients with abdominal aortic aneurysms. J Vasc Surg. 2014;59:202–9.

    Article  PubMed  Google Scholar 

  53. Holderfield MT, Hughes CCW. Crosstalk between vascular endothelial growth factor, notch, and transforming growth factor-β in vascular morphogenesis. Circ Res. 2008;102:637.

    Article  CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

NMM performed the experiments, analyzed the data, and wrote the manuscript; RBD performed the experiments; IB performed the experiments and analyzed the data; MF contributed in sample collection and writing the manuscript; JMA conceived and designed the study, analyzed the data and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank all the participants for their contribution and participation.

Funding was provided by Sociedad Española de Cardiología (Grant Beca de Investigación Clínica 2016).

Competing interests

The authors declare that they have no competing interests.

Availability of data and supporting materials

All relevant data are within the paper and its supporting information files.

Ethics approval and consent to participate

Study approved by the local ethical commitee (Hospital Universitari de Sant Joan, Reus); each patient signed a written informed consent.

Funding

This work was supported by a grant for clinical investigation from the Sociedad Española de Cardiología.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Josep M. Alegret.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Martínez-Micaelo, N., Beltrán-Debón, R., Baiges, I. et al. Specific circulating microRNA signature of bicuspid aortic valve disease. J Transl Med 15, 76 (2017). https://doi.org/10.1186/s12967-017-1176-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-017-1176-x

Keywords