key: cord-0911270-7010x3ff
authors: Wilson, Julie C.; Kealy, David; James, Sally R.; Plowman, Tobias; Newling, Katherine; Jagger, Christopher; Filbey, Kara; Mann, Elizabeth R.; Konkel, Joanne E.; Menon, Madhvi; Knight, Sean B.; Simpson, Angela; Prihartadi, Aliya; Forshaw, Greg; Todd, Neil; Yates, David R.A.; Grainger, John R.; Hussell, Tracy; Kaye, Paul M.; Signoret, Nathalie; Lagos, Dimitris
title: Integrated miRNA/cytokine/chemokine profiling reveals severity-associated step changes and principal correlates of fatality in COVID-19
date: 2021-12-20
journal: iScience
DOI: 10.1016/j.isci.2021.103672
sha: 0fd0d065152afb1091a638a2794cb0f537f8d436
doc_id: 911270
cord_uid: 7010x3ff

Inflammatory cytokines and chemokines (CC) drive COVID-19 pathology. Yet, patients with similar circulating CC levels present with different disease severity. Here, we determined 171 microRNAomes from 58 hospitalised COVID-19 patients (Cohort 1) and levels of 25 cytokines and chemokines (CC) in the same samples. Combining microRNA (miRNA) and CC measurements allowed for discrimination of severe cases with greater accuracy than using miRNA or CC levels alone. Severity group-specific associations between miRNAs and COVID-19-associated CC (e.g. IL6, CCL20) or clinical hallmarks of COVID-19 (e.g. neutrophilia, hypoalbuminemia) separated patients with similar CC levels but different disease severity. Analysis of an independent cohort of 108 patients from a different centre (Cohort 2) demonstrated feasibility of CC/miRNA profiling in leftover hospital blood samples with similar severe disease CC and miRNA profiles, and revealed CCL20, IL6, IL10, and miR-451a as key correlates of fatal COVID-19. These findings highlight that systemic miRNA/CC networks underpin severe COVID-19.

The COVID-19 pandemic has caused more than 5.2 million deaths (as of 1 st December 2021, supplemental oxygen support while others develop critical disease requiring invasive ventilation in intensive care units (ICUs) . Age, comorbidities, host genetics, and the type of immune response are all thought to play a decisive role in shaping the trajectory of the infection (Esai Selvan, 2020).

There is a precarious balance between protective immunity and immune-driven pathology in COVID-19, with numerous studies identifying systemic and tissue correlates of severe and fatal disease. Systemic clinical and immune correlates of severity include elevated CRP and D-dimer levels, increased numbers of neutrophils, lymphopenia, T cell exhaustion, increased number of proliferating plasmablasts and extrafollicular B cell activation, and generation of platelet-sequestering non-classical monocytes (Chen et al., 2020a; Diao et al., 2020; Schulte-Schrepping et al., 2020; Terpos et al., 2020; Vanderbeke et al., 2021; Woodruff et al., 2020; Wu et al., 2020) . Consistent with the acute respiratory distress syndrome (ARDS) presentation of severe COVID-19, multiple studies have also reported increased circulating levels of proinflammatory cytokines and chemokines (CC) associated with COVID-19 severity including IL6, IL10, CCL2, CXCL10, and GMCSF (Consortium, 2021; Del Valle et al., 2020; Laing et al., It has been proposed that miRNAs might be involved in controlling inflammation during COVID-19 (Aslani et al., 2021; Gasparello et al., 2021; Houshmandfar et al., 2021; Matarese et al., 2020; Plowman and Lagos, 2021) . Theoretical CC/miRNA regulatory networks have been proposed but not validated (Khokhar et al., 2022) , whereas the potential of miRNA-like therapeutics has also been suggested (Ahn et al., 2021; Dash et al., 2021; Hum et al., 2021; Park et al., 2021) . Of note, initial studies have identified potential changes in miRNA expression in COVID-19 and signatures from targeted testing of 41 miRNAs have been linked to ICU treatment and mortality in a Spanish cohort (de Gonzalo-Calvo et al., 2021) .

Furthermore, it has been reported that suppression of the hsa-miR-320 family is linked to COVID-19 severity (Duecker et al., 2021) . Correlations between specific miRNAs and response to treatment (Sabbatinelli et al., 2021) or D-dimer levels and risk of thrombosis in COVID-19 patients have also been proposed (Gambardella et al., 2021b) . In addition, hsa-miR-2392 has been proposed both as a biomarker and therapeutic target of SARS-CoV-2 infection (McDonald et al., 2021) .

We hypothesised that combining measurements of circulating CC and miRNAs would reveal qualitative changes associated with increased COVID-19 severity, leading to insight into determinants of COVID-19 outcome. We aimed to explore whether integrating miRNA with CC profiles would reveal a set of more accurate correlates of COVID-19 severity and outcome, and to identify severity-specific correlations of miRNAs with cytokine/chemokine and clinical parameter profiles in COVID-19. To achieve this, we measured CC levels as potential initiators of pathogenic signalling and circulating miRNAs as indicators of intracellular gene networks operating in combination with or in response to these CC. We first performed cytokine, chemokine and miRNA profiling in 171 blood plasma samples from 58 hospitalised patients with mild, moderate, or severe COVID-19 (Cohort 1). We found that integrating clinical parameters with systemic correlates of immune response and tissue status (cytokines, chemokines, and miRNAs) revealed step changes occurring between clinical severity stages in COVID-19 and novel correlates of severe COVID-19. These severity group-specific J o u r n a l P r e -p r o o f qualitative features of COVID-19 were unique to each severity group and could only be observed by combining miRNA and CC profiling. Based on these findings, we analysed 166 samples from a second cohort of 108 hospitalised individuals using blood samples that were leftover (or waste) from routine clinical testing following hospital admission (Cohort 2).

Critically, we were able to detect CC and miRNAs in these samples with similar severedisease associated factors. In this cohort, we were also able to search for principal systemic correlates of outcome (death or discharge). Again, this revealed that combinations of CC with miRNAs (e.g. CCL20, IL6, IL10, hsa-miR-451a) were better at predicting outcome than CC or miRNAs alone. Overall, our findings provide novel insight into trajectories underpinning severe disease in symptomatic COVID-19 patients and support that integrating miRNA and cytokine/chemokine profiling can shed light into determinants of disease severity as well as haematological and immunopathological features of COVID-19. Importantly, we demonstrate that miRNA/CC measurements can be reliably performed in leftover blood samples that can be easily accessed in any hospital setting.

We performed circulating miRNA profiling as a nucleic-acid-based approach to discovering peripheral correlates of disease in COVID-19 (Figure S1A). We analysed samples from a first cohort (Cohort 1) of 58 hospitalised patients recruited between April and November 2020 in four hospitals in Manchester, UK as part of the CIRCO cohort (Table S1 ) (Mann et al., 2020) .

These samples were collected for research purposes. Patients were categorised in three groups of disease severity, mild (N=20), moderate (N=21), and severe (N=17), based on the requirement for respiratory support as described before (Mann et al., 2020) . Clinical measurements and comorbidity and symptom data were also available ( Table S2, S3) . For several patients, follow-up plasma samples were analysed, collected at different time-points (typically taken every 1-2 days, Table S4 ).

J o u r n a l P r e -p r o o f 6 Expression of miRNAs was determined by Nanostring profiling and 280 miRNAs were detected in at least 10% of the samples (see Materials and Methods, Table S5-S6) . First, we identified miRNAs that were differentially expressed in samples from individuals with severe disease vs the rest of the cohort (moderate and mild disease). When all samples were included in the analysis (N=171 samples) we identified 32 significantly differentially expressed (DE) miRNAs (adjusted p < 0.1; Figure 1A and Table S5 ). Down-regulated miRNAs in severe disease samples (25 in total, Table S5 ) included hsa-miR-323a-3p, hsa-miR-451a, hsa-miR-188-5p, and hsa-miR-432-5p, while up-regulated miRNAs (8 in total) included hsa-miR-4454+7975 (the mature sequence of hsa-miR-7975 differs only by one base from hsa-miR-4454), hsa-miR-let-7a-5p, hsa-miR-let-7f-5p, hsa-miR-630, hsa-miR-24-3p, and hsa-miR-29b-3p ( Figure 1B , Table S5 ). We repeated the analysis including only the first available sample for each patient (N=58 samples). Although, differences in miRNA expression did not reach significance after correcting for multiple testing, we identified six miRNAs which were DE in severe samples at p<0.05 and absolute log2 fold change (absLFC)>1 ( Figure 1C , D and Table   S5 ). These were hsa-miR-4454+7975, hsa-miR-1285-5p, hsa-miR-320e, and hsa-miR-29b-3p, which were all up-regulated, and hsa-miR-451a and hsa-miR-144-3p, which were downregulated. hsa-miR-4454+7975, hsa-miR-29b-3p and hsa-miR-451a were identified in both analyses (including all samples or first samples). For these DE miRNAs we did not observe any significant changes in expression over the duration of hospitalisation (not shown). Having determined miRNA signatures associated with severe COVID-19 in hospitalised individuals, we performed functional association analysis using the TAM 2.0 tool (Li et al., 2018) . Using the list of DE expressed miRNAs in severe cases (corresponding to Figure 1A) , enriched functional terms for up-regulated miRNAs included cell death, innate immunity, and cell cycle ( Figure 1E ). For down-regulated DE miRNAs, enriched terms included Regulation of Stem Cell, Vascular inflammation, and Chemotaxis ( Figure 1E ). This indicated terms that were likely to reflect the immune response to infection (e.g. Innate immunity) but also terms that could reflect tissue pathology (e.g. cell death). This was also the case, when performing functional term enrichment analysis using the list of all DE miRNAs based on first available samples only J o u r n a l P r e -p r o o f (corresponding to Figure 1B) , with terms referring to wound healing and extracellular matrix remodeling being most enriched ( Figure S1B ). We note that statistical significance was lower due to the smaller number of DE miRNAs. We then explored functional associations of top predicted targets of up-regulated miRNAs in severe samples using the Targetscan tool (Agarwal et al., 2015) (Figure S1C) . Notably enriched gene ontology process associations included terms involving extracellular matrix remodeling and responses to hypoxia consistent with severe COVID-19 pathology ( Figure 1F) . Overall, our analyses revealed a miRNA signature associated with severe COVID-19, which was reflective of tissue pathology and damage.

Progression from mild to moderate and severe COVID-19 is characterised by a global reduction in circulating miRNA levels.

Next, we explored miRNA signatures of mild COVID-19 disease within hospitalised individuals. When performing comparisons, either including all samples or just first available samples, the most notable finding was a decrease in circulating miRNA levels for the vast majority of detectable miRNAs in individuals with severe and moderate disease when compared to individuals with mild disease (Figure 2A-D) . When including all samples, mild cases were characterised by significantly higher levels (adjusted p<0.1) of miRNAs including hsa-miR-495-3p, hsa-miR-150-5p, hsa-miR-432-5p, hsa-miR-376a-3p, hsa-miR-146a-5p, and hsa-miR-451a ( Figure 2B , Table S6 ). When including only the first available sample for each patient differences were more modest (p<0.05, absLFC>1), and DE miRNAs included hsa-miR-548g-3p, hsa-miR-150-5p, hsa-miR-590-5p, hsa-miR-518b, hsa-miR-363-3p, and hsa-miR-495-3p ( Figure 2C -D, Table S6 ). Functional association analysis, using TAM 2.0, indicated an enrichment in miRNAs contributing to immune processes amongst DE miRNAs in mild COVID-19 (corresponding to Figure 2A) , and more generic cellular processes when considering only the first samples ( Figure 2E ). Exploring gene ontology terms amongst the top predicted targets of up-regulated miRNAs (Figure 2A) , we observed enrichment in functional terms associated with mRNA and small RNA metabolism ( Figure 2F ). This could J o u r n a l P r e -p r o o f reflect the type response against SARS-CoV-2, an RNA virus. Alternatively, a change in RNA metabolism, including processing of small RNAs, could be associated with the higher circulating miRNA levels in patients with mild disease. Overall, these analyses revealed a miRNA signature associated with mild COVID-19, reflecting elements of both the immune response to SARS-CoV-2 infection but also tissue pathology.

Cytokine and chemokine profiling reveal quantitative and qualitative differences between COVID-19 severity groups.

In addition to miRNA profiling, we measured cytokine (Figure 3A , B) and chemokine ( Figure   3C , D) levels in the same plasma samples. When including all samples as independent variables we observed significant (adjusted p <0.05) up-regulation of most tested cytokines in severe samples ( Figure 3A , Table S7 ). However, when taking only the first available samples for each patient only up-regulation of IL-6 reached statistical significance when comparing severe samples to moderate and mild ( Figure 3B ). With regards to circulating chemokine levels, CCL3, CXCL9, CCL20, and CXCL1 were significantly up-regulated in samples from individuals with severe COVID-19 when compared against mild and moderate cases ( Figure   3C ). Repeating the analysis to only include first samples, showed again increased levels of CXCL9 and CCL20 in severe disease, but also increased levels of CXCL10 and a downregulation of CCL17 in patients with severe COVID-19 ( Figure 3D ). Comparison of mild samples against the rest of the cohort (moderate and severe) revealed significantly lower levels of CXCL9, CCL20, IL6, IL10, IFNγ, and IL4 in mild samples, whereas levels of CCL5 were significantly higher (Table S7) . Within the studied sampling time-frame we did not observe any statistically significant patterns of regulation for DE CC over time (not shown).

We performed correlation analysis between individual CC, taking each sample as an independent observation ( Figure S2A ) or only including the first available sample for each patient ( Figure S2B ). However, we also performed analyses using the average levels for each patient. This approach has the benefit of capturing a representative and stable cytokine profile per patient, by including multiple sampling points for most patients (Table S4 ) but avoiding over-representation of some patients due to a higher number of collected samples. The abovementioned lack of significant changes in CC levels over the studied period meant that by taking an average of all measurements we were less likely to misrepresent the levels of any specific protein in the correlation matrix. Following this approach, we found that IL2, IL4, IL17A, IL17F and TNF were positively correlated in all three severity groups ( Figure 3E, Figure S2C ), as were IL5 and IFNγ but to a lesser extent particularly in the moderate group ( Figure 3E , Figure   S2A , S2B). Notably, IL10 was also part of this correlation cluster but only for patients with mild disease. Furthermore, CCL5 and CXCL1 were positively correlated in all groups. CXCL11 and CXCL5 were positively correlated with these, but only in mild and moderate disease groups respectively. CXCL10 and CXCL9 also correlated in all groups, although the correlation was stronger in severe cases ( Figure 3E , Figure S2D ). IL9 was positively correlated with CCL3, CCL4, and CXCL8 in severe samples but to a lesser extent, or not at all in moderate and mild samples, respectively ( Figure 3E ). CXCL1 and IL17F were negatively correlated only in moderate and severe samples, the effect being stronger in severe. The same was observed for CXCL10 and CCL17. Notably, we found a strong qualitative difference between the three severity groups with regards to correlation between CCL20 and IFNγ, which were not correlated in mild, negatively correlated in moderate, and positively correlated in severe ( Figure 3F ). To identify relationships between severity groups and high or low levels of cytokines or chemokines we first performed correspondence analysis with count data obtained by assigning values as low, medium or high for each severity group (see Methods). This showed that, elevated IL6 (p<0.001), IL9 (p<0.05), IL10 (p<0.01), and CXCL1 (p<0.05) were significantly associated with the severe disease group, while high values of CCL5 were associated with mild samples ( Figure S2E) . Overall, the CC profiling presented here and correlation analyses per severity group demonstrated that in addition to quantitative, there were also qualitative changes between severity groups.

To gain more insight into differences between severity groups we performed partial least squares regression (PLSR) analysis on the miRNA and CC profiles together (including all samples as independent observations). We found separation between samples from patients with severe disease and those with mild or moderate disease ( Figure 4A and Figure S3A ).

Using PLS models, we performed leave-one-patient-out (L-O-P-O) classification to assess the ability to predict severity based on our dataset. As discrimination between mild and moderate samples was poor, we repeated the analysis with just two classes, severe vs mild/moderate. The analysis was performed using only CC data, all miRNAs, only DE miRNAs, or CC data together with DE miRNAs. We found that combining the CC levels with DE miRNAs resulted in the highest overall accuracy (87%) in severity classification ( Figure 4B ).

To prevent individual patients with multiple samples dominating the analysis, we calculated correlations using average values (one per patient). Notably, correlations between miRNAs and CC calculated across the whole cohort were generally weak ( Figure S3B ). However, when performing the analysis within severity groups stronger and distinct correlations were discovered ( Figure S3C and Table S8 ). The number of these correlations was notably higher in the severe disease group ( Figure S3C , miRNAs and CC that show at least one significant correlation with correlation coefficient higher than 0.4). We observed severity group-specific features. CXCL1 showed positive correlations with a high number of miRNAs in severe disease compared to mild and moderate. CXCL1 also clustered with CXCL9 and CXCL10 in severe disease but not in moderate. IL6 showed several positive correlations with miRNAs and clustered with IL10, CCL20, GMCSF, IL22, and IFNγ in severe samples, and this was not the case for moderate and mild groups (with the exemption of CCL20). Interestingly, IL17F

showed predominantly negative correlations with miRNAs specifically in the severe disease group ( Figure S3C) . Strong, severity-group specific correlations between CCL17 and miR-92a-3p (mild), CXCL1 and miR-1283 (moderate), IL-17F and miR-425-5p (severe), and CXCL9 and miR-1246 (severe) were also observed ( Figure 4C ). It is of note that, although J o u r n a l P r e -p r o o f these correlations were severity group-specific, none of these miRNAs were significantly differentially expressed in severe or mild samples (Table S5-S6) .

For differentially expressed CC in severe or mild samples we determined the five miRNAs with the highest absolute correlation coefficients per severity group (Table S8 ). In addition to severity group-specific top miRNA correlates of plasma CC levels, this analysis revealed a striking qualitative difference for the severe group, namely the increased number of miRNAs featuring as top correlating with multiple CC (Figure 4D ). In severe patients, miR-30d-5p was in the top 5 most correlated miRNAs for 9 CC. Furthermore, hsa-miR-425-5p and hsa-miR-19a-3p feature as top correlates of 8 CC, with hsa-miR-19b-3p featuring a further 6 times (Table S8 ). There were no miRNAs in the top 5 lists of more than 6 CC in mild or moderate cases. Conversely just 27 miRNAs appeared only once as one of the top 5 strongest correlations in severe cases, the number increasing to 39 and 34 for moderate and mild cases, respectively ( Figure 4D ). Of note, we observed distinct top correlating miRNA for each severity group for all tested CC, including for example IL6 and CCL20, two cytokines strongly associated with COVID-19 severity in our data ( Figure 4E ). In addition to different miRNAs, we found several cytokines with opposing correlation trends. For example, all top strongest miRNA correlations for CCL17 were negative in the severe group but positive in the mild and moderate, the reverse being the case for CXCL10. Similarly, the mild group was characterised by all top CCL5 correlations being negative in contrast to moderate and severe groups ( Figure   4F ). Conversely, miRNAs showed group-specific correlation patterns. For example, hsa-miR-24-3p, which was upregulated in patients with severe disease in cohort 1 (Table S3) , showed a negative correlation with TNF and IL4 in the severe cases group only (Table S8) .

Interestingly, the above cross-correlation analysis with miRNA expression (Figure S3C) separated in a severity-dependent manner the IL2, IL4, IL17A, IL17F, TNF cluster found in Figure 3E (when only CC expression was considered). This suggested that co-expression of these pro-inflammatory cytokines might be associated with distinct immunological and/or J o u r n a l P r e -p r o o f tissue pathology effects in each subgroup. We further explored this by focusing on miRNAs that contributed to enrichment of the Cell Death functional association amongst DE miRNAs in severe samples ( Figure 1E ). These miRNAs were miR-24-3p, Let-7f-5p, miR-630, miR-29b-3p, and Let-7a-5p. We reasoned that a set of chemokines or cytokines that would correlate with these miRNAs as a group could potentially be primary candidates to be effectors of lung tissue damage. We used canonical correlation analysis (CCA) to relate a linear combination of the miRNAs associated with cell death with a linear combination of the CC that showed the highest correlations individually with at least one of these miRNAs in severe samples ( Table S9 ). The coefficients in the linear combinations with maximal correlation, obtained using samples from only 12 of the 17 patients in the severe group, are given in Figure   4G . For this data, the correlation for this first canonical pair was 0.8 (p = 4.79x10 -6 ), identifying IL4 and TNF as the predominant contributors to the correlation with Cell Death miRNAs, with minor contributions from IL17F, CCL3, IFNγ, IL9, and IL10. We note here that collinearity due to correlations within these CC may make individual contributions difficult to estimate. Applying the coefficients to the data for samples from the remaining 5 severe patients showed that the model generalises to unseen data with the correlation over all severe samples being 0.76. The canonical variates are plotted for both training and test data in Figure 4G . The above findings identified systemic signatures that were both quantitatively and qualitatively different between COVID-19 severity groups.

patients.

We then explored how these could potentially underpin clinical features of COVID-19 by correlating levels of CC with clinical measurements. Clinical measurements were acquired only at the time of first blood sampling, therefore for this analysis we used CC levels from the first available plasma sample for each participant.

In severe samples, we found a cluster including IL6, IFNγ, CCL2, and CCL20 positively correlating with ALT, CRP, saturated O2, and creatinine and negatively correlating with HCT ( Figure 5A) . A second cluster including TNF, IL17A, IL17F, IL9, IL2, and IL4 correlated positively with several cell populations including neutrophils and eosinophils but largely negatively with monocytes and lymphocytes. This clustering profile was specific to severe cases. Several strong and group-specific correlations were observed. IL9 correlated positively and strongly with lymphocytes in moderate cases, modestly in mild, and not at all in severe ( Figure 5A ). Other group-specific correlations included IL17F with BMI and CCL11 with temperature in mild, CCL11 with MCV and CXCL10 with WCC in moderate, and CCL5 with potassium and CXCL9 with ALP in severe ( Figure 5B) . Interestingly, we noted that cytokines widely associated with COVID-19 severity such IL6, CXCL10, and TNF showed multiple but modest correlations with clinical parameters in all groups, likely reflecting their multifunctional role in disease development.

We also investigated significant correlations between levels of individual miRNAs and clinical measurements. The clinical parameters separated in different clusters between groups and the number of strong correlations (absolute value > 0.4) was higher in moderate cases ( Figure   6A -C). Of note clinical parameters, previously associated with COVID-19 severity, separated in different clusters between groups. For example, albumin and globulin (Feketea and Vlacha, 2020 ) clustered together in mild and moderate samples but not in severe (Figure 6A -C) and the top miRNA correlates of albumin and globulin were distinct between severity groups ( Figure 6D ). In addition, the top miRNA correlates of lymphocyte numbers showed positive correlations in the moderate and severe groups but negative correlations in the mild, whilst we also observed severity group-specific correlations between miRNAs and neutrophils ( Figure 6D) . Interestingly, for CRP ( Figure 6D ) and haemoglobin (Figure S4) , we found switches in the type of correlation and identity of miRNAs between mild and moderate cases, whereas only modest correlations were observed in the severe group ( Figure 6C) . Similar group-specific correlations were observed for all tested clinical parameters ( Figure S4 ). In all groups we observed multiple positive correlations between platelets and individual miRNAs J o u r n a l P r e -p r o o f ( Figure S4 ). However, we did not find any miRNAs that correlated with platelets in all groups.

Eleven miRNAs strongly correlated with platelets in at least two groups. Overall, integration of cytokine, chemokine, and miRNA profiles with clinical parameters identified commonalities but also significant differences between disease severity groups in hospitalised individuals.

Application of CC/miRNA profiling using leftover blood samples from routine clinical testing highlights factors associated with fatal COVID-19.

To test the potential for wider application of miRNA/CC profiling in clinical practice, we analysed an independent COVID-19 patient cohort including 166 leftover blood samples from 108 hospitalised individuals (66 mild, 21 moderate, 21 severe) between January and March 2020 (cohort 2). These samples were collected for clinical assessment and diagnostic purposes (e.g. complete blood counts). Importantly, samples had to be stored at room temperature (typically 12-48h) before plasma separation, to be available in case they were required for further clinical tests. Samples were collected at admission for all patients, with some patients having follow up samples (Table S10) . In this cohort, we measured levels of miRNAs identified as differentially expressed in cohort 1 by qRTPCR, including hsa-miR-4454 (we chose this over hsa-miR-7975 as it is more abundant based on publicly available data in www.miRBase.org), hsa-miR-29b-3p, hsa-miR-451a, hsa-miR-630, hsa-miR-146a-5p, hsa-miR-495-3p, hsa-miR-320e, and hsa-miR-144-3p. We also measured levels of hsa-miR-21-5p as a control miRNA that did not show differential expression between groups in cohort 1 (Figures 1-2, and Tables S5-6) . We opted to measure miRNA levels by qRTPCR as an orthogonal technique to Nanostring profiling, appropriate for candidate miRNA measurement. Figure 2 were also measured. Strikingly, despite storage of these samples at room temperature before plasma isolation, we observed notable similarities in differential expression. As in cohort 1, levels of CCL2, CCL20, IL6, IL9, IL10, IFNγ, IL17A, IL17F, IL4, IL22, hsa-miR-4454, and miR-630 were significantly (adjusted p<0.05) different in severe cases (Figure 7A-C) . It is of note that several cytokine and chemokine receptors are predicted targets of the two miRNAs found to be associated with severe disease in both J o u r n a l P r e -p r o o f cohorts. Predicted hsa- miR-4454 targets include IL10, IL1R1, IL13RA1, IL18BP, IL10RB, and   CCR4, and predicted hsa-miR-630 targets include IL7, IL17D, IL23R, IL15RA, CCR4, CXCL13, CXCL14, CCR2, CCL11, and CCR5. We observed similar core correlations between CC (i.e. correlations occurring in all groups) as in cohort 1, including CXCL10 with CXCL9 (and CXCL11 in cohort 2), IL4, with IL71A, IL17F and TNF (and IL2 to a lesser extent), CXCL1 with CXCL5, and IL6 with IL10 and IL22 (Figure S5A-C) . Similarly to cohort 1, CXCL10 and CCL17 showed a negative correlation only in severe and moderate cases. We note that we did not observe a strong correlation between CCL20 and IFNγ in severe cases, although we still observed the strongest negative correlation between CCL20 and IFNγ in moderate cases as was the case with cohort 1 (compare Figure S5B and Figure 3F ). We did not perform correlation analyses between miRNAs and CC, as we only measured levels of 9 preselected miRNAs rather than the whole miRNome.

The higher number of patients in the cohort 2 allowed us to test CC and miRNA correlates of COVID-19 outcome (discharge or death). Twenty-four patients were deceased and 84 discharged (Table S10). We found that fatal COVID-19 was associated with higher levels of IL6, IL10, IL22, CXCL10, CCL2, and CCL20 and lower levels of miR-495-3p, miR-451a, and miR-29b-3p (Figure 7D-F) . To identify main systemic correlates of fatality, we used forwards stepwise logistic regression to construct models connecting miRNAs and CC with outcome.

We started with a null model and allowed any of the variables found to be statistically significant individually (Figure 7A -F) to be added if the model improved. This resulted in a model including CCL20, IL6, IL10 and hsa-miR-451a ( Figure 7G, Table S11 ), which was compared to a model with only CCL20 (the most statistically significantly associated factor with fatal COVID-19, adjusted p=1.02x10 9 ). Although classification accuracies were comparable between the two models, both the residual deviance (the lack of fit to the perfect model) and the Akaike information criterion (AIC; a measure of prediction error) were smaller for this model, indicating that the CCL20/IL6/IL10/hsa-miR-451a model fitted the data better and suggesting a potential role for these CC and miRNA in fatal COVID-19. A third model that J o u r n a l P r e -p r o o f considered all variables without strong outliers (CXCL10, CCL17, CCL2, CCL3, CXCL9,   CXCL5, CCL20, CXCL11, CCL4, IL6, IL10, TNF, IL17F, IL22, hsa -miR-4454, hsa-miR-29b-3p, hsa-miR-451a, hsa-miR-630, hsa-miR-146a-5p, hsa-miR-495-3p, and hsa-miR-144-3p) resulted in inclusion of 11 variables, including three miRNAs: CCL17, CCL2, CXCL5, CCL20, IL6, IL10, IL17F, IL22, hsa-miR-4454, hsa-miR-495-3p, hsa-miR-451a (Figure 7G-H, Table   S11 ). This 3miRNA/11CC model showed the highest accuracy (79-90%) and lowest deviance (79.204) and AIC (103.2). Taken together, the second cohort analyses demonstrated the feasibility of miRNA and CC profiling in leftover blood samples, identified CC and miRNAs that were statistically elevated in severe cases in both cohorts, and allowed us to build a model highlighting CC and miRNA correlates of fatal COVID-19.

Our study represents the largest circulating miRNA profiling study in COVID-19, including whole miRNAome and multiplex CC analysis of 337 samples from 166 patients from 4 hospitals and two different geographical locations in the UK. Using miRNA profiling to deconvolute CC abundance resulted in identifying sets of miRNAs that correlate in a severityspecific manner with CC already associated with COVID-19 severity such as IL6 and CXCL10 (Coperchini et al., 2021) . Overall, we show that determining correlations between miRNAs and these CC can robustly identify patients with severe disease in a manner not possible through analysis of CC alone. Importantly, these miRNAs reflect intracellular gene networks (e.g. hypoxia, apoptosis, etc) providing mechanistic insight into the contextdependent function of CC in severe COVID-19. It is of note that several of the severity groupspecific miRNA-CC or miRNA-clinical parameter correlations we report involve miRNAs that are not significantly differentially expressed in severe or mild samples. This indicates the emergence of disease-associated novel networks rather than just quantitative increase of interactions operating in a continuum during progression from mild to moderate to severe disease. Interestingly, although such qualitative changes are more numerous and profound when combining miRNA and CC profiling, a small number can be already determined when J o u r n a l P r e -p r o o f exploring CC levels alone. In Cohort 1, IFNγ and CCL20 showed a strong positive correlation only in patients with severe disease. Of note, IFNγ is an essential trigger of antimicrobial pulmonary responses, but increased IFNγ levels in COVID-19 blood and airways samples are associated with pathogenesis and pulmonary tissue damage (Chua et al., 2020; Karki et al., 2021; Montalvo Villalba et al., 2020) , possibly through a detrimental effect of IFNγ on epithelial cells (Heuberger et al., 2021) . On the other hand, CCL20 is secreted by the bronchial epithelium and important for lung immunity, but it can lead to pathogenic responses (Ito et al., 2011) . Therefore, the strongly correlated IFNγ and CCL20 may mark a circulating readout of severe respiratory disease complications.

Our approach can be used as a foundation for stratifying patients and understanding response to CC-targeting treatments such as anti-IL6R blocking antibodies (Galvan-Roman et al., 2021) . Indeed, small-scale studies have proposed candidate miRNA predictors of response to IL6R blockade (Sabbatinelli et al., 2021) . In this respect, testing the applicability of our findings from Cohort 1 when analysing leftover (or waste) blood samples (Cohort 2), which are obtained in every hospital as part of routine care, has wide-reaching implications in clinical practice. Identifying CC and miRNAs with altered levels in patients with severe COVID-19 in both cohorts is notable, as in addition to significant differences in method and time of sampling (e.g immediate processing and freezing, Ficoll-based separation, collection within 1-2 days of admission for Cohort 1 vs storage at room temperature for 12-24h, plasma isolation by centrifugation, and collection at hospital admission in Cohort 2), Cohort 1 was predominantly recruited during the first wave of the pandemic but Cohort 2 was recruited during the second wave of the pandemic in the UK. Therefore, commonalities between the two cohorts span different dominant SARS-CoV-2 variants and changes in clinical practice that occurred between the first two waves of the pandemic. Differences between the two cohorts should be interpreted cautiously as they can be a result of technical or immune-virological differences between the two cohorts.

A correlation between circulating and diseased tissue profiles has been extensively reported in cancer, and autoimmune and infectious diseases (Correia et al., 2017; Sohel, 2020; Wittmann and Jack, 2010; Zeng et al., 2014) . Our analyses in Cohort 1 show that both circulating miRNAs and their predicted targets are indeed reflective of processes occurring at the tissue level (e.g. cell death, response to hypoxia, extracellular matrix remodelling). This provides a powerful tool in further understanding main activities of multifunctional CC. For example, we identify a significant linkage between IL4 and TNF and cell death-associated miRNAs, including hsa-miR-24-3p, which has previously been linked to cerebrovascular effects (Gambardella et al., 2021a; Mone et al., 2021) . These associations do not imply causality, as changes in circulating CC levels can also be the consequence of tissue damage.

Interestingly, hsa-miR-24-3p was also proposed as a potential component of computationally predicted miRNA/CC regulatory networks operating during COVID-19 (Khokhar et al., 2022) , and was with hsa-miR-335-5p one of the only two miRNAs from this study to be highlighted in our analyses. In addition, a smaller study compared 10 COVID-19 patients with 10 healthy controls and identified 55 DE miRNAs in COVID-19 patients (Farr et al., 2021) . Seven of these miRNAs (hsa-let-7a-5p, hsa-let-7f-5p, hsa-miR-150-5p, hsa-miR-423-5p, hsa-miR-432-5p, hsa-miR-451a, hsa-miR-651-5p) were also found to be associated with severity in our study, as they were DE in severe or mild patients. Another study reported down-regulation of hsa-miR-320a,b,c family members when comparing 10 severe cases to 8 healthy controls (Duecker et al., 2021) , something we did not observe in our study. This discrepancy could be due to the small size of the previous study, the fact that the previous study compared COVID-19 patients to healthy controls using a different miRNA detection platform (Nanostring vs small RNA sequencing), or differential expression of family members between cohorts.

A strength of our analysis in Cohort 1, was the ability to perform correlation analysis with clinical laboratory measurements. In severe samples, two CC clusters were observed one of which (TNF, IL17A, IL17F, IL9, IL2, and IL4) correlated positively with neutrophils and largely negatively with monocytes and lymphocytes, providing insight into inflammatory profiles associated with lymphopenia and neutrophilia occurring in severe COVID-19 Terpos et al., 2020) . Of note, several miRNAs correlate with platelets in our analysis.

Although we cannot exclude the possibility that this could be partly due to platelet contamination or lysis, amongst a list of 11 miRNAs previously reported to strongly correlate with platelet contamination (Mitchell et al., 2016) only two (hsa-miR-93-5p and hsa-miR-17-5p) correlated with platelets in more than one severity group in our study. Alternatively, as platelet-secreted miRNAs reflect platelet activation and function (Czajka et al., 2021) , severity group-specific platelet-associated miRNA signatures could provide novel insight into the observed thrombocytopenia and platelet aggregation phenotype in severe COVID-19 Rampotas and Pavord, 2020) .

Although Cohort 1 allowed us to perform unbiased miRNA and CC profiling in high-quality clinical samples, its smaller size meant that it was not possible to attempt analyses of outcome.

Cohort 2 however, allowed us to perform this analysis. The two cohorts could not be merged due to the different methodologies used to measure miRNA levels. Others have reported a link between elevated IL6 and IL10 levels to fatality (Li et al., 2021) providing external validation of our data collected from leftover blood samples. However, these factors alone are not sufficient for prediction of outcome (Liu et al., 2021b) . Our study also revealed CCL20 as key correlate of fatal COVID-19. Increased circulating CCL20 levels in blood as well as respiratory tract samples have been reported in severe COVID-19 patients and linked to ARDS (Chua et al., 2020; Hue et al., 2020; Saris et al., 2021) , but the profound link to fatal COVID-19 observed here is less appreciated. Importantly, combining miRNA with CC profiling leads to improved and more accurate models of correlates of fatal COVID-19 including hsa-miR-451a, hsa-miR-495-3p, and hsa-miR-4454. Of note hsa-miR-451a was also linked to COVID-19 in the above-mentioned comparison of COVID-19 patients to healthy controls (Farr et al., 2021) . These findings provide a solid foundation for other studies further validating fatal COVID-19-associated CC and miRNAs.

Overall, we propose that integrated miRNA/CC profiling can determine key systemic changes that reflect the step-wise progression of COVID-19. Other methodologies such as whole blood RNA sequencing, single cell RNA sequencing, high-resolution proteomics, or autoantibody J o u r n a l P r e -p r o o f profiling also support the existence of distinct disease progression trajectories characterised by unique quantitative and qualitative features (Bergamaschi et al., 2021; Bost et al., 2020; Liu et al., 2021a; Stephenson et al., 2021; Wang et al., 2021) . Identifying qualitative differences between different COVID-19 severity manifestations and quantitative correlates of potentially fatal disease in symptomatic individuals and investigating their mechanistic relevance can provide important clues for development of future COVID-19 therapies. Using exclusively circulating measurements and easy-to-obtain samples to dissect these differences means that our work already provides a set of candidate biomarkers that can be used by the clinical and research community to study pathogenesis and treatment of COVID-19. We propose that targeted tandem miRNA-CC profiling provides a flexible non-invasive toolkit to help mechanistic studies, patient stratification, and response-to-treatment monitoring for future and current COVID-19 therapeutics.

One limitation of our study is that its size does not allow us to perform extensive analysis of all possible confounders (Tables S1-3) and risk analysis to identify independent predictors of COVID-19 severity or outcome. As such future studies in independent cohorts will need to further validate our results. Another limitation of our study is that only patients from the North of England were recruited. Although the severity-associated CC profiles identified here do not significantly differ to those reported in previous studies (Consortium, 2021; Del Valle et al., 2020; Laing et al., 2020; Mann et al., 2020; Thwaites et al., 2021; Zhao et al., 2020) , this geographic limitation should be highlighted. In terms of other biases associated with observational cohort studies (Sedgwick, 2014) , we note that although the number of follow-up samples was not different between severity groups (Table S4) , the duration varied. To counter this bias, we performed some analyses including only the first available sample or the mean value from all samples for each participant. Furthermore, although we cannot exclude selection biases, such effects are likely to be reduced through including only hospitalised COVID-19 patients (all exposed and positive to SARS-CoV-2; no healthy controls) and J o u r n a l P r e -p r o o f studying two cohorts, one of which recruited prospectively. Lastly, as all our measurements were taken from plasma the cellular source of detected CC and miRNAs is not known, therefore it is not possible to assess whether there are direct interactions between the detected miRNAs and CC.

Volcano plot comparing miRNA expression in severe against mild and moderate samples including only the first available samples for each patient. Dotted lines correspond to log2 fold change (LFC) greater than one and Mann Whitney p < 0.05.

Violin plots of DE expressed miRNAs in severe samples when only the first sample per patient is included (p<0.05 and absolute LFC>1).

Significantly over-represented functional terms within up-regulated (shown in black) and

down-regulated (shown in grey) miRNAs in severe cases (shown in red in Figure 1A ), including all samples as independent variables.

F. Significantly over-represented gene ontology process terms amongst predicted targets of miRNAs up-regulated in severe cases (red dots in the top right quadrant in Figure 1A ). A. Violin plots of measured cytokines in mild (green), moderate (blue), and severe (magenta) samples, including all samples. Stars indicate significance for the comparison of severe against mild and moderate as follows: * adjusted p < 0.05, ** adjusted p < 0.01, *** adjusted p < 0.005.

As in A, but analysis including only first available samples per patient C. As in A, but for circulating chemokines D. As in C, but analysis including only first available samples per patient E. Heatmaps showing Spearman correlation between cytokines and chemokines for mild, moderate, and severe groups separately. The average value over all available measurements for each patient was used in the calculation of correlations.

F. Plots of IFNγ against CCL20 levels in mild (green), moderate (blue), and severe (magenta) patients. Correlation coefficients (c.c.) are shown above each scatter plot.

A. Score plots for the first two latent variables from PLSR obtained using cytokine and chemokine data together with DE miRNA measurements. The plot on the left was obtained using all three severity levels (mild, moderate, and severe) whilst that on the right was obtained using only two classes (severe and mild/moderate). B. As in A, but for chemokine levels.

C. As in A, but for miRNA levels.

Violin plots of measured cytokines in samples from discharged (white) and deceased (grey) patients, in the leftover blood sample cohort. Stars indicate significance for the comparison of severe against mild and moderate as follows: * adjusted p < 0.05, ** adjusted p < 0.01, *** adjusted p < 0.005. E. As in D, but for chemokine levels.

F. As in D, but for miRNA levels.

G. Results of forward regression models including accuracy for deceased and discharged patients, residual deviance, and AIC. Results are shown for the CCL20-only model, the CCL20/IL6/IL10/miR-451a model, and the 3 miRNAs/11 CC models (see Table S11 ). 

Requests for resources, datasets, and protocols, should be directed to and will be fulfilled by the lead contact, Dimitris Lagos (Dimitris.lagos@york.ac.uk).

This study did not generate new unique reagents.

 Nanostring miRNA profiling data have been deposited at GEO and are publicly available as of the date of publication. Accession numbers are listed in the key resources table.

 This study does not report original code

Clinical study design and sample collection J o u r n a l P r e -p r o o f

The clinical recruitment and study design for cohort 1 has been reported elsewhere (Mann et al., 2020) . Briefly, individuals with COVID-19 were recruited from 3 hospitals in the Greater Samples were collected as soon after admission as possible and at one-or two-day intervals thereafter (Supplemental Table S4 ).

For cohort 2, the leftover blood cohort, individuals with COVID-19 were recruited at York District Hospital, York and Scarborough NHS Trust between January and March 2020 (Supplemental Table S10 ). Patients were recruited at admission. Ethics approval was obtained from Yorkshire & The Humber -Leeds West Research Ethics Committee (REC reference 19/YH/0394 with approved SA002 amendment of IRAS project 269597). Samples were first used for complete blood counts and stored for 12-48h at room temperature before released for plasma isolation.

For both cohorts, clinical information was extracted from written/electronic medical records.

Only patients who tested positive for SARS-CoV-2 by reverse transcriptase polymerase chain reaction (RT-PCR) on nasopharyngeal/oropharyngeal swabs or sputum were included.

Disease severity was scored based on degree of respiratory failure (Mann et al., 2020) , as below:

Mild: less than 3L/min or less than 28% supplemental oxygen and managed in a ward.

Moderate: less than 10L/min or less than 60% supplemental oxygen, managed in a ward environment; chronic NIV or CPAP (at home) or acute NIV for individuals with COPD.

Severe: any of more than 10L/min or less than 60% supplemental oxygen, or use of acute NIV, or managed in ICU with invasive ventilation.

Where severity of disease changed during admission, the highest disease severity score was selected for classification.

For cohort 1, whole venous blood was collected in tubes containing EDTA (Starstedt) and diluted 1:1 with PBS. Peripheral blood mononuclear cells (PBMCs) were isolated by density gradient centrifugation using Ficoll-Paque Plus (GE Healthcare) and 50 ml SepMate tubes (STEMCELL technologies) according to the manufacturer's protocol. Blood plasma, at the top of the gradient was collected and stored at -80 o C until usage. For cohort 2, whole blood was collected in tubes containing EDTA (Starstedt). Blood plasma was collected by centrifugation (1500g for 10 min at 4 degrees).

For Nanostring analyses (Cohort 1), RNA was extracted from 800μl plasma the Norgen Plasma/Serum RNA Purification Midi kit (Norgen, Cat #56100). A spike-in non-human miRNA (osa-miR-414) was added to each sample (5μl from a 200 pM per sample) for normalisation.

For miRNA profiling we used the Nanostring nCounter™ platform and the human v3 miRNA panel was used containing probes for 827 human miRNAs. Extraction of raw reads and analysis of the individual datasets was carried out using the nSolver™ Analysis Software. For each miRNA the mean plus one standard deviation of negative controls from all counts was subtracted. A minimum value of 1 to all counts. All miRNAs which were expressed above this threshold in fewer than 10% of samples were removed, leaving 280 miRNAs for further analysis. The spike-in count (osa-miR414) was used to normalise miRNA counts. Differential expression analysis on these normalised values was performed by performing Mann Whitney J o u r n a l P r e -p r o o f U test, followed by a false discovery rate correction using the Benjamini Hochberg correction to generate adjusted p values.

For cohort 2, RNA was extracted from 200μl of undiluted plasma using the Norgen Plasma/Serum RNA Purification Mini kit (Norgen, Cat #55000) after addition of the spike RNA at the same concentration as above. Targeted analysis of 9 miRNAs was performed using miRNA Taqman assays (Cat #4427975). Reverse transcription was performed using the Taqman MicroRNA Reverse Transcription Kit (Cat #4366596), by multiplexing primers for all 10 miRNAs (9 miRNAs+spike) after 1 in 6 dilution and using 5μl RNA as template. Quantitative PCR was performed using a StepOnePlus Real Time PCR System (ThermoFisher) and using Taqman Universal PCR Master Mix (Cat #4305719) and standard cycling conditions as per manufacturer's instructions (Taqman miRNA assays).

With the exemption of GM-CSF, for cytokine and chemokine profiling bead-based multiplex flow cytometry assays were used. Cytokines and chemokines levels were quantified from patients' plasma fractions using three sets of pre-defined LEGENDplex panels: Human Th-2 cytokine Panel version 2 (IL2, IL4, IL5, IL6, IL9, IL10, IL13, IL17A, IL17F, IL22, IFNγ and TNF), Proinflammatory chemokine Panel 1 (CCL2, CCL3, CCL4, CCL11, CCL17 and CCL20; CXCL1, CXCL5, CXCL8, CXCL9, CXCL10 and CXCL11) plus a separate CCL5 1-plex kit (BioLegend, San Diego, USA) according to manufacturer's instructions. GM-CSF was detected and quantified using a sandwich ABTS ELISA kit (PeproTech EC, Ltd, UK). Statistical significance in differences in cytokine and chemokine expression between groups were assessed using Mann Whitney U test, followed by Benjamini Hochberg correction to generate adjusted p values.

Data integration and statistical methods J o u r n a l P r e -p r o o f Partial least squares regression (PLSR) was performed using the R package 'pls' with the response variable encoded as 1.0, 2.0 and 3.0 for mild, moderate and severe observations respectively. Leave-one-patient-out (L-O-P-O) classification was carried out by assigning the class corresponding to the closest integer to the output response. For each patient, a model was built without using the data for any of that patient's samples and the results are given for the prediction of the observations left out for each patient in turn. As results showed mild and moderate patients could not be distinguished, the L-O-P-O classification was repeated using just two classes, with mild and moderate observations combined as one class.

For correspondence analysis (CA), min-max scaling was applied to the data for each chemokine and cytokine so that each variable had a minimum of zero and maximum of 3.0 over all observations. With values less than 1.0 considered as 'low', values between 1.0 and 2.0 as 'medium' and values > 2.0 considered as 'high', count data for each severity group was obtained for each variable. To reduce the number of variables, medium counts were excluded from CA, which was performed using the corresp() function from the 'MASS' package in R. All count data was used in chi-squared tests to assess the significance of any association. Balloon plots were created using only counts for either high values or low values with the function ballonplot() from the R package 'gplot'.

The R package 'CCA' was used for canonical correlation analysis (CCA) to find linear combinations of chemokines and cytokines with maximal correlation to a linear combination of the up-regulated miRNAs associated with cell death. To reduce the likelihood of overfitting, we only included cytokines and chemokines with individual correlation to any of these miRNAs > 0.4 in absolute value. We developed a model using the 43 samples from 12 patients (chosen randomly from the 17 available observations) with severe disease, using Wilks' Lambda with Rao's F-approximation to assess the significance of correlations. We used the data for the remaining 5 severe patients as an independent test set.

For all correlations between individual variables, Spearman's rank-order correlation was used.

Where there were missing values for clinical data, correlations were calculated over all patients for which values were available. Heatmaps showing multiple correlations were J o u r n a l P r e -p r o o f created using the package 'pheatmap' in R with thresholds for the minimum absolute values shown chosen to limit the size of the heatmaps.

Stepwise logistic regression was also carried out in R using the functions 'glm' and 'step'.

Models were compared using the residual deviance (the lack of fit to the perfect model) and

Akaike Information Criterion (AIC; a measure of the prediction error). Classification results were obtained by assigning the class with the highest probability.

Predicting effective microRNA target sites in mammalian mRNAs

AGO CLIP-based imputation of potent siRNA sequences targeting SARS-CoV-2 with antifibrotic miRNA-like activity

Cytokine storm in the pathophysiology of COVID-19: Possible functional disturbances of miRNAs

Longitudinal analysis reveals that delayed bystander CD8+ T cell activation and early immune pathology distinguish severe COVID-19 from mild disease

Host-Viral Infection Maps Reveal Signatures of Severe COVID-19 Patients

Clinical and immunological features of severe and moderate coronavirus disease 2019

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Serum Cytokine and Chemokine Profile in Relation to the Severity of Coronavirus Disease 2019 in China

COVID-19 severity correlates with airway epithelium-immune cell interactions identified by single-cell analysis

A blood atlas of COVID-19 defines hallmarks of disease severity and specificity

Interleukin-6, CXCL10 and Infiltrating Macrophages in COVID-19-Related Cytokine Storm: Not One for All But All for One! Front Immunol 12

Circulating microRNAs as Potential Biomarkers of Infectious Disease

MicroRNA as Potential Biomarkers of Platelet Function on Antiplatelet Therapy: A Review

Therapeutic Significance of microRNA-Mediated Regulation of PARP-1 in SARS-CoV-2 Infection

Circulating microRNA profiles predict the severity of COVID-19 in hospitalized patients

An inflammatory cytokine signature predicts COVID-19 severity and survival

Reduction and Functional Exhaustion of T Cells in Patients With Coronavirus Disease 2019 (COVID-19). Front Immunol 11

The MiR-320 Family Is Strongly Downregulated in Patients with COVID-19 Induced Severe Respiratory Failure

Risk factors for death from COVID-19

Altered microRNA expression in COVID-19 patients enables identification of SARS-CoV-2 infection

The Diagnostic Significance of Usual Biochemical Parameters in Coronavirus Disease 19 (COVID-19): Albumin to Globulin Ratio and CRP to Albumin Ratio

IL-6 serum levels predict severity and response to tocilizumab in COVID-19: An observational study

Role of endothelial miR-24 in COVID-19 cerebrovascular events

Exosomal microRNAs Drive Thrombosis in COVID-19

In vitro induction of interleukin-8 by SARS-CoV-2 Spike protein is inhibited in bronchial epithelial IB3-1 cells by a miR-93-5p agomiR

Lymphopenia and neutrophilia at admission predicts severity and mortality in patients with COVID-19: a meta-analysis

Epithelial response to IFN-gamma promotes SARS-CoV-2 infection

miRNA-223 as a regulator of inflammation and NLRP3 inflammasome, the main fragments in the puzzle of immunopathogenesis of different inflammatory diseases and COVID-19

Uncontrolled Innate and Impaired Adaptive Immune Responses in Patients with COVID-19 Acute Respiratory Distress Syndrome

MicroRNA Mimics or Inhibitors as Antiviral Therapeutic Approaches Against COVID-19

CCR6 as a mediator of immunity in the lung and gut

Synergism of TNF-alpha and IFNgamma Triggers Inflammatory Cell Death, Tissue Damage, and Mortality in SARS-CoV-2 Infection and Cytokine Shock Syndromes

MicroRNAs based regulation of cytokine regulating immune expressed genes and their transcription factors in COVID-19

A dynamic COVID-19 immune signature includes associations with poor prognosis

TAM 2.0: tool for MicroRNA set analysis

Dynamic changes in serum IL-6, IL-8, and IL-10 predict the outcome of ICU patients with severe COVID-19

Association between IL-6 and severe disease and mortality in COVID-19 disease: a systematic review and meta-analysis

Longitudinal analyses reveal immunological misfiring in severe COVID-19

Longitudinal immune profiling reveals key myeloid signatures associated with COVID-19

miR-98 Regulates TMPRSS2 Expression in Human Endothelial Cells: Key Implications for COVID-19

Role of miR-2392 in driving SARS-CoV-2 infection

Platelets confound the measurement of extracellular miRNA in archived plasma

Targets the Transmembrane Glycoprotein Neuropilin-1 in Human Brain Microvascular Endothelial Cells. Noncoding

Interferon gamma, TGF-beta1 and RANTES expression in upper airway samples from SARS-CoV-2 infected patients

Potential Therapeutic Effect of Micrornas in Extracellular Vesicles from Mesenchymal Stem Cells against SARS-CoV-2

Non-Coding RNAs in COVID-19: Emerging Insights and Current Questions

Platelet aggregates, a marker of severe COVID-19 disease

Decreased serum levels of the inflammaging marker miR-146a are associated with clinical non-response to tocilizumab in COVID-19 patients

T cells and CCL20 in the lungs of mechanically ventilated patients with COVID-19

Severe COVID-19 Is Marked by a Dysregulated Myeloid Cell Compartment

Bias in observational study designs: prospective cohort studies

Circulating microRNAs as biomarkers in cancer diagnosis

Single-cell multi-omics analysis of the immune response in COVID-19

Hematological findings and complications of COVID-19

Inflammatory profiles across the spectrum of disease reveal a distinct role for GM-CSF in severe COVID-19

Monocyte-driven atypical cytokine storm and aberrant neutrophil activation as key mediators of COVID-19 disease severity

Diverse functional autoantibodies in patients with COVID-19

Serum microRNAs as powerful cancer biomarkers

Extrafollicular B cell responses correlate with neutralizing antibodies and morbidity in COVID-19

A new coronavirus associated with human respiratory disease in China

The study was funded by the UKRI MRC/NIHR award to the UK Coronavirus Immunology