Skip to main content

Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification

Abstract

Background

Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored.

Methods

Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold.

Results

Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12–3.50, P-value = 4.13 × 10−15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99–2.49, P-value = 5.70 × 10−46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72–0.74).

Conclusions

Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

Background

Lung cancer is a multifactorial disease with high incidence and mortality [1, 2]. Environmental exposures including tobacco smoking [3,4,5,6], occupational exposures [7, 8], and air pollution [9,10,11], as well as heritable genetics, contribute to lung cancer risk [12]. Unfortunately, a majority of lung cancer cases are diagnosed at a late disease stage with a poor 5-year survival rate of less than 5%. There is, therefore, an urgent and unmet need to detect lung cancer early when prevention or earlier intervention is possible [1]. Polygenic risk scores (PRSs) have emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS) [13,14,15,16]. Lung cancer PRSs and their potential clinical utility have been explored in European and Chinese populations [2, 17,18,19,20]. Lung cancer PRSs may provide additional information for personalized risk stratification and prediction beyond non-genetic risk factors [17,18,19,20].

As PRS prediction moves towards clinical implementation in personalized medicine, accurate and unbiased PRS predictions for any single individual are needed. Unfortunately, to date, the predictive accuracy of PRS has been mostly evaluated at the population level using cohort-level metrics of prediction R2 and area under the curve (AUC), with its precision at a single individual level remaining largely unexplored [21, 22]. Moreover, whether individual-level PRS instability would influence the subsequent clinical utilization of PRS-based risk stratification and prediction is also of great interest to discover [21,22,23,24,25]. In recent studies, the individual level uncertainty in PRS estimation on various traits including height and body mass index, and diseases including breast cancer, hypertension, and dementia have been assessed in the British population using UK Biobank data [21, 22]. Another study showed that post-traumatic stress disorder and type 2 diabetes PRSs estimated among populations from different ancestries have very modest correlations at the individual level [24]. To our knowledge, there has been no study investigating the stability of lung cancer PRS for individuals and its potential impact on subsequent prediction for downstream clinical applications.

Using data from the International Lung Cancer Consortium (ILCCO) [16], we estimated lung cancer PRSs and constructed corresponding confidence/credible interval (CI) for each individual using two statistical approaches: (1) the weighted sum of 16 GWAS-derived significant single nucleotide polymorphisms (SNP) loci that have been validated in European descent population and the CI through bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes) [21]. We further evaluated the impact of individual-level PRS uncertainty on PRS-based ranking, risk stratification, and prediction in conjunction with non-genetic risk factors of age, gender, and smoking history. Our research shows that the uncertainty of lung cancer PRS at the individual level greatly impacts the subsequent performance of individual risk stratification and prediction, highlighting the importance of cautious clinical interpretation and implementation in precision medicine.

Methods

Study population

We conducted our study in 30,060 participants of European ancestry (17,166 lung cancer cases and 12,894 controls) from 25 lung cancer OncoArray studies of ILCCO including both population-based and hospital-based case-control studies. The basic characteristics of these studies are summarized in Table S1 (Fig. 1, Additional file 1: Table S1). The eligibility criteria of studies to be included in ILCCO were that they had a study protocol for subject recruitment and a structured questionnaire for baseline lifestyle information. A detailed description of each study was described in the consortium flagship paper [16, 26]. Baseline demographics of age, gender, self-reported race/ethnicity, smoking history, and lung cancer histology information were collected and adjusted for in multivariable logistic regression analyses.

Fig. 1
figure 1

Overview of the study. The study was conducted in 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium (ILCCO). The lung cancer PRSs and corresponding confidence/credible interval were constructed using two statistical approaches for each individual—(1) the weighted sum of 16 GWAS-derived significant SNP loci that have been validated in European descent population and the confidence interval through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the credible interval through posteriors sampling (PRS-Bayes). The individual-level PRS uncertainty was characterized and the impact on subsequent risk stratification and prediction were evaluated

Genotyping, imputation, and quality control

Genotyping was performed at the Center for Inherited Disease Research using the OncoArray platform with 533,631 SNPs from Illumina and imputation was conducted based on 1000G phase 3 reference panel [27]. Standard QC procedures were applied to keep variants with genotype calling rate > 95%, minor allele frequency > 1%, and deviation from Hardy-Weinberg equilibrium with P-value > 10−10 in controls. A total of 5,097,871 variants remained in our statistical analyses. To adjust for subtle population stratification, we performed a principal component analysis using PLINK v1.9 [28] and adopted the top 10 principal components (PCs) following reference [29].

PRS estimation and confidence/credible interval (CI) construction

PRS-16-CV was calculated as the weighted sum of 16 GWAS-derived common SNP loci that have been validated in European descent populations [18]. Effect sizes were estimated through five-fold cross-validation after adjusting for age, gender, smoking status, and 10 PCs. Confidence intervals (CI) for PRS-16-CV were constructed using 1000 bootstrap samples within each fold. Therefore, each individual would have their own individual-level PRS distribution characterized by the PRS-16-CV bootstrapped mean and CI for downstream risk stratification and prediction. Cross-validation and bootstrapping were implemented using R v4.1.0. In addition, we calculated the PRS-16 based on the same 16 SNP loci with effect sizes directly extracted from the GWAS catalog and previous literature (Additional file 1: Table S2-3), and no confidence interval was provided for the PRS-16 [16, 18, 30,31,32,33].

On the other hand, PRS-Bayes was estimated using a Bayesian approach leveraging the LDpred2 framework, and PRS-Bayes credible interval was constructed via posterior distribution sampling using the method proposed by Ding et al. [21, 34]. More specifically, we utilized the same training data in PRS-16-CV to estimate a posterior distribution of effect sizes, given the observed genotype and lung cancer status. For each variant, a sample of effect size estimates was drawn from the posterior distribution of the causal effect size using Markov chain Monte Carlo. A credible interval of the PRS-Bayes estimator was constructed by aggregating the number of effect allele copies weighted by each of the drawn estimates.

Characterization of PRS uncertainty at the individual level

We characterized the individual level uncertainty in PRS estimation with standard deviation (s.d.) and a pre-specified confidence/credible level of CI. With a pre-specified confidence/credible level p, we derived an individual p-level CI of PRS by obtaining the empirical (1-p)/2 and (1+p)/2 quantiles from the individual level PRS distribution.

Impact of individual-level PRS uncertainty on PRS-based risk stratification

Having obtained PRSs and corresponding CI, we stratified individuals into different PRS risk subgroups based on the relationship between their own PRS CI and the population level threshold (Fig. 2). Individuals with their PRS CI above a pre-specified population-level threshold t at the upper tail (e.g., t = 90th percentile) were classified as certainly high genetic risk and similarly for individuals with PRS CI below the population level threshold t at the lower tail (e.g., t = 10th percentile) as certainly low genetic risk. Individuals whose CI covered the population level threshold were considered uncertain. As a comparison, we classified individuals based on their estimated PRS mean and population threshold without taking individual-level certainty into account.

Fig. 2
figure 2

Risk stratification based on PRS confidence/credible interval (CI). The large distribution illustrates the PRS distribution at the population level, and the four small ones refer to individual PRS distributions for participants with different PRS-based risks of lung cancer. The dashed horizontal lines indicate the population level thresholds for risk stratification. Individuals with their PRS CI above a pre-specified population-level threshold t at the upper tail (e.g., t = 90th percentile) were classified as certainly high genetic risk, and similarly for individuals with PRS CI below the population level threshold t at the lower tail (e.g., t = 10th percentile) as certainly low genetic risk. Individuals whose CI covered the population level threshold were considered uncertain

To evaluate the degree of consistency for PRS-based risk stratification, we counted the number of overlapped individuals that are concordantly identified as the same PRS-based risk by different PRS estimators. Sensitivity analyses were conducted by changing the confidence/credible level p and risk threshold t at the population level. In addition to the two PRS approaches, we also assessed the concordance of PRS-based risk stratification among 22 lung cancer PRS in the PGS catalog.

Impact of individual level PRS uncertainty on relative risk and risk prediction

We applied multivariable logistic regression to discover the effect of PRS-based risk subgroups on lung cancer risk controlling for age, gender, and smoking history in the individuals who were identified with certainty. Stratified analyses by gender, smoking status, and lung cancer histology were similarly conducted. Lung cancer risk prediction models were constructed by integrating both PRS risk groups and other non-genetic baseline covariates, such as age, gender, and smoking history. The prediction model performance was evaluated using five-fold cross-validation based on the metric of AUC. To investigate the potential impact of individual-level uncertainty on risk prediction, we constructed a risk prediction model within the individuals identified with certainty (certainly high vs. low risk) and compared the model performance constructed in the high (n = 3006) and low (n = 3006) risk deciles identified by PRS-16-CV mean without taking individual level uncertainty into account under the exact same model specification.

Results

Table 1 summarizes the baseline characteristics of 30,060 participants that were included in the study. Ever smokers were significantly enriched in lung cancer cases with a significantly longer median smoking pack-years (39 pack-years) compared to controls (13 pack-years). Among lung cancer patients, 76.4% of the cases were non-small cell lung cancer (NSCLC) with an enrichment of lung adenocarcinoma (38.1%). Three PRSs (PRS-16, PRS-16-CV, and PRS-Bayes) of lung cancer risk were calculated for each individual and statistically higher PRS mean (mPRS) were consistently observed in lung cancer patients compared to controls (mPRS-16 in lung cancer cases: 1.26, controls: 1.21, P-value < 2.2e−16; mPRS-16-CV in lung cancer cases: 1.21, controls: 1.16, P-value < 2.2e−16; mPRS-Bayes in lung cancer cases: − 0.04, controls: − 0.11, P-value < 2.2e−16), suggesting a potentially higher genetic risk in lung cancer patients (Additional file 2: Fig. S1).

Table 1 Baseline demographic information of the study population

Characterization of lung cancer PRS uncertainty at the individual level

The variability of lung cancer PRS at the individual level was evaluated by leveraging the CIs constructed in PRS-16-CV and PRS-Bayes. A considerable variation of PRS point estimates was observed for both methods. On average, a larger standard deviation (s.d.) of individual-level PRS distribution was observed in PRS-Bayes (mean s.d.: 0.88, 95% CI of s.d.: 0.68–1.11) compared to PRS-16-CV (mean s.d.: 0.12, 95% CI of s.d.: 0.09–0.15), suggesting a larger variability using PRS-Bayes at the individual level. For illustration purposes, we show the individual-level PRS distribution for both methods in 100 individuals (Fig. 3).

Fig. 3
figure 3

Individual-level distribution of PRS-16-CV and PRS-Bayes. A Individual-level PRS distributions obtained from PRS-16-CV. B Individual-level PRS distributions obtained from PRS-Bayes. For illustration purposes, here we only show the individual-level PRS distributions for 100 participants

Impact of Individual-level PRS uncertainty on lung cancer risk stratification

All individuals were stratified into PRS deciles and at a confidence/credible level of 95%. Using PRS-16-CV, only 25.0% (751/3006) of individuals with PRS point estimates in the lowest decile of PRS and 16.8% (505/3006) in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individual (Table 2). We further conducted sensitivity analysis by changing the population level threshold t and the confidence/credible level p. We varied the range of the confidence/credible level p from 0 to 100% and fixed the threshold at t = 10th or 5th percentile for the low-risk population and at t = 90th or 95th percentile for the high-risk population. The proportion of certainty was negatively correlated with the confidence/credible level p for both the high-risk and low-risk classification (Additional file 2: Fig. S2). A consistently lower percentage of certain classifications was observed in PRS-Bayes than in PRS-16-CV across all confidence/credible levels. Similar relationships between the proportion and the confidence/credible level p were observed across subgroups of gender, histology, and smoking status (Additional file 2: Fig. S3). The proportion of certainty decreases as more stringent threshold t and confidence/credible level p are specified.

Table 2 PRS-based rankings identified by PRS-16-CV and PRS-Bayes

Furthermore, the individual-level uncertainty greatly impacts PRS-based rankings. Each individual would have a distribution of the rankings for each PRS point estimate obtained from both PRS-16-CV and PRS-Bayes. We calculated the mean and range of the rankings for each individual. Substantial variability was also observed in the rankings identified by both methods for each PRS decile (Table 2). A wider range of rankings was observed for PRS-Bayes compared to using PRS-16-CV. The minimal range of the rankings was observed in the lowest decile (mean ranking 7th, range: 0–66th) identified by PRS-16-CV, and individuals above the 90th percentile of PRS can be anywhere from the 25th–100th percentile. In contrast, individuals in each PRS decile can be anywhere from the 0th–100th using PRS-Bayes when their CI is taken into consideration.

To answer the question of whether individuals would be commonly assigned into the same genetic risk categories, we counted the number of overlapped individuals that were commonly identified using different PRS estimators of PRS-16, PRS-16-CV, and PRS-Bayes (Table 3). Overall, the degree of overlap decreases as the population level threshold increases. PRS-16 and PRS-16-CV identified a fair number of overlapped individuals as the same SNP loci were utilized. Two thousand four hundred seventy (82%) individuals were concordantly classified as high risk (> 90th percentile) using PRS-16 and PRS-16-CV, while PRS-Bayes agreed with either PRS-16-CV or PRS-16 on 19% of high-risk (> 90th percentile) stratification. In addition, we assessed the concordance of 22 lung cancer PGSs available in the PGS catalog, and similar modest correlations in PRS-based risk stratification were observed using PGS estimate (Additional file 2: Fig. S4-7). Taken together, substantial disagreement was observed for lung cancer PRS at the individual level using different PRS estimators.

Table 3 Overlap of commonly identified individuals using different PRS estimators

Impact of individual-level uncertainty on relative risk of PRS deciles on lung cancer

Next, we evaluated its impact on lung cancer risk in individuals identified by PRS-16-CV taking individual-level uncertainty into account. In contrast to risk stratification based on PRS mean, only individuals in the top and lowest deciles were able to be identified with certainty when taking variance in PRS estimates into account; hence, the relative risk was only evaluated in the two PRS-based risk subgroups. An increased effect size for PRS CI (OR = 2.73, 95% CI: 2.12–3.50, P-value = 4.13 × 10−15) was compared to using PRS mean (OR = 2.23, 95% CI: 1.99–2.49, P-value = 5.70 × 10−46) (Table 4). Similar improvement was observed in stratified analyses by gender, lung cancer histology, and smoking status (Table 5). The largest increase in relative risk of lung cancer from OR = 2.63 to OR = 4.22 was observed in never smokers, suggesting a potentially larger genetic contribution to lung cancer and thus more susceptible to individual-level PRS uncertainty.

Table 4 Impact of individual-level uncertainty on relative risk of PRS deciles on lung cancer risk
Table 5 Impact of individual-level uncertainty on the relative risk of PRS deciles on lung cancer risk in stratified analyses by gender, lung cancer histology, and smoking status

Impact of individual-level PRS uncertainty on lung cancer risk prediction

Finally, we evaluated the impact of the individual-level PRS uncertainty on lung cancer risk prediction by including the two PRS-based risk subgroups identified by PRS-16-CV. Prediction models were constructed based on PRS-based risk subgroups and non-genetic risk factors of age, gender, and smoking history. An improved prediction model performance was consistently observed in all models when considering individual-level uncertainty compared to using PRS mean solely (Table 6). The discriminative performance improved as we gradually incorporated non-genetic predictors into the model in a sequence of age, gender, and smoking information. The highest AUC (0.73, 95% CI = 0.72–0.74) was achieved by the model whose predictors were a combination of PRS-based risk categories and non-genetic predictors of age, gender, and detailed smoking pack years. Taking individual-level uncertainty of genetic risk into account in conjunction with non-genetic risk factors improves lung cancer risk prediction.

Table 6 Impact of individual-level uncertainty on lung cancer risk prediction performance

Discussion

Our work explored lung cancer PRS uncertainty for individuals in populations with European ancestry using both GWAS-derived SNPs and LDpred2 and demonstrated that individual-level PRS uncertainty greatly impacts PRS-based rankings, risk stratification, and ultimately prediction of lung cancer risk. The substantial recent interest in translating lung cancer PRS and other predictive tools like deep learning models from low radiation dose chest computed tomography (LDCT) for future lung cancer risk prediction necessitates a careful assessment of individual-level uncertainty to truly accomplish personalized risk assessment [35, 36]. Taking individual-level uncertainty into PRS-based risk stratification and prediction provides more accurate risk stratification and improves the ability to identify high-risk subjects and recommend LDCT screening with more certainty. It is imperative to develop more stable PRS and account for uncertainty at the individual level in risk stratification to help ensure reliable clinical applications of PRS in real-world settings.

A critical concern of PRS application is delivering inaccurate risk estimates at the individual level, and wrongly categorizing an individual as low or high genetic risk based on unstable PRS estimates. The downstream effect could be inappropriate or even contradictory medical advice or clinical decisions [25]. Our study shows substantial disagreement in risk categorization using different lung cancer PRSs (PRS-16, PRS-16-CV, and PRS-Bayes), suggesting that individuals that were identified as very high risk by one PRS method may not be classified as such by another. Similar modest correlations were observed for another 22 lung cancer PRS in the PGS catalog. This issue is not unique to lung cancer, as large discordances have also been found in breast cancer, hypertension, and dementia using approaches to calculate PRS in a white British population [22]. More complex traits or diseases and different ancestries between the discovery and target population of GWAS may lead to even more profound inconsistency at the individual level [24]. PRS needs to be reliable and reproducible if it is going to inform personalized decision-making in clinical settings.

Comparing the two PRS generative methods, a much larger variability was observed in PRS-Bayes, resulting in no individuals can be identified with sufficient certainty. Ding et al. also found only a limited proportion of individuals are classified as high risk with certainty across 13 traits in UK Biobank [21]. PRS-Bayes includes a much greater number of non-zero weight SNPs (~2000 SNPs) with small effect sizes and this may explain a higher proportion of SNP-based heritability compared to PRS-16-CV (16 SNPs). In the meanwhile, the numerous SNPs with small effect sizes may contribute to a higher variance in PRS estimates. PRS-Bayes may perform better in terms of lung cancer risk prediction at the population level as the accuracy is determined by the proportion of phenotypic variance explained by variants included and the improved population prediction error. On the other hand, PRS-16-CV is relatively parsimonious and contains only those that have been robustly validated SNPs and have a minor individual-level variance. Additionally, we constructed and evaluated the individual-level PRS uncertainty and population-level prediction accuracy using nine experimentally validated SNPs that are associated with lung cancer risk [37] (Additional file 1: Table S4). More individuals can be identified with certainty while similar prediction accuracy was achieved (Additional file 1: Table S5-6). This may suggest that experimentally validated fine-mapped variants may be more likely to be biologically causal variants and thus have stable effects and PRS risk stratification performance at the individual level for lung cancer prediction given the similar population-level prediction accuracy. A trade-off between these methods may be manifested by the compromise between prediction performance at the population level and stability at the individual level, and this needs to be considered cautiously when developing new PRS methods tailoring for different purposes and applications [21].

Given the individual-level PRS uncertainty, it is imperative to cautiously interpret and implement it in clinical applications. From our analyses, using the PRS-based risk subgroup alone resulted in the lowest AUC, and including other non-genetic factors of age, gender, and smoking history largely improved the lung cancer risk prediction performance. PRS stand-alone is imperfect and the use of PRS as a covariate or effect modifier in conjunction with non-genetic risk factors may inform better outcomes [22, 38]. Despite modest improvement in prediction accuracy using PRS-based risk subgroups taking individual-level uncertainty into account, likely due to the limited sample sizes of eligible individuals, it attains both PRS stability and predictive accuracy given the data. Also, it reflects the situation that it is challenging for the current PRS approaches to balance reproducible PRS for individuals with sufficient certainty and high prediction accuracy at the population level. It is crucial and pressing to develop guidelines in constructing PRS to minimize the extent to which individuals could be provided with imprecise or contradictory clinical advice and intervention as well as PRS reporting standards in light of the patient and primary care providers’ perspective [25].

There are several improvements and future directions for the present study. First, more statistical approaches to construct lung cancer PRS, e.g., regularization-based methods, can be included [39]. Second, we did not assess the potential interactions between lung cancer PRS and non-genetic risk factors when taking individual-level uncertainty into account, and we would expect more uncertainty arising from environmental factors as well as from the gene-environment interactions in predicting the risk of lung cancer [40]. Further research is needed to assess the impact on performance when using commonly applied LD clumping methods for SNP selection, particularly with varying parameters. Lastly, our study was conducted in a population of European ancestry; thus, the generalizability of the results in other populations of different genetic ancestries is a concern given the discrepancy of genetic structures for causal variants, allele frequencies, and LD patterns across ancestries. To realize the equitable and reliable potential of PRS, it is necessary to carry out larger, multi-ancestry GWAS and further investigate the uncertainty at the individual level in PRS.

Conclusions

Our study characterized the uncertainty of lung cancer PRS at the individual level and evaluated its impact on subsequent PRS-based ranking, risk stratification, and prediction in populations with European ancestry. It is imperative to develop reproducible PRS, reliable PRS-based clinical applications as well as guidelines to report and communicate PRS to patients and their primary physicians.

Availability of data and materials

Genotype and phenotype data used in this work are available from dbGaP (https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs001273.v4.p2; phs001273.v3.p2; )[2].

References

  1. Center., A.C.S.C.S Key statistics for lung cancer. 2022 [cited 2022 March 12].

    Google Scholar 

  2. Byun J, et al. Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer. Nat Genet. 2022;54(8):1167–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Cancer, I.A.f.R.o. IARC monographs on the evaluation of carcinogenic risks to humans. Vol. 83. Lyon: IARC; 2004. p. 1452.

    Google Scholar 

  4. Doll R, et al. Mortality in relation to smoking: 50 years' observations on male British doctors. BMJ. 2004;328(7455):1519.

    Article  PubMed  PubMed Central  Google Scholar 

  5. O'Keeffe LM, et al. Smoking as a risk factor for lung cancer in women and men: a systematic review and meta-analysis. BMJ Open. 2018;8(10):e021611.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Wang X, et al. Association between smoking history and tumor mutation burden in advanced non-small cell lung cancer. Cancer Res. 2021;81(9):2566–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Nawrot TS, et al. Association of total cancer and lung cancer with environmental exposure to cadmium: the meta-analytical evidence. Cancer Causes Control. 2015;26(9):1281–8.

    Article  PubMed  Google Scholar 

  8. van der Bij S, et al. Lung cancer risk at low cumulative asbestos exposure: meta-regression of the exposure-response relationship. Cancer Causes Control. 2013;24(1):1–12.

    Article  PubMed  Google Scholar 

  9. Raaschou-Nielsen O, et al. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European Study of Cohorts for Air Pollution Effects (ESCAPE). Lancet Oncol. 2013;14(9):813–22.

    Article  PubMed  Google Scholar 

  10. Lissowska J, et al. Lung cancer and indoor pollution from heating and cooking with solid fuels: the IARC international multicentre case-control study in Eastern/Central Europe and the United Kingdom. Am J Epidemiol. 2005;162(4):326–33.

    Article  PubMed  Google Scholar 

  11. Blechter B, et al. Sub-multiplicative interaction between polygenic risk score and household coal use in relation to lung adenocarcinoma among never-smoking women in Asia. Environ Int. 2021;147:105975.

    Article  PubMed  Google Scholar 

  12. Malhotra J, et al. Risk factors for lung cancer worldwide. Eur Respir J. 2016;48(3):889–902.

    Article  PubMed  Google Scholar 

  13. Bosse Y, Amos CI. A decade of GWAS results in lung cancer. Cancer Epidemiol Biomarkers Prev. 2018;27(4):363–79.

    Article  PubMed  Google Scholar 

  14. Wang Y, et al. Deciphering associations for lung cancer risk through imputation and analysis of 12,316 cases and 16,831 controls. Eur J Hum Genet. 2015;23(12):1723–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Timofeeva MN, et al. Influence of common genetic variation on lung cancer risk: meta-analysis of 14 900 cases and 29 485 controls. Hum Mol Genet. 2012;21(22):4980–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. McKay JD, et al. Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat Genet. 2017;49(7):1126–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Hung RJ, et al. Assessing lung cancer absolute risk trajectory based on a polygenic risk model. Cancer Res. 2021;81(6):1607–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Jia G, et al. Evaluating the utility of polygenic risk scores in identifying high-risk individuals for eight common cancers. JNCI Cancer Spectr. 2020;4(3):pkaa021.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Dai J, et al. Identification of risk loci and a polygenic risk score for lung cancer: a large-scale prospective cohort study in Chinese populations. Lancet Respir Med. 2019;7(10):881–91.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Jia G, et al. Incorporating both genetic and tobacco smoking data to identify high-risk smokers for lung cancer screening. Carcinogenesis. 2021;42(6):874–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Ding Y, et al. Large uncertainty in individual polygenic risk score estimation impacts PRS-based risk stratification. Nat Genet. 2022;54(1):30–9.

    Article  CAS  PubMed  Google Scholar 

  22. Clifton L, et al. Assessing agreement between different polygenic risk scores in the UK Biobank. Sci Rep. 2022;12(1):12812.

    Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

  23. Torkamani A, Wineinger NE, Topol EJ. The personal and clinical utility of polygenic risk scores. Nat Rev Genet. 2018;19(9):581–90.

    Article  CAS  PubMed  Google Scholar 

  24. Schultz LM, et al. Stability of polygenic scores across discovery genome-wide association studies. HGG Adv. 2022;3(2):100091.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Lewis ACF, et al. Patient and provider perspectives on polygenic risk scores: implications for clinical reporting and utilization. Genome Med. 2022;14(1):114.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Amos CI, et al. The OncoArray Consortium: a network for understanding the genetic architecture of common cancers. Cancer Epidemiol Biomarkers Prev. 2017;26(1):126–35.

    Article  PubMed  Google Scholar 

  27. Genomes Project Consortium, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74.

    Article  ADS  Google Scholar 

  28. Chang CC, et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Marees AT, et al. A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int J Methods Psychiatr Res. 2018;27(2):e1608.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Wang Y, et al. Common 5p15.33 and 6p21.33 variants influence lung cancer risk. Nat Genet. 2008;40(12):1407–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Qin N, et al. Comprehensive functional annotation of susceptibility variants identifies genetic heterogeneity between lung adenocarcinoma and squamous cell carcinoma. Front Med. 2021;15(2):275–91.

    Article  PubMed  Google Scholar 

  32. Sakaue S, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53(10):1415–24.

    Article  CAS  PubMed  Google Scholar 

  33. Wang Y, et al. Rare variants of large effect in BRCA2 and CHEK2 affect risk of lung cancer. Nat Genet. 2014;46(7):736–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Prive F, Arbel J, Vilhjalmsson BJ. LDpred2: better, faster, stronger. Bioinformatics. 2020;36:5424–543.

    Article  CAS  PubMed Central  Google Scholar 

  35. Mikhael PG, et al. Sybil: a validated deep learning model to predict future lung cancer risk from a single low-dose chest computed tomography. J Clin Oncol. 2023;41(12):2191–200.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Silvestri GA, Jett JR. The intersection of lung cancer screening, radiomics, and artificial intelligence: can one scan really predict the future development of lung cancer? J Clin Oncol. 2023;41(12):2141–3.

    Article  PubMed  Google Scholar 

  37. Long E, et al. Functional studies of lung cancer GWAS beyond association. Hum Mol Genet. 2022;31(R1):R22–36.

    Article  MathSciNet  PubMed  PubMed Central  Google Scholar 

  38. Huntley C, et al. Utility of polygenic risk scores in UK cancer screening: a modelling analysis. Lancet Oncol. 2023;24(6):658–68.

    Article  PubMed  Google Scholar 

  39. Speed D, Balding DJ. MultiBLUP: improved SNP-based prediction for complex traits. Genome Res. 2014;24(9):1550–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wang H, et al. Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank. Sci Adv. 2019;5(8):eaaw3538.

    Article  PubMed  PubMed Central  ADS  Google Scholar 

Download references

Acknowledgements

We thank all the participants from the International Lung Cancer Consortium for sharing their data.

Funding

XW, ZZ, and DCC are supported by 5U01CA209414 from the National Cancer Institute. XL is supported by R35-CA197449 and U19-CA203654 from the National Cancer Institute, and U01-HG009088 and U01HG012064 from NHGRI.

Author information

Authors and Affiliations

Authors

Contributions

XW, ZZ, and DCC conceived the study; XW, ZZ, YD, and TC carried out the analyses. The International Lung Cancer Consortium coordinated genotyping and sequencing data acquisition. XL and DCC supervised the whole work and the statistical analyses. XW and ZZ drafted the manuscript and generated the figures. XW, ZZ, YD, and TC made critical revisions. All authors read and approved the final manuscript.

Corresponding author

Correspondence to David C. Christiani.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the institutional review board of the Massachusetts General Brigham Health Care System (MGB) protocol number: IRB Protocol No: 1999P004935. All participants provided written informed consent to participate in the study. The study was performed in accordance with the principles of the Helsinki Declaration.

Consent for publication

Not applicable.

Competing interests

M. C. Aldrich reports funding from NIH and GO2 foundation, and she serves an advisor to Guardant Health. While M. Johansson and J. D. McKay are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy, or views of the International Agency for Research on Cancer/World Health Organization. The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

Participating samples and their origin sites in the ILCCO lung cancer OncoArray project. Table S2 Information of 16 GWAS-derived lung cancer SNPs included in PRS-16 and PRS-16-CV. Table S3 Information of 19 GWAS-derived lung cancer SNPs that have been validated in Caucasians from prior studies. Table S4 Fine-mapped lung cancer risk variants that have been experimentally validated. Table S5 Rankings of individuals identified by three PRS models. Table S6 AUC of predicting lung cancer in the confident individuals.

Additional file 2: Figure S1

The population-level distribution of lung cancer risk PRSs. Figure S2 PRS CI-based stratification uncertainty. Figure S3. PRS-based stratification uncertainty in subgroups. Figure S4-S7 Proportions of concordant individuals by 22 lung cancer PRSs in the PGS catalog based on different population thresholds.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, X., Zhang, Z., Ding, Y. et al. Impact of individual level uncertainty of lung cancer polygenic risk score (PRS) on risk stratification. Genome Med 16, 22 (2024). https://doi.org/10.1186/s13073-024-01298-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13073-024-01298-4

Keywords