Integrated analysis of COVID-19 multi-omics data for eQTLs reveals genetic mechanisms underlying disease severity

The global pandemic caused by the SARS-CoV-2 virus provided an unprecedented opportunity to investigate genetic factors influencing the disease severity of the viral infection. Despite a plethora of recent research on both SARS-CoV-2 and COVID-19, few have taken a systems biology approach to address individual-level variation, especially based on non-European populations. Accordingly, we analyzed multi-omics data generated at three timepoints from 193 Korean COVID-19 patients with mild or severe symptoms, composed of whole genome sequencing, blood-based single-cell RNA-sequencing (2.15M cells), 195 cytokine profiles, and human leukocyte antigen (HLA) allele data. We identified expression quantitative trait loci (eQTLs), disease severity interacting eQTLs (n = 388), and disease progression interacting eQTLs (n = 945) for various cell types. We elucidated a complex regulatory mechanism involving HLA genes and their targets, and identified genetic determinants of cytokine levels. Finally, we show how regulation of ieQTLs is established by upstream transcription factors (TFs), illustrating complex regulation of the IGFBP7 ieQTL by a combined action of two TFs, which is potentially important in conferring differential severity. This study illuminates an efficient molecular interrogation framework that can be applied toward understanding infectious disease progression in individuals of different genotypes.