Thoracic Aortic Aneurysm (TAA) is characterized by the dilation of the aorta and is fatal if not diagnosed and treated appropriately. The underlying genetic mechanisms have not been completely delineated, so better knowledge of the physiopathology of TAAs is needed to improve detection and therapy. MicroRNAs (miRNAs) regulate gene expression post-transcriptionally and are known to be involved in cardiovascular diseases (CVDs). The current study aimed to identify miRNAs that can be used as possible biomarkers for the early diagnosis of patients with ascending TAAs (ATAAs). MiRNA expression was profiled by NanoString nCounter technology using 12 samples including tissue and pre- and post-surgical plasma from ATAA patients. Four miRNAs were selected and further validated by real time polymerase chain reaction (RT-PCR) in 22 plasma samples from which three miRNAs (hsa-miR140-5p, hsa-miR-191-5p and hsa-miR-214-3p) showed significant expression level differences between the two types of plasma samples. Further analyses of the corresponding predicted target genes by these miRNAs, revealed two genes (Myotubularin-related protein 4 (MTMR4) and Phosphatase 1 catalytic subunit β (PPP1CB)) whose expression was inversely correlated with the expression of their respective miRNAs. Overall, in this pilot study, we identified three miRNAs that might serve as potential biomarkers and therapeutic targets in ATAA.

Thoracic aortic aneurysm (TAA) may be defined as a localized or diffuse dilatation of the aorta to at least a 50 percent increase in diameter compared with the expected normal diameter. The aortic root, ascending aorta, aortic arch, or descending aorta may be affected and are classified accordingly [1]. Sixty percent of TAAs involve the aortic root and/or ascending aorta, 40% involve the descending aorta, 10% involve the arch, and 10% involve the thoracoabdominal aorta [2]. The etiology, natural history, and treatment of thoracic aneurysms differ for each of these segments. Ascending TAAs (ATAAs) most often result from cystic medial degeneration, which appears histologically as smooth muscle cell dropout and elastic fiber degeneration. Medial degeneration leads to weakening of the aortic wall, which in turn results in aortic dilatation and aneurysm formation [2]. Cystic medial degeneration occurs normally to some extent with aging, but the process is accelerated by hypertension [3]. ATAAs can either be sporadic or genetically inherited. In younger patients, it is more often associated with a genetic predisposition that can be familial [4] or related to defined genetic disorders such as Marfan syndrome.

The natural history of TAA is one of the slow expansion with a progressive increase in the risk of aortic dissection (AD) at larger aortic sizes. The slow expansion of TAA means that most patients with TAA are asymptomatic and many patients will succumb to other disease processes without being aware that TAA is present. Aneurysms that do produce symptoms are typically very large and at risk for rupture, which is associated with high rates of morbidity and mortality. Death from aneurysmal rupture is one of the fifteen leading causes of deaths in most series. According to the Centers for Disease Control and Prevention, the incidence of ATAA is estimated to be approximately 10 per 100000 person-years [5].

Despite their diverse etiologies, all TAAs are prone to two lethal acute events: AD and aortic rupture. In AD, blood is diverted from its usual location within the lumen of the aorta into a false lumen within the media through a tear in the intima. AD is a catastrophic event in which the layers of the aortic wall split apart. An aortic rupture includes a complete tear or transection through all three layers of the aorta, resulting in massive internal hemorrhage [6]. Eleftheriades et al., 2005 [7], have described the detection of TAAs as essential due to the lethality of the disease, as well as difficult as approximately 95% of TAAs patients are asymptomatic before presented with a devastating fatal or near-fatal dissection or rupture. The mortality is said to be 1% per hour once AD has occurred. Symptomatic patients or those that have complications related to the aneurysm should undergo surgical repair. The mortality of elective surgical repair of ascending aortic aneurysms in large centers is 3–5% [2]. Conservative management of asymptomatic TAA aims to lessen stress on the aorta and limit further aortic expansion. Asymptomatic patients who do not meet criteria for repair also require ongoing aneurysm observation [1,8,9]. Screening for asymptomatic TAAs before an event occurs is of great clinical importance. However, there are no standard screening programs for TAAs. Consequently, TAAs are commonly detected incidentally by contrast-enhanced CT scanning and MR angiography. It is critically important to detect the aneurysm before it reaches dangerous diameters, as surgical treatment is safest and most effective on an elective basis. These detection modalities are not safe or cost-effective as screening methods for the general population. Therefore, better screening methods are warranted for TAAs which are non-invasive, cheap, safe, and widely available. A biomarker is a molecular indicator of a biological process or state and in clinical application, biomarkers are reflective of the pathogenic activity of a disease. A biomarker must be able to measure objectively and quantitatively. The quantity of the biomarker can correlate with the type, risk, or severity of disease. Therefore, ideal biomarkers should be able to predict aneurysm formation well before dilatation, can be detected via imaging techniques, and may also provide information on how close an aneurysm is to dissection or rupture. While a number of biomarkers of aortic aneurysms have been identified and shown promise, they have not led to any clinically utilized testing programs.

MicroRNAs (miRNAs) are a group of highly conserved, small noncoding RNAs 21–25 nucleotides long, which usually negatively modulate gene expression at the post-transcriptional level by incomplete or complete complementary binding to target sequences within the 3′ untranslated regions (UTRs) of mRNA [10]. Therefore, an increased expression of a specific miRNA leads to decreased expression of its target genes and conversely [11].

MiRNAs have been discovered to be important regulators for several diseases, therefore, they can also be significant biomarkers in the detection of diseases, since they have several characteristics of an ideal biomarker [12]. MiRNAs are present in several body fluids which can be collected without or via minimally invasive methods, such as blood, saliva, and urine. They are stable in blood plasma and can withstand extreme conditions such as high pH and they can be stored for a long time under appropriate conditions without their biological integrity being affected [13]. They are differentially expressed in different tissues making them specific for each tissue, which is another significant characteristic for a biomarker [14]. In several studies, it was identified that the expression of miRNA in diseased samples compared with healthy controls is different, therefore these small molecules can be used for the detection of a condition or the progression of the disease [15,16]. The aberrant expression of miRNAs has also been linked to cardiovascular diseases (CVDs) and is likely that they also play a key role in aneurysm formation [17–20]. Studies have already investigated and identified miRNAs which are differentially expressed in aneurysmal sample [18,20–25]. These results show promise for miRNA-based biomarkers that can be used to detect and screen for TAAs. Additionally, by studying both the relevant genes and miRNAs in conjunction with the intermodulation between them may clarify the molecular mechanism of the disease and provide reliable molecular targets for detection also leading to prevention and treatment.

The current study was divided into two phases. During the primary phase, the aim was to identify miRNAs with significantly different expression levels between plasma samples from patients, pre- and post-operatively. The secondary phase, aimed to determine whether the predicted target genes of the miRNAs detected, showed a specific pattern of expression in the two types of plasma samples and speculate that their expression changes could be associated with TAAs.

Sample collection

Samples were collected from 11 male and 4 female patients with sporadic ATAAs undergoing surgery for the repair of the aneurysm at The American Medical Centre, Cyprus. Briefly, the ascending sort extending from the sino-tubular junction to the proximal aortic arch was resected. The resected aorta was replaced with a synthetic graft. Any bicuspid or tricuspid valve with moderate to severe aortic valve regurgitation or stenosis were replaced. In cases of AD with aortic valve insufficiency, we proceeded with resuspension of the valve commissures and aortic valve repair.

The tissue samples (T) were taken from the widest region of the ascending aorta at the time of surgery as well as blood samples before (BB) and 7–8 days after (BA) the surgery, for each participant. The criteria used for the selection of these samples was that all patients presented with an aortic diameter ≥ 50 mm and did not have any transfusions during or immediately after the surgery. Patients’ characteristics including, aortic valve type, presence of hypertension, diabetes or hyperlipidemia were recorded (Table 1). Upon initial evaluation, the patients did not have any manifestation of any phenotypic findings (i.e. iris flocculi, livedo reticularis, aortic valve abnormalities) hinting toward connective tissue disorders or aortic pathology. The excised tissue samples were rinsed in saline solution cut and snap-frozen in liquid nitrogen and stored at −80°C. The blood samples were collected in EDTA tubes and steps were followed in order to reduce the issue of hemolysis. When the blood samples were collected, the needle was removed from the syringe and the blood slowly expelled into the EDTA tube. Once the samples were collected, within 20–30 min they were centrifuged and the plasma collected. This involved a double centrifugation at 35 rpm for 10 min, where the plasma was initially removed carefully to prevent carryover of any cell debris and placed into a clean 2-ml tube and then re-centrifuged for further 10 min and carefully transferred to another clean 2-ml tube. The plasma samples were stored in −80°C until needed. The presentstudy was approved by the Cyprus National Bioethics Committee and each patient gave their signed consent.

Table 1
Patients characteristics
CharacteristicsPrimary phaseSecondary phase
Number (%)
Gender (Male) 4 (100%) 7/11 (63.6%) 
Average age 73.2 ± 9.4 62.5 ± 7 
Tricuspid aortic valve 4 (100%) 6/11 (54.5%) 
Hypertension 0 (0%) 4/11 (36.4%) 
Diabetes 0 (0%) 1/11 (9.1%) 
Hyperlipidemia 0 (0%) 4/11 (36.4%) 
CharacteristicsPrimary phaseSecondary phase
Number (%)
Gender (Male) 4 (100%) 7/11 (63.6%) 
Average age 73.2 ± 9.4 62.5 ± 7 
Tricuspid aortic valve 4 (100%) 6/11 (54.5%) 
Hypertension 0 (0%) 4/11 (36.4%) 
Diabetes 0 (0%) 1/11 (9.1%) 
Hyperlipidemia 0 (0%) 4/11 (36.4%) 

Primary phase: investigation of differentially expressed miRNAs between plasma samples

For the primary phase, miRNA expression analysis, a total of 12 samples were used. These were from four male patients with ATAAs undergoing surgery, including the aortic tissue sample and the plasma samples (BB and BA) from each individual.

Total RNA extraction

For the plasma samples total RNA was extracted using Plasma/Serum RNA Purification Midi Kit (Norgen Biotek Corp.) and concentrated and purified further using RNA Clean-up and Micro-Elute Kit (Norgen Biotek Corp.), following the manufacturer’s instructions. From the four tissue samples, a small piece of approximately 5 mg was used where it was initially homogenized and then sonicated to further disrupt the cells of the tissue. Total RNA was isolated, from each sample, using miRCURY RNA Isolation Kit-Tissue (Exiqon), in accordance with the manufacturer’s instructions. The A260/280 ratio and the A260/230 ratios for the purified RNA was measured on a Nanodrop 8000 spectrophotometer (Thermo Fisher Scientific). These data are summarized in Supplementary Table S1. To monitor for hemolysis, an RNA sample extracted from a hemolyzed blood sample was used as a control. The level of the erythrocyte specific hsa-miR-451 was compared between this control sample with each of the extracted plasma samples prior to any experimental work undertaken. This hemolysis indicator allowed for the elimination of samples with high expression levels of hsa-miR-451 [26].

miRNA expression assay

The multiplexed NanoString nCounter Human v3 miRNA expression assay (NanoString Technologies, Seattle, WA, U.S.A.) was used to profile 798 human miRNAs. Samples from four patients were initially used, which included T, BB and (BA) samples after the repair operation. The assay was performed according to the manufacturer’s protocol. Briefly, up to 100 ng of total RNA was used as input material, with 3 μl of the threefold-diluted sample (in the tissue samples only). A specific DNA tag was ligated on to the 3′ end of each mature miRNA, providing an exclusive identification for each miRNA species in the sample. The tagging was performed in a multiplexed ligation reaction utilizing reverse complementary bridge oligonucleotides to dispose the ligation of each miRNA to its designated tag. All hybridization reactions were incubated at 64°C for 18 h. Excess tags were then removed, and the resulting material was hybridized to a panel of fluorescently labeled, bar-coded reporter probes specific to the miRNA of interest. Abundances of miRNAs were quantified with the nCounter Prep Station via counting individual fluorescent barcodes and quantifying target miRNA molecules present in each sample. The count number for each sample was detected and the mean value for each sample type (BB, BA and T) was calculated. Comparisons were made between tissue samples vs plasma from before the operation (T vs BB), tissue samples vs plasma samples from after the operation (T vs BA) and plasma samples before the operation vs plasma sample after the operation (BB vs BA).

Statistical analysis

All statistical analysis was conducted using R software version 3.2. NanoString data pre-processing and normalization followed by differential expression analysis was performed using DEseq2 package [27] and in-house scripts [28] identifying the mean values from the patients’ samples. The mean value of negative controls was set as the lower threshold for each sample; only miRNAs with at least 50% of their values above the lower threshold were included in the analysis. The false discovery rate method was used to control multiple observations, and the adjusted P-value was used for further analysis. Fold-change difference were calculated for each comparison and the P-value was determined where P-value <0.05 was considered significant [28–30]. The miRNAs which showed significant expression level differences between BB and BA samples were selected for further investigation. MirPath software from Diana tool [31,32], a prediction algorithm designed using databases of experimentally verified miRNA targets, was used in order to identify pathways which were predicted to be significantly (P-value <0.05) associated with the eight miRNAs.

Validation study

The miRNAs that were predicted to be significantly associated with pathways correlated with ATAAs were used for validation purposes on an additional 22 plasma samples (11 BB and 11 BA) by real time polymerase chain reaction (RT-PCR; Table 1).

cDNA synthesis and RT-PCR

For cDNA synthesis, qScript microRNA cDNA Synthesis Kit (Quanta Biosciences) was used following manufacturer’s instructions. For RT-PCR, mi-real-time Eva Green Master mix (ROX+) (Metabion) was used. For this reaction, 5 μl of the master mix were used, with 1 μl of PCR primer and 4 μl of diluted cDNA (3:20). The primer sequences for the miRNAs were available online by Quantabio (http://www.quantabio.com/products/microrna-profiling) and ordered from Metabion. Table 2 shows the sequences selected for each miRNA.

Table 2
Primer sequences of miRNAs selected for validation
miRNAPrimer sequence
hsa-miR-140-5p CGCAGTGGTTTTACCCTATG 
hsa-miR-191-5p AACGGAATCCCAAAAGCAG 
hsa-miR-214-3p GCACAGCAGGCACAGACAG 
hsa-let-7i-5p CGTTCTGAGGTAGTAGTTTGTGCT 
hsa-miR-103a-3p (control) CAGCATTGTACAGGGCTATGAA 
miRNAPrimer sequence
hsa-miR-140-5p CGCAGTGGTTTTACCCTATG 
hsa-miR-191-5p AACGGAATCCCAAAAGCAG 
hsa-miR-214-3p GCACAGCAGGCACAGACAG 
hsa-let-7i-5p CGTTCTGAGGTAGTAGTTTGTGCT 
hsa-miR-103a-3p (control) CAGCATTGTACAGGGCTATGAA 

The amplification conditions included: denaturation at 95°C (10 min), amplification at 95°C (15 sec) and extension at 60°C (1 min) for 40 cycles using the 7500 Real-Time PCR System (Applied Biosystems). The miRNA, hsa-miR103a-3p was used as an endogenous control for all the reactions. In addition, non-enzyme controls (NEC) and non-template controls (NTC) were used where, in NEC, reverse transcriptase was not added during synthesis of cDNA and in NTC, no cDNA template was added, instead it was replaced by water. The RT-PCR were performed in triplicates for each sample.

Data analysis

For the analysis of the RT-PCR data, normalized Ct values (ΔCt) were calculated for each sample according to the mean threshold of the miRNA of interest and the endogenous control. The ΔCt of each miRNA in BB samples were compared with the ΔCt of BA samples and Wilcoxon signed-rank test was performed where a P- value was considered significant at <0.05.

Secondary phase: targeted genes for the validated miRNAs

Targeted genes

A total of 21 prediction tools, including BcmicrO (https://tools4mirs.org/software/target_prediction/bcmicro/) [33], BiTargeting (https://www.cs.bgu.ac.il/∼vaksler/BiTargeting.htm) [34], CoMeTa (http://cometa.tigem.it/) [35], Cupid (http://cupidtool.sourceforge.net/) [36], DIANA (http://diana.imis.athena-innovation.gr/DianaTools/index.php, ElMMo3 (http://www.mirz.unibas.ch/ElMMo3/index.php) [37], MBStar (https://www.isical.ac.in/∼bioinfo_miu/MBStar/MBStar_download20.htm) [38] microrna.org (http://www.microrna.org/microrna/home.do) [39], MirAncesTar (https://www.cs.mcgill.ca/∼blanchem/mirancestar/) [40,41], miRcode (http://www.mircode.org/) [42], mirCoX (http://210.212.254.116/mircox/pages/index.php) [43], miRDB (http://www.mirdb.org/) [44], MirMAP (https://mirmap.ezlab.org/) [45], MirTar (http://mirtar.mbc.nctu.edu.tw/human/) [46], miRTar2GO (http://www.mirtar2go.org/) [47], Mirza-G (http://www.clipz.unibas.ch/mirzag/) [48], MultiMiTar (https://www.isical.ac.in/∼bioinfo_miu/multimitar.htm) [49], PicTar (https://pictar.mdc-berlin.de/) [50], RepTar [51], RNAhybrid (https://bibiserv.cebitec.uni-bielefeld.de/rnahybrid) [52] and TargetRank (http://hollywood.mit.edu/targetrank/) [53], were used in order to identify predicted target genes for the miRNAs which showed significant differences between BB and BA samples. The genes selected were those detected to be predicted targets of two or more of the validated miRNAs. In order to identify the most promising target genes, a comprehensive literature search was performed in PubMed using specific key terms; ‘gene name’ and/or ‘aneurysm’, ‘cardiovascular disease’, ‘aorta’ and ‘miRNA of interest’. In addition, NCBI-Gene (https://www.ncbi.nlm.nih.gov/gene) and UniProt databases (https://www.uniprot.org/) were searched in order to identify any possible association of these proteins with aneurysms [54].

cDNA synthesis and RT-PCR

Using the same 22 samples as with the miRNA validation phase, the expression levels of the mRNAs from the predicted nine target genes chosen were examined. Primers were designed based on the UTR of the mRNAs using TargetScan Human and Ensembl genome browser. cDNA synthesis and RT-PCR were performed using qScript microRNA cDNA Synthesis Kit (Quanta Biosciences) and EvaGreen master mix (ROX+), as described earlier. As an internal control, 18S rRNA was used. Table 3 shows the mRNA primer sequence for the target genes.

Table 3
mRNA primer sequence
Sequence nameSequence
CCND2-Forward AGTTATTGCTGGTGCAAAGA 
CCND2-Reverse GCCTCTAAGAAGTGGAGAGG 
ABLIM1-Forward AGCTCTGTCGAATTCATGG 
ABLIM1-Reverse CTTTGAGTCACACACACTCG 
CRKL- Forward TGTCTTAGTGTTTTAGTGACTAGGG 
CRKL-Reverse AGATGATGGATGAACCACAC 
MTMR4-Forward TTAGGGGAACAGTTCTAGGG 
MTMR4-Reverse TAGGAAAAGTCAGGTGCAAA 
HEY1-Forward TTGGTCTGTTTTTCTCCTTG 
HEY1-Reverse CAACGTTTTGCTCAGAAATA 
FNDC5-Forward ACAGGATAAAGGGCAGATGT 
FNDC5-Reverse TGCTCTAAGTGGATCAGAGG 
NFIA-Forward CATGCGACTTCAAGAAGGTTT 
NFIA-Reverse TGTGTCTTAATGCACTCTACTTTATCC 
PPP1CB-Forward TCATTTTCAGCACAGTGCAA 
PPP1CB-Reverse CAGTCAGGCTTCATATCAATTAGAC 
TJP1-Forward TCGTCACTGCTTAACTTCACA 
TJP1-Reverse TAGGCCAGGGCCATAGTAAA 
Sequence nameSequence
CCND2-Forward AGTTATTGCTGGTGCAAAGA 
CCND2-Reverse GCCTCTAAGAAGTGGAGAGG 
ABLIM1-Forward AGCTCTGTCGAATTCATGG 
ABLIM1-Reverse CTTTGAGTCACACACACTCG 
CRKL- Forward TGTCTTAGTGTTTTAGTGACTAGGG 
CRKL-Reverse AGATGATGGATGAACCACAC 
MTMR4-Forward TTAGGGGAACAGTTCTAGGG 
MTMR4-Reverse TAGGAAAAGTCAGGTGCAAA 
HEY1-Forward TTGGTCTGTTTTTCTCCTTG 
HEY1-Reverse CAACGTTTTGCTCAGAAATA 
FNDC5-Forward ACAGGATAAAGGGCAGATGT 
FNDC5-Reverse TGCTCTAAGTGGATCAGAGG 
NFIA-Forward CATGCGACTTCAAGAAGGTTT 
NFIA-Reverse TGTGTCTTAATGCACTCTACTTTATCC 
PPP1CB-Forward TCATTTTCAGCACAGTGCAA 
PPP1CB-Reverse CAGTCAGGCTTCATATCAATTAGAC 
TJP1-Forward TCGTCACTGCTTAACTTCACA 
TJP1-Reverse TAGGCCAGGGCCATAGTAAA 

Data analysis

The expression levels of the mRNAs were analyzed according to their average ΔCt as described earlier and they were further analyzed according to their fold-change expression differences between BB and BA samples. In cases were the fold-change was 1.5 or above, the difference was considered significant. Validated target genes of deregulated miRNAs were determined using MirTarbase (http://mirtarbase.mbc.nctu.edu.tw/php/index.php), which is a database of manually curated miRNA target gene interactions, obtained from the literature [55]. Data underlying the present study are available at the NCBI Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE151378.

Primary phase: investigation of differentially expressed miRNAs between plasma samples

Using the gene expression profile obtained from the Nanostring data, a total of 80 miRNAs were detected in eight plasma samples and they were also confirmed in four tissue samples. The detection of the circulating levels of the miRNAs may be dependent on the diseased aortic tissue samples either by originating or being enriched in this particular tissue. Eight of these miRNAs, showed significant difference (P≤0.001) in their expression levels between BB and BA plasma samples and in all cases the expression was lower in BB samples while the levels were much higher in the tissue samples (Table 4).

Table 4
miRNAs showing significant difference between BB and BA samples
miRNAMean BBMean BAP-valueMean T
Hsa-let-7c-5p 0.22 5.47 0.007 1444.88 
Hsa-let-7i-5p 3.30 0.002 748.18 
Hsa-miR-140-5p 0.32 3.92 0.001 258.81 
Hsa-miR-191-5p 0.18 2.89 0.021 498.97 
Hsa-miR-214-3p 2.26 0.007 48.27 
Hsa-miR-223-3p 0.64 15.71 0.001 96.89 
Hsa-miR-23a-3p 0.98 8.83 0.015 12306.94 
Hsa-miR-30b-5p 2.51 0.001 8.69 
miRNAMean BBMean BAP-valueMean T
Hsa-let-7c-5p 0.22 5.47 0.007 1444.88 
Hsa-let-7i-5p 3.30 0.002 748.18 
Hsa-miR-140-5p 0.32 3.92 0.001 258.81 
Hsa-miR-191-5p 0.18 2.89 0.021 498.97 
Hsa-miR-214-3p 2.26 0.007 48.27 
Hsa-miR-223-3p 0.64 15.71 0.001 96.89 
Hsa-miR-23a-3p 0.98 8.83 0.015 12306.94 
Hsa-miR-30b-5p 2.51 0.001 8.69 

Plasma samples before (BB) and after (BA) repair surgery as well as the removed aortic sample (T).

he hsa-miR-223-3p presented the highest expression difference between BB and BA samples while the hsa-miR-191-5p, showed the lowest significance between the two sample types. For three, miRNAs (hsa-let-7i-5p, hsa-miR-214-3p and hsa-miR-30b-5p), the mean BB values were zero and this could be due to the extremely low expression, therefore they did not reach the threshold and could not be detected.

miRNAs and ATAA pathways

In order to investigate any possible correlation between the eight miRNAs and ATAA pathways, the MirPath software was used. Three miRNAs (hsa-let-7i-5p, hsa-miR-140-5p, and hsa-miR-214-3p) were identified that significantly interact with relevant pathways, namely the Extracellular Matrix (ECM)–Receptor Interaction Pathway and one (hsa-miR-191-5p) with the Focal Adhesion Pathway (Table 5).

Table 5
miRNAs and predicted ATAA pathways
miRNAsPathwayP-value
Hsa-let-7i-5p ECM–Receptor Interaction 5.16 e-7 
Hsa-miR-140-5p ECM–Receptor Interaction 0.00017 
Hsa-miR-191-5p Focal Adhesion 0.049 
Hsa-miR-214-3p ECM–Receptor Interaction 1.26 e-14 
miRNAsPathwayP-value
Hsa-let-7i-5p ECM–Receptor Interaction 5.16 e-7 
Hsa-miR-140-5p ECM–Receptor Interaction 0.00017 
Hsa-miR-191-5p Focal Adhesion 0.049 
Hsa-miR-214-3p ECM–Receptor Interaction 1.26 e-14 

miRNA validation

The above four miRNAs chosen for validation in 22 plasma samples, showed lower expression in BB compared with BA samples, further supporting the results obtained using the nCounter miRNA Expression Assay. However, only hsa-miR-140-5p, hsa-miR-191-5p and hsa-miR-214-3p showed statistically different expression levels between the two plasma sample types with P-values, 0.0004, 0.007, and 2.83e−6, respectively (Figure 1). A flow chart summarizing the primary phase of the study is shown in Figure 2A.

Differentially expressed plasma miRNAs from ATAA patients

Figure 1
Differentially expressed plasma miRNAs from ATAA patients

Relative levels of hsa-let-7i-5p, hsa-miR-140-5p, hsa-miR-191-5p, and hsa-miR-214-3p were validated in plasma samples from 11 patients collect before (BB) and after (BA) repair surgery. The following three, hsa-miR-140-5p (P-value = 0.0004), hsa-miR-191-5p (P-value = 0.007), and hsa-miR-214-3p (P-value = 2.83e−6) showed statistically significant differences. The expression threshold (Ct) values of each miRNA were normalized to hsa-miR-103a-3p. *P<0.05; NS, non-significant.

Figure 1
Differentially expressed plasma miRNAs from ATAA patients

Relative levels of hsa-let-7i-5p, hsa-miR-140-5p, hsa-miR-191-5p, and hsa-miR-214-3p were validated in plasma samples from 11 patients collect before (BB) and after (BA) repair surgery. The following three, hsa-miR-140-5p (P-value = 0.0004), hsa-miR-191-5p (P-value = 0.007), and hsa-miR-214-3p (P-value = 2.83e−6) showed statistically significant differences. The expression threshold (Ct) values of each miRNA were normalized to hsa-miR-103a-3p. *P<0.05; NS, non-significant.

Close modal

Summary of experimental design and analysis

Figure 2
Summary of experimental design and analysis

(A) Primary phase of experimental work which included the investigation of differentially expressed miRNAs between plasma samples from patients collected before (BB) and after (BA) repair surgery. (B) Secondary phase involving the investigation of the target genes for the validated miRNAs

Figure 2
Summary of experimental design and analysis

(A) Primary phase of experimental work which included the investigation of differentially expressed miRNAs between plasma samples from patients collected before (BB) and after (BA) repair surgery. (B) Secondary phase involving the investigation of the target genes for the validated miRNAs

Close modal

Secondary phase: targeted genes for the validated miRNAs

Targeted genes

The three miRNAs which were detected to be significantly different between BB and BA samples were further investigated using 21 prediction tools, identifying 59 genes as predicted targets for at least two of these three miRNAs. For this selection, at least 15 tools presented a consistent result for all miRNA-interacted target genes. After an in-depth literature search, which was performed for each of these genes, nine were selected since they were previously identified and could possibly be associated with CVDs including ATAAs (Table 6).

Table 6
The selected predicted target genes
Gene namemiRNA
ABLIM1 Hsa-miR-140-5p
Hsa-miR-214-3p 
CCND2 Hsa-mir-191-5p
Hsa-miR-214-3p 
CRKL Hsa-mir-191-5p
Hsa-miR-214-3p
Hsa-miR-140-5p 
FNDC5 Hsa-miR-140-5p
Hsa-miR-214-3p 
HEY1 Hsa-miR-140-5p
Hsa-miR-191-5p 
MTMR4 Hsa-miR-140-5p
Hsa-miR-214-3p 
NFIA Hsa-miR-191-5p
Hsa-miR-214-3p 
PPP1CB Hsa-miR-140-5p
Hsa-miR-191-5p 
TJP1 Hsa-miR-140-5p
Hsa-miR-191-5p 
Gene namemiRNA
ABLIM1 Hsa-miR-140-5p
Hsa-miR-214-3p 
CCND2 Hsa-mir-191-5p
Hsa-miR-214-3p 
CRKL Hsa-mir-191-5p
Hsa-miR-214-3p
Hsa-miR-140-5p 
FNDC5 Hsa-miR-140-5p
Hsa-miR-214-3p 
HEY1 Hsa-miR-140-5p
Hsa-miR-191-5p 
MTMR4 Hsa-miR-140-5p
Hsa-miR-214-3p 
NFIA Hsa-miR-191-5p
Hsa-miR-214-3p 
PPP1CB Hsa-miR-140-5p
Hsa-miR-191-5p 
TJP1 Hsa-miR-140-5p
Hsa-miR-191-5p 

Hsa-miR-140-5p is predicted to target seven (ABLIM1, CRKL, FNDC5, HEY1, MTMR4, PPP1CB, TJP1) of these genes and hsa-miR-191-5p (CCND2, CRKL, HEY1, NFIA, PPP1CB, TJP1) as well as hsa-miR-214-3p (ABLIM1, CCND2, CRKL, FNDC5, MTMR4, NFIA) to target six from the total of nine genes. CRKL is the only gene which is predicted to be targeted by all three validated miRNAs.

Based on our results, all nine mRNAs examined (CCND2, CRKL, HEY1, ABLIM1, FNDC5, TJP1, MTMR4, NFIA and PPP1CB) showed higher expression in BB samples compared with BA samples with Myotubularin-related protein 4 (MTMR4) and Phosphatase 1 catalytic subunit β (PPP1CB) showing statistically significant differences with P-values 0.02 and 0.045, respectively (Figure 3). Further research using the miRTarBase showed strong experimental evidence associating CCND2 with hsa-miR-191 and CRKL with hsa-miR-214-3p. Figure 2B represents a flow chart summarizing secondary phase of the study.

Differentially expressed genes in plasma from ATAA patients
Figure 3
Differentially expressed genes in plasma from ATAA patients

Average ΔCt values of expressed genes (ABLIM1, CCND2, CRKL, FNDC5, HEY1, MTMR4, NFIA, PPP1CB and TJP1) in plasma samples from 11 patients collect before (BB) and after (BA) repair surgery. The values are lower in the BB compared with BA samples and so indicating higher expression levels in BB samples. The MTMR4 and the PPP1CB reached a statistical significance.

Figure 3
Differentially expressed genes in plasma from ATAA patients

Average ΔCt values of expressed genes (ABLIM1, CCND2, CRKL, FNDC5, HEY1, MTMR4, NFIA, PPP1CB and TJP1) in plasma samples from 11 patients collect before (BB) and after (BA) repair surgery. The values are lower in the BB compared with BA samples and so indicating higher expression levels in BB samples. The MTMR4 and the PPP1CB reached a statistical significance.

Close modal

Although a number of miRNAs have been associated with ATAAs [18,19,23,56], no biomarker(s) exists for the non-invasive, early detection or screening for the disease.

In the present study, from a total of 80 miRNAs identified, eight (hsa-let-7c-5p, hsa-let-7i-5p, hsa-miR-140-5p, hsa-miR-191-5p, hsa-miR-214-3p, hsa-miR-223-3p, hsa-miR-23a-3p, hsa-miR-30b-5p), were identified with statistically significant different expression levels (P<0.05) between plasma samples collected before (BB) and after (BA) repair surgery from patients with ATAAs. From these eight miRNAs, four miRNAs, were selected based on their predicted interaction with pathways associated with TAAs, and were further validated in a larger sample size of 22 plasma samples showing the same pattern of expression, lower expression in BB compared with BA samples. However, only three (hsa-miR-140-5p, hsa-miR-190-5p and hsa-miR-214-3p) out of four miRNAs showed significant differences at the statistical level between BB and BA samples. According to the literature search performed by the authors, and to the best of their knowledge, as yet neither of these three miRNAs have been previously described as biomarkers associated with ATAAs.

These miRNAs have been previously described as potential biomarkers in a number of cardiovascular as well as other diseases. More specifically, hsa-miR-140-5p, was detected to be down-regulated in heart failure [57] as well as in carcinomas like Wilm’s tumor [58] and non-small-cell lung cancer [59] and overexpressed in children with obesity [60].

The expression levels of hsa-miR-214-3p were previously identified to be up-regulated in serum samples from patients with atrial fibrillation [61] whereas in patients with sinonasal inverted papilloma [62] and Parkinson [63] it was down-regulated. In addition, hsa-miR-191-5p was also down-regulated in numerous diseases including, acute myocardial infraction [64], Duchenne muscular dystrophy [65], renal cell carcinoma [66] and significantly overexpressed in attention-deficit/hyperactivity disorder [67]. Furthermore, a study by Licholai and colleagues (2016) identified [23] an overexpression of hsa-miR-191-5p in aneurysmal tissue when compared with normal aortic tissue. MiRNAs participate in various cell processes by targeting several genes that participate in a wide range of cell signaling pathways. More specifically, a single miRNA may affect expression of several target genes and may induce changes in various processes and pathways presenting additional mechanisms by which a disease may be induced [68,69]. Therefore, miRNAs are involved in many pathological processes that are present in different diseases and they may show a cell-specific regulation. The value of any suggested biomarker could be possible targeted therapeutics through the common pathways that are involved in the different diseases.

For the selection of genes that were targets of at least two out of the three validated miRNAs, a total of 59 were predicted as targets. A literature search investigated the possible association of these genes with CVD including aortic aneurysms. This was followed by gene expression analysis revealing that the mRNA levels of two genes (MTMR4, and PPP1CB) were statistically significant (P=0.03 and P=0.045, respectively) as well as being inversely correlated with the corresponding miRNA expression levels in the plasma samples.

MTMR4 (target for hsa-miR-140-5p and hsa-miR-214-3p) and PPP1CB (target for hsa-miR-140-5p and hsa-miR-191-5p) indicate a possible association with Transforming Growth Factor (TGF)-β signaling pathway, which in term plays a significant role in TAAs [70,71]. In particular, MTMR4 and PPP1C interact and de-phosphorylate SMAD2/3 complex and TGF-β receptor of this pathway, respectively and thus control its effects downstream. TGF-β signaling pathway targets and controls the activation of different genes [72], where a dysregulation or an overexpression of this pathway has been associated with syndromes known to cause aortic aneurysms [72–74].

Furthermore, it has been associated with matrix metalloproteinase (MMP) expression which plays a role in collagen degradation, (collagen being a major component of the aortic wall), suggesting that an overexpression of this pathway could lead to the up-regulation of MMPs and thus TAA formation [75–77]. The corresponding target genes for the validated miRNAs could be further investigated for any possible diagnostic and/or therapeutic value.

To further confirm these results, a larger validation study including more BB and BA samples should be performed to assess the sensitivity and specificity of the selected miRNAs that could be used as biomarkers for the early identification of the disease. In addition, the collection of plasma samples at least 1 month after the surgery could be more representative as normal-state condition compared with those samples that were collected 5–7 days post-surgically. Control plasma samples will need to be included in the extended study.

Overall, the three miRNAs described in this pilot study may prove to be promising biomarkers associated with ATAAs. More studies are warranted for further exploration of their contribution to ATAA pathogenesis and their potential role as clinical biomarkers.

The authors declare that there are no competing interests associated with the manuscript.

This work was supported by The Cyprus School of Molecular Medicine; and The Cyprus Institute of Neurology and Genetics.

Conceived and designed the experiments: Marios A. Cariolou, Noam Shomron and Evy Bashiardes. Performed the experiments: Areti Moushi and Nir Pillar. Analyzed the data: Areti Moushi and Nir Pillar. Prepared the manuscript: Evy Bashiardes, Areti Moushi and Anna Keravnou. All authors reviewed and approved the final version of the manuscript.

The authors would like to thank all the patients who consented and participated in the present study.

AD

aortic dissection

ATAA

ascending thoracic aortic aneurysm

BA

blood samples collected after surgery

BB

blood samples collected before surgery

CVD

cardiovascular disease

miRNA

microRNA

MTMR4

myotubularin-related protein 4

NEC

non-enzyme control

NTC

non-template control

PPP1CB

Phosphatase 1 catalytic subunit β

RT-PCR

real time polymerase chain reaction

T

tissue sample

TAA

thoracic aortic aneurysm

TGF-β

transforming growth factor-β

UTR

untranslated region

1.
Hiratzka
L.F.
,
Bakris
G.L.
,
Beckman
J.A.
,
Bersin
R.M.
,
Carr
V.F.
,
Casey
D.E.
Jr
et al.
(
2010
)
2010 ACCF/AHA/AATS/ACR/ASA/SCA/SCAI/SIR/STS/SVM guidelines for the diagnosis and management of patients with Thoracic Aortic Disease: a report of the American College of Cardiology Foundation/American Heart Association Task Force on Practice Guidelines, American Association for Thoracic Surgery, American College of Radiology, American Stroke Association, Society of Cardiovascular Anesthesiologists, Society for Cardiovascular Angiography and Interventions, Society of Interventional Radiology, Society of Thoracic Surgeons, and Society for Vascular Medicine
.
Circulation
121
,
e266
e369
[PubMed]
2.
Isselbacher
E.M.
(
2005
)
Thoracic and abdominal aortic aneurysms
.
Circulation
111
,
816
828
[PubMed]
3.
Guo
D.C.
,
Papke
C.L.
,
He
R.
and
Milewicz
D.M.
(
2006
)
Pathogenesis of thoracic and abdominal aortic aneurysms
.
Ann. N.Y. Acad. Sci.
1085
,
339
352
[PubMed]
4.
Coady
M.A.
,
Rizzo
J.A.
,
Hammond
G.L.
,
Kopf
G.S.
and
Elefteriades
J.A.
(
1999
)
Surgical intervention criteria for thoracic aortic aneurysms: a study of growth rates and complications
.
Ann. Thorac. Surg.
67
,
1922
1926
,
[PubMed]
5.
System CfDCaPNVS
(
2017
)
LCWK1. Deaths, Percent of Total Deaths, and Death Rates for the 15 Leading Causes of Death in 5-year Age Groups, by Race and Sex: United States, 1999-2015
.
6.
Sidloff
D.A.
,
Stather
P.W.
,
Choke
E.
,
Bown
M.J.
and
Sayers
R.D.
(
2014
)
A systematic review and meta-analysis of the association between markers of hemostasis and abdominal aortic aneurysm presence and size
.
J. Vasc. Surg.
59
,
528.e4
535.e4
7.
Elefteriades
J.A.
(
2005
)
Beating a sudden killer
.
Sci. Am.
293
,
64
71
[PubMed]
8.
Hiratzka
L.F.
,
Creager
M.A.
,
Isselbacher
E.M.
,
Svensson
L.G.
,
Nishimura
R.A.
,
Bonow
R.O.
et al.
(
2016
)
Surgery for aortic dilatation in patients with bicuspid aortic valves: a statement of clarification from the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines
.
J. Am. Coll. Cardiol.
67
,
724
731
[PubMed]
9.
Svensson
L.G.
,
Adams
D.H.
,
Bonow
R.O.
,
Kouchoukos
N.T.
,
Miller
D.C.
,
O'Gara
P.T.
et al.
(
2013
)
Aortic valve and ascending aorta guidelines for management and quality measures
.
Ann. Thorac. Surg.
95
,
S1
S66
[PubMed]
10.
Chitwood
D.H.
and
Timmermans
M.C.
(
2007
)
Target mimics modulate miRNAs
.
Nat. Genet.
39
,
935
936
[PubMed]
11.
Houard
X.
,
Rouzet
F.
,
Touat
Z.
,
Philippe
M.
,
Dominguez
M.
,
Fontaine
V.
et al.
(
2007
)
Topology of the fibrinolytic system within the mural thrombus of human abdominal aortic aneurysms
.
J. Pathol.
212
,
20
28
[PubMed]
12.
van Rooij
E.
,
Purcell
A.L.
and
Levin
A.A.
(
2012
)
Developing microRNA therapeutics
.
Circ. Res.
110
,
496
507
[PubMed]
13.
Bronze-da-Rocha
E.
(
2014
)
MicroRNAs expression profiles in cardiovascular diseases
.
Biomed Res. Int.
2014
,
985408
[PubMed]
14.
Jin
H.
,
Li
C.
,
Ge
H.
,
Jiang
Y.
and
Li
Y.
(
2013
)
Circulating microRNA: a novel potential biomarker for early diagnosis of intracranial aneurysm rupture a case control study
.
J. Transl. Med.
11
,
296
[PubMed]
15.
Asaga
S.
,
Kuo
C.
,
Nguyen
T.
,
Terpenning
M.
,
Giuliano
A.E.
and
Hoon
D.S.
(
2011
)
Direct serum assay for microRNA-21 concentrations in early and advanced breast cancer
.
Clin. Chem.
57
,
84
91
[PubMed]
16.
Fichtlscherer
S.
,
De Rosa
S.
,
Fox
H.
,
Schwietz
T.
,
Fischer
A.
,
Liebetrau
C.
et al.
(
2010
)
Circulating microRNAs in patients with coronary artery disease
.
Circ. Res.
107
,
677
684
[PubMed]
17.
Boon
R.A.
,
Seeger
T.
,
Heydt
S.
,
Fischer
A.
,
Hergenreider
E.
,
Horrevoets
A.J.
et al.
(
2011
)
MicroRNA-29 in aortic dilation: implications for aneurysm formation
.
Circ. Res.
109
,
1115
1119
[PubMed]
18.
Jones
J.A.
,
Stroud
R.E.
,
O'Quinn
E.C.
,
Black
L.E.
,
Barth
J.L.
,
Elefteriades
J.A.
et al.
(
2011
)
Selective microRNA suppression in human thoracic aneurysms: relationship of miR-29a to aortic size and proteolytic induction
.
Circ. Cardiovasc. Genet.
4
,
605
613
[PubMed]
19.
Li
Y.
and
Maegdefessel
L.
(
2017
)
Non-coding RNA contribution to thoracic and abdominal aortic aneurysm disease development and progression
.
Front. Physiol.
8
,
429
[PubMed]
20.
Liao
M.
,
Zou
S.
,
Weng
J.
,
Hou
L.
,
Yang
L.
,
Zhao
Z.
et al.
(
2011
)
A microRNA profile comparison between thoracic aortic dissection and normal thoracic aorta indicates the potential role of microRNAs in contributing to thoracic aortic dissection pathogenesis
.
J. Vasc. Surg.
53
,
1341.e3
1349.e3
21.
Ikonomidis
J.S.
,
Ivey
C.R.
,
Wheeler
J.B.
,
Akerman
A.W.
,
Rice
A.
,
Patel
R.K.
et al.
(
2013
)
Plasma biomarkers for distinguishing etiologic subtypes of thoracic aortic aneurysm disease
.
J. Thorac. Cardiovasc. Surg.
145
,
1326
1333
[PubMed]
22.
Jondeau
G.
and
Boileau
C.
(
2012
)
Genetics of thoracic aortic aneurysms
.
Curr. Atheroscler. Rep.
14
,
219
226
[PubMed]
23.
Licholai
S.
,
Blaż
M.
,
Kapelak
B.
and
Sanak
M.
(
2016
)
Unbiased profile of microRNA expression in ascending aortic aneurysm tissue appoints molecular pathways contributing to the pathology
.
Ann. Thorac. Surg.
102
,
1245
1252
[PubMed]
24.
Moushi
A.
,
Michailidou
K.
,
Soteriou
M.
,
Cariolou
M.
and
Bashiardes
E.
(
2018
)
MicroRNAs as possible biomarkers for screening of aortic aneurysms: a systematic review and validation study
.
Biomarkers
23
,
253
264
25.
Patuzzo
C.
,
Pasquali
A.
,
Malerba
G.
,
Trabetti
E.
,
Pignatti
P.
,
Tessari
M.
et al.
(
2012
)
A preliminary microRNA analysis of non syndromic thoracic aortic aneurysms
.
Balkan J. Med. Genet.
15
,
51
55
[PubMed]
26.
Pritchard
C.C.
,
Kroh
E.
,
Wood
B.
,
Arroyo
J.D.
,
Dougherty
K.J.
,
Miyaji
M.M.
et al.
(
2012
)
Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies
.
Cancer Prev. Res.
5
,
492
497
[PubMed]
27.
Love
M.I.
,
Huber
W.
and
Anders
S.
(
2014
)
Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2
.
Genome Biol.
15
,
550
[PubMed]
28.
Pillar
N.
,
Haguel
D.
,
Grad
M.
,
Shapira
G.
,
Yoffe
L.
and
Shomron
N.
(
2019
)
Characterization of microRNA and gene expression profiles following ricin intoxication
.
Toxins
11
,
250
[PubMed]
29.
Amarilyo
G.
,
Pillar
N.
,
Ben-Zvi
I.
,
Weissglas-Volkov
D.
,
Zalcman
J.
,
Harel
L.
et al.
(
2018
)
Analysis of microRNAs in familial Mediterranean fever
.
PLoS ONE
13
,
e0197829
[PubMed]
30.
Pillar
N.
,
Bairey
O.
,
Goldschmidt
N.
,
Fellig
Y.
,
Rosenblat
Y.
,
Shehtman
I.
et al.
(
2017
)
MicroRNAs as predictors for CNS relapse of systemic diffuse large B-cell lymphoma
.
Oncotarget
8
,
86020
86030
[PubMed]
31.
Vlachos
I.S.
and
Hatzigeorgiou
A.G.
(
2017
)
Functional analysis of miRNAs using the DIANA tools online suite
.
Methods Mol. Biol.
1517
,
25
50
32.
Vlachos
I.S.
,
Zagganas
K.
,
Paraskevopoulou
M.D.
,
Georgakilas
G.
,
Karagkouni
D.
,
Vergoulis
T.
et al.
(
2015
)
DIANA-miRPath v3.0: deciphering microRNA function with experimental support
.
Nucleic Acids Res.
43
,
W460
W466
[PubMed]
33.
Yue
D.
,
Guo
M.
,
Chen
Y.
and
Huang
Y.
(
2012
)
A Bayesian decision fusion approach for microRNA target prediction
.
BMC Genomics
13
,
S13
[PubMed]
34.
Veksler-Lublinsky
I.
,
Shemer-Avni
Y.
,
Kedem
K.
and
Ziv-Ukelson
M.
(
2010
)
Gene bi-targeting by viral and human miRNAs
.
BMC Bioinformatics
11
,
249
[PubMed]
35.
Gennarino
V.A.
,
D'Angelo
G.
,
Dharmalingam
G.
,
Fernandez
S.
,
Russolillo
G.
,
Sanges
R.
et al.
(
2012
)
Identification of microRNA-regulated gene networks by expression analysis of target genes
.
Genome Res.
22
,
1163
1172
[PubMed]
36.
Chiu
H.S.
,
Llobet-Navas
D.
,
Yang
X.
,
Chung
W.J.
,
Ambesi-Impiombato
A.
,
Iyer
A.
et al.
(
2015
)
Cupid: simultaneous reconstruction of microRNA-target and ceRNA networks
.
Genome Res.
25
,
257
267
[PubMed]
37.
Gaidatzis
D.
,
van Nimwegen
E.
,
Hausser
J.
and
Zavolan
M.
(
2007
)
Inference of miRNA targets using evolutionary conservation and pathway analysis
.
BMC Bioinformatics
8
,
69
[PubMed]
38.
Bandyopadhyay
S.
,
Ghosh
D.
,
Mitra
R.
and
Zhao
Z.
(
2015
)
MBSTAR: multiple instance learning for predicting specific functional binding sites in microRNA targets
.
Sci. Rep.
5
,
8004
[PubMed]
39.
Betel
D.
,
Wilson
M.
,
Gabow
A.
,
Marks
D.S.
and
Sander
C.
(
2008
)
The microRNA.org resource: targets and expression
.
Nucleic Acids Res.
36
,
D149
D153
[PubMed]
40.
Leclercq
M.
,
Diallo
A.B.
and
Blanchette
M.
(
2013
)
Computational prediction of the localization of microRNAs within their pre-miRNA
.
Nucleic Acids Res.
41
,
7200
7211
[PubMed]
41.
Leclercq
M.
,
Diallo
A.B.
and
Blanchette
M.
(
2017
)
Prediction of human miRNA target genes using computationally reconstructed ancestral mammalian sequences
.
Nucleic Acids Res.
45
,
556
566
[PubMed]
42.
Jeggari
A.
,
Marks
D.S.
and
Larsson
E.
(
2012
)
miRcode: a map of putative microRNA target sites in the long non-coding transcriptome
.
Bioinformatics
28
,
2062
2063
[PubMed]
43.
Giles
C.B.
,
Girija-Devi
R.
,
Dozmorov
M.G.
and
Wren
J.D.
(
2013
)
mirCoX: a database of miRNA-mRNA expression correlations derived from RNA-seq meta-analysis
.
BMC Bioinformatics
14
,
S17
[PubMed]
44.
Chen
Y.
and
Wang
X.
(
2020
)
miRDB: an online database for prediction of functional microRNA targets
.
Nucleic Acids Res.
48
,
D127
D131
[PubMed]
45.
Vejnar
C.E.
,
Blum
M.
and
Zdobnov
E.M.
(
2013
)
miRmap web: Comprehensive microRNA target prediction online
.
Nucleic Acids Res.
41
,
W165
W168
[PubMed]
46.
Hsu
J.B.
,
Chiu
C.M.
,
Hsu
S.D.
,
Huang
W.Y.
,
Chien
C.H.
,
Lee
T.Y.
et al.
(
2011
)
miRTar: an integrated system for identifying miRNA-target interactions in human
.
BMC Bioinformatics
12
,
300
[PubMed]
47.
Ahadi
A.
,
Sablok
G.
and
Hutvagner
G.
(
2017
)
miRTar2GO: a novel rule-based model learning method for cell line specific microRNA target prediction that integrates Ago2 CLIP-Seq and validated microRNA-target interaction data
.
Nucleic Acids Res.
45
,
e42
[PubMed]
48.
Gumienny
R.
and
Zavolan
M.
(
2015
)
Accurate transcriptome-wide prediction of microRNA targets and small interfering RNA off-targets with MIRZA-G
.
Nucleic Acids Res.
43
,
1380
1391
[PubMed]
49.
Mitra
R.
and
Bandyopadhyay
S.
(
2011
)
MultiMiTar: a novel multi objective optimization based miRNA-target prediction method
.
PLoS ONE
6
,
e24583
[PubMed]
50.
Chen
K.
and
Rajewsky
N.
(
2006
)
Natural selection on human microRNA binding sites inferred from SNP data
.
Nat. Genet.
38
,
1452
1456
[PubMed]
51.
Elefant
N.
,
Berger
A.
,
Shein
H.
,
Hofree
M.
,
Margalit
H.
and
Altuvia
Y.
(
2011
)
RepTar: a database of predicted cellular targets of host and viral miRNAs
.
Nucleic Acids Res.
39
,
D188
D194
[PubMed]
52.
Krüger
J.
and
Rehmsmeier
M.
(
2006
)
RNAhybrid: microRNA target prediction easy, fast and flexible
.
Nucleic Acids Res.
34
,
W451
W454
[PubMed]
53.
Nielsen
C.B.
,
Shomron
N.
,
Sandberg
R.
,
Hornstein
E.
,
Kitzman
J.
and
Burge
C.B.
(
2007
)
Determinants of targeting by endogenous and exogenous microRNAs and siRNAs
.
RNA
13
,
1894
1910
54.
UniProt Consortium T.
(
2018
)
UniProt: the universal protein knowledgebase
.
Nucleic Acids Res.
46
,
2699
[PubMed]
55.
Chou
C.H.
,
Shrestha
S.
,
Yang
C.D.
,
Chang
N.W.
,
Lin
Y.L.
,
Liao
K.W.
et al.
(
2018
)
miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions
.
Nucleic Acids Res.
46
,
D296
D302
[PubMed]
56.
Wang
X.J.
,
Huang
B.
,
Yang
Y.M.
,
Zhang
L.
,
Su
W.J.
,
Tian
L.
et al.
(
2015
)
Differential expression of microRNAs in aortic tissue and plasma in patients with acute aortic dissection
.
J. Geriatr. Cardiol.
12
,
655
661
[PubMed]
57.
Marques
F.Z.
,
Vizi
D.
,
Khammy
O.
,
Mariani
J.A.
and
Kaye
D.M.
(
2016
)
The transcardiac gradient of cardio-microRNAs in the failing heart
.
Eur. J. Heart Fail.
18
,
1000
1008
[PubMed]
58.
Liu
Z.
,
He
F.
,
OuYang
S.
,
Li
Y.
,
Ma
F.
,
Chang
H.
et al.
(
2019
)
miR-140-5p could suppress tumor proliferation and progression by targeting TGFBRI/SMAD2/3 and IGF-1R/AKT signaling pathways in Wilms’ tumor
.
BMC Cancer
19
,
405
[PubMed]
59.
Wang
K.
,
Chen
M.
and
Wu
W.
(
2017
)
Analysis of microRNA (miRNA) expression profiles reveals 11 key biomarkers associated with non-small cell lung cancer
.
World J. Surg. Oncol.
15
,
175
[PubMed]
60.
Oses
M.
,
Margareto Sanchez
J.
,
Portillo
M.P.
,
Aguilera
C.M.
and
Labayen
I.
(
2019
)
Circulating miRNAs as biomarkers of obesity and obesity-associated comorbidities in children and adolescents: a systematic review
.
Nutrients
11
,
[PubMed]
61.
Natsume
Y.
,
Oaku
K.
,
Takahashi
K.
,
Nakamura
W.
,
Oono
A.
,
Hamada
S.
et al.
(
2018
)
Combined analysis of human and experimental murine samples identified novel circulating microRNAs as biomarkers for atrial fibrillation
.
Circ. J.
82
,
965
973
[PubMed]
62.
Teng
Y.
,
Li
Y.
,
Lin
Z.
,
Gao
Y.
,
Cao
X.
,
Lou
X.
et al.
(
2018
)
Analysis of miRNA expression profiling identifies miR-214-3p as a novel biomarker in sinonasal inverted papilloma
.
Epigenomics
10
,
1541
1553
[PubMed]
63.
Schulz
J.
,
Takousis
P.
,
Wohlers
I.
,
Itua
I.O.G.
,
Dobricic
V.
,
Rücker
G.
et al.
(
2019
)
Meta-analyses identify differentially expressed micrornas in Parkinson’s disease
.
Ann. Neurol.
85
,
835
851
[PubMed]
64.
Hsu
A.
,
Chen
S.J.
,
Chang
Y.S.
,
Chen
H.C.
and
Chu
P.H.
(
2014
)
Systemic approach to identify serum microRNAs as potential biomarkers for acute myocardial infarction
.
Biomed Res. Int.
2014
,
418628
[PubMed]
65.
Mousa
N.O.
,
Abdellatif
A.
,
Fahmy
N.
,
Zada
S.
,
El-Fawal
H.
and
Osman
A.
(
2020
)
Circulating MicroRNAs in Duchenne muscular dystrophy
.
Clin. Neurol. Neurosurg.
189
,
105634
[PubMed]
66.
Chen
P.
,
Pan
X.
,
Zhao
L.
,
Jin
L.
,
Lin
C.
,
Quan
J.
et al.
(
2018
)
MicroRNA-191-5p exerts a tumor suppressive role in renal cell carcinoma
.
Exp. Ther. Med.
15
,
1686
1693
[PubMed]
67.
Sánchez-Mora
C.
,
Soler Artigas
M.
,
Garcia-Martínez
I.
,
Pagerols
M.
,
Rovira
P.
,
Richarte
V.
et al.
(
2019
)
Epigenetic signature for attention-deficit/hyperactivity disorder: identification of miR-26b-5p, miR-185-5p, and miR-191-5p as potential biomarkers in peripheral blood mononuclear cells
.
Neuropsychopharmacology
44
,
890
897
[PubMed]
68.
Boileau
A.
,
Lino Cardenas
C.L.
,
Courtois
A.
,
Zhang
L.
,
Rodosthenous
R.S.
,
Das
S.
et al.
(
2019
)
MiR-574-5p: a circulating marker of thoracic aortic aneurysm
.
Int. J. Mol. Sci.
20
,
3924
[PubMed]
69.
Gasiulė
S.
,
Stankevičius
V.
,
Patamsytė
V.
,
Ražanskas
R.
,
Žukovas
G.
,
Kapustina
Ž
et al.
(
2019
)
Tissue-specific miRNAs regulate the development of thoracic aortic aneurysm: the emerging role of KLF4 network
.
J. Clin. Med.
8
,
70.
Korrodi-Gregório
L.
,
Silva
J.V.
,
Santos-Sousa
L.
,
Freitas
M.J.
,
Felgueiras
J.
and
Fardilha
M.
(
2014
)
TGF-β cascade regulation by PPP1 and its interactors -impact on prostate cancer development and therapy
.
J. Cell. Mol. Med.
18
,
555
567
[PubMed]
71.
Yu
J.
,
Pan
L.
,
Qin
X.
,
Chen
H.
,
Xu
Y.
,
Chen
Y.
et al.
(
2010
)
MTMR4 attenuates transforming growth factor beta (TGFbeta) signaling by dephosphorylating R-Smads in endosomes
.
J. Biol. Chem.
285
,
8454
8462
[PubMed]
72.
Jones
J.A.
,
Spinale
F.G.
and
Ikonomidis
J.S.
(
2009
)
Transforming growth factor-beta signaling in thoracic aortic aneurysm development: a paradox in pathogenesis
.
J. Vasc. Res.
46
,
119
137
[PubMed]
73.
Loeys
B.L.
,
Chen
J.
,
Neptune
E.R.
,
Judge
D.P.
,
Podowski
M.
,
Holm
T.
et al.
(
2005
)
A syndrome of altered cardiovascular, craniofacial, neurocognitive and skeletal development caused by mutations in TGFBR1 or TGFBR2
.
Nat. Genet.
37
,
275
281
[PubMed]
74.
Yuan
S.M.
,
Ma
H.H.
,
Zhang
R.S.
and
Jing
H.
(
2011
)
Transforming growth factor-β signaling pathway in Marfan’s syndrome: a preliminary histopathological study
.
VASA Zeitschrift Fur Gefasskrankheiten
40
,
369
374
[PubMed]
75.
Jabłońska-Trypuć
A.
,
Matejczyk
M.
and
Rosochacki
S.
(
2016
)
Matrix metalloproteinases (MMPs), the main extracellular matrix (ECM) enzymes in collagen degradation, as a target for anticancer drugs
.
J. Enzyme Inhib. Med. Chem.
31
,
177
183
[PubMed]
76.
Konrad
L.
,
Scheiber
J.A.
,
Schwarz
L.
,
Schrader
A.J.
and
Hofmann
R.
(
2009
)
TGF-beta1 and TGF-beta2 strongly enhance the secretion of plasminogen activator inhibitor-1 and matrix metalloproteinase-9 of the human prostate cancer cell line PC-3
.
Regul. Pept.
155
,
28
32
[PubMed]
77.
Krstić
J.
,
Trivanović
D.
,
Mojsilović
S.
and
Santibanez
J.F.
(
2015
)
Transforming growth factor-beta and oxidative stress interplay: implications in tumorigenesis and cancer progression
.
Oxid. Med. Cell. Longev.
2015
,
654594
[PubMed]

Author notes

*

These authors contributed equally to this work.

This is an open access article published by Portland Press Limited on behalf of the Biochemical Society and distributed under the Creative Commons Attribution License 4.0 (CC BY).

Supplementary data