Leveraging fine-mapping and multipopulation training data to improve cross-population polygenic risk scores

Asthma-Genomic Advances Toward Risk Prediction

Asthma is a common complex airway disease whose prediction of disease risk and most severe outcomes is crucial in clinical practice for adequate clinical management. This review discusses the latest findings in asthma genomics and current obstacles faced in moving forward to translational medicine. While genome-wide association studies have provided valuable insights into the genetic basis of asthma, there are challenges that must be addressed to improve disease prediction, such as the need for diverse representation, the functional characterization of genetic variants identified, variant selection for genetic testing, and refining prediction models using polygenic risk scores.

Pharmacogenomic scores in psychiatry: systematic review of current evidence

In the past two decades, significant progress has been made in the development of polygenic scores (PGSs). One specific application of PGSs is the development and potential use of pharmacogenomic- scores (PGx-scores) to identify patients who can benefit from a specific medication or are likely to experience side effects. This systematic review comprehensively evaluates published PGx-score studies in psychiatry and provides insights into their potential clinical use and avenues for future development. A systematic literature search was conducted across PubMed, EMBASE, and Web of Science databases until 22 August 2023. This review included fifty-three primary studies, of which the majority (69.8%) were conducted using samples of European ancestry. We found that over 90% of PGx-scores in psychiatry have been developed based on psychiatric and medical diagnoses or trait variants, rather than pharmacogenomic variants. Among these PGx-scores, the polygenic score for schizophrenia (PGSSCZ) has been most extensively studied in relation to its impact on treatment outcomes (32 publications). Twenty (62.5%) of these studies suggest that individuals with higher PGSSCZ have negative outcomes from psychotropic treatment - poorer treatment response, higher rates of treatment resistance, more antipsychotic-induced side effects, or more psychiatric hospitalizations, while the remaining studies did not find significant associations. Although PGx-scores alone accounted for at best 5.6% of the variance in treatment outcomes (in schizophrenia treatment resistance), together with clinical variables they explained up to 13.7% (in bipolar lithium response), suggesting that clinical translation might be achieved by including PGx-scores in multivariable models. In conclusion, our literature review found that there are still very few studies developing PGx-scores using pharmacogenomic variants. Research with larger and diverse populations is required to develop clinically relevant PGx-scores, using biology-informed and multi-phenotypic polygenic scoring approaches, as well as by integrating clinical variables with these scores to facilitate their translation to psychiatric practice.

Sources of gene expression variation in a globally diverse human cohort

Genetic variation that influences gene expression and splicing is a key source of phenotypic diversity1-5. Although invaluable, studies investigating these links in humans have been strongly biased towards participants of European ancestries, which constrains generalizability and hinders evolutionary research. Here to address these limitations, we developed MAGE, an open-access RNA sequencing dataset of lymphoblastoid cell lines from 731 individuals from the 1000 Genomes Project6, spread across 5 continental groups and 26 populations. Most variation in gene expression (92%) and splicing (95%) was distributed within versus between populations, which mirrored the variation in DNA sequence. We mapped associations between genetic variants and expression and splicing of nearby genes (cis-expression quantitative trait loci (eQTLs) and cis-splicing QTLs (sQTLs), respectively). We identified more than 15,000 putatively causal eQTLs and more than 16,000 putatively causal sQTLs that are enriched for relevant epigenomic signatures. These include 1,310 eQTLs and 1,657 sQTLs that are largely private to underrepresented populations. Our data further indicate that the magnitude and direction of causal eQTL effects are highly consistent across populations. Moreover, the apparent 'population-specific' effects observed in previous studies were largely driven by low resolution or additional independent eQTLs of the same genes that were not detected. Together, our study expands our understanding of human gene expression diversity and provides an inclusive resource for studying the evolution and function of human genomes.

Fish oil supplementation modifies the associations between genetically predicted and observed concentrations of blood lipids: a cross-sectional gene-diet interaction study in UK Biobank

BACKGROUND

Dyslipidemia is a well-known risk factor for cardiovascular disease, the leading cause of mortality worldwide. Although habitual intake of fish oil is associated with cardioprotective effects through triglyceride reduction, the interactions of fish oil with the genetic predisposition to dysregulated lipids remain elusive.

OBJECTIVES

We examined whether fish oil supplementation modifies the association between genetically predicted and observed concentrations of total cholesterol, low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, and triglycerides.

METHODS

A total of 441,985 participants with complete genetic and phenotypic data from the UK Biobank were included. Polygenic scores (PGS) of the 4 lipids were calculated in participants of diverse ancestries. For each lipid, multivariable linear regression models were used to assess if fish oil supplementation modified the association between PGS and the observed circulating concentration, with adjustment for relevant covariates.

RESULTS

Fish oil supplementation attenuates the associations between genetically predicted and observed circulating concentrations of total cholesterol, LDL cholesterol, and triglycerides while accentuating the corresponding association for HDL cholesterol among 424,090 participants of European ancestry. Consistent significant findings were obtained using PGS calculated based on multiple genome-wide association studies or alternative PGS methods. For triglycerides, each standard deviation (SD) increment in PGS is associated with 0.254 [95% confidence interval (CI): 0.248, 0.259] SD increase in the observed concentration among European-ancestry participants who reported fish oil usage. In contrast, a stronger association was observed in nonusers (0.267; 95% CI: 0.263, 0.270). Consistently, we showed that fish oil significantly attenuates the association between genetically predicted and observed concentrations of triglycerides in African-ancestry participants.

CONCLUSIONS

Fish oil supplementation attenuates the association between genetically predicted and observed circulating concentrations of total cholesterol, LDL cholesterol, and triglycerides while accentuating the corresponding association for HDL cholesterol in individuals of European ancestry. Further research is needed to understand the clinical implications of these findings.

Real-time dynamic polygenic prediction for streaming data

Polygenic risk scores (PRSs) are promising tools for advancing precision medicine. However, existing PRS construction methods rely on static summary statistics derived from genome-wide association studies (GWASs), which are often updated at lengthy intervals. As genetic data and health outcomes are continuously being generated at an ever-increasing pace, the current PRS training and deployment paradigm is suboptimal in maximizing the prediction accuracy of PRSs for incoming patients in healthcare settings. Here, we introduce real-time PRS-CS (rtPRS-CS), which enables online, dynamic refinement and calibration of PRS as each new sample is collected, without the need to perform intermediate GWASs. Through extensive simulation studies, we evaluate the performance of rtPRS-CS across various genetic architectures and training sample sizes. Leveraging quantitative traits from the Mass General Brigham Biobank and UK Biobank, we show that rtPRS-CS can integrate massive streaming data to enhance PRS prediction over time. We further apply rtPRS-CS to 22 schizophrenia cohorts in 7 Asian regions, demonstrating the clinical utility of rtPRS-CS in dynamically predicting and stratifying disease risk across diverse genetic ancestries.

Autoimmune hypothyroidism GWAS reveals independent autoimmune and thyroid-specific contributions and an inverse relation with cancer risk

The high prevalence of autoimmune hypothyroidism (AIHT) - more than 5% in human populations - provides a unique opportunity to unlock the most complete picture to date of genetic loci that underlie systemic and organ-specific autoimmunity. Using a meta-analysis of 81,718 AIHT cases in FinnGen and the UK Biobank, we dissect associations along axes of thyroid dysfunction and autoimmunity. This largest-to-date scan of hypothyroidism identifies 418 independent associations (p < 5×10- 8), more than half of which have not previously been documented in thyroid disease. In 48 of these, a protein-coding variant is the lead SNP or is highly correlated (r2 > 0.95) with the lead SNP at the locus, including low-frequency coding variants at LAG3, ZAP70, TG, TNFSF11, IRF3, S1PR4, HABP2, ZNF429 as well as established variants at ADCY7, IFIH1 and TYK2. The variants at LAG3 (P67T), ZAP70 (T155M), and TG (Q655X) are highly enriched in Finland and functional experiments in T-cells demonstrate that the ZAP70:T155M allele reduces T-cell activation. By employing a large-scale scan of non-thyroid autoimmunity and a published meta-analysis of TSH levels, we use a Bayesian classifier to dissect the associated loci into distinct groupings and from this estimate, a significant proportion are involved in systemic (i.e., general to multiple autoimmune conditions) autoimmunity (34%) and another subset in thyroid-specific dysfunction (17%). By comparing these association results further to other common disease endpoints, we identify a noteworthy overlap with skin cancer, with 10% of AIHT loci showing a consistent but opposite pattern of association where alleles that increase the risk of hypothyroidism have protective effects for skin cancer. The association results, including genes encoding checkpoint inhibitors and other genes affecting protein levels of PD1, bolster the causal role of natural variation in autoimmunity influencing cancer outcomes.

Influence of Polygenic Background on the Clinical Presentation of Familial Hypercholesterolemia

BACKGROUND

Heterozygous familial hypercholesterolemia (FH) is among the most common genetic conditions worldwide that affects ≈ 1 in 300 individuals. FH is characterized by increased levels of low-density lipoprotein cholesterol (LDL-C) and increased risk of coronary artery disease (CAD), but there is a wide spectrum of severity within the FH population. This variability in expression is incompletely explained by known risk factors. We hypothesized that genome-wide genetic influences, as represented by polygenic risk scores (PRSs) for cardiometabolic traits, would influence the phenotypic severity of FH.

METHODS

We studied individuals with clinically diagnosed FH (n=1123) from the FH Canada National Registry, as well as individuals with genetically identified FH from the UK Biobank (n=723). For all individuals, we used genome-wide gene array data to calculate PRSs for CAD, LDL-C, lipoprotein(a), and other cardiometabolic traits. We compared the distribution of PRSs in individuals with clinically diagnosed FH, genetically diagnosed FH, and non-FH controls and examined the association of the PRSs with the risk of atherosclerotic cardiovascular disease.

RESULTS

Individuals with clinically diagnosed FH had higher levels of LDL-C, and the incidence of atherosclerotic cardiovascular disease was higher in individuals with clinically diagnosed compared with genetically identified FH. Individuals with clinically diagnosed FH displayed enrichment for higher PRSs for CAD, LDL-C, and lipoprotein(a) but not for other cardiometabolic risk factors. The CAD PRS was associated with a risk of atherosclerotic cardiovascular disease among individuals with an FH-causing genetic variant.

CONCLUSIONS

Genetic background, as expressed by genome-wide PRSs for CAD, LDL-C, and lipoprotein(a), influences the phenotypic severity of FH, expanding our understanding of the determinants that contribute to the variable expressivity of FH. A PRS for CAD may aid in risk prediction among individuals with FH.

Calibrated prediction intervals for polygenic scores across diverse contexts

Polygenic scores (PGS) have emerged as the tool of choice for genomic prediction in a wide range of fields. We show that PGS performance varies broadly across contexts and biobanks. Contexts such as age, sex and income can impact PGS accuracy with similar magnitudes as genetic ancestry. Here we introduce an approach (CalPred) that models all contexts jointly to produce prediction intervals that vary across contexts to achieve calibration (include the trait with 90% probability), whereas existing methods are miscalibrated. In analyses of 72 traits across large and diverse biobanks (All of Us and UK Biobank), we find that prediction intervals required adjustment by up to 80% for quantitative traits. For disease traits, PGS-based predictions were miscalibrated across socioeconomic contexts such as annual household income levels, further highlighting the need of accounting for context information in PGS-based prediction across diverse populations.

Developing an optimal stratification model for colorectal cancer screening and reducing racial disparities in multi-center population-based studies

BACKGROUND

Early detection of colorectal neoplasms can reduce the colorectal cancer (CRC) burden by timely intervention for high-risk individuals. However, effective risk prediction models are lacking for personalized CRC early screening in East Asian (EAS) population. We aimed to develop, validate, and optimize a comprehensive risk prediction model across all stages of the dynamic adenoma-carcinoma sequence in EAS population.

METHODS

To develop precision risk-stratification and intervention strategies, we developed three trans-ancestry PRSs targeting colorectal neoplasms: (1) using 148 previously identified CRC risk loci (PRS148); (2) SNPs selection from large-scale meta-analysis data by clumping and thresholding (PRS183); (3) PRS-CSx, a Bayesian approach for genome-wide risk prediction (PRSGenomewide). Then, the performance of each PRS was assessed and validated in two independent cross-sectional screening sets, including 4600 patients with advanced colorectal neoplasm, 4495 patients with non-advanced adenoma, and 21,199 normal individuals from the ZJCRC (Zhejiang colorectal cancer set; EAS) and PLCO (the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial; European, EUR) studies. The optimal PRS was further incorporated with lifestyle factors to stratify individual risk and ultimately tested in the PLCO and UK Biobank prospective cohorts, totaling 350,013 participants.

RESULTS

Three trans-ancestry PRSs achieved moderately improved predictive performance in EAS compared to EUR populations. Remarkably, the PRSs effectively facilitated a thorough risk assessment across all stages of the dynamic adenoma-carcinoma sequence. Among these models, PRS183 demonstrated the optimal discriminatory ability in both EAS and EUR validation datasets, particularly for individuals at risk of colorectal neoplasms. Using two large-scale and independent prospective cohorts, we further confirmed a significant dose-response effect of PRS183 on incident colorectal neoplasms. Incorporating PRS183 with lifestyle factors into a comprehensive strategy improves risk stratification and discriminatory accuracy compared to using PRS or lifestyle factors separately. This comprehensive risk-stratified model shows potential in addressing missed diagnoses in screening tests (best NPV = 0.93), while moderately reducing unnecessary screening (best PPV = 0.32).

CONCLUSIONS

Our comprehensive risk-stratified model in population-based CRC screening trials represents a promising advancement in personalized risk assessment, facilitating tailored CRC screening in the EAS population. This approach enhances the transferability of PRSs across ancestries and thereby helps address health disparity.

Germline prediction of immune checkpoint inhibitor discontinuation for immune-related adverse events

Introduction

Immune checkpoint inhibitors (ICIs) can yield remarkable clinical responses in subsets of patients with solid tumors but can also often lead to immune-related adverse events (irAEs). Predictive features of clinically severe irAEs leading to cessation of ICIs have yet to be established. Using data from 1,327 patients with lung cancer treated with ICIs between 2009 and 2022 at four academic medical centers, we evaluated the association of a germline polygenic risk score for autoimmune disease and discontinuation of ICIs due to irAEs.

Methods

Using Cox proportional hazards model, we assessed the association between a polygenic risk score for autoimmune disease (PRSAD) and cessation of ICI therapy due to irAEs. All models were adjusted for age at diagnosis, sex, lung cancer histology, type of therapy, recruiting center, and the first 5 principal components. To further understand the differential effects of type of therapy and disease stage on the association between PRSAD and cessation of ICI due to irAEs, we conducted stratified logistic regression analysis by type of ICI therapy and disease stage.

Results

We found an association between PRSAD and ICI cessation due to irAEs (HR per SD = 1.18, 95% CI = 1.02 - 1.37, P = 0.03). This association was particularly strong in patients who had ICI cessation due to irAEs within three months of therapy initiation (HR per SD = 1.38, 95% CI = 1.08 - 1.78, P = 0.01). Individuals in the top 20th percentile of PRSAD had 7.2% ICI discontinuation for irAEs by three months, compared to 3.9% discontinuation by three months among patients in the bottom 80th percentile (log-rank P = 0.02). In addition, among patients who received combination PD-1/PD-L1 and CTLA-4 inhibitor therapy, PRSAD had an OR per SD of 1.86 (95% CI = 1.08 - 3.51, P = 0.04).

Conclusions

We demonstrate an association between a polygenic risk score for autoimmune disease and early ICI discontinuation for irAEs, particularly among patients treated with combination ICI therapy. Our results suggest that germline genetics may be used as an adjunctive tool for risk stratification around ICI clinical decision-making in solid tumor oncology.

shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

Optimizing clinico-genomic disease prediction across ancestries: a machine learning strategy with Pareto improvement

BACKGROUND

Accurate prediction of an individual's predisposition to diseases is vital for preventive medicine and early intervention. Various statistical and machine learning models have been developed for disease prediction using clinico-genomic data. However, the accuracy of clinico-genomic prediction of diseases may vary significantly across ancestry groups due to their unequal representation in clinical genomic datasets.

METHODS

We introduced a deep transfer learning approach to improve the performance of clinico-genomic prediction models for data-disadvantaged ancestry groups. We conducted machine learning experiments on multi-ancestral genomic datasets of lung cancer, prostate cancer, and Alzheimer's disease, as well as on synthetic datasets with built-in data inequality and distribution shifts across ancestry groups.

RESULTS

Deep transfer learning significantly improved disease prediction accuracy for data-disadvantaged populations in our multi-ancestral machine learning experiments. In contrast, transfer learning based on linear frameworks did not achieve comparable improvements for these data-disadvantaged populations.

CONCLUSIONS

This study shows that deep transfer learning can enhance fairness in multi-ancestral machine learning by improving prediction accuracy for data-disadvantaged populations without compromising prediction accuracy for other populations, thus providing a Pareto improvement towards equitable clinico-genomic prediction of diseases.

The acceptability and clinical impact of using polygenic scores for risk-estimation of common cancers in primary care: a systematic review

BACKGROUND

Polygenic scores (PGS) have been developed for cancer risk-estimation and show potential as tools to prompt earlier referral for high-risk individuals and aid risk-stratification within cancer screening programmes. This review explores the potential for using PGS to identify individuals at risk of the most common cancers seen in primary care.

METHODS

Two electronic databases were searched up until November 2023 to identify quantitative, qualitative, and mixed methods studies that reported on the acceptability and clinical impact of using PGS to identify individuals at highest risk of breast, prostate, colorectal and lung cancer in primary care. The Mixed Methods Appraisal Tool (MMAT) was used to assess the quality of included studies and a narrative synthesis was used to analyse data.

RESULTS

A total of 190 papers were identified, 18 of which were eligible for inclusion. A cancer risk-assessment tool incorporating PGS was acceptable to the general practice population and their healthcare providers but major challenges to implementation were identified, including lack of evidence for PGS in non-European ancestry and a need for healthcare provider education in genomic medicine. A PGS cancer risk-assessment had relatively limited impact on psychosocial outcomes and health behaviours. However, for prostate cancer, potential applications for its use in primary care were shown.

CONCLUSIONS

Cancer risk assessment incorporating PGS in primary care is acceptable to patients and healthcare providers but there is a paucity of research exploring clinical impact. Few studies were identified, and more research is required before clinical implementation of PGS can be recommended.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization

BACKGROUND

Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups.

METHODS

We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSPT) and ancestry-based continuous shrinkage priors (PRSCSx) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176,988 individuals across 9 diverse cohorts.

RESULTS

Multi-ancestry PRSPT and PRSCSx outperformed ancestry-specific PRSPT and PRSCSx across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, PRSPTmult and PRSCSxmult) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. PRSPTmult demonstrated the strongest association with CHD in individuals of South Asian ancestry and European ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian ancestry (1.56 [1.50-1.61]), Hispanic/Latino ancestry (1.38 [1.24-1.54]), and African ancestry (1.16 [1.11-1.21]). PRSCSxmult showed the strongest associations in South Asian ancestry (2.67 [2.38-3.00]) and European ancestry (1.65 [1.59-1.71]), lower in East Asian ancestry (1.59 [1.54-1.64]), Hispanic/Latino ancestry (1.51 [1.35-1.69]), and the lowest in African ancestry (1.20 [1.15-1.26]).

CONCLUSIONS

The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African ancestry. This highlights the need for larger genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Functional dissection of complex and molecular trait variants at single nucleotide resolution

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries

We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using ∼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.

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.

Haplotype function score improves biological interpretation and cross-ancestry polygenic prediction of human complex traits

We propose a new framework for human genetic association studies: at each locus, a deep learning model (in this study, Sei) is used to calculate the functional genomic activity score for two haplotypes per individual. This score, defined as the Haplotype Function Score (HFS), replaces the original genotype in association studies. Applying the HFS framework to 14 complex traits in the UK Biobank, we identified 3619 independent HFS-trait associations with a significance of p < 5 × 10-8. Fine-mapping revealed 2699 causal associations, corresponding to a median increase of 63 causal findings per trait compared with single-nucleotide polymorphism (SNP)-based analysis. HFS-based enrichment analysis uncovered 727 pathway-trait associations and 153 tissue-trait associations with strong biological interpretability, including 'circadian pathway-chronotype' and 'arachidonic acid-intelligence'. Lastly, we applied least absolute shrinkage and selection operator (LASSO) regression to integrate HFS prediction score with SNP-based polygenic risk scores, which showed an improvement of 16.1-39.8% in cross-ancestry polygenic prediction. We concluded that HFS is a promising strategy for understanding the genetic basis of human complex traits.

Distinct explanations underlie gene-environment interactions in the UK Biobank

The role of gene-environment (GxE) interaction in disease and complex trait architectures is widely hypothesized, but currently unknown. Here, we apply three statistical approaches to quantify and distinguish three different types of GxE interaction for a given trait and E variable. First, we detect locus-specific GxE interaction by testing for genetic correlation r g < 1 across E bins. Second, we detect genome-wide effects of the E variable on genetic variance by leveraging polygenic risk scores (PRS) to test for significant PRSxE in a regression of phenotypes on PRS, E, and PRSxE, together with differences in SNP-heritability across E bins. Third, we detect genome-wide proportional amplification of genetic and environmental effects as a function of the E variable by testing for significant PRSxE with no differences in SNP-heritability across E bins. Simulations show that these approaches achieve high sensitivity and specificity in distinguishing these three GxE scenarios. We applied our framework to 33 UK Biobank traits (25 quantitative traits and 8 diseases; average N = 325 K ) and 10 E variables spanning lifestyle, diet, and other environmental exposures. First, we identified 19 trait-E pairs with r g significantly < 1 (FDR<5%) (average r g = 0.95 ); for example, white blood cell count had r g = 0.95 (s.e. 0.01) between smokers and non-smokers. Second, we identified 28 trait-E pairs with significant PRSxE and significant SNP-heritability differences across E bins; for example, BMI had a significant PRSxE for physical activity (P=4.6e-5) with 5% larger SNP-heritability in the largest versus smallest quintiles of physical activity (P=7e-4). Third, we identified 15 trait-E pairs with significant PRSxE with no SNP-heritability differences across E bins; for example, waist-hip ratio adjusted for BMI had a significant PRSxE effect for time spent watching television (P=5e-3) with no SNP-heritability differences. Across the three scenarios, 8 of the trait-E pairs involved disease traits, whose interpretation is complicated by scale effects. Analyses using biological sex as the E variable produced additional significant findings in each of the three scenarios. Overall, we infer a significant contribution of GxE and GxSex effects to complex trait and disease variance.

Evaluating and improving health equity and fairness of polygenic scores

Polygenic scores (PGSs) are quantitative metrics for predicting phenotypic values, such as human height or disease status. Some PGS methods require only summary statistics of a relevant genome-wide association study (GWAS) for their score. One such method is Lassosum, which inherits the model selection advantages of Lasso to select a meaningful subset of the GWAS single-nucleotide polymorphisms as predictors from their association statistics. However, even efficient scores like Lassosum, when derived from European-based GWASs, are poor predictors of phenotype for subjects of non-European ancestry; that is, they have limited portability to other ancestries. To increase the portability of Lassosum, when GWAS information and estimates of linkage disequilibrium are available for both ancestries, we propose Joint-Lassosum (JLS). In the simulation settings we explore, JLS provides more accurate PGSs compared to other methods, especially when measured in terms of fairness. In analyses of UK Biobank data, JLS was computationally more efficient but slightly less accurate than a Bayesian comparator, SDPRX. Like all PGS methods, JLS requires selection of predictors, which are determined by data-driven tuning parameters. We describe a new approach to selecting tuning parameters and note its relevance for model selection for any PGS. We also draw connections to the literature on algorithmic fairness and discuss how JLS can help mitigate fairness-related harms that might result from the use of PGSs in clinical settings. While no PGS method is likely to be universally portable, due to the diversity of human populations and unequal information content of GWASs for different ancestries, JLS is an effective approach for enhancing portability and reducing predictive bias.

MUSSEL: Enhanced Bayesian polygenic risk prediction leveraging information across multiple ancestry groups

Polygenic risk scores (PRSs) are now showing promising predictive performance on a wide variety of complex traits and diseases, but there exists a substantial performance gap across populations. We propose MUSSEL, a method for ancestry-specific polygenic prediction that borrows information in summary statistics from genome-wide association studies (GWASs) across multiple ancestry groups via Bayesian hierarchical modeling and ensemble learning. In our simulation studies and data analyses across four distinct studies, totaling 5.7 million participants with a substantial ancestral diversity, MUSSEL shows promising performance compared to alternatives. For example, MUSSEL has an average gain in prediction R2 across 11 continuous traits of 40.2% and 49.3% compared to PRS-CSx and CT-SLEB, respectively, in the African ancestry population. The best-performing method, however, varies by GWAS sample size, target ancestry, trait architecture, and linkage disequilibrium reference samples; thus, ultimately a combination of methods may be needed to generate the most robust PRSs across diverse populations.

Integrative polygenic risk score improves the prediction accuracy of complex traits and diseases

Polygenic risk scores (PRSs) are an emerging tool to predict the clinical phenotypes and outcomes of individuals. We propose PRSmix, a framework that leverages the PRS corpus of a target trait to improve prediction accuracy, and PRSmix+, which incorporates genetically correlated traits to better capture the human genetic architecture for 47 and 32 diseases/traits in European and South Asian ancestries, respectively. PRSmix demonstrated a mean prediction accuracy improvement of 1.20-fold (95% confidence interval [CI], [1.10; 1.3]; p = 9.17 × 10-5) and 1.19-fold (95% CI, [1.11; 1.27]; p = 1.92 × 10-6), and PRSmix+ improved the prediction accuracy by 1.72-fold (95% CI, [1.40; 2.04]; p = 7.58 × 10-6) and 1.42-fold (95% CI, [1.25; 1.59]; p = 8.01 × 10-7) in European and South Asian ancestries, respectively. Compared to the previously cross-trait-combination methods with scores from pre-defined correlated traits, we demonstrated that our method improved prediction accuracy for coronary artery disease up to 3.27-fold (95% CI, [2.1; 4.44]; p value after false discovery rate (FDR) correction = 2.6 × 10-4). Our method provides a comprehensive framework to benchmark and leverage the combined power of PRS for maximal performance in a desired target population.

Immunogenetics of lithium response and psychiatric phenotypes in patients with bipolar disorder

The link between bipolar disorder (BP) and immune dysfunction remains controversial. While epidemiological studies have long suggested an association, recent research has found only limited evidence of such a relationship. To clarify this, we performed an exploratory study of the contributions of immune-relevant genetic factors to the response to lithium (Li) treatment and the clinical presentation of BP. First, we assessed the association of a large collection of immune-related genes (4925) with Li response, defined by the Retrospective Assessment of the Lithium Response Phenotype Scale (Alda scale), and clinical characteristics in patients with BP from the International Consortium on Lithium Genetics (ConLi+Gen, N = 2374). Second, we calculated here previously published polygenic scores (PGSs) for immune-related traits and evaluated their associations with Li response and clinical features. Overall, we observed relatively weak associations (p < 1 × 10-4) with BP phenotypes within immune-related genes. Network and functional enrichment analyses of the top findings from the association analyses of Li response variables showed an overrepresentation of pathways participating in cell adhesion and intercellular communication. These appeared to converge on the well-known Li-induced inhibition of GSK-3β. Association analyses of age-at-onset, number of mood episodes, and presence of psychosis, substance abuse and/or suicidal ideation suggested modest contributions of genes such as RTN4, XKR4, NRXN1, NRG1/3 and GRK5 to disease characteristics. PGS analyses returned weak associations (p < 0.05) between inflammation markers and the studied BP phenotypes. Our results suggest a modest relationship between immunity and clinical features in BP. More research is needed to assess the potential therapeutic relevance.

Polygenic Scores and Preclinical Cardiovascular Disease in Individuals With HIV: Insights From the REPRIEVE Trial

BACKGROUND

Coronary artery disease (CAD) is a leading cause of death among the 38.4 million people with HIV globally. The extent to which cardiovascular polygenic risk scores (PRSs) derived in non-HIV populations generalize to people with HIV is not well understood.

METHODS AND RESULTS

PRSs for CAD (GPSMult) and lipid traits were calculated in a global cohort of people with HIV treated with antiretroviral therapy with low-to-moderate atherosclerotic cardiovascular disease risk enrolled in REPRIEVE (Randomized Trial to Prevent Vascular Events in HIV). The PRSs were associated with baseline lipid traits in 4495 genotyped participants, and with subclinical CAD in a subset of 662 who underwent coronary computed tomography angiography. Among participants who underwent coronary computed tomography angiography (mean age, 50.9 [SD, 5.8] years; 16.1% women; 41.8% African, 57.3% European, 1.1% Asian), GPSMult was associated with plaque presence with odds ratio (OR) per SD in GPSMult of 1.42 (95% CI, 1.20-1.68; P=3.8×10-5), stenosis >50% (OR, 2.39 [95% CI, 1.48-3.85]; P=3.4×10-4), and noncalcified/vulnerable plaque (OR, 1.45 [95% CI, 1.23-1.72]; P=9.6×10-6). Effects were consistent in subgroups of age, sex, 10-year atherosclerotic cardiovascular disease risk, ancestry, and CD4 count. Adding GPSMult to established risk factors increased the C-statistic for predicting plaque presence from 0.718 to 0.734 (P=0.02). Furthermore, a PRS for low-density lipoprotein cholesterol was associated with plaque presence with OR of 1.21 (95% CI, 1.01-1.44; P=0.04), and partially calcified plaque with OR of 1.21 (95% CI, 1.01-1.45; P=0.04) per SD.

CONCLUSIONS

Among people with HIV treated with antiretroviral therapy without documented atherosclerotic cardiovascular disease and at low-to-moderate calculated risk in REPRIEVE, an externally developed CAD PRS was predictive of subclinical atherosclerosis. PRS for low-density lipoprotein cholesterol was also associated with subclinical atherosclerosis, supporting a role for low-density lipoprotein cholesterol in HIV-associated CAD.

REGISTRATION

URL: https://www.reprievetrial.org; Unique identifier: NCT02344290.

Bayesian approach to assessing population differences in genetic risk of disease with application to prostate cancer

Population differences in risk of disease are common, but the potential genetic basis for these differences is not well understood. A standard approach is to compare genetic risk across populations by testing for mean differences in polygenic scores, but existing studies that use this approach do not account for statistical noise in effect estimates (i.e., the GWAS betas) that arise due to the finite sample size of GWAS training data. Here, we show using Bayesian polygenic score methods that the level of uncertainty in estimates of genetic risk differences across populations is highly dependent on the GWAS training sample size, the polygenicity (number of causal variants), and genetic distance (FST) between the populations considered. We derive a Wald test for formally assessing the difference in genetic risk across populations, which we show to have calibrated type 1 error rates under a simplified assumption that all SNPs are independent, which we achieve in practise using linkage disequilibrium (LD) pruning. We further provide closed-form expressions for assessing the uncertainty in estimates of relative genetic risk across populations under the special case of an infinitesimal genetic architecture. We suggest that for many complex traits and diseases, particularly those with more polygenic architectures, current GWAS sample sizes are insufficient to detect moderate differences in genetic risk across populations, though more substantial differences in relative genetic risk (relative risk > 1.5) can be detected. We show that conventional approaches that do not account for sampling error from the training sample, such as using a simple t-test, have very high type 1 error rates. When applying our approach to prostate cancer, we demonstrate a higher genetic risk in African Ancestry men, with lower risk in men of European followed by East Asian ancestry.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis

BACKGROUND

'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056).

RESULTS

We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null.

CONCLUSION

Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Genetic and molecular architecture of complex traits

Human genetics has emerged as one of the most dynamic areas of biology, with a broadening societal impact. In this review, we discuss recent achievements, ongoing efforts, and future challenges in the field. Advances in technology, statistical methods, and the growing scale of research efforts have all provided many insights into the processes that have given rise to the current patterns of genetic variation. Vast maps of genetic associations with human traits and diseases have allowed characterization of their genetic architecture. Finally, studies of molecular and cellular effects of genetic variants have provided insights into biological processes underlying disease. Many outstanding questions remain, but the field is well poised for groundbreaking discoveries as it increases the use of genetic data to understand both the history of our species and its applications to improve human health.

Recent advances in polygenic scores: translation, equitability, methods and FAIR tools

Polygenic scores (PGS) can be used for risk stratification by quantifying individuals' genetic predisposition to disease, and many potentially clinically useful applications have been proposed. Here, we review the latest potential benefits of PGS in the clinic and challenges to implementation. PGS could augment risk stratification through combined use with traditional risk factors (demographics, disease-specific risk factors, family history, etc.), to support diagnostic pathways, to predict groups with therapeutic benefits, and to increase the efficiency of clinical trials. However, there exist challenges to maximizing the clinical utility of PGS, including FAIR (Findable, Accessible, Interoperable, and Reusable) use and standardized sharing of the genomic data needed to develop and recalculate PGS, the equitable performance of PGS across populations and ancestries, the generation of robust and reproducible PGS calculations, and the responsible communication and interpretation of results. We outline how these challenges may be overcome analytically and with more diverse data as well as highlight sustained community efforts to achieve equitable, impactful, and responsible use of PGS in healthcare.

Fine-mapping genomic loci refines bipolar disorder risk genes

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Distinguishing different psychiatric disorders using DDx-PRS

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Improving polygenic risk prediction in admixed populations by explicitly modeling ancestral-differential effects via GAUDI

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women's Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

A Stacking Framework for Polygenic Risk Prediction in Admixed Individuals

Polygenic risk scores (PRS) are summaries of an individual's personalized genetic risk for a trait or disease. However, PRS often perform poorly for phenotype prediction when the ancestry of the target population does not match the population in which GWAS effect sizes were estimated. For many populations this can be addressed by performing GWAS in the target population. However, admixed individuals (whose genomes can be traced to multiple ancestral populations) lie on an ancestry continuum and are not easily represented as a discrete population. Here, we propose slaPRS (stacking local ancestry PRS), which incorporates multiple ancestry GWAS to alleviate the ancestry dependence of PRS in admixed samples. slaPRS uses ensemble learning (stacking) to combine local population specific PRS in regions across the genome. We compare slaPRS to single population PRS and a method that combines single population PRS globally. In simulations, slaPRS outperformed existing approaches and reduced the ancestry dependence of PRS in African Americans. In lipid traits from African British individuals (UK Biobank), slaPRS again improved on single population PRS while performing comparably to the globally combined PRS. slaPRS provides a data-driven and flexible framework to incorporate multiple population-specific GWAS and local ancestry in samples of admixed ancestry.

Cross-population applications of genomics to understand the risk of multifactorial traits involving inflammation and immunity

The interplay between genetic and environmental factors plays a significant role in interindividual variation in immune and inflammatory responses. The availability of high-throughput low-cost genotyping and next-generation sequencing has revolutionized our ability to identify human genetic variation and understand how this varies within and between populations, and the relationship with disease. In this review, we explore the potential of genomics for patient benefit, specifically in the diagnosis, prognosis and treatment of inflammatory and immune-related diseases. We summarize the knowledge arising from genetic and functional genomic approaches, and the opportunity for personalized medicine. The review covers applications in infectious diseases, rare immunodeficiencies and autoimmune diseases, illustrating advances in diagnosis and understanding risk including use of polygenic risk scores. We further explore the application for patient stratification and drug target prioritization. The review highlights a key challenge to the field arising from the lack of sufficient representation of genetically diverse populations in genomic studies. This currently limits the clinical utility of genetic-based diagnostic and risk-based applications in non-Caucasian populations. We highlight current genome projects, initiatives and biobanks from diverse populations and how this is being used to improve healthcare globally by improving our understanding of genetic susceptibility to diseases and regional pathogens such as malaria and tuberculosis. Future directions and opportunities for personalized medicine and wider application of genomics in health care are described, for the benefit of individual patients and populations worldwide.

Highly parameterized polygenic scores tend to overfit to population stratification via random effects

Polygenic scores (PGSs), increasingly used in clinical settings, frequently include many genetic variants, with performance typically peaking at thousands of variants. Such highly parameterized PGSs often include variants that do not pass a genome-wide significance threshold. We propose a mathematical perspective that renders the effects of many of these non-significant variants random rather than causal, with the randomness capturing population structure. We devise methods to assess variant effect randomness and population stratification bias. Applying these methods to 141 traits from the UK Biobank, we find that, for many PGSs, the effects of non-significant variants are considerably random, with the extent of randomness associated with the degree of overfitting to population structure of the discovery cohort. Our findings explain why highly parameterized PGSs simultaneously have superior cohort-specific performance and limited generalizability, suggesting the critical need for variant randomness tests in PGS evaluation. Supporting code and a dashboard are available at https://github.com/songlab-cal/StratPGS.

Evaluating the Efficacy of Type 2 Diabetes Polygenic Risk Scores in an Independent European Population

Numerous type 2 diabetes (T2D) polygenic risk scores (PGSs) have been developed to predict individuals' predisposition to the disease. An independent assessment and verification of the best-performing PGS are warranted to allow for a rapid application of developed models. To date, only 3% of T2D PGSs have been evaluated. In this study, we assessed all (n = 102) presently published T2D PGSs in an independent cohort of 3718 individuals, which has not been included in the construction or fine-tuning of any T2D PGS so far. We further chose the best-performing PGS, assessed its performance across major population principal component analysis (PCA) clusters, and compared it with newly developed population-specific T2D PGS. Our findings revealed that 88% of the published PGSs were significantly associated with T2D; however, their performance was lower than what had been previously reported. We found a positive association of PGS improvement over the years (p-value = 8.01 × 10-4 with PGS002771 currently showing the best discriminatory power (area under the receiver operating characteristic (AUROC) = 0.669) and PGS003443 exhibiting the strongest association PGS003443 (odds ratio (OR) = 1.899). Further investigation revealed no difference in PGS performance across major population PCA clusters and when compared with newly developed population-specific PGS. Our findings revealed a positive trend in T2D PGS performance, consistently identifying high-T2D-risk individuals in an independent European population.

PGS-Depot: a comprehensive resource for polygenic scores constructed by summary statistics based methods

Polygenic score (PGS) is an important tool for the genetic prediction of complex traits. However, there are currently no resources providing comprehensive PGSs computed from published summary statistics, and it is difficult to implement and run different PGS methods due to the complexity of their pipelines and parameter settings. To address these issues, we introduce a new resource called PGS-Depot containing the most comprehensive set of publicly available disease-related GWAS summary statistics. PGS-Depot includes 5585 high quality summary statistics (1933 quantitative and 3652 binary trait statistics) curated from 1564 traits in European and East Asian populations. A standardized best-practice pipeline is used to implement 11 summary statistics-based PGS methods, each with different model assumptions and estimation procedures. The prediction performance of each method can be compared for both in- and cross-ancestry populations, and users can also submit their own summary statistics to obtain custom PGS with the available methods. Other features include searching for PGSs by trait name, publication, cohort information, population, or the MeSH ontology tree and searching for trait descriptions with the experimental factor ontology (EFO). All scores, SNP effect sizes and summary statistics can be downloaded via FTP. PGS-Depot is freely available at http://www.pgsdepot.net.

Improving fine-mapping by modeling infinitesimal effects

Fine-mapping aims to identify causal genetic variants for phenotypes. Bayesian fine-mapping algorithms (for example, SuSiE, FINEMAP, ABF and COJO-ABF) are widely used, but assessing posterior probability calibration remains challenging in real data, where model misspecification probably exists, and true causal variants are unknown. We introduce replication failure rate (RFR), a metric to assess fine-mapping consistency by downsampling. SuSiE, FINEMAP and COJO-ABF show high RFR, indicating potential overconfidence in their output. Simulations reveal that nonsparse genetic architecture can lead to miscalibration, while imputation noise, nonuniform distribution of causal variants and quality control filters have minimal impact. Here we present SuSiE-inf and FINEMAP-inf, fine-mapping methods modeling infinitesimal effects alongside fewer larger causal effects. Our methods show improved calibration, RFR and functional enrichment, competitive recall and computational efficiency. Notably, using our methods' posterior effect sizes substantially increases polygenic risk score accuracy over SuSiE and FINEMAP. Our work improves causal variant identification for complex traits, a fundamental goal of human genetics.

Machine Learning Strategies for Improved Phenotype Prediction in Underrepresented Populations

Precision medicine models often perform better for populations of European ancestry due to the over-representation of this group in the genomic datasets and large-scale biobanks from which the models are constructed. As a result, prediction models may misrepresent or provide less accurate treatment recommendations for underrepresented populations, contributing to health disparities. This study introduces an adaptable machine learning toolkit that integrates multiple existing methodologies and novel techniques to enhance the prediction accuracy for underrepresented populations in genomic datasets. By leveraging machine learning techniques, including gradient boosting and automated methods, coupled with novel population-conditional re-sampling techniques, our method significantly improves the phenotypic prediction from single nucleotide polymorphism (SNP) data for diverse populations. We evaluate our approach using the UK Biobank, which is composed primarily of British individuals with European ancestry, and a minority representation of groups with Asian and African ancestry. Performance metrics demonstrate substantial improvements in phenotype prediction for underrepresented groups, achieving prediction accuracy comparable to that of the majority group. This approach represents a significant step towards improving prediction accuracy amidst current dataset diversity challenges. By integrating a tailored pipeline, our approach fosters more equitable validity and utility of statistical genetics methods, paving the way for more inclusive models and outcomes.

Principles and methods for transferring polygenic risk scores across global populations

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Utility of Polygenic Scores for Differentiating Diabetes Diagnosis Among Patients With Atypical Phenotypes of Diabetes

CONTEXT

Misclassification of diabetes type occurs in people with atypical presentations of type 1 diabetes (T1D) or type 2 diabetes (T2D). Although current clinical guidelines suggest clinical variables and treatment response as ways to help differentiate diabetes type, they remain insufficient for people with atypical presentations.

OBJECTIVE

This work aimed to assess the clinical utility of 2 polygenic scores (PGSs) in differentiating between T1D and T2D.

METHODS

Patients diagnosed with diabetes in the UK Biobank were studied (N = 41 787), including 464 (1%) and 15 923 (38%) who met the criteria for classic T1D and T2D, respectively, and 25 400 (61%) atypical diabetes. The validity of 2 published PGSs for T1D (PGST1D) and T2D (PGST2D) in differentiating classic T1D or T2D was assessed using C statistic. The utility of genetic probability for T1D based on PGSs (GenProb-T1D) was evaluated in atypical diabetes patients.

RESULTS

The joint performance of PGST1D and PGST2D for differentiating classic T1D or T2D was outstanding (C statistic = 0.91), significantly higher than that of PGST1D alone (0.88) and PGST2D alone (0.70), both P less than .001. Using an optimal cutoff of GenProb-T1D, 23% of patients with atypical diabetes had a higher probability of T1D and its validity was independently supported by clinical presentations that are characteristic of T1D.

CONCLUSION

PGST1D and PGST2D can be used to discriminate classic T1D and T2D and have potential clinical utility for differentiating these 2 types of diseases among patients with atypical diabetes.

Pervasive correlations between causal disease effects of proximal SNPs vary with functional annotations and implicate stabilizing selection

The genetic architecture of human diseases and complex traits has been extensively studied, but little is known about the relationship of causal disease effect sizes between proximal SNPs, which have largely been assumed to be independent. We introduce a new method, LD SNP-pair effect correlation regression (LDSPEC), to estimate the correlation of causal disease effect sizes of derived alleles between proximal SNPs, depending on their allele frequencies, LD, and functional annotations; LDSPEC produced robust estimates in simulations across various genetic architectures. We applied LDSPEC to 70 diseases and complex traits from the UK Biobank (average N=306K), meta-analyzing results across diseases/traits. We detected significantly nonzero effect correlations for proximal SNP pairs (e.g., -0.37±0.09 for low-frequency positive-LD 0-100bp SNP pairs) that decayed with distance (e.g., -0.07±0.01 for low-frequency positive-LD 1-10kb), varied with allele frequency (e.g., -0.15±0.04 for common positive-LD 0-100bp), and varied with LD between SNPs (e.g., +0.12±0.05 for common negative-LD 0-100bp) (because we consider derived alleles, positive-LD and negative-LD SNP pairs may yield very different results). We further determined that SNP pairs with shared functions had stronger effect correlations that spanned longer genomic distances, e.g., -0.37±0.08 for low-frequency positive-LD same-gene promoter SNP pairs (average genomic distance of 47kb (due to alternative splicing)) and -0.32±0.04 for low-frequency positive-LD H3K27ac 0-1kb SNP pairs. Consequently, SNP-heritability estimates were substantially smaller than estimates of the sum of causal effect size variances across all SNPs (ratio of 0.87±0.02 across diseases/traits), particularly for certain functional annotations (e.g., 0.78±0.01 for common Super enhancer SNPs)-even though these quantities are widely assumed to be equal. We recapitulated our findings via forward simulations with an evolutionary model involving stabilizing selection, implicating the action of linkage masking, whereby haplotypes containing linked SNPs with opposite effects on disease have reduced effects on fitness and escape negative selection.

Pervasive correlations between causal disease effects of proximal SNPs vary with functional annotations and implicate stabilizing selection

The genetic architecture of human diseases and complex traits has been extensively studied, but little is known about the relationship of causal disease effect sizes between proximal SNPs, which have largely been assumed to be independent. We introduce a new method, LD SNP-pair effect correlation regression (LDSPEC), to estimate the correlation of causal disease effect sizes of derived alleles between proximal SNPs, depending on their allele frequencies, LD, and functional annotations; LDSPEC produced robust estimates in simulations across various genetic architectures. We applied LDSPEC to 70 diseases and complex traits from the UK Biobank (average N=306K), meta-analyzing results across diseases/traits. We detected significantly nonzero effect correlations for proximal SNP pairs (e.g., -0.37±0.09 for low-frequency positive-LD 0-100bp SNP pairs) that decayed with distance (e.g., -0.07±0.01 for low-frequency positive-LD 1-10kb), varied with allele frequency (e.g., -0.15±0.04 for common positive-LD 0-100bp), and varied with LD between SNPs (e.g., +0.12±0.05 for common negative-LD 0-100bp) (because we consider derived alleles, positive-LD and negative-LD SNP pairs may yield very different results). We further determined that SNP pairs with shared functions had stronger effect correlations that spanned longer genomic distances, e.g., -0.37±0.08 for low-frequency positive-LD same-gene promoter SNP pairs (average genomic distance of 47kb (due to alternative splicing)) and -0.32±0.04 for low-frequency positive-LD H3K27ac 0-1kb SNP pairs. Consequently, SNP-heritability estimates were substantially smaller than estimates of the sum of causal effect size variances across all SNPs (ratio of 0.87±0.02 across diseases/traits), particularly for certain functional annotations (e.g., 0.78±0.01 for common Super enhancer SNPs)-even though these quantities are widely assumed to be equal. We recapitulated our findings via forward simulations with an evolutionary model involving stabilizing selection, implicating the action of linkage masking, whereby haplotypes containing linked SNPs with opposite effects on disease have reduced effects on fitness and escape negative selection.

Phenotype integration improves power and preserves specificity in biobank-based genetic studies of major depressive disorder

Biobanks often contain several phenotypes relevant to diseases such as major depressive disorder (MDD), with partly distinct genetic architectures. Researchers face complex tradeoffs between shallow (large sample size, low specificity/sensitivity) and deep (small sample size, high specificity/sensitivity) phenotypes, and the optimal choices are often unclear. Here we propose to integrate these phenotypes to combine the benefits of each. We use phenotype imputation to integrate information across hundreds of MDD-relevant phenotypes, which significantly increases genome-wide association study (GWAS) power and polygenic risk score (PRS) prediction accuracy of the deepest available MDD phenotype in UK Biobank, LifetimeMDD. We demonstrate that imputation preserves specificity in its genetic architecture using a novel PRS-based pleiotropy metric. We further find that integration via summary statistics also enhances GWAS power and PRS predictions, but can introduce nonspecific genetic effects depending on input. Our work provides a simple and scalable approach to improve genetic studies in large biobanks by integrating shallow and deep phenotypes.

Can polygenic risk scores help explain disease prevalence differences around the world? A worldwide investigation

Complex disorders are caused by a combination of genetic, environmental and lifestyle factors, and their prevalence can vary greatly across different populations. The extent to which genetic risk, as identified by Genome Wide Association Study (GWAS), correlates to disease prevalence in different populations has not been investigated systematically. Here, we studied 14 different complex disorders and explored whether polygenic risk scores (PRS) based on current GWAS correlate to disease prevalence within Europe and around the world. A clear variation in GWAS-based genetic risk was observed based on ancestry and we identified populations that have a higher genetic liability for developing certain disorders. We found that for four out of the 14 studied disorders, PRS significantly correlates to disease prevalence within Europe. We also found significant correlations between worldwide disease prevalence and PRS for eight of the studied disorders with Multiple Sclerosis genetic risk having the highest correlation to disease prevalence. Based on current GWAS results, the across population differences in genetic risk for certain disorders can potentially be used to understand differences in disease prevalence and identify populations with the highest genetic liability. The study highlights both the limitations of PRS based on current GWAS but also the fact that in some cases, PRS may already have high predictive power. This could be due to the genetic architecture of specific disorders or increased GWAS power in some cases.

Ancestry-specific polygenic risk scores are risk enhancers for clinical cardiovascular disease assessments

Clinical implementation of new prediction models requires evaluation of their utility in a broad range of intended use populations. Here we develop and validate ancestry-specific Polygenic Risk Scores (PRSs) for Coronary Artery Disease (CAD) using 29,389 individuals from diverse cohorts and genetic ancestry groups. The CAD PRSs outperform published scores with an average Odds Ratio per Standard Deviation of 1.57 (SD = 0.14) and identify between 12% and 24% of individuals with high genetic risk. Using this risk factor to reclassify borderline or intermediate 10 year Atherosclerotic Cardiovascular Disease (ASCVD) risk improves assessments for both CAD (Net Reclassification Improvement (NRI) = 13.14% (95% CI 9.23-17.06%)) and ASCVD (NRI = 10.70 (95% CI 7.35-14.05)) in an independent cohort of 9,691 individuals. Our analyses demonstrate that using PRSs as Risk Enhancers improves ASCVD risk assessments outlining an approach for guiding ASCVD prevention with genetic information.

Power of inclusion: Enhancing polygenic prediction with admixed individuals

Admixed individuals offer unique opportunities for addressing limited transferability in polygenic scores (PGSs), given the substantial trans-ancestry genetic correlation in many complex traits. However, they are rarely considered in PGS training, given the challenges in representing ancestry-matched linkage-disequilibrium reference panels for admixed individuals. Here we present inclusive PGS (iPGS), which captures ancestry-shared genetic effects by finding the exact solution for penalized regression on individual-level data and is thus naturally applicable to admixed individuals. We validate our approach in a simulation study across 33 configurations with varying heritability, polygenicity, and ancestry composition in the training set. When iPGS is applied to n = 237,055 ancestry-diverse individuals in the UK Biobank, it shows the greatest improvements in Africans by 48.9% on average across 60 quantitative traits and up to 50-fold improvements for some traits (neutrophil count, R2 = 0.058) over the baseline model trained on the same number of European individuals. When we allowed iPGS to use n = 284,661 individuals, we observed an average improvement of 60.8% for African, 11.6% for South Asian, 7.3% for non-British White, 4.8% for White British, and 17.8% for the other individuals. We further developed iPGS+refit to jointly model the ancestry-shared and -dependent genetic effects when heterogeneous genetic associations were present. For neutrophil count, for example, iPGS+refit showed the highest predictive performance in the African group (R2 = 0.115), which exceeds the best predictive performance for the White British group (R2 = 0.090 in the iPGS model), even though only 1.49% of individuals used in the iPGS training are of African ancestry. Our results indicate the power of including diverse individuals for developing more equitable PGS models.

The GenoVA study: Equitable implementation of a pragmatic randomized trial of polygenic-risk scoring in primary care

Polygenic risk scores (PRSs) hold promise for disease risk assessment and prevention. The Genomic Medicine at Veterans Affairs (GenoVA) Study is addressing three main challenges to the clinical implementation of PRSs in preventive care: defining and determining their clinical utility, implementing them in time-constrained primary care settings, and countering their potential to exacerbate healthcare disparities. The study processes used to test patients, report their PRS results to them and their primary care providers (PCPs), and promote the use of those results in clinical decision-making are modeled on common practices in primary care. The following diseases were chosen for their prevalence and familiarity to PCPs: coronary artery disease; type 2 diabetes; atrial fibrillation; and breast, colorectal, and prostate cancers. A randomized clinical trial (RCT) design and primary outcome of time-to-new-diagnosis of a target disease bring methodological rigor to the question of the clinical utility of PRS implementation. The study's pragmatic RCT design enhances its relevance to how PRS might reasonably be implemented in primary care. Steps the study has taken to promote health equity include the thoughtful handling of genetic ancestry in PRS construction and reporting and enhanced recruitment strategies to address underrepresentation in research participation. To date, enhanced recruitment efforts have been both necessary and successful: participants of underrepresented race and ethnicity groups have been less likely to enroll in the study than expected but ultimately achieved proportional representation through targeted efforts. The GenoVA Study experience to date offers insights for evaluating the clinical utility of equitable PRS implementation in adult primary care.

Polygenic risk scores for disease risk prediction in Africa: current challenges and future directions

Early identification of genetic risk factors for complex diseases can enable timely interventions and prevent serious outcomes, including mortality. While the genetics underlying many Mendelian diseases have been elucidated, it is harder to predict risk for complex diseases arising from the combined effects of many genetic variants with smaller individual effects on disease aetiology. Polygenic risk scores (PRS), which combine multiple contributing variants to predict disease risk, have the potential to influence the implementation for precision medicine. However, the majority of existing PRS were developed from European data with limited transferability to African populations. Notably, African populations have diverse genetic backgrounds, and a genomic architecture with smaller haplotype blocks compared to European genomes. Subsequently, growing evidence shows that using large-scale African ancestry cohorts as discovery for PRS development may generate more generalizable findings. Here, we (1) discuss the factors contributing to the poor transferability of PRS in African populations, (2) showcase the novel Africa genomic datasets for PRS development, (3) explore the potential clinical utility of PRS in African populations, and (4) provide insight into the future of PRS in Africa.

Machine Learning Strategies for Improved Phenotype Prediction in Underrepresented Populations

Precision medicine models often perform better for populations of European ancestry due to the over-representation of this group in the genomic datasets and large-scale biobanks from which the models are constructed. As a result, prediction models may misrepresent or provide less accurate treatment recommendations for underrepresented populations, contributing to health disparities. This study introduces an adaptable machine learning toolkit that integrates multiple existing methodologies and novel techniques to enhance the prediction accuracy for underrepresented populations in genomic datasets. By leveraging machine learning techniques, including gradient boosting and automated methods, coupled with novel population-conditional re-sampling techniques, our method significantly improves the phenotypic prediction from single nucleotide polymorphism (SNP) data for diverse populations. We evaluate our approach using the UK Biobank, which is composed primarily of British individuals with European ancestry, and a minority representation of groups with Asian and African ancestry. Performance metrics demonstrate substantial improvements in phenotype prediction for underrepresented groups, achieving prediction accuracy comparable to that of the majority group. This approach represents a significant step towards improving prediction accuracy amidst current dataset diversity challenges. By integrating a tailored pipeline, our approach fosters more equitable validity and utility of statistical genetics methods, paving the way for more inclusive models and outcomes.