Background: A personalized approach to occupational medicine allows specialists to 1
prevent professional hazards such as stress-related depression and anxiety in extreme work environments. Objective: we aim to detect genetic markers of low resilience to stress. Methods: The study cohort included 97 elite athletes and 167 special forces personnel. The research team collected buccal mucosa samples and examined psychological status with the Hospital Anxiety and Depression Scale (HADS). We assessed 35 variants within selected genes that are most often associated with low resilience to stress, anxiety, and depression. Fisher’s Exact test was used to determine nonrandom associations between scores in the HADS scale and the genetic variants. We also trained machine learning models to predict score values from genotyping findings and ranked genetic biomarkers according to their predictive power. Results: High-risk depression profiles included C/T genotype of MTHFR C677T and A/C variant of MTHFR A1298C. Susceptibility to anxiety was associated with several polymorphisms regulating neuroactive substances, immune response, and coagulation. The ML models accurately detect depression or anxiety levels with MAE/ROV of 17.69±1.35 and 17.86±2.09% respectively. Conclusions: The study findings justify a polymorphic nature of anxiety and confirm the immune system's involvement in regulating stress response.