The final model of highest predictive quality (Hosmer-Lemeshow test: χ2(7) = 4.016, p = 0.778; Nagelkerke R2 = 75.2%, correctly classified 83.3% of cases) showed that females carrying both major alleles of TMPRSS2 rs2070788 were less likely to develop severe COVID-19. Here, TMPRSS2 is linked to COVID-19.