To identify genes involved in anxiety/fear traits, we analyzed the gene expression profile in the amygdala of genetically heterogeneous NIH-HS rats. The NIH-HS rat stock has revealed to be a unique genetic resource for the fine mapping of Quantitative Trait Loci (QTLs) to very small genomic regions, due to the high amount of genetic recombinants accumulated along more than 50 breeding generations, and for the same reason it can be expected that those genetically heterogeneous rats should be especially useful for studying differential gene expression as a function of anxiety-(or other)-related traits. We selected high- and low-anxious NIH-HS rats differing in their number of avoidances in a single 50-trial session of the two-way active avoidance task. Rats were also tested in unconditioned anxiety tests (e.g. elevated zero-maze). Three weeks after behavioural testing, the amygdala was dissected and prepared for the microarray study. There appeared 6 significantly down-regulated and 28 up-regulated genes (fold-change > |2|, FDR < 0.05) between the low- and high-anxious groups, with central nervous system-related functions. Regression analyses (stepwise) revealed that differential expression of some genes could be predictive of anxiety/fear responses. Among those genes for which the present results suggest a link with individual differences in trait anxiety, six relevant genes were examined with qRT-PCR, four of which (Ucn3, Tacr3, H2-M9 and Arr3) were validated. Remarkably, some of them are characterized by sharing known functions related with hormonal HPA-axis responses to (and/or modulation of) stress, anxiety or fear, and putative involvement in related neurobehavioural functions. The results confirm the usefulness of NIH-HS rats as a good animal model for research on the neurogenetic basis of anxiety and fear, while suggesting the involvement of some neuropeptide/neuroendocrine pathways on the development of differential anxiety profiles.