The severity of SARS-CoV-2 infection has been related to uncontrolled inflammatory innate responses and impaired adaptive immune responses mostly due to exhausted T lymphocytes and lymphopenia. In this work we have characterized the nature of the lymphopenia and demonstrate a set of factors that hinder the effective control of virus infection and the activation and arming of effector cytotoxic T CD8 cells and showing signatures defining a high-risk population. We performed immune profiling of the T helper (Th) CD4+ and T CD8+ cell compartments in peripheral blood of 144 COVID-19 patients using multiparametric flow cytometry analysis. On the one hand, there was a consistent lymphopenia with an overrepresentation of non-functional T cells, with an increased percentage of naive Th cells (CD45RA+, CXCR3-, CCR4-, CCR6-, CCR10-) and persistently low frequency of markers associated with Th1, Th17, and Th1/Th17 memory-effector T cells compared to healthy donors. On the other hand, the most profound alteration affected the Th1 subset, which may explain the poor T cells responses and the persistent blood virus load. Finally, the decrease in Th1 cells may also explain the low frequency of CD4+ and CD8+ T cells that express the HLA-DR and CD38 activation markers observed in numerous patients who showed minimal or no lymphocyte activation response. We also identified the percentage of HLA-DR+CD4+ T cells, PD-1+CD+4/CD8+ T cells in blood, and the neutrophil/lymphocyte ratio as useful factors for predicting critical illness and fatal outcome in patients with confirmed COVID-19.

Critical illness in COVID-19 is an extreme and clinically homogeneous disease phenotype that we have previously shown1 to be highly efficient for discovery of genetic associations2. Despite the advanced stage of illness at presentation, we have shown that host genetics in patients who are critically ill with COVID-19 can identify immunomodulatory therapies with strong beneficial effects in this group3. Here we analyse 24,202 cases of COVID-19 with critical illness comprising a combination of microarray genotype and whole-genome sequencing data from cases of critical illness in the international GenOMICC (11,440 cases) study, combined with other studies recruiting hospitalized patients with a strong focus on severe and critical disease: ISARIC4C (676 cases) and the SCOURGE consortium (5,934 cases). To put these results in the context of existing work, we conduct a meta-analysis of the new GenOMICC genome-wide association study (GWAS) results with previously published data. We find 49 genome-wide significant associations, of which 16 have not been reported previously. To investigate the therapeutic implications of these findings, we infer the structural consequences of protein-coding variants, and combine our GWAS results with gene expression data using a monocyte transcriptome-wide association study (TWAS) model, as well as gene and protein expression using Mendelian randomization. We identify potentially druggable targets in multiple systems, including inflammatory signalling (JAK1), monocyte–macrophage activation and endothelial permeability (PDE4A), immunometabolism (SLC2A5 and AK5), and host factors required for viral entry and replication (TMPRSS2 and RAB2A).

Here we describe the results of a genome-wide study conducted in 11 939 COVID-19 positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (p < 5x10−8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (p = 1.3x10−22 and p = 8.1x10−12, respectively), and for variants in 9q21.32 near TLE1 only among females (p = 4.4x10−8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (p = 2.7x10−8) and ARHGAP33 (p = 1.3x10−8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, p = 4.1x10−8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥ 60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

Human leukocyte antigen (HLA) plays an important role in immune responses to infections, especially in the development of acquired immunity. Given the high degree of polymorphisms that HLA molecules present, some will be more or less effective in controlling SARS-CoV-2 infection. We wanted to analyze whether certain polymorphisms may be involved in the protection or susceptibility to COVID-19. To do this, we studied the polymorphisms in HLA class I (HLA-A, -B and -C) and II (HLA-DRB1 and HLA-DQB1) molecules in 450 patients who required hospitalization for COVID-19, creating one of the largest HLA-typed patient cohort to date. Our results show that there is no relationship between HLA polymorphisms or haplotypes and susceptibility or protection to COVID-19. Our results may contribute to resolve the contradictory data on the role of HLA polymorphisms in COVID-19 infection.

The vaccines designed against the SARS-CoV-2 coronavirus are based on the spike (S) protein. Processing of the S protein by antigen-presenting cells (APC) and its subsequent presentation to T cells is an essential part of the development of a humoral response. HLA-class II alleles are considered immune response genes because their codified molecules, expressed on the surface of APCs (macrophages, dendritic, and B cells) present antigenic peptides to T cell via their T cell receptor (TCR). The HLA-class II genes are highly polymorphic, regulating what specific peptides induce follicular helper T cells (TFH) and promote B lymphocyte differentiation into plasma or memory B cells. This work hypothesizes that the presence of certain HLA-class II alleles could be associated with the intensity of the humoral response (amount, length) to the SARS-CoV2 mRNA 1273 vaccine. We have studied the relationship between the HLA-class II typing of 87 health workers and the level of antibodies produced 30 days after vaccination. We show a possible association between the HLA-DRB1* 07:01 allele and the HLA-DRB1*07:01~DQA1*02:01~DQB1*02:02 haplotype to a higher production of antibodies 30 days after the administration of the second dose of mRNA-1273.