Regarding subcortical regions, the GO analysis revealed overrepresentation of 23 genes involved in the biological process of neuronal development (Table 2), including NeuroD2. Apart from a high correlation with the BPND of the 5-HT1AR (ρ = 0.81), NeuroD2 also showed a similar SD (SD = 2.13) as the BPND of the 5-HT1AR. Alterations in NeuroD2 have been linked to schizophrenia, where patients with NeuroD2 polymorphisms showed reduced cognitive functions (Spellmann et al. 2017). The gene discussed is NEUROD2; the disease is schizophrenia.