Thrombosis-related circulating miR-16-5p is associated with disease severity in patients hospitalised for COVID-19.

SARS-CoV-2 tropism for the ACE2 receptor, along with the multifaceted inflammatory reaction, is likely to drive the generalized hypercoagulable and thrombotic state seen in patients with COVID-19. Using the original bioinformatic workflow and network medicine approaches we reanalysed four coronavirus-related expression datasets and performed co-expression analysis focused on thrombosis and ACE2 related genes. We identified microRNAs (miRNAs) which play role in ACE2-related thrombosis in coronavirus infection and further, we validated the expressions of precisely selected miRNAs-related to thrombosis (miR-16-5p, miR-27a-3p, let-7b-5p and miR-155-5p) in 79 hospitalized COVID-19 patients and 32 healthy volunteers by qRT-PCR. Consequently, we aimed to unravel whether bioinformatic prioritization could guide selection of miRNAs with a potential of diagnostic and prognostic biomarkers associated with disease severity in patients hospitalized for COVID-19. In bioinformatic analysis, we identified EGFR, HSP90AA1, APP, TP53, PTEN, UBC, FN1, ELAVL1 and CALM1 as regulatory genes which could play a pivotal role in COVID-19 related thrombosis. We also found miR-16-5p, miR-27a-3p, let-7b-5p and miR-155-5p as regulators in the coagulation and thrombosis process. In silico predictions were further confirmed in patients hospitalized for COVID-19. The expression levels of miR-16-5p and let-7b in COVID-19 patients were lower at baseline, 7-days and 21-day after admission compared to the healthy controls (p < 0.0001 for all time points for both miRNAs). The expression levels of miR-27a-3p and miR-155-5p in COVID-19 patients were higher at day 21 compared to the healthy controls (p = 0.007 and p < 0.001, respectively). A low baseline miR-16-5p expression presents predictive utility in assessment of the hospital length of stay or death in follow-up as a composite endpoint (AUC:0.810, 95% CI, 0.71-0.91, p < 0.0001) and low baseline expression of miR-16-5p and diabetes mellitus are independent predictors of increased length of stay or death according to a multivariate analysis (OR: 9.417; 95% CI, 2.647-33.506; p = 0.0005 and OR: 6.257; 95% CI, 1.049-37.316; p = 0.044, respectively). This study enabled us to better characterize changes in gene expression and signalling pathways related to hypercoagulable and thrombotic conditions in COVID-19. In this study we identified and validated miRNAs which could serve as novel, thrombosis-related predictive biomarkers of the COVID-19 complications, and can be used for early stratification of patients and prediction of severity of infection development in an individual.Abbreviations: ACE2, angiotensin-converting enzyme 2AF, atrial fibrillationAPP, Amyloid Beta Precursor ProteinaPTT, activated partial thromboplastin timeAUC, Area under the curveAß, amyloid betaBMI, body mass indexCAD, coronary artery diseaseCALM1, Calmodulin 1 geneCaM, calmodulinCCND1, Cyclin D1CI, confidence intervalCOPD, chronic obstructive pulmonary diseaseCOVID-19, Coronavirus disease 2019CRP, C-reactive proteinCV, CardiovascularCVDs, cardiovascular diseasesDE, differentially expressedDM, diabetes mellitusEGFR, Epithelial growth factor receptorELAVL1, ELAV Like RNA Binding Protein 1FLNA, Filamin AFN1, Fibronectin 1GEO, Gene Expression OmnibushiPSC-CMs, Human induced pluripotent stem cell-derived cardiomyocytesHSP90AA1, Heat Shock Protein 90 Alpha Family Class A Member 1Hsp90α, heat shock protein 90αICU, intensive care unitIL, interleukinIQR, interquartile rangelncRNAs, long non-coding RNAsMI, myocardial infarctionMiRNA, MiR, microRNAmRNA, messenger RNAncRNA, non-coding RNANERI, network-medicine based integrative approachNF-kB, nuclear factor kappa-light-chain-enhancer of activated B cellsNPV, negative predictive valueNXF, nuclear export factorPBMCs, Peripheral blood mononuclear cellsPCT, procalcitoninPPI, Protein-protein interactionsPPV, positive predictive valuePTEN, phosphatase and tensin homologqPCR, quantitative polymerase chain reactionROC, receiver operating characteristicSARS-CoV-2, severe acute respiratory syndrome coronavirus 2SD, standard deviationTLR4, Toll-like receptor 4TM, thrombomodulinTP53, Tumour protein P53UBC, Ubiquitin CWBC, white blood cells.

NERI: network-medicine based integrative approach for disease gene prioritization by relative importance.

BACKGROUND: Complex diseases are characterized as being polygenic and multifactorial, so this poses a challenge regarding the search for genes related to them. With the advent of high-throughput technologies for genome sequencing, gene expression measurements (transcriptome), and protein-protein interactions, complex diseases have been sistematically investigated. Particularly, Protein-Protein Interaction (PPI) networks have been used to prioritize genes related to complex diseases according to its topological features. However, PPI networks are affected by ascertainment bias, in which more studied proteins tend to have more connections, degrading the results quality. Additionally, methods using only PPI networks can provide only static and non-specific results, since the topologies of these networks are not specific of a given disease. RESULTS: The goal of this work is to develop a methodology that integrates PPI networks with disease specific data sources, such as GWAS and gene expression, to find genes more specific of a given complex disease. After the integration of PPI networks and gene expression data, the resulting network is used to connect genes related to the disease through the shortest paths that have the greatest concordance between their gene expressions. Both case and control expression data are used separately and, at the end, the most altered genes between the two conditions are selected. To evaluate the method, schizophrenia was adopted as case study. CONCLUSION: Results show that the proposed method successfully retrieves differentially coexpressed genes in two conditions, while avoiding the bias from literature. Moreover we were able to achieve a greater concordance in the selection of important genes from different microarray studies of the same disease and to produce a more specific gene set related to the studied disease.

Alteration of circulating ACE2-network related microRNAs in patients with COVID-19.

Angiotensin converting enzyme 2 (ACE2) serves as the primary receptor for the SARS-CoV-2 virus and has implications for the functioning of the cardiovascular system. Based on our previously published bioinformatic analysis, in this study we aimed to analyze the diagnostic and predictive utility of miRNAs (miR-10b-5p, miR-124-3p, miR-200b-3p, miR-26b-5p, miR-302c-5p) identified as top regulators of ACE2 network with potential to affect cardiomyocytes and cardiovascular system in patients with COVID-19. The expression of miRNAs was determined through qRT-PCR in a cohort of 79 hospitalized COVID-19 patients as well as 32 healthy volunteers. Blood samples and clinical data of COVID-19 patients were collected at admission, 7-days and 21-days after admission. We also performed SHAP analysis of clinical data and miRNAs target predictions and advanced enrichment analyses. Low expression of miR-200b-3p at the seventh day of admission is indicative of predictive value in determining the length of hospital stay and/or the likelihood of mortality, as shown in ROC curve analysis with an AUC of 0.730 and a p-value of 0.002. MiR-26b-5p expression levels in COVID-19 patients were lower at the baseline, 7 and 21-days of admission compared to the healthy controls (P < 0.0001). Similarly, miR-10b-5p expression levels were lower at the baseline and 21-days post admission (P = 0.001). The opposite situation was observed in miR-124-3p and miR-302c-5p. Enrichment analysis showed influence of analyzed miRNAs on IL-2 signaling pathway and multiple cardiovascular diseases through COVID-19-related targets. Moreover, the COVID-19-related genes regulated by miR-200b-3p were linked to T cell protein tyrosine phosphatase and the HIF-1 transcriptional activity in hypoxia. Analysis focused on COVID-19 associated genes showed that all analyzed miRNAs are strongly affecting disease pathways related to CVDs which could be explained by their strong interaction with the ACE2 network.

ACE2 Interaction Networks in COVID-19: A Physiological Framework for Prediction of Outcome in Patients with Cardiovascular Risk Factors.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (coronavirus disease 2019; COVID-19) is associated with adverse outcomes in patients with cardiovascular disease (CVD). The aim of the study was to characterize the interaction between SARS-CoV-2 and Angiotensin-Converting Enzyme 2 (ACE2) functional networks with a focus on CVD. METHODS: Using the network medicine approach and publicly available datasets, we investigated ACE2 tissue expression and described ACE2 interaction networks that could be affected by SARS-CoV-2 infection in the heart, lungs and nervous system. We compared them with changes in ACE-2 networks following SARS-CoV-2 infection by analyzing public data of human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs). This analysis was performed using the Network by Relative Importance (NERI) algorithm, which integrates protein-protein interaction with co-expression networks. We also performed miRNA-target predictions to identify which miRNAs regulate ACE2-related networks and could play a role in the COVID19 outcome. Finally, we performed enrichment analysis for identifying the main COVID-19 risk groups. RESULTS: We found similar ACE2 expression confidence levels in respiratory and cardiovascular systems, supporting that heart tissue is a potential target of SARS-CoV-2. Analysis of ACE2 interaction networks in infected hiPSC-CMs identified multiple hub genes with corrupted signaling which can be responsible for cardiovascular symptoms. The most affected genes were EGFR (Epidermal Growth Factor Receptor), FN1 (Fibronectin 1), TP53, HSP90AA1, and APP (Amyloid Beta Precursor Protein), while the most affected interactions were associated with MAST2 and CALM1 (Calmodulin 1). Enrichment analysis revealed multiple diseases associated with the interaction networks of ACE2, especially cancerous diseases, obesity, hypertensive disease, Alzheimer's disease, non-insulin-dependent diabetes mellitus, and congestive heart failure. Among affected ACE2-network components connected with the SARS-Cov-2 interactome, we identified AGT (Angiotensinogen), CAT (Catalase), DPP4 (Dipeptidyl Peptidase 4), CCL2 (C-C Motif Chemokine Ligand 2), TFRC (Transferrin Receptor) and CAV1 (Caveolin-1), associated with cardiovascular risk factors. We described for the first time miRNAs which were common regulators of ACE2 networks and virus-related proteins in all analyzed datasets. The top miRNAs regulating ACE2 networks were miR-27a-3p, miR-26b-5p, miR-10b-5p, miR-302c-5p, hsa-miR-587, hsa-miR-1305, hsa-miR-200b-3p, hsa-miR-124-3p, and hsa-miR-16-5p. CONCLUSION: Our study provides a complete mechanistic framework for investigating the ACE2 network which was validated by expression data. This framework predicted risk groups, including the established ones, thus providing reliable novel information regarding the complexity of signaling pathways affected by SARS-CoV-2. It also identified miRNAs that could be used in personalized diagnosis in COVID-19.