In model 1 of the multinomial regression using the PHQ-9 as a continuous variable, there was a significant effect of depression symptoms on glycaemic control (B = 0.03, p = 0.039), after adjusting for variables which had a significant association with PHQ-9 as potential confounders (namely, gender, BMI, prescribed insulin, having past history of depression and hypertension, and not following any physical activity programme (see Table 2). This evidence concerns the gene INS and depressive disorder.