Polygenic scoring accuracy varies across the genetic ancestry continuum.

Polygenic scores (PGSs) have limited portability across different groupings of individuals (for example, by genetic ancestries and/or social determinants of health), preventing their equitable use1-3. PGS portability has typically been assessed using a single aggregate population-level statistic (for example, R2)4, ignoring inter-individual variation within the population. Here, using a large and diverse Los Angeles biobank5 (ATLAS, n = 36,778) along with the UK Biobank6 (UKBB, n = 487,409), we show that PGS accuracy decreases individual-to-individual along the continuum of genetic ancestries7 in all considered populations, even within traditionally labelled 'homogeneous' genetic ancestries. The decreasing trend is well captured by a continuous measure of genetic distance (GD) from the PGS training data: Pearson correlation of -0.95 between GD and PGS accuracy averaged across 84 traits. When applying PGS models trained on individuals labelled as white British in the UKBB to individuals with European ancestries in ATLAS, individuals in the furthest GD decile have 14% lower accuracy relative to the closest decile; notably, the closest GD decile of individuals with Hispanic Latino American ancestries show similar PGS performance to the furthest GD decile of individuals with European ancestries. GD is significantly correlated with PGS estimates themselves for 82 of 84 traits, further emphasizing the importance of incorporating the continuum of genetic ancestries in PGS interpretation. Our results highlight the need to move away from discrete genetic ancestry clusters towards the continuum of genetic ancestries when considering PGSs.

Towards fair and clinically relevant polygenic predictions.

Lennon et al. recently proposed a clinical polygenic score (PGS) pipeline as part of the Electronic Medical Records and Genomics (eMERGE) network initiative. In this spotlight article we discuss the broader context for the use of PGS in preventive medicine and highlight key limitations and challenges facing their inclusion in prediction models.

Genetic liability estimated from large-scale family data improves genetic prediction, risk score profiling, and gene mapping for major depression.

Large biobank samples provide an opportunity to integrate broad phenotyping, familial records, and molecular genetics data to study complex traits and diseases. We introduce Pearson-Aitken Family Genetic Risk Scores (PA-FGRS), a method for estimating disease liability from patterns of diagnoses in extended, age-censored genealogical records. We then apply the method to study a paradigmatic complex disorder, major depressive disorder (MDD), using the iPSYCH2015 case-cohort study of 30,949 MDD cases, 39,655 random population controls, and more than 2 million relatives. We show that combining PA-FGRS liabilities estimated from family records with molecular genotypes of probands improves three lines of inquiry. Incorporating PA-FGRS liabilities improves classification of MDD over and above polygenic scores, identifies robust genetic contributions to clinical heterogeneity in MDD associated with comorbidity, recurrence, and severity and can improve the power of genome-wide association studies. Our method is flexible and easy to use, and our study approaches are generalizable to other datasets and other complex traits and diseases.

Inferring disease architecture and predictive ability with LDpred2-auto.

LDpred2 is a widely used Bayesian method for building polygenic scores (PGSs). LDpred2-auto can infer the two parameters from the LDpred model, the SNP heritability h2 and polygenicity p, so that it does not require an additional validation dataset to choose best-performing parameters. The main aim of this paper is to properly validate the use of LDpred2-auto for inferring multiple genetic parameters. Here, we present a new version of LDpred2-auto that adds an optional third parameter α to its model, for modeling negative selection. We then validate the inference of these three parameters (or two, when using the previous model). We also show that LDpred2-auto provides per-variant probabilities of being causal that are well calibrated and can therefore be used for fine-mapping purposes. We also introduce a formula to infer the out-of-sample predictive performance r2 of the resulting PGS directly from the Gibbs sampler of LDpred2-auto. Finally, we extend the set of HapMap3 variants recommended to use with LDpred2 with 37% more variants to improve the coverage of this set, and we show that this new set of variants captures 12% more heritability and provides 6% more predictive performance, on average, in UK Biobank analyses.

Portability of 245 polygenic scores when derived from the UK Biobank and applied to 9 ancestry groups from the same cohort.

The low portability of polygenic scores (PGSs) across global populations is a major concern that must be addressed before PGSs can be used for everyone in the clinic. Indeed, prediction accuracy has been shown to decay as a function of the genetic distance between the training and test cohorts. However, such cohorts differ not only in their genetic distance but also in their geographical distance and their data collection and assaying, conflating multiple factors. In this study, we examine the extent to which PGSs are transferable between ancestries by deriving polygenic scores for 245 curated traits from the UK Biobank data and applying them in nine ancestry groups from the same cohort. By restricting both training and testing to the UK Biobank data, we reduce the risk of environmental and genotyping confounding from using different cohorts. We define the nine ancestry groups at a sub-continental level, based on a simple, robust, and effective method that we introduce here. We then apply two different predictive methods to derive polygenic scores for all 245 phenotypes and show a systematic and dramatic reduction in portability of PGSs trained using Northwestern European individuals and applied to nine ancestry groups. These analyses demonstrate that prediction already drops off within European ancestries and reduces globally in proportion to genetic distance. Altogether, our study provides unique and robust insights into the PGS portability problem.

Accounting for age of onset and family history improves power in genome-wide association studies.

Genome-wide association studies (GWASs) have revolutionized human genetics, allowing researchers to identify thousands of disease-related genes and possible drug targets. However, case-control status does not account for the fact that not all controls may have lived through their period of risk for the disorder of interest. This can be quantified by examining the age-of-onset distribution and the age of the controls or the age of onset for cases. The age-of-onset distribution may also depend on information such as sex and birth year. In addition, family history is not routinely included in the assessment of control status. Here, we present LT-FH++, an extension of the liability threshold model conditioned on family history (LT-FH), which jointly accounts for age of onset and sex as well as family history. Using simulations, we show that, when family history and the age-of-onset distribution are available, the proposed approach yields statistically significant power gains over LT-FH and large power gains over genome-wide association study by proxy (GWAX). We applied our method to four psychiatric disorders available in the iPSYCH data and to mortality in the UK Biobank and found 20 genome-wide significant associations with LT-FH++, compared to ten for LT-FH and eight for a standard case-control GWAS. As more genetic data with linked electronic health records become available to researchers, we expect methods that account for additional health information, such as LT-FH++, to become even more beneficial.

Deep integrative models for large-scale human genomics.

Polygenic risk scores (PRSs) are expected to play a critical role in precision medicine. Currently, PRS predictors are generally based on linear models using summary statistics, and more recently individual-level data. However, these predictors mainly capture additive relationships and are limited in data modalities they can use. We developed a deep learning framework (EIR) for PRS prediction which includes a model, genome-local-net (GLN), specifically designed for large-scale genomics data. The framework supports multi-task learning, automatic integration of other clinical and biochemical data, and model explainability. When applied to individual-level data from the UK Biobank, the GLN model demonstrated a competitive performance compared to established neural network architectures, particularly for certain traits, showcasing its potential in modeling complex genetic relationships. Furthermore, the GLN model outperformed linear PRS methods for Type 1 Diabetes, likely due to modeling non-additive genetic effects and epistasis. This was supported by our identification of widespread non-additive genetic effects and epistasis in the context of T1D. Finally, we constructed PRS models that integrated genotype, blood, urine, and anthropometric data and found that this improved performance for 93% of the 290 diseases and disorders considered. EIR is available at https://github.com/arnor-sigurdsson/EIR.

Leveraging both individual-level genetic data and GWAS summary statistics increases polygenic prediction.

The accuracy of polygenic risk scores (PRSs) to predict complex diseases increases with the training sample size. PRSs are generally derived based on summary statistics from large meta-analyses of multiple genome-wide association studies (GWASs). However, it is now common for researchers to have access to large individual-level data as well, such as the UK Biobank data. To the best of our knowledge, it has not yet been explored how best to combine both types of data (summary statistics and individual-level data) to optimize polygenic prediction. The most widely used approach to combine data is the meta-analysis of GWAS summary statistics (meta-GWAS), but we show that it does not always provide the most accurate PRS. Through simulations and using 12 real case-control and quantitative traits from both iPSYCH and UK Biobank along with external GWAS summary statistics, we compare meta-GWAS with two alternative data-combining approaches, stacked clumping and thresholding (SCT) and meta-PRS. We find that, when large individual-level data are available, the linear combination of PRSs (meta-PRS) is both a simple alternative to meta-GWAS and often more accurate.

Exploring the Genetic Risk of Childhood Daytime Urinary Incontinence: A Genome-Wide Association Study.

PURPOSE: Childhood incontinence is stigmatized and underprioritized, and a basic understanding of its pathogenesis is missing. Our goal was to identify risk-conferring genetic variants in daytime urinary incontinence (DUI). MATERIALS AND METHODS: We conducted a genome-wide association study in the Danish iPSYCH2015 cohort. Cases (3024) were identified through DUI diagnosis codes and redeemed prescriptions for DUI medication in patients aged 5 to 20 years. Controls (30,240), selected from the same sample, were matched to cases on sex and psychiatric diagnoses, if any, and down-sampled to a 1:10 case:control ratio. Replication was performed in the Icelandic deCODE cohort (5475 cases/287,773 controls). Single-nucleotide polymorphism heritability was calculated using the genome-based restricted maximum likelihood method. Cross-trait genetic correlation was estimated using linkage disequilibrium score regression. Polygenic risk scores generated with LDpred2-auto and BOLT-LMM were assessed for association. RESULTS: Variants on chromosome 6 (rs12210989, odds ratio [OR] 1.24, 95% CI 1.17-1.32, P = 3.21 × 10-12) and 20 (rs4809801, OR 1.18, 95% CI 1.11-1.25, P = 3.66 × 10-8) reached genome-wide significance and implicated the PRDM13 and RIPOR3 genes. Chromosome 6 findings were replicated (P = .024, OR 1.09, 95% CI 1.01-1.16). Liability scale heritability ranged from 10.20% (95% CI 6.40%-14.00%) to 15.30% (95% CI 9.66%-20.94%). DUI and nocturnal enuresis showed positive genetic correlation (rg = 1.28 ± 0.38, P = .0007). DUI was associated with attention-deficit/hyperactivity disorder (OR 1.098, 95% CI 1.046-1.152, P < .0001) and BMI (OR 1.129, 95% CI 1.081-1.178, P < .0001) polygenic risk. CONCLUSIONS: Common genetic variants contribute to the risk of childhood DUI, and genes important in neuronal development and detrusor smooth muscle activity were implicated. These findings may help guide identification of new treatment targets.

Polygenic liabilities and treatment trajectories in early-onset depression: a Danish register-based study.

BACKGROUND: The clinical course of major depressive disorder (MDD) is heterogeneous, and early-onset MDD often has a more severe and complex clinical course. Our goal was to determine whether polygenic scores (PGSs) for psychiatric disorders are associated with treatment trajectories in early-onset MDD treated in secondary care. METHODS: Data were drawn from the iPSYCH2015 sample, which includes all individuals born in Denmark between 1981 and 2008 who were treated in secondary care for depression between 1995 and 2015. We selected unrelated individuals of European ancestry with an MDD diagnosis between ages 10-25 (N = 10577). Seven-year trajectories of hospital contacts for depression were modeled using Latent Class Growth Analysis. Associations between PGS for MDD, bipolar disorder, schizophrenia, ADHD, and anorexia and trajectories of MDD contacts were modeled using multinomial logistic regressions. RESULTS: We identified four trajectory patterns: brief contact (65%), prolonged initial contact (20%), later re-entry (8%), and persistent contact (7%). Relative to the brief contact trajectory, higher PGS for ADHD was associated with a decreased odds of membership in the prolonged initial contact (odds ratio = 1.06, 95% confidence interval = 1.01-1.11) and persistent contact (1.12, 1.03-1.21) trajectories, while PGS-AN was associated with increased odds of membership in the persistent contact trajectory (1.12, 1.03-1.21). CONCLUSIONS: We found significant associations between polygenic liabilities for psychiatric disorders and treatment trajectories in patients with secondary-treated early-onset MDD. These findings help elucidate the relationship between a patient's genetics and their clinical course; however, the effect sizes are small and therefore unlikely to have predictive value in clinical settings.

Interplay of polygenic liability with birth-related, somatic, and psychosocial factors in anorexia nervosa risk: a nationwide study.

BACKGROUND: Although several types of risk factors for anorexia nervosa (AN) have been identified, including birth-related factors, somatic, and psychosocial risk factors, their interplay with genetic susceptibility remains unclear. Genetic and epidemiological interplay in AN risk were examined using data from Danish nationwide registers. AN polygenic risk score (PRS) and risk factor associations, confounding from AN PRS and/or parental psychiatric history on the association between the risk factors and AN risk, and interactions between AN PRS and each level of target risk factor on AN risk were estimated. METHODS: Participants were individuals born in Denmark between 1981 and 2008 including nationwide-representative data from the iPSYCH2015, and Danish AN cases from the Anorexia Nervosa Genetics Initiative and Eating Disorder Genetics Initiative cohorts. A total of 7003 individuals with AN and 45 229 individuals without a registered AN diagnosis were included. We included 22 AN risk factors from Danish registers. RESULTS: Risk factors showing association with PRS for AN included urbanicity, parental ages, genitourinary tract infection, and parental socioeconomic factors. Risk factors showed the expected association to AN risk, and this association was only slightly attenuated when adjusted for parental history of psychiatric disorders or/and for the AN PRS. The interaction analyses revealed a differential effect of AN PRS according to the level of the following risk factors: sex, maternal age, genitourinary tract infection, C-section, parental socioeconomic factors and psychiatric history. CONCLUSIONS: Our findings provide evidence for interactions between AN PRS and certain risk-factors, illustrating potential diverse risk pathways to AN diagnosis.

Shared familial risk for type 2 diabetes mellitus and psychiatric disorders: a nationwide multigenerational genetics study.

BACKGROUND: Psychiatric disorders and type 2 diabetes mellitus (T2DM) are heritable, polygenic, and often comorbid conditions, yet knowledge about their potential shared familial risk is lacking. We used family designs and T2DM polygenic risk score (T2DM-PRS) to investigate the genetic associations between psychiatric disorders and T2DM. METHODS: We linked 659 906 individuals born in Denmark 1990-2000 to their parents, grandparents, and aunts/uncles using population-based registers. We compared rates of T2DM in relatives of children with and without a diagnosis of any or one of 11 specific psychiatric disorders, including neuropsychiatric and neurodevelopmental disorders, using Cox regression. In a genotyped sample (iPSYCH2015) of individuals born 1981-2008 (n = 134 403), we used logistic regression to estimate associations between a T2DM-PRS and these psychiatric disorders. RESULTS: Among 5 235 300 relative pairs, relatives of individuals with a psychiatric disorder had an increased risk for T2DM with stronger associations for closer relatives (parents:hazard ratio = 1.38, 95% confidence interval 1.35-1.42; grandparents: 1.14, 1.13-1.15; and aunts/uncles: 1.19, 1.16-1.22). In the genetic sample, one standard deviation increase in T2DM-PRS was associated with an increased risk for any psychiatric disorder (odds ratio = 1.11, 1.08-1.14). Both familial T2DM and T2DM-PRS were significantly associated with seven of 11 psychiatric disorders, most strongly with attention-deficit/hyperactivity disorder and conduct disorder, and inversely with anorexia nervosa. CONCLUSIONS: Our findings of familial co-aggregation and higher T2DM polygenic liability associated with psychiatric disorders point toward shared familial risk. This suggests that part of the comorbidity is explained by shared familial risks. The underlying mechanisms still remain largely unknown and the contributions of genetics and environment need further investigation.

Phenomewide Association Study of Health Outcomes Associated With the Genetic Correlates of 25 Hydroxyvitamin D Concentration and Vitamin D Binding Protein Concentration.

While it is known that vitamin D deficiency is associated with adverse bone outcomes, it remains unclear whether low vitamin D status may increase the risk of a wider range of health outcomes. We had the opportunity to explore the association between common genetic variants associated with both 25 hydroxyvitamin D (25OHD) and the vitamin D binding protein (DBP, encoded by the GC gene) with a comprehensive range of health disorders and laboratory tests in a large academic medical center. We used summary statistics for 25OHD and DBP to generate polygenic scores (PGS) for 66,482 participants with primarily European ancestry and 13,285 participants with primarily African ancestry from the Vanderbilt University Medical Center Biobank (BioVU). We examined the predictive properties of PGS25OHD, and two scores related to DBP concentration with respect to 1322 health-related phenotypes and 315 laboratory-measured phenotypes from electronic health records. In those with European ancestry: (a) the PGS25OHD and PGSDBP scores, and individual SNPs rs4588 and rs7041 were associated with both 25OHD concentration and 1,25 dihydroxyvitamin D concentrations; (b) higher PGS25OHD was associated with decreased concentrations of triglycerides and cholesterol, and reduced risks of vitamin D deficiency, disorders of lipid metabolism, and diabetes. In general, the findings for the African ancestry group were consistent with findings from the European ancestry analyses. Our study confirms the utility of PGS and two key variants within the GC gene (rs4588 and rs7041) to predict the risk of vitamin D deficiency in clinical settings and highlights the shared biology between vitamin D-related genetic pathways a range of health outcomes.

Multitrait GWAS to connect disease variants and biological mechanisms.

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Association between Mental Disorders and Subsequent Medical Conditions.

BACKGROUND: Persons with mental disorders are at a higher risk than the general population for the subsequent development of certain medical conditions. METHODS: We used a population-based cohort from Danish national registries that included data on more than 5.9 million persons born in Denmark from 1900 through 2015 and followed them from 2000 through 2016, for a total of 83.9 million person-years. We assessed 10 broad types of mental disorders and 9 broad categories of medical conditions (which encompassed 31 specific conditions). We used Cox regression models to calculate overall hazard ratios and time-dependent hazard ratios for pairs of mental disorders and medical conditions, after adjustment for age, sex, calendar time, and previous mental disorders. Absolute risks were estimated with the use of competing-risks survival analyses. RESULTS: A total of 698,874 of 5,940,299 persons (11.8%) were identified as having a mental disorder. The median age of the total population was 32.1 years at entry into the cohort and 48.7 years at the time of the last follow-up. Persons with a mental disorder had a higher risk than those without such disorders with respect to 76 of 90 pairs of mental disorders and medical conditions. The median hazard ratio for an association between a mental disorder and a medical condition was 1.37. The lowest hazard ratio was 0.82 for organic mental disorders and the broad category of cancer (95% confidence interval [CI], 0.80 to 0.84), and the highest was 3.62 for eating disorders and urogenital conditions (95% CI, 3.11 to 4.22). Several specific pairs showed a reduced risk (e.g., schizophrenia and musculoskeletal conditions). Risks varied according to the time since the diagnosis of a mental disorder. The absolute risk of a medical condition within 15 years after a mental disorder was diagnosed varied from 0.6% for a urogenital condition among persons with a developmental disorder to 54.1% for a circulatory disorder among those with an organic mental disorder. CONCLUSIONS: Most mental disorders were associated with an increased risk of a subsequent medical condition; hazard ratios ranged from 0.82 to 3.62 and varied according to the time since the diagnosis of the mental disorder. (Funded by the Danish National Research Foundation and others; COMO-GMC ClinicalTrials.gov number, NCT03847753.).

Genetic liability to posttraumatic stress disorder and its association with postpartum depression.

BACKGROUND: Childbirth may be a traumatic experience and vulnerability to posttraumatic stress disorder (PTSD) may increase the risk of postpartum depression (PPD). We investigated whether genetic vulnerability to PTSD as measured by polygenic score (PGS) increases the risk of PPD and whether a predisposition to PTSD in PPD cases exceeds that of major depressive disorder (MDD) outside the postpartum period. METHODS: This case-control study included participants from the iPSYCH2015, a case-cohort of all singletons born in Denmark between 1981 and 2008. Restricting to women born between 1981 and 1997 and excluding women with a first diagnosis other than depression (N = 22 613), 333 were identified with PPD. For each PPD case, 999 representing the background population and 993 with MDD outside the postpartum were matched by calendar year at birth, cohort selection, and age. PTSD PGS was calculated from summary statistics from the Psychiatric Genomics Consortium with LDpred2-auto. Odds ratios (ORs) were estimated using conditional logistic regression adjusted for parental psychiatric history and country of origin, PGS for MDD and age at first birth, and the first 10 principal components. RESULTS: The PTSD PGS was significantly associated with PPD (OR 1.42, 95% CI 1.20-1.68 per standard deviation increase in PTSD PGS) compared to healthy female controls. Genetic PTSD vulnerability in PPD cases did not exceed that of matched female depression cases outside the postpartum period (OR 1.10, 95% CI 0.94-1.30 per standard deviation increase). CONCLUSIONS: Genetic vulnerability to PTSD increased the risk of PPD but did not differ between PPD cases and women with depression at other times.

Polygenic liability, stressful life events and risk for secondary-treated depression in early life: a nationwide register-based case-cohort study.

BACKGROUND: In this study, we examined the relationship between polygenic liability for depression and number of stressful life events (SLEs) as risk factors for early-onset depression treated in inpatient, outpatient or emergency room settings at psychiatric hospitals in Denmark. METHODS: Data were drawn from the iPSYCH2012 case-cohort sample, a population-based sample of individuals born in Denmark between 1981 and 2005. The sample included 18 532 individuals who were diagnosed with depression by a psychiatrist by age 31 years, and a comparison group of 20 184 individuals. Information on SLEs was obtained from nationwide registers and operationalized as a time-varying count variable. Hazard ratios and cumulative incidence rates were estimated using Cox regressions. RESULTS: Risk for depression increased by 35% with each standard deviation increase in polygenic liability (p < 0.0001), and 36% (p < 0.0001) with each additional SLE. There was a small interaction between polygenic liability and SLEs (ß = -0.04, p = 0.0009). The probability of being diagnosed with depression in a hospital-based setting between ages 15 and 31 years ranged from 1.5% among males in the lowest quartile of polygenic liability with 0 events by age 15, to 18.8% among females in the highest quartile of polygenic liability with 4+ events by age 15. CONCLUSIONS: These findings suggest that although there is minimal interaction between polygenic liability and SLEs as risk factors for hospital-treated depression, combining information on these two important risk factors could potentially be useful for identifying high-risk individuals.

Genome-wide association study of febrile seizures implicates fever response and neuronal excitability genes.

Febrile seizures represent the most common type of pathological brain activity in young children and are influenced by genetic, environmental and developmental factors. In a minority of cases, febrile seizures precede later development of epilepsy. We conducted a genome-wide association study of febrile seizures in 7635 cases and 83 966 controls identifying and replicating seven new loci, all with P < 5 × 10-10. Variants at two loci were functionally related to altered expression of the fever response genes PTGER3 and IL10, and four other loci harboured genes (BSN, ERC2, GABRG2, HERC1) influencing neuronal excitability by regulating neurotransmitter release and binding, vesicular transport or membrane trafficking at the synapse. Four previously reported loci (SCN1A, SCN2A, ANO3 and 12q21.33) were all confirmed. Collectively, the seven novel and four previously reported loci explained 2.8% of the variance in liability to febrile seizures, and the single nucleotide polymorphism heritability based on all common autosomal single nucleotide polymorphisms was 10.8%. GABRG2, SCN1A and SCN2A are well-established epilepsy genes and, overall, we found positive genetic correlations with epilepsies (rg = 0.39, P = 1.68 × 10-4). Further, we found that higher polygenic risk scores for febrile seizures were associated with epilepsy and with history of hospital admission for febrile seizures. Finally, we found that polygenic risk of febrile seizures was lower in febrile seizure patients with neuropsychiatric disease compared to febrile seizure patients in a general population sample. In conclusion, this largest genetic investigation of febrile seizures to date implicates central fever response genes as well as genes affecting neuronal excitability, including several known epilepsy genes. Further functional and genetic studies based on these findings will provide important insights into the complex pathophysiological processes of seizures with and without fever.

Risk factors for anorexia nervosa: A population-based investigation of sex differences in polygenic risk and early life exposures.

OBJECTIVE: To examine sex differences in risk factors for anorexia nervosa (AN). METHOD: This population-based study involved 44,743 individuals (6,239 AN cases including 5,818 females and 421 males, and 38,504 controls including 18,818 females and 19,686 males) born in Denmark between May 1981 and December 2009. Follow-up began on the individual's sixth birthday and ended at AN diagnosis, emigration, death, or December 31, 2016, whichever occurred first. Exposures included socioeconomic status (SES), pregnancy, birth, and early childhood factors based on data from Danish registers, and psychiatric and metabolic polygenic risk scores (PRS) based on genetic data. Hazard ratios were estimated using weighted Cox proportional hazards models stratified by sex (assigned at birth), with AN diagnosis as the outcome. RESULTS: The effects of early life exposures and PRS on AN risk were comparable between females and males. Although we observed some differences in the magnitude and direction of effects, there were no significant interactions between sex and SES, pregnancy, birth, or early childhood exposures. The effects of most PRS on AN risk were highly similar between the sexes. We observed significant sex-specific effects of parental psychiatric history and body mass index PRS, though these effects did not survive corrections for multiple comparisons. CONCLUSIONS: Risk factors for AN are comparable between females and males. Collaboration across countries with large registers is needed to further investigate sex-specific effects of genetic, biological, and environmental exposures on AN risk, including exposures in later childhood and adolescence as well as the additive effects of exposures. PUBLIC SIGNIFICANCE: Sex differences in the prevalence and clinical presentation of AN warrant examination of sex-specific risk factors. This population-based study indicates that the effects of polygenic risk and early life exposures on AN risk are comparable between females and males. Collaboration between countries with large registers is needed to further investigate sex-specific AN risk factors and improve early identification of AN.

Making the Most of Clumping and Thresholding for Polygenic Scores.

Polygenic prediction has the potential to contribute to precision medicine. Clumping and thresholding (C+T) is a widely used method to derive polygenic scores. When using C+T, several p value thresholds are tested to maximize predictive ability of the derived polygenic scores. Along with this p value threshold, we propose to tune three other hyper-parameters for C+T. We implement an efficient way to derive thousands of different C+T scores corresponding to a grid over four hyper-parameters. For example, it takes a few hours to derive 123K different C+T scores for 300K individuals and 1M variants using 16 physical cores. We find that optimizing over these four hyper-parameters improves the predictive performance of C+T in both simulations and real data applications as compared to tuning only the p value threshold. A particularly large increase can be noted when predicting depression status, from an AUC of 0.557 (95% CI: [0.544-0.569]) when tuning only the p value threshold to an AUC of 0.592 (95% CI: [0.580-0.604]) when tuning all four hyper-parameters we propose for C+T. We further propose stacked clumping and thresholding (SCT), a polygenic score that results from stacking all derived C+T scores. Instead of choosing one set of hyper-parameters that maximizes prediction in some training set, SCT learns an optimal linear combination of all C+T scores by using an efficient penalized regression. We apply SCT to eight different case-control diseases in the UK biobank data and find that SCT substantially improves prediction accuracy with an average AUC increase of 0.035 over standard C+T.