Proteomic analysis of circulating immune cells identifies novel cellular phenotypes associated with COVID-19 severity

Certain serum proteins, including CRP and D-dimer, have prognostic value in patients with SARS-CoV-2 infection. Nonetheless, these factors are non-specific, and provide limited mechanistic insight into the peripheral blood mononuclear cell (PBMC) populations which drive the pathogenesis of severe COVID-19. To identify novel cellular phenotypes associated with disease progression, we here describe a comprehensive, unbiased analysis of the total and plasma membrane proteomes of PBMCs from a cohort of 40 unvaccinated individuals with SARS-CoV-2 infection, spanning the whole spectrum of disease severity. Combined with RNA-seq and flow cytometry data from the same donors, we define a comprehensive multi-omic profile for each severity level, revealing cumulative immune cell dysregulation in progressive disease. In particular, the cell surface proteins CEACAMs1, 6 and 8, CD177, CD63 and CD89 are strongly associated with severe COVID-19, corresponding to the emergence of atypical CD3+CD4+CD177+ and CD16+CEACAM1/6/8+ mononuclear cells. Utilisation of these markers may facilitate real-time patient assessment by flow cytometry, and identify immune cell populations that could be targeted to ameliorate immunopathology.

Mapping interindividual dynamics of innate immune response at single-cell resolution

Common genetic variants modulate the cellular response to viruses and are implicated in a range of immune pathologies, including infectious and autoimmune diseases. The transcriptional antiviral response is known to vary between infected cells from a single individual, yet how genetic variants across individuals modulate the antiviral response (and its cell-to-cell variability) is not well understood. Here, we triggered the antiviral response in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-seq. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), the first statistical approach designed to identify dynamic eQTLs across a transcriptional trajectory of cell populations, without aggregating single-cell data into pseudo-bulk. This allows us to uncover the underlying architecture and variability of antiviral response across responding cells, and to identify more than two thousands eQTLs modulating the dynamic changes during this response. Many of these eQTLs colocalise with risk loci identified in GWAS of infectious and autoimmune diseases. As a case study, we focus on a COVID-19 susceptibility locus, colocalised with the antiviral OAS1 splicing QTL. We validated it in blood cells from a patient cohort and in the infected nasal cells of a patient with the risk allele, demonstrating the utility of GASPACHO to fine-map and functionally characterise a genetic locus. In summary, our novel analytical approach provides a new framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

The impact of hypoxia on B cells in COVID-19

Prominent early features of COVID-19 include severe, often clinically silent, hypoxia and a pronounced reduction in B cells, the latter important in defence against SARS-CoV-2. This brought to mind the phenotype of mice with VHL-deficient B cells, in which Hypoxia-Inducible Factors are constitutively active, suggesting hypoxia might drive B cell abnormalities in COVID-19. We demonstrated the breadth of early and persistent defects in B cell subsets in moderate/severe COVID-19, including reduced marginal zone-like, memory and transitional B cells, changes we also observed in B cell VHL-deficient mice. This was corroborated by hypoxia-related transcriptional changes in COVID-19 patients, and by similar B cell abnormalities in mice kept in hypoxic conditions, including reduced marginal zone and germinal center B cells. Thus hypoxia might contribute to B cell pathology in COVID-19, and in other hypoxic states. Through this mechanism it may impact on COVID-19 outcome, and be remediable through early oxygen therapy.

The local and systemic response to SARS-CoV-2 infection in children and adults

While a substantial proportion of adults infected with SARS-CoV-2 progress to develop severe disease, children rarely manifest respiratory complications. Therefore, understanding differences in the local and systemic response to SARS-CoV-2 infection between children and adults may provide important clues about the pathogenesis of SARS-CoV-2 infection. To address this, we first generated a healthy reference multi-omics single cell data set from children (n=30) in whom we have profiled triple matched samples: nasal and tracheal brushings and PBMCs, where we track the developmental changes for 42 airway and 31 blood cell populations from infancy, through childhood to adolescence. This has revealed the presence of naive B and T lymphocytes in neonates and infants with a unique gene expression signature bearing hallmarks of innate immunity. We then contrast the healthy reference with equivalent data from severe paediatric and adult COVID-19 patients (total n=27), from the same three types of samples: upper and lower airways and blood. We found striking differences: children with COVID-19 as opposed to adults had a higher proportion of innate lymphoid and non-clonally expanded naive T cells in peripheral blood, and a limited interferon-response signature. In the airway epithelium, we found the highest viral load in goblet and ciliated cells and describe a novel inflammatory epithelial cell population. These cells represent a transitional regenerative state between secretory and ciliated cells; they were found in healthy children and were enriched in paediatric and adult COVID-19 patients. Epithelial cells display an antiviral and neutrophil-recruiting gene signature that is weaker in severe paediatric versus adult COVID-19. Our matched blood and airway samples allowed us to study the spatial dynamics of infection. Lastly, we provide a user-friendly interface for this data1 as a highly granular reference for the study of immune responses in airways and blood in children.

Factor V is an immune inhibitor that is expressed at increased levels in leukocytes of patients with severe Covid-19

Severe Covid-19 is associated with elevated plasma Factor V (FV) and increased risk of thromboembolism. We report that neutrophils, T regulatory cells (Tregs), and monocytes from patients with severe Covid-19 express FV, and expression correlates with T cell lymphopenia. In vitro full length FV, but not FV activated by thrombin cleavage, suppresses T cell proliferation. Increased and prolonged FV expression by cells of the innate and adaptive immune systems may contribute to lymphopenia in severe Covid-19. Activation by thrombin destroys the immunosuppressive properties of FV. Anticoagulation in Covid-19 patients may have the unintended consequence of suppressing the adaptive immune system.

Early immune pathology and persistent dysregulation characterise severe COVID-19

In a study of 207 SARS-CoV2-infected individuals with a range of severities followed over 12 weeks from symptom onset, we demonstrate that an early robust bystander CD8 T cell immune response, without systemic inflammation, is characteristic of asymptomatic or mild disease. Those presenting to hospital had delayed bystander responses and systemic inflammation already evident at around symptom onset. Such early evidence of inflammation suggests immunopathology may be inevitable in some individuals, or that preventative intervention might be needed before symptom onset. Viral load does not correlate with the development of this pathological response, but does with its subsequent severity. Immune recovery is complex, with profound persistent cellular abnormalities correlating with a change in the nature of the inflammatory response, where signatures characteristic of increased oxidative phosphorylation and reactive-oxygen species-associated inflammation replace those driven by TNF and IL-6. These late immunometabolic inflammatory changes and unresolved immune defects may have clinical implications.

The cellular immune response to COVID-19 deciphered by single cell multi-omics across three UK centres

The COVID-19 pandemic, caused by SARS coronavirus 2 (SARS-CoV-2), has resulted in excess morbidity and mortality as well as economic decline. To characterise the systemic host immune response to SARS-CoV-2, we performed single-cell RNA-sequencing coupled with analysis of cell surface proteins, providing molecular profiling of over 800,000 peripheral blood mononuclear cells from a cohort of 130 patients with COVID-19. Our cohort, from three UK centres, spans the spectrum of clinical presentations and disease severities ranging from asymptomatic to critical. Three control groups were included: healthy volunteers, patients suffering from a non-COVID-19 severe respiratory illness and healthy individuals administered with intravenous lipopolysaccharide to model an acute inflammatory response. Full single cell transcriptomes coupled with quantification of 188 cell surface proteins, and T and B lymphocyte antigen receptor repertoires have provided several insights into COVID-19: 1. a new non-classical monocyte state that sequesters platelets and replenishes the alveolar macrophage pool; 2. platelet activation accompanied by early priming towards megakaryopoiesis in immature haematopoietic stem/progenitor cells and expansion of megakaryocyte-primed progenitors; 3. increased clonally expanded CD8+ effector:effector memory T cells, and proliferating CD4+ and CD8+ T cells in patients with more severe disease; and 4. relative increase of IgA plasmablasts in asymptomatic stages that switches to expansion of IgG plasmablasts and plasma cells, accompanied with higher incidence of BCR sharing, as disease severity increases. All data and analysis results are available for interrogation and data mining through an intuitive web portal. Together, these data detail the cellular processes present in peripheral blood during an acute immune response to COVID-19, and serve as a template for multi-omic single cell data integration across multiple centers to rapidly build powerful resources to help combat diseases such as COVID-19.