Multi-ancestry meta-analyses of lung cancer in the Million Veteran Program reveal novel risk loci and elucidate smoking-independent genetic risk

Bryan R. Gorman, Sun-Gou Ji, Michael Francis, Anoop K. Sendamarai, Yunling Shi, Poornima Devineni, Uma Saxena, Elizabeth Partan, Andrea K. DeVito, Jinyoung Byun, Younghun Han, Xiangjun Xiao, Don D. Sin, Wim Timens, Jennifer Moser, Sumitra Muralidhar, Rachel Ramoni, Rayjean J. Hung, James D. McKay, Yohan Bossé, Ryan Sun, Christopher I. Amos, VA Million Veteran Program, Saiju Pyarajan
doi: https://doi.org/10.1101/2024.04.25.24306313

Abstract

Lung cancer remains the leading cause of cancer mortality, despite declines in smoking rates. Previous lung cancer genome-wide association studies (GWAS) have identified numerous loci, but separating the genetic risks of lung cancer and smoking behavioral susceptibility remains challenging. We performed multi-ancestry GWAS meta-analyses of lung cancer using the Million Veteran Program (MVP) cohort and a previous study of European-ancestry individuals, comprising 42,102 cases and 181,270 controls, followed by replication in an independent cohort of 19,404 cases and 17,378 controls. We further performed conditional meta-analyses on cigarettes per day and identified two novel, replicated loci, including the 19p13.11 pleiotropic cancer locus in LUSC. Overall, we report twelve novel risk loci for overall lung cancer, lung adenocarcinoma (LUAD), and squamous cell lung carcinoma (LUSC), nine of which were externally replicated. Finally, we performed phenome-wide association studies (PheWAS) on polygenic risk scores (PRS) for lung cancer, with and without conditioning on smoking. The unconditioned lung cancer PRS was associated with smoking status in controls, illustrating reduced predictive utility in non-smokers. Additionally, our PRS demonstrates smoking-independent pleiotropy of lung cancer risk across neoplasms and metabolic traits.

Introduction

Lung cancer remains the leading cause of overall cancer mortality, as the most prevalent cancer type in men, and the second highest in women after breast cancer1–3. Despite declines in smoking rates in the US since the 1980s4, tobacco use is currently implicated in upwards of 80% of lung cancer diagnoses1. Even in those who have never smoked, nor had meaningful exposure to environmental carcinogens1,5, there exists a heritable risk component of lung cancer conferred by genetic factors6–8. Differentiating the mutations which directly predispose an individual to lung cancer from those whose effect is mediated through environmental components remains challenging.

Genome-wide association studies (GWAS) have identified lung cancer risk variants associated with oncogenic processes such as immune response7, cell cycle regulation9, and those affecting DNA damage response and genomic stability8. Several lung cancer GWAS have also reported strong effects of genes such as CHRNA nicotine receptor genes which putatively increase the risk of lung cancer through behavioral predisposition towards smoking5. Characteristic molecular markers and genetic risk factors in smokers and never-smokers have been identified10,11, though fewer variants have been found in GWAS performed exclusively in never-smokers12.

Lung cancer has a heterogeneous genetic architecture across ancestral groups13,14. In the two most well-studied ancestries, European (EA) and East Asian (EAS), the majority of genome-wide significant loci are not shared15,16; this is in agreement with molecular studies showing differences in tumor characteristics between EA and EAS17. Smaller African ancestry (AA) cohorts have replicated known loci from EA or EAS8,18, though no AA-specific GWAS loci have been reported.

In this study, we examined lung cancer genetic variation in EA as well as in the largest AA cohort to-date. Our discovery analysis is performed in an older cohort of mostly male US veterans in the Department of Veterans Affairs Million Veteran Program (MVP)19. Lung cancer incidence is approximately twice as high in men than in women2, and additionally MVP contains a large number of cigarette smokers, positioning this biobank as particularly valuable for this analysis. We performed GWAS in overall cases of lung cancer as well as two non-small cell lung cancer (NSCLC) subtypes, adenocarcinoma (LUAD) and squamous cell lung carcinoma (LUSC).

Results

Genome-wide association studies for lung cancer

We performed a GWAS on overall lung cancer within EA participants in MVP (10,398 lung cancer cases and 62,708 controls; Supplementary Data 1), followed by a meta-analysis with the EA International Lung Cancer Consortium OncoArray study (ILCCO; McKay et al., 2017)7, for a total of 39,781 cases and 119,158 controls (Supplementary Fig. 1). The EA meta-analysis for overall lung cancer identified 26 conditionally independent SNPs within 17 genome-wide significant loci (P<5×10−8; Supplementary Fig. 2a; Supplementary Data 2). All 12 loci reported by ILCCO7 were confirmed, with consistent direction of effect on all single nucleotide polymorphisms (SNPs) with P<1×10−5, as well as high correlation of effect sizes and allele frequency (Supplementary Fig. 3). Of the 17 genome-wide significant loci for overall lung cancer, four were novel with respect to the broader literature: neuronal growth regulator LSAMP, WNT signaling regulator NMUR2, DNA damage repair protein XCL2, and hedgehog signaling regulator TULP3, (Table 1; Supplementary Fig. 4a-d).

View this table:
Table 1:

Novel genome-wide significant loci and their respective index variants associated with lung cancer risk in European-ancestry meta-analyses from MVP and ILCCO7 cohorts, MVP African ancestry, multi-ancestry meta-analyses, and in European-ancestry meta-analyses after conditioning on cigarettes per day. LUAD, adenocarcinoma; LUSC, squamous cell carcinoma; EA, effect allele; NEA, non-effect allele; EAF, effect allele frequency in the given population; OR (95% CI), odds ratio and 95% confidence interval.

Further association tests stratified by cancer subtypes LUAD and LUSC in MVP EA (Supplementary Fig. 2bc; Supplementary Data 3-4) replicated associations reported by ILCCO7 (Supplementary Fig. 3) and identified additional loci. Two novel EA meta-analysis loci were identified for LUAD, proto-oncogene MYC and Wnt signaling inhibitor TLE3 (Table 1; Supplementary Fig. 4e-h). For LUSC, we identified one novel locus at 10q24.31 near NFκB inhibitor CHUK and BLOC1S2. Across all subtypes for EA meta-analysis index variants, the MVP cohort had associations with P<0.05 in all but one in overall lung cancer, five in LUAD, including approximately nominal significance at rs67824503 (MYC; P=0.057), and one in LUSC (Supplementary Data 2-4).

We investigated expression quantitative trait loci (eQTL) relationships between top SNPs from the EA meta-analysis across all lung cancer GWAS in GTEx v8 Lung20 and the Lung eQTL Consortium21 (Supplementary Data 2-4). This analysis showed that the LUSC index SNP rs36229791 on 10q24.31 was associated with the mRNA expression levels of BLOC1S2 (Fig. 1a-d), consistent with previous TWAS22. BLOC1S2 is an oncogene whose gene product is associated with centrosome function; centrosomal abnormalities have previously been observed in vitro in LUSC23,24.

Figure 1.
Figure 1. Highlighted novel GWAS loci.

a-d) The meta-analysis of squamous cell lung carcinoma (LUSC) in European ancestry (EA) identifies a novel locus at 10q24.31. a) Odds ratios for rs36229791 in LUSC compared to lung adenocarcinoma (LUAD) and overall lung cancer. b) BLOC1S2 expression varies by genotype at rs36229791. c) BLOC1S2 eQTL t statistic vs LUSC z statistic. d) Regional association plot showing SNP significance and genes around lead SNP rs36229791. e) The African ancestry GWAS highlights a putatively novel locus on chr12 at LINC00944. The risk allele has effectively 0% frequency in EA.

We improved our variant selection by fine-mapping and estimating credible sets of candidate causal variants in EA meta analysis using sum of single effects (SuSiE)25,26 modeling. For overall lung cancer, LUAD, and LUSC, we identified 23, 23, and 9 high quality credible sets, respectively, containing 370, 246, and 192 total SNPs (Supplementary Data 5).

GWAS in AA

We analyzed overall lung cancer risk in 2,438 cases and 62,112 controls of African ancestry (AA), the largest AA GWAS discovery cohort to date (Supplementary Fig. 5a). Two loci reached genome-wide significance in our discovery scan: 15q25, replicating the association in CHRNA5 for AA populations reported by an earlier GWAS18, and a putative novel locus at 12q23 with index SNP rs78994068 (Table 1; Fig. 1e). We further performed GWAS in AA within LUAD and LUSC subtypes but found no genome-wide significant associations (Supplementary Fig. 5b-c).

The putative AA locus at 12q23 is driven by six SNPs in high linkage disequilibrium (LD; R2>0.8) found in long non-coding RNAs LINC00943 and LINC00944 (Fig. 1e). These imputed SNPs all had odds ratios (ORs) close to 2, with 1.3% frequency in AA and 0% in EA, consistent with gnomAD v3. LINC00944 is highly expressed in immune cells and blood, and enriched in T cell pathways in lung tissue and cancer27–30. We fine-mapped this locus to define a 95% credible set (Supplementary Data 6), and annotated the functional consequence of the variants using the Variant Effect Predictor (VEP)31. Two variants, rs78994068 and rs115962601, were in a known enhancer regulatory region (ENSR00000974920) and thus may involve regulatory changes. However, this locus was directionally consistent but not significant in our AA replication cohort (discussed below); therefore, larger-scale AA analyses are needed to confirm this finding.

GWAS multi-ancestry meta-analysis

We conducted a fixed-effect inverse variance-weighted multi-ancestry meta-analysis, combining the EA meta-analysis and the MVP AA GWAS for overall lung cancer, LUAD, and LUSC (Supplementary Data 7-9; Supplementary Fig. 6a-c). This analysis identified two additional novel genome-wide significant loci in overall lung cancer (Table 1; Supplementary Fig. 4i-j): ubiquitin ligase JADE2, previously associated with smoking initiation32, and RNA polymerase-associated RPAP3. Neither of these novel multi-ancestry meta-analysis loci were reported in a recent multi-ancestry analysis by Byun et al.8 that included fewer AA and more Asian ancestry samples, indicating the value our larger AA sample provided for novel discovery. All genome-wide significant EA meta-analysis associations reached genome-wide significance in the multi-ancestry meta-analysis except rs11855650 (TLE3) in LUAD (P=6.19×10−8). We additionally performed random effects meta-analyses using the Han-Eskin method (RE2)33, and observed similar P-values to the fixed effect meta-analysis, with all index variants PRE2<5×10−8 (Supplementary Data 7-9).

Polygenic risk scoring

To gain an understanding of the penetrance and pleiotropy of lung cancer risk, we constructed PRSs based on the ILCCO summary statistics7 for every EA subject in MVP. As expected, the PRS was highly associated with both lung cancer risk as well as smoking behavior (Supplementary Fig. 7a-b). Even after removing individuals with any history of lung cancer risk to prevent enrichment of risk factors and comorbidities, the association with smoking behavior remained, suggesting that the PRS is partially capturing genetic smoking behavioral risk factors (Supplementary Fig. 7c). In all groups, individuals at the top decile of the PRS were at significantly higher risk of lung cancer than those in the lowest decile.

Multi-trait conditional analysis for smoking status

Despite adjusting for smoking status, both in MVP EA and ILCCO7, a significant genetic correlation was observed between all subsets of lung cancer GWAS and a recently published GWAS of smoking behaviors34 (Fig. 2a, Supplementary Data 10). In order to remove all residual effects of smoking on lung cancer susceptibility, we conducted a multi-trait-based conditional and joint analysis (mtCOJO)35,36, conditioning on a GWAS for cigarettes per day34, which was the smoking trait most strongly correlated with overall lung cancer and subtype GWAS from the EA meta-analysis. Because lung cancer case selection also preferentially selects smokers, conventional adjustment for smoking may inadvertently cause selection bias, which functions as a collider to induce biased genetic effects37. mtCOJO is considered more robust to potential collider bias than conventional covariate adjustment35,36. The total observed-scale SNP-heritability38 of lung cancer risk decreased substantially after conditioning on cigarettes per day, from 5.4% to 3.1% in overall LC, from 6.7% to 5.5% in LUAD, and from 5.8% to 3.8% in LUSC (Fig. 2b; Supplementary Data 11).

Figure 2.
Figure 2. Association of lung cancer GWAS with smoking behaviors.

a) Genetic correlations (with 95% confidence interval) between the lung cancer GWAS and smoking behaviors, including smoking initiation, cigarettes per day, smoking cessation, and age of initiation. b) SNP heritability for the meta-analysis and conditional meta-analysis. The heritability decreases in the conditional analysis for overall lung cancer as well as both subtypes, suggesting that some portion of the heritability of lung cancer is due to smoking behavior. c) Polygenic risk scores (PRS) based on standard lung cancer GWAS (blue) performs worse in never-smokers than former or current smokers, while conditioning on smoking behavior (orange) results in similar performance.

Significant loci from the conditional analyses are shown in Supplementary Fig. 8-9 and Supplementary Data 12-14. As expected, the statistical significance of loci harboring smoking-related genes (e.g., CHRNA5, CYP2A6, CHRNA4) dropped to below genome-wide significance after conditioning (Fig. 3). Conversely, five signals (four loci) became significant only after conditioning, including novel signals at MMS22L in overall lung cancer and 19p13.1 (ABHD8) in LUSC. MMS22L is a novel GWAS signal but was previously identified as overexpressed in lung cancer in genome-wide gene expression scan39. These may represent biological lung cancer signals partially masked by countervailing genetic effects on smoking behavior. We performed fine-mapping to identify candidate causal variants in the conditioned EA meta-analysis summary statistics, and for overall lung cancer, LUAD, and LUSC, we identified 11, 15, and 6 high quality credible sets, respectively, containing a total of 243, 277, and 78 SNPs (Supplementary Data 5).

Figure 3.
Figure 3. Forest plot of genome-wide significant associations.

Within each cancer subtype, changes in effect size and significance are shown before and after conditioning on cigarettes per day. Novel loci are indicated by an asterisk after the gene name (*). Loci that became significant after conditioning (P<5×10−8) are in red.

We constructed PRS based on mtCOJO-conditioned ILCCO summary statistics7 to directly compare the predictive performance of PRS derived from the conditioned and non-conditioned GWAS in MVP EA. While the PRS based on the non-conditioned overall lung cancer GWAS exhibited reduced performance in never-smokers compared to ever-smokers, the PRS based on the conditional analysis resulted in similar performance across smoking status (Fig. 2c; Supplementary Data 15).

Replication of novel variants in OncoArray and combined meta-analysis

We queried the OncoArray Consortium Lung Study (OncoArray) as an external non-overlapping replication dataset for our significant GWAS signals (Supplementary Data 16-17). For GWAS in EA meta-analysis for overall lung cancer, LUAD, and LUSC, we replicated five of seven novel loci (P<0.01) in an OncoArray European ancestry cohort: XCL2 and TLE3 in overall lung cancer, MYC and TLE3 in LUAD, and BLOC1S2 in LUSC. The novel African ancestry association for overall lung cancer at LINC00944 was not replicated. We meta-analyzed OncoArray European and African ancestry participants to replicate our multi-ancestry meta-analysis signals for overall lung cancer at RPAP3 (P=0.0044) and JADE2 which bordered on nominal significance (rs329122; P=0.053). For the two novel loci which were identified in EA meta-analysis conditioned on cigarettes per day, we included smoking as a covariate for association analysis in the OncoArray European ancestry cohort. These association signals were replicated for overall lung cancer at MMS22L (P=0.006) and LUSC at ABHD8 (P=0.003). In a variant-level replication of 137 conditionally independent discovery associations which fell within ≤1 Mb of a previously reported lung cancer GWAS signal, 134 had P<0.05 in OncoArray, and 42 had P<5×10−8 (Supplementary Data 18).

We then performed a combined meta-analysis of our discovery results with OncoArray replication results (Supplementary Data 18). We considered a conservative threshold of P=4.17×10−9 (P=5×10−8/12 total GWAS analyses) to be significant, which was met by 9 of the 12 loci. Because rs329122 in JADE2 achieved the more conservative significance threshold (P=3.69×10−9), and has also been associated with smoking behavior32 and identified as a splicing-related variant associated with lung cancer40, we considered this locus to be replicated. In the combined meta-analysis we observed similar P-values in fixed effects and random effects (RE2) models.

Next, for all previously reported lung cancer and subtype loci in this study, we identified lung cancer associations from GWAS Catalog which fell within the same loci as our index variants (Supplementary Data 19). We confirmed two loci that previously had been reported only in a recent genome-wide association by proxy (GWAx) of lung cancer41: CENPC (rs75675343) in overall lung cancer in the EA meta-analysis (P=2.40×10−8) and the multi ancestry meta-analysis, and TP53BP1 in overall lung cancer in the multi-ancestry meta-analysis (rs9920763; P=1.63×10−8). Our multi-ancestry meta-analysis for overall lung cancer also confirmed a recently reported locus at 4q32.2 (NAF1)15 in East Asian ancestry.

Multi-trait analysis with breast cancer

At 19p13.1, a known pleiotropic cancer locus42,43, the index SNP of LUSC conditioned on smoking (rs61494113) sits in a gene-rich region where a recent fine-mapping effort of breast cancer risk loci44 proposed two independent associations, one affecting the regulation of ABHD8 and MRPL34, and another causing a coding mutation in ANKLE1. Here, we used the increased power provided by a multi-trait analysis of GWAS (MTAG)45 of LUSC and estrogen receptor negative (ER−) breast cancer46 to disentangle the complex relationships between cancer risk and the genes in this locus (Fig. 4a). Overexpression of ABHD8 has been shown to significantly reduce cell migration42,43. Similar odds ratios at rs61494113 were observed across LUSC and breast cancer, and MTAG enhanced the GWAS signal at this locus (Fig. 4b).

Figure 4.
Figure 4. Significant locus after conditioning on smoking behavior, 19p13.11, has pleiotropic associations with ER-negative breast cancer.

a) Regional association plot of the 19p13.11 multi-trait analysis of GWAS (MTAG) locus. b) Odds ratios for lead SNP rs61494113 in squamous cell lung carcinoma (LUSC), before and after conditioning, and MTAG analysis, compared to lung adenocarcinoma and overall lung cancer. c) ABHD8 expression varies by genotype at rs61494113. d) ABHD8 eQTL t statistic vs LUSC z statistic; red X’s indicate the 95% credible set.

We used the coloc-SuSiE method47 to assess colocalized associations between pairs of credible sets in this locus underlying the risk of LUSC and ER− breast cancer, allowing for multiple causal signals. We found evidence for a shared causal signal between credible sets in the LUSC conditional meta-analysis and ER− breast cancer (97.7% posterior probability; Supplementary Data 20). The index SNPs for the credible sets of LUSC conditioned on smoking and ER− breast cancer (rs61494113 and rs56069439, respectively) have r2=0.99.

The eQTL effect of ABHD8 was replicated in multiple tissues of GTEx v8, including Lung (Fig. 4c). Interestingly, the group of SNPs in the LUSC-BC credible set did not have the most significant eQTL effect, suggesting a complex relationship between the multiple causal variants at the locus and gene expression (Fig. 4d). For instance, a recent splice variant analysis48 implicated splicing of BABAM1 (a BRCA1-interacting protein) as a culprit of the associations observed in 19p13.1. Consistent with previous reports42,43, the cancer risk-increasing haplotype was correlated with increased expression of ABHD8 and alternative splicing of BABAM1. However, there was no overlap between the 95% eQTL credible sets of ABHD8 and BABAM1, and neither of the credible sets included rs61494113.

Phenome-wide association study

Finally, to investigate the pleiotropy of lung cancer genetic risk in the absence of the overwhelming effect of smoking behavior, we performed PheWAS in MVP using the PRS scores constructed from the ILCCO summary statistics7 for overall lung cancer, both based on the standard GWAS (“unconditioned PRS”; Fig. 5a; Supplementary Data 21) and the GWAS conditioned on cigarettes per day using mtCOJO (“conditioned PRS”; Fig. 5b; Supplementary Data 22). Each PRS was tested for association with 1,772 phecode-based phenotypes. Overall, 240 phenotypes were associated with the unconditioned PRS and 112 were associated with the conditioned PRS at a Bonferroni-corrected significance threshold (P<0.05/1,772). Although lung cancer remained a top association with the conditioned PRS, the association with tobacco use disorder was greatly reduced, from an OR associated with a standard deviation increase in the PRS of 1.151 [1.142-1.160] (P=2.32×10−237) in the unconditioned PRS to OR=1.046 [1.038-1.053] (P=1.05×10−32) in the conditioned PRS. However, the effect on alcohol use disorder was only modestly attenuated between the unconditioned (OR=1.098 [1.089-1.108]; P=1.05×10−87) and conditioned LC (OR=1.078 [1.069-1.088], P=4.41×10−60) PRSs. Whether a role for alcohol in lung cancer exists independently of smoking is controversial49,50; this analysis suggests that may be the case. Other putatively smoking-related associations, such as chronic obstructive pulmonary disease, pneumonia, and peripheral vascular disease were greatly diminished with the conditioned PRS. Mood disorders, depression, and post-traumatic stress disorder, were also significantly associated with the unconditioned PRS but no longer significantly associated with the conditioned PRS, reflecting neuropsychiatric correlates of smoking behavior.

Figure 5.
Figure 5. Phenome-wide association study (PheWAS) of polygenic risk scores (PRS) of lung cancer and lung cancer conditioned on cigarettes per day.

a) PheWAS of PRS on lung cancer is mostly confounded with smoking associations. b) PheWAS of the conditional meta-analysis PRS shows associations with skin cancer and metabolic traits.

Intriguingly, a category of metabolic traits that were not associated with the unconditioned PRS were highly associated with the conditioned PRS and in a negative effect direction. We observed protective associations of the conditioned PRS with metabolic traits such as type 2 diabetes (OR=0.945 [0.938-0.952], P=9.46×10−52) and obesity (OR=0.952 [0.945-0.959], P=2.48×10−41). Neither were associated with the unconditioned PRS (OR=1.006 [0.999-1.014]; P=0.092, and OR=1.005 [0.998-1.012]; P=0.183, respectively). Other traits in this category included sleep apnea and hyperlipidemia. These findings are consistent with prior observational findings of an inverse relationship between BMI and lung cancer51 and illustrate the extent to which smoking may be a major confounder of this relationship.

Finally, we observed strong associations of the lung cancer PRS with skin cancer and related traits, such as actinic keratitis. In basal cell carcinoma, the OR increased from 1.087 [1.072-1.102] (P=6.06×10−32) with the unconditioned PRS to 1.105 [1.090-1.120] (P=1.82×10−47) with the conditioned PRS. As a sensitivity analysis, we tested the strength of this association after removing the TERT locus, which is prominently associated with both traits. Doing so only modestly reduced the effect of the conditioned PRS to OR=1.092 [1.077-1.107] (P=4.08×10−36). Thus, our results are consistent with a genome-wide genetic correlation between lung cancer and basal cell carcinoma that is strengthened when the effect of smoking is removed. Overall, our results suggest that the biology underlying lung cancer risk may be partially masked by the residual genetic load of smoking.

Discussion

We identified novel lung cancer-associated loci in a new cohort of EA and AA participants, including the largest AA cohort analyzed to-date. We also show that, despite studies on the genetic basis of lung cancer risk taking smoking status into account, the effects of smoking continue to obfuscate our understanding of lung cancer genetics. In particular, we report two novel loci, at MMS22L (overall) and ABHD8 (LUSC), which may be partially masked by countervailing genetic effects on smoking. Our replication analysis which adjusted for smoking pack-years confirmed these loci. Additionally, our analyses demonstrated that PRSs for lung cancer contain large uncorrected genetic loading for smoking behavioral factors. Our results indicate that controlling for these factors can improve risk assessment models, potentially improving lung cancer screening even for non-smokers. Finally, our phenomic scans comparing PRSs derived from GWAS with and without genomic conditioning on smoking showed divergent associations across numerous traits, especially metabolic phenotypes.

The increased sample size in this study enabled the interpretation of multiple causal variants underlying the gene-rich ADHL8-BABAM1 region, synthesizing prior observations into a clearer understanding of this locus. Our other novel loci strengthen established lung cancer mechanisms. We identify for the first time a susceptibility locus at MYC, a well-known oncogene and master immune regulator. XCL2 is involved in cellular response to inflammatory cytokines52. LSAMP is a tumor suppressor gene in osteosarcoma53, and 3q13.31 homozygous deletions have been implicated in tumorigenesis54. TLE3 is a transcriptional corepressor involved in tumorigenesis and immune function55. The transcription factor TULP3 has been implicated in pancreatic ductal adenocarcinoma and colorectal cancer56. XCL2, NMUR2, and TULP3 may also be related to cancer progression via G-protein-coupled receptor (GPCR) signaling pathways57. JADE2 expression has been experimentally linked to NSCLC58, and has been identified in GWAS of smoking behavior34. Finally, DNA damage repair mechanisms emerge, including RPAP3, an RNA polymerase that may be involved in DNA damage repair regulation59, and MMS22L which repairs double strand breaks60.

Although smoking is the major risk factor for lung cancer, it is important to clearly disentangle the effect of smoking to fully understand the complex genetic and environmental causes of lung cancer. Our approach enables the development of new polygenic scores, which can improve precision medicine applications for lung cancer in both smokers and nonsmokers.

Author contributions statement

Drafted the manuscript: B.R.G., M.F., S.-G. J., A.K.S., E.P., A.K.D., S.P.

Acquired the data: B.R.G., S.-G. J., A.K.S., Y.S., P.D., U.S., D.D.S., W.T., J.M., S.M., R.R., R.J.H., J.D.M., Y.B., C.I.A., S.P.

Analyzed the data: B.R.G., S.-G. J., M.F., A.K.S., Y.S., P.D., U.S., Y.B., R.S.

Critically revised the manuscript for important intellectual content: all authors.

Subject terms and techniques

Biological sciences > Cancer > Lung cancer

Biological sciences > Genetics > Genetic association study > Genome-wide association studies

Data Availability

The full summary level association data from the individual population analyses in MVP will be available upon publication via the dbGaP study accession number phs001672.

Competing interests

S.-G.J. is an employee and shareholder of BridgeBio Pharma. The other authors declare no competing interests.

Methods

Cohort definition

Patients were identified from MVP participants19 utilizing clinical information available through the United States Department of Veterans Affairs (VA) Corporate Data Warehouse (CDW) with ICD codes for primary lung cancer. Occurrences of the ICD-9 codes 162.3, 162.4, 162.5, 162.8, and 162.9 or the ICD-10 codes C34.10, C34.11, C34.12, C34.2, C34.30, C34.31, C34.32, C34.80, C34.81, C34.82, C34.90, C34.91, and C34.92 were used in case identification. Patients with secondary lung cancer were excluded from the cohort using ICD-9/10 codes 197.x, C78.00, C78.01, and C78.02. Additional patients were identified in the VA Cancer Registry using ICD-O site, including lung/bronchus, other respiratory system or intrathoracic organs, or trachea. The Cancer Registry was also used to determine the lung cancer subtypes LUAD and LUSC among cases.

Preliminary totals of 18,633 and 10,845 patients with MVP participation were identified from the VA CDW and Cancer Registry, respectively. A combined cohort of 20,631 unique patients was generated for further analysis. The cohort was predominantly male (∼95%) with a median age of 64–68 for sub-cohorts, depending on ancestry assignments and cancer subtypes. The cohort was curated further to remove any participant with missing data. The final cohorts are described in Supplementary Data 1.

Once patients were identified from VA’s CDW and Cancer Registry, cases were used to gather records related to age, sex, smoking status, and ancestry. Smoking status included former, current, and never, based on the MVP survey at the time of enrollment and on electronic medical records. Ancestry was defined using a machine learning algorithm that harmonizes self-reported ethnicity and genetic ancestry (HARE)61. All analyses described here were performed on patients of EA or AA ancestry in ancestry-stratified cohorts. Additionally, the cohorts were further stratified by lung cancer subtypes for analysis. Matched controls were selected based on age, gender, smoking status, and HARE assignments. Age was binned into 5-year intervals for this purpose.

Array genotyping, genotype quality control, and principal component analysis

Genotyping and quality control were conducted as described previously62. Briefly, we removed all samples with excess heterozygosity (F statistic<-0.1), excess relatedness (kinship coefficient≥0.1 with 7 or more MVP samples), and samples with call rates <98.5%. Additional samples with a mismatch between self-reported sex and genetic sex were removed.

Principal component (PC) analysis was conducted using PLINK 2.063 (v2.00a3LM), on a pruned set of SNPs (window size 1Mb, step size 80, r2<0.1, minor allele frequency (MAF)<0.01, Hardy-Weinberg equilibrium P<1×10−10, missingness rate<10%) within European ancestry (EA) and African ancestry (AA) on unrelated individuals, where unrelated individuals were defined as greater than third-degree relatives as previously described62. PCs were then projected onto related individuals in EA.

Imputation

Prior to imputation, a within-cohort pre-phasing procedure was applied across the whole cohort by chromosome using Eagle264. Imputation was then conducted on pre-phased genotypes using Minimac465 and the 1000 Genomes Phase 3 (v5) reference panel66 in 20Mb chunks and 3Mb flanking regions. Quality of imputation (Minimac Rsq or INFO) was then re-computed in EA and AA separately to be used as filters for respective GWAS. Imputed loci reaching genome-wide significance were tested for deviation from Hardy-Weinberg equilibrium (HWE) in 61,538 EA controls (Supplementary Data 23). Of the 93 conditionally independent SNPs across the GWAS analyses, 6 SNPs had a significant (P < 1×10−6) HWE signal; unsurprisingly, the strongest HWE signal was from SNPs in the Major Histocompatibility Complex region. However, none of the 12 novel loci reported in Table 1 significantly deviated from HWE.

Association analyses

For the EA lung cancer overall and subtype GWAS, we performed standard logistic regression using PLINK 2.0 (v2.00a2LM)63 with a matched control design. EA GWAS was performed in unrelated individuals, defined as greater than third-degree relatives. For the AA lung cancer overall and subtype analyses, because the case numbers were smaller, we performed a mixed-model logistic regression using REGENIE (v1.0.6.7)67; REGENIE applies a whole genome regression model to control for relatedness and population structure, and includes a Firth correction to control for bias in rare SNPs as well as case-control imbalance. GWAS covariates for each ancestry included age, age-squared, sex, and smoking status as a categorical variable (current, former, never), and the first ten principal components. Participants with missing smoking status (n=786) were removed.

EA meta-analysis

We performed inverse-variance weighted meta-analyses of MVP-EA summary statistics and summary statistics previously reported by ILCCO7 using METAL (v20100505)68 with scheme STDERR. Significant inflation across GWAS and meta-analyses was not observed (all genomic control values (λ) for GWAS in this study ≤1.15). Only variants present in both studies were meta-analyzed. We further performed a sensitivity analysis using the Han-Eskin random effects model (RE2) in METASOFT v2.0.1 33.

Lung eQTL consortium

The lung tissues used for eQTL analyses were from human subjects who underwent lung surgery at three academic sites: Laval University, University of British Columbia (UBC), and University of Groningen. Genotyping was carried out using the Illumina Human1M-Duo BeadChip. Expression profiling was performed using an Affymetrix custom array (see GEO platform GPL10379). Only samples that passed genotyping and gene expression quality controls were considered for eQTL analysis, leaving sample sizes of 409 for Laval, 287 for UBC, and 342 for Groningen. Within each set, genotypes were imputed in each cohort with the Michigan Imputation Server65 using the Haplotype Reference Consortium69 version 1 (HRC.r1-1) data as a reference set, and gene expression values were adjusted for age, sex, and smoking status. Normalized gene expression values from each set were then combined with ComBat70. eQTLs were calculated using a linear regression model and additive genotype effects as implemented in the Matrix eQTL package in R71. Cis-eQTLs were defined by a 2 Mb window, i.e., 1 Mb distance on either side of lung cancer-associated SNPs. Pre-computed lung eQTLs were also obtained from the Genotype-Tissue Expression (GTEx) Portal20. Lung eQTLs in GTEx (version 8) are based on 515 individuals and calculated using FastQTL72.

Fine-mapping

We performed Bayesian fine-mapping the genome-wide significant loci from EA meta-analysis and AA using the FinnGen fine-mapping pipeline73 (https://github.com/FINNGEN/finemapping-pipeline) and SuSiE25,26. Pairwise SNP correlations were calculated directly from imputed dosages on European-ancestry MVP samples from this analysis using LDSTORE 2.073. The maximum number of allowed causal SNPs at each locus was set to 10. Fine-mapping regions which overlapped the major histocompatibility complex (MHC; chr6:25,000,000-34,000,000) were excluded. High quality credible sets were defined as those with minimum r2<0.5 between variants. The functional consequences of the AA credible set variants were annotated using the Variant Effect Predictor (VEP)31.

Replication analysis

External replication was performed for all genome-wide significant associations in overall lung cancer, LUAD, and LUSC in OncoArray Consortium Lung Study (OncoArray)8,74. Replication for genome-wide significant multi-ancestry associations was performed in a fixed effects meta-analysis of OncoArray CEU Europeans for significant EA meta-analysis associations, and in an YRI AA meta-analysis composed of 5 studies8 for significant MVP AA associations. Meta-analysis associations from this study were replicated against a meta-analysis of these OncoArray groups. To replicate significant variants from EA analysis conditioned on smoking, pack-years was additionally included as a covariate in replication cohorts. There was no participant overlap between the replication cohorts and the ILCCO study7 used in the discovery scan. Covariates included the first five genetic principal components and participant study sites. Proxy SNPs were used to replicate known associations at rs75675343 (rs2318539/4:67831628:C:A; R2 =1) and rs4586884 (rs4435699/4:164019500:C:G; R2 =0.999).

Multi-ancestry meta-analysis

A multi-ancestry meta-analysis of MVP EA and AA cohorts with summary statistics previously reported by ILCCO7 was conducted in METAL68 using an inverse variance-weighted fixed effects scheme. Only variants present in two or more cohorts were meta-analyzed. Index variants were defined using the two-stage “clumping” procedure implemented in the Functional Mapping and Annotation (FUMA) platform75. In this process, genome-wide significant variants are collapsed into LD blocks (r2>0.6) and subsequently re-clumped to yield approximately independent (r2<0.1) signals; adjacent signals separated by <250kb are ligated to form independent loci. Novel variants are defined as meta-analysis index variants located >1Mb from previously reported lung cancer associations. We additionally performed a sensitivity analysis using the random effects model (RE2) in METASOFT v2.0.133.

Polygenic risk score (PRS) calculation

We used PRS-CS76 to generate effect size estimates under a Bayesian shrinkage framework, and then used PLINK 2.0 (v2.00a3LM)63 to linearly combine weights into a risk score using a global shrinkage prior of 1×10−4, which is recommended for less polygenic traits. Finally, scores were normalized to a mean of 0 and a standard deviation of 1.

Multi-trait analyses

In order to remove all residual effects of smoking on lung cancer susceptibility, we conducted a multi-trait meta-analysis35 conditioned on cigarettes per day, which was shown to be most significantly correlated with all lung cancer GWAS34. The meta-analysis was performed on the EA meta-analysis summary statistics using mtCOJO, part of the GCTA software package77. An LD reference was constructed from 50,000 MVP EA samples.

Multi-trait analysis of GWAS (MTAG)45 (v0.9.0) was applied using genome-wide LUSC summary statistics after conditioning on cigarettes per day, and estrogen receptor negative (ER−) breast cancer summary statistics46 which were munged using LDSC (v1.01)38. Single causal variant colocalization between LUSC conditioned on cigarettes per day and ER-breast cancer was performed using Coloc (R; version 4)78 for variants at ABHD8 (chr19: 17,350,000 to 17,475,000). A posterior probability > 0.9 for Hypothesis 4 (both traits are associated and share a single causal variant) was used as the criteria for colocalization.

Heritability and genetic correlations

Linkage Disequilibrium score regression (LDSC) v1.0.1 was used to calculate observed-scale SNP-heritability38 using lung cancer and subtypes summary statistics, before and after conditioning on cigarettes per day. Pairwise genetic correlations were estimated between lung cancer and subtypes from MVP, ILCCO7, and EA meta-analysis, and four smoking traits (smoking initiation, cigarettes per day, smoking cessation, and age of initiation)34.

Conditional and joint SNP analysis

To find independently associated genome-wide significant SNPs at each locus in a stepwise fashion, we used GCTA-COJO using the --cojo-slct option. An LD reference was constructed from 50,000 MVP EA samples. Variants with MAF<0.01 in the COJO reference panel were not included in identification of independent signals. LDTrait79 was queried to identify previously published significant GWAS variants within 1Mb of our index variants in all populations. Novel loci were defined as those at which the index variant was not within ±500_kb of previously reported genome-wide significant lead SNPs for lung cancer or its subtypes in any ancestry.

Phenome-wide association study (PheWAS)

We conducted a PheWAS of electronic health record-derived phenotypes and lab results in EA subjects using either the normalized PRS as the predictor or independently associated genome-wide significant SNPs. Comparison of unconditioned PRS PheWAS and conditioned PRS PheWAS were based on ILCCO summary statistics7 and used MVP EA as the out-of-sample test set. Associations were tested using the R PheWAS package80 version 0.1 with QC procedures described previously81. Control and sex-based exclusion criteria were applied.

Supplementary Figure captions

Supplementary Fig. 1. Study overview. Genome-wide association studies were performed in Million Veteran Program (MVP) European and African ancestry (AA) cohorts for overall lung cancer, adenocarcinoma, and squamous cell carcinoma. MVP and International Lung Cancer Consortium OncoArray (ILCCO) European cohorts were meta-analyzed, and further meta-analyzed with AA for multi-ancestry meta-analysis. Multi-trait conditional meta-analysis was performed on EA using average cigarettes per day from Liu et al. (2019). Replication and combined meta-analysis was performed using external OncoArray cohorts.

Supplementary Fig. 2. Manhattan plots and quantile-quantile (QQ) plots for European meta-analyses. Manhattan and QQ plots are shown for a) overall lung cancer; b) lung adenocarcinoma (LUAD); and c) squamous cell lung carcinoma (LUSC). Cytoband positions for significant loci are noted in each Manhattan plot; putatively novel loci identified in this study are in red; externally replicated novel loci are indicated by a box. Genomic control (λ) values, LDSC intercepts, and sample sizes are inset in QQ plots.

Supplementary Fig. 3. Effect allele frequency concordance between International Lung Cancer Consortium OncoArray (ILCCO) and Million Veteran Program European ancestry (EA) GWAS. (a-c) Effect allele frequency concordance for all variants tested in both studies with P<1×10−5 in ILCCO for a) overall lung cancer, b) lung adenocarcinoma, and c) squamous cell lung carcinoma. Points are styled based on significance level in MVP. (d-f) Effect size concordance for genome-wide significant variants in d) overall lung cancer, e) lung adenocarcinoma, and f) squamous cell lung carcinoma. One-to-one concordance is shown as a dashed line. Index variants from the EA meta-analysis between ILCCO and MVP are annotated by locus. Novel significant loci after meta-analysis are annotated in red.

Supplementary Fig. 4. Genome-wide significant novel lung cancer loci. Forest plots (left) and regional Manhattan plots (right) for novel loci from European meta-analysis: a) XCL2, b) LSAMP, c) NMUR2, d) TUPL3, e) MYC, f) TLE3, and g) BLOC1S2; and from multi-ancestry meta-analysis: h) JADE2; i) RPAP3.

Supplementary Fig. 5. Manhattan plots and quantile-quantile (QQ) plots for MVP African ancestry. Manhattan and QQ plots are shown for a) African ancestry overall lung cancer; b) lung adenocarcinoma (LUAD); and c) squamous cell lung carcinoma (LUSC). Cytoband positions for significant loci are noted in each Manhattan plot; putatively novel loci identified in this study are in red. Genomic control (λ) values and sample sizes are inset in QQ plots.

Supplementary Fig. 6. Manhattan plots and quantile-quantile (QQ) plots for multi-ancestry meta-analyses. Manhattan and QQ plots are shown for a) the multi-ancestry meta-analysis in overall lung cancer; b) lung adenocarcinoma (LUAD); and c) squamous cell lung carcinoma (LUSC). Cytoband positions for significant loci are noted in each Manhattan plot; novel loci not identified in the European meta-analysis are in red; externally replicated novel loci are indicated by a box. Genomic control (λ) values and sample sizes are inset in QQ plots.

Supplementary Fig. 7. Association of the lung cancer polygenic risk score (PRS) with lung cancer by smoking status. a) Association of the lung cancer PRS with overall lung cancer risk. The risk of lung cancer reached an odds ratio (OR) of 2.51 (95% confidence interval: 1.80, 3.51) in the top decile. b) Association of the lung cancer PRS with lung cancer risk in never-smokers. Among never-smokers, lung cancer risk reached an OR of 2.67 (2.40, 2.98) in the top decile. c) Association of the lung cancer PRS with lung cancer risk in ever-smokers with no history of lung cancer. The top PRS decile was associated with an OR of 1.25 (1.18, 1.32).

Supplementary Fig. 8. Manhattan plots and quantile-quantile (QQ) plots for European meta-analyses conditioned on cigarettes per day. Manhattan and QQ plots for a) overall lung cancer conditioned on cigarettes per day; b) lung adenocarcinoma (LUAD) conditioned on cigarettes per day; and c) squamous cell lung carcinoma (LUSC) conditioned on cigarettes per day. Cytoband positions for significant loci are noted in each Manhattan plot; novel loci not identified in the European meta-analysis are in red; externally replicated novel loci are indicated by a box. Genomic control (λ) values, LDSC intercepts, and sample sizes are inset in QQ plots.

Supplementary Fig. 9. Novel loci for overall lung cancer and squamous cell carcinoma conditioned on smoking. Forest plots (left) and regional Manhattan plots (right) for novel loci identified in the European meta-analysis conditioned on cigarettes per day: a) MMS22L in overall lung cancer and b) ABHD8 in squamous cell lung cancer.

Acknowledgements

This work was supported by award #MVP000 from the United States Department of Veterans Affairs (VA) Million Veteran Program. The contents of this publication are the sole responsibility of the authors and do not necessarily represent the views of VA or the United States Government. Where authors 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. Full consortium acknowledgements for MVP and the ILCCO OncoArray study7 are provided in Supplementary Information.

References

  1. 1.↵
    Schabath, M. B. & Cote, M. L. Cancer Progress and Priorities: Lung Cancer. Cancer Epidemiol. Biomarkers Prev. 28, 1563–1579 (2019).
    OpenUrlFREE Full Text
  2. 2.↵
    Leiter, A., Veluswamy, R. R. & Wisnivesky, J. P. The global burden of lung cancer: current status and future trends. Nat. Rev. Clin. Oncol. 20, 624–639 (2023).
    OpenUrlCrossRef
  3. 3.↵
    Sung, H. et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 71, 209–249 (2021).
    OpenUrlCrossRefPubMed
  4. 4.↵
    Siegel, R. L., Miller, K. D., Fuchs, H. E. & Jemal, A. Cancer statistics, 2022. CA Cancer J. Clin. 72, 7–33 (2022).
    OpenUrlCrossRefPubMed
  5. 5.↵
    Bossé, Y. & Amos, C. I. A Decade of GWAS Results in Lung Cancer. Cancer Epidemiol. Biomarkers Prev. 27, 363–379 (2018).
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    Timofeeva, M. N. et al. Influence of common genetic variation on lung cancer risk: meta-analysis of 14 900 cases and 29 485 controls. Hum. Mol. Genet. 21, 4980–4995 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    McKay, J. D. et al. Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat. Genet. 49, 1126–1132 (2017).
    OpenUrlPubMed
  8. 8.↵
    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. 54, 1167–1177 (2022).
    OpenUrl
  9. 9.↵
    Wang, Y. et al. SNP rs17079281 decreases lung cancer risk through creating an YY1-binding site to suppress DCBLD1 expression. Oncogene 39, 4092–4102 (2020).
    OpenUrlCrossRef
  10. 10.↵
    Zhang, T. et al. Genomic and evolutionary classification of lung cancer in never smokers. Nat. Genet. 53, 1348–1359 (2021).
    OpenUrlCrossRefPubMed
  11. 11.↵
    Govindan, R. et al. Genomic landscape of non-small cell lung cancer in smokers and never-smokers. Cell 150, 1121–1134 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.↵
    Wang, Z. et al. Meta-analysis of genome-wide association studies identifies multiple lung cancer susceptibility loci in never-smoking Asian women. Hum. Mol. Genet. 25, 620–629 (2016).
    OpenUrlCrossRefPubMed
  13. 13.↵
    Schabath, M. B., Cress, D. & Munoz-Antonia, T. Racial and Ethnic Differences in the Epidemiology and Genomics of Lung Cancer. Cancer Control 23, 338–346 (2016).
    OpenUrlCrossRef
  14. 14.↵
    Long, E., Patel, H., Byun, J., Amos, C. I. & Choi, J. Functional studies of lung cancer GWAS beyond association. Hum. Mol. Genet. 31, R22–R36 (2022).
    OpenUrl
  15. 15.↵
    Shi, J. et al. Genome-wide association study of lung adenocarcinoma in East Asia and comparison with a European population. Nat. Commun. 14, 3043 (2023).
    OpenUrl
  16. 16.↵
    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 7, 881–891 (2019).
    OpenUrl
  17. 17.↵
    Nahar, R. et al. Elucidating the genomic architecture of Asian EGFR-mutant lung adenocarcinoma through multi-region exome sequencing. Nat. Commun. 9, 216 (2018).
    OpenUrlCrossRef
  18. 18.↵
    Zanetti, K. A. et al. Genome-wide association study confirms lung cancer susceptibility loci on chromosomes 5p15 and 15q25 in an African-American population. Lung Cancer 98, 33–42 (2016).
    OpenUrlCrossRef
  19. 19.↵
    Gaziano, J. M. et al. Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J. Clin. Epidemiol. 70, 214–223 (2016).
    OpenUrlCrossRefPubMed
  20. 20.↵
    GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 369, 1318–1330 (2020).
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    Hao, K. et al. Lung eQTLs to help reveal the molecular underpinnings of asthma. PLoS Genet. 8, e1003029 (2012).
    OpenUrlCrossRefPubMed
  22. 22.↵
    Bossé, Y. et al. Transcriptome-wide association study reveals candidate causal genes for lung cancer. Int. J. Cancer 146, 1862–1878 (2020).
    OpenUrl
  23. 23.↵
    Koutsami, M. K. et al. Centrosome abnormalities are frequently observed in non-small-cell lung cancer and are associated with aneuploidy and cyclin E overexpression. J. Pathol. 209, 512–521 (2006).
    OpenUrlCrossRefPubMed
  24. 24.↵
    Chan, J. Y. A clinical overview of centrosome amplification in human cancers. Int. J. Biol. Sci. 7, 1122–1144 (2011).
    OpenUrlCrossRefPubMed
  25. 25.↵
    Wang, G., Sarkar, A., Carbonetto, P. & Stephens, M. A simple new approach to variable selection in regression, with application to genetic fine mapping. J. R. Stat. Soc. Series B Stat. Methodol. 82, 1273–1300 (2020).
    OpenUrl
  26. 26.↵
    Zou, Y., Carbonetto, P., Wang, G. & Stephens, M. Fine-mapping from summary data with the ‘Sum of Single Effects’ model. PLoS Genet. 18, e1010299 (2022).
    OpenUrl
  27. 27.↵
    de Goede, O. M. et al. Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease. Cell 184, 2633–2648.e19 (2021).
    OpenUrl
  28. 28.
    Li, Y. et al. Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers. Nat. Commun. 11, 1000 (2020).
    OpenUrlCrossRefPubMed
  29. 29.
    de Santiago, P. R. et al. Immune-related IncRNA LINC00944 responds to variations in ADAR1 levels and it is associated with breast cancer prognosis. Life Sci. 268, 118956 (2021).
    OpenUrl
  30. 30.↵
    Chen, D. et al. Genome-wide analysis of long noncoding RNA (lncRNA) expression in colorectal cancer tissues from patients with liver metastasis. Cancer Med. 5, 1629–1639 (2016).
    OpenUrl
  31. 31.↵
    McLaren, W. et al. The Ensembl Variant Effect Predictor. Genome Biol. 17, 122 (2016).
    OpenUrlCrossRefPubMed
  32. 32.↵
    Saunders, G. R. B. et al. Genetic diversity fuels gene discovery for tobacco and alcohol use. Nature 612, 720–724 (2022).
    OpenUrlCrossRef
  33. 33.↵
    Han, B. & Eskin, E. Random-effects model aimed at discovering associations in meta-analysis of genome-wide association studies. Am. J. Hum. Genet. 88, 586–598 (2011).
    OpenUrlCrossRefPubMed
  34. 34.↵
    Liu, M. et al. Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat. Genet. 51, 237–244 (2019).
    OpenUrlCrossRefPubMed
  35. 35.↵
    Zhu, Z. et al. Causal associations between risk factors and common diseases inferred from GWAS summary data. Nat. Commun. 9, 1–12 (2018).
    OpenUrlCrossRefPubMed
  36. 36.↵
    Xue, A. et al. Genome-wide analyses of behavioural traits are subject to bias by misreports and longitudinal changes. Nat. Commun. 12, 20211 (2021).
    OpenUrl
  37. 37.↵
    Munafò, M. R., Tilling, K., Taylor, A. E., Evans, D. M. & Davey Smith, G. Collider scope: when selection bias can substantially influence observed associations. Int. J. Epidemiol. 47, 226–235 (2018).
    OpenUrlCrossRefPubMed
  38. 38.↵
    Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295 (2015).
    OpenUrlCrossRefPubMed
  39. 39.↵
    Nguyen, M.-H., Ueda, K., Nakamura, Y. & Daigo, Y. Identification of a novel oncogene, MMS22L, involved in lung and esophageal carcinogenesis. Int. J. Oncol. 41, 1285–1296 (2012).
    OpenUrlPubMed
  40. 40.↵
    Yang, W. et al. Deciphering associations between three RNA splicing-related genetic variants and lung cancer risk. NPJ Precis Oncol 6, 48 (2022).
    OpenUrl
  41. 41.↵
    Gabriel, A. A. G. et al. Genetic Analysis of Lung Cancer and the Germline Impact on Somatic Mutation Burden. J. Natl. Cancer Inst. 114, 1159–1166 (2022).
    OpenUrl
  42. 42.↵
    Lawrenson, K. et al. Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus. Nat. Commun. 7, 12675 (2016).
    OpenUrlCrossRefPubMed
  43. 43.↵
    Lesseur, C. et al. Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers. PLoS Genet. 17, e1009254 (2021).
    OpenUrlCrossRef
  44. 44.↵
    Fachal, L. et al. Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes. Nat. Genet. 52, 56–73 (2020).
    OpenUrlCrossRef
  45. 45.↵
    Turley, P. et al. Multi-trait analysis of genome-wide association summary statistics using MTAG. Nat. Genet. 50, 229–237 (2018).
    OpenUrlCrossRefPubMed
  46. 46.↵
    Michailidou, K. et al. Association analysis identifies 65 new breast cancer risk loci. Nature 551, 92–94 (2017).
    OpenUrlCrossRefPubMed
  47. 47.↵
    Wallace, C. A more accurate method for colocalisation analysis allowing for multiple causal variants. PLoS Genet. 17, e1009440 (2021).
    OpenUrlCrossRefPubMed
  48. 48.↵
    Gusev, A. et al. A transcriptome-wide association study of high-grade serous epithelial ovarian cancer identifies new susceptibility genes and splice variants. Nat. Genet. 51, 815–823 (2019).
    OpenUrlCrossRef
  49. 49.↵
    Brenner, D. R. et al. Alcohol consumption and lung cancer risk: A pooled analysis from the International Lung Cancer Consortium and the SYNERGY study. Cancer Epidemiol. 58, 25–32 (2019).
    OpenUrl
  50. 50.↵
    Larsson, S. C. et al. Smoking, alcohol consumption, and cancer: A mendelian randomisation study in UK Biobank and international genetic consortia participants. PLoS Med. 17, e1003178 (2020).
    OpenUrlPubMed
  51. 51.↵
    Petrelli, F. et al. Association of Obesity With Survival Outcomes in Patients With Cancer: A Systematic Review and Meta-analysis. JAMA Netw Open 4, e213520 (2021).
    OpenUrl
  52. 52.↵
    Lan, T., Chen, L. & Wei, X. Inflammatory Cytokines in Cancer: Comprehensive Understanding and Clinical Progress in Gene Therapy. Cells 10, (2021).
  53. 53.↵
    Kresse, S. H. et al. LSAMP, a novel candidate tumor suppressor gene in human osteosarcomas, identified by array comparative genomic hybridization. Genes Chromosomes Cancer 48, 679–693 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    Xie, J. et al. Copy number analysis identifies tumor suppressive lncRNAs in human osteosarcoma. Int. J. Oncol. 50, 863–872 (2017).
    OpenUrl
  55. 55.↵
    Yu, G. et al. Roles of transducin-like enhancer of split (TLE) family proteins in tumorigenesis and immune regulation. Front Cell Dev Biol 10, 1010639 (2022).
    OpenUrl
  56. 56.↵
    Sartor, I. T. S., Recamonde-Mendoza, M. & Ashton-Prolla, P. TULP3: A potential biomarker in colorectal cancer? PLoS One 14, e0210762 (2019).
    OpenUrl
  57. 57.↵
    Chaudhary, P. K. & Kim, S. An Insight into GPCR and G-Proteins as Cancer Drivers. Cells 10, (2021).
  58. 58.↵
    Murphy, C. et al. An Analysis of JADE2 in Non-Small Cell Lung Cancer (NSCLC). Biomedicines 11, (2023).
  59. 59.↵
    Ni, L. et al. RPAP3 interacts with Reptin to regulate UV-induced phosphorylation of H2AX and DNA damage. J. Cell. Biochem. 106, 920–928 (2009).
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    Saredi, G. et al. H4K20me0 marks post-replicative chromatin and recruits the TONSL– MMS22L DNA repair complex. Nature 534, 714–718 (2016).
    OpenUrlCrossRefPubMed
  61. 61.↵
    Fang, H. et al. Harmonizing Genetic Ancestry and Self-identified Race/Ethnicity in Genome-wide Association Studies. Am. J. Hum. Genet. 105, 763–772 (2019).
    OpenUrlCrossRef
  62. 62.↵
    Hunter-Zinck, H. et al. Genotyping Array Design and Data Quality Control in the Million Veteran Program. Am. J. Hum. Genet. 106, 535–548 (2020).
    OpenUrlCrossRef
  63. 63.↵
    Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7 (2015).
    OpenUrlCrossRefPubMed
  64. 64.↵
    Loh, P.-R. et al. Reference-based phasing using the Haplotype Reference Consortium panel. Nat. Genet. 48, 1443–1448 (2016).
    OpenUrlCrossRefPubMed
  65. 65.↵
    Das, S. et al. Next-generation genotype imputation service and methods. Nat. Genet. 48, 1284–1287 (2016).
    OpenUrlCrossRefPubMed
  66. 66.↵
    1000 Genomes Project Consortium et al. A global reference for human genetic variation. Nature 526, 68–74 (2015).
    OpenUrlCrossRefPubMed
  67. 67.↵
    Mbatchou, J. et al. Computationally efficient whole-genome regression for quantitative and binary traits. Nat. Genet. 53, 1097–1103 (2021).
    OpenUrlCrossRefPubMed
  68. 68.↵
    Willer, C. J., Li, Y. & Abecasis, G. R. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26, 2190–2191 (2010).
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    McCarthy, S. et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat. Genet. 48, 1279–1283 (2016).
    OpenUrlCrossRefPubMed
  70. 70.↵
    Johnson, W. E., Li, C. & Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8, 118–127 (2007).
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    Shabalin, A. A. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics 28, 1353–1358 (2012).
    OpenUrlCrossRefPubMedWeb of Science
  72. 72.↵
    Ongen, H., Buil, A., Brown, A. A., Dermitzakis, E. T. & Delaneau, O. Fast and efficient QTL mapper for thousands of molecular phenotypes. Bioinformatics 32, 1479–1485 (2016).
    OpenUrlCrossRefPubMed
  73. 73.↵
    Benner, C. et al. FINEMAP: efficient variable selection using summary data from genome-wide association studies. Bioinformatics 32, 1493–1501 (2016).
    OpenUrlCrossRefPubMed
  74. 74.↵
    Amos, C. I. et al. The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers. Cancer Epidemiol. Biomarkers Prev. 26, 126–135 (2017).
    OpenUrlAbstract/FREE Full Text
  75. 75.↵
    Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826 (2017).
    OpenUrlCrossRefPubMed
  76. 76.↵
    Ge, T., Chen, C.-Y., Ni, Y., Feng, Y.-C. A. & Smoller, J. W. Polygenic prediction via Bayesian regression and continuous shrinkage priors. Nat. Commun. 10, 1776 (2019).
    OpenUrl
  77. 77.↵
    Yang, J., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82 (2011).
    OpenUrlCrossRefPubMed
  78. 78.↵
    Giambartolomei, C. et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 10, e1004383 (2014).
    OpenUrlCrossRefPubMed
  79. 79.↵
    Machiela, M. J. & Chanock, S. J. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 31, 3555–3557 (2015).
    OpenUrlCrossRefPubMed
  80. 80.↵
    Carroll, R. J., Bastarache, L. & Denny, J. C. R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment. Bioinformatics 30, 2375–2376 (2014).
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    Klarin, D. et al. Genetics of blood lipids among ∼300,000 multi-ethnic participants of the Million Veteran Program. Nat. Genet. 50, 1514–1523 (2018).
    OpenUrlCrossRefPubMed
Back to top
Posted April 26, 2024.
Download PDF

Supplementary Material

Data/Code
Email
Multi-ancestry meta-analyses of lung cancer in the Million Veteran Program reveal novel risk loci and elucidate smoking-independent genetic risk
Bryan R. Gorman, Sun-Gou Ji, Michael Francis, Anoop K. Sendamarai, Yunling Shi, Poornima Devineni, Uma Saxena, Elizabeth Partan, Andrea K. DeVito, Jinyoung Byun, Younghun Han, Xiangjun Xiao, Don D. Sin, Wim Timens, Jennifer Moser, Sumitra Muralidhar, Rachel Ramoni, Rayjean J. Hung, James D. McKay, Yohan Bossé, Ryan Sun, Christopher I. Amos, VA Million Veteran Program, Saiju Pyarajan
medRxiv 2024.04.25.24306313; doi: https://doi.org/10.1101/2024.04.25.24306313
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area