Early life stress may precipitate psychopathology, at least in part, by influencing amygdala function. Converging evidence across species suggests that links between childhood stress and amygdala function may be dependent upon hypothalamic-pituitary-adrenal (HPA) axis function. Using data from college-attending non-Hispanic European-Americans (n=308) who completed the Duke Neurogenetics Study, we examined whether early life stress (ELS) and HPA axis genetic variation interact to predict threat-related amygdala function as well as psychopathology symptoms. A biologically-informed multilocus profile score (BIMPS) captured HPA axis genetic variation (FKBP5 rs1360780, CRHR1 rs110402; NR3C2 rs5522/rs4635799) previously associated with its function (higher BIMPS are reflective of higher HPA axis activity). BOLD fMRI data were acquired while participants completed an emotional face matching task. ELS and depression and anxiety symptoms were measured using the childhood trauma questionnaire and the mood and anxiety symptom questionnaire, respectively. The interaction between HPA axis BIMPS and ELS was associated with right amygdala reactivity to threat-related stimuli, after accounting for multiple testing (empirical-p=0.016). Among individuals with higher BIMPS (i.e., the upper 21.4%), ELS was positively coupled with threat-related amygdala reactivity, which was absent among those with average or low BIMPS. Further, higher BIMPS were associated with greater self-reported anxious arousal, though there was no evidence that amygdala function mediated this relationship. Polygenic variation linked to HPA axis function may moderate the effects of early life stress on threat-related amygdala function and confer risk for anxiety symptomatology. However, what, if any, neural mechanisms may mediate the relationship between HPA axis BIMPS and anxiety symptomatology remains unclear.