AbstractEmotional stress is a leading risk factor in the development of neuropsychiatric disorders possibly via immune activation. P2X7 receptors promote neuroinflammation, and research suggests a relationship between chromosome region 12q2431, in which the P2X7R gene is located, and development of mood disorders, however, few studies concentrate on its association with anxiety. Our aim was to investigate the effects of P2RX7 variation in interaction with early childhood traumas and recent stressors on anxiety. 1752 participants completed questionnaires assessing childhood adversities and recent negative life events, provided data on anxiety using the Brief Symptom Inventory, and were genotyped for 681 SNPs in the P2RX7 gene, 335 of which passed quality control and were entered into linear regression models followed by a linkage disequilibrium-based clumping procedure to identify clumps of SNPs with a significant main or interaction effect. We identified a significant clump with top SNP rs67881993 and containing a set of 29SNPs that are in high LD, which significantly interacted with early childhood traumas but not with recent stress conveying a protective effect against increased anxiety in those exposed to early adversities. Our study demonstrated that P2RX7 variants interact with distal and more etiological stressors in influencing the severity of anxiety symptoms, supporting previous scarce results and demonstrating its role in moderating the effects of stress.