Methods for modeling gene-environment interplay using polygenic risk scores

Polygenic risk scores (PRS) are increasingly recognized as pivotal tools for quantifying disease risk through the aggregation of multiple genetic variants. As sample sizes in genome-wide association studies (GWAS) continue to expand and PRS become more powerful, they are set to play a key role in translational research and personalized medicine. Understanding the interplay of PRS with environmental factors is critical for interpreting and applying PRS in diverse contexts. This interplay manifests in two forms: PRS-by-environment interaction (PRS × E) and gene-environment correlation (rGE). However, despite the growing application and importance of PRS, there are limited guidelines for performing PRS × E interaction analyses while controlling for rGE, which can lead to inconsistencies across studies and misinterpretation of results. Here we provide a review of different methods for performing PRSxE interaction in various epidemiological study designs, propose recommendations for best-practice, and discuss future challenges.