Calculator for the potential benefit of embryo selection experiments using Polygenic Risk Scores
This repository contains R code used to calculate the possible gain in phenotype value achievable in embryo selection experiments under Gaussian models for the genetic architecture. The repository accompanies the paper [1]. Please cite this paper if you're using the package.
clone the repository into a directory of your choice, and start an R session within this directory.
After running: G <- gain.moments(10, 0.3, 'exact', 1)
G$E.G will contain the expected gain (in units of standard deviations) in selecting the top quantitative trait out of n=10 embryos, for selection based on a polygenic score which explaines r2ps=0.3 of the trait's variance. G$Var.G will contain the variance of this gain.
For any questions, please contact Shai Carmi (shai.carmi@mail.huji.ac.il) or Or Zuk (or.zuk@mail.huji.ac.il)
[1] Screening Human Embryos for Polygenic Traits has Limited Utility
E. Karavani, O. Zuk, D. Zeevi, G. Atzmon, N. Barzilai, N.C. Stefanis, A. Hatzimanolis, N. Smyrnis, D. Avramopoulos, L. Kruglyak, M. Lam, T. Lencz and S. Carmi , Cell 179(6), p. 1424-1435 (2019)
https://www.sciencedirect.com/science/article/pii/S0092867419312103
[2] Utility of polygenic embryo screening for disease depends on the selection strategy
T. Lencz, D. Backenroth, E. Granot-Hershkovitz, A. Green, K. Gettler, J.H. Cho, O. Weissbrod, O. Zuk and S. Carmi
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8510582/