Through feature selection processes, LASSO regression identified 21 candidate biomarkers linked to CKM–depression comorbidity (Fig. 2a), while the Boruta algorithm prioritized 54 variables (Fig. 2b). This evidence concerns the gene CKM and major depressive disorder.