Then, to investigate if plasma levels of VIP were associated with functional connectivity in the selected regions, we performed a multiple regression analysis between VIP and functional connectivity values between bilateral amygdala seed regions and all target regions including the hippocampus, parahippocampal gyri, and lateral and medial OFC regions while controlling for both HADS anxiety and depression. This evidence concerns the gene VIP and depressive symptom measurement.