Deconvoluting complex correlates of COVID-19 severity with a multi-omic pandemic tracking strategy

COVID-19 severity, rather than sex or age, predicts SARS-CoV-2 kinetics, and SARS-CoV-2 viral load from lower respiratory tract specimens may predict severe disease days before clinical deterioration for COVID-19 patients.

Metabolic biomarkers measured from single blood test can identify apparently healthy people at high susceptibility for developing severe pneumonia, and may also be useful for preventive COVID-19 screening.

Viral dynamics of SARS-CoV-2 in saliva from infected patients

Article Figures and data Abstract eLife digest Introduction Results Discussion Materials and methods Data availability References Decision letter Author response Article and author information Metrics Abstract End-stage kidney disease (ESKD) patients are at high risk of severe COVID-19. We measured 436 circulating proteins in serial blood samples from hospitalised and non-hospitalised ESKD patients with COVID-19 (n = 256 samples from 55 patients). Comparison to 51 non-infected patients revealed 221 differentially expressed proteins, with consistent results in a separate subcohort of 46 COVID-19 patients. Two hundred and three proteins were associated with clinical severity, including IL6, markers of monocyte recruitment (e.g. CCL2, CCL7), neutrophil activation (e.g. proteinase-3), and epithelial injury (e.g. KRT19). Machine-learning identified predictors of severity including IL18BP, CTSD, GDF15, and KRT19. Survival analysis with joint models revealed 69 predictors of death. Longitudinal modelling with linear mixed models uncovered 32 proteins displaying different temporal profiles in severe versus non-severe disease, including integrins and adhesion molecules. These data implicate epithelial damage, innate immune activation, and leucocyte–endothelial interactions in the pathology of severe COVID-19 and provide a resource for identifying drug targets. eLife digest COVID-19 varies from a mild illness in some people to fatal disease in others. Patients with severe disease tend to be older and have underlying medical problems. People with kidney failure have a particularly high risk of developing severe or fatal COVID-19. Patients with severe COVID-19 have high levels of inflammation, causing damage to tissues around the body. Many drugs that target inflammation have already been developed for other diseases. Therefore, to repurpose existing drugs or design new treatments, it is important to determine which proteins drive inflammation in COVID-19. Here, Gisby, Clarke, Medjeral-Thomas et al. measured 436 proteins in the blood of patients with kidney failure and compared the levels between patients who had COVID-19 to those who did not. This revealed that patients with COVID-19 had increased levels of hundreds of proteins involved in inflammation and tissue injury. Using a combination of statistical and machine learning analyses, Gisby et al. probed the data for proteins that might predict a more severe disease progression. In total, over 200 proteins were linked to disease severity, and 69 with increased risk of death. Tracking how levels of blood proteins changed over time revealed further differences between mild and severe disease. Comparing this data with a similar study of COVID-19 in people without kidney failure showed many similarities. This suggests that the findings may apply to COVID-19 patients more generally. Identifying the proteins that are a cause of severe COVID-19 – rather than just correlated with it – is an important next step that could help to select new drugs for severe COVID-19. Introduction Coronavirus disease (COVID-19), caused by the SARS-CoV-2 virus, displays wide clinical heterogeneity from asymptomatic to fatal disease. Patients with severe disease exhibit marked inflammatory responses and immunopathology. The mechanisms underlying this remain incompletely characterised, and the key molecular mediators are yet to be determined. The first treatment shown to reduce mortality from COVID-19 in randomised trials was dexamethasone (Horby et al., 2020), a corticosteroid that has broad non-specific effects on the immune system. Even with corticosteroid treatment, mortality in severe COVID-19 remains significant. There is a wide armamentarium of existing drugs that target inflammation more selectively, providing potential repurposing opportunities for the treatment of COVID-19. Recently, the REMAP-CAP trial has demonstrated efficacy of anti-IL6 receptor blockade in patients admitted to intensive care units with severe disease (Gordon et al., 2021). In order to select the most promising agents for future trials, we urgently need to better understand the molecular drivers of severe disease. Proteins are the effector molecules of biology and the targets of most drugs. Therefore, proteomic profiling to identify the key mediators of severe disease provides a valuable tool for identifying and prioritising potential drug targets (Suhre et al., 2021). Risk factors for severe or fatal COVID-19 include age, male sex, non-European ancestry, obesity, diabetes mellitus, cardiovascular disease, and immunosuppression (Williamson et al., 2020). End-stage kidney disease (ESKD) is one of the strongest risk factors for severe COVID-19 (estimated hazard ratio for death 3.69) (Williamson et al., 2020), and ESKD patients hospitalised with COVID-19 have a mortality of approximately 30% (Docherty et al., 2020; Corbett et al., 2020; Ng et al., 2020; Valeri et al., 2020). ESKD patients have a high prevalence of vascular and cardiometabolic disease (e.g. hypertension, ischaemic heart disease, diabetes), either as a result of the underlying cause of their renal disease or as a consequence of renal failure. In addition, ESKD results in both relative immunosuppression and chronic low-grade inflammation, which may impact viral defence and the host inflammatory response. Here we performed proteomic profiling of serial blood samples of ESKD patients with COVID-19, leveraging the unique opportunity for longitudinal sampling in both the outpatient and inpatient settings afforded by a large multi-ethnic haemodialysis cohort (Figure 1a). These data revealed 221 proteins that are dysregulated in COVID-19 versus matched non-infected ESKD patients. Using linear mixed models, joint models, and machine learning, we identified proteins that are markers of COVID-19 severity and risk of death. Finally, we characterised the temporal dynamics of the blood proteomic response during COVID-19 infection in ESKD patients, uncovering 32 proteins that display altered trajectories in patients with severe versus non-severe disease. Figure 1 with 1 supplement see all Download asset Open asset Study design. (a) Schematic representing a summary of the patient cohorts, sampling, and the major analyses. Blue and red stick figures represent outpatients and hospitalised patients, respectively. (b) Timing of serial blood sampling in relation to clinical course of COVID-19 (subcohort A). Black asterisks indicate when samples were obtained. Three patients were already in hospital prior to COVID-19 diagnosis (indicated by red bars). Results We recruited 55 ESKD patients with COVID-19 (subcohort A; Table 1). All patients were receiving haemodialysis prior to acquiring COVID-19. Blood samples were taken as soon as feasible following COVID-19 diagnosis. At time of initial sample, 30 patients were outpatients attending haemodialysis sessions and 25 were hospitalised inpatients (see Materials and methods, Figure 1). Following the initial blood sample, serial sampling was performed for 51/55 patients. We also recruited 51 non-infected haemodialysis patients as ESKD controls, mirroring the age, sex, and ethnicity distribution of the COVID-19 cases (Figure 1—figure supplement 1a–c). We used the Olink proteomics platform to measure 436 proteins (Supplementary file 1a) in 256 plasma samples from the COVID-19 patients and the 51 control samples. The proteins measured consisted of five multiplex 'panels' focussed on proteins relevant to immuno-inflammation, cardiovascular, and cardiometabolic disease. The 436 proteins assayed showed strong enrichment for immune-related proteins (Supplementary file 1b). Table 1 Characteristics of subcohort A. COVID-19-positive ESKD patients (n = 55)ESKD controls (n = 51) OverallPeak severity mild or moderate (n = 28)Peak severity severe or critical (n = 27)Age Median (IQR)72.2 62.5–77.373.4 65.5–76.468.5 61.8–78.870.1 62.2–75.1Sex M F39 (70.9%) 16 (29.1%)18 (64.3%) 10 (35.7%)21 (77.8%) 6 (22.2%)36 (70.6%) 15 (29.4%)Ethnicity White Black South Asian Asian (other) Other16 (29.1%) 8 (14.5%) 18 (32.7%) 4 (7.3%) 9 (16.4%)5 (17.9%) 5 (17.9%) 10 (35.7%) 1 (3.6%) 7 (25.0%)11 (40.7%) 3 (11.1%) 8 (29.6%) 3 (11.1%) 2 (7.4%)13 (25.5%) 8 (15.7%) 20 (39.2%) 3 (5.9%) 7 (13.7%)Diabetes34 (61.8%)*16 (57.1%)18 (66.7%)24 (47.1%)*Current smoker1 (1.8%)1 (3.6%)00ESKD cause DN Genetic GN HTN/vascular Other Unknown29 (52.7%) 1 (1.8%) 3 (5.5%) 5 (9.1%) 8 (14.5%) 9 (16.4%)14 (50.0%) 1 (3.6%) 1 (3.6%) 3 (10.7%) 5 (17.9%) 4 (14.3%)15 (55.6%) 0 2 (7.4%) 2 (7.4%) 3 (11.1%) 5 (18.5%)20 (39.2%) 1 (2.0%) 9 (17.6%) 7 (13.7%) 4 (7.8%) 10 (19.6%)Hospitalisation due to COVID-19†33 (60%)6 (21.4%)27 (100%)N/AFatal COVID-199 (16.3%)0 (0%)9 (33.3%)N/A DN = diabetic nephropathy. GN = glomerulonephritis. HTN = hypertension. IQR = inter-quartile range. 'South Asian' represents individuals with Indian, Pakistani, or Bangladeshi ancestry. Subsets defined according to peak WHO severity over the course of the illness. N/A = not applicable. *One patient had type 1 diabetes, the remainder type 2. †3 patients were hospitalised prior to COVID-19 diagnosis. 8 patients diagnosed with COVID-19 as outpatients subsequently deteriorated were hospitalised. In addition, we performed the Olink proteomic assays in 52 serum samples from a separate set of 46 COVID-19-positive ESKD patients (subcohort B) and 11 serum samples from ESKD COVID-19-negative controls (a subset of the controls described above). For the large majority of patients in subcohort B, only a single timepoint was available. A higher proportion of these patients (41/46, 89%) were hospitalised and had severe disease (Table 2) than in subcohort A (Figure 1, Table 1). Table 2 Characteristics of subcohort B. COVID-19-positive ESKD patients (n = 46)COVID-19-negative ESKD controls (n = 11)*Age Median (IQR)64.3 60.3–73.071.6 (61.7–73.9)Sex M F32 (69.6%) 14 (30.4%)8 (72.3%) 3 (27.3%)Ethnicity White Black South Asian Asian (other) Other11 (23.9%) 8 (17.4%) 12 (26.1%) 7 (15.2%) 8 (17.4%)3 (27.3%) 3 (27.3%) 3 (27.3%) 0 2 (18.2%)Diabetes29 (63.0%)6 (54.5%)Current smoker2 (4.3%)0 (%)ESKD cause DN Genetic GN HTN/vascular Other Unknown19 (41.3%) 1 (2.2%) 7 (15.2%) 3 (6.5%) 3 (6.5%) 13 (28.3%)5 (45.5%) 0 1 (9.1%) 1 (9.1%) 2 (18.2%) 2 (18.2%)Hospitalisation due to COVID-1941 (89.1%)N/ASevere or critical COVID-1933 (71.7%)N/AFatal COVID-199 (19.6%)N/A DN = diabetic nephropathy. GN = glomerulonephritis. HTN = hypertension. IQR = inter-quartile range. 'South Asian' represents individuals with Indian, Pakistani, or Bangladeshi ancestry. Subsets defined according to peak WHO severity over the course of the illness. N/A = not applicable. *These 11 controls are a subset of the control patients used in subcohort A. Proteomic differences between COVID-19-positive and -negative ESKD patients Principal component analysis (PCA) of proteomic data from subcohort A demonstrated differences between samples from COVID-19-positive cases and controls, although the two groups did not separate into discrete clusters (Figure 2a,b). To examine the effects of COVID-19 on the plasma proteome, we performed a differential expression analysis in subcohort A between COVID-19 cases (n = 256 samples passing quality control [QC] from 55 patients) and non-infected ESKD controls (n = 51) using linear mixed models, which account for serial samples from the same individual (see Materials and methods). This revealed 221 proteins associated with COVID-19 (5% false discovery rate, FDR); the vast majority were upregulated, with only 40 downregulated (Figure 3a, Supplementary file 1c). In order to provide a succinct and standardised nomenclature, we report proteins by the symbols of the genes encoding them (see Supplementary file 1a for a mapping of symbols to full protein names). The most strongly upregulated proteins (in terms of fold change) were DDX58, CCL7, IL6, CXCL11, KRT19, and CXCL10, and the most strongly downregulated were SERPINA5, CCL16, FABP2, PON3, ITGA11, and MMP12 (Figure 3—figure supplement 1). Notably, many of the upregulated proteins were chemotaxins. Figure 2 with 2 supplements see all Download asset Open asset Principal component analysis. PC = principal component. Each point represents a sample. Colouring indicates COVID-19 status. The directions and relative sizes of the six largest PC loadings are plotted as arrows (middle column). (a, b) Subcohort A. Due to serial sampling, there are multiple samples for most patients. The proportion of variance explained in subcohort A by each PC is shown in parentheses on the axis labels. (c, d) Subcohort B. Samples are projected into the PCA coordinates from subcohort A. Figure 3 with 4 supplements see all Download asset Open asset Identification of dysregulated proteins. (a) Proteins upregulated (red) or downregulated (blue) in COVID-19-positive patients versus COVID-19-negative ESKD patients n = 256 plasma samples from 55 COVID-19-positive patients, versus n = 51 ESKD controls (one sample per control patient). (b) Proteins associated with disease severity associations of protein levels against WHO severity score at the time of sampling. Linear gradient indicates the effect size. A positive effect size (red) indicates that an increase in protein level is associated with increasing disease severity and a negative gradient (blue) the opposite. n = 256 plasma samples from 55 COVID-19-positive patients. For (a, b), p-values from linear mixed models after Benjamini–Hochberg adjustment; significance threshold = 5% FDR; dark-grey = non-significant. (c) Heatmap showing protein levels for selected proteins with strong associations with severity. Each column represents a sample (n = 256 COVID-19 samples and 51 non-infected samples). Each row represents a protein. Proteins are annotated using the symbol of their encoding gene. For the purposes of legibility, not all significantly associated proteins are shown; the heatmap is limited to the 17% most up- or downregulated proteins (by effect size) of those with a significant association. Proteins are ordered by hierarchical clustering. Samples are ordered by WHO severity at the time of blood sample ('Severity'). 'Overall course' indicates the peak WHO severity over the course of the illness. We observed that a high proportion of the measured proteins were associated with COVID-19. Given the highly targeted nature of the Olink panels that we used (enriched for immune and inflammation-related proteins), this was not surprising. Nevertheless, to ensure that the Benjamini–Hochberg adjustment of p-values was controlling the FDR at the 5% level, we performed two additional analyses (see Materials and methods). First, we estimated the FDR using an alternative method (the plug-in procedure ; Hastie et al., 2001); this confirmed appropriate FDR control. Second, we used permutation to estimate the distribution of the number of proteins expected to be declared significant under the null hypothesis of no association between any proteins and COVID-19. This showed that the probability of observing the number of differentially abundant proteins we identified was highly unlikely under the null (empirical p<1×10−5; Figure 3—figure supplement 2). Although our COVID-19-negative controls were well matched in terms of age, sex, and ethnicity (Figure 1—figure supplement 1a–c), perfect matching of comorbidities was not feasible in the context of the healthcare emergency at the time of patient recruitment. There was a higher prevalence of diabetes in the COVID-19 cases compared to the controls (61.8% versus 47.1%, respectively; Table 1). To evaluate whether differing rates of diabetes had impacted the proteins identified as differentially abundant between cases and controls, we performed a sensitivity analysis adding diabetes as an additional covariate in the linear mixed model. This did not materially affect our findings; estimated effect sizes and –log10 p-values from models with and without the inclusion of diabetes were highly correlated (Pearson r > 0.99, and r = 0.95, respectively; Figure 3—figure supplement 3a,b). Full results from both models are shown in Supplementary file 1c. Similarly, there were also differences in the underlying cause of ESKD in cases compared to controls (Table 1). We therefore performed a further sensitivity analysis adjusting for underlying cause of renal failure. This did not make any meaningful difference to our results (Figure 3—figure supplement 3c,d, Supplementary file 1c). We also considered the possibility that timing of haemodialysis might affect the plasma proteome. To minimise the impact of this, all samples were taken prior to haemodialysis. For the large majority (86.6%) of samples, the most recent haemodialysis was between 48 and 72 hr prior to blood draw. This consistency in timing of blood sampling reduces the potential for impact of this issue. Nevertheless, to evaluate whether timing of haemodialysis might have impacted our results, we performed a sensitivity analysis including time from last haemodialysis as a covariate. Our results were not materially affected by this, with −log10 p-values and estimated effect sizes very highly correlated with those obtained without inclusion of this covariate (Pearson r > 0.99 for effect size estimates and for −log10 p-values; Figure 3—figure supplement 4a,b, Supplementary file 1c). We used the smaller subcohort B (n = 52 serum samples from 46 patients with COVID-19; see Materials and methods) for validation. We first projected the data from subcohort B into the PCA space of subcohort A to examine how well the separation of cases and controls in the PCA space replicated (see Materials and methods). This revealed clearer separation of infected and non-infected patients than in subcohort A (Figure 2c,d), perhaps reflecting the higher proportion of hospitalised patients (41 of 46 patients) in subcohort B (Table 2). We next performed differential abundance analysis in subcohort B and found 201 proteins that were dysregulated in cases versus controls (5% FDR) (Supplementary file 1c). Of the 221 differentially abundant proteins from subcohort A, 150 (69.7%) were also identified in subcohort B at 5% FDR (Figure 4a). Effect sizes in each dataset showed a strong correlation (r = 0.80, Figure 4b). This demonstrates that our findings are highly reproducible despite differences in sample sizes and blood materials (plasma versus serum in subcohorts A and B, respectively). Figure 4 with 1 supplement see all Download asset Open asset Validation. (a) Overlap between the significant associations in the differential abundance analysis between ESKD patients with and without COVID-19 in subcohorts A and B. 5% FDR was used as the significance threshold in both analyses. (b) Comparison of estimated effect sizes for all 436 proteins in the differential abundance analyses (COVID-19 positive versus negative) in subcohort A and B. Each point represents a protein. Pearson's r is shown. Differential abundance analyses were performed using linear mixed models. Subcohort A analysis (plasma samples): 256 samples from 55 COVID-19 patients versus 51 non-infected patient samples (single time-point). Subcohort B (serum samples): 52 samples from 55 COVID-19 patients and 11 non-infected patient samples (single timepoint). Proteins associated with COVID-19 severity Examination of the principal components plot labelling samples by clinical severity at the time of sampling (defined by WHO severity scores, graded as mild, moderate, severe, or critical) demonstrated a gradient of COVID-19 severity, best captured by principal components 1 and 3 (Figure 2—figure supplement 1a). To determine the proteomic effects of COVID-19 severity, we tested for associations between proteins and WHO severity score at the time of blood sampling, using linear mixed models with severity encoded as an ordinal predictor (see Materials and methods). This analysis revealed 203 proteins associated with severity (Figure 3b, Supplementary file 1d). The majority of these were upregulated in more severe disease, with only 42 downregulated. A sensitivity analysis adjusting for time since last haemodialysis made no significant impact on our results (Figure 3—figure supplement 4c,d, Supplementary file 1d). Consistent with previous reports, we found that severe COVID-19 was characterised by elevated IL6. In addition, we observed a signature of upregulated monocyte chemokines (e.g. CCL2, CCL7, CXCL10), neutrophil activation and degranulation (e.g. PRTN3, MPO), and epithelial injury (e.g. KRT19, AREG, PSIP1, GRN). (Figures 3b,c and 5). SERPINA5 and leptin showed the greatest downregulation as COVID-19 severity increased (Figure 3b,c). Figure 5 Download asset Open asset proteins strongly associated with COVID-19 severity. showing distribution of plasma protein levels according to COVID-19 at the time of blood draw. indicate and inter-quartile range. n = 256 samples from 55 COVID-19 patients and 51 samples from non-infected patients. WHO severity indicates the clinical severity score of the patient at the time the sample was n = moderate n = severe n = critical n = 15 samples. monocyte markers of epithelial injury. two neutrophil and IL6. We next how the COVID-19 severity protein signature to the proteins that are differentially abundant between cases and The majority of proteins were also identified as differentially abundant in the COVID-19-positive versus -negative analysis (Figure fold for proteins in COVID-19 versus non-infected patients were correlated with effect sizes in the severity that the proteins most upregulated in cases versus controls also to the greatest in severe disease (Figure there were some (e.g. that were strongly associated with severity, not differentially expressed in infected versus non-infected patients (Figure Figure 6 Download asset Open asset Comparison of proteins differentially expressed in COVID-19 with those associated with clinical severity. (a) Overlap between the proteins significantly differentially expressed in COVID-19 (n = 256 COVID-19 samples and 51 non-infected versus those associated with severity n = 256 (subcohort A). 5% FDR was used as the significant in both analyses. (b) Comparison of effect sizes for each protein in the COVID-19-positive versus -negative analysis and severity analysis Each point represents a protein. Pearson's r is shown. (c) of proteins associated with severity, not significantly differentially abundant in the of all cases versus showing distribution of plasma protein levels according to COVID-19 at the time of blood draw. indicate and inter-quartile range. n = 256 samples from 55 COVID-19 patients and 51 samples from non-infected patients. WHO severity indicates the clinical severity score of the patient at the time the sample was n = moderate n = severe n = critical n = 15 samples. learning to predict COVID-19 severity PCA revealed that some samples from patients who had mild or moderate disease at the time of sampling with samples from patients with severe disease (Figure 2—figure supplement 1a). Examination of the same PCA plot labelling samples according to the clinical course by peak WHO severity score over the of the (Figure 2—figure supplement revealed that these samples from individuals who subsequently developed severe or critical disease. This that molecular may clinical To evaluate this we used learning to whether the proteomic signature of the first blood sample for each patient in our dataset could identify whether the patient either had severe COVID-19 at the time of sampling or severe disease in the differential expression analyses each protein of all proteins in the Using we a on the first sample for each COVID-19 patient to predict the clinical defined by peak WHO severity. For the purposes of this we clinical course into either WHO or The method in peak severity. using only predictors and clinical the method in peak severity. clinical proteins did not compared to using proteomic predictors that the information in the clinical predictors is captured at the proteomic we not that proteomic profiling is to clinical for risk during this the selected by the proteins of We therefore the to identify key proteins by (see Materials and methods, Supplementary file The most important proteins for the of or future severe disease were IL18BP, GDF15, KRT19, and (Figure is that this as a key of severe disease. Figure 7 with 1 supplement see all Download asset Open asset of severe COVID-19 and death. (a) The 12 most important proteins for clinical course (defined by peak COVID-19 WHO using a is important for it is to in many of and be to the have a The all was used as the (b) Proteins that are significant predictors of death n = 256 samples from 55 COVID-19-positive patients, of Risk estimates are from a joint model. indicate For proteins with a positive risk a higher to a high risk of and for proteins with negative Proteins associated with risk of death of 55 patients in subcohort A We therefore to identify proteins associated with risk of death. To the nature of protein for of we joint models, which linear mixed models and models et al., (see Materials and methods). This analysis identified proteins for which increased was associated with increased risk of death (Figure Supplementary file including and and 25 proteins for which increased was associated with risk of including and with clinical A number of clinical have associations with COVID-19 (e.g. elevated inflammatory and et al., 2020). We therefore compared our proteomic data from COVID-19 patients at each timepoint to clinical using linear mixed models (see Materials and methods). We found associations between plasma proteins and all clinical (Figure Supplementary file Many of these proteins were also markers of severity (e.g. IL6, KRT19, and were strongly associated with and and Of CCL7, a monocyte that was also identified as an important of severity by the was associated with monocyte and inflammatory neutrophil was associated with which IL6, and and with the and Figure 8 Download asset Open asset of clinical markers with plasma proteins. Proteins that are (red) or (blue) associated with clinical (5% p-values from differential abundance analysis using linear mixed models after Benjamini–Hochberg = non-significant. Two associations were found for shown – see Supplementary file Longitudinal analysis proteins with temporal profiles according to severity The immune response to infection is and therefore provide only the serial sampling in our dataset (Figure we the temporal of each protein and whether or not any protein trajectories in patients with a versus clinical This was a

End-stage kidney disease (ESKD) patients are at high risk of severe COVID-19. We measured 436 circulating proteins in serial blood samples from hospitalised and non-hospitalised ESKD patients with COVID-19 (n = 256 samples from 55 patients). Comparison to 51 non-infected patients revealed 221 differentially expressed proteins, with consistent results in a separate subcohort of 46 COVID-19 patients. Two hundred and three proteins were associated with clinical severity, including IL6, markers of monocyte recruitment (e.g. CCL2, CCL7), neutrophil activation (e.g. proteinase-3), and epithelial injury (e.g. KRT19). Machine-learning identified predictors of severity including IL18BP, CTSD, GDF15, and KRT19. Survival analysis with joint models revealed 69 predictors of death. Longitudinal modelling with linear mixed models uncovered 32 proteins displaying different temporal profiles in severe versus non-severe disease, including integrins and adhesion molecules. These data implicate epithelial damage, innate immune activation, and leucocyte–endothelial interactions in the pathology of severe COVID-19 and provide a resource for identifying drug targets.

COVID-19 in solid organ transplant recipients: Initial report from the US epicenter.

Metagenomic next-generation sequencing (mNGS) offers an agnostic approach for emerging pathogen detection directly from clinical specimens. In contrast to targeted methods, mNGS also provides valuable information on the composition of the microbiome and might uncover coinfections that may associate with disease progression and impact prognosis. To evaluate the use of mNGS for detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and/or other infecting pathogens, we applied direct Oxford Nanopore long-read third-generation metatranscriptomic and metagenomic sequencing. Nasopharyngeal (NP) swab specimens from 50 patients under investigation for CoV disease 2019 (COVID-19) were sequenced, and the data were analyzed by the CosmosID bioinformatics platform. Further, we characterized coinfections and the microbiome associated with a four-point severity index. SARS-CoV-2 was identified in 77.5% (31/40) of samples positive by RT-PCR, correlating with lower cycle threshold (Ct) values and fewer days from symptom onset. At the time of sampling, possible bacterial or viral coinfections were detected in 12.5% of SARS-CoV-2-positive specimens. A decrease in microbial diversity was observed among COVID-19-confirmed patients (Shannon diversity index, P = 0.0082; Chao richness estimate, P = 0.0097; Simpson diversity index, P = 0.018), and differences in microbial communities were linked to disease severity (P = 0.022). Furthermore, statistically significant shifts in the microbiome were identified among SARS-CoV-2-positive and -negative patients, in the latter of whom a higher abundance of Propionibacteriaceae (P = 0.028) and a reduction in the abundance of Corynebacterium accolens (P = 0.025) were observed. Our study corroborates the growing evidence that increased SARS-CoV-2 RNA detection from NP swabs is associated with the early stages rather than the severity of COVID-19. Further, we demonstrate that SARS-CoV-2 causes a significant change in the respiratory microbiome. This work illustrates the utility of mNGS for the detection of SARS-CoV-2, for diagnosing coinfections without viral target enrichment or amplification, and for the analysis of the respiratory microbiome.IMPORTANCE SARS-CoV-2 has presented a rapidly accelerating global public health crisis. The ability to detect and analyze viral RNA from minimally invasive patient specimens is critical to the public health response. Metagenomic next-generation sequencing (mNGS) offers an opportunity to detect SARS-CoV-2 from nasopharyngeal (NP) swabs. This approach also provides information on the composition of the respiratory microbiome and its relationship to coinfections or the presence of other organisms that may impact SARS-CoV-2 disease progression and prognosis. Here, using direct Oxford Nanopore long-read third-generation metatranscriptomic and metagenomic sequencing of NP swab specimens from 50 patients under investigation for COVID-19, we detected SARS-CoV-2 sequences by applying the CosmosID bioinformatics platform. Further, we characterized coinfections and detected a decrease in the diversity of the microbiomes in these patients. Statistically significant shifts in the microbiome were identified among COVID-19-positive and -negative patients, in the latter of whom a higher abundance of Propionibacteriaceae and a reduction in the abundance of Corynebacterium accolens were observed. Our study also corroborates the growing evidence that increased SARS-CoV-2 RNA detection from NP swabs is associated with the early stages of disease rather than with severity of disease. This work illustrates the utility of mNGS for the detection and analysis of SARS-CoV-2 from NP swabs without viral target enrichment or amplification and for the analysis of the respiratory microbiome.

Severe acute coronavirus disease 2019 (COVID-19) is uncommon in children; however, its development can lead to longer-term health problems. Understanding factors associated with clinical deterioration in paediatric patients is therefore of public health relevance. Early triage enables closer monitoring, tailored counselling, and timely escalation of care. International studies have associated obesity, chronic illness, and certain sociodemographic factors with worse outcomes. However, robust real-world datasets remain sparse, particularly in Germany, where statutory surveillance often fails to capture detailed clinical data. To address this gap, we analysed 731 polymerase-chain-reaction (PCR)-confirmed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in children aged ≤ 15 years reported to the Cologne and Augsburg health departments between January and December 2021. • Primary: Identify independent factors associated with a severe acute COVID-19 course in this population-based study with defined sampling framepaediatric cohort. • Secondary (exploratory): Estimate the prevalence of long COVID–compatible symptoms (persisting >4 weeks) and explore how their frequency changes across acute COVID-19 severity strata. This analysis is based on the CoCo-Fakt cross-sectional study, which captured 731 PCR-confirmed SARS-CoV-2 infections in children aged ≤15 years between January and December 2021. The parent questionnaire captured sociodemographic variables, body mass index (BMI), chronic illnesses, quarantine details, acute COVID-19 severity (asymptomatic to severe), and persisting symptoms (>4 weeks). We then used multivariable logistic regression to examine whether age, sex, socioeconomic status, migration background, BMI, and chronic illness were independently associated with a severe disease course. For exploratory comparisons between children with and without long COVID, we applied t-tests for continuous and Fisher’s exact or chi-square tests for categorical variables. Among the included participants, 67 (9.6%) experienced a severe disease course of COVID-19. In multivariable analysis, chronic illness emerged as the only independent factor associated with severe COVID-19, conferring an almost sixfold higher odds of severe COVID-19 (OR 5.90, 95% CI 2.98–11.68). BMI showed a positive trend but was not statistically significant. Exploratory analyses indicated associations with older age, chronic illness, and increasing acute disease severity. The most frequently reported symptoms were fatigue, sleep disturbances, and problems concentrating. Chronic illness was consistently associated with a higher likelihood of severe acute COVID-19 in children. Given the cross-sectional design, exploratory, parent-reported data suggest older age and more severe acute disease courses may increase the likelihood of persisting symptoms, which warrants confirmation in prospective cohorts. Therefore, children with more severe acute disease or chronic illness may benefit from tailored follow-up to better understand potential long-term impairment, including potential progression to myalgic encephalomyelitis/chronic fatigue syndrome.

Association of Rhesus factor blood type with risk of SARS-CoV-2 infection and COVID-19 severity.

TO THE EDITOR: We read with great interest the paper by Deng et al. ((1)). Obesity increases the risk of coronavirus disease 2019 (COVID-19) morbidity and mortality ((2, 3)). However, Deng et al. suggested that not simply obesity but also visceral adiposity is an independent risk factor for COVID-19 complications in young adults ((1)). Interestingly, visceral fat could serve as a reservoir for the virus and amplify the inflammatory response ((4, 5)). Deng et al. reported computed tomography (CT) data of ectopic fat depots, such as liver fat and epicardial adipose tissue (EAT), in young patients with COVID-19 ((1)). EAT, the visceral fat depot of the heart, has been suggested to play a role in COVID-19 myocardial inflammation ((6-8)). Hence, we retrospectively analyzed EAT from CT scans of patients who were admitted for COVID-19. We collected data from 41 patients with laboratory-confirmed COVID-19 infection who were admitted at the Policlinico San Donato, San Donato Milanese, Milan, University of Milan, Italy, between April 1 and April 9, 2020. A confirmed case of COVID-19 was defined by a positive result on a reverse transcriptase-polymerase chain reaction assay of a specimen collected on a nasopharyngeal swab. Chest CT scan was performed on admission day 1 in patients with suspected or confirmed COVID-19 infection to evaluate the presence of pulmonary embolism. EAT measurement was retrospectively obtained from each CT scan and analyzed according to the clinical and radiological criteria defining COVID-19 severity. EAT and subcutaneous adipose tissue (SAT) density was defined as mean attenuation expressed in Hounsfield units (HU). Patients' features on admission are reported in Table 1. More than half (54%) of the patients presented clinical and CT signs of pulmonary embolism. Almost two-thirds (26) of the patients had no coronary calcium content (CAC) score, and only two patients had a severe CAC score (> 400). Overall, EAT HU was significantly greater than SAT HU (−95 vs. −118 HU; P < 0.01); mean EAT thickness was 5.5 mm. We then compared EAT attenuation between the four different groups of patients according to COVID-19 severity. EAT attenuation significantly increased with increasing COVID-19 severity, whereas SAT attenuation did not substantially change, as depicted in Figure 1. Patients with severe and critical COVID-19 had significantly greater EAT attenuation than those presenting with mild and moderate COVID-19 (P ≤ 0.01 for all the comparisons; 95% CI: −99 to −69 HU). Results were substantially similar when EAT HU was calculated with or without contrast. CT-measured EAT thickness was similar among the groups of COVID-19 severity. EAT HU was significant correlated (r = −0.45; P < 0.05) with high-sensitivity troponin T levels, while there was no significant correlation with interleukin-6 levels (r = 0.05; P = 0.76). Moreover, EAT HU was significantly correlated with peripheral oxygen saturation (r = −054; P < 0.05) and body temperature (r = −0.43; P < 0.05). EAT attenuation reflects inflammatory changes within the fat depot ((9)). In our analysis, EAT showed imaging signs of increased inflammation in patients with severe and critical COVID-19. CT-measured EAT attenuation was consensually greater with higher degree of COVID-19 severity. Remarkably, EAT attenuation was similar to that observed in coronary artery disease despite most of these patients with COVID-19 having no history of coronary artery disease and no CAC. Deng et al. ((1)) found not quite statistically significant differences in CT-EAT attenuation between patients with moderate and severe COVID-19. However, in our study, patients were older and more critical than those evaluated by Deng et al. ((1)). Hence, the lower prevalence of patients with severe COVID-19 in the study by Deng et al. may explain the lack of statistically significant difference in CT-EAT attenuation. Also, elderly participants were more likely to have myocardial inflammation and injury ((1)). Notably, we found that SAT attenuation, unlike EAT, did not progress with the severity of COVID-19. This could be somehow consistent with Deng et al. who found no difference in SAT thickness between moderate and severe patients. Although our analysis presents with some limitations, cardiac imaging, particularly CT-measured EAT attenuation, could play a diagnostic and prognostic role in patients with COVID-19 with obesity ((10)). The authors declared no conflict of interest.

IntroductionThere is considerable variability in symptoms and severity of COVID-19 among patients infected by the SARS-CoV-2 virus. Linking host and virus genome sequence information to antibody response and biological information...

Serial simultaneously self-swabbed samples from multiple sites show similarly decreasing SARS-CoV-2 loads in COVID-19 cases of differing clinical severity

Predicting Severe Outcomes in COVID-19

Vascular mechanisms and manifestations of COVID-19

Early studies have reported various electrolyte abnormalities at admission in patients who progress to the severe form of coronavirus disease 2019 (COVID-19). As electrolyte imbalance may not only impact patient care, but provide insight into the pathophysiology of COVID-19, we aimed to analyse all early data reported on electrolytes in COVID-19 patients with and without severe form. An electronic search of Medline (PubMed interface), Scopus and Web of Science was performed for articles comparing electrolytes (sodium, potassium, chloride and calcium) between COVID-19 patients with and without severe disease. A pooled analysis was performed to estimate the weighted mean difference (WMD) with 95% confidence interval. Five studies with a total sample size of 1415 COVID-19 patients. Sodium was significantly lower in patients with severe COVID-19 (WMD: -0.91 mmol/L [95% CI: -1.33 to -0.50 mmol/L]). Similarly, potassium was also significantly lower in COVID-19 patients with severe disease (WMD: -0.12 mmol/L [95% CI: -0.18 to -0.07 mmol/L], I2=33%). For chloride, no statistical differences were observed between patients with severe and non-severe COVID-19 (WMD: 0.30 mmol/L [95% CI: -0.41 to 1.01 mmol/L]). For calcium, a statistically significant lower concentration was noted in patients with severe COVID-19 (WMD: -0.20 mmol/L [95% CI: -0.25 to -0.20 mmol/L]). This pooled analysis confirms that COVID-19 severity is associated with lower serum concentrations of sodium, potassium and calcium. We recommend electrolytes be measured at initial presentation and serially monitored during hospitalization in order to establish timely and appropriate corrective actions.