Fine mapping of chromosome 5p15.33 based on a targeted deep sequencing and high density genotyping identifies novel lung cancer susceptibility loci

Candidate pathway analysis of surfactant proteins identifies CTSH and SFTA2 that influences lung cancer risk

Pulmonary surfactant is a lipoprotein synthesized and secreted by alveolar type II cells in lung. We evaluated the associations between 200,139 single nucleotide polymorphisms (SNPs) of 40 surfactant-related genes and lung cancer risk using genotyped data from two independent lung cancer genome-wide association studies. Discovery data included 18,082 cases and 13,780 controls of European ancestry. Replication data included 1,914 cases and 3,065 controls of European descent. Using multivariate logistic regression, we found novel SNPs in surfactant-related genes CTSH [rs34577742 C > T, odds ratio (OR) = 0.90, 95% confidence interval (CI) = 0.89-0.93, P = 7.64 × 10-9] and SFTA2 (rs3095153 G > A, OR = 1.16, 95% CI = 1.10-1.21, P = 1.27 × 10-9) associated with overall lung cancer in the discovery data and validated in an independent replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.80-0.96, P = 5.76 × 10-3) and SFTA2 (rs3095153 G > A, OR = 1.14, 95% CI = 1.01-1.28, P = 3.25 × 10-2). Among ever smokers, we found SNPs in CTSH (rs34577742 C > T, OR = 0.89, 95% CI = 0.85-0.92, P = 1.94 × 10-7) and SFTA2 (rs3095152 G > A, OR = 1.20, 95% CI = 1.14-1.27, P = 4.25 × 10-11) associated with overall lung cancer in the discovery data and validated in the replication data-CTSH (rs34577742 C > T, OR = 0.88, 95% CI = 0.79-0.97, P = 1.64 × 10-2) and SFTA2 (rs3095152 G > A, OR = 1.15, 95% CI = 1.01-1.30, P = 3.81 × 10-2). Subsequent transcriptome-wide association study using expression weights from a lung expression quantitative trait loci study revealed genes most strongly associated with lung cancer are CTSH (PTWAS = 2.44 × 10-4) and SFTA2 (PTWAS = 2.32 × 10-6).

Ethnicity-specific association between TERT rs2736100 (A > C) polymorphism and lung cancer risk: a comprehensive meta-analysis

The rs2736100 (A > C) polymorphism of the second intron of Telomerase reverse transcriptase (TERT) has been confirmed to be closely associated with the risk of Lung cancer (LC), but there is still no unified conclusion on the results of its association with LC. This study included Genome-wide association studies (GWAS) and case-control studies reported so far on this association between TERT rs2736100 polymorphism and LC to clarify such a correlation with LC and the differences in it between different ethnicities and different types of LC. Relevant literatures published before May 7, 2022 on 'TERT rs2736100 polymorphism and LC susceptibility' in PubMed, EMbase, CENTRAL, MEDLINE databases were searched through the Internet, and data were extracted. Statistical analysis of data was performed in Revman5.3 software, including drawing forest diagrams, drawing funnel diagrams and so on. Sensitivity and publication bias analysis were performed in Stata 12.0 software. The C allele of TERT rs2736100 was associated with the risk of LC (Overall population: [OR] = 1.21, 95%CI [1.17, 1.25]; Caucasians: [OR] = 1.11, 95%CI [1.06, 1.17]; Asians: [OR] = 1.26, 95%CI [1.21, 1.30]), and Asians had a higher risk of LC than Caucasians (C vs. A: Caucasians: [OR] = 1.11 /Asians: [OR]) = 1.26). The other gene models also showed similar results. The results of stratified analysis of LC patients showed that the C allele was associated with the risk of Non-small-cell lung carcinoma (NSCLC) and Lung adenocarcinoma (LUAD), and the risk of NSCLC and LUAD in Asians was higher than that in Caucasians. The C allele was associated with the risk of Lung squamous cell carcinoma (LUSC) and Small cell lung carcinoma(SCLC) in Asians but not in Caucasians. NSCLC patients ([OR] = 1.27) had a stronger correlation than SCLC patients ([OR] = 1.03), and LUAD patients ([OR] = 1.32) had a stronger correlation than LUSC patients ([OR] = 1.09).In addition, the C allele of TERT rs2736100 was associated with the risk of LC, NSCLC and LUAD in both smoking groups and non-smoking groups, and the risk of LC in non-smokers of different ethnic groups was higher than that in smokers. In the Asians, non-smoking women were more at risk of developing LUAD. The C allele of TERT rs2736100 is a risk factor for LC, NSCLC, and LUAD in different ethnic groups, and the Asian population is at a more common risk. The C allele is a risk factor for LUSC and SCLC in Asians but not in Caucasians. And smoking is not the most critical factor that causes variation in TERT rs2736100 to increase the risk of most LC (NSCLC, LUAD). Therefore, LC is a multi-etiological disease caused by a combination of genetic, environmental and lifestyle factors.

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

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).

Fanconi anemia and dyskeratosis congenita/telomere biology disorders: Two inherited bone marrow failure syndromes with genomic instability

Inherited bone marrow failure syndromes (IBMFS) are a complex and heterogeneous group of genetic diseases. To date, at least 13 IBMFS have been characterized. Their pathophysiology is associated with germline pathogenic variants in genes that affect hematopoiesis. A couple of these diseases also have genomic instability, Fanconi anemia due to DNA damage repair deficiency and dyskeratosis congenita/telomere biology disorders as a result of an alteration in telomere maintenance. Patients can have extramedullary manifestations, including cancer and functional or structural physical abnormalities. Furthermore, the phenotypic spectrum varies from cryptic features to patients with significantly evident manifestations. These diseases require a high index of suspicion and should be considered in any patient with abnormal hematopoiesis, even if extramedullary manifestations are not evident. This review describes the disrupted cellular processes that lead to the affected maintenance of the genome structure, contrasting the dysmorphological and oncological phenotypes of Fanconi anemia and dyskeratosis congenita/telomere biology disorders. Through a dysmorphological analysis, we describe the phenotypic features that allow to make the differential diagnosis and the early identification of patients, even before the onset of hematological or oncological manifestations. From the oncological perspective, we analyzed the spectrum and risks of cancers in patients and carriers.

Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10-8); rs79942605 (OR = 2.17, P = 2.81×10-8) and rs208908 (OR = 0.70, P = 4.54×10-8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

Functional studies of lung cancer GWAS beyond association

Fourteen years after the first genome-wide association study (GWAS) of lung cancer was published, approximately forty-five genomic loci have now been significantly associated with lung cancer risk. While functional characterization was performed for several of these loci, a comprehensive summary of current molecular understanding of lung cancer risk has been lacking. Further, many novel computational and experimental tools now became available to accelerate the functional assessment of disease-associated variants, moving beyond locus-by-locus approaches. In this review, we first highlight the heterogeneity of lung cancer GWAS findings across histological subtypes, ancestries, and smoking status, which poses unique challenges to follow-up studies. We then summarize the published lung cancer post-GWAS studies for each risk-associated locus to assess the current understanding of biological mechanisms beyond the initial statistical association. We further summarize strategies for GWAS functional follow-up studies considering cutting-edge functional genomics tools and providing a catalog of available resources relevant to lung cancer. Overall, we aim to highlight the importance of integrating computational and experimental approaches to draw biological insights from the lung cancer GWAS results beyond association.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Assessing Lung Cancer Absolute Risk Trajectory Based on a Polygenic Risk Model

Lung cancer is the leading cause of cancer-related death globally. An improved risk stratification strategy can increase efficiency of low-dose CT (LDCT) screening. Here we assessed whether individual's genetic background has clinical utility for risk stratification in the context of LDCT screening. On the basis of 13,119 patients with lung cancer and 10,008 controls with European ancestry in the International Lung Cancer Consortium, we constructed a polygenic risk score (PRS) via 10-fold cross-validation with regularized penalized regression. The performance of risk model integrating PRS, including calibration and ability to discriminate, was assessed using UK Biobank data (N = 335,931). Absolute risk was estimated on the basis of age-specific lung cancer incidence and all-cause mortality as competing risk. To evaluate its potential clinical utility, the PRS distribution was simulated in the National Lung Screening Trial (N = 50,772 participants). The lung cancer ORs for individuals at the top decile of the PRS distribution versus those at bottom 10% was 2.39 [95% confidence interval (CI) = 1.92-3.00; P = 1.80 × 10-14] in the validation set (P trend = 5.26 × 10-20). The OR per SD of PRS increase was 1.26 (95% CI = 1.20-1.32; P = 9.69 × 10-23) for overall lung cancer risk in the validation set. When considering absolute risks, individuals at different PRS deciles showed differential trajectories of 5-year and cumulative absolute risk. The age reaching the LDCT screening recommendation threshold can vary by 4 to 8 years, depending on the individual's genetic background, smoking status, and family history. Collectively, these results suggest that individual's genetic background may inform the optimal lung cancer LDCT screening strategy. SIGNIFICANCE: Three large-scale datasets reveal that, after accounting for risk factors, an individual's genetics can affect their lung cancer risk trajectory, thus may inform the optimal timing for LDCT screening.

Mild catalytic defects of tert rs61748181 polymorphism affect the clinical presentation of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a disorder of accelerated lung aging. Multiple pieces of evidence support that the aging biomarker short telomeres, which can be caused by mutations in telomerase reverse transcriptase (TERT), contribute to COPD pathogenesis. We hypothesized that short telomere risk-associated single nucleotide polymorphisms (SNPs) in TERT, while not able to drive COPD development, nonetheless modify the disease presentation. We set out to test the SNP carrying status in a longitudinal study of smokers with COPD and found that rapid decline of FEV1 in lung function was associated with the minor allele of rs61748181 (adjusted odds ratio 2.49, p = 0.038). Biochemical evaluation of ex vivo engineered human cell models revealed that primary cells expressing the minor allele of rs61748181 had suboptimal telomere length maintenance due to reduced telomerase catalytic activity, despite having comparable cell growth kinetics as WT-TERT expressing cells. This ex vivo observation translated clinically in that shorter telomeres were found in minor allele carriers in a sub-population of COPD patients with non-declining lung function, over the 5-year period of the longitudinal study. Collectively, our data suggest that functional TERT SNPs with mild catalytic defects are nonetheless implicated in the clinical presentation of COPD.

Rare deleterious germline variants and risk of lung cancer

Recent studies suggest that rare variants exhibit stronger effect sizes and might play a crucial role in the etiology of lung cancers (LC). Whole exome plus targeted sequencing of germline DNA was performed on 1045 LC cases and 885 controls in the discovery set. To unveil the inherited causal variants, we focused on rare and predicted deleterious variants and small indels enriched in cases or controls. Promising candidates were further validated in a series of 26,803 LCs and 555,107 controls. During discovery, we identified 25 rare deleterious variants associated with LC susceptibility, including 13 reported in ClinVar. Of the five validated candidates, we discovered two pathogenic variants in known LC susceptibility loci, ATM p.V2716A (Odds Ratio [OR] 19.55, 95%CI 5.04-75.6) and MPZL2 p.I24M frameshift deletion (OR 3.88, 95%CI 1.71-8.8); and three in novel LC susceptibility genes, POMC c.*28delT at 3' UTR (OR 4.33, 95%CI 2.03-9.24), STAU2 p.N364M frameshift deletion (OR 4.48, 95%CI 1.73-11.55), and MLNR p.Q334V frameshift deletion (OR 2.69, 95%CI 1.33-5.43). The potential cancer-promoting role of selected candidate genes and variants was further supported by endogenous DNA damage assays. Our analyses led to the identification of new rare deleterious variants with LC susceptibility. However, in-depth mechanistic studies are still needed to evaluate the pathogenic effects of these specific alleles.

Genomic analyses of flow-sorted Hodgkin Reed-Sternberg cells reveal complementary mechanisms of immune evasion

Classical Hodgkin lymphoma (cHL) is composed of rare malignant Hodgkin Reed-Sternberg (HRS) cells within an extensive, but ineffective, inflammatory/immune cell infiltrate. HRS cells exhibit near-universal somatic copy gains of chromosome 9p/9p24.1, which increase expression of the programmed cell death protein 1 (PD-1) ligands. To define genetic mechanisms of response and resistance to PD-1 blockade and identify complementary treatment targets, we performed whole-exome sequencing of flow cytometry-sorted HRS cells from 23 excisional biopsies of newly diagnosed cHLs, including 8 Epstein-Barr virus-positive (EBV+) tumors. We identified significantly mutated cancer candidate genes (CCGs) as well as somatic copy number alterations and structural variations and characterized their contribution to disease-defining immune evasion mechanisms and nuclear factor κB (NF-κB), JAK/STAT, and PI3K signaling pathways. EBV- cHLs had a higher prevalence of genetic alterations in the NF-κB and major histocompatibility complex class I antigen presentation pathways. In this young cHL cohort (median age, 26 years), we identified a predominant mutational signature of spontaneous deamination of cytosine- phosphate-guanines ("Aging"), in addition to apolipoprotein B mRNA editing catalytic polypeptide-like, activation-induced cytidine deaminase, and microsatellite instability (MSI)-associated hypermutation. In particular, the mutational burden in EBV- cHLs was among the highest reported, similar to that of carcinogen-induced tumors. Together, the overall high mutational burden, MSI-associated hypermutation, and newly identified genetic alterations represent additional potential bases for the efficacy of PD-1 blockade in cHL. Of note, recurrent cHL alterations, including B2M, TNFAIP3, STAT6, GNA13, and XPO1 mutations and 2p/2p15, 6p21.32, 6q23.3, and 9p/9p24.1 copy number alterations, were also identified in >20% of primary mediastinal B-cell lymphomas, highlighting shared pathogenetic mechanisms in these diseases.

Investigation of Leukocyte Telomere Length and Genetic Variants in Chromosome 5p15.33 as Prognostic Markers in Lung Cancer

BACKGROUND

Lung cancer remains the leading cause of cancer mortality with relatively few prognostic biomarkers. We investigated associations with overall survival for telomere length (TL) and genetic variation in chromosome 5p15.33, an established telomere maintenance locus.

METHODS

Leukocyte TL was measured after diagnosis in 807 patients with non-small cell lung cancer (NSCLC) from the Princess Margaret Cancer Center in Toronto and assessed prospectively in 767 NSCLC cases from the Copenhagen City Heart Study and the Copenhagen General Population Study. Associations with all-cause mortality were tested for 723 variants in 5p15.33, genotyped in 4,672 NSCLC cases.

RESULTS

Short telomeres (≤10th percentile) were associated with poor prognosis for adenocarcinoma in both populations: TL measured 6 months after diagnosis [HR = 1.65; 95% confidence intervals (CI), 1.04-2.64] and for those diagnosed within 5 years after blood sampling (HR = 2.42; 95% CI, 1.37-4.28). Short TL was associated with mortality in never smokers with NSCLC (HR = 10.29; 95% CI, 1.86-56.86) and adenocarcinoma (HR = 11.31; 95% CI, 1.96-65.24). Analyses in 5p15.33 identified statistically significant prognostic associations for rs56266421-G in LPCAT1 (HR = 1.86; 95% CI, 1.38-2.52; P = 4.5 × 10^-5) in stage I-IIIA NSCLC, and for the SLC6A3 gene with OS in females with NSCLC (P = 1.6 × 10^-3).

CONCLUSIONS

Our findings support the potential clinical utility of TL, particularly for adenocarcinoma patients, while associations in chromosome 5p15.33 warrant further exploration.

IMPACT

This is the largest lung cancer study of leukocyte TL and OS, and the first to examine the impact of the timing of TL measurement. Our findings suggest that extremely short telomeres are indicative of poor prognosis in NSCLC.

Protein-altering germline mutations implicate novel genes related to lung cancer development

Few germline mutations are known to affect lung cancer risk. We performed analyses of rare variants from 39,146 individuals of European ancestry and investigated gene expression levels in 7,773 samples. We find a large-effect association with an ATM L2307F (rs56009889) mutation in adenocarcinoma for discovery (adjusted Odds Ratio = 8.82, P = 1.18 × 10-15) and replication (adjusted OR = 2.93, P = 2.22 × 10-3) that is more pronounced in females (adjusted OR = 6.81 and 3.19 and for discovery and replication). We observe an excess loss of heterozygosity in lung tumors among ATM L2307F allele carriers. L2307F is more frequent (4%) among Ashkenazi Jewish populations. We also observe an association in discovery (adjusted OR = 2.61, P = 7.98 × 10-22) and replication datasets (adjusted OR = 1.55, P = 0.06) with a loss-of-function mutation, Q4X (rs150665432) of an uncharacterized gene, KIAA0930. Our findings implicate germline genetic variants in ATM with lung cancer susceptibility and suggest KIAA0930 as a novel candidate gene for lung cancer risk.

TERT-rs33963617 and CLPTM1L-rs77518573 reduce the risk of non-small cell lung cancer in Chinese population

BACKGROUND

Genome-wide association studies (GWAS) have identified 5p15.33 as a susceptible locus for lung cancer. However, for non-small cell lung cancer (NSCLC), low-frequency risk variants in this region have not been systematically studied. We intended to explore the associations between low-frequency variants on 5p15.33 and NSCLC using a next-generation sequencing based approach in this study.

METHODS

We have acquisited the variation spectrum of 400 NSCLC patients on 5p15.33 by sequencing the targeted region before. Candidate variants were primarily selected by restricting the minor allele frequency (MAF 1-5%) and then by comparing their frequency in 400 NSCLC patients with 1008 East Asians from The genome Aggregation Database (gnomAD). The associations between candidate variants and NSCLC were discovered and replicated in two case-control sets: discovery stage with 960 cases and 916 controls, and replication stage with 1596 cases and 1614 controls in total.

RESULTS

Five low-frequency variants were selected as our candidates and subsequent association analyses showed that 2 polymorphisms were significantly associated with risk of NSCLC, including rs33963617 (OR = 0.63, 95% CI: 0.53-0.76, P = 3.80 × 10^-7) in TERT and rs77518573 (OR = 0.73, 95% CI: 0.63-0.84, P = 2.00 × 10^-5) in upstream of CLPTM1L. When stratified by histologic subtype, a significant association was only investigated in adenocarcinoma for rs77518573. We also observed an obvious cumulative effect of the two significant variants.

CONCLUSIONS

We newly identified two NSCLC related variants on chromosome 5p15.33. Both TERT-rs33963617 and CLPTM1L-rs77518573 conferred reduced risk for NSCLC in Chinese Han population.

Lung Cancer Risk in Never-Smokers of European Descent is Associated With Genetic Variation in the 5p15.33 TERT-CLPTM1Ll Region

INTRODUCTION

Inherited susceptibility to lung cancer risk in never-smokers is poorly understood. The major reason for this gap in knowledge is that this disease is relatively uncommon (except in Asians), making it difficult to assemble an adequate study sample. In this study we conducted a genome-wide association study on the largest, to date, set of European-descent never-smokers with lung cancer.

METHODS

We conducted a two-phase (discovery and replication) genome-wide association study in never-smokers of European descent. We further augmented the sample by performing a meta-analysis with never-smokers from the recent OncoArray study, which resulted in a total of 3636 cases and 6295 controls. We also compare our findings with those in smokers with lung cancer.

RESULTS

We detected three genome-wide statistically significant single nucleotide polymorphisms rs31490 (odds ratio [OR]: 0.769, 95% confidence interval [CI]: 0.722-0.820; p value 5.31 × 10-16), rs380286 (OR: 0.770, 95% CI: 0.723-0.820; p value 4.32 × 10-16), and rs4975616 (OR: 0.778, 95% CI: 0.730-0.829; p value 1.04 × 10-14). All three mapped to Chromosome 5 CLPTM1L-TERT region, previously shown to be associated with lung cancer risk in smokers and in never-smoker Asian women, and risk of other cancers including breast, ovarian, colorectal, and prostate.

CONCLUSIONS

We found that genetic susceptibility to lung cancer in never-smokers is associated to genetic variants with pan-cancer risk effects. The comparison with smokers shows that top variants previously shown to be associated with lung cancer risk only confer risk in the presence of tobacco exposure, underscoring the importance of gene-environment interactions in the etiology of this disease.

Fine Mapping in Chromosome 3q28 Identified Two Variants Associated with Lung Cancer Risk in Asian Population

Genome-wide association studies (GWASs) have consistently identified chromosome 3q28 as a lung cancer susceptibility region. To further characterize the potential genetic mechanism of the variants in this region, we conducted a fine-mapping study on chromosome 3q28 region. We performed a target resequencing in 200 lung cancer cases and 300 controls in the screening and followed by validation in multi-ethnic lung cancer GWASs with 12,843 cases and 12,639 controls. For our identified novel variants, we conducted expression quantitative trait loci (eQTL) analysis to reveal the potential target genes. Two susceptibility variants were identified (rs4396880: G>A, OR = 0.35, 95%CI: 0.20-0.62, P = 3.01×10^-4; and rs3856776: C>T, OR = 2.05, 95%CI: 1.32-3.18, P = 1.49×10^-3) and further supported in Asian population (rs4396880: OR = 0.88, P = 7.43×10^-6; and rs3856776: OR =1.17, P = 1.64×10^-4). The eQTL analysis showed the A allele of rs4396880 was significantly associated with higher mRNA expression of TP63 (P = 1.70×10^-4) in lung tissues, while rs3856776 showed significant association with the expression of LEPREL1-AS1 (P = 6.90×10^-3), which was the antisense RNA of LEPREL1 and could suppress the translation of LEPREL1. Notably, LEPREL1 was aberrantly downregulated (P = 2.54×10^-18) in lung tumor tissues based on TCGA database. In conclusion, this is the first fine-mapping analysis of 3q28 region in Han Chinese, and we found two variants associated with lung cancer susceptibility in Asian population. What's more, rs3856776 was newly identified and might modulate lung cancer susceptibility by suppressing the function of LEPREL1.

Mendelian Randomization and mediation analysis of leukocyte telomere length and risk of lung and head and neck cancers

BACKGROUND

Evidence from observational studies of telomere length (TL) has been conflicting regarding its direction of association with cancer risk. We investigated the causal relevance of TL for lung and head and neck cancers using Mendelian Randomization (MR) and mediation analyses.

METHODS

We developed a novel genetic instrument for TL in chromosome 5p15.33, using variants identified through deep-sequencing, that were genotyped in 2051 cancer-free subjects. Next, we conducted an MR analysis of lung (16 396 cases, 13 013 controls) and head and neck cancer (4415 cases, 5013 controls) using eight genetic instruments for TL. Lastly, the 5p15.33 instrument and distinct 5p15.33 lung cancer risk loci were evaluated using two-sample mediation analysis, to quantify their direct and indirect, telomere-mediated, effects.

RESULTS

The multi-allelic 5p15.33 instrument explained 1.49-2.00% of TL variation in our data (p = 2.6 × 10-9). The MR analysis estimated that a 1000 base-pair increase in TL increases risk of lung cancer [odds ratio (OR) = 1.41, 95% confidence interval (CI): 1.20-1.65] and lung adenocarcinoma (OR = 1.92, 95% CI: 1.51-2.22), but not squamous lung carcinoma (OR = 1.04, 95% CI: 0.83-1.29) or head and neck cancers (OR = 0.90, 95% CI: 0.70-1.05). Mediation analysis of the 5p15.33 instrument indicated an absence of direct effects on lung cancer risk (OR = 1.00, 95% CI: 0.95-1.04). Analysis of distinct 5p15.33 susceptibility variants estimated that TL mediates up to 40% of the observed associations with lung cancer risk.

CONCLUSIONS

Our findings support a causal role for long telomeres in lung cancer aetiology, particularly for adenocarcinoma, and demonstrate that telomere maintenance partially mediates the lung cancer susceptibility conferred by 5p15.33 loci.

Genetic interaction analysis among oncogenesis-related genes revealed novel genes and networks in lung cancer development

The development of cancer is driven by the accumulation of many oncogenesis-related genetic alterations and tumorigenesis is triggered by complex networks of involved genes rather than independent actions. To explore the epistasis existing among oncogenesis-related genes in lung cancer development, we conducted pairwise genetic interaction analyses among 35,031 SNPs from 2027 oncogenesis-related genes. The genotypes from three independent genome-wide association studies including a total of 24,037 lung cancer patients and 20,401 healthy controls with Caucasian ancestry were analyzed in the study. Using a two-stage study design including discovery and replication studies, and stringent Bonferroni correction for multiple statistical analysis, we identified significant genetic interactions between SNPs in RGL1:RAD51B (OR=0.44, p value=3.27x10-11 in overall lung cancer and OR=0.41, p value=9.71x10-11 in non-small cell lung cancer), SYNE1:RNF43 (OR=0.73, p value=1.01x10-12 in adenocarcinoma) and FHIT:TSPAN8 (OR=1.82, p value=7.62x10-11 in squamous cell carcinoma) in our analysis. None of these genes have been identified from previous main effect association studies in lung cancer. Further eQTL gene expression analysis in lung tissues provided information supporting the functional role of the identified epistasis in lung tumorigenesis. Gene set enrichment analysis revealed potential pathways and gene networks underlying molecular mechanisms in overall lung cancer as well as histology subtypes development. Our results provide evidence that genetic interactions between oncogenesis-related genes play an important role in lung tumorigenesis and epistasis analysis, combined with functional annotation, provides a valuable tool for uncovering functional novel susceptibility genes that contribute to lung cancer development by interacting with other modifier genes.

Fine mapping of MHC region in lung cancer highlights independent susceptibility loci by ethnicity

The basis for associations between lung cancer and major histocompatibility complex genes is not completely understood. Here the authors further consider genetic variation within the MHC region in lung cancer patients and identify independent associations within HLA genes that explain MHC lung cancer associations in Europeans and Asian populations.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA–tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

Positional integration of lung adenocarcinoma susceptibility loci with primary human alveolar epithelial cell epigenomes

AIM

To identify functional lung adenocarcinoma (LUAD) risk SNPs.

MATERIALS & METHODS

Eighteen validated LUAD risk SNPs (p ≤ 5 × 10^-8) and 930 SNPs in high linkage disequilibrium (r² > 0.5) were integrated with epigenomic information from primary human alveolar epithelial cells. Enhancer-associated SNPs likely affecting transcription factor-binding sites were predicted. Three SNPs were functionally investigated using luciferase assays, expression quantitative trait loci and cancer-specific expression.

RESULTS

Forty-seven SNPs mapped to putative enhancers; 11 located to open chromatin. Of these, seven altered predicted transcription factor-binding motifs. Rs6942067 showed allele-specific luciferase expression and expression quantitative trait loci analysis indicates that it influences expression of DCBLD1, a gene that encodes an unknown membrane protein and is overexpressed in LUAD.

CONCLUSION

Integration of candidate LUAD risk SNPS with epigenomic marks from normal alveolar epithelium identified numerous candidate functional LUAD risk SNPs including rs6942067, which appears to affect DCBLD1 expression. Data deposition: Data are provided in GEO record GSE84273.

Multi-Omics Analysis Reveals a HIF Network and Hub Gene EPAS1 Associated with Lung Adenocarcinoma

Recent technological advancements have permitted high-throughput measurement of the human genome, epigenome, metabolome, transcriptome, and proteome at the population level. We hypothesized that subsets of genes identified from omic studies might have closely related biological functions and thus might interact directly at the network level. Therefore, we conducted an integrative analysis of multi-omic datasets of non-small cell lung cancer (NSCLC) to search for association patterns beyond the genome and transcriptome. A large, complex, and robust gene network containing well-known lung cancer-related genes, including EGFR and TERT, was identified from combined gene lists for lung adenocarcinoma. Members of the hypoxia-inducible factor (HIF) gene family were at the center of this network. Subsequent sequencing of network hub genes within a subset of samples from the Transdisciplinary Research in Cancer of the Lung-International Lung Cancer Consortium (TRICL-ILCCO) consortium revealed a SNP (rs12614710) in EPAS1 associated with NSCLC that reached genome-wide significance (OR = 1.50; 95% CI: 1.31-1.72; p = 7.75 × 10^-9). Using imputed data, we found that this SNP remained significant in the entire TRICL-ILCCO consortium (p = .03). Additional functional studies are warranted to better understand interrelationships among genetic polymorphisms, DNA methylation status, and EPAS1 expression.

Fine mapping in TERT-CLPTM1L region identified three independent lung cancer susceptibility signals: A large-scale multi-ethnic population study

Genome-wide association studies (GWAS) and fine mapping studies have identified multiple lung cancer susceptibility variants in TERT-CLPTM1L region. However, it is still unclear about the relationship between these risk variants and the independent lung cancer risk signals in this region. Therefore, we evaluated the independent susceptibility signals for lung cancer and explored the potential functional variants in this region. Sequential conditional analysis was used to detect the independent susceptibility loci based on four lung cancer GWAS datasets with 12 843 lung cases and 12 639 controls. Comprehensively functional annotations were performed for each independent signal. Three independent susceptibility signals were identified in multi-ethnic population. For the first signal, rs2736100 showed the most significant association with lung cancer risk (C > A, OR = 0.82, 95%CI: 0.79-0.85, P = 1.98 × 10^-25 ). Rs36019446 was the top-ranked site (A > G, OR = 0.88, 95%CI: 0.84-0.92, P = 1.74 × 10^-9 ) in the second signal. For the third signal, rs326048 was the leading SNP (A > G, OR = 0.91, 95%CI: 0.87-0.95, P = 1.38 × 10^-5 ). The following subgroup analysis found the same three loci among Asian population. Further, we compared the difference between various subgroup populations. Functional annotations revealed that rs2736100, rs27996 (r²  = 0.85 with rs36019446) and rs326049 (r²  = 0.73 with rs326048) could be potential functional variants in these three risk signals, respectively. In conclusion, although multiple variants have been found associated with lung cancer risk in TERT-CLPTM1L region, our findings indicated that there are three independent lung cancer susceptibility signals in this region.

Genome-wide interaction study of smoking behavior and non-small cell lung cancer risk in Caucasian population

Non-small cell lung cancer is the most common type of lung cancer. Both environmental and genetic risk factors contribute to lung carcinogenesis. We conducted a genome-wide interaction analysis between single nucleotide polymorphisms (SNPs) and smoking status (never- versus ever-smokers) in a European-descent population. We adopted a two-step analysis strategy in the discovery stage: we first conducted a case-only interaction analysis to assess the relationship between SNPs and smoking behavior using 13336 non-small cell lung cancer cases. Candidate SNPs with P-value <0.001 were further analyzed using a standard case-control interaction analysis including 13970 controls. The significant SNPs with P-value <3.5 × 10-5 (correcting for multiple tests) from the case-control analysis in the discovery stage were further validated using an independent replication dataset comprising 5377 controls and 3054 non-small cell lung cancer cases. We further stratified the analysis by histological subtypes. Two novel SNPs, rs6441286 and rs17723637, were identified for overall lung cancer risk. The interaction odds ratio and meta-analysis P-value for these two SNPs were 1.24 with 6.96 × 10-7 and 1.37 with 3.49 × 10-7, respectively. In addition, interaction of smoking with rs4751674 was identified in squamous cell lung carcinoma with an odds ratio of 0.58 and P-value of 8.12 × 10-7. This study is by far the largest genome-wide SNP-smoking interaction analysis reported for lung cancer. The three identified novel SNPs provide potential candidate biomarkers for lung cancer risk screening and intervention. The results from our study reinforce that gene-smoking interactions play important roles in the etiology of lung cancer and account for part of the missing heritability of this disease.

Fine mapping of chromosome 5p15.33 identifies novel lung cancer susceptibility loci in Han Chinese

Genome-wide association studies in European and Asian populations have consistently identified chromosome 5p15.33 as a lung cancer susceptibility region. To investigate further the genetic architecture of common variants in this region, we conducted a two-stage fine-mapping analysis discovered by targeted resequencing of 200 cases and 300 controls individually, and validated in multiethnic lung cancer Genome wide association studies (GWASs) with 12,843 cases and 12,639 controls. Two independent variants were identified in approximate conditional analysis with GCTA and consistently validated in lung cancer GWASs in both Asian and European populations. These were rs10054203 in TERT (resequencing: OR = 1.69, p = 2.70 × 10^-4 ; validation: OR = 1.34, p = 2.10 × 10^-23 for Asian, and OR = 1.09, p = 6.00 × 10^-3 for European), and rs397640 in CLPTM1L (resequencing: OR = 0.37, p = 1.19 × 10^-4 ; validation: OR = 0.75, p = 5.89 × 10^-8 for Asian, and OR = 0.90, p = 2.40 × 10^-2 for European). Expression quantitative trait loci analysis showed the risk allele (C) of rs10054203 was significantly associated with lower mRNA expression of CTD-2245Ef15.3 (p = 0.019) and Tubulin Polymerization-Promoting Protein (TPPP, p = 0.031) in 167 lung tissues. In conclusion, in this largest and first resequencing-based fine-mapping analysis of 5p15.33 region in Han Chinese, we identified two novel variants associated with lung cancer susceptibility. Further validation studies and functional work is required to confirm the roles of the newly discovered variants.

Association between GWAS-identified lung adenocarcinoma susceptibility loci and EGFR mutations in never-smoking Asian women, and comparison with findings from Western populations

To evaluate associations by EGFR mutation status for lung adenocarcinoma risk among never-smoking Asian women, we conducted a meta-analysis of 11 loci previously identified in genome-wide association studies (GWAS). Genotyping in an additional 10,780 never-smoking cases and 10,938 never-smoking controls from Asia confirmed associations with eight known single nucleotide polymorphisms (SNPs). Two new signals were observed at genome-wide significance (P < 5 × 10-8), namely, rs7216064 (17q24.3, BPTF), for overall lung adenocarcinoma risk, and rs3817963 (6p21.3, BTNL2) which is specific to cases with EGFR mutations. In further sub-analyses by EGFR status, rs9387478 (ROS1/DCBLD1) and rs2179920 (HLA-DPB1) showed stronger estimated associations in EGFR-positive compared to EGFR-negative cases. Comparison of the overall associations with published results in Western populations revealed that the majority of these findings were distinct, underscoring the importance of distinct contributing factors for smoking and non-smoking lung cancer. Our results extend the catalogue of regions associated with lung adenocarcinoma in non-smoking Asian women and highlight the importance of how the germline could inform risk for specific tumour mutation patterns, which could have important translational implications.

The OncoArray Consortium: A Network for Understanding the Genetic Architecture of Common Cancers

BACKGROUND

Common cancers develop through a multistep process often including inherited susceptibility. Collaboration among multiple institutions, and funding from multiple sources, has allowed the development of an inexpensive genotyping microarray, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits.

METHODS

The OncoArray can be genotyped using a novel technology developed by Illumina to facilitate efficient genotyping. The consortium developed standard approaches for selecting SNPs for study, for quality control of markers, and for ancestry analysis. The array was genotyped at selected sites and with prespecified replicate samples to permit evaluation of genotyping accuracy among centers and by ethnic background.

RESULTS

The OncoArray consortium genotyped 447,705 samples. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis.

CONCLUSIONS

Results from these analyses will enable researchers to identify new susceptibility loci, perform fine-mapping of new or known loci associated with either single or multiple cancers, assess the degree of overlap in cancer causation and pleiotropic effects of loci that have been identified for disease-specific risk, and jointly model genetic, environmental, and lifestyle-related exposures.

IMPACT

Ongoing analyses will shed light on etiology and risk assessment for many types of cancer. Cancer Epidemiol Biomarkers Prev; 26(1); 126-35. ©2016 AACR.

Three new pancreatic cancer susceptibility signals identified on chromosomes 1q32.1, 5p15.33 and 8q24.21

Genome-wide association studies (GWAS) have identified common pancreatic cancer susceptibility variants at 13 chromosomal loci in individuals of European descent. To identify new susceptibility variants, we performed imputation based on 1000 Genomes (1000G) Project data and association analysis using 5,107 case and 8,845 control subjects from 27 cohort and case-control studies that participated in the PanScan I-III GWAS. This analysis, in combination with a two-staged replication in an additional 6,076 case and 7,555 control subjects from the PANcreatic Disease ReseArch (PANDoRA) and Pancreatic Cancer Case-Control (PanC4) Consortia uncovered 3 new pancreatic cancer risk signals marked by single nucleotide polymorphisms (SNPs) rs2816938 at chromosome 1q32.1 (per allele odds ratio (OR) = 1.20, P = 4.88x10 -15), rs10094872 at 8q24.21 (OR = 1.15, P = 3.22x10 -9) and rs35226131 at 5p15.33 (OR = 0.71, P = 1.70x10 -8). These SNPs represent independent risk variants at previously identified pancreatic cancer risk loci on chr1q32.1 ( NR5A2), chr8q24.21 ( MYC) and chr5p15.33 ( CLPTM1L- TERT) as per analyses conditioned on previously reported susceptibility variants. We assessed expression of candidate genes at the three risk loci in histologically normal ( n = 10) and tumor ( n = 8) derived pancreatic tissue samples and observed a marked reduction of NR5A2 expression (chr1q32.1) in the tumors (fold change -7.6, P = 5.7x10 -8). This finding was validated in a second set of paired ( n = 20) histologically normal and tumor derived pancreatic tissue samples (average fold change for three NR5A2 isoforms -31.3 to -95.7, P = 7.5x10 -4-2.0x10 -3). Our study has identified new susceptibility variants independently conferring pancreatic cancer risk that merit functional follow-up to identify target genes and explain the underlying biology.

The causal relevance of body mass index in different histological types of lung cancer: A Mendelian randomization study

Body mass index (BMI) is inversely associated with lung cancer risk in observational studies, even though it increases the risk of several other cancers, which could indicate confounding by tobacco smoking or reverse causality. We used the two-sample Mendelian randomization (MR) approach to circumvent these limitations of observational epidemiology by constructing a genetic instrument for BMI, based on results from the GIANT consortium, which was evaluated in relation to lung cancer risk using GWAS results on 16,572 lung cancer cases and 21,480 controls. Results were stratified by histological subtype, smoking status and sex. An increase of one standard deviation (SD) in BMI (4.65 Kg/m(2)) raised the risk for lung cancer overall (OR = 1.13; P = 0.10). This was driven by associations with squamous cell (SQ) carcinoma (OR = 1.45; P = 1.2 × 10(-3)) and small cell (SC) carcinoma (OR = 1.81; P = 0.01). An inverse trend was seen for adenocarcinoma (AD) (OR = 0.82; P = 0.06). In stratified analyses, a 1 SD increase in BMI was inversely associated with overall lung cancer in never smokers (OR = 0.50; P = 0.02). These results indicate that higher BMI may increase the risk of certain types of lung cancer, in particular SQ and SC carcinoma.