Prognostic characteristics inform risk stratification in intensive care unit (ICU) patients with coronavirus disease 2019 (COVID-19). We obtained blood samples (n = 474) from hospitalized COVID-19 patients (n = 123), non-COVID-19 ICU sepsis patients (n = 25) and healthy controls (n = 30). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was detected in plasma or serum (RNAemia) of COVID-19 ICU patients when neutralizing antibody response was low. RNAemia is associated with higher 28-day ICU mortality (hazard ratio [HR], 1.84 [95% CI, 1.22–2.77] adjusted for age and sex). RNAemia is comparable in performance to the best protein predictors. Mannose binding lectin 2 and pentraxin-3 (PTX3), two activators of the complement pathway of the innate immune system, are positively associated with mortality. Machine learning identified ‘Age, RNAemia’ and ‘Age, PTX3’ as the best binary signatures associated with 28-day ICU mortality. In longitudinal comparisons, COVID-19 ICU patients have a distinct proteomic trajectory associated with mortality, with recovery of many liver-derived proteins indicating survival. Finally, proteins of the complement system and galectin-3-binding protein (LGALS3BP) are identified as interaction partners of SARS-CoV-2 spike glycoprotein. LGALS3BP overexpression inhibits spike-pseudoparticle uptake and spike-induced cell-cell fusion in vitro.

Reverse transcription (RT) of RNA templates containing RNA modifications leads to synthesis of cDNA containing information on the modification in the form of misincorporation, arrest, or nucleotide skipping events. A compilation of such events from multiple cDNAs represents an RT-signature that is typical for a given modification, but, as we show here, depends also on the reverse transcriptase enzyme. A comparison of 13 different enzymes revealed a range of RT-signatures, with individual enzymes exhibiting average arrest rates between 20 and 75%, as well as average misincorporation rates between 30 and 75% in the read-through cDNA. Using RT-signatures from individual enzymes to train a random forest model as a machine learning regimen for prediction of modifications, we found strongly variegated success rates for the prediction of methylated purines, as exemplified with N1-methyladenosine (m1A). Among the 13 enzymes, a correlation was found between read length, misincorporation, and prediction success. Inversely, low average read length was correlated to high arrest rate and lower prediction success. The three most successful polymerases were then applied to the characterization of RT-signatures of other methylated purines. Guanosines featuring methyl groups on the Watson-Crick face were identified with high confidence, but discrimination between m1G and m22G was only partially successful. In summary, the results suggest that, given sufficient coverage and a set of specifically optimized reaction conditions for reverse transcription, all RNA modifications that impede Watson-Crick bonds can be distinguished by their RT-signature.

Aims Coronavirus disease 2019 (COVID-19) can lead to multiorgan damage. MicroRNAs (miRNAs) in blood reflect cell activation and tissue injury. We aimed to determine the association of circulating miRNAs with COVID-19 severity and 28-day intensive care unit (ICU) mortality. Methods and results We performed RNA-Seq in plasma of healthy controls (n = 11), non-severe (n = 18) and severe (n = 18) COVID-19 patients and selected 14 miRNAs according to cell- and tissue origin for measurement by reverse transcription quantitative polymerase chain reaction (RT-qPCR) in a separate cohort of mild (n = 6), moderate (n = 39) and severe (n = 16) patients. Candidates were then measured by RT-qPCR in longitudinal samples of ICU COVID-19 patients (n = 240 samples from n = 65 patients). 60 miRNAs, including platelet-, endothelial-, hepatocyte- and cardiomyocyte-derived miRNAs, were differentially expressed depending on severity, with increased miR-133a and reduced miR-122 also being associated with 28-day mortality. We leveraged mass spectrometry-based proteomics data for corresponding protein trajectories. Myocyte-derived (myomiR) miR-133a was inversely associated with neutrophil counts and positively with proteins related to neutrophil degranulation, such as myeloperoxidase. In contrast, levels of hepatocyte-derived miR-122 correlated to liver parameters and to liver-derived positive (inverse association) and negative acute phase proteins (positive association). Finally, we compared miRNAs to established markers of COVID-19 severity and outcome, i.e. SARS-CoV-2 RNAemia, age, BMI, D-dimer and troponin. Whilst RNAemia, age and troponin were better predictors of mortality, miR-133a and miR-122 showed superior classification performance for severity. In binary and triplet combinations, miRNAs improved classification performance of established markers for severity and mortality. Conclusion Circulating miRNAs of different tissue origin, including several known cardiometabolic biomarkers, rise with COVID-19 severity. MyomiR miR-133a and liver-derived miR-122 also relate to 28-day mortality. MiR-133a reflects inflammation-induced myocyte damage, whilst miR-122 reflects the hepatic acute phase response. Translational perspective Adding biomarkers to conventional scores for illness severity and mortality could improve prognostic performance in COVID-19 patients. Circulating miRNAs are emerging as promising biomarkers with tissue specific origins but have only sparsely been investigated in COVID-19. We quantified circulating miRNAs of different tissue origin in COVID-19 patients, identifying several miRNAs of the cardiometabolic system to be associated with severity. Myocyte-derived miR-133a and liver-derived miR-122 also associated with mortality. Through longitudinal proteomics measurements, we related myomiR miR-133a release to neutrophil activation and miR-122 release to the hepatic acute phase response. Our findings highlight key pathophysiological changes and provide first evidence on the performance of miRNA biomarkers in COVID-19.