1. Epidemiology and Global Health

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome

  1. Ricardo Cortez Cardoso Penha  Is a corresponding author
  2. Karl Smith-Byrne
  3. Joshua R Atkins
  4. Philip C Haycock
  5. Siddhartha Kar
  6. Veryan Codd
  7. Nilesh J Samani
  8. Christopher Nelson
  9. Maja Milojevic
  10. Aurélie AG Gabriel
  11. Christopher Amos
  12. Paul Brennan
  13. Rayjean J Hung
  14. Linda Kachuri
  15. James D Mckay  Is a corresponding author
  1. Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), France
  2. Cancer Epidemiology Unit, University of Oxford, United Kingdom
  3. MRC Integrative Epidemiology Unit, Bristol Population Health Science Institute, Bristol Medical School (PHS), United Kingdom
  4. Department of Cardiovascular Sciences, University of Leicester, United Kingdom
  5. NIHR Leicester Biomedical Research Centre, Glenfield Hospital, United Kingdom
  6. Ludwig Lausanne Branch, Faculty of Biology and Medicine, Switzerland
  7. Institute for Clinical and Translational Research, Baylor College of Medicine, United States
  8. Lunenfeld-Tanenbaum Research Institute, Sinai Health, Canada
  9. Departament of Epidemiology and Population Health, Stanford University, United States
https://doi.org/10.7554/eLife.83118
Share this article
Cite this article
  1. Ricardo Cortez Cardoso Penha
  2. Karl Smith-Byrne
  3. Joshua R Atkins
  4. Philip C Haycock
  5. Siddhartha Kar
  6. Veryan Codd
  7. Nilesh J Samani
  8. Christopher Nelson
  9. Maja Milojevic
  10. Aurélie AG Gabriel
  11. Christopher Amos
  12. Paul Brennan
  13. Rayjean J Hung
  14. Linda Kachuri
  15. James D Mckay
(2023)
Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome
eLife 12:e83118.
https://doi.org/10.7554/eLife.83118
  2. Download BibTeX
  3. Download .RIS

Abstract

Background:

Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours.

Methods:

We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343).

Results:

Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity.

Conclusions:

This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas.

Funding:

Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17‐022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

Editor's evaluation

This study is of interest to epidemiologists and geneticists studying the association between telomere length and lung cancer risk. This work provides useful insight into risk factors for lung cancer. Overall, the results of this study are solid, as the genetic instrument used here is better powered and the battery of Mendelian randomization analysis makes this broad set of results convincing compared to previous work.

https://doi.org/10.7554/eLife.83118.sa0

Introduction

Telomeres are a complex of repetitive TTAGGG sequences and nucleoproteins located at the end of chromosomes and have an essential role in sustaining cell proliferation and preserving genome integrity (de Lange, 2009). Telomere length progressively shortens with age in proliferative somatic cells due to incomplete telomeric regions replication (Watson, 1972) and low activity of the telomerase TERT in adult cells. The shortening of the telomere length results in cell cycle arrest, cellular senescence, and apoptosis in somatic cells (Harley et al., 1990). The maintenance of telomere length, which allows cancer cells to escape the telomere-mediated cell death pathways, is one feature related to the hallmarks of cancer (Hanahan and Weinberg, 2011).

Telomere length appears to vary between individuals and has been studied in relation to many diseases. In observational studies, telomere length is measured as the average length of telomeric sequences in a given tissue (Montpetit et al., 2014). Telomere length appears correlated across tissue types (Demanelis et al., 2020), and as such, leukocyte telomere length (LTL) is generally measured in epidemiologic studies as a proxy for telomere length in other tissues. Recently, LTL has been measured in 472,174 individuals from the UK Biobank (UKBB; Codd et al., 2021), and LTL was associated with multiple biomedical traits (i.e. pulmonary and cardiovascular diseases, haematological traits, lymphomas, kidney cancer, and other cancer types). Genetic analysis of LTL also revealed 138 genetic loci linked to LTL across a variety of different genes involved in telomere biology and DNA repair (Codd et al., 2021).

In the context of lung cancer, genetic variants at several loci have been associated with both LTL and lung cancer risk, including variants near the TERT, TERC, OBFC1, and RTEL1 genes, fundamental to telomere length maintenance (McKay et al., 2008; Wang et al., 2008; Rafnar et al., 2009; Kachuri et al., 2016; McKay et al., 2017). The effects of the telomere-related variants appear more relevant to lung adenocarcinoma risk than other histologic subtypes (McKay et al., 2017; Landi et al., 2009). Accordingly, a causal relationship between LTL and susceptibility to lung cancer was observed using Mendelian randomisation (MR) approaches (Zhang et al., 2015; Haycock et al., 2017; Kachuri et al., 2019) as well as in observational studies that have associated directly measured telomere length with risk of lung cancer (Sanchez-Espiridion et al., 2014; Zhang et al., 2017).

The aim of the current work was to investigate the relationship between genetically predicted LTL and lung cancer, including lung cancer histological subtypes and smoking status. To this end, we conducted genome-wide correlations, MR, and colocalisation analyses to explore the relationship between LTL and lung cancer. We additionally undertook polygenic risk score (PRS) analysis using the LTL genetic instrument to explore the influence of LTL on the demographic, clinical, and molecular features of lung adenocarcinoma tumours.

Materials and methods

Reporting guidelines

Request a detailed protocol

The current study has been reported according to the STROBE-MR guidelines (Reporting Standards Document).

Data

Genome-wide association studies (GWASs) summary statistics for lung cancer (29,239 cases and 56,450 controls) and stratified by histological subtype (squamous cell carcinoma, small-cell carcinoma, and adenocarcinoma) and smoking status (ever and never smokers) were obtained from the International Lung Cancer Consortium (ILCCO; McKay et al., 2017). All analyses of LTL requiring summary statistics used results from a GWAS of LTL in 464,716 individuals of European ancestry from the UKBB (Codd et al., 2021). Downstream analyses considered additional lung cancer risk factors, such as lung function and cigarette smoking. We obtained GWAS summary statistics for forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) from a published UKBB analysis (Kachuri et al., 2020). For smoking behaviour traits, we used results from the GWAS and Sequencing Consortium of Alcohol and Nicotine use (GSCAN) consortium meta-analysis of cigarettes per day (continuous), smoking initiation (ever versus never), smoking cessation (successfully quit versus continuing), and age at smoking initiation (continuous; Liu et al., 2019) excluding the UKBB participants. For the obesity-related traits (continuous), we used the results from the UKBB and GIANT meta-analysis of BMI (Pulit et al., 2019) and waist-to-hip ratio (WHR; Pulit et al., 2019), or OpenGWAS data using UKBB participants (Elsworth et al., 2020) for high-density lipoprotein (HDL), triglycerides, and systolic and diastolic blood pressure. For the alcohol behaviour trait, we obtained the results from GSCAN phase 2 of drinks per week (continuous; Saunders et al., 2022). Colocalisation analyses of gene expression used lung tissue expression quantitative trait loci (eQTL) summary statistics from the Genotype-Tissue Expression (GTEx) data version 8.

Analyses of molecular phenotypes were performed using 343 lung adenocarcinoma samples of European ancestry from The Cancer Genome Atlas (TCGA) cohort with germline profile, RNA-sequencing, and epidemiological data available. Genotyping and imputation of germline variants have been described elsewhere (Gabriel et al., 2022). The total somatic mutation burden of TCGA samples was obtained from Ellrott et al., 2018, and DNA mutational signatures were extracted and attributed, as previously described (Gabriel et al., 2022). RNA-sequencing data were obtained from TCGA data portal using TCGA biolinks package in R (version 2.22.3; Colaprico et al., 2016). Telomere length measurement by whole-genome sequencing (WGS-measured TL, 655 samples across cancer sites) was retrieved from Barthel et al., 2017.

Tumour genomic characteristics were defined by the analyses of the TCGA data, including gene expression-based scores of telomerase activity (Barthel et al., 2017) and cellular proliferation (Thorsson et al., 2018), as well as the observed frequency of somatic homologous recombination-related events (represented as a homologous recombination repair deficiency score), and the average number of somatic copy number alteration within the tumours (Knijnenburg et al., 2018).

Linkage disequilibrium score regression

Request a detailed protocol

Genetic correlations across traits were calculated using linkage disequilibrium score regression (LDSC) by the LDSC package (v1.0.0; Bulik-Sullivan et al., 2015). Linkage disequilibrium (LD) scores were generated on the 1000 Genomes Project Phase 3 reference panel with the HLA region excluded as provided by the package due to long range LD patterns. The genome-wide correlations that passed Bonferroni correction (adjusted p-values<0.05) were considered statically significant.

Mendelian randomisation

Request a detailed protocol

MR is a method for interrogating relationships between putative risk factors and health outcomes by using genetic variants associated with the exposure of interest, typically obtained from GWAS, as instrumental variables. Assuming that fundamental MR assumptions are satisfied, this approach can be said to identify unbiased causal estimates. The genetic instrument for LTL was defined as the set of 144 genetic variants that were genome-wide significant (p<5x10–08) but not in linkage disequilibrium with each other (r2<0.01) and restricted to common genetic variation (minor allele frequency >1%) in European populations. Proxy variants in LD (r2>0.8, whenever possible) were chosen when a genetic variant was not available in the lung cancer GWAS. Primary MR analyses were conducted using the inverse-variance method with multiplicative random-effects (Yavorska and Burgess, 2017). Sensitivity analyses to horizontal pleiotropy and other violations of MR assumptions were performed using other MR estimation methods, such as weighted median, MR-Egger, contamination mixture model, MR-PRESSO, and MR-RAPS (Yavorska and Burgess, 2017; Sanderson et al., 2022). Multivariable MR (MVMR) methods included the inverse-variance weighted, MR-Egger, and least absolute shrinkage and selection operator (LASSO)-based methods (Yavorska and Burgess, 2017; Sanderson et al., 2019).

Colocalisation methods

Request a detailed protocol

Unlike MR, where the goal is to assess the evidence for a causal effect of an exposure on an outcome, colocalisation is agnostic with respect to direction of effect and only assesses the probability that the two traits are affected by the same genetic variants at a given locus. Colocalisation can be viewed as a complementary approach for evaluating MR assumptions within specific genes or regions since strong evidence of colocalisation indicates overlap in genetic mechanisms affecting LTL and lung cancer. We used COLOC (v5.1.0; Wallace, 2020) to estimate the posterior probability for two traits sharing the same causal variant (PP4) in a 150 kb LD window, with PP4 >0.70 corresponding to strong evidence of colocalisation, as previously suggested (Wallace, 2020; Lopes et al., 2022). Priors chosen for the colocalisation analyses were p1=10–3, p2=10–4, and p12=10–5, or approximately, a 75% prior belief that a signal will only be observed in the LTL GWAS and less than 0.01% prior belief in favour of colocalisation between the two traits at a given locus (Giambartolomei et al., 2014). Conditioning and masking colocalisation methods were also used as they may identify putative shared causal variants in the presence of multiple causal variants present in a defined LD window (Wallace, 2021). We present the average PP4 from all methods as our posterior belief in favour of colocalisation between LTL and lung cancer risk. Multi-trait colocalisation based on a clustering algorithm was also performed using HyPrColoc (v1.0) to identify shared genetic signals with other lung cancer-related traits (Foley et al., 2021).

Principal component analysis based on RNA-sequencing data

Request a detailed protocol

Read counts of RNA-sequencing data were normalised within (GC-content and gene length) and between (sequencing depth) lane procedures by EDASeq R package (version 2.28.0; Risso et al., 2011) and excluding low read counts. Principal component analysis was applied using singular value decomposition method, after excluding extreme outliers. Pathway analyses were conducted using Gene Set Enrichment Analysis software (GSEA, version 4.2.3; Subramanian et al., 2005) on gene annotations from Gene Ontology database. Pathway analyses were restricted to the top 500 genes positively and negatively correlated with each principal component that passed multiple-testing correction (Bonferroni-adjusted p-value<0.05 for 74,465 tests), which is the maximum number of genes supported by the online software.

The PRS for LTL was composed of the same 144 variants used in the MR analysis and was computed as the sum of the individual’s beta-weighted genotypes using PRSice-2 software (Choi and O’Reilly, 2019). Associations were estimated per SD increase in the PRS, which was normalised to have a mean of zero across lung adenocarcinoma samples of European ancestry within the TCGA cohort. The associations between the eigenvalues of the gene expression principal components (outcome) and demographic, clinical, and genomic features related to genome stability (predictors derived from TCGA published papers and TCGA data portal, except for the DNA mutational signatures [Gabriel et al., 2022]) were calculated using a multivariate linear regression model.

Inferring PC2 gene expression signature based on RNA-sequencing data

Request a detailed protocol

The TCGA-lung adenocarcinoma (LUAD) tumour samples were split into training (70%, N=255) and validation (30%, N=108) datasets. Principal components analysis based on RNA-sequencing data was performed to summarise the gene expression profiles of lung adenocarcinoma tumours into five principal components in the training and validation datasets separately, as previously described (see ‘Principal component analysis based on RNA-sequencing data’ section in methods). Subsequently, we applied the partial least squares-based method called rigid transformation (Hubert and Branden, 2003) to align the first five principal components in both datasets. This method compares embeddings, low-dimensional representations in both datasets, performing a slightly rotation of principal components in order to translate and match them in both training and validation datasets. To select the most informative genes of PC2 in the training dataset, RNA levels of the genes correlated with PC2 (N=3914 out of 14,893 genes, FDR <0.05 for 14,893 tests) were selected as variables for the LASSO regression models. The LASSO tune parameters were chosen by resampling the training dataset (1000 bootstraps: root mean of SE=0.12 ± 0.0004, Lambda = 10 x 10–10, r2=0.99 ± 0.00007) using the tidymodels metapackage in R (v1.0.0; wrapper of glmnet). The 10 genes selected by the LASSO model were used to infer the gene expression signature of PC2 by adding up the scaled values of the log-normalised read counts of each gene multiplied by the respective LASSO regression coefficients. For validation purpose, the inferred PC2 signature was calculated in the validation dataset and compared with the observed principal components. After validation in the subset of the TCGA-LUAD cohort, the inferred PC2 signature was calculated in TCGA-LUSC dataset to compare differences between lung cancer histological subtypes.

Results

Genome-wide genetic correlations

The design of the study is represented in Figure 1—figure supplement 1. We first assessed the shared genetic basis of telomere length, lung cancer risk, and other putative lung cancer risk factors, such as smoking behaviours (age start smoking, smoking cessation, smoking initiation, and cigarettes per day) and lung function (FEV1 and FVC) using genome-wide correlations (Figure 1A). There was little evidence for genetic correlations by LDSC between LTL variants and lung cancer (rg = −0.01, p=0.88) or when stratified by histologic subtypes (Figure 1A). Increasing LTL was genetically correlated with older age at smoking initiation (rg = 0.13, p=3.0 × 10–3), and negatively correlated with smoking cessation: (rg = −0.21, p=6.9 × 10–09), smoking initiation (rg = −0.16, p=1.3 × 10–10), and cigarettes per day (rg = −0.19, p=2.1 × 10–08). Longer LTL was genetically correlated with improved lung function, as indicated by increasing values of FEV1 (rg = 0.09, p=5.1 × 10–07) and FVC (rg = 0.09, p=1.1 × 10–05). To better understand the absence of genome-wide correlations between LTL and lung cancer, we visualised the Z-scores for each trait for approximately 1.2 million variants included in the LDSC analyses (Figure 1B). A subgroup of variants associated with longer LTL was correlated with increased lung adenocarcinoma risk, while the subgroup of smoking-behaviour associated variants, which also conferred an increased risk of lung adenocarcinoma, tended to have lower LTL.

Figure 1 with 1 supplement see all
Download asset Open asset
Genetic correlations between leukocyte telomere length (LTL) and lung cancer (LC) related traits.

(A) Heatmap representing the genetic correlation analyses (rg) for LTL across LC, histological subtypes (lung adenocarcinoma [ADE], squamous cell carcinoma [SQC], and small-cell carcinoma [SCC]), smoking propensity (cigarettes per day [CPD], smoking cessation [SmkCes], Smoking initiation [SmkInit], and age of smoking initiation [AgeSmk]), and lung function related (forced vital capacity [FVC] and forced expiratory volume [FEV1]) traits. The black star indicates correlations that passed Bonferroni correction (p<4x10–04). Heritability (h2) as the proportion of the phenotypic variance caused by SNPs. (B) Plot of Z-scores (ADE versus LTL), restricting to the Hapmap SNPs (~1.2 million) but excluding HLA region. Genome-wide significant SNPs (p<5x10–08) for each trait were coloured (CPD in red, SmkInit in dark red, LTL in dark blue, AgeSmk in blue, SmkCes in lightblue, and not genome-wide hits for LTL or any other selected trait in white). Linear regression line was coloured in red.

MR analyses

From the 490 genetic instruments associated with LTL at genome-wide significance (p<5x10–08), 144 LTL genetic instruments, that explained ~3.5% of the variance in LTL, and were in low-linkage disequilibrium (r2<0.01) were used in MR analysis (Supplementary file 1a). As a sensitivity analysis, a PRS composed of these genetic instruments was associated with TL estimated from WGS in blood samples across TCGA cohorts (Beta = 0.03, 95%CI = 0.01–0.05, p=0.001) but was not associated with TL in tumour material from the same patients (Figure 2—figure supplement 1).

MR analyses demonstrated that longer genetically predicted LTL was associated with increased lung cancer risk (OR = 1.62, 95% CI = 1.44–1.84, p=9.91 × 10–15) (Figure 2; Supplementary file 1). Longer LTL conferred the largest increase in risk for lung adenocarcinoma tumours (OR = 2.43, 95% CI = 2.02–2.92, p=3.76 × 10–21), but there was limited evidence of a causal relationship for other histologic subtypes, such as squamous cell carcinoma (OR = 1.00, 95% CI = 0.84–1.19, p=0.98) and small-cell carcinoma (OR = 1.13, 95% CI = 0.87–1.45, p=0.34; Figure 2, Supplementary file 1). However, our study was underpowered to detect an association between lung small-cell carcinoma and LTL at OR of 1.13 and considering alpha type-1 error rate of 5% (Figure 2—figure supplement 1). When stratifying the analyses by smoking status, LTL was associated with lung cancer risk in both never (OR = 2.02, 95% CI = 1.45–2.83, p=3.78 × 10–05) and ever smokers (OR = 1.54, 95% CI = 1.34–1.76, p=7.75 × 10–10; Figure 2, Supplementary file 1). Evidence for negative pleiotropy (supplementary file 1c) and heterogeneity (supplementary file 1d) was observed for all lung cancer outcomes except for squamous cell carcinoma. However, a significant association for LTL and lung cancer risk was found for methods robust to the significant directional pleiotropy (MR-Egger: lung cancer [OR = 2.35, p=3.37 × 10–13]; lung adenocarcinoma [OR = 4.48, p=7.30 × 10–17]; never smokers [OR = 6.84, p=2.07 × 10–10]; supplementary file 1b). Leave-one-out analyses detected only one outlier, rs7705526 in TERT, resulting in >10% change in MR effect size for associated lung cancer subtypes (supplementary file 1e). MVMR analyses considering instruments related to LTL and WHR, HDL, total triglycerides, systolic blood pressure, smoking, and alcohol intake, as well as multiple traits combined, suggested that the association between LTL and lung adenocarcinoma risk is independent of smoking propensity, obesity-related, and alcohol intake-related traits (supplementary file 1f, 1g).

Figure 2 with 1 supplement see all
Download asset Open asset
Genetically predicted leukocyte telomere length (LTL) association with lung cancer.

Lung cancer (by histology or by smoking status) risk associations with the LTL instrument from the inverse-variance-weighted MR analyses are expressed as OR per SD increase in genetically predicted LTL. Statistically significant associations with p-values<0.05 (red square). Heterogeneity is estimated by the statistic I2, tau variance of subgroups (τ2), and p-values for Cochran’s Q heterogeneity measure.

Colocalisation analyses

We investigated whether there was evidence of shared genetic signals between LTL and lung adenocarcinoma at loci centred on the 144 genetic instruments used in MR analyses using colocalisation methods (Figure 3A, supplementary file 1h). Loci with evidence of colocalisation between LTL and lung adenocarcinoma tended to be near genes that encode telomerase subunits and its associated complex, including genetic variants at TERT (5p15.33; rs116539972, rs7705526, rs61748181, rs71593392, and rs140648021), TERC (3q26.2; rs12638862 and rs146546514), and OBFC1 (10q24.33; rs9419958 and rs139122544). Several colocalised loci mapped to genes that have not been previously linked to lung cancer risk at genome-wide significant level: MPHOSPH6 (16q23.3; rs2303262), PRPF6 (20q13.33; rs80150989), and POLI (18q21.2; rs2276182). Other telomere maintenance genes showed limited evidence of colocalisation with lung adenocarcinoma (i.e. TERF1 and PIF1). For instance, while the RTEL1 locus (20q13.33: rs117238689, rs115610405, rs35640778, and rs35902944) harboured variants associated with both LTL and lung adenocarcinoma (Figure 3—figure supplement 1), these signals appeared to be distinct and independent of each other (avg_PP3=0.999, avg_PP4=0.001; Figure 3A, supplementary file 1h).

Figure 3 with 1 supplement see all
Download asset Open asset
Colocalisation analyses for the genetic loci defined by the 144 leukocyte telomere length (LTL) variants.

(A) Distribution of the average posterior probability for shared genetic loci between LTL and lung adenocarcinoma, highlighting in orange the telomere maintenance loci that colocalised (avg_PP4≥0.70) and in blue the ones where there was limited evidence for colocalisation (avg_PP4<0.70). Dashed red line represents the arbitrary avg_PP4 cutoff of 0.70. Representative stack plots for the multi-trait colocalisation results within (B) MPHOSPH6 and (C) OBFC1 loci, centred on a 150 kb LD window of rs2303262 and rs9419958 variants, respectively. Left Y-axis represents the –log10(p-values) of the association in the respective genome-wide association study for a given trait. The right Y-axis represents the recombination rate for the genetic loci. The X-axis represents the chromosome position. SNPs are coloured by the linkage disequilibrium correlation threshold (r2) with the query labelled SNP in European population. Sentinel SNPs within the defined LD window were labelled in each trait.

We further evaluated whether the loci colocalised between LTL and lung adenocarcinoma also shared genetic signals with other traits related to lung cancer susceptibility (supplementary file 1i). Multi-trait analyses at the 16q23.3 locus colocalised rs2303262 with MPHOSPH6 expression in lung tissue, FVC and FEV1, but not with any of the traits related to smoking behaviour (p=0.72; Figure 3B, supplementary file 1i). We additionally identified evidence of colocalisation (p=0.74) between lung adenocarcinoma, LTL, and gene expression in lung epithelial cells for two variants at the OBFC1 locus: rs139122544 and rs9419958 (Figure 3C, supplementary file 1i).

Genetically predicted LTL association with tumour features

We investigated the impact of genetically predicted LTL on lung adenocarcinoma tumour features by estimating molecular expression patterns within 343 lung adenocarcinomas tumours using principal component analysis in RNA-sequencing data. The first five components explained ~54% of the observed variance in the RNA-sequencing data (Figure 4, Figure 4—figure supplement 1). To explore the biological meaning of the five components, we performed pathway analyses for the top 500 genes with the highest loadings in each component (supplementary file 1j, supplementary file 1k). Overall, the genes correlated with each component tended to be enriched for specific cell signaling pathways (PC1: RNA processing; PC2: cell-cycle; PC3: metabolic processes; PC4: immune response; PC5: cellular response to stress and DNA damage; false discovery rate <5%; supplementary file 1l).

Figure 4 with 1 supplement see all
Download asset Open asset
Associations between molecular expression patterns of lung adenocarcinoma tumours, LTL PRS, and The Cancer Genome Atlas (TCGA) features.

(A) LTL PRS association with the first five principal components based on RNA-sequencing data of lung adenocarcinomas tumours (n=343). Results are expressed as beta estimate per SD increase in genetically predicted LTL. Linear regression model adjusted by sex, age, smoking status, and PC1-5 (genetic ancestry) covariates. Statistically significant associations with p-values<0.05 (red square). (B) Heatmap representing the correlations among PC2 and selected molecular features related to telomere length canonical roles. LTL = leukocyte telomere length; PRS = polygenic risk score; PC = principal component; TMB = tumour total mutation burden; HRD = homologous recombination deficiency, SBS (single base substitution DNA mutational signatures). SBS1 and SBS5 are DNA mutational signatures associated with age-related processes, and SBS4 is associated with tobacco smoking exposure. X-shaped marker to cross correlations with p-value>0.05.

We then tested the association between the PRS composed of the 144 genetic instruments selected for MR analysis and the five components of gene expression within lung adenocarcinoma tumours (Figure 4A). The LTL PRS was positively associated with the second component (PC2) of tumour expression (Beta = 0.17, 95% CI = 0.12–0.19, p=1.0 × 10–3; Figure 4A). In multivariate analysis, higher values of PC2 tended to be associated with patients older at diagnosis (p=0.001), female (p=0.005), being never smokers (p=0.04), and diagnosed with early-stage tumours (p=0.002; Table 1). PC2 was also highly correlated with gene expression-based measure of cell proliferation and several genomic features related to genomic stability (Figure 4B). In multivariate analysis, higher values of PC2 were associated with reduced tumour proliferation (p=3.7 × 10–30), lower somatic copy number alternations (p=1.6 × 10–05), and higher tumour telomerase activity scores (p=1.6 × 10–5). Multivariate analysis also indicated that LTL PRS remained an independent predictor of PC2 when considering these genomic features (p=0.009; Table 1). It is noteworthy only nominal associations between LTL PRS and above-mentioned features and none remained statistically significant after correction for multiple testing (supplementary file 1m). We next inferred the gene expression signature of PC2, based on 10 genes informative of this component selected by the LASSO regression models, in both lung adenocarcinoma (TCGA-LUAD) and squamous cell carcinoma (TCGA-LUSC) cohorts (Figure 5, Figure 5—figure supplement 1). The association between LTL PRS and inferred PC2 was observed in TCGA-LUAD (p=0.001) but not in TCGA-LUSC (p=0.729) cases (Figure 5A). The inferred PC2 signature levels were higher in TCGA-LUAD than in TCGA-LUSC (Figure 5B), while higher proliferation rate (Figure 5C, p=1.45 × 10–141) and TERT activity (Figure 5D, p=1.36 × 10–20) were observed in TCGA-LUSC than in TCGA-LUAD cases. Of note, the low RNA levels of the telomere-related genes (less than five read counts), such as TERT and TERC, in both TCGA-LUAD and TCGA-LUSC tumour samples limited the direct comparison of these genes between these cohorts.

Figure 5 with 1 supplement see all
Download asset Open asset
Comparing inferred PC2 gene expression signature by lung cancer histological subtypes.

(A) Leukocyte telomere length (LTL) polygenic risk score (PRS) association with the 10-gene expression signature of PC2 in lung adenocarcinoma (The Cancer Genome Atlas [TCGA]-LUAD, N=343) and squamous cell carcinoma (TCGA-LUSC, N=338) cases from TCGA dataset. Results are expressed as beta estimate per SD increase in genetically predicted LTL. Linear regression model adjusted by sex, age, and PC1-5 (genetic ancestry) covariates, PC2 signature as outcome. Statistically significant associations with -values<0.05. Values per SD of (B) PC2, (C) proliferation score, and (D) telomerase/TERT activity gene expression signatures by lung cancer histological subtypes (TCGA-LUAD and TCGA-LUSC). p-Values derived from Student’s t-tests.

Table 1
Association between PC2 (outcome) and lung adenocarcinoma tumour features in univariate and multivariate models (n=343).
Non-molecular features
PredictorsUnivariate modelMultivariate model
OR/Beta (SE)p-valueOR/Beta (SE)p-value
Age at diagnosis*0.17±0.050.0010.17±0.050.001
Gender (male)†0.73±0.110.0050.74±0.110.005
Smoking status (ever)†0.67±0.160.0130.72±0.150.035
Tumour stage (late)†0.67±0.130.0020.67±0.130.002
Molecular features
PredictorsUnivariate modelMultivariate model
Beta (SE)p-valueBeta (SE)p-value
LTL PRS ‡0.17±0.050.0010.10±0.040.009
Telomerase activity–0.37±0.059.34E-130.25±0.061.32E-05
Proliferation–0.69±0.043.30E-46–0.80±0.063.66E-30
Copy number alteration–0.41±0.056.36E-16–0.23±0.051.62E-05
Homologous recombination deficiency–0.4±0.058.32E-150.12±0.060.048
Tumour total mutation burden–0.28±0.051.37E-07–0.09±0.240.695
SBS1–0.18±0.050.0010.01±0.050.827
SBS4–0.24±0.056.36E-060.04±0.180.814
SBS5–0.27±0.054.84E-070.03±0.090.770

  1. SBS (single base substitution DNA mutational signatures); LTL = leukocyte telomere length; PC = principal component; PRS = polygenic risk score.

  2. *

    age of diagnosis represented as beta estimate per 1 unit of SD.

  3. †

    OR per 1 unit of SD.

  4. ‡

    LTL PRS is adjusted by first five PC of genetic ancestry in the univariate model.

Discussion

The maintenance of telomere length is one of the hallmarks of cancer, being critical for cell proliferation and genome integrity (Hanahan and Weinberg, 2011). Individual differences in telomere length, measured either directly or indirectly by germline determinates, have been linked with multiple diseases, including cancer susceptibility (Codd et al., 2021). The measurement of LTL within the UKBB has provided a resource for the development of a more powerful set of genetic instruments that capture a greater proportion of variation in LTL compared to previous studies (Codd et al., 2021). We applied genetic determinants of LTL to the largest GWAS of lung cancer to further characterise the role of telomere maintenance in lung cancer aetiology.

Using an MR analysis framework, we confirmed the previously reported relationship between genetically predicted longer LTL and increased risk of lung cancer. Our expanded genetic instrument detected systematic negative pleiotropy, which has not been observed in previous MR studies (Zhang et al., 2015; Haycock et al., 2017; Kachuri et al., 2019). Correcting for this pervasive directional bias resulted in substantially larger effects of LTL on risk of lung adenocarcinoma and lung cancer in never smokers, implying that LTL may even be more important to these phenotypes than previously estimated (Zhang et al., 2015; Haycock et al., 2017; Kachuri et al., 2019). Our observations in never smokers were also supported by multivariate MR analyses where adjustment for smoking did not attenuate the effect of LTL on lung cancer susceptibility. We used MVMR analysis (Sanderson et al., 2019) to assess the potential that factors, such as BMI, smoking, alcohol use, and other obesity related factors, may account for the association between LTL and lung cancer. While the influence of alternative unknown potential confounding factors cannot be excluded, the association of LTL with lung cancer risk was materially unaltered following adjustment for a range of potential confounders considered in our MVMR analyses.

Colocalisation analyses for the variants selected as the MR genetic instrument highlighted shared genetic signals between LTL and lung adenocarcinoma, including loci near genes related to telomere length maintenance (TERT, TERC, and OBFC1) and three genetic loci not previously linked to lung cancer susceptibility (POLI, PRPF6, and MPHOSPH6). The lung cancer risk allele of the MPHOSPH6 sentinel variant (rs2303262) was associated with longer LTL, reduced pulmonary function, and increased MPHOSPH6 gene expression in lung tissue. MPHOSPH6 encodes an enzyme associated with the RNA exosome complex where it modulates RNA binding activity. PRPF6 in 20q13.33 is involved in androgen binding and has been shown to promote colon tumour growth via preferential splicing of genes involved in proliferation (Adler et al., 2014). POLI is a member of the Y-family of DNA damage-tolerant polymerases involved in translesion synthesis (McIntyre, 2020). As part of its role in DNA repair and replication stress, POLI interacts with TP53 to bypass barriers during DNA replication, which may confer a pro-survival effect to stem cells and cancer cells (Guo et al., 2021).

Despite the limited evidence of colocalisation between several loci important for telomere maintenance and lung cancer, the MR analyses restricting to the non-colocalised LTL SNPs pointed out that these loci might lie in the same causal pathway of LTL and lung cancer, highlighting the heterogeneity in the genetic effects of LTL loci.

We additionally identified the relationship between genetic determinants of LTL and a specific gene expression component in lung adenocarcinoma tumours. The aspect of this component associated with longer LTL was also associated with demographic and clinical features, such as never smoking, female, and early-stage tumours compared with other lung adenocarcinoma patients, which is an interesting but not completely understood lung cancer strata. This expression component also tended to be related to genomic features related to genomically stable tumours and strikingly associated with cell proliferation score, implying that this component might be a proxy for this feature. These results appear consistent with the canonical role of telomere length in preserving genome stability and cell proliferation (de Lange, 2009).

A plausible explanation for why long LTL was associated with an increased risk of lung cancer might be that individuals with longer telomeres have lower rates of telomere attrition compared to individuals with shorter telomeres. Given a very large population of histologically normal cells, even a very small difference in telomere attrition would change the probability that a given cell is able to escape the telomere-mediated cell death pathways (Aviv et al., 2017). Such inter-individual differences could suffice to explain the modest lung cancer risk observed in our MR analyses. However, it is not clear why longer TL would be more relevant to lung adenocarcinoma compared to other lung cancer subtypes. A suggestion may come from our observation that longer LTL is related to genomically stable lung tumours (such as lung adenocarcinomas in never smokers and tumours with lower proliferation rates) but not genomically unstable lung tumours (such as heavy smoking related, highly proliferating lung squamous carcinomas). One possible hypothesis is that histologic normal cells exposed to highly genotoxic compounds, such as tobacco smoking, might require an intrinsic activation of telomere length maintenance at early stages of carcinogenesis that would allow them to survive, and therefore, genetic differences in telomere length are less relevant in these cells. By contrast, in more genomically stable lung tumours, where TL attrition rate is more modest, the hypothesis related to differences in TL length may be more relevant and potentially explain the heterogeneity in genetic effects between lung tumours. Alternately, we also note that the cell of origin may also differ, with lung adenocarcinoma postulated to be mostly derived from alveolar type 2 cells, the squamous cell carcinoma is from bronchiolar epithelium cells (Sainz de Aja et al., 2021), possibly suggesting that LTL might be more relevant to the former.

One surprising finding from the genetic analysis was that despite the robust and large effects of LTL on lung cancer risk observed in the MR analyses, the genetic correlation between LTL and lung cancer was effectively null. The LDSC approach considers genetic variants across the entire genome, whereas the MR approach preferentially selects variants based on their association with LTL, restricting to those that achieved genome-wide significance. One possibility for the lack of genetic correlation between LTL and lung cancer is that genetic variants may differ in the direction that they influence these traits, and we used the smoking behaviour traits to exemplify that. For example, the subgroup of genetic variants noted at genome-wide significance from LTL studies was associated with increased LTL and lung cancer risk. However, the subgroup related to smoking behaviours which, in turn, are linked with increased lung cancer risk, tends to decrease LTL. If such opposing effects were widespread across the genome, it could account for the lack of genetic correlation between LTL and lung cancer estimated by LDSC and highlights the complex nature of the genetic variants that determine LTL and lung cancer risk.

Some limitations of this study should be acknowledged. A limited sample size might have limited the detection of an association between lung small-cell carcinoma and LTL. Our colocalisation approach is generally more conservative and may fail to accurately determine the posterior probability for shared genetic signals in the presence of multiple independent associations in a given locus (Hukku et al., 2021), which may be a reasonable explanation for the lack of colocalisation observed at RTEL1 locus, and we stress that many of the variants that are COLOC negative are likely to be associated with lung cancer. Furthermore, the relatively small sample size of the lung adenocarcinoma cohort from TCGA may have reduced the power of our study, and larger cohorts of expression profiles tumours will be necessary to validate and explore some of our findings. The potential limitations such as collider bias within the lung adenocarcinoma case only study design should also be considered.

In conclusion, we describe an association between long genetically predicted LTL and lung cancer risk, which provides insights into how telomere length influences the genetic basis of lung cancer aetiology, including never smoking and female lung adenocarcinoma, which is an enigmatic subset of lung cancer. By using a novel framework to explore the biological implications of genetically complex traits, we unravel one gene expression component, highly correlated with proliferation rate score and other genomic stability-related features, associated with LTL in lung adenocarcinoma tumours. As far as we are aware, this is the first time an association between a PRS related to an aetiological factor, such as telomere length, and a particular expression component in the lung tumours is reported. Our findings suggest that lung adenocarcinoma patients with longer LTL might have more genomically stable tumours than the ones with shorter LTL, shedding some light on telomere biology in those tumours.

Appendix 1

Additional methods: Datasets

UK Biobank

The UK Biobank (UKBB) is a large cohort with more than 500,000 individuals from the United Kingdom who were recruited between 2006 and 2010, as previously described (Bycroft et al., 2018). At the time of recruitment, individuals provided electronic signed consent, answered questions regarding socio-demographic, lifestyle, health-related factors, and physical measures. The genetic data were obtained using the UK BiLEVE and the UKBB Axiom arrays. UKBB has ethical approval from the National Information Governance Board for Health and Social Care and the North-West Multicentre Research Ethics Committee (Ref: 11/NW/0382).

The latest genome-wide association study (GWAS) summary statistics of leukocyte telomere length (LTL; https://figshare.com/s/caa99dc0f76d62990195) were performed on 464,716 UKBB participants of European ancestry with 19.4 million imputed variants ( minor allele frequency≥ 1% and imputation quality score ≥0.3), adjusting by age, sex, array, and the first 10 principal components. LTL was measured from UKBB baseline samples, and values were log-transformed and Z-standardised.

The forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) GWAS summary statistics were conducted on 372,750 and 370,638 individuals of European genetic ancestry, respectively, who had spirometry data and sample passed quality control of genetic data (Hardy–Weinberg equilibrium at p>1 × 10−5, MAF >0.5%, and imputation quality score ≥0.3). Linear regression models were applied to standardised Z-scores for FEV1 and FVC, adjusting by age, sex, genotyping array, and 15 principal components. These GWAS summary statistics are available by the corresponding authors of the original article upon reasonable request (see Kachuri et al., 2020).

International lung cancer consortium

The International Lung Cancer Consortium (ILCCO) is an international group of lung cancer researchers for sharing comparable data from ongoing lung cancer case-control and cohort studies. The lung cancer GWAS summary statistics used in the current work is derived from a collaborative effort between Transdisciplinary Research of Cancer in Lung of the ILCCO (TRICL-ILCCO) and the Lung Cancer Cohort Consortium (LC3). All participants in these studies signed an informed consent, approved by the local internal review board or ethics committee.

The lung cancer GWAS summary statistic (29,266 cases and 56,450 controls of European descent; McKay et al., 2017) is meta-analysis of the combined imputed genotypes from OncoArray series (14,803 cases and 12,262 controls) and previous lung cancer GWAS (IARC, MDACC, SLRI, ICR, Harvard, NCI, Germany, and deCODE studies; 44,188 controls and 14,436 cases). The meta-analysis only considered imputed variants with MAF >1%, imputation quality score ≥0.3, and with little evidence for heterogeneity in effect between the studies (Cochran’s Q statistic: p>0.05), adjusting by sex, age, and principal components. Additionally, analyses were conducted by lung cancer histological subtype (adenocarcinoma [11,273 cases and 55,483 controls], squamous cell carcinoma [7426 cases and 55,627 controls], and small-cell carcinoma [2664 cases and 21,444 controls]). Of note, histological subtype information was not available for all studies.

The GWAS and Sequencing Consortium of Alcohol and Nicotine use consortium

The GWAS and Sequencing Consortium of Alcohol and Nicotine use dataset is derived from an international genetic association meta-analysis consortium for tobacco smoking and alcohol consumption traits across 30 studies and up to 1.2 million individuals (Liu et al., 2019). For the GWAS summary statistics used in the current study, UKBB participants were excluded from the meta-analysis, and only genetic variants present in more than 19 studies, MAF >1%, and imputation quality score >0.3 were considered. The smoking traits were defined as follow: cigarettes per day (CPD; N=337,334; continuous variable defined as average number of cigarettes smoked per day, either as a current smoker or former smoker), smoking initiation (SmkInit; N=1,232,091; categorical variable defined as ever smokers versus never smokers), smoking cessation (SmkCes; N=547,219; categorical variable defined as current smokers versus former smokers), and age of initiation (AgeSmk; N=341,427; continuous variable defined as the age at which an individual first became a regular smoker). Age, sex, and genetic principal components were used as covariates in all analyses. GWAS summary statistics were downloaded from https://conservancy.umn.edu/handle/11299/201564.

The Genotype-Tissue Expression (GTEx) project data is a project that aims to investigate the effects of genetic variants in transcriptome and their impact on regulatory mechanisms of a traits and diseases (GTEx Consortium, 2020). The GTEx data contain 15,201 RNA-sequencing samples from 49 tissues of 838 postmortem donors. For the current study, we accessed version 8 of the GTEx program to obtain the effect estimates and risk alleles of eQTL genetic variants (in normal lung tissue) around the queried SNPs for the multi-trait colocalisation analyses.

The Cancer Genome Atlas

The Cancer Genome Atlas (TCGA) is currently the biggest resource of cancer genomics, with more than 20,000 primary cancer and matched normal samples molecularly characterised across 33 cancer types with available demographic and clinical data.

For the current study, the access for the germline data of TCGA dataset was obtained via TCGA project #2731 via dbGAP (phs000178.v11.p8). Genotyping and imputation from TCGA data were performed in-house and described elsewhere (Gabriel et al., 2022). Briefly, germline variants with low genotyping call rate (<97%), MAF <1%, strong deviation (p-value<10–8) from the Hardy–Weinberg equilibrium, and allele frequency differing more than 20% from 1000 Genomes (1000 G) within European ancestry group were removed. Phasing was performed using Eagle (v2.4.1; Gabriel and Lipinski, 2021) with the 1000 genomes phase 3 reference panel, and Minimac4 (v1.0.1; Gabriel and Lipinski, 2021) was used to perform the imputation (window size of 500 kb; Das et al., 2016). Samples with discrepancies between genetic inferred and self-reported sex as well as with high inter-sample relatedness (PI_HAT >0.185%) were removed.

Of the 480 lung adenocarcinoma cases with genotype data, 364 cases had available tumour RNA-sequencing data (more details in Data availability section). We excluded 20 cases due to missing information for smoking status (ever smokers versus never smokers) from the downstream analyses. After summarising the normalised log-transformed read counts for all genes in two dimensions, one sample was removed due to extreme gene expression profile in comparison with all the lung adenocarcinoma tumours. Thus, 343 lung adenocarcinoma tumour samples were selected for the principal component analysis and association with LTL polygenic risk score (PRS), molecular, clinic, and demographic variables of TCGA dataset.

Data availability

Lung cancer GWAS summary statistics obtained from ILCCO can be accessed by the database of Genotypes and Phenotypes (dbGAP) under accession phs000876.v1.p1. The GWAS summary statistics for tobacco-smoking behaviors (GSCAN: https://conservancy.umn.edu/handle/11299/201564), LTL (https://figshare.com/s/caa99dc0f76d62990195), and GTEx version 8 (downloaded via GTEx google cloud resource) are publicly available. Germline data of TCGA cohorts were accessed by dbGAP under accession number phs000178.v11.p8 and project application #2731. RNA-sequencing data from TCGA cohorts were retrieved from GDC open-access data portal (https://portal.gdc.cancer.gov/) using TCGAbiolinks package in R. TCGA-related data are publicly available as described in the data section. The code for LDSC analysis is available at: https://github.com/bulik/ldsc/wiki/Heritability-andGenetic-Correlation. The codes used in this study for two-sample MR, colocalisation, multi-trait colocalisation, and principal component analyses are available at https://github.com/ricardocortezcardoso/Telomere_Length_Code, (Penha, 2023a, copy archived at swh:1:rev:f365df300919c46bb99a96b4040d90576fc878e2). Plots were created using R packages ‘meta’ (v5.5, forest plots), ‘corrplot’ (v0.92, correlation matrix), and ‘ggplot2’ (v3.3.6). The R package to generate stackplots for visualisation of the multi-trait colocalisation results is available at https://github.com/jrs95/gassocplot, (Penha, 2023b, copy archived at swh:1:rev:ae6a59dff2e43d39eead3d483af1d50f151c3d5b).

References

    1. Haycock PC
    2. Burgess S
    3. Nounu A
    4. Zheng J
    5. Okoli GN
    6. Bowden J
    7. Wade KH
    8. Timpson NJ
    9. Evans DM
    10. Willeit P
    11. Aviv A
    12. Gaunt TR
    13. Hemani G
    14. Mangino M
    15. Ellis HP
    16. Kurian KM
    17. Pooley KA
    18. Eeles RA
    19. Lee JE
    20. Fang S
    21. Chen WV
    22. Law MH
    23. Bowdler LM
    24. Iles MM
    25. Yang Q
    26. Worrall BB
    27. Markus HS
    28. Hung RJ
    29. Amos CI
    30. Spurdle AB
    31. Thompson DJ
    32. O’Mara TA
    33. Wolpin B
    34. Amundadottir L
    35. Stolzenberg-Solomon R
    36. Trichopoulou A
    37. Onland-Moret NC
    38. Lund E
    39. Duell EJ
    40. Canzian F
    41. Severi G
    42. Overvad K
    43. Gunter MJ
    44. Tumino R
    45. Svenson U
    46. van Rij A
    47. Baas AF
    48. Bown MJ
    49. Samani NJ
    50. van t’Hof FNG
    51. Tromp G
    52. Jones GT
    53. Kuivaniemi H
    54. Elmore JR
    55. Johansson M
    56. Mckay J
    57. Scelo G
    58. Carreras-Torres R
    59. Gaborieau V
    60. Brennan P
    61. Bracci PM
    62. Neale RE
    63. Olson SH
    64. Gallinger S
    65. Li D
    66. Petersen GM
    67. Risch HA
    68. Klein AP
    69. Han J
    70. Abnet CC
    71. Freedman ND
    72. Taylor PR
    73. Maris JM
    74. Aben KK
    75. Kiemeney LA
    76. Vermeulen SH
    77. Wiencke JK
    78. Walsh KM
    79. Wrensch M
    80. Rice T
    81. Turnbull C
    82. Litchfield K
    83. Paternoster L
    84. Standl M
    85. Abecasis GR
    86. SanGiovanni JP
    87. Li Y
    88. Mijatovic V
    89. Sapkota Y
    90. Low S-K
    91. Zondervan KT
    92. Montgomery GW
    93. Nyholt DR
    94. van Heel DA
    95. Hunt K
    96. Arking DE
    97. Ashar FN
    98. Sotoodehnia N
    99. Woo D
    100. Rosand J
    101. Comeau ME
    102. Brown WM
    103. Silverman EK
    104. Hokanson JE
    105. Cho MH
    106. Hui J
    107. Ferreira MA
    108. Thompson PJ
    109. Morrison AC
    110. Felix JF
    111. Smith NL
    112. Christiano AM
    113. Petukhova L
    114. Betz RC
    115. Fan X
    116. Zhang X
    117. Zhu C
    118. Langefeld CD
    119. Thompson SD
    120. Wang F
    121. Lin X
    122. Schwartz DA
    123. Fingerlin T
    124. Rotter JI
    125. Cotch MF
    126. Jensen RA
    127. Munz M
    128. Dommisch H
    129. Schaefer AS
    130. Han F
    131. Ollila HM
    132. Hillary RP
    133. Albagha O
    134. Ralston SH
    135. Zeng C
    136. Zheng W
    137. Shu X-O
    138. Reis A
    139. Uebe S
    140. Hüffmeier U
    141. Kawamura Y
    142. Otowa T
    143. Sasaki T
    144. Hibberd ML
    145. Davila S
    146. Xie G
    147. Siminovitch K
    148. Bei J-X
    149. Zeng Y-X
    150. Försti A
    151. Chen B
    152. Landi S
    153. Franke A
    154. Fischer A
    155. Ellinghaus D
    156. Flores C
    157. Noth I
    158. Ma S-F
    159. Foo JN
    160. Liu J
    161. Kim J-W
    162. Cox DG
    163. Delattre O
    164. Mirabeau O
    165. Skibola CF
    166. Tang CS
    167. Garcia-Barcelo M
    168. Chang K-P
    169. Su W-H
    170. Chang Y-S
    171. Martin NG
    172. Gordon S
    173. Wade TD
    174. Lee C
    175. Kubo M
    176. Cha P-C
    177. Nakamura Y
    178. Levy D
    179. Kimura M
    180. Hwang S-J
    181. Hunt S
    182. Spector T
    183. Soranzo N
    184. Manichaikul AW
    185. Barr RG
    186. Kahali B
    187. Speliotes E
    188. Yerges-Armstrong LM
    189. Cheng C-Y
    190. Jonas JB
    191. Wong TY
    192. Fogh I
    193. Lin K
    194. Powell JF
    195. Rice K
    196. Relton CL
    197. Martin RM
    198. Davey Smith G
    199. Telomeres Mendelian Randomization Collaboration
    (2017) Association between telomere length and risk of cancer and non-neoplastic diseases: a mendelian randomization study
    JAMA Oncology 3:636–651.
    https://doi.org/10.1001/jamaoncol.2016.5945
    1. Liu M
    2. Jiang Y
    3. Wedow R
    4. Li Y
    5. Brazel DM
    6. Chen F
    7. Datta G
    8. Davila-Velderrain J
    9. McGuire D
    10. Tian C
    11. Zhan X
    12. Choquet H
    13. Docherty AR
    14. Faul JD
    15. Foerster JR
    16. Fritsche LG
    17. Gabrielsen ME
    18. Gordon SD
    19. Haessler J
    20. Hottenga J-J
    21. Huang H
    22. Jang S-K
    23. Jansen PR
    24. Ling Y
    25. Mägi R
    26. Matoba N
    27. McMahon G
    28. Mulas A
    29. Orrù V
    30. Palviainen T
    31. Pandit A
    32. Reginsson GW
    33. Skogholt AH
    34. Smith JA
    35. Taylor AE
    36. Turman C
    37. Willemsen G
    38. Young H
    39. Young KA
    40. Zajac GJM
    41. Zhao W
    42. Zhou W
    43. Bjornsdottir G
    44. Boardman JD
    45. Boehnke M
    46. Boomsma DI
    47. Chen C
    48. Cucca F
    49. Davies GE
    50. Eaton CB
    51. Ehringer MA
    52. Esko T
    53. Fiorillo E
    54. Gillespie NA
    55. Gudbjartsson DF
    56. Haller T
    57. Harris KM
    58. Heath AC
    59. Hewitt JK
    60. Hickie IB
    61. Hokanson JE
    62. Hopfer CJ
    63. Hunter DJ
    64. Iacono WG
    65. Johnson EO
    66. Kamatani Y
    67. Kardia SLR
    68. Keller MC
    69. Kellis M
    70. Kooperberg C
    71. Kraft P
    72. Krauter KS
    73. Laakso M
    74. Lind PA
    75. Loukola A
    76. Lutz SM
    77. Madden PAF
    78. Martin NG
    79. McGue M
    80. McQueen MB
    81. Medland SE
    82. Metspalu A
    83. Mohlke KL
    84. Nielsen JB
    85. Okada Y
    86. Peters U
    87. Polderman TJC
    88. Posthuma D
    89. Reiner AP
    90. Rice JP
    91. Rimm E
    92. Rose RJ
    93. Runarsdottir V
    94. Stallings MC
    95. Stančáková A
    96. Stefansson H
    97. Thai KK
    98. Tindle HA
    99. Tyrfingsson T
    100. Wall TL
    101. Weir DR
    102. Weisner C
    103. Whitfield JB
    104. Winsvold BS
    105. Yin J
    106. Zuccolo L
    107. Bierut LJ
    108. Hveem K
    109. Lee JJ
    110. Munafò MR
    111. Saccone NL
    112. Willer CJ
    113. Cornelis MC
    114. David SP
    115. Hinds DA
    116. Jorgenson E
    117. Kaprio J
    118. Stitzel JA
    119. Stefansson K
    120. Thorgeirsson TE
    121. Abecasis G
    122. Liu DJ
    123. Vrieze S
    124. 23andMe Research Team
    125. HUNT All-In Psychiatry
    (2019) Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use
    Nature Genetics 51:237–244.
    https://doi.org/10.1038/s41588-018-0307-5
    1. McKay JD
    2. Hung RJ
    3. Han Y
    4. Zong X
    5. Carreras-Torres R
    6. Christiani DC
    7. Caporaso NE
    8. Johansson M
    9. Xiao X
    10. Li Y
    11. Byun J
    12. Dunning A
    13. Pooley KA
    14. Qian DC
    15. Ji X
    16. Liu G
    17. Timofeeva MN
    18. Bojesen SE
    19. Wu X
    20. Le Marchand L
    21. Albanes D
    22. Bickeböller H
    23. Aldrich MC
    24. Bush WS
    25. Tardon A
    26. Rennert G
    27. Teare MD
    28. Field JK
    29. Kiemeney LA
    30. Lazarus P
    31. Haugen A
    32. Lam S
    33. Schabath MB
    34. Andrew AS
    35. Shen H
    36. Hong Y-C
    37. Yuan J-M
    38. Bertazzi PA
    39. Pesatori AC
    40. Ye Y
    41. Diao N
    42. Su L
    43. Zhang R
    44. Brhane Y
    45. Leighl N
    46. Johansen JS
    47. Mellemgaard A
    48. Saliba W
    49. Haiman CA
    50. Wilkens LR
    51. Fernandez-Somoano A
    52. Fernandez-Tardon G
    53. van der Heijden HFM
    54. Kim JH
    55. Dai J
    56. Hu Z
    57. Davies MPA
    58. Marcus MW
    59. Brunnström H
    60. Manjer J
    61. Melander O
    62. Muller DC
    63. Overvad K
    64. Trichopoulou A
    65. Tumino R
    66. Doherty JA
    67. Barnett MP
    68. Chen C
    69. Goodman GE
    70. Cox A
    71. Taylor F
    72. Woll P
    73. Brüske I
    74. Wichmann H-E
    75. Manz J
    76. Muley TR
    77. Risch A
    78. Rosenberger A
    79. Grankvist K
    80. Johansson M
    81. Shepherd FA
    82. Tsao M-S
    83. Arnold SM
    84. Haura EB
    85. Bolca C
    86. Holcatova I
    87. Janout V
    88. Kontic M
    89. Lissowska J
    90. Mukeria A
    91. Ognjanovic S
    92. Orlowski TM
    93. Scelo G
    94. Swiatkowska B
    95. Zaridze D
    96. Bakke P
    97. Skaug V
    98. Zienolddiny S
    99. Duell EJ
    100. Butler LM
    101. Koh W-P
    102. Gao Y-T
    103. Houlston RS
    104. McLaughlin J
    105. Stevens VL
    106. Joubert P
    107. Lamontagne M
    108. Nickle DC
    109. Obeidat M
    110. Timens W
    111. Zhu B
    112. Song L
    113. Kachuri L
    114. Artigas MS
    115. Tobin MD
    116. Wain LV
    117. Rafnar T
    118. Thorgeirsson TE
    119. Reginsson GW
    120. Stefansson K
    121. Hancock DB
    122. Bierut LJ
    123. Spitz MR
    124. Gaddis NC
    125. Lutz SM
    126. Gu F
    127. Johnson EO
    128. Kamal A
    129. Pikielny C
    130. Zhu D
    131. Lindströem S
    132. Jiang X
    133. Tyndale RF
    134. Chenevix-Trench G
    135. Beesley J
    136. Bossé Y
    137. Chanock S
    138. Brennan P
    139. Landi MT
    140. Amos CI
    141. SpiroMeta Consortium
    (2017) Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes
    Nature Genetics 49:1126–1132.
    https://doi.org/10.1038/ng.3892
    1. Saunders GRB
    2. Wang X
    3. Chen F
    4. Jang S-K
    5. Liu M
    6. Wang C
    7. Gao S
    8. Jiang Y
    9. Khunsriraksakul C
    10. Otto JM
    11. Addison C
    12. Akiyama M
    13. Albert CM
    14. Aliev F
    15. Alonso A
    16. Arnett DK
    17. Ashley-Koch AE
    18. Ashrani AA
    19. Barnes KC
    20. Barr RG
    21. Bartz TM
    22. Becker DM
    23. Bielak LF
    24. Benjamin EJ
    25. Bis JC
    26. Bjornsdottir G
    27. Blangero J
    28. Bleecker ER
    29. Boardman JD
    30. Boerwinkle E
    31. Boomsma DI
    32. Boorgula MP
    33. Bowden DW
    34. Brody JA
    35. Cade BE
    36. Chasman DI
    37. Chavan S
    38. Chen Y-DI
    39. Chen Z
    40. Cheng I
    41. Cho MH
    42. Choquet H
    43. Cole JW
    44. Cornelis MC
    45. Cucca F
    46. Curran JE
    47. de Andrade M
    48. Dick DM
    49. Docherty AR
    50. Duggirala R
    51. Eaton CB
    52. Ehringer MA
    53. Esko T
    54. Faul JD
    55. Fernandes Silva L
    56. Fiorillo E
    57. Fornage M
    58. Freedman BI
    59. Gabrielsen ME
    60. Garrett ME
    61. Gharib SA
    62. Gieger C
    63. Gillespie N
    64. Glahn DC
    65. Gordon SD
    66. Gu CC
    67. Gu D
    68. Gudbjartsson DF
    69. Guo X
    70. Haessler J
    71. Hall ME
    72. Haller T
    73. Harris KM
    74. He J
    75. Herd P
    76. Hewitt JK
    77. Hickie I
    78. Hidalgo B
    79. Hokanson JE
    80. Hopfer C
    81. Hottenga J
    82. Hou L
    83. Huang H
    84. Hung Y-J
    85. Hunter DJ
    86. Hveem K
    87. Hwang S-J
    88. Hwu C-M
    89. Iacono W
    90. Irvin MR
    91. Jee YH
    92. Johnson EO
    93. Joo YY
    94. Jorgenson E
    95. Justice AE
    96. Kamatani Y
    97. Kaplan RC
    98. Kaprio J
    99. Kardia SLR
    100. Keller MC
    101. Kelly TN
    102. Kooperberg C
    103. Korhonen T
    104. Kraft P
    105. Krauter K
    106. Kuusisto J
    107. Laakso M
    108. Lasky-Su J
    109. Lee W-J
    110. Lee JJ
    111. Levy D
    112. Li L
    113. Li K
    114. Li Y
    115. Lin K
    116. Lind PA
    117. Liu C
    118. Lloyd-Jones DM
    119. Lutz SM
    120. Ma J
    121. Mägi R
    122. Manichaikul A
    123. Martin NG
    124. Mathur R
    125. Matoba N
    126. McArdle PF
    127. McGue M
    128. McQueen MB
    129. Medland SE
    130. Metspalu A
    131. Meyers DA
    132. Millwood IY
    133. Mitchell BD
    134. Mohlke KL
    135. Moll M
    136. Montasser ME
    137. Morrison AC
    138. Mulas A
    139. Nielsen JB
    140. North KE
    141. Oelsner EC
    142. Okada Y
    143. Orrù V
    144. Palmer ND
    145. Palviainen T
    146. Pandit A
    147. Park SL
    148. Peters U
    149. Peters A
    150. Peyser PA
    151. Polderman TJC
    152. Rafaels N
    153. Redline S
    154. Reed RM
    155. Reiner AP
    156. Rice JP
    157. Rich SS
    158. Richmond NE
    159. Roan C
    160. Rotter JI
    161. Rueschman MN
    162. Runarsdottir V
    163. Saccone NL
    164. Schwartz DA
    165. Shadyab AH
    166. Shi J
    167. Shringarpure SS
    168. Sicinski K
    169. Skogholt AH
    170. Smith JA
    171. Smith NL
    172. Sotoodehnia N
    173. Stallings MC
    174. Stefansson H
    175. Stefansson K
    176. Stitzel JA
    177. Sun X
    178. Syed M
    179. Tal-Singer R
    180. Taylor AE
    181. Taylor KD
    182. Telen MJ
    183. Thai KK
    184. Tiwari H
    185. Turman C
    186. Tyrfingsson T
    187. Wall TL
    188. Walters RG
    189. Weir DR
    190. Weiss ST
    191. White WB
    192. Whitfield JB
    193. Wiggins KL
    194. Willemsen G
    195. Willer CJ
    196. Winsvold BS
    197. Xu H
    198. Yanek LR
    199. Yin J
    200. Young KL
    201. Young KA
    202. Yu B
    203. Zhao W
    204. Zhou W
    205. Zöllner S
    206. Zuccolo L
    207. Batini C
    208. Bergen AW
    209. Bierut LJ
    210. David SP
    211. Gagliano Taliun SA
    212. Hancock DB
    213. Jiang B
    214. Munafò MR
    215. Thorgeirsson TE
    216. Liu DJ
    217. Vrieze S
    218. 23andMe Research Team
    219. Biobank Japan Project
    (2022) Genetic diversity fuels gene discovery for tobacco and alcohol use
    Nature 612:720–724.
    https://doi.org/10.1038/s41586-022-05477-4

Article and author information

Author details

  1. Ricardo Cortez Cardoso Penha

    Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), Lyon, France
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    For correspondence
    cortezr@iarc.who.int
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2847-1993

  2. Karl Smith-Byrne

    Cancer Epidemiology Unit, University of Oxford, Oxford, United Kingdom
    Contribution
    Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared

  3. Joshua R Atkins

    Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), Lyon, France
    Contribution
    Formal analysis, Methodology, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared

  4. Philip C Haycock

    MRC Integrative Epidemiology Unit, Bristol Population Health Science Institute, Bristol Medical School (PHS), Bristol, United Kingdom
    Contribution
    Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5001-3350

  5. Siddhartha Kar

    MRC Integrative Epidemiology Unit, Bristol Population Health Science Institute, Bristol Medical School (PHS), Bristol, United Kingdom
    Contribution
    Writing – review and editing
    Competing interests
    No competing interests declared

  6. Veryan Codd

    1. Department of Cardiovascular Sciences, University of Leicester, Leicester, United Kingdom
    2. NIHR Leicester Biomedical Research Centre, Glenfield Hospital, Leicester, United Kingdom
    Contribution
    Data curation, Writing – review and editing
    Competing interests
    No competing interests declared

  7. Nilesh J Samani

    1. Department of Cardiovascular Sciences, University of Leicester, Leicester, United Kingdom
    2. NIHR Leicester Biomedical Research Centre, Glenfield Hospital, Leicester, United Kingdom
    Contribution
    Data curation, Writing – review and editing
    Competing interests
    No competing interests declared

  8. Christopher Nelson

    1. Department of Cardiovascular Sciences, University of Leicester, Leicester, United Kingdom
    2. NIHR Leicester Biomedical Research Centre, Glenfield Hospital, Leicester, United Kingdom
    Contribution
    Data curation, Writing – review and editing
    Competing interests
    No competing interests declared

  9. Maja Milojevic

    Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), Lyon, France
    Contribution
    Data curation, Writing – review and editing
    Competing interests
    No competing interests declared

  10. Aurélie AG Gabriel

    Ludwig Lausanne Branch, Faculty of Biology and Medicine, Lausanne, Switzerland
    Contribution
    Data curation, Writing – review and editing
    Competing interests
    No competing interests declared

  11. Christopher Amos

    Institute for Clinical and Translational Research, Baylor College of Medicine, Houston, United States
    Contribution
    Funding acquisition, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared

  12. Paul Brennan

    Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), Lyon, France
    Contribution
    Writing – review and editing
    Competing interests
    No competing interests declared

  13. Rayjean J Hung

    Lunenfeld-Tanenbaum Research Institute, Sinai Health, Toronto, Canada
    Contribution
    Writing – review and editing
    Competing interests
    No competing interests declared

  14. Linda Kachuri

    Departament of Epidemiology and Population Health, Stanford University, Stanford, United States
    Contribution
    Formal analysis, Investigation, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared

  15. James D Mckay

    Genomic Epidemiology branch, International Agency for Research on Cancer/World Health Organization (IARC/WHO), Lyon, France
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    mckayj@iarc.who.int
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1787-3874

Funding

Terry Fox Foundation (IARC Postdoctoral Fellowship)

  • Ricardo Cortez Cardoso Penha

National Institutes of Health (K99CA246076)

  • Linda Kachuri

Cancer Prevention and Research Institute of Texas (RR170048)

  • Christopher Amos

Cancer Research UK (K99CA246076)

  • Philip C Haycock

Cancer Research UK (C18281/A29019)

  • Philip C Haycock

Institut National Du Cancer (GeniLuc 2017–1-TABAC-03-CIRC-1 - [TABAC 17‐022])

  • Philip C Haycock

National Cancer Institute (INTEGRAL NIH 5U19CA203654-03)

  • Philip C Haycock

Agence Nationale pour la Recherche (ANR-10-INBS-09)

  • Philip C Haycock

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We would like to acknowledge the TCGA Research Network (https://www.cancer.gov/tcga) and the contribution of specimen donors and research groups involved in this resource. We also would like to acknowledge the ILCCO consortium, the participants of the UK biobank and GTEx project and the supporting bodies (https://commonfund.nih.gov/GTEx), specimen donors, and research groups. This work was supported by the Institut National du Cancer (INCa) (GeniLuc 2017–1-TABAC-03-CIRC-1 - [TABAC 17‐022]), NIH/NCI, INTEGRAL NIH 5U19CA203654-03, Cancer Research UK (grant number C18281/A29019), the France Génomique National infrastructure, funded as part of the « Investissements d’Avenir » program managed by the Agence Nationale pour la Recherche (contract ANR-10-INBS-09). Christopher Amos is a Research Scholar of the Cancer Prevention Institute of Texas ( RR170048). The work of Ricardo Cortez Cardoso Penha reported in this paper was undertaken during the tenure of an IARC Postdoctoral Fellowship at the International Agency for Research on Cancer. Linda Kachuri is supported by funding from the National Institutes of Health (K99CA246076). Philip Haycock is supported by funding from Cancer Research UK (C18281/A29019).

Ethics

Human subjects: Mendelian Randomization analyses did not require ethical approval as it used secondary, genome-wide association data from studies that obtained informed consent from all participants and ethical approval from review boards and/or ethics committees. Individual level data, accessed under accession number phs000178.v11.p8 and project application #2731, was subject to the ethnical policies that govern the Cancer Genome Atlas. These policies can be found at https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga/history/policies.

Copyright

© 2023, Cortez Cardoso Penha et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,609
    views
  • 212
    downloads
  • 11
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Citations by DOI

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ricardo Cortez Cardoso Penha
  2. Karl Smith-Byrne
  3. Joshua R Atkins
  4. Philip C Haycock
  5. Siddhartha Kar
  6. Veryan Codd
  7. Nilesh J Samani
  8. Christopher Nelson
  9. Maja Milojevic
  10. Aurélie AG Gabriel
  11. Christopher Amos
  12. Paul Brennan
  13. Rayjean J Hung
  14. Linda Kachuri
  15. James D Mckay
(2023)
Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome
eLife 12:e83118.
https://doi.org/10.7554/eLife.83118

Share this article

https://doi.org/10.7554/eLife.83118

Categories and tags

Research organism