Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge1–5. Here we conducted a genome-wide association study (GWAS) involving 2,393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3,289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target.

Consecutive mRNA vaccinations against SARS-CoV-2 reinforced both innate and adaptive immune responses. However, it remains unclear whether the enhanced innate immune responses are mediated by epigenetic regulation and, if so, whether these effects persist. Using mass cytometry, RNA-seq, and ATAC-seq, we show that BNT162b2 mRNA vaccination upregulated antiviral and IFN-stimulated gene expression in monocytes with greater effects after the second vaccination than those after the first vaccination. Transcription factor-binding motif analysis also revealed enriched IFN regulatory factors and PU.1 motifs in accessible chromatin regions. Importantly, although consecutive BNT162b2 mRNA vaccinations boosted innate immune responses and caused epigenetic changes in isolated monocytes, we showed that these effects occur only transiently and disappear 4 weeks after the second vaccination. Furthermore, single-cell RNA sequencing analysis revealed that a similar gene signature was impaired in the monocytes of unvaccinated COVID-19 patients with acute respiratory distress syndrome. These results reinforce the importance of the innate immune response in the determination of COVID-19 severity but indicate that, unlike adaptive immunity, innate immunity is not unexpectedly sustained even after consecutive vaccination. This study, which focuses on innate immmune memory, may provide novel insights into the vaccine development against infectious diseases.

Sex-biased humoral immune responses to COVID-19 patients have been observed, but the cellular basis for this is not understood. Using single-cell proteomics by mass cytometry, we find disrupted regulation of humoral immunity in COVID-19 patients, with a sex-biased loss of circulating follicular regulatory T cells (cTfr) at a significantly greater rate in male patients. In addition, a male sex-associated cellular network of T-peripheral helper, plasma blasts, proliferating and extrafollicular/atypical CD11c + memory B cells was strongly positively correlated with neutralizing antibody concentrations and negatively correlated with cTfr frequency. These results suggest that sex-specific differences to the balance of cTfr and a network of extrafollicular antibody production-associated cell types may be a key factor in the altered humoral immune responses between male and female COVID-19 patients.

Coronavirus disease 2019 (COVID-19) is a recently-emerged infectious disease that has caused millions of deaths, where comprehensive understanding of disease mechanisms is still unestablished. In particular, studies of gene expression dynamics and regulation landscape in COVID-19 infected individuals are limited. Here, we report on a thorough analysis of whole blood RNA-seq data from 465 genotyped samples from the Japan COVID-19 Task Force, including 359 severe and 106 non-severe COVID-19 cases. We discover 1169 putative causal expression quantitative trait loci (eQTLs) including 34 possible colocalizations with biobank fine-mapping results of hematopoietic traits in a Japanese population, 1549 putative causal splice QTLs (sQTLs; e.g. two independent sQTLs at TOR1AIP1), as well as biologically interpretable trans-eQTL examples (e.g., REST and STING1), all fine-mapped at single variant resolution. We perform differential gene expression analysis to elucidate 198 genes with increased expression in severe COVID-19 cases and enriched for innate immune-related functions. Finally, we evaluate the limited but non-zero effect of COVID-19 phenotype on eQTL discovery, and highlight the presence of COVID-19 severity-interaction eQTLs (ieQTLs; e.g., CLEC4C and MYBL2). Our study provides a comprehensive catalog of whole blood regulatory variants in Japanese, as well as a reference for transcriptional landscapes in response to COVID-19 infection.