Next Article in Journal
Hepatitis B Virus and Tuberculosis Are Associated with Increased Noncommunicable Disease Risk among Treatment-Naïve People with HIV: Opportunities for Prevention, Early Detection and Management of Comorbidities in Sierra Leone
Previous Article in Journal
Procedural Outcome Following Stent-Assisted Coiling for Wide-Necked Aneurysms Using Three Different Stent Models: A Single-Center Experience
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCM
Volume 11
Issue 12
10.3390/jcm11123470
jcm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

miRNA Regulatory Networks Associated with Peripheral Vascular Diseases

by
Daniel P. Zalewski
1,*,†,
Karol P. Ruszel
2,†,
Andrzej Stępniewski
3,
Dariusz Gałkowski
4,
Marcin Feldo
5,‡,
Janusz Kocki
2,‡ and
Anna Bogucka-Kocka
1,‡
1
Chair and Department of Biology and Genetics, Medical University of Lublin, 4a Chodźki St., 20-093 Lublin, Poland
2
Department of Clinical Genetics, Chair of Medical Genetics, Medical University of Lublin, 11 Radziwiłłowska St., 20-080 Lublin, Poland
3
Ecotech Complex Analytical and Programme Centre for Advanced Environmentally Friendly Technologies, University of Marie Curie-Skłodowska, 39 Głęboka St., 20-612 Lublin, Poland
4
Department of Pathology and Laboratory Medicine, Rutgers-Robert Wood Johnson Medical School, One Robert Wood Johnson Place, New Brunswick, NJ 08903-0019, USA
5
Chair and Department of Vascular Surgery and Angiology, Medical University of Lublin, 11 Staszica St., 20-081 Lublin, Poland
*
Author to whom correspondence should be addressed.
†
These authors contributed to this paper as first authors.
‡
These authors contributed to this paper as senior authors.
J. Clin. Med. 2022, 11(12), 3470; https://doi.org/10.3390/jcm11123470
Submission received: 29 March 2022 / Revised: 13 May 2022 / Accepted: 14 June 2022 / Published: 16 June 2022
(This article belongs to the Section Clinical Laboratory Medicine)
Browse Figures
Versions Notes

Abstract

:
A growing body of evidence indicates a crucial role of miRNA regulatory function in a variety of mechanisms that contribute to the development of diseases. In our previous work, alterations in miRNA expression levels and targeted genes were shown in peripheral blood mononuclear cells (PBMCs) from patients with lower extremity artery disease (LEAD), abdominal aortic aneurysm (AAA), and chronic venous disease (CVD) in comparison with healthy controls. In this paper, previously obtained miRNA expression profiles were compared between the LEAD, AAA, and CVD groups to find either similarities or differences within the studied diseases. Differentially expressed miRNAs were identified using the DESeq2 method implemented in the R programming software. Pairwise comparisons (LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD) were performed and revealed 10, 8, and 17 differentially expressed miRNA transcripts, respectively. The functional analysis of the obtained miRNAs was conducted using the miRNet 2.0 online tool and disclosed associations with inflammation and cellular differentiation, motility, and death. The miRNet 2.0 tool was also used to identify regulatory interactions between dysregulated miRNAs and target genes in patients with LEAD, AAA, and CVD. The presented research provides new information about similarities and differences in the miRNA-dependent regulatory mechanisms involved in the pathogenesis of LEAD, AAA, and CVD.

1. Introduction

Permanent progress in the development of research methods allows the elucidation of molecular mechanisms involved in the pathogenesis of human diseases. One of the most intensively studied aspects in this field is the role of the microRNA (miRNA) regulatory function in the processes underlying pathological conditions.
miRNAs constitute a group of non-coding and single-stranded RNA molecules with a length of 18–25 nucleotides. The fundamental role of miRNA is post-transcriptional regulation of gene expression through direct binding to messenger RNA (mRNA) strands, causing mRNA degradation or inhibition of translation [1,2]. Many previous studies have indicated altered expression of miRNAs in patients with various diseases, including cancer, autoimmunology, infections, preeclampsia, and cardiovascular diseases, indicating a potential utility of miRNAs in diagnosis and therapy [2,3,4,5].
Our research group is focused on the effect of alterations in the miRNA–gene regulatory network on the pathogenesis of three vascular diseases: lower extremity artery disease (LEAD), abdominal aortic aneurysm (AAA), and chronic venous disease (CVD). LEAD is caused by the chronic development of atherosclerotic plaques in the arteries of the lower extremities, leading to ischemic events, including intermittent claudication and critical limb ischemia [6,7,8,9,10]. The development of AAA is related to a focal dilatation of the abdominal aorta measuring 50% greater than the normal proximal segment or >3 cm in maximum diameter. The progression of AAA can lead to aneurysm rupture, an event that reaches 60–80% of mortality [11,12,13,14,15]. In turn, CVD is defined as a syndrome of chronic morphological and functional abnormalities in the venous circulation of the lower extremities, resulting from venous valve incompetence, venous reflux, and venous hypertension [16,17,18,19,20,21,22,23,24]. Despite the different clinical images of these diseases, they share many common pathological conditions, including endothelial dysfunction, inflammatory response, and vascular remodeling [6,7,11,12,16,17,20]. Furthermore, the high prevalence, multifactorial character, and burden of severe life-threatening consequences of these diseases made them major global health problems [6,8,9,10,13,14,15,18,19,22,24]. Therefore, intensive research is desired to propose new and more effective strategies for the diagnosis and therapy of LEAD, AAA, and CVD.
Many studies have indicated an important role of miRNAs in vascular pathology [25,26,27,28] and have pointed to certain dysregulated miRNAs as candidates for biomarkers of LEAD [29,30], AAA [31], and CVD [32,33,34]. Special attention is paid to the identification of circulatory miRNA biomarkers (from whole blood, blood cells, plasma, or serum) with potential utility as diagnostic, prognostic, and therapeutic targets for these diseases [35,36,37,38,39,40].
In our previous studies, we presented an integrated analysis of miRNA and gene expression in peripheral blood mononuclear cells (PBMCs) of patients with LEAD, AAA, and CVD in comparison with control subjects to identify alterations in miRNA regulatory function associated with these diseases [40,41,42]. However, it is unknown whether the miRNA expression profiles differ between the analyzed diseases. Therefore, in this paper, we present the results of the next step of our study, which includes a pairwise comparative analysis among the LEAD, AAA, and CVD groups, together with the functional analysis and identification of potential targets of dysregulated miRNAs. Identification of similarities and differences in the expression of the miRNAs in the regulatory network could provide a deeper insight into the vascular pathology underlying these diseases. Identification of disease-specific biomarkers could potentially enable the development of new methods for diagnosis and treatment.

2. Materials and Methods

In our previous studies, miRNA expression profiles of PBMCs were compared between patients with vascular diseases (LEAD, AAA, and CVD) and healthy controls [40,41,42]. For a deeper characterization of the molecular mechanisms involved in the pathogenesis of LEAD, AAA, and CVD, we performed and presented in this paper the next step of the study, where miRNA expression profiles were compared between the studied diseases using the following comparisons: LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD.

2.1. Study Participants

The study was carried out in accordance with the Declaration of Helsinki and after the approval of the Ethics Committee of the Medical University of Lublin (decision No. KE-0254/341/2015). The study group consisted of 40 patients with LEAD, 28 patients with AAA, and 34 patients with CVD. The clinical examination of the patients was carried out at the Independent Public Clinical Hospital No. 1 in Lublin. Informed and signed consent was obtained from all study subjects.
Patients with LEAD were evaluated using peripheral pulse examination, ankle-brachial index test, treadmill test, angiography, and duplex color flow ultrasound scanning. Individuals diagnosed with intermittent ischemia-associated claudication without critical ischemia (Trans-Atlantic Inter-Society Consensus score B or C and Rutherford category 2 or 3) and with atherosclerotic lesions localized in the femoral, iliac, or popliteal arteries were included. Patients with AAA were diagnosed using duplex ultrasonography and contrast-enhanced spiral computed tomography with volume-rendered reconstructions. AAA patients with a maximum diameter of the aneurysm ranging from 5.6 to 7.8 cm (mean = 6.39, standard deviation = 0.633) were included. CVD patients were examined using tourniquet test, auscultation, and duplex ultrasound scanning, and those who had symptoms classified according to CEAP as varicose veins (C2) of superficial veins (As) with primary etiology (Ep) and reflux pathophysiology (Pr) were included. Established exclusion criteria common to all studied diseases include type 1 diabetes mellitus, previous vascular surgery, inflammatory diseases, and pregnancy. A detailed presentation of the clinical characteristics of the participants and a full description of the inclusion and exclusion criteria are provided in our previous articles [40,41,42].

2.2. miRNA Expression Datasets

The miRNA expression datasets were generated by RNA sequencing of the PBMC samples obtained from the study participants as described in our previous papers [40,41,42]. Briefly, PBMC specimens were isolated from whole blood samples using density gradient centrifugation with Gradisol L reagent (Aqua-Med, Łódź, Poland). A diversity of white blood cell subpopulations in the studied groups was evaluated using laboratory analysis of blood morphology (Figure S1). Small RNA fractions were isolated from PBMC samples using the MirVana microRNA Isolation Kit (Ambion, Austin, TX, USA). The quantity and quality of the obtained small RNA samples were assessed using the Agilent 2100 Bioanalyzer (Agilent Small RNA Kit, Agilent Technologies, Santa Clara, CA, USA). Small RNA samples were subjected to preparation of miRNA libraries using the Ion Total RNA-Seq Kit v2, Magnetic Bead Cleanup Module kit, and Ion Xpress RNA-Seq Barcode 01-16 Kit (Life Technologies, Carlsbad, CA, USA). Libraries were sequenced on the Ion 540 chips (Life Technologies) using the Ion S5 XL System (Thermo Fisher Scientific, Waltham, MA, USA). Raw sequences were aligned to 2792 human miRNAs from miRBase v21 (http://www.mirbase.org) using the Torrent Suite Software v5.0.4. and the Ion Torrent Small RNA Plugin v5.0.5r3 plugin (Thermo Fisher Scientific, Waltham, MA, USA). The parameters describing miRNA libraries and the primary results of the sequencing data analysis are provided in the supplementary materials provided together with our previous publications [40,41,42].

2.3. Data Analysis

Data analysis was performed using the R environment (version 4.2.0, https://www.r-project.org) and the appropriate packages according to the corresponding reference manuals.
The statistical significance of the differences in demographic and clinical parameters between the LEAD, AAA, and CVD groups was examined using the Kruskal–Wallis rank sum test for continuous variables (kruskal.test function in R) and the two-sided Fisher’s exact test for categorical variables (fisher.test function in R).
All other statistical procedures applied to expression datasets and subsequent bioinformatic analysis were previously described in detail in [40,41,42]. The uniformity of the expression data was confirmed on the boxplot of Cook’s distances of miRNAs across all analyzed samples (Figure S2). The quality of differential expression analysis results was assessed on MA plots and histograms of the p-values obtained for each comparison (Figures S3–S5). Differential expression analysis of expression datasets was performed on biological replicates using the DESeq2 method implemented in the DESeq2 1.36.0 package [43] (https://bioconductor.org/packages/release/bioc/html/DESeq2.html (accessed on 28 March 2022)).
DESeq2 analysis was performed on expression data filtered out of genes with a mean of reads number lower than 1. The analysis was carried out in multivariate mode, using disease status as a variable of interest and the following characteristics as covariates: age; body mass index; morphology test results (hematocrit and counts of monocytes, eosinophils, basophils, and platelets); creatinine levels; sex; smoking (never and former smokers vs. current smokers); hypertension; type 2 diabetes; coronary artery disease; and medication with statins, acetylsalicylic acid, clopidogrel, beta-adrenergic blockers, angiotensin-converting enzyme inhibitors, Ca2+ channel blockers, fibrates, micronized diosmin, and preparation containing hesperidin, Ruscus aculeatus, and vitamin C. Differentially expressed miRNAs received from the DESeq2 analysis with p-value below 0.05, absolute log2 fold change > 1, and mean of normalized counts >10 were selected.
Venn diagrams and heatmap with hierarchical clustering were created using VennDiagram 1.7.3 (https://cran.r-project.org/web/packages/VennDiagram/index.html (accessed on 28 March 2022)) [44] and pheatmap 1.0.12 (https://cran.r-project.org/web/packages/pheatmap/index.html (accessed on 28 March 2022)) packages, respectively.
A receiver operating characteristic (ROC) analysis implemented in the pROC 1.18.0 package [45] (https://cran.r-project.org/web/packages/pROC/index.html (accessed on 28 March 2022)) was used to evaluate the classification performance of selected miRNAs. For ROC analysis and visualizations, filtered expression data were transformed using regularized log normalization (the rlog function in DESeq2 package).
The biological role of selected miRNAs was explored using the miRNet 2.0 online platform [46] (https://www.mirnet.ca) using a hypergeometric test applied to functional terms of the ‘miRNA-Function’ category. The miRNet 2.0 database was also used to identify interactions between miRNAs selected in this study and genes selected in the previous study [47]. Obtained interactions formed a regulatory network, which was presented using Cytoscape v3.9.1 software (https://cytoscape.org/) [48].
The identification of relationships between the characteristics of the study participants and the expression of selected genes was performed using the Spearman rank correlation test implemented in the Hmisc 4.7-0. package (https://cran.r-project.org/web/packages/Hmisc/index.html (accessed on 28 March 2022)) as well as a two-sided Mann–Whitney U test implemented in wilcox.test function in R. Multivariate linear regression models were constructed using the lm() base function in R, and obtained results were summarized using sjPlot 2.8.10 package (https://cran.r-project.org/web/packages/sjPlot/index.html (accessed on 28 March 2022)).

3. Results

3.1. The Study Group Characteristics

The study group included 40 patients with LEAD, 28 patients with AAA, and 34 patients with CVD. Information on the characteristics of the participants was provided separately for each disease in our previous papers [40,41,42], and here it is summarized in Table 1. Statistically significant differences (p < 0.05) in some characteristics were disclosed between the LEAD, AAA, and CVD groups, especially in age, hypertension status, blood cell counts, creatinine serum level, and medication (Table 1). The observed differences were difficult to avoid during the construction of the study population because they reflect different risk factors, comorbidities, and medications related to particular diseases. The relationships between these factors and the expression of selected miRNAs are investigated in Section 3.5.

3.2. Identification of Differentially Expressed miRNAs between the LEAD, AAA, and CVD Groups

The miRNA expression profiles were compared between the studied diseases using the following comparisons: LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD. The multivariate DESeq2 method was applied to select the most relevant differentially expressed miRNAs using disease status as a variable of interest and characteristics differentiating the study groups with a statistical significance of p < 0.05 (see Table 1). Application of established selection criteria allowed the selection of 10 miRNA transcripts (8 upregulated and 2 downregulated) from the LEAD vs. AAA comparison, 8 miRNA transcripts (4 upregulated and 4 downregulated) from the LEAD vs. CVD comparison, and 17 miRNA transcripts (9 upregulated and 8 downregulated) from the AAA vs. CVD comparisons (Table 2). The detailed results of the differential expression analysis concerning all miRNA transcripts are provided in Tables S1–S3.
To identify similarities and differences between the sets of miRNA transcripts selected from the performed comparisons, the sets of 10, 8, and 17 miRNA transcripts were compared in a Venn diagram (Figure 1A). Six miRNAs were common for the LEAD vs. CVD and AAA vs. CVD comparisons, and one miRNA, hsa-miR-135b-5p, was common for the LEAD vs. CVD and AAA vs. CVD comparisons (Figure 1A). The expression of selected miRNA transcripts was visualized on a heatmap with hierarchical clustering (Figure 1B).
The classification performance of selected miRNA transcripts was further evaluated using ROC analysis, which shows moderate classification ability (areas under the ROC curves ranged from 0.824 to 0.450) (Table 2). The detailed results of the ROC analysis are provided in Figures S6–S8 and Table S4.

3.3. Selected miRNAs Are Closely Related to Vascular Pathology

The biological role of dysregulated miRNAs was explored using the miRNet 2.0 online platform (https://www.mirnet.ca). This tool requires mature miRNA IDs; therefore, the names of the miRNA transcripts were transformed into the corresponding mature miRNA IDs according to miRBase 22.1 (http://www.mirbase.org). Thus, functional analysis was performed for 7 unique mature miRNAs given by 10 miRNA transcripts selected from the LEAD vs. AAA comparison, 8 unique mature miRNAs given by 8 miRNA transcripts selected from the LEAD vs. CVD comparison, and 16 unique mature miRNAs given by 17 miRNA transcripts selected from the AAA vs. CVD comparison (Table 2). The functional analysis was performed using a hypergeometric test applied to functional terms of the ‘miRNA-Function’ category in the miRNet 2.0 tool. The top 15 most enriched functional terms (with the lowest p-value of enrichment) and associated miRNAs are presented in Figure 2 and Table S5. Analyzed miRNAs were associated with processes closely related to vascular pathology, including inflammation and cell differentiation, motility, and death (Figure 2).

3.4. Selected miRNAs Regulate Genes Differentially Expressed between the LEAD, AAA, and CVD Groups

To construct a miRNA–gene regulatory network for selected miRNAs, the results of this study were integrated with the results of our previous paper, in which the whole transcriptome expression profiles of LEAD, AAA, and CVD patients were analyzed and the sets of 21, 58, and 10 differentially expressed genes were revealed from the LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively [47]. The altered expression of these genes could be an effect of the dysregulation of miRNAs found in the current study. To examine this effect, the sets of miRNAs selected in this study were subjected to the miRNet 2.0 tool to identify the targeted genes, and then the received sets of targets were searched for differentially expressed genes found in our previous work.
The identification of targets for seven miRNAs selected from the LEAD vs. AAA comparison revealed a set of 7359 genes, containing 3 genes (GIT2, UFM1, and YBX1) previously shown to be dysregulated between the LEAD and AAA groups. The identification of targets for eight miRNAs selected from the LEAD vs. CVD comparison disclosed a group of 5980 genes, containing 8 genes (ARL6IP1, C1orf216, DNAH1, EIF3C, FAM167A, HECTD4, PSME1, and TSC2) previously shown to be dysregulated between the LEAD and CVD groups. The identification of targets for 16 miRNAs selected from the AAA vs. CVD comparison disclosed a set of 8006 genes, containing 2 genes (MALT1 and STMN3) previously shown to be dysregulated between the AAA and CVD groups. Detailed information regarding identified interactions is provided in Table S6.
The interactions found between miRNAs selected in this study and genes reported in the previous study were visualized on regulatory networks, including four unique miRNA–gene pairs between two miRNAs and three genes resulting from the LEAD vs. AAA comparison, eight unique miRNA–gene pairs between three miRNAs and eight genes resulting from the LEAD vs. CVD comparison, and three miRNA–gene pairs between three miRNAs and two genes resulting from the AAA vs. CVD comparison (Figure 3).

3.5. Some of the Selected miRNAs Are Associated with the Clinical and Demographic Characteristics of the Study Groups

A correlation analysis was performed between the expression data of the miRNA transcripts selected in this study and such patients’ characteristics as age, body mass index, and morphology test results (hematocrit and counts of red blood cells, leukocytes, neutrophils, lymphocytes, monocytes, eosinophils, basophils, and platelets), as well as hemoglobin, creatinine, and urea levels. There were no strong or moderate correlations, and only a weak negative correlation (0.4 < |R| < 0.5) was determined between age and expression of has-miR-31-3p (R = −0.42, Benjamini–Hochberg FDR = 0.00026). All correlation results are presented in Table S7.
The two-sided Mann–Whitney U test was used to identify relationships between the expression of selected miRNA transcripts and the following categorical variables: sex; smoking (never and former smokers vs. current smokers); hypertension; type 2 diabetes; coronary artery disease; myocardial infarction; and medication with statins, acetylsalicylic acid, clopidogrel, beta-adrenergic blockers, angiotensin-converting enzyme inhibitors, Ca2+ channel blockers, fibrates, and micronized diosmin, as well as preparation with hesperidin, Ruscus aculeatus, and vitamin C (refer to Table 1). The relationships revealed with a statistical significance of p < 0.05 corrected by the Benjamini–Hochberg FDR are presented in Table 3. Six miRNAs (hsa-miR-6503-3p, -4433b-5p, 4433a-3p, -199b-5p, -21-3p, and -1261) were related to sex, four miRNAs (hsa-miR-664a-3p, -199b-5p, -6503-3p, and -21-3p) to medication with acetylsalicylic acid status, and two miRNAs (hsa-miR-664a-3p and -6503-3p) to the medication with diosmin (Table 3, Figure S9). All associated miRNAs were selected from the AAA vs. CVD comparison.
The results of the correlation analysis and Mann–Whitney U test showed that certain characteristics could confound differential expression of selected miRNAs. To explore these associations in more detail, multivariate linear regression models were constructed using the expression of particular miRNAs as a dependent variable and all 29 analyzed characteristics as independent variables. The entire results are provided in Table S9. The adjusted coefficients of determination obtained for the performed models were generally low and ranged from −0.059 to 0.351. Furthermore, a low number of variables (maximum 2) in particular models was found to contribute to the explanation of miRNA expression with a statistical significance of p < 0.01. These variables, together with associated miRNAs, are presented in Table 4. The characteristics that appeared the most frequently as related to the analyzed miRNAs were monocyte counts and levels of creatinine and urea (Table 4).
Obtained results indicate that the dysregulations of some miRNAs selected in this study could be affected by differences in demographic and clinical characteristics between the LEAD, AAA, and CVD groups. Although the found relationships seem rather weak, further investigations with larger and more balanced populations are needed to validate these findings.

4. Discussion

The elucidation of the molecular mechanisms regulating the pathological conditions responsible for the development of diseases could provide valuable information that could be useful in developing new methods for the management of these diseases. Many recent studies have shown that miRNAs are promising targets for the diagnosis and treatment of various conditions, including LEAD [49,50,51,52,53,54], AAA [28,55,56,57], and CVD [58,59]. In particular, many investigations were carried out to evaluate the utility of proangiogenic properties of miRNAs in the therapy of peripheral atherosclerosis [52]. For example, adipose tissue mesenchymal stem cells transfected with miRNA-126 exhibit a promising therapeutic potential in the diabetic mouse model with critical limb ischemia [49]. The proangiogenic properties of miR-548j-5p have been demonstrated to exert a beneficial effect on hindlimb ischemia in mouse models [50]. Many miRNAs have also been identified to be involved in various types of cardiomyocyte death, making them interesting therapeutic targets in ischemic events such as acute myocardial infarction and heart failure [60].
Our study contributes to the efforts by investigating differences and similarities in miRNA regulatory networks in PBMCs between these diseases. As a result, 10, 8, and 17 differentially expressed miRNA transcripts were selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively (Figure 1, Table 2).
In our study, a higher expression of hsa-miR-135b-5p was found in the LEAD group compared to patients with AAA and CVD using next-generation sequencing. This result is consistent with a previous study, which demonstrated that this miRNA is upregulated in serum samples from patients with coronary artery disease compared to healthy donors. Further analyses showed that hsa-miR-135b-5p promotes endothelial and vascular smooth muscle cell proliferation by suppressing MEF2C expression [61]. Our results suggest that this mechanism could be especially pronounced in LEAD in contrast to AAA and CVD; however, further studies are required to validate this conclusion.
Six miRNAs were found in our study to be dysregulated in arterial diseases (LEAD and AAA) in comparison with CVD (Figure 1). Among them, hsa-miR-99a-5p and hsa-miR-99a-3p were downregulated, and lower expression of both miRNAs was previously associated with the regulation of cell proliferation. hsa-miR-99a-5p was shown to target HOXA1, which promotes the proliferation, migration, and invasion of human aortic smooth muscle cells [62]. In turn, hsa-miR-99a-3p was downregulated in head and neck squamous cell carcinoma, and VEGFA was reported to be the main target gene for hsa-miR-99a-3p [63]. Therefore, our results showed that the enhancement of cell proliferation resulting from a lower expression of these miRNAs seems to be a common mechanism involved in both LEAD and AAA. This hypothesis could be supported by the upregulation of another miRNA in LEAD and AAA, hsa-miR-196a-5p, which exerts pro-oncogenic properties in various types of cancers [64,65,66].
Furthermore, a higher expression of this miRNA presented in our study was associated with its target, TSC2 (TSC complex subunit 2), which had a lower expression level in the LEAD group compared to CVD patients (Figure 3). This could indicate a higher proliferation status of LEAD since TSC2 is a known tumor suppressor acting mainly by inhibiting mTOR signaling [67,68].
Upregulation of hsa-miR-196a-5p was also linked to its other target, MALT1 (MALT1 paracaspase), which was downregulated in AAA vs. CVD group (Figure 3). Malt1 belongs to CARMA/Bcl10/MALT1 signalosomes activating NF-κB signaling in endothelial and vascular smooth muscle cells and thus promoting vascular inflammation [69,70]. Therefore, this pro-inflammatory mechanism appeared to be more active in CVD than LEAD, and more research is needed to confirm this conclusion.
Our study showed that the increase in expression of hsa-miR-212-3p is accompanied by an elevated level of its potential target, EIF3C (eukaryotic translation initiation factor 3 subunit C), in LEAD patients compared to CVD subjects (Figure 3). Increased expression of EIF3C, both in cells and exosomes, has previously been shown to promote the development and progression of hepatocellular carcinoma by improving cell proliferation and angiogenesis [71,72]. These processes are also involved in the pathology of atherosclerosis and could be promoted during this disease by upregulation of EIF3C.
In our study, a higher level of hsa-miR-124-3p was found in the LEAD group compared to AAA patients. This finding is in line with a previous study, where the expression of this miRNA was increased either in the ischemic tissue of the hindlimb ischemia model or in hypoxic human umbilical vein endothelial cells, where it exerted an anti-angiogenesis effect by targeting STAT3 [73]. Furthermore, our results showed that the target for this miRNA, GIT2, had a higher expression in the AAA group compared to LEAD patients (Figure 3). GIT2 is a key regulator of aging [74]; therefore, it could be assumed that AAA development is more closely associated with GIT2-related aging than LEAD.
Another gene found in our study as a potential target for has-miR-124-3p is YBX1 (Figure 3). Similarly to GIT2, the expression of YBX1 was found to be higher in AAA patients compared to the LEAD group. This gene has previously been reported to play an important role in aneurysm development through inhibition of proliferation and aggravating apoptosis of smooth muscle cells via upregulation of p21 [75].
These findings suggest that altered expression of hsa-miR-124-3p could change the character of vascular cells between the pro-apoptotic (in AAA) and pro-proliferative state (in LEAD). This hypothesis requires validation in further studies.
In our study, dysregulations in miRNA regulatory networks associated with LEAD, AAA, and CVD were identified to at least partially explain similarities and differences in the molecular mechanisms underlying the pathology of these diseases. However, we would like to emphasize that this work has a descriptive character, and all conclusions presented in Section 4 are hypotheses that should be validated in further studies. Further investigations are needed to evaluate the distribution of miRNA expression in subpopulations of PBMCs and to validate the presented miRNA–gene interactions. The translation potential of the results obtained in PBMCs to other blood components, including serum or plasma, should also be evaluated. It should also be emphasized that the presented study concerns diseases in a specific stage of development. Patients with LEAD suffered intermittent claudication, but without events of critical limb ischemia in the medical history. In the AAA group, patients with aneurysms sized 5.6 to 7.8 cm were included, which are typically classified as large aneurysms, in contrast to small aneurysms characterized by a diameter below 5.5 cm. Regarding CVD, patients with superficial varicose veins with primary etiology and reflux pathophysiology were qualified, without more advanced symptoms, such as edema, skin changes, and ulcers. Furthermore, due to the inability to use qPCR validation and an independent cohort to validate differentially expressed miRNAs in this study, the extended methodology was applied, including DESeq2 and ROC analyses, with multiple thresholds for miRNA selection (p-value < 0.05, absolute log2 fold change > 2, and mean of normalized counts > 10). Although the unification of the cut-off criteria for the selection of differentially expressed miRNAs from performed comparisons guaranteed the comparability of the obtained results between the studied groups, it could lead to the loss of some relevant miRNAs. Therefore, further studies focused on genes and miRNAs regulating particular processes closely related to the studied diseases (e.g., angiogenesis or inflammation) should be performed to reveal possible all miRNAs and genes involved in the development of LEAD, AAA, and CVD. Potential bias of our results could also be introduced by demographic and clinical differences between the studied groups. In particular, CVD patients are significantly younger than subjects with other diseases, and hypertensive patients are more represented in the LEAD group than in other groups. There are also clear differences in the medications taken by the studied patients (Table 1). To address this bias, we performed multivariate DESeq2 analysis, and the relationships between such characteristics and dysregulated miRNAs selected in this study were investigated. Certain miRNAs were found to be associated with age, sex, monocyte levels, medication, and other features (Table 3 and Table 4). Therefore, the change in expression of these miRNAs could be affected by these characteristics, and more studies should be performed in larger and more balanced cohorts to validate our results.
Our study demonstrated that investigations of miRNA and whole transcriptome expression profiles provide abundant information about molecular mechanisms that contribute to vascular diseases. The integrated analysis of miRNA and gene expression allowed us to identify dysregulations in the miRNA regulatory network in LEAD, AAA, and CVD. However, more explanatory studies are required to verify the hypotheses drawn from the obtained results and to evaluate the potential diagnostic and therapeutic utility of dysregulated miRNAs and genes presented in this study.

5. Conclusions

The presented work demonstrated similarities and differences in the dysregulations of the miRNA regulatory network in LEAD, AAA, and CVD. The identified miRNAs and their potential targets revealed abundant information about vascular pathology, providing many hypotheses on how dysregulated miRNA function in that PBMCs can contribute to the development of LEAD, AAA, and CVD. Therefore, this study indicates new research paths for further validation studies and provides new candidates for transcriptomic biomarkers of LEAD, AAA, and CVD with potential utility in differentiating these diseases.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm11123470/s1, Figure S1: Proportions of white blood cell subpopulations in the groups of patients with LEAD (lower extremity artery disease), AAA (abdominal aortic aneurysm), and CVD (chronic venous disease), Figure S2: Boxplot presenting log10 of Cook’s distances of miRNAs across samples, Figure S3: Quality control plots generated for results of differential miRNA expression analysis performed by DESeq2 package between 40 LEAD and 34 CVD patients, Figure S4: Quality control plots generated for results of differential miRNA expression analysis performed by DESeq2 package between 40 LEAD and 28 AAA patients, Figure S5: Quality control plots generated for results of differential miRNA expression analysis performed by DESeq2 package between 28 AAA and 34 CVD patients, Figure S6: Receiver operating characteristic (ROC) analysis plots performed for 10 miRNA transcripts selected from LEAD vs. AAA comparison, Figure S7: Receiver operating characteristic (ROC) analysis plots performed for 8 miRNA transcripts selected from LEAD vs. CVD comparison, Figure S8: Receiver operating characteristic (ROC) analysis plots performed for 17 miRNA transcripts selected from AAA vs. CVD comparison, Figure S9: Boxplots presenting expression of miRNAs in relation to their statistically significantly associated categorical characteristics of study population, as resulted from two-sided Mann–Whitney U test, Table S1: Differentially expressed miRNAs resulting from LEAD vs. AAA comparison using multivariate DESeq2 analysis, Table S2: Differentially expressed miRNAs resulting from LEAD vs. CVD comparison using multivariate DESeq2 analysis, Table S3: Differentially expressed miRNAs resulting from AAA vs. CVD comparison using multivariate DESeq2 analysis, Table S4: Results of receiver operating characteristic (ROC) analysis performed for 10, 8, and 17 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively, Table S5: Results of functional analysis performed using miRNet 2.0 website tool (https://www.mirnet.ca) on 7, 8, and 16 miRNAs selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively, Table S6: Interactions identified between miRNAs selected in this study from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons and genes selected in our previous study from corresponding comparisons, Table S7: Results of correlation analysis between characteristics of studied groups and expression of 10, 8, and 17 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively, Table S8: Results of two-sided Mann–Whitney U test performed to investigate differences between categorical characteristics of study subjects and expression of 10, 8, and 17 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively, Table S9: Results of multivariate regression analysis between characteristics of studied groups (independent variables) and expression of 10, 8, and 17 miRNA transcripts (dependent variables) selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively.

Author Contributions

Conceptualization, M.F., J.K., and A.B.-K.; data curation, D.P.Z. and A.S.; formal analysis, D.P.Z., K.P.R., and A.B.-K.; funding acquisition, A.B.-K.; investigation, D.P.Z. and K.P.R.; methodology, D.P.Z., K.P.R., M.F., J.K., and A.B.-K.; project administration, A.S. and A.B.-K.; resources, A.S., M.F., J.K., and A.B.-K.; software, D.P.Z. and A.S.; supervision, A.S., M.F., J.K., and A.B.-K.; validation, D.P.Z., K.P.R., M.F., and A.B.-K.; visualization, D.P.Z.; writing—original draft, D.P.Z.; writing—review and editing, K.P.R., D.G., M.F., and A.B.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by statutory funds of the Medical University of Lublin, number DS43 to A.B.-K.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of the Medical University of Lublin (protocol code KE-0254/341/2015, 17 December 2015).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data generated for this study are openly available in the FigShare repository at https://doi.org/10.6084/m9.figshare.19446851.v2.

Acknowledgments

The research was performed using the equipment purchased within the project “The equipment of innovative laboratories doing research on new medicines used in the therapy of civilization and neoplastic diseases” within the Operational Program Development of Eastern Poland 2007–2013, Priority Axis I Modern Economy, Operations I.3 Innovation Promotion, and presented results were partially obtained using the apparatus located in the Center ECOTECH-COMPLEX in Lublin, purchased under the project “ECOTECH-COMPLEX Man, Environment, Production” project No. POIG.02.01.00-06-212/09 funded from the European Regional Development Fund 2007-2013 by the National Center for Research and Development in Poland.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, D.; Sung, Y.M.; Park, J.; Kim, S.; Kim, J.; Park, J.; Ha, H.; Bae, J.Y.; Kim, S.; Baek, D. General Rules for Functional MicroRNA Targeting. Nat. Genet. 2016, 48, 1517–1526. [Google Scholar] [CrossRef] [PubMed]
  2. Saliminejad, K.; Khorshid, H.R.K.; Fard, S.S.; Ghaffari, S.H. An Overview of MicroRNAs: Biology, Functions, Therapeutics, and Analysis Methods. J. Cell. Physiol. 2019, 234, 5451–5465. [Google Scholar] [CrossRef] [PubMed]
  3. Rupaimoole, R.; Slack, F.J. MicroRNA Therapeutics: Towards a New Era for the Management of Cancer and Other Diseases. Nat. Rev. Drug Discov. 2017, 16, 203–222. [Google Scholar] [CrossRef] [PubMed]
  4. Kalayinia, S.; Arjmand, F.; Maleki, M.; Malakootian, M.; Singh, C.P. MicroRNAs: Roles in Cardiovascular Development and Disease. Cardiovasc. Pathol. 2021, 50, 107296. [Google Scholar] [CrossRef] [PubMed]
  5. Munjas, J.; Sopić, M.; Stefanović, A.; Košir, R.; Ninić, A.; Joksić, I.; Antonić, T.; Spasojević-Kalimanovska, V.; Prosenc Zmrzljak, U. Non-Coding RNAs in Preeclampsia—Molecular Mechanisms and Diagnostic Potential. Int. J. Mol. Sci. 2021, 22, 10652. [Google Scholar] [CrossRef] [PubMed]
  6. Aboyans, V.; Ricco, J.-B.; Bartelink, M.-L.E.L.; Björck, M.; Brodmann, M.; Cohnert, T.; Collet, J.-P.; Czerny, M.; De Carlo, M.; Debus, S.; et al. 2017 ESC Guidelines on the Diagnosis and Treatment of Peripheral Arterial Diseases, in Collaboration with the European Society for Vascular Surgery (ESVS). Eur. Heart J. 2018, 39, 763–816. [Google Scholar] [CrossRef] [Green Version]
  7. Hamburg, N.M.; Creager, M.A. Pathophysiology of Intermittent Claudication in Peripheral Artery Disease. Circ. J. 2017, 81, 281–289. [Google Scholar] [CrossRef] [Green Version]
  8. Song, P.; Rudan, D.; Zhu, Y.; Fowkes, F.J.I.; Rahimi, K.; Fowkes, F.G.R.; Rudan, I. Global, Regional, and National Prevalence and Risk Factors for Peripheral Artery Disease in 2015: An Updated Systematic Review and Analysis. Lancet Glob. Health 2019, 7, e1020–e1030. [Google Scholar] [CrossRef] [Green Version]
  9. Shu, J.; Santulli, G. Update on Peripheral Artery Disease: Epidemiology and Evidence-Based Facts. Atherosclerosis 2018, 275, 379–381. [Google Scholar] [CrossRef]
  10. Nehler, M.R.; Duval, S.; Diao, L.; Annex, B.H.; Hiatt, W.R.; Rogers, K.; Zakharyan, A.; Hirsch, A.T. Epidemiology of Peripheral Arterial Disease and Critical Limb Ischemia in an Insured National Population. J. Vasc. Surg. 2014, 60, 686–695.e2. [Google Scholar] [CrossRef] [Green Version]
  11. Chaikof, E.L.; Dalman, R.L.; Eskandari, M.K.; Jackson, B.M.; Lee, W.A.; Mansour, M.A.; Mastracci, T.M.; Mell, M.; Murad, M.H.; Nguyen, L.L.; et al. The Society for Vascular Surgery Practice Guidelines on the Care of Patients with an Abdominal Aortic Aneurysm. J. Vasc. Surg. 2018, 67, 2–77.e2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lattanzi, S. Abdominal Aortic Aneurysms: Pathophysiology and Clinical Issues. J. Intern. Med. 2020, 288, 376–378. [Google Scholar] [CrossRef] [PubMed]
  13. Kent, K.C.; Zwolak, R.M.; Egorova, N.N.; Riles, T.S.; Manganaro, A.; Moskowitz, A.J.; Gelijns, A.C.; Greco, G. Analysis of Risk Factors for Abdominal Aortic Aneurysm in a Cohort of More than 3 Million Individuals. J. Vasc. Surg. 2010, 52, 539–548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ahmed, R.; Ghoorah, K.; Kunadian, V. Abdominal Aortic Aneurysms and Risk Factors for Adverse Events. Cardiol. Rev. 2016, 24, 88–93. [Google Scholar] [CrossRef]
  15. Murali Krishna, S.; Morton, S.K.; Li, J.; Golledge, J. Risk Factors and Mouse Models of Abdominal Aortic Aneurysm Rupture. Int. J. Mol. Sci. 2020, 21, 7250. [Google Scholar] [CrossRef]
  16. Eberhardt, R.T.; Raffetto, J.D. Chronic Venous Insufficiency. Circulation 2014, 130, 333–346. [Google Scholar] [CrossRef]
  17. Ligi, D.; Croce, L.; Mannello, F. Chronic Venous Disorders: The Dangerous, the Good, and the Diverse. Int. J. Mol. Sci. 2018, 19, 2544. [Google Scholar] [CrossRef] [Green Version]
  18. Nicolaides, A.; Kakkos, S.; Baekgaard, N.; Comerota, A.; de Maeseneer, M.; Eklof, B.; Giannoukas, A.D.; Lugli, M.; Maleti, O.; Myers, K.; et al. Management of Chronic Venous Disorders of the Lower Limbs. Guidelines According to Scientific Evidence. Part I. Int. Angiol. 2018, 37, 181–254. [Google Scholar] [CrossRef]
  19. Nicolaides, A.; Kakkos, S.; Baekgaard, N.; Comerota, A.; de Maeseneer, M.; Eklof, B.; Giannoukas, A.; Lugli, M.; Maleti, O.; Mansilha, A.; et al. Management of Chronic Venous Disorders of the Lower Limbs. Guidelines According to Scientific Evidence. Part II. Int. Angiol. 2020, 39, 175–240. [Google Scholar] [CrossRef]
  20. Mansilha, A.; Sousa, J. Pathophysiological Mechanisms of Chronic Venous Disease and Implications for Venoactive Drug Therapy. Int. J. Mol. Sci. 2018, 19, 1669. [Google Scholar] [CrossRef] [Green Version]
  21. Lurie, F.; Passman, M.; Meisner, M.; Dalsing, M.; Masuda, E.; Welch, H.; Bush, R.L.; Blebea, J.; Carpentier, P.H.; De Maeseneer, M.; et al. The 2020 Update of the CEAP Classification System and Reporting Standards. J. Vasc. Surg. Venous Lymphat. Disord. 2020, 8, 342–352. [Google Scholar] [CrossRef] [PubMed]
  22. Vuylsteke, M.E.; Colman, R.; Thomis, S.; Guillaume, G.; Van Quickenborne, D.; Staelens, I. An Epidemiological Survey of Venous Disease Among General Practitioner Attendees in Different Geographical Regions on the Globe: The Final Results of the Vein Consult Program. Angiology 2018, 69, 779–785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Conte, M.S.; Bradbury, A.W.; Kolh, P.; White, J.V.; Dick, F.; Fitridge, R.; Mills, J.L.; Ricco, J.-B.; Suresh, K.R.; Murad, M.H.; et al. Global Vascular Guidelines on the Management of Chronic Limb-Threatening Ischemia. Eur. J. Vasc. Endovasc. Surg. 2019, 58, S1–S109.e33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wittens, C.; Davies, A.H.; Bækgaard, N.; Broholm, R.; Cavezzi, A.; Chastanet, S.; de Wolf, M.; Eggen, C.; Giannoukas, A.; Gohel, M.; et al. Editor’s Choice—Management of Chronic Venous Disease: Clinical Practice Guidelines of the European Society for Vascular Surgery (ESVS). Eur. J. Vasc. Endovasc. Surg. 2015, 49, 678–737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Qin, S.; Zhang, C. MicroRNAs in Vascular Disease. J. Cardiovasc. Pharmacol. 2011, 57, 8–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Fernández-Hernando, C.; Suárez, Y. MicroRNAs in Endothelial Cell Homeostasis and Vascular Disease. Curr. Opin. Hematol. 2018, 25, 227–236. [Google Scholar] [CrossRef]
  27. Ono, K.; Horie, T.; Baba, O.; Kimura, M.; Tsuji, S.; Rodriguez, R.R.; Miyagawa, S.; Kimura, T. Functional Non-Coding RNAs in Vascular Diseases. FEBS J. 2020, 288, 6315–6330. [Google Scholar] [CrossRef]
  28. Kumar, S.; Williams, D.; Sur, S.; Wang, J.-Y.; Jo, H. Role of Flow-Sensitive MicroRNAs and Long Noncoding RNAs in Vascular Dysfunction and Atherosclerosis. Vascul. Pharmacol. 2019, 114, 76–92. [Google Scholar] [CrossRef]
  29. Vogiatzi, G.; Oikonomou, E.; Deftereos, S.; Siasos, G.; Tousoulis, D. Peripheral Artery Disease: A Micro-RNA-Related Condition? Curr. Opin. Pharmacol. 2018, 39, 105–112. [Google Scholar] [CrossRef]
  30. Hamburg, N.M.; Leeper, N.J. Therapeutic Potential of Modulating MicroRNA in Peripheral Artery Disease. Curr. Vasc. Pharmacol. 2015, 13, 316–323. [Google Scholar] [CrossRef]
  31. Iyer, V.; Rowbotham, S.; Biros, E.; Bingley, J.; Golledge, J. A Systematic Review Investigating the Association of MicroRNAs with Human Abdominal Aortic Aneurysms. Atherosclerosis 2017, 261, 78–89. [Google Scholar] [CrossRef] [PubMed]
  32. Cui, C.; Liu, G.; Huang, Y.; Lu, X.; Lu, M.; Huang, X.; Li, W.; Jiang, M. MicroRNA Profiling in Great Saphenous Vein Tissues of Patients with Chronic Venous Insufficiency. Tohoku J. Exp. Med. 2012, 228, 341–350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Huang, X.; Liu, Z.; Shen, L.; Jin, Y.; Xu, G.; Zhang, Z.; Fang, C.; Guan, W.; Liu, C. Augmentation of MiR-202 in Varicose Veins Modulates Phenotypic Transition of Vascular Smooth Muscle Cells by Targeting Proliferator–Activated Receptor-γ Coactivator-1α. J. Cell. Biochem. 2019, 120, 10031–10042. [Google Scholar] [CrossRef] [PubMed]
  34. Anwar, M.A.; Adesina-Georgiadis, K.N.; Spagou, K.; Vorkas, P.A.; Li, J.V.; Shalhoub, J.; Holmes, E.; Davies, A.H. A Comprehensive Characterisation of the Metabolic Profile of Varicose Veins; Implications in Elaborating Plausible Cellular Pathways for Disease Pathogenesis. Sci. Rep. 2017, 7, 2989. [Google Scholar] [CrossRef]
  35. Saenz-Pipaon, G.; Martinez-Aguilar, E.; Orbe, J.; González Miqueo, A.; Fernandez-Alonso, L.; Paramo, J.A.; Roncal, C. The Role of Circulating Biomarkers in Peripheral Arterial Disease. Int. J. Mol. Sci. 2021, 22, 3601. [Google Scholar] [CrossRef]
  36. Barbalata, T.; Moraru, O.E.; Stancu, C.S.; Devaux, Y.; Simionescu, M.; Sima, A.V.; Niculescu, L.S. Increased MiR-142 Levels in Plasma and Atherosclerotic Plaques from Peripheral Artery Disease Patients with Post-Surgery Cardiovascular Events. Int. J. Mol. Sci. 2020, 21, 9600. [Google Scholar] [CrossRef]
  37. Plana, E.; Gálvez, L.; Medina, P.; Navarro, S.; Fornés-Ferrer, V.; Panadero, J.; Miralles, M. Identification of Novel MicroRNA Profiles Dysregulated in Plasma and Tissue of Abdominal Aortic Aneurysm Patients. Int. J. Mol. Sci. 2020, 21, 4600. [Google Scholar] [CrossRef]
  38. Torres-Do Rego, A.; Barrientos, M.; Ortega-Hernández, A.; Modrego, J.; Gómez-Gordo, R.; Álvarez-Sala, L.A.; Cachofeiro, V.; Gómez-Garre, D. Identification of a Plasma Microrna Signature as Biomarker of Subaneurysmal Aortic Dilation in Patients with High Cardiovascular Risk. J. Clin. Med. 2020, 9, 2783. [Google Scholar] [CrossRef]
  39. Nana, P.; Dakis, K.; Brodis, A.; Spanos, K.; Kouvelos, G. Circulating Biomarkers for the Prediction of Abdominal Aortic Aneurysm Growth. J. Clin. Med. 2021, 10, 1718. [Google Scholar] [CrossRef]
  40. Zalewski, D.P.; Ruszel, K.P.; Stępniewski, A.; Gałkowski, D.; Bogucki, J.; Komsta, Ł.; Kołodziej, P.; Chmiel, P.; Zubilewicz, T.; Feldo, M.; et al. Dysregulations of MicroRNA and Gene Expression in Chronic Venous Disease. J. Clin. Med. 2020, 9, 1251. [Google Scholar] [CrossRef]
  41. Bogucka-Kocka, A.; Zalewski, D.P.; Ruszel, K.P.; Stępniewski, A.; Gałkowski, D.; Bogucki, J.; Komsta, Ł.; Kołodziej, P.; Zubilewicz, T.; Feldo, M.; et al. Dysregulation of MicroRNA Regulatory Network in Lower Extremities Arterial Disease. Front. Genet. 2019, 10, 1200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Zalewski, D.P.; Ruszel, K.P.; Stępniewski, A.; Gałkowski, D.; Bogucki, J.; Komsta, Ł.; Kołodziej, P.; Chmiel, P.; Zubilewicz, T.; Feldo, M.; et al. Dysregulation of MicroRNA Modulatory Network in Abdominal Aortic Aneurysm. J. Clin. Med. 2020, 9, 1974. [Google Scholar] [CrossRef] [PubMed]
  43. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Chen, H.; Boutros, P.C. VennDiagram: A Package for the Generation of Highly-Customizable Venn and Euler Diagrams in R. BMC Bioinform. 2011, 12, 35. [Google Scholar] [CrossRef] [Green Version]
  45. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.-C.; Müller, M. PROC: An Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef]
  46. Chang, L.; Zhou, G.; Soufan, O.; Xia, J. MiRNet 2.0: Network-Based Visual Analytics for MiRNA Functional Analysis and Systems Biology. Nucleic Acids Res. 2020, 48, W244–W251. [Google Scholar] [CrossRef]
  47. Zalewski, D.P.; Ruszel, K.P.; Stępniewski, A.; Gałkowski, D.; Bogucki, J.; Kołodziej, P.; Szymańska, J.; Płachno, B.J.; Zubilewicz, T.; Feldo, M.; et al. Identification of Transcriptomic Differences between Lower Extremities Arterial Disease, Abdominal Aortic Aneurysm and Chronic Venous Disease in Peripheral Blood Mononuclear Cells Specimens. Int. J. Mol. Sci. 2021, 22, 3200. [Google Scholar] [CrossRef]
  48. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  49. Asadi-Yousefabad, S.-L.; Nammian, P.; Tabei, S.M.B.; Daneshi, S.; Fallahi, J.; Razban, V.; Sheikhha, M.H. Angiogenesis in Diabetic Mouse Model with Critical Limb Ischemia; Cell and Gene Therapy. Microvasc. Res. 2022, 141, 104339. [Google Scholar] [CrossRef]
  50. Lee, C.-Y.; Lin, S.-J.; Wu, T.-C. MiR-548j-5p Regulates Angiogenesis in Peripheral Artery Disease. Sci. Rep. 2022, 12, 838. [Google Scholar] [CrossRef]
  51. Lamin, V.; Verry, J.; Eigner-Bybee, I.; Fuqua, J.D.; Wong, T.; Lira, V.A.; Dokun, A.O. Modulation of MiR-29a and ADAM12 Reduces Post-Ischemic Skeletal Muscle Injury and Improves Perfusion Recovery and Skeletal Muscle Function in a Mouse Model of Type 2 Diabetes and Peripheral Artery Disease. Int. J. Mol. Sci. 2022, 23, 429. [Google Scholar] [CrossRef] [PubMed]
  52. Theofilis, P.; Oikonomou, E.; Vogiatzi, G.; Antonopoulos, A.S.; Siasos, G.; Iliopoulos, D.C.; Perrea, D.; Tsioufis, C.; Tousoulis, D. The Impact of Proangiogenic MicroRNA Modulation on Blood Flow Recovery Following Hind Limb Ischemia. A Systematic Review and Meta-Analysis of Animal Studies. Vascul. Pharmacol. 2021, 141, 106906. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, M.; Zhang, W.; Zhang, L.; Wang, L.; Li, J.; Shu, C.; Li, X. Roles of MicroRNAs in Peripheral Artery In-Stent Restenosis after Endovascular Treatment. BioMed Res. Int. 2021, 2021, e9935671. [Google Scholar] [CrossRef] [PubMed]
  54. Xu, H.; Ni, Y.-Q.; Liu, Y.-S. Mechanisms of Action of MiRNAs and LncRNAs in Extracellular Vesicle in Atherosclerosis. Front. Cardiovasc. Med. 2021, 8, 733985. [Google Scholar] [CrossRef] [PubMed]
  55. Yang, Z.; Zhang, L.; Liu, Y.; Zeng, W.; Wang, K. Potency of MiR-144-3p in Promoting Abdominal Aortic Aneurysm Progression in Mice Correlates with Apoptosis of Smooth Muscle Cells. Vascul. Pharmacol. 2022, 142, 106901. [Google Scholar] [CrossRef] [PubMed]
  56. Si, X.; Chen, Q.; Zhang, J.; Zhou, W.; Chen, L.; Chen, J.; Deng, N.; Li, W.; Liu, D.; Wang, L.; et al. MicroRNA-23b Prevents Aortic Aneurysm Formation by Inhibiting Smooth Muscle Cell Phenotypic Switching via FoxO4 Suppression. Life Sci. 2022, 288, 119092. [Google Scholar] [CrossRef]
  57. Huang, T.; Liu, S.; Liu, R.; Pan, B.; Wang, W. Inhibition of MiR-188-5p Suppresses Progression of Experimental Abdominal Aortic Aneurysms. J. Cardiovasc. Pharmacol. 2021, 77, 107–114. [Google Scholar] [CrossRef]
  58. Cao, Y.; Cao, Z.; Wang, W.; Jie, X.; Li, L.; Cao, Y.; Cao, Z.; Wang, W.; Jie, X.; Li, L. MicroRNA-199a-5p Regulates FOXC2 to Control Human Vascular Smooth Muscle Cell Phenotypic Switch. Mol. Med. Rep. 2021, 24, 627. [Google Scholar] [CrossRef]
  59. Wu, J.; Li, X.; Li, D.; Ren, X.; Li, Y.; Herter, E.K.; Qian, M.; Toma, M.-A.; Wintler, A.-M.; Sérézal, I.G.; et al. MicroRNA-34 Family Enhances Wound Inflammation by Targeting LGR4. J. Investig. Dermatol. 2020, 140, 465–476.e11. [Google Scholar] [CrossRef]
  60. Kansakar, U.; Varzideh, F.; Mone, P.; Jankauskas, S.S.; Santulli, G. Functional Role of MicroRNAs in Regulating Cardiomyocyte Death. Cells 2022, 11, 983. [Google Scholar] [CrossRef]
  61. Xu, Z.; Han, Y.; Liu, J.; Jiang, F.; Hu, H.; Wang, Y.; Liu, Q.; Gong, Y.; Li, X. MiR-135b-5p and MiR-499a-3p Promote Cell Proliferation and Migration in Atherosclerosis by Directly Targeting MEF2C. Sci. Rep. 2015, 5, 12276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Han, Z.; Guan, Y.; Liu, B.; Lin, Y.; Yan, Y.; Wang, H.; Wang, H.; Jing, B. MicroRNA-99a-5p Alleviates Atherosclerosis via Regulating Homeobox A1. Life Sci. 2019, 232, 116664. [Google Scholar] [CrossRef] [PubMed]
  63. Wei, G.-G.; Guo, W.-P.; Tang, Z.-Y.; Li, S.-H.; Wu, H.-Y.; Zhang, L.-C. Expression Level and Prospective Mechanism of MiRNA-99a-3p in Head and Neck Squamous Cell Carcinoma Based on MiRNA-Chip and MiRNA-Sequencing Data in 1, 167 Cases. Pathol. Res. Pract. 2019, 215, 963–976. [Google Scholar] [CrossRef] [PubMed]
  64. Pan, Y.; Shu, X.; Sun, L.; Yu, L.; Sun, L.; Yang, Z.; Ran, Y. MiR-196a-5p Modulates Gastric Cancer Stem Cell Characteristics by Targeting Smad4. Int. J. Oncol. 2017, 50, 1965–1976. [Google Scholar] [CrossRef] [Green Version]
  65. Xin, H.; Wang, C.; Liu, Z. MiR-196a-5p Promotes Metastasis of Colorectal Cancer via Targeting IκBα. BMC Cancer 2019, 19, 30. [Google Scholar] [CrossRef]
  66. Li, X.; Chen, Z.; Ni, Y.; Bian, C.; Huang, J.; Chen, L.; Xie, X.; Wang, J. Tumor-Associated Macrophages Secret Exosomal MiR-155 and MiR-196a-5p to Promote Metastasis of Non-Small-Cell Lung Cancer. Transl Lung Cancer Res. 2021, 10, 1338–1354. [Google Scholar] [CrossRef]
  67. Huang, J.; Manning, B.D. A Complex Interplay between Akt, TSC2 and the Two MTOR Complexes. Biochem. Soc. Trans. 2009, 37, 217–222. [Google Scholar] [CrossRef] [Green Version]
  68. Motoshima, H.; Goldstein, B.J.; Igata, M.; Araki, E. AMPK and Cell Proliferation—AMPK as a Therapeutic Target for Atherosclerosis and Cancer. J. Physiol. 2006, 574, 63–71. [Google Scholar] [CrossRef]
  69. McAllister-Lucas, L.M.; Jin, X.; Gu, S.; Siu, K.; McDonnell, S.; Ruland, J.; Delekta, P.C.; Van Beek, M.; Lucas, P.C. The CARMA3-Bcl10-MALT1 Signalosome Promotes Angiotensin II-Dependent Vascular Inflammation and Atherogenesis. J. Biol. Chem. 2010, 285, 25880–25884. [Google Scholar] [CrossRef] [Green Version]
  70. Birch, C.A.; Molinar-Inglis, O.; Trejo, J. Subcellular Hot Spots of GPCR Signaling Promote Vascular Inflammation. Curr. Opin. Endocr. Metab. Res. 2021, 16, 37–42. [Google Scholar] [CrossRef]
  71. Lee, H.-Y.; Chen, C.-K.; Ho, C.-M.; Lee, S.-S.; Chang, C.-Y.; Chen, K.-J.; Jou, Y.-S. EIF3C-Enhanced Exosome Secretion Promotes Angiogenesis and Tumorigenesis of Human Hepatocellular Carcinoma. Oncotarget 2018, 9, 13193–13205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Li, T.; Li, S.; Chen, D.; Chen, B.; Yu, T.; Zhao, F.; Wang, Q.; Yao, M.; Huang, S.; Chen, Z.; et al. Transcriptomic Analyses of RNA-Binding Proteins Reveal EIF3c Promotes Cell Proliferation in Hepatocellular Carcinoma. Cancer Sci. 2017, 108, 877–885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Shi, Y.; Xu, X.; Luan, P.; Kou, W.; Li, M.; Yu, Q.; Zhuang, J.; Xu, Y.; Peng, W.; Jian, W. MiR-124-3p Regulates Angiogenesis in Peripheral Arterial Disease by Targeting STAT3. Mol. Med. Rep. 2020, 22, 4890–4898. [Google Scholar] [CrossRef] [PubMed]
  74. Van Gastel, J.; Boddaert, J.; Jushaj, A.; Premont, R.T.; Luttrell, L.M.; Janssens, J.; Martin, B.; Maudsley, S. GIT2—A Keystone in Ageing and Age-Related Disease. Ageing Res. Rev. 2018, 43, 46–63. [Google Scholar] [CrossRef]
  75. He, X.; Wang, S.; Li, M.; Zhong, L.; Zheng, H.; Sun, Y.; Lai, Y.; Chen, X.; Wei, G.; Si, X.; et al. Long Noncoding RNA GAS5 Induces Abdominal Aortic Aneurysm Formation by Promoting Smooth Muscle Apoptosis. Theranostics 2019, 9, 5558–5576. [Google Scholar] [CrossRef]
Jcm 11 03470 g001 550
Figure 1. The comparison of the sets of differentially expressed miRNA transcripts selected from pairwise comparative analysis performed among the LEAD, AAA, and CVD groups using DESeq2 method. (A) The comparison of the sets of miRNA transcripts selected from LEAD vs. CVD, LEAD vs. AAA, and AAA vs. CVD comparisons using the following selection criteria: p-value < 0.05, absolute log2 fold change > 1, and mean of normalized counts > 10. The number in the intersection fields of the Venn diagram represents the number of miRNA transcripts common to the comparisons. The frames below the Venn diagram include six miRNAs common for AAA vs. CVD and LEAD vs. CVD comparisons and one miRNA common for LEAD vs. CVD and LEAD vs. AAA comparisons. (B) Heatmap of the expression of 28 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons and included in the Venn diagram in panel A. Hierarchical clustering was performed using the average method applied to Euclidean distances. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Figure 1. The comparison of the sets of differentially expressed miRNA transcripts selected from pairwise comparative analysis performed among the LEAD, AAA, and CVD groups using DESeq2 method. (A) The comparison of the sets of miRNA transcripts selected from LEAD vs. CVD, LEAD vs. AAA, and AAA vs. CVD comparisons using the following selection criteria: p-value < 0.05, absolute log2 fold change > 1, and mean of normalized counts > 10. The number in the intersection fields of the Venn diagram represents the number of miRNA transcripts common to the comparisons. The frames below the Venn diagram include six miRNAs common for AAA vs. CVD and LEAD vs. CVD comparisons and one miRNA common for LEAD vs. CVD and LEAD vs. AAA comparisons. (B) Heatmap of the expression of 28 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons and included in the Venn diagram in panel A. Hierarchical clustering was performed using the average method applied to Euclidean distances. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Jcm 11 03470 g001
Jcm 11 03470 g002 550
Figure 2. The functional network of miRNAs selected from pairwise comparisons among the LEAD, AAA, and CVD groups. (A) Network for miRNAs selected from the LEAD vs. AAA comparison. (B) Network for miRNAs selected from the LEAD vs. CVD comparison. (C) Network for miRNAs selected from the AAA vs. CVD comparison. Each panel presents the top 15 most enriched functional terms (with the lowest p-value of enrichment) of ‘miRNA Function’ category in the miRNet 2.0 tool as well as associated miRNAs. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Figure 2. The functional network of miRNAs selected from pairwise comparisons among the LEAD, AAA, and CVD groups. (A) Network for miRNAs selected from the LEAD vs. AAA comparison. (B) Network for miRNAs selected from the LEAD vs. CVD comparison. (C) Network for miRNAs selected from the AAA vs. CVD comparison. Each panel presents the top 15 most enriched functional terms (with the lowest p-value of enrichment) of ‘miRNA Function’ category in the miRNet 2.0 tool as well as associated miRNAs. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Jcm 11 03470 g002
Jcm 11 03470 g003 550
Figure 3. The regulatory network containing miRNAs selected in the current study and genes selected from a previous study [47] regarding (A) LEAD vs. AAA comparison, (B) AAA vs. CVD comparison, and (C) LEAD vs. CVD comparison. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Figure 3. The regulatory network containing miRNAs selected in the current study and genes selected from a previous study [47] regarding (A) LEAD vs. AAA comparison, (B) AAA vs. CVD comparison, and (C) LEAD vs. CVD comparison. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease.
Jcm 11 03470 g003
Table
Table 1. Demographical and clinical information regarding the study subjects. The table presents extended data gathered from our previous studies [40,41,42].
Table 1. Demographical and clinical information regarding the study subjects. The table presents extended data gathered from our previous studies [40,41,42].
CharacteristicPatients with LEAD
(n = 40)		Patients with AAA
(n = 28)		Patients with CVD (n = 34)p1
Age57.6 ± 9.82
56.5 (43–71)		66.4 ± 4.52
67 (57–76)		44.1 ± 10.07
42 (27–78)		9.966 × 10−12
Sex males/females35 (87.5%)/5 (12.5%)25 (89.3%)/
3 (10.7%)		17 (50%)/
17 (50%)		2.180 × 10−4
Body mass index (BMI)27.2 ± 2.62
27.3 (21.9–31.6)		25.1 ± 3.30
25.2 (18.0–31.3)		23.9 ± 2.35
24.0 (20.1–28.8)		1.670 × 10−5
Risk factors and cardiovascular comorbidities
Smoking never/former/current0 (0%)/18 (45%)/
22 (55%)		9 (32.1%)/
9 (32.1%)/10 (35.7%)		16 (47%)/13 (38%)/
5 (14.7%)		2.369 × 10−5
Diabetes type 25 (12.5%)6 (21.4%)0 (0%)9.309 × 10−3
Hypertension36 (90%)19 (67.9%)0 (0%)1.066 × 10−17
Coronary artery disease (CAD)11 (27.5%)7 (25%)0 (0%)8.508 × 10−4
Myocardial infarction8 (20%)3 (10.7%)0 (0%)0.011
Stroke/transient ischemic attack2 (5%)1 (3.6%)0 (0%)0.493
Hematological and biochemical blood parameters
Red blood cells (M/µL)4.74 ± 0.30
4.75 (4.11–5.18)		4.94 ± 0.21
4.99 (4.56–5.50)		4.78 ± 0.34
4.83 (3.78–5.50)		0.017
White blood cells (K/µL)5.49 ± 0.69
5.27 (4.45–6.89)		5.66 ± 0.70
5.78 (4.44–6.90)		5.76 ± 0.71
5.78 (4.67–6.90)		0.361
Lymphocytes (K/µL)3.04 ± 0.54
3.03 (2.01–3.99)		2.99 ± 0.44
3.02 (2.05–3.78)		2.99 ± 0.45
3.00 (2.04–3.99)		0.771
Monocytes (K/µL)0.47 ± 0.15
0.41 (0.22–0.87)		0.39 ± 0.12
0.33 (0.26–0.67)		0.33 ± 0.07
0.32 (0.19–0.56)		9.727 × 10−6
Neutrophils (K/µL)4.21 ± 0.47
4.07 (3.51–5.21)		4.17 ± 0.43
4.20 (3.33–4.78)		4.17 ± 0.62
4.31 (2.12–4.99)		0.809
Eosinophils (K/µL)0.21 ± 0.09
0.20 (0.10–0.56)		0.18 ± 0.05
0.18 (0.11–0.31)		0.16 ± 0.04
0.17 (0.08–0.23)		6.072 × 10−3
Basophils (K/µL)0.10 ± 0.03
0.1 (0.07–0.19)		0.11 ± 0.02
0.11 (0.08–0.20)		0.10 ± 0.03
0.10 (0.02–0.20)		3.955 × 10−3
Platelets (K/µL)309.3 ± 75.7
295.5 (179–561)		419.9 ± 124.0
404.5 (211–756)		368.26 ± 104.15
365.5 (211–756)		4.643 × 10−5
Hemoglobin (g/dL)14.12 ± 0.52
14.01 (12.99–14.99)		14.02 ± 0.51
13.99 (13.34–15.00)		14.08 ± 0.48
14.06 (13.30–14.99)		0.692
Hematocrit (%)41.3 ± 1.42
41.5 (38.4–43.8)		40.8 ± 1.30
41 (38–43)		39.9 ± 1.63
40 (37–43)		6.775 × 10−4
Creatinine (mmol/L)78.7 ± 12.6
79 (56–99)		54.2 ± 11.5
51.5 (39–87)		58.2 ± 14.20
56 (37–93)		2.086 × 10−10
Urea (mmol/L)4.69 ± 0.83
4.79 (2.99–6.02)		4.66 ± 0.67
4.69 (3.45–5.88)		5.09 ± 0.86
4.94 (3.45–6.87)		0.119
Medication
Statins34 (85%)13 (46.4%)0 (0%)2.090 × 10−13
Acetylsalicylic acid40 (100%)27 (96.4%)0 (0%)7.961 × 10−26
Clopidogrel8 (20%)3 (10.7%)0 (0%)6.382 × 10−3
Beta-adrenergic blockers27 (67.5%)16 (57.1%)0 (0%)4.427 × 10−12
Angiotensin-converting enzyme inhibitor20 (50%)4 (14.3%)0 (0%)6.266 × 10−7
Ca2+ channel blockers11 (27.5%)2 (7.14%)0 (0%)5.431 × 10−4
Fibrates5 (12.5%)2 (7.14%)0 (0%)0.092
Metformin2 (5%)3 (10.7%)0 (0%)0.148
Gliclazide4 (10%)4 (14.3%)0 (0%)0.055
Micronized diosmin0 (0%)0 (0%)24 (70.6%)1.012 × 10−15
Preparation with hesperidin, Ruscus aculeatus, and vitamin C0 (0%)0 (0%)15 (44.1%)4.550 × 10−8
1 Statistical significance of differences between lower extremity artery disease (LEAD), abdominal aortic aneurysm (AAA), and chronic venous disease (CVD) groups, calculated using Kruskal–Wallis rank sum test for continuous-type variables (age, BMI, hematological and biochemical blood parameters) and two-sided Fisher’s exact test for categorical-type variables (sex, risk factors, cardiovascular comorbidities, and medication). Data for continuous-type variables were expressed as mean ± SD and median (range). Data for categorical-type variables were expressed as numbers (percentage).
Table
Table 2. Differential expression parameters of 10, 8, and 17 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively.
Table 2. Differential expression parameters of 10, 8, and 17 miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons, respectively.
No.miRNA TranscriptmiRNA ID 1pFold ChangeROC-AUC
LEAD vs. AAA—Upregulated miRNAs
1.hsa-mir-124-2_hsa-miR-124-3phsa-miR-124-3p6.715 × 10−47.4410.585
2.hsa-mir-124-1_hsa-miR-124-3phsa-miR-124-3p1.665 × 10−36.0150.550
3.hsa-mir-124-3_hsa-miR-124-3phsa-miR-124-3p1.946 × 10−36.4550.562
4.hsa-mir-7641-2_hsa-miR-7641hsa-miR-76417.019 × 10−33.3980.646
5.hsa-mir-205_hsa-miR-205-5phsa-miR-205-5p9.943 × 10−33.4140.558
6.hsa-mir-31_hsa-miR-31-3phsa-miR-31-3p1.066 × 10−22.0400.749
7.hsa-mir-135b_hsa-miR-135b-5phsa-miR-135b-5p1.121 × 10−22.5650.656
8.hsa-mir-7641-1_hsa-miR-7641hsa-miR-76412.886 × 10−22.2310.629
LEAD vs. AAA—Downregulated miRNAs
9.hsa-mir-4531_hsa-miR-4531hsa-miR-45318.039 × 10−40.4390.521
10.hsa-mir-1303_hsa-miR-1303hsa-miR-13039.764 × 10−40.3860.629
LEAD vs. CVD—Upregulated miRNAs
1.hsa-mir-196a-1_hsa-miR-196a-5phsa-miR-196a-5p7.766 × 10−314.1060.603
2.hsa-mir-98_hsa-miR-98-5phsa-miR-98-5p1.599 × 10−22.0860.646
3.hsa-mir-212_hsa-miR-212-3phsa-miR-212-3p3.163 × 10−22.3390.558
4.hsa-mir-135b_hsa-miR-135b-5phsa-miR-135b-5p3.597 × 10−27.9390.524
LEAD vs. CVD—Downregulated miRNAs
5.hsa-mir-99a_hsa-miR-99a-5phsa-miR-99a-5p5.074 × 10−30.2110.626
6.hsa-mir-99a_hsa-miR-99a-3phsa-miR-99a-3p5.534 × 10−30.1130.607
7.hsa-mir-125b-2_hsa-miR-125b-5phsa-miR-125b-5p1.317 × 10−20.2590.659
8.hsa-mir-193b_hsa-miR-193b-3phsa-miR-193b-3p4.348 × 10−20.2630.527
AAA vs. CVD—Upregulated miRNAs
1.hsa-mir-98_hsa-miR-98-5phsa-miR-98-5p5.135 × 10−32.2270.674
2.hsa-mir-196a-1_hsa-miR-196a-5phsa-miR-196a-5p8.853 × 10−311.4380.516
3.hsa-mir-199b_hsa-miR-199b-5phsa-miR-199b-5p1.133 × 10−22.6080.750
4.hsa-mir-7-1_hsa-miR-7-5phsa-miR-7-5p1.496 × 10−22.0790.601
5.hsa-mir-1261_hsa-miR-1261hsa-miR-12611.766 × 10−22.2560.611
6.hsa-mir-21_hsa-miR-21-3phsa-miR-21-3p3.115 × 10−22.5670.798
7.hsa-mir-99b_hsa-miR-99b-3phsa-miR-99b-3p3.642 × 10−23.6520.564
8.hsa-mir-6503_hsa-miR-6503-3phsa-miR-6503-3p4.350 × 10−22.3130.733
9.hsa-mir-550a-3_hsa-miR-550a-3phsa-miR-550a-3p4.373 × 10−22.0260.693
AAA vs. CVD—Downregulated miRNAs
10.hsa-mir-664a_hsa-miR-664a-3phsa-miR-664a-3p1.218 × 10−40.4990.824
11.hsa-mir-125b-2_hsa-miR-125b-5phsa-miR-125b-5p8.773 × 10−40.1820.619
12.hsa-mir-99a_hsa-miR-99a-5phsa-miR-99a-5p1.189 × 10−30.1850.608
13.hsa-mir-99a_hsa-miR-99a-3phsa-miR-99a-3p1.251 × 10−30.0920.640
14.hsa-mir-125b-1_hsa-miR-125b-5phsa-miR-125b-5p6.648 × 10−30.2470.598
15.hsa-mir-4433b_hsa-miR-4433b-5phsa-miR-4433b-5p2.297 × 10−20.1970.544
16.hsa-mir-4433_hsa-miR-4433-3phsa-miR-4433a-3p2.297 × 10−20.1970.544
17.hsa-mir-193b_hsa-miR-193b-3phsa-miR-193b-3p4.728 × 10−20.2920.450
1 According to miRBase 22.1 (http://www.mirbase.org). The table presents p and fold change values received from the DESeq2 analysis and areas under ROC curves (ROC-AUC) received from the ROC analysis. miRNA transcripts were divided into upregulated and downregulated groups within each comparison and ordered according to increasing p-value. AAA—abdominal aortic aneurysm, CVD—chronic venous disease, LEAD—lower extremity artery disease, ROC—receiver operating characteristic.
Table
Table 3. Relationships between categorical characteristics of study subjects and expression of miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons. The table presents statistically significant relationships (p < 0.05, Benjamini–Hochberg FDR correction) obtained from two-sided Mann–Whitney U test (the entire results are provided in Table S8).
Table 3. Relationships between categorical characteristics of study subjects and expression of miRNA transcripts selected from LEAD vs. AAA, LEAD vs. CVD, and AAA vs. CVD comparisons. The table presents statistically significant relationships (p < 0.05, Benjamini–Hochberg FDR correction) obtained from two-sided Mann–Whitney U test (the entire results are provided in Table S8).
CharacteristicmiRNA TranscriptmiRNA ID 1p
Sexhsa-mir-6503_hsa-miR-6503-3phsa-miR-6503-3p3.02 × 10−3
hsa-mir-4433b_hsa-miR-4433b-5phsa-miR-4433b-5p1.17 × 10−2
hsa-mir-4433_hsa-miR-4433-3phsa-miR-4433a-3p1.17 × 10−2
hsa-mir-199b_hsa-miR-199b-5phsa-miR-199b-5p1.17 × 10−2
hsa-mir-21_hsa-miR-21-3phsa-miR-21-3p1.48 × 10−2
hsa-mir-1261_hsa-miR-1261hsa-miR-12612.21 × 10−2
Acetylsalicylic acid medicationhsa-mir-664a_hsa-miR-664a-3phsa-miR-664a-3p3.02 × 10−3
hsa-mir-199b_hsa-miR-199b-5phsa-miR-199b-5p9.42 × 10−3
hsa-mir-6503_hsa-miR-6503-3phsa-miR-6503-3p1.17 × 10−2
hsa-mir-21_hsa-miR-21-3phsa-miR-21-3p1.48 × 10−2
Diosmin medicationhsa-mir-664a_hsa-miR-664a-3phsa-miR-664a-3p3.42 × 10−2
hsa-mir-6503_hsa-miR-6503-3phsa-miR-6503-3p4.04 × 10−2
1 According to miRBase 22.1 (http://www.mirbase.org).
Table
Table 4. Results of the multivariate linear regression performed using the expression of particular miRNA transcripts as a dependent variable and the characteristics of the studied subjects as independent variables. The table presents regression coefficients and associated p-values of characteristics resulting in p < 0.01. The entire results of the multivariate regression analysis are provided in Table S9.
Table 4. Results of the multivariate linear regression performed using the expression of particular miRNA transcripts as a dependent variable and the characteristics of the studied subjects as independent variables. The table presents regression coefficients and associated p-values of characteristics resulting in p < 0.01. The entire results of the multivariate regression analysis are provided in Table S9.
CharacteristicmiRNA TranscriptmiRNA ID 1Regression Coefficientp
Monocyte countshsa-mir-99b_hsa-miR-99b-3p 4hsa-miR-99b-3p−1.386.07 × 10−3
hsa-mir-4433b_hsa-miR-4433b-5p 4hsa-miR-4433b-5p−2.126.03 × 10−3
hsa-mir-4433_hsa-miR-4433-3p 4hsa-miR-4433a-3p−2.126.03 × 10−3
hsa-mir-31_hsa-miR-31-3p 2hsa-miR-31-3p2.232.01 × 10−3
hsa-mir-7-1_hsa-miR-7-5p 4hsa-miR-7-5p0.921.44 × 10−3
Creatininehsa-mir-125b-1_hsa-miR-125b-5p 4hsa-miR-125b-5p−0.012.02 × 10−3
hsa-mir-125b-2_hsa-miR-125b-5p 3,4hsa-miR-125b-5p −0.016.17 × 10−3
Ureahsa-mir-125b-1_hsa-miR-125b-5p 4hsa-miR-125b-5p0.205.24 × 10−3
hsa-mir-125b-2_hsa-miR-125b-5p 3,4hsa-miR-125b-5p0.214.55 × 10−3
Sexhsa-mir-7641-2_hsa-miR-7641 2hsa-miR-76410.724.40 × 10−3
Myocardial infarctionhsa-mir-7-1_hsa-miR-7-5p 4hsa-miR-7-5p−0.575.99 × 10−3
1 According to miRBase 22.1 (http://www.mirbase.org), 2 miRNA transcripts selected from LEAD vs. AAA comparison, 3 miRNA transcripts selected from LEAD vs. CVD comparison, 4 miRNA transcripts selected from AAA vs. CVD comparison, p—statistical significance, LEAD—lower extremity artery disease, AAA—abdominal aortic aneurysm, CVD—chronic venous disease.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zalewski, D.P.; Ruszel, K.P.; Stępniewski, A.; Gałkowski, D.; Feldo, M.; Kocki, J.; Bogucka-Kocka, A. miRNA Regulatory Networks Associated with Peripheral Vascular Diseases. J. Clin. Med. 2022, 11, 3470. https://doi.org/10.3390/jcm11123470

AMA Style

Zalewski DP, Ruszel KP, Stępniewski A, Gałkowski D, Feldo M, Kocki J, Bogucka-Kocka A. miRNA Regulatory Networks Associated with Peripheral Vascular Diseases. Journal of Clinical Medicine. 2022; 11(12):3470. https://doi.org/10.3390/jcm11123470

Chicago/Turabian Style

Zalewski, Daniel P., Karol P. Ruszel, Andrzej Stępniewski, Dariusz Gałkowski, Marcin Feldo, Janusz Kocki, and Anna Bogucka-Kocka. 2022. "miRNA Regulatory Networks Associated with Peripheral Vascular Diseases" Journal of Clinical Medicine 11, no. 12: 3470. https://doi.org/10.3390/jcm11123470

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop