the GWAS summary statistics. PGS were calculated from the imputed genotype dosages to the 1000 Genomes (Phase 3 Release 5) reference panel. We fitted linear mixed models, which controlled for relatedness using a Genetic Relatedness Matrix (GRM) and covariates sex, age, two dummy variables for the GWAS array used, and the first five genetic principal components. The parameters of the model were estimated using GCTA 1.9 [36]. The linear @@model was as follows:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$CD{\rm{symptom}}\,{\rm{score}} = {\rm{intercept}} + {\rm{Covariates}} \ast b + c \ast PGS + G$$\end{document}CDsymptomscore=intercept+Covariates*b+c*PGS+G
