Integration of questionnaire-based risk factors improves polygenic risk scores for human coronary heart disease and type 2 diabetes

A unified framework for estimating country-specific cumulative incidence for 18 diseases stratified by polygenic risk

Polygenic scores (PGSs) offer the ability to predict genetic risk for complex diseases across the life course; a key benefit over short-term prediction models. To produce risk estimates relevant to clinical and public health decision-making, it is important to account for varying effects due to age and sex. Here, we develop a novel framework to estimate country-, age-, and sex-specific estimates of cumulative incidence stratified by PGS for 18 high-burden diseases. We integrate PGS associations from seven studies in four countries (N = 1,197,129) with disease incidences from the Global Burden of Disease. PGS has a significant sex-specific effect for asthma, hip osteoarthritis, gout, coronary heart disease and type 2 diabetes (T2D), with all but T2D exhibiting a larger effect in men. PGS has a larger effect in younger individuals for 13 diseases, with effects decreasing linearly with age. We show for breast cancer that, relative to individuals in the bottom 20% of polygenic risk, the top 5% attain an absolute risk for screening eligibility 16.3 years earlier. Our framework increases the generalizability of results from biobank studies and the accuracy of absolute risk estimates by appropriately accounting for age- and sex-specific PGS effects. Our results highlight the potential of PGS as a screening tool which may assist in the early prevention of common diseases.

Polygenic risk score predicting susceptibility and outcome of benign prostatic hyperplasia in the Han Chinese

BACKGROUND

Given the high prevalence of BPH among elderly men, pinpointing those at elevated risk can aid in early intervention and effective management. This study aimed to explore that polygenic risk score (PRS) is effective in predicting benign prostatic hyperplasia (BPH) incidence, prognosis and risk of operation in Han Chinese.

METHODS

A retrospective cohort study included 12,474 male participants (6,237 with BPH and 6,237 non-BPH controls) from the Taiwan Precision Medicine Initiative (TPMI). Genotyping was performed using the Affymetrix Genome-Wide TWB 2.0 SNP Array. PRS was calculated using PGS001865, comprising 1,712 single nucleotide polymorphisms. Logistic regression models assessed the association between PRS and BPH incidence, adjusting for age and prostate-specific antigen (PSA) levels. The study also examined the relationship between PSA, prostate volume, and response to 5-α-reductase inhibitor (5ARI) treatment, as well as the association between PRS and the risk of TURP.

RESULTS

Individuals in the highest PRS quartile (Q4) had a significantly higher risk of BPH compared to the lowest quartile (Q1) (OR = 1.51, 95% CI = 1.274-1.783, p < 0.0001), after adjusting for PSA level. The Q4 group exhibited larger prostate volumes and a smaller volume reduction after 5ARI treatment. The Q1 group had a lower cumulative TURP probability at 3, 5, and 10 years compared to the Q4 group. PRS Q4 was an independent risk factor for TURP.

CONCLUSIONS

In this Han Chinese cohort, higher PRS was associated with an increased susceptibility to BPH, larger prostate volumes, poorer response to 5ARI treatment, and a higher risk of TURP. Larger prospective studies with longer follow-up are warranted to further validate these findings.

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.

High-Resolution Genotyping of Formalin-Fixed Tissue Accurately Estimates Polygenic Risk Scores in Human Diseases

Formalin-fixed paraffin-embedded (FFPE) tissues stored in biobanks and pathology archives are a vast but underutilized source for molecular studies on different diseases. Beyond being the "gold standard" for preservation of diagnostic human tissues, FFPE samples retain similar genetic information as matching blood samples, which could make FFPE samples an ideal resource for genomic analysis. However, research on this resource has been hindered by the perception that DNA extracted from FFPE samples is of poor quality. Here, we show that germline disease-predisposing variants and polygenic risk scores (PRS) can be identified from FFPE normal tissue (FFPE-NT) DNA with high accuracy. We optimized the performance of FFPE-NT DNA on a genome-wide array containing 657,675 variants. Via a series of testing and validation phases, we established a protocol for FFPE-NT genotyping with results comparable with blood genotyping. The median call rate of FFPE-NT samples in the validation phase was 99.85% (range 98.26%-99.94%) and median concordance with matching blood samples was 99.79% (range 98.85%-99.9%). We also demonstrated that a rare pathogenic PALB2 genetic variant predisposing to cancer can be correctly identified in FFPE-NT samples. We further imputed the FFPE-NT genotype data and calculated the FFPE-NT genome-wide PRS in 3 diseases and 4 disease risk variables. In all cases, FFPE-NT and matching blood PRS were highly concordant (all Pearson's r > 0.95). The ability to precisely genotype FFPE-NT on a genome-wide array enables translational genomics applications of archived FFPE-NT samples with the possibility to link to corresponding phenotypes and longitudinal health data.

Incorporating polygenic risk into the Leicester Risk Assessment score for 10-year risk prediction of type 2 diabetes

AIMS

We evaluated whether incorporating information on ethnic background and polygenic risk enhanced the Leicester Risk Assessment (LRA) score for predicting 10-year risk of type 2 diabetes.

METHODS

The sample included 202,529 UK Biobank participants aged 40-69 years. We computed the LRA score, and developed two new risk scores using training data (80% sample): LRArev, which incorporated additional information on ethnic background, and LRAprs, which incorporated polygenic risk for type 2 diabetes. We assessed discriminative and reclassification performance in a test set (20% sample). Type 2 diabetes was ascertained using primary care, hospital inpatient and death registry records.

RESULTS

Over 10 years, 7,476 participants developed type 2 diabetes. The Harrell's C indexes were 0.796 (95% Confidence Interval [CI] 0.785, 0.806), 0.802 (95% CI 0.792, 0.813), and 0.829 (95% CI 0.820, 0.839) for the LRA, LRArev and LRAprs scores, respectively. The LRAprs score significantly improved the overall reclassification compared to the LRA (net reclassification index [NRI] = 0.033, 95% CI 0.015, 0.049) and LRArev (NRI = 0.040, 95% CI 0.024, 0.055) scores.

CONCLUSIONS

Polygenic risk moderately improved the performance of the existing LRA score for 10-year risk prediction of type 2 diabetes.

Integration of risk factor polygenic risk score with disease polygenic risk score for disease prediction

Polygenic risk score (PRS) is useful for capturing an individual's genetic susceptibility. However, previous studies have not fully exploited the potential of the risk factor PRS (RFPRS) for disease prediction. We explored the potential of integrating disease-related RFPRSs with disease PRS to enhance disease prediction performance. We constructed 112 RFPRSs and analyzed the association of RFPRSs with diseases to identify disease-related RFPRSs in 700 diseases, using the UK Biobank dataset. We uncovered 6157 statistically significant associations between 247 diseases and 109 RFPRSs. We estimated the disease PRSs of 70 diseases that exhibited statistically significant heritability, to generate RFDiseasemetaPRS-a combined PRS integrating RFPRSs and disease PRS-and compare the prediction performance metrics between RFDiseasemetaPRS and disease PRS. RFDiseasemetaPRS showed better performance for Nagelkerke's pseudo-R2, odds ratio (OR) per 1 SD, net reclassification improvement (NRI) values and difference of R2 considered by variance of R2 in 31 out of 70 diseases. Additionally, we assessed risk classification between two models by examining OR between the top 10% and remaining 90% individuals for the 31 diseases; RFDiseasemetaPRS exhibited better R2, NRI and OR than disease PRS. These findings highlight the importance of utilizing RFDiseasemetaPRS, which can provide personalized healthcare and tailored prevention strategies.

A linear weighted combination of polygenic scores for a broad range of traits improves prediction of coronary heart disease

Polygenic scores (PGS) for coronary heart disease (CHD) are constructed using GWAS summary statistics for CHD. However, pleiotropy is pervasive in biology and disease-associated variants often share etiologic pathways with multiple traits. Therefore, incorporating GWAS summary statistics of additional traits could improve the performance of PGS for CHD. Using lasso regression models, we developed two multi-PGS for CHD: 1) multiPGSCHD, utilizing GWAS summary statistics for CHD, its risk factors, and other ASCVD as training data and the UK Biobank for tuning, and 2) extendedPGSCHD, using existing PGS for a broader range of traits in the PGS Catalog as training data and the Atherosclerosis Risk in Communities Study (ARIC) cohort for tuning. We evaluated the performance of multiPGSCHD and extendedPGSCHD in the Mayo Clinic Biobank, an independent cohort of 43,578 adults of European ancestry which included 4,479 CHD cases and 39,099 controls. In the Mayo Clinic Biobank, a 1 SD increase in multiPGSCHD and extendedPGSCHD was associated with a 1.66-fold (95% CI: 1.60-1.71) and 1.70-fold (95% CI: 1.64-1.76) increased odds of CHD, respectively, in models that included age, sex, and 10 PCs, whereas an already published PGS for CHD (CHD_PRSCS) increased the odds by 1.50 (95% CI: 1.45-1.56). In the highest deciles of extendedPGSCHD, multiPGSCHD, and CHD_PRSCS, 18.4%, 17.5%, and 16.3% of patients had CHD, respectively.

Integration of Biomarker Polygenic Risk Score Improves Prediction of Coronary Heart Disease

There are several established biomarkers for coronary heart disease (CHD), including blood pressure, cholesterol, and lipoproteins. It is of high interest to determine how a combined polygenic risk score (PRS) of CHD-associated biomarkers (BioPRS) can further improve genetic prediction of CHD. We developed CHDBioPRS, combining BioPRS with PRS of CHD in the UK Biobank and tested it on FinnGen. We found that BioPRS was clearly predictive of CHD and that CHDBioPRS improved the standard CHD PRS. The largest effect was observed with early onset cases in FinnGen, with HRs above 2 per standard deviation of CHDBioPRS.

Minimal improvement in coronary artery disease risk prediction in Chinese population using polygenic risk scores: evidence from the China Kadoorie Biobank

BACKGROUND

Several studies have reported that polygenic risk scores (PRSs) can enhance risk prediction of coronary artery disease (CAD) in European populations. However, research on this topic is far from sufficient in non-European countries, including China. We aimed to evaluate the potential of PRS for predicting CAD for primary prevention in the Chinese population.

METHODS

Participants with genome-wide genotypic data from the China Kadoorie Biobank were divided into training ( n = 28,490) and testing sets ( n = 72,150). Ten previously developed PRSs were evaluated, and new ones were developed using clumping and thresholding or LDpred method. The PRS showing the strongest association with CAD in the training set was selected to further evaluate its effects on improving the traditional CAD risk-prediction model in the testing set. Genetic risk was computed by summing the product of the weights and allele dosages across genome-wide single-nucleotide polymorphisms. Prediction of the 10-year first CAD events was assessed using hazard ratios (HRs) and measures of model discrimination, calibration, and net reclassification improvement (NRI). Hard CAD (nonfatal I21-I23 and fatal I20-I25) and soft CAD (all fatal or nonfatal I20-I25) were analyzed separately.

RESULTS

In the testing set, 1214 hard and 7201 soft CAD cases were documented during a mean follow-up of 11.2 years. The HR per standard deviation of the optimal PRS was 1.26 (95% CI:1.19-1.33) for hard CAD. Based on a traditional CAD risk prediction model containing only non-laboratory-based information, the addition of PRS for hard CAD increased Harrell's C index by 0.001 (-0.001 to 0.003) in women and 0.003 (0.001 to 0.005) in men. Among the different high-risk thresholds ranging from 1% to 10%, the highest categorical NRI was 3.2% (95% CI: 0.4-6.0%) at a high-risk threshold of 10.0% in women. The association of the PRS with soft CAD was much weaker than with hard CAD, leading to minimal or no improvement in the soft CAD model.

CONCLUSIONS

In this Chinese population sample, the current PRSs minimally changed risk discrimination and offered little improvement in risk stratification for soft CAD. Therefore, this may not be suitable for promoting genetic screening in the general Chinese population to improve CAD risk prediction.

Genome-wide polygenic risk score for type 2 diabetes in Indian population

Genome-wide polygenic risk scores (PRS) for lifestyle disorders, like Type 2 Diabetes (T2D), are useful in identifying at-risk individuals early on in life, and to guide them towards healthier lifestyles. The current study was aimed at developing PRS for the Indian population using imputed genotype data from UK Biobank and testing the developed PRS on data from GenomegaDB of Indians living in India. 959 T2D cases and 2,818 controls were selected from Indian participants of UK Biobank to develop the PRS. Summary statistics available for South Asians, from the DIAMANTE consortium, were used to weigh genetic variants. LDpred2 algorithm was used to adjust the effect of linkage disequilibrium among the variants. The association of PRS with T2D, after adjusting for age, sex and top ten genetic principal components, was found to be very significant (AUC = 0.7953, OR = 2.9856 [95% CI: 2.7044-3.2961]). When participants were divided into four PRS quartile groups, the odds of developing T2D increased sequentially with the higher PRS groups. The highest PRS group (top 25%) showed 5.79 fold increased risk compared to the rest of the participants (75%). The PRS derived using the same set of variants was found to be significantly associated with T2D in the test dataset of 445 Indians (AUC = 0.7781, OR = 1.6656 [95%CI = 0.6127-4.5278]). Our study demonstrates a framework to derive Indian-specific PRS for T2D. The accuracy of the derived PRS shows it's potential to be used as a prognostic metric to stratify individuals, and to recommend personalized preventive strategies.

A multi-ancestry polygenic risk score improves risk prediction for coronary artery disease

Identification of individuals at highest risk of coronary artery disease (CAD)-ideally before onset-remains an important public health need. Prior studies have developed genome-wide polygenic scores to enable risk stratification, reflecting the substantial inherited component to CAD risk. Here we develop a new and significantly improved polygenic score for CAD, termed GPSMult, that incorporates genome-wide association data across five ancestries for CAD (>269,000 cases and >1,178,000 controls) and ten CAD risk factors. GPSMult strongly associated with prevalent CAD (odds ratio per standard deviation 2.14, 95% confidence interval 2.10-2.19, P < 0.001) in UK Biobank participants of European ancestry, identifying 20.0% of the population with 3-fold increased risk and conversely 13.9% with 3-fold decreased risk as compared with those in the middle quintile. GPSMult was also associated with incident CAD events (hazard ratio per standard deviation 1.73, 95% confidence interval 1.70-1.76, P < 0.001), identifying 3% of healthy individuals with risk of future CAD events equivalent to those with existing disease and significantly improving risk discrimination and reclassification. Across multiethnic, external validation datasets inclusive of 33,096, 124,467, 16,433 and 16,874 participants of African, European, Hispanic and South Asian ancestry, respectively, GPSMult demonstrated increased strength of associations across all ancestries and outperformed all available previously published CAD polygenic scores. These data contribute a new GPSMult for CAD to the field and provide a generalizable framework for how large-scale integration of genetic association data for CAD and related traits from diverse populations can meaningfully improve polygenic risk prediction.

Genomic Innovation in Early Life Cardiovascular Disease Prevention and Treatment

Cardiovascular disease (CVD) is a leading cause of morbidity and mortality globally. Although CVD events do not typically manifest until older adulthood, CVD develops gradually across the life-course, beginning with the elevation of risk factors observed as early as childhood or adolescence and the emergence of subclinical disease that can occur in young adulthood or midlife. Genomic background, which is determined at zygote formation, is among the earliest risk factors for CVD. With major advances in molecular technology, including the emergence of gene-editing techniques, along with deep whole-genome sequencing and high-throughput array-based genotyping, scientists now have the opportunity to not only discover genomic mechanisms underlying CVD but use this knowledge for the life-course prevention and treatment of these conditions. The current review focuses on innovations in the field of genomics and their applications to monogenic and polygenic CVD prevention and treatment. With respect to monogenic CVD, we discuss how the emergence of whole-genome sequencing technology has accelerated the discovery of disease-causing variants, allowing comprehensive screening and early, aggressive CVD mitigation strategies in patients and their families. We further describe advances in gene editing technology, which might soon make possible cures for CVD conditions once thought untreatable. In relation to polygenic CVD, we focus on recent innovations that leverage findings of genome-wide association studies to identify druggable gene targets and develop predictive genomic models of disease, which are already facilitating breakthroughs in the life-course treatment and prevention of CVD. Gaps in current research and future directions of genomics studies are also discussed. In aggregate, we hope to underline the value of leveraging genomics and broader multiomics information for characterizing CVD conditions, work which promises to expand precision approaches for the life-course prevention and treatment of CVD.

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

Polygenic risk scores (PRS) are an emerging tool to predict the clinical phenotypes and outcomes of individuals. Validation and transferability of existing PRS across independent datasets and diverse ancestries are limited, which hinders the practical utility and exacerbates health disparities. We propose PRSmix, a framework that evaluates and 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. We applied PRSmix to 47 and 32 diseases/traits in European and South Asian ancestries, respectively. PRSmix demonstrated a mean prediction accuracy improvement of 1.20-fold (95% CI: [1.10; 1.3]; P-value = 9.17 × 10-5) and 1.19-fold (95% CI: [1.11; 1.27]; P-value = 1.92 × 10-6), and PRSmix+ improved the prediction accuracy by 1.72-fold (95% CI: [1.40; 2.04]; P-value = 7.58 × 10-6) and 1.42-fold (95% CI: [1.25; 1.59]; P-value = 8.01 × 10-7) in European and South Asian ancestries, respectively. Compared to the previously established cross-trait-combination method with scores from pre-defined correlated traits, we demonstrated that our method can improve prediction accuracy for coronary artery disease up to 3.27-fold (95% CI: [2.1; 4.44]; P-value after 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.

Integrating Common Risk Factors with Polygenic Scores Improves the Prediction of Type 2 Diabetes

We tested associations between 13 established genetic variants and type 2 diabetes (T2D) in 1371 study participants from the Volga-Ural region of the Eurasian continent, and evaluated the predictive ability of the model containing polygenic scores for the variants associated with T2D in our dataset, alone and in combination with other risk factors such as age and sex. Using logistic regression analysis, we found associations with T2D for the CCL20 rs6749704 (OR = 1.68, P_FDR = 3.40 × 10^-5), CCR5 rs333 (OR = 1.99, P_FDR = 0.033), ADIPOQ rs17366743 (OR = 3.17, P_FDR = 2.64 × 10^-4), TCF7L2 rs114758349 (OR = 1.77, P_FDR = 9.37 × 10^-5), and CCL2 rs1024611 (OR = 1.38, P_FDR = 0.033) polymorphisms. We showed that the most informative prognostic model included weighted polygenic scores for these five loci, and non-genetic factors such as age and sex (AUC 85.8%, 95%CI 83.7-87.8%). Compared to the model containing only non-genetic parameters, adding the polygenic score for the five T2D-associated loci showed improved net reclassification (NRI = 37.62%, 1.39 × 10^-6). Inclusion of all 13 tested SNPs to the model with age and sex did not improve the predictive ability compared to the model containing five T2D-associated variants (NRI = -17.86, p = 0.093). The five variants associated with T2D in people from the Volga-Ural region are linked to inflammation (CCR5, CCL2, CCL20) and glucose metabolism regulation (TCF7L, ADIPOQ2). Further studies in independent groups of T2D patients should validate the prognostic value of the model and elucidate the molecular mechanisms of the disease development.

Clinicians' Perceptions towards Precision Medicine Tools for Cardiovascular Disease Risk Stratification in South Africa

Cardiovascular diseases (CVDs) are a leading cause of mortality and morbidity in South Africa. Risk stratification is the preferred approach to disease prevention, but identifying patients at high risk for CVD remains challenging. Assessing genetic risk could improve stratification and inform a clinically relevant precision medicine (PM) approach. Clinicians are critical to PM adoption, thus, this study explores practicing clinicians’ perceptions of PM-based CVD risk stratification in South Africa’s public health setting. Practicing clinicians (n = 109) at four teaching hospitals in Johannesburg, South Africa, completed an electronic self-administered survey. The effect of demographic and professional characteristics on PM-based CVD risk stratification perceptions was assessed. Fewer than 25% of respondents used clinical genetic testing, and 14% had formal genetics training. 78% had a low mean knowledge score, with higher scores associated with genetic training (p < 0.0005) and research involvement (p < 0.05). Despite limited knowledge and resources, 84% perceived PM approaches positively. 57% felt confident in applying the PM-based approach, with those already undertaking CVD risk stratification more confident (p < 0.001). High cost and limited access to genetics services are key barriers. Integrating genetic information into established clinical tools will likely increase confidence in using PM approaches. Addressing the genetics training gap and investment into the country’s genomics capacity is needed to advance PM in South Africa.

Genome-wide risk prediction of common diseases across ancestries in one million people

Polygenic risk scores (PRS) measure genetic disease susceptibility by combining risk effects across the genome. For coronary artery disease (CAD), type 2 diabetes (T2D), and breast and prostate cancer, we performed cross-ancestry evaluation of genome-wide PRSs in six biobanks in Europe, the United States, and Asia. We studied transferability of these highly polygenic, genome-wide PRSs across global ancestries, within European populations with different health-care systems, and local population substructures in a population isolate. All four PRSs had similar accuracy across European and Asian populations, with poorer transferability in the smaller group of individuals of African ancestry. The PRSs had highly similar effect sizes in different populations of European ancestry, and in early- and late-settlement regions with different recent population bottlenecks in Finland. Comparing genome-wide PRSs to PRSs containing a smaller number of variants, the highly polygenic, genome-wide PRSs generally displayed higher effect sizes and better transferability across global ancestries. Our findings indicate that in the populations investigated, the current genome-wide polygenic scores for common diseases have potential for clinical utility within different health-care settings for individuals of European ancestry, but that the utility in individuals of African ancestry is currently much lower.