Specific circulating microRNA signature of bicuspid aortic valve disease
Journal of Translational Medicine volume 15, Article number: 76 (2017)
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.
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.
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  and the risk of aortic dissection and thus often requires prophylactic aortic surgery [2–6]. 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 . 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 , vulnerable coronary artery disease , heart failure , myocardial infarction , abdominal aortic aneurysm  and coronary artery disease . 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.
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 , on a transesophageal echocardiogram  or on a cardiac magnetic resonance . 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.
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.
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 , 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 . The Benjamini and Hochberg method was used to adjust the p values for multiple testing and to control the false discovery rate .
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) , 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) , miRTarBase (2016 version; mirtarbase.mbc.nctu.edu.tw) , DIANA-microT-CDS (v3; microrna.gr/microT-CDS)  and miRDB (2016 version; mirdb.org)  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)  and WebGestalt (webgestalt.org) , respectively.
Cytoscape version 3.0.1  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.
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).
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.
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).
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.
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.
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).
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.  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.  determined that the down-regulation of miR-141 in stenotic bicuspid leaflets regulates BMP-2-mediated valve calcification. Wu et al.  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 . 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 [32–34]. 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 . 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 . 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 . 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 . Moreover, Forte et al. , 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 , 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  and is considered a central player in liver biology and in the regulation of cholesterol and fatty acid metabolism . Interestingly, Beaumont et al.  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 . 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 . Furthermore, the expression of miR-130a in the aorta is correlated with vascular remodelling in spontaneously hypertensive rats (SHRs) . Regarding miR-486, Holliday et al.  demonstrated that the expression of miR-486 in human aortic valvular endothelial cells was shear dependent, and Navon et al.  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 [48–51] 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 . 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 . 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.
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.
tricuspid aortic valve
tricuspid aortic valve with non-dilated aorta
bicuspid aortic valve
bicuspid aortic valve with non-dilated aorta
wall systolic stress
inhibitors of matrix metalloproteinases
- TGF-β1 :
transforming growth factor-beta
vascular endothelial growth factor A
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.
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.
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.
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.
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.
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.
Elefteriades JA. Thoracic aortic aneurysm: reading the enemy’s playbook. Yale J Biol Med. 2008;81:175–86.
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.
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.
Goettsch C, Hutcheson JD, Aikawa E. MicroRNA in cardiovascular calcification: focus on targets and extracellular vesicle delivery mechanisms. Circ Res. 2013;112:1073–84.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Agarwal V, Bell GW, Nam J-W, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
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.
Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43:D146–52.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013;41:W77–83.
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.
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.
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.
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.
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.
Siu SC, Silversides CK. Bicuspid aortic valve disease. J Am Coll Cardiol. 2010;55:2789–800.
Abdulkareem N, Smelt J, Jahangiri M. Bicuspid aortic valve aortopathy: genetics, pathophysiology and medical therapy. Interact Cardiovasc Thorac Surg. 2013;17:554–9.
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.
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.
Heldin CH, Miyazono K, ten Dijke P. TGF-beta signalling from cell membrane to nucleus through SMAD proteins. Nature. 1997;390:465–71.
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.
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.
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.
Jopling C. Liver-specific microRNA-122: biogenesis and function. RNA Biol. 2012;9:137–42.
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.
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.
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.
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.
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.
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.
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.
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.
Leng R, Zha L, Tang L. MiR-718 represses VEGF and inhibits ovarian cancer cell progression. FEBS Lett. 2014;588:2078–86.
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.
Shi J, Liu H, Wang H, Kong X. MicroRNA expression signature in degenerative aortic stenosis. Biomed Res Int. 2016;2016:4682172.
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.
Holderfield MT, Hughes CCW. Crosstalk between vascular endothelial growth factor, notch, and transforming growth factor-β in vascular morphogenesis. Circ Res. 2008;102:637.
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.
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).
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.
This work was supported by a grant for clinical investigation from the Sociedad Española de Cardiología.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
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