Circulating Plasma miRNA and Clinical/Hemodynamic Characteristics Provide Additional Predictive Information About Acute Pulmonary Thromboembolism, Chronic Thromboembolic Pulmonary Hypertension and Idiopathic Pulmonary Hypertension

Idiopathic pulmonary artery hypertension (IPAH), chronic thromboembolic pulmonary hypertension (CTEPH), and acute pulmonary embolism (APTE) are life-threatening cardiopulmonary diseases without specific surgical or medical treatment. Although APTE, CTEPH and IPAH are different pulmonary vascular diseases in terms of clinical presentation, prevalence, pathophysiology and prognosis, the identification of their circulating microRNA (miRNAs) might help in recognizing differences in their outcome evolution and clinical forms. The aim of this study was to describe the APTE, CTEPH, and IPAH-associated miRNAs and to predict their target genes. The target genes of the key differentially expressed miRNAs were analyzed, and functional enrichment analyses were carried out. The miRNAs were detected using RT-PCR. Finally, we incorporated plasma circulating miRNAs in baseline and clinical characteristics of the patients to detect differences between APTE and CTEPH in time of evolution, and differences between CTEPH and IPAH in diseases form. We found five top circulating plasma miRNAs in common with APTE, CTEPH and IPAH assembled in one conglomerate. Among them, miR-let-7i-5p expression was upregulated in APTE and IPAH, while miRNA-320a was upregulated in CTEP and IPAH. The network construction for target genes showed 11 genes regulated by let-7i-5p and 20 genes regulated by miR-320a, all of them regulators of pulmonary arterial adventitial fibroblasts, pulmonary artery endothelial cell, and pulmonary artery smooth muscle cells. AR (androgen receptor), a target gene of hsa-let-7i-5p and has-miR-320a, was enriched in pathways in cancer, whereas PRKCA (Protein Kinase C Alpha), also a target gene of hsa-let-7i-5p and has-miR-320a, was enriched in KEGG pathways, such as pathways in cancer, glioma, and PI3K-Akt signaling pathway. We inferred that CTEPH might be the consequence of abnormal remodeling in APTE, while unbalance between the hyperproliferative and apoptosis-resistant phenotype of pulmonary arterial adventitial fibroblasts, pulmonary artery endothelial cell and pulmonary artery smooth muscle cells in pulmonary artery confer differences in IPAH and CTEPH diseases form. We concluded that the incorporation of plasma circulating let-7i-5p and miRNA-320a in baseline and clinical characteristics of the patients reinforces differences between APTE and CTEPH in outcome evolution, as well as differences between CTEPH and IPAH in diseases form.


INTRODUCTION
Acute pulmonary thromboembolism (APTE), chronic thromboembolic pulmonary hypertension (CTEPH) and idiopathic pulmonary artery hypertension (IPAH) are pulmonary vascular diseases, that can reveal a pathophysiological continuum. On the other hand, they are rather different pulmonary vascular diseases in terms of clinical presentation, prevalence, pathophysiology, and prognosis.
Acute pulmonary thromboembolism is characterized by acute obstruction of the pulmonary artery usually leading to pulmonary artery hypertension (PH), being the third most common cause of acute cardiovascular death (Douma et al., 2010;Goldhaber et al., 2012;Konstantinides et al., 2014). However, due to the non-specific clinical symptoms and signs, the diagnosis of APTE is often missed or delayed (Kostadima and Zakynthinos, 2007). CTEPH, a complication of APTE and the fourth types of PH, is characterized by the remodeling of the pulmonary arteries by fibro-organized tissue Connell et al., 2015(Simonneau et al., 2018. The incidence of CTEPH in APTE survivors is about 3-9%, and 2 years survival in untreated patients with a mean pulmonary artery pressure greater than 50 mmHg was as low as 10% (Simonneau et al., 2018). However, the early recognition of CTEPH progression is difficult due to the gradual beginning and absence of effective biomarker Miao et al., 2017). IPAH is a rare but life-threatening cardiopulmonary disease of unknown etiology, with progressive loss of quality of life, showing a 5 years overall survival rate of 50% (Berger et al., 2011;Simonneau et al., 2018;Koudstaal et al., 2020). IPAH characterizes by progressively increased pulmonary artery pressure and pulmonary vascular resistance, and up to now, there is no specific treatment (Santos-Ferreira et al., 2020;Zhang et al., 2020). Thus, comprehensive analyses of potential mechanisms with possible implications in time of evolution and disease form in APTE, CTEPH and IPAH are urgently needed. Interestingly, they share common pathophysiology changes, including fibrinolysis, pulmonary artery endothelial cell (PAEC) proliferation; pulmonary artery smooth muscle cell (PASMC) proliferation, migration, and contraction; inflammation; and pulmonary arterial adventitial fibroblasts (PAAF) proliferation, activation, and migration. Therefore, there should be some differences between these three pulmonary vascular diseases at molecular level, especially microRNAs (miRNAs).
MicroRNAs (miRNAs) are small, non-coding endogenous RNA molecules consisting of approximately 20 nucleotides that suppress gene expression post-transcriptionally by binding to the ''seed sequences'' in 39 untranslated regions (UTRs) mRNAs (Chen et al., 2008;Kosaka et al., 2010). Identified as stable cell-free miRNAs in serum or plasma, circulating miRNAs are discharged passively and selectively to blood by various cells, and may play as receiver or messenger in cell communication (Kosaka et al., 2010). Through illness, miRNAs in the sick cells are released into the circulation, and the circulating miRNA profile is adopted with the disease characteristics (Chen et al., 2008). Therefore, circulating miRNAs have been widely explored as potential blood-based biomarkers for individual cardiovascular diseases diagnosis. Furthermore, it is also reasonable to assume that identification of molecular pathways and axon guidance enriched by differentially expressed miRNAs target genes and interactions between them provide a molecular basis for the unbalance between hyperproliferative and apoptosis-resistant phenotype particularly in the three principal pulmonary artery cell involved including PAEC, PASMC, and PAAF (Graham et al., 2010;Chen and Jiang, 2018;Yuan, 2019;Santos-Ferreira et al., 2020). Thus, we hypothesize that analyzing different circulating plasma miRNAs in APTE, CTEPH and IPAH, we can describe the early APTE evolution to CTEPH, as well as irreversible fibrocellular narrowing of pulmonary arteries in diseases form, including CTEPH and IPAH.
This study aimed to describe the APTE, CTEPH, and IPAHassociated miRNAs. miRNAs differentially expressed in APTE, CTEPH, and IPAH samples, compared with healthy samples were identified, and the target genes were predicted. The target genes of the key differentially expressed miRNAs were analyzed, and functional enrichment analyses were carried out. The miRNAs were detected using RT-PCR. Finally, we incorporated plasma circulating miRNAs in baseline and clinical characteristics of the patients to detect differences between APTE and CTEPH in time of evolution, as well as differences between CTEPH and IPAH in diseases form.  (Simonneau et al., 2019). Patients with CTEPH were inoperable and managed with targeted IPAH therapy according to the 2018 guideline (Simonneau et al., 2019). Table 1 summarizes the baseline and clinical characteristics of our NGS cohort and healthy subjects.

Blood Plasma Collection and RNA Isolation
Blood samples (6 ml each) were collected into EDTA anticoagulant tubes and mixed by upside-down rotation ten times. After a twostep centrifugation process (1,100 g at 4°C for 10 min, followed by 1,300 g at 4°C for 20 min), the plasma was transferred into rnase/ DNase-free tubes and stored at −80°C until further analysis. We obtained RNA plasma using the miRNeasy Serum/Plasma kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions. The RNA quantity and integrity were determined by Qubit ® 3.0 Fluorometer (Invitrogen, Life Technologies, CA, United States) and the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States), respectively.

Library Preparation and High-Throughput Sequencing
Small RNA library was built with TruSeq Small RNA Library Prep kit (Illumina Inc, San Diego, CA, United States) following the manufacturer's recommendations. We transcribed the RNA to cDNA and amplified it using a common primer and a primer containing the barcode index sequences. This cDNA was separated by PAGE 6%. Libraries was recovered via purification and quality checked on Bioanalyzer 2100 (Agilent Technologies). Library concentration was confirmed (>2 nmol/L) by qPCR using the KAPA Library Quantification kit (KAPA Biosystems). Indexed small RNA libraries were multiplexed in equimolar amounts, denatured, and loaded on flow-cell lanes for sequencing. All samples were sequenced at the Institute of Biotechnology (IBTEC-UNESP, Botucatu, Brazil) using the Illumina NextSeq 500 according to the manufacturer's instructions consisting of 36bp single-end reads in 50 cycles.

Screening of Differentially Expressed miRNAs
Small RNA sequencing raw data was translated into a raw FASTQ sequence one. Then, we filtered the FASTQ format reads to remove rows with low quality via adaptor sequences and sequences inclusion with a quality score <20, using Fastx-Toolkit (http://hannonlab.cshl.edu/fastxtoolkit/). We aligned the filtered reads against the reference human genome, version hg19, using Bowtie v2.0.6. Then, the reads from each individual were merged to identify mature miRNAs using miRDeep2.0.0.8. In detail, the reads were mapped and quantified using the miRDeep2 mapper. pl module and quantifier. pl module, respectively, with standard parameters. Mature miRNAs were downloaded from mirbase v21 (http://www.mirbase.org). The differentially expressed miRNAs (DEmiRNAs) between the samples were normalized and identified by the Bioconductor/R DESeq2 package using ten counts per miRNA as a detection threshold. The criteria for microRNA expression level measuring was CPM (counts per million) value, and we adjusted p values to control the false discovery rate (FDR). We identified significant DEmiRNAs using log2 (fold change) >1 and adjusted FDR.

Functional Enrichment Analysis for Differentially Expressed miRNAs
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to assess the functional annotation of the DEmiRNAs and their enrichment in different biological pathways, respectively. Functional annotation and miRNAs pathways were performed using the DIANA-miRPath v3.0 database according to the published instructions (Vlachos et al., 2015). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to assess the functional annotation of the DEmiRNAs and their enrichment in different biological pathways, respectively. The functional annotation and miRNAs pathways were performed using the DIANA-miRPath v3.0 database according to the published instructions (Vlachos et al., 2015).

Construction of the Differentially Expressed miRNAs-Gene Regulatory Network
In order to predict the DEmiRNA target genes (Hamberg et al., 2016), we used the miRTargetLink (www.ccb.uni-saarland.de/ mirtargetlink) to predict the DEmiRNA target genes (Hamberg et al., 2016). Then, we calculated and visualized the mirtarbase and miRanda interaction network information (Enright et al., 2003;Chou et al., 2016). Thus, we only considered the validated targets with strong evidence. The Venn diagram indicating the intersected targets genes was generated by a Draw Venn Diagram online tool (http://bioinformatics.psb.ugent.be/webtools/Venn/).

Quantitative Reverse-Transcription PCR Confirmation
For quantitative real-time reverse transcription PCR (qRT-PCR) confirmation, we selected the miRNAs with the highest significant differences in expression revealed by next-generation sequencing experiment. Then, we performed the quantification of miRNA expression using miScript miRNA PCR primer assays and miScript SdYBRGreen PCR kit (Qiagen, Valencia, CA, United States) on a StepOnePlus ™ Real-Time PCR System (Applied Biosystems, Foster, United States) following the manufacturer's instructions. Finally, we calculated each miRNA levels using the comparative cycle threshold method (ddCT method) relative to the mean of miR-16 and miR-483. All samples were performed in triplicate.

Statistical Analysis
The chi-square or Fisher's exact test were performed to assess differences in clinical variables. In addition, General Linear Model was applied to evaluate the associations between miRNAs and clinically related parameters. The statistical software program IBM SPSS (version 22; Armonk, NY, United States of America) performed all analyses. p value < 0.05 was considered as statistically significant.

Screening of Differentially Expressed miRNA
Overall, 299 known mature miRNAs were identified in our NGS cohort according to the QC inclusion criteria ( Supplementary  Table S2; Figures 1A,B). Among these, a total of 21 differentially expressed miRNAs were obtained according to the cut-off criteria (p adj <0.05, |log 2 (foldchange)| >1) (Supplementary Table S2).
The heat map of these differentially expressed miRNAs is shown in Figure 2.

Detection of miRNAs Using RT-PCR
The expression of let-7i-5p and miRNA-320a in APTE samples was significantly higher than that of the control samples (p 0.01 and p 0.02, respectively; Figures 6A,B). In contrast, the expression of let-7i-5p and miRNA-320a was significantly lower in CTEPH compared to the control samples (p 0.02 and p 0.01, respectively; Figures 6C,D). IPAH samples presented significant higher expression of let-7i-5p than control samples (p < 0.01; Figure 6E) and significant lower expression of miRNA-320a compared to the control samples (p 0.02; Figure 6F). Table 1 summarizes the baseline and clinical characteristics of patients and healthy subjects. Patients and healthy subjects presented a similar distribution of age, gender, race, and tobacco history. Deep venous thrombosis was detected in six out of 47 (13%) patients. In patients who underwent right cardiac catheterization, the level of mean pulmonary artery pressure (mPAP) in IPAH was significantly higher than CTEPH. Equally significant was the increased systolic pulmonary artery pressure (sPAP) in IPAH compared to CTEPH. Furthermore, pulmonary vascular resistance (PVR) was significantly increased in IPAH compared to patients with CTEPH. Oxygen saturation and heart rate were lower in IPAH patients than CTEPH groups. The WHO functional class major than two was significantly associated with CTEPH patients. Tricuspid annular plane systolic excursion (TAPSE) done by echocardiography, was reduced in patients with IPAH when compared to APTE. Brain natriuretic peptide (BNP) or its stable inactive pro-FIGURE 2 | Circulating miRNAs signature in APTE, CTEPH, PAH patients, and healthy controls. Hierarchical clustering and heatmap diagram of the optimum circulating miRNAs from APTE (n 7), CTEPH (n 8), PAH (n 5), and healthy controls (n 2). Green and red colors represent increased and decreased expression levels (transformed into z-scores), respectively, when compared with the average level of each miRNA. Columns and lines represent samples and miRNA, respectively. APTE, Acute Pulmonary Thromboembolism; CTEPH, Chronic Thromboembolic Pulmonary Hypertension; PAH, Pulmonary artery Hypertension; and CTR, Controls.

Baseline and Clinical Characteristics of Patients and Healthy Subjects
Frontiers in Pharmacology | www.frontiersin.org May 2021 | Volume 12 | Article 648769 hormone (NT-proBNP) was higher in CTEPH than IPAH patients. Patients with APTE presented a significantly increased level of CRP and D-dimer compared to CTEPH and PAH.

Association Between Circulating Plasma miRNA Expression, Clinical and Baseline Characteristics of the Patients
Differentially expressed miRNAs were associated with baseline and clinical characteristics of the patients, as shown in Table 3. Expression of let-7i-5p > 2.30 and miRNA-320a >2.58 was associated with APTE, high level of C-reactive protein and D-dymer. In contrast, miRNA-320a expression ≤2.58 was associated with CTEPH (p 0.04) and IPAH (p 0.03).

Rationale of the Study
We described five highly expressed circulating plasma miRNAs present in APTE, CTEPH and IPAH assembled in one conglomerate. APTE, CTEPH and IPAH are rather different pulmonary vascular diseases in terms of clinical presentation, prevalence, pathophysiology, and prognosis, making it difficult to extrapolate the results among them. However, incorporating miRNA panel into baseline and clinical characteristics of the patients, we provide additional information about outcome evolution (APTE and CTEPH) and the disease form (CTEPH and IPAH), which may have implications in life-threatening cardiopulmonary disease, and specific surgical or medical treatment. To our knowledge it is an innovative study suggesting a differential spectrum of evolution and form of the disease among APTE, CTEPH and IPAH. Among the miRNA 5-signature, let-7i-5p and miR-320a were the most expressed and functionally enriched by cancer, glioma, and PI3K-Akt signaling pathways. The network construction for target genes showed 11 genes regulated by let-7i-5p and 20 genes regulated by miR-320a. These genes regulated by let-7i-5p and miR-320a are modulators of PAAF, PAEC, and PASMC, respectively involved in proliferation, vasoconstriction, inflammation, and DNA damage. These target genes regulated by differentially expressed miRNAs were mainly enriched in GO terms and KEGG pathways. AR (androgen receptor), a target gene of hsa-let-7i-5p and has-miR-320a, was enriched in pathways in cancer, whereas PRKCA (Protein kinase C Alpha), also a target gene of hsa-let-7i-5p and has-miR-320a, was enriched in KEGG pathways, such as pathways in cancer, glioma, and PI3K-Akt signaling pathway.
Taken together, the interaction between genes and their miRNA may play important roles in the time of evolution and disease form in pulmonary hypertensive and thromboembolic diseases, as discussed below.

Outcome Evolution in Pulmonary Thromboembolism (APTE and CTEPH)
Our results demonstrated that miR-let-7i-5p expression was upregulated in APTE and downregulated in CTEPH. Furthermore, these upregulated miRNAs in APTE were significantly enriched in cancer, glioma, and PI3K-Akt signaling pathways, whereas the downregulated miRNAs observed in CTEPH were mainly enriched in cancer and apelin signaling pathways.
Yang's group showed that 24 miRNAs were upregulated, and 22 miRNAs were downregulated (including miR-let-7i-5p) in CTEPH compared with APTE (Miao et al., 2017). Likewise, Guo FIGURE 3 | Heatmap and box plots graphs representing the five miRNAs-let-7i-5p, miR-320a, miR-320b-1, miR-320b-2, and miR-1291 significantly associated with the APTE group compared to the CTEPH group after testing adjustment by FDR (Benjamini-Hochberg). Red and blue colors represent increased and decreased expression levels (transformed into z-scores), respectively, when compared with the average level of each miRNA. APTE, acute pulmonary thromboembolism; CTEPH, thromboembolic pulmonary hypertension; PAH, idiopathic pulmonary hypertension; and CTR, Controls. et al. also found that miR-let-7i-5p was downregulated in CTEPH patients when compared to healthy controls, as well as miR-140-3p, miR-93, miR-22, miR-106b; and that miR-1260, miR-602, miR-129-5p, miR-1908 and miR-483-5p were upregulated . As a tumor suppressor or tumor promoter, miR-let-7b has been found to halt cell proliferation, adhesion, and invasion by targeting Protein Code Gene (PKA1), Diaphanous Related Formin 2 (DIAPH2), Radixin (RDX), RAS, MYC, and High Mobility Group A2 (HMGA2) genes and their respective proteins. Let-7b expression is also differently regulated in breast cancer and the data reveal its possible role in DNA repair capacity during breast carcinogenesis (Hu et al., 2013). Similarly, miR-let-7b was shown to be significantly downregulated in osteosarcoma tissues and cell lines and the functional studies revealed that the antitumor effect of miR-let-7b was probably due to targeting and suppressing Insulin Like Growth Factor 1 Receptor (IGF1R) expression (Zhang et al., 2019). Wang et al. demonstrated that miR-let-7b binds the 3′-UTR of Angiopoietin (ANGPT2) to induce migration and tube formation of Human Umbelical Vein Endothelial Cells (HUVECs), and Human Cardiac Progenitor Cells (hCPFs)-exosome transport to endothelial cells expressing relatively low amounts of miR-let-7b, promoting angiogenesis through upregulation of Angiopoietin 2 (ANGPT2) (Wang et al., 2015). MiR-let-7i-5p has been reported as pivotal in regulating cell proliferation and migration. In fact, investigations have demonstrated that PASMCs excessively proliferate and migrate to the intima layer of pulmonary artery after APTE, remodeling the pulmonary vasculature and increasing pulmonary vascular resistance (Douma et al., 2010;Goldhaber et al., 2012;Konstantinides et al., 2014). Thus, we hypothesize that upregulation of miR-let-7i-5p may be intimately associated with excessive proliferation and migration of PASMCs following APTE which may evolve to CTEPH. Ham et al. showed that the levels of miR-let-7i-5p were upregulated in APTE compared with CTEPH, thus suppressing apoptosis via caspase in HUVECs (Ham et al., 2015). In this way, we have also shown that 11 genes were a direct target of miR-let-7i-5p, which upregulation modulates the proliferation of human PASMCs. Among them, Interleukin13 (IL13) upregulation induces PASMCs in the pulmonary artery through signal transducer and activator of transcription STAT6 or through phosphatidylinositol 3-kinase (PI3K), STAT3 and mitogen activated protein kinase (MAPK) pathways (Yuan, 2019). Oxidant, arginase 2 (ARG2) and hypoxia-inducible factor 1α (HIF1) are also involved in the proliferation of PASMCs (Yuan, 2019). The experimental Schistosomiasisassociated PAH potentially depends on upregulation of IL-13, which can activate TGF-β1 via phospho-STAT6 signaling, or via increased matrix metalloproteinase-inducing muscular hypertrophy and periadventitial fibrosis in pulmonary arteries (Graham et al., 2010). The human eukaryotic initiation factor EIF2C1 gene downregulates protein synthesis in response to amino acid starvation, hypoxia, and viral infection but upregulated specific stress response proteins. EIF2C1 expression is increased in smooth muscle cells in the vessel wall and interstitial tissue (Chen and Jiang, 2018). Toll-Like Receptor four Signaling (TLR4) has a prominent functional impact on platelet activity, hemostasis, and thrombosis (Vallance et al., 2017). Bone morphogenetic protein-4 (BMP4) acts downstream of HIF-1 and mediates hypoxia-induced upregulation of transient receptor potential canonical (TRPC), leading to increased basal [Ca2+]i in pulmonary arterial smooth muscle cells (PASMCs), promoting chronic hypoxiainduced pulmonary hypertension pathogenesis (Wang et al., 2015).
The up-regulation in miR-let-7i-5p observed only in APTE and the downregulation of miR-let-7i-5p observed only in CTEPH reinforces the abnormal remodeling between APTE and CTEPH. These differences may be associated with time of evolution, since differences were observed in cardiac index, C-reactive protein and D-dimer.

Disease Form in Pulmonary Artery Hypertension (CTEPH and IPAH)
IPAH is an isolated small-vessel disease comprising vasoconstriction, remodeling, and thrombosis of distal vessels. On the other hand, CTEPH is characterized by large-vessel occlusion with a concomitant small-vessel disease following thromboembolism event. The pulmonary microvascular disease as a common denominator of IPAH and CTEPH probably plays a major role in pathophysiology (Berger et al., 2011;Koudstaal et al., 2020).
In the current study, we carried out NGS analysis and detection of miRNAs to understand the key miRNAs associated with disease form in CTEPH and IPAH. We found that let-7i-5p and miR-320a were down-regulated in CTEPH, while let-7i-5p was up-regulated in IPAH. Furthermore, 20 target genes were coregulated by miR-320a and other miRNAs. Androgen receptor (AR), a target gene of hsa-miR-320a, was enriched in cancer pathways. Protein kinase C Alpha (PRKCA), also a target gene of hsa-let-7i-5p and has-miR-320a, was enriched in 14 of 18 KEGG pathways, such as in cancer, glioma, and PI3K-Akt signaling pathway.
Recent studies have shown that miRNA dysregulation plays an important role in the hyperproliferative and apoptosis-resistant phenotype of pulmonary vascular cells in PH, including PAAF, PAEC, and PASMC respectively involved in proliferation, vasoconstriction, inflammation, and DNA damage (Santos-Ferreira et al., 2020). Previous studies have also uncovered that dysregulation of miRNAs such as miR-204 and miR-21 was associated with the pathobiology of IPAH (Courboulin et al., 2011;Rhodes et al., 2013;Bienertova-Vasku et al., 2015), suggesting that miRNAs may configure novel therapeutic targets. In this line, the downregulation of miR-let-7i-5p has also been observed elsewhere (Wang et al., 2013). While miR-let-7i-5p was downregulated in CTEPH, it was upregulated in IPAH and this inverse association may be crucial to differentiate both disease forms. As a tumor promoter, miR-let-7b has been found to enhance cell proliferation, adhesion, and invasion by targeting the genes of PKA1, DIAPH2, RDX, RAS, MYC, and HMGA2 and respective proteins, creating an unbalance between hyperproliferative and apoptosis-resistant phenotype of PAAF, PAEC and PASMC in IPAH. Moreover, in the present study, we demonstrated that 20 genes were a direct target of miR-320a. It has been reported these miRNA regulates cell proliferation, differentiation, and invasion in human colon cancer (Sun et al., 2012), and erythroid differentiation (Mittal et al., 2013) as well as bladder carcinoma (Shang et al., 2014) and when upregulated, miRNA-320a act toward proliferation of human PASMCs. Among the 20 target genes, a family of water-selective membrane channels promoting endothelial cell migration and angiogenesis (aquaporin-AQP) was found to be overexpressed in PAH patients, indicating its potential role in the treatment of PH (Saadoun et al., 2005;Schuoler et al., 2017;Gräf et al., 2018).
Considering the above evidence, we inferred that unbalance between the hyperproliferative and apoptosis-resistant phenotype of PAAF, PAEC and PASMC in pulmonary artery confer differences in IPAH and CTEPH about the diseases form. The incorporation of miRNA-320a and the unbalance of miR-let-7i-5p in baseline added by clinical characteristics signatures reinforce differences between CTEPH and IPAH.

CONCLUSION
We concluded that circulating plasma miRNA and clinical/ hemodynamic characteristics provide additional predictive information about the evolution outcome and the form of pulmonary thromboembolic and hypertensive diseases.

DATA AVAILABILITY STATEMENT
We have completed the deposition of our BioProject database in a public, community-supported repository, under the number PRJNA690230 which can be accessible with the following link: http://www.ncbi.nlm.nih.gov/bioproject/690230.  Abbreviations: APTE, acute pulmonary thromboembolism; CPTEH, chronic pulmonary thromboembolism with pulmonary hypertension; IPAH, idiopathic pulmonary hypertension; mean pulmonary arterial pression sPAH, systolic pulmonary arterial pressure; PVR, pulmonary vascular resistance; TAPSE, tricuspid annular plane systolic pressure; CI, cardiac index; CRP, C-reactive protein. General Linear Model was carried out to analyze the effect of several dependent variables on clinical features of the subjects. β β coefficient. P P value. a FDR statistical test (statistical method for correcting the p value for multiple tests) was carried out to compare independent variables between two groups.
Bolded values are statically significant.