When osteocalcin was introduced as the dependent variable, both COVID-19 diagnosis (β = −0.375, p = 0.005) and the presence of stress hyperglycemia (β = −0.344, p = 0.009) were retained by the model (R2 = 0.347, F = 12.224, p < 0.001). This evidence concerns the gene BGLAP and COVID-19.