Understanding inter-individual differences in stress response requires the explanation of genetic influences at multiple phenotypic levels, including complex behaviours and the metabolic responses of brain regions to emotional stimuli. Neuropeptide Y (NPY) is anxiolytic 1,2 and its release is induced by stress 3 . NPY is abundantly expressed in regions of the limbic system that are implicated in arousal and in the assignment of emotional valences to stimuli and memories [4][5][6] . Here we show that haplotype-driven NPY expression predicts brain responses to emotional and stress challenges and also inversely correlates with trait anxiety. NPY haplotypes predicted levels of NPY messenger RNA in postmortem brain and lymphoblasts, and levels of plasma NPY. Lower haplotype-driven NPY expression predicted higher emotion-induced activation of the amygdala, as well as diminished resiliency as assessed by pain/stress-induced activations of endogenous opioid neurotransmission in various brain regions. A single nucleotide polymorphism (SNP rs16147) located in the promoter region alters NPY expression in vitro and seems to account for more than half of the variation in expression in vivo. These convergent findings are consistent with the function of NPY as an anxiolytic peptide and help to explain inter-individual variation in resiliency to stress, a risk factor for many diseases.Reprints and permissions information is available at www.nature.com/reprints. Correspondence and requests for materials should be addressed to D.G. (E-mail: davidgoldman@mail.nih.gov). * These authors contributed equally to this work. † Present address: Innovation Centre China, AstraZeneca Global R&D, Shanghai 201203, China. Supplementary Fig. 1b). Five haplotypes (H1-H5) account for 93.8% of chromosomes in this block (Fig. 1a).We observed haplotype-driven NPY mRNA expression in postmortem brain (US Caucasians, Miami sample) by detecting the differential expression of alleles at single nucleotide polymorphism (SNP) rs5574 C/T, selected because of its high frequency and location in the transcript. Of these 28 samples, chosen because all were heterozygous for rs5574, 16 (57%) showed differential allele expression at an allele ratio of more than 1.2, in either direction. H1 and H4 were low-expression haplotypes, H2 was high, H3 was intermediate and H5 was unclassified because only two H1/H5 heterozygous brains were available (Fig. 1b). This expression-based functional classification is consistent with a cladistically based clustering of haplotypes, indicating that expression variation is linked to gene ancestry (Fig. 1a). The effects on expression of the more common H1, H2 and H3 haplotypes were verified in 47 lymphoblastoid cell lines derived from healthy Finnish men ( Fig. 1c) representing the six common diplotypes (72% of all diplotypes). On the basis of lymphoblast NPY mRNA levels, the expression value for each haplotype was calculated by regression analysis. Expression values for the six common diplotypes were well predicted under a co...