Abstract The present study examines whether neuroticism is predicted by genetic vulnerability, summarized as polygenic risk score for neuroticism (PRS _N ), in interaction with bullying, parental bonding, and childhood adversity. Data were derived from a general population adolescent and young adult twin cohort. The final sample consisted of 202 monozygotic and 436 dizygotic twins and 319 twin pairs. The Short Eysenck Personality questionnaire was used to measure neuroticism. PRS _N was trained on the results from the Genetics of Personality Consortium (GPC) and United Kingdom Biobank (UKB) cohorts, yielding two different PRS _N . Multilevel mixed-effects models were used to analyze the main and interacting associations of PRS _N , childhood adversity, bullying, and parental bonding style with neuroticism. We found no evidence of gene–environment correlation. PRS _N thresholds of .005 and .2 were chosen, based on GPC and UKB datasets, respectively. After correction for confounders, all the individual variables were associated with the expression of neuroticism: both PRS _N from GPC and UKB, childhood adversity, maternal bonding, paternal bonding, and bullying in primary school and secondary school. However, the results indicated no evidence for gene–environment interaction in this cohort. These results suggest that genetic vulnerability on the one hand and negative life events (childhood adversity and bullying) and positive life events (optimal parental bonding) on the other represent noninteracting pathways to neuroticism.