Background Somatic syndrome is one of the remarkably prevalent issues in primary health care and subspecialty settings. We aimed to elucidate multidimensional associations between somatic symptoms with major mental problems and personality traits in the framework of the quantile regression model with a Bayesian approach. Methods A total of 4763 employees at Isfahan University of Medical Sciences and Health Services in Isfahan province, Iran, filled out four validated questionnaires including Hospital Anxiety and Depression Scale (HADS), NEO Questionnaire, General Health Questionnaire (GHQ) and PHQ-15 for somatic symptom severity. In addition, Functional Gastrointestinal Disorders (FGIDs) were determined using Rome IV criteria. Exploratory Factor Analysis (EFA) and Bayesian regularized quantile regression with adaptive LASSO penalization were applied for reduced dimension of somatic symptoms and variable selection and parameter estimation, respectively. Results The 25 major somatic symptoms were grouped into four factors including general, upper gastrointestinal, lower gastrointestinal and respiratory by EFA. Stress, depression, and anxiety had significant effects on all of the four extracted factors. The effect of anxiety in each four extracted factors was more than stress and depression. Neuroticism and agreeableness had significant effects on all of the four extracted factors, generally ( p < 0.05). Conclusions Given the high prevalence of somatic symptoms and psychosomatic complaints in correlation with the diverse range of mental co-morbidities, developing more detailed diagnostic tools and methods is crucial; nonetheless, it seems that providing better interdisciplinary approaches in general medical practice is groundwork.

Knowledge and awareness of the healthcare professionals regarding the sexual problems should be increased. Management of sexual dysfunction in a holistic approach in the primary care services might improve the wellness and quality of life of the women.