A large-scale genome-wide gene-gene interaction study of lung cancer susceptibility in Europeans with a trans-ethnic validation in Asians

Introduction

Although genome-wide association studies have been conducted to investigate genetic variation of lung tumorigenesis, little is known about gene-gene (G×G) interactions that may influence the risk of non-small cell lung cancer (NSCLC).

Methods

Leveraging a total of 445,221 European-descent participants from the International Lung Cancer Consortium OncoArray project, Transdisciplinary Research in Cancer of the Lung and UK Biobank, we performed a large-scale genome-wide G×G interaction study on European NSCLC risk via a series of analyses. First, we used BiForce to evaluate and rank over 58 billion G×G interactions from 340,958 SNPs. Then, the top interactions were further tested by demographically adjusted logistic regression models. Finally, we used the selected interactions to build lung cancer screening models of NSCLC, separately, for never and ever smokers.

Results

With the Bonferroni correction, we identified eight statistically significant pairs of SNPs, which predominantly appeared in the 6p21.32 and 5p15.33 regions (e.g., rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.17, P = 6.57×10-13; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.17, P = 2.43×10-13; rs2858859HLA-DQA1 and rs9275572HLA-DQA2, ORinteraction = 1.15, P = 2.84×10-13; rs2853668TERT and rs62329694CLPTM1L, ORinteraction = 0.73, P = 2.70×10-13). Notably, even with much genetic heterogeneity across ethnicities, three pairs of SNPs in the 6p21.32 region identified from the European-ancestry population remained significant among an Asian population from the Nanjing Medical University Global Screening Array project (rs521828C6orf10 and rs204999PRRT1, ORinteraction = 1.13, P = 0.008; rs3135369BTNL2 and rs2858859HLA-DQA1, ORinteraction = 1.11, P = 5.23×10-4; rs3135369BTNL2 and rs9271300HLA-DQA1, ORinteraction = 0.89, P = 0.006). The interaction empowered polygenetic risk score (iPRS) that integrated classical polygenetic risk score and G×G information score was remarkable in lung cancer risk stratification.

Conclusions

Significant G×G interactions were identified and enriched in the 5p15.33 and 6p21.32 regions, which may enhance lung cancer screening models.