We conducted whole-exome sequencing in 233 hospitalized COVID-19 patients and identified four PID gene (UNC13D, AP3B1, RNF168, DHX58) variants were significantly enriched in COVID-19 patients experiencing severe cytokine storms. Here, UNC13D is linked to COVID-19.