Search

Abstract The alarming increase in global rates of metabolic diseases (MetDs) and their association with cancer risk renders them a considerable burden on our society. The interplay of environmental and genetic factors in causing MetDs may be reflected in DNA methylation patterns, particularly at non-canonical (non-B) DNA structures, such as G-quadruplexes (G4s) or R-loops. To gain insight into the mechanisms of MetD progression, we focused on DNA methylation and functional analyses on intragenic regions of two MetD risk genes, the glucokinase (GCK) exon 7 and the transmembrane 6 superfamily 2 (TM6SF2) intron 2-exon 3 boundary, which harbor non-B DNA motifs for G4s and R-loops. Pyrosequencing of 148 blood samples from a nested cohort study revealed significant differential methylation in GCK and TM6SF2 in MetD patients versus healthy controls. Furthermore, these regions harbor hypervariable and differentially methylated CpGs also in hepatocellular carcinoma versus normal tissue samples from The Cancer Genome Atlas (TCGA). Permanganate/S1 nuclease footprinting with direct adapter ligation (PDAL-Seq), native polyacrylamide DNA gel electrophoresis and circular dichroism (CD) spectroscopy revealed the formation of G4 structures in these regions and demonstrated that their topology and stability is affected by DNA methylation. Detailed analyses including histone marks, chromatin conformation capture data, and luciferase reporter assays, highlighted the cell-type specific regulatory function of the target regions. Based on our analyses, we hypothesize that changes in DNA methylation lead to topological changes, especially in GCK exon 7, and cause the activation of alternative regulatory elements or potentially play a role in alternative splicing. Our analyses provide a new view on the mechanisms underlying the progression of MetDs and their link to hepatocellular carcinomas, unveiling non-B DNA structures as important key players already in early disease stages.

Abstract Background Although polygenic risk score (PRS) has emerged as a promising tool for predicting cancer risk from genome-wide association studies (GWAS), the individual-level accuracy of lung cancer PRS and the extent to which its impact on subsequent clinical applications remains largely unexplored. Methods Lung cancer PRSs and confidence/credible interval (CI) were constructed using two statistical approaches for each individual: (1) the weighted sum of 16 GWAS-derived significant SNP loci and the CI through the bootstrapping method (PRS-16-CV) and (2) LDpred2 and the CI through posteriors sampling (PRS-Bayes), among 17,166 lung cancer cases and 12,894 controls with European ancestry from the International Lung Cancer Consortium. Individuals were classified into different genetic risk subgroups based on the relationship between their own PRS mean/PRS CI and the population level threshold. Results Considerable variances in PRS point estimates at the individual level were observed for both methods, with an average standard deviation (s.d.) of 0.12 for PRS-16-CV and a much larger s.d. of 0.88 for PRS-Bayes. Using PRS-16-CV, only 25.0% of individuals with PRS point estimates in the lowest decile of PRS and 16.8% in the highest decile have their entire 95% CI fully contained in the lowest and highest decile, respectively, while PRS-Bayes was unable to find any eligible individuals. Only 19% of the individuals were concordantly identified as having high genetic risk (> 90th percentile) using the two PRS estimators. An increased relative risk of lung cancer comparing the highest PRS percentile to the lowest was observed when taking the CI into account (OR = 2.73, 95% CI: 2.12–3.50, P-value = 4.13 × 10−15) compared to using PRS-16-CV mean (OR = 2.23, 95% CI: 1.99–2.49, P-value = 5.70 × 10−46). Improved risk prediction performance with higher AUC was consistently observed in individuals identified by PRS-16-CV CI, and the best performance was achieved by incorporating age, gender, and detailed smoking pack-years (AUC: 0.73, 95% CI = 0.72–0.74). Conclusions Lung cancer PRS estimates using different methods have modest correlations at the individual level, highlighting the importance of considering individual-level uncertainty when evaluating the practical utility of PRS.

Abstract Background Epigenetic alterations are a near-universal feature of human malignancy and have been detected in malignant cells as well as in easily accessible specimens such as blood and urine. These findings offer promising applications in cancer detection, subtyping, and treatment monitoring. However, much of the current evidence is based on findings in retrospective studies and may reflect epigenetic patterns that have already been influenced by the onset of the disease. Methods Studying breast cancer, we established genome-scale DNA methylation profiles of prospectively collected buffy coat samples (n = 702) from a case–control study nested within the EPIC-Heidelberg cohort using reduced representation bisulphite sequencing (RRBS). Results We observed cancer-specific DNA methylation events in buffy coat samples. Increased DNA methylation in genomic regions associated with SURF6 and REXO1/CTB31O20.3 was linked to the length of time to diagnosis in the prospectively collected buffy coat DNA from individuals who subsequently developed breast cancer. Using machine learning methods, we piloted a DNA methylation-based classifier that predicted case–control status in a held-out validation set with 76.5% accuracy, in some cases up to 15 years before clinical diagnosis of the disease. Conclusions Taken together, our findings suggest a model of gradual accumulation of cancer-associated DNA methylation patterns in peripheral blood, which may be detected long before clinical manifestation of cancer. Such changes may provide useful markers for risk stratification and, ultimately, personalized cancer prevention.

Summary: Background: Tumour-promoting inflammation is a “hallmark” of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear. Methods: We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis. Findings: We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10–1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20–1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53–0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88–0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48–4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk. Interpretation: Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated. Funding: Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Abstract Lung adenocarcinoma is the most common type of lung cancer. Known risk variants explain only a small fraction of lung adenocarcinoma heritability. Here, we conducted a two-stage genome-wide association study of lung adenocarcinoma of East Asian ancestry (21,658 cases and 150,676 controls; 54.5% never-smokers) and identified 12 novel susceptibility variants, bringing the total number to 28 at 25 independent loci. Transcriptome-wide association analyses together with colocalization studies using a Taiwanese lung expression quantitative trait loci dataset (n = 115) identified novel candidate genes, including FADS1 at 11q12 and ELF5 at 11p13. In a multi-ancestry meta-analysis of East Asian and European studies, four loci were identified at 2p11, 4q32, 16q23, and 18q12. At the same time, most of our findings in East Asian populations showed no evidence of association in European populations. In our studies drawn from East Asian populations, a polygenic risk score based on the 25 loci had a stronger association in never-smokers vs. individuals with a history of smoking (Pinteraction = 0.0058). These findings provide new insights into the etiology of lung adenocarcinoma in individuals from East Asian populations, which could be important in developing translational applications.

Cystic fibrosis (CF) is a monogenic disease, characterized by massive chronic lung inflammation. The observed variability in clinical phenotypes in monozygotic CF twins is likely associated with the extent of inflammation. This study sought to investigate inflammation-related aberrant DNA methylation in CF twins and to determine to what extent acquired methylation changes may be associated with lung cancer. Blood-based genome-wide DNA methylation analysis was performed to compare the DNA methylomes of monozygotic twins, from the European CF Twin and Sibling Study with various degrees of disease severity. Putatively inflammation-related and differentially methylated positions were selected from a large lung cancer case-control study and investigated in blood by targeted bisulphite next-generation-sequencing. An inflammation-related locus located in the Plakophilin-3 (PKP3) gene was functionally analysed regarding promoter and enhancer activity in presence and absence of methylation using luciferase reporter assays. We confirmed in a unique cohort that monozygotic twins, even if clinically discordant, have only minor differences in global DNA methylation patterns and blood cell composition. Further, we determined the most differentially methylated positions, a high proportion of which are blood cell-type-specific, whereas others may be acquired and thus have potential relevance in the context of inflammation as lung cancer risk factors. We identified a sequence in the gene body of PKP3 which is hypermethylated in blood from CF twins with severe phenotype and highly variably methylated in lung cancer patients and controls, independent of known clinical parameters, and showed that this region exhibits methylation-dependent promoter activity in lung epithelial cells.

Abstract Limited efforts have been made in assessing the effect of genome-wide profiling of RNA splicing-related variation on lung cancer risk. In the present study, we first identified RNA splicing-related genetic variants linked to lung cancer in a genome-wide profiling analysis and then conducted a two-stage (discovery and replication) association study in populations of European ancestry. Discovery and validation were conducted sequentially with a total of 29,266 cases and 56,450 controls from both the Transdisciplinary Research in Cancer of the Lung and the International Lung Cancer Consortium as well as the OncoArray database. For those variants identified as significant in the two datasets, we further performed stratified analyses by smoking status and histological type and investigated their effects on gene expression and potential regulatory mechanisms. We identified three genetic variants significantly associated with lung cancer risk: rs329118 in JADE2 (P = 8.80E−09), rs2285521 in GGA2 (P = 4.43E−08), and rs198459 in MYRF (P = 1.60E−06). The combined effects of all three SNPs were more evident in lung squamous cell carcinomas (P = 1.81E−08, P = 6.21E−08, and P = 7.93E−04, respectively) than in lung adenocarcinomas and in ever smokers (P = 9.80E−05, P = 2.70E−04, and P = 2.90E−05, respectively) than in never smokers. Gene expression quantitative trait analysis suggested a role for the SNPs in regulating transcriptional expression of the corresponding target genes. In conclusion, we report that three RNA splicing-related genetic variants contribute to lung cancer susceptibility in European populations. However, additional validation is needed, and specific splicing mechanisms of the target genes underlying the observed associations also warrants further exploration.

Genome-wide association studies (GWAS) have identified SNPs linked with lung cancer risk. Our aim was to discover the genes, non-coding RNAs, and regulatory elements within GWAS-identified risk regions that are deregulated in non-small cell lung carcinoma (NSCLC) to identify novel, clinically targetable genes and mechanisms in carcinogenesis. A targeted bisulphite-sequencing approach was used to comprehensively investigate DNA methylation changes occurring within lung cancer risk regions in 17 NSCLC and adjacent normal tissue pairs. We report differences in differentially methylated regions between adenocarcinoma and squamous cell carcinoma. Among the minimal regions found to be differentially methylated in at least 50% of the patients, 7 candidates were replicated in 2 independent cohorts (n = 27 and n = 87) and the potential of 6 as methylation-dependent regulatory elements was confirmed by functional assays. This study contributes to understanding the pathways implicated in lung cancer initiation and progression, and provides new potential targets for cancer treatment.

Abstract Background Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility. Aim To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent. Methods Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups. Results No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR = 1.20; 95% CI 1.13–1.27; p = 5.6 × 10–10) and never smokers (e.g., maker rs1133683 Axin2; OR = 1.27; 95% CI 1.19–1.35; p = 1.0 × 10–12). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants. Conclusions The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers.

Genetic activation of hedgehog/glioma‐associated oncogene homolog (HH/GLI) signaling causes basal cell carcinoma (BCC), a very frequent nonmelanoma skin cancer. Small molecule targeting of the essential HH effector Smoothened (SMO) has proven an effective therapy of BCC, though the frequent development of drug resistance poses major challenges to anti‐HH treatments. In light of recent breakthroughs in cancer immunotherapy, we analyzed the possible immunosuppressive mechanisms in HH/GLI‐induced BCC in detail. Using a genetic mouse model of BCC, we identified profound differences in the infiltration of BCC lesions with cells of the adaptive and innate immune system. Epidermal activation of Hh/Gli signaling led to an accumulation of immunosuppressive regulatory T cells, and to an increased expression of immune checkpoint molecules including programmed death (PD)‐1/PD‐ligand 1. Anti‐PD‐1 monotherapy, however, did not reduce tumor growth, presumably due to the lack of immunogenic mutations in common BCC mouse models, as shown by whole‐exome sequencing. BCC lesions also displayed a marked infiltration with neutrophils, the depletion of which unexpectedly promoted BCC growth. The study provides a comprehensive survey of and novel insights into the immune status of murine BCC and serves as a basis for the design of efficacious rational combination treatments. This study also underlines the need for predictive immunogenic mouse models of BCC to evaluate the efficacy of immunotherapeutic strategies in vivo.

In lung cancer, relatively few germline mutations are known to impact risk. Here the authors looked at rare variants in 39,146 individuals and find novel germline mutations associated with risk, as well as implicating ATM and a new candidate gene for lung cancer risk.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Exposure to a prenatal high-fat (HF) diet leads to an impaired metabolic phenotype in mouse offspring. The underlying mechanisms, however, are not yet fully understood. Therefore, this study investigated whether the impaired metabolic phenotype may be mediated through altered hepatic DNA methylation and gene expression. We showed that exposure to a prenatal HF diet altered the offspring's hepatic gene expression of pathways involved in lipid synthesis and uptake (SREBP), oxidative stress response [nuclear factor (erythroid-derived 2)-like 2 (Nrf2)], and cell proliferation. The downregulation of the SREBP pathway related to previously reported decreased hepatic lipid uptake and postprandial hypertriglyceridemia in the offspring exposed to the prenatal HF diet. The upregulation of the Nrf2 pathway was associated with increased oxidative stress levels in offspring livers. The prenatal HF diet also induced hypermethylation of transcription factor (TF) binding sites upstream of lipin 1 (Lpin1), a gene involved in lipid metabolism. Furthermore, DNA methylation of Lpin1 TF binding sites correlated with mRNA expression of Lpin1. These findings suggest that the effect of a prenatal HF diet on the adult offspring's metabolic phenotype are regulated by changes in hepatic gene expression and DNA methylation.

Similarities in cancers can be studied to interrogate their etiology. Here, the authors use genome-wide association study summary statistics from six cancer types based on 296,215 cases and 301,319 controls of European ancestry, showing that solid tumours arising from different tissues share a degree of common germline genetic basis.

Squamous cell carcinomas (SqCC) of the aerodigestive tract have similar etiological risk factors. Although genetic risk variants for individual cancers have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. To identify novel and pleotropic SqCC risk variants, we performed a meta-analysis of GWAS data on lung SqCC (LuSqCC), oro/pharyngeal SqCC (OSqCC), laryngeal SqCC (LaSqCC) and esophageal SqCC (ESqCC) cancers, totaling 13,887 cases and 61,961 controls of European ancestry. We identified one novel genome-wide significant (Pmeta

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

The chromosome 15q25.1 locus is a leading susceptibility region for lung cancer. Here, the authors interrogate three GWAS cohorts with 42,901 individuals to investigate potential pathological pathways such as gated channel activity and neuroactive ligand receptor interaction in lung cancer etiology.

Abstract The T-cell protein tyrosine phosphatase (TCPTP) pathway consists of signaling events mediated by TCPTP. Mutations and genetic variants of some genes in the TCPTP pathway are associated with lung cancer risk and survival. In the present study, we first investigated associations of 5,162 single nucleotide polymorphisms (SNPs) in 43 genes of this TCPTP pathway with lung cancer risk by using summary data of six published genome-wide association studies (GWAS) of 12,160 cases and 16,838 controls. We identified 11 independent SNPs in eight genes after correction for multiple comparisons by a false discovery rate

Abstract Hypomethylating agents (HMAs) have been widely used over the last decade, approved for use in myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML) and acute myeloid leukemia (AML). The proposed central mechanism of action of HMAs, is the reversal of aberrant methylation in tumor cells, thus reactivating CpG-island promoters and leading to (re)expression of tumor suppressor genes. Recent investigations into the mode of action of azacitidine (AZA) and decitabine (DAC) have revealed new molecular mechanisms that impinge on tumor immunity via induction of an interferon response, through activation of endogenous retroviral elements (ERVs) that are normally epigenetically silenced. Although the global demethylation of DNA by HMAs can induce anti-tumor effects, it can also upregulate the expression of inhibitory immune checkpoint receptors and their ligands, resulting in secondary resistance to HMAs. Recent studies have, however, suggested that this could be exploited to prime or (re)sensitize tumors to immune checkpoint inhibitor therapies. In recent years, immune checkpoints have been targeted by novel therapies, with the aim of (re)activating the host immune system to specifically eliminate malignant cells. Antibodies blocking checkpoint receptors have been FDA-approved for some solid tumors and a plethora of clinical trials testing these and other checkpoint inhibitors are under way. This review will discuss AZA and DAC novel mechanisms of action resulting from the re-expression of pathologically hypermethylated promoters of gene sets that are related to interferon signaling, antigen presentation and inflammation. We also review new insights into the molecular mechanisms of action of transient, low-dose HMAs on various tumor types and discuss the potential of new treatment options and combinations.

Background: Recent meta-analyses show that individuals with high risk variants in CHRNA5 on chromosome 15q25 are likely to develop lung cancer earlier than those with low-risk genotypes. The same high-risk genetic variants also predict nicotine dependence and delayed smoking cessation. It is unclear whether smoking cessation confers the same benefits in terms of lung cancer risk reduction for those who possess CHRNA5 risk variants versus those who do not. Methods: Meta-analyses examined the association between smoking cessation and lung cancer risk in 15 studies of individuals with European ancestry who possessed varying rs16969968 genotypes (N = 12,690 ever smokers, including 6988 cases of lung cancer and 5702 controls) in the International Lung Cancer Consortium. Results: Smoking cessation (former vs. current smokers) was associated with a lower likelihood of lung cancer (OR = 0.48, 95%CI = 0.30–0.75, p = 0.0015). Among lung cancer patients, smoking cessation was associated with a 7-year delay in median age of lung cancer diagnosis (HR = 0.68, 95%CI = 0.61–0.77, p = 4.9 ∗ 10–10). The CHRNA5 rs16969968 risk genotype (AA) was associated with increased risk and earlier diagnosis for lung cancer, but the beneficial effects of smoking cessation were very similar in those with and without the risk genotype. Conclusion: We demonstrate that quitting smoking is highly beneficial in reducing lung cancer risks for smokers regardless of their CHRNA5 rs16969968 genetic risk status. Smokers with high-risk CHRNA5 genotypes, on average, can largely eliminate their elevated genetic risk for lung cancer by quitting smoking- cutting their risk of lung cancer in half and delaying its onset by 7 years for those who develop it. These results: 1) underscore the potential value of smoking cessation for all smokers, 2) suggest that CHRNA5 rs16969968 genotype affects lung cancer diagnosis through its effects on smoking, and 3) have potential value for framing preventive interventions for those who smoke.

INTRODUCTION:Gene-set analysis (GSA) is an approach using the results of single-marker genome-wide association studies when investigating pathways as a whole with respect to the genetic basis of a disease. METHODS:We performed a meta-analysis of seven GSAs for lung cancer, applying the method META-GSA. Overall, the information taken from 11,365 cases and 22,505 controls from within the TRICL/ILCCO consortia was used to investigate a total of 234 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. RESULTS:META-GSA reveals the systemic lupus erythematosus KEGG pathway hsa05322, driven by the gene region 6p21-22, as also implicated in lung cancer (p = 0.0306). This gene region is known to be associated with squamous cell lung carcinoma. The most important genes driving the significance of this pathway belong to the genomic areas HIST1-H4L, -1BN, -2BN, -H2AK, -H4K and C2/C4A/C4B. Within these areas, the markers most significantly associated with LC are rs13194781 (located within HIST12BN) and rs1270942 (located between C2 and C4A). CONCLUSIONS:We have discovered a pathway currently marked as specific to systemic lupus erythematosus as being significantly implicated in lung cancer. The gene region 6p21-22 in this pathway appears to be more extensively associated with lung cancer than previously assumed. Given wide-stretched linkage disequilibrium to the area APOM/BAG6/MSH5, there is currently simply not enough information or evidence to conclude whether the potential pleiotropy of lung cancer and systemic lupus erythematosus is spurious, biological, or mediated. Further research into this pathway and gene region will be necessary.

BACKGROUND:Assessing the relationship between lung cancer and metabolic conditions is challenging because of the confounding effect of tobacco. Mendelian randomization (MR), or the use of genetic instrumental variables to assess causality, may help to identify the metabolic drivers of lung cancer. METHODS AND FINDINGS:We identified genetic instruments for potential metabolic risk factors and evaluated these in relation to risk using 29,266 lung cancer cases (including 11,273 adenocarcinomas, 7,426 squamous cell and 2,664 small cell cases) and 56,450 controls. The MR risk analysis suggested a causal effect of body mass index (BMI) on lung cancer risk for two of the three major histological subtypes, with evidence of a risk increase for squamous cell carcinoma (odds ratio (OR) [95% confidence interval (CI)] = 1.20 [1.01-1.43] and for small cell lung cancer (OR [95%CI] = 1.52 [1.15-2.00]) for each standard deviation (SD) increase in BMI [4.6 kg/m2]), but not for adenocarcinoma (OR [95%CI] = 0.93 [0.79-1.08]) (Pheterogeneity = 4.3x10-3). Additional analysis using a genetic instrument for BMI showed that each SD increase in BMI increased cigarette consumption by 1.27 cigarettes per day (P = 2.1x10-3), providing novel evidence that a genetic susceptibility to obesity influences smoking patterns. There was also evidence that low-density lipoprotein cholesterol was inversely associated with lung cancer overall risk (OR [95%CI] = 0.90 [0.84-0.97] per SD of 38 mg/dl), while fasting insulin was positively associated (OR [95%CI] = 1.63 [1.25-2.13] per SD of 44.4 pmol/l). Sensitivity analyses including a weighted-median approach and MR-Egger test did not detect other pleiotropic effects biasing the main results. CONCLUSIONS:Our results are consistent with a causal role of fasting insulin and low-density lipoprotein cholesterol in lung cancer etiology, as well as for BMI in squamous cell and small cell carcinoma. The latter relation may be mediated by a previously unrecognized effect of obesity on smoking behavior.

BackgroundObservational studies examining associations between adult height and risk of colorectal, prostate, and lung cancers have generated mixed results. We conducted meta-analyses using data from prospective cohort studies and further carried out Mendelian randomization analyses, using height-associated genetic variants identified in a genome-wide association study (GWAS), to evaluate the association of adult height with these cancers.Methods and findingsA systematic review of prospective studies was conducted using the PubMed, Embase, and Web of Science databases. Using meta-analyses, results obtained from 62 studies were summarized for the association of a 10-cm increase in height with cancer risk. Mendelian randomization analyses were conducted using summary statistics obtained for 423 genetic variants identified from a recent GWAS of adult height and from a cancer genetics consortium study of multiple cancers that included 47,800 cases and 81,353 controls. For a 10-cm increase in height, the summary relative risks derived from the meta-analyses of prospective studies were 1.12 (95% CI 1.10, 1.15), 1.07 (95% CI 1.05, 1.10), and 1.06 (95% CI 1.02, 1.11) for colorectal, prostate, and lung cancers, respectively. Mendelian randomization analyses showed increased risks of colorectal (odds ratio [OR] = 1.58, 95% CI 1.14, 2.18) and lung cancer (OR = 1.10, 95% CI 1.00, 1.22) associated with each 10-cm increase in genetically predicted height. No association was observed for prostate cancer (OR = 1.03, 95% CI 0.92, 1.15). Our meta-analysis was limited to published studies. The sample size for the Mendelian randomization analysis of colorectal cancer was relatively small, thus affecting the precision of the point estimate.ConclusionsOur study provides evidence for a potential causal association of adult height with the risk of colorectal and lung cancers and suggests that certain genetic factors and biological pathways affecting adult height may also affect the risk of these cancers.

Gene-set analysis (GSA) methods are used as complementary approaches to genome-wide association studies (GWASs). The single marker association estimates of a predefined set of genes are either contrasted with those of all remaining genes or with a null non-associated background. To pool the p-values from several GSAs, it is important to take into account the concordance of the observed patterns resulting from single marker association point estimates across any given gene set. Here we propose an enhanced version of Fisher's inverse χ2-method META-GSA, however weighting each study to account for imperfect correlation between association patterns.We investigated the performance of META-GSA by simulating GWASs with 500 cases and 500 controls at 100 diallelic markers in 20 different scenarios, simulating different relative risks between 1 and 1.5 in gene sets of 10 genes. Wilcoxon's rank sum test was applied as GSA for each study. We found that META-GSA has greater power to discover truly associated gene sets than simple pooling of the p-values, by e.g. 59% versus 37%, when the true relative risk for 5 of 10 genes was assume to be 1.5. Under the null hypothesis of no difference in the true association pattern between the gene set of interest and the set of remaining genes, the results of both approaches are almost uncorrelated. We recommend not relying on p-values alone when combining the results of independent GSAs.We applied META-GSA to pool the results of four case-control GWASs of lung cancer risk (Central European Study and Toronto/Lunenfeld-Tanenbaum Research Institute Study; German Lung Cancer Study and MD Anderson Cancer Center Study), which had already been analyzed separately with four different GSA methods (EASE; SLAT, mSUMSTAT and GenGen). This application revealed the pathway GO0015291 "transmembrane transporter activity" as significantly enriched with associated genes (GSA-method: EASE, p = 0.0315 corrected for multiple testing). Similar results were found for GO0015464 "acetylcholine receptor activity" but only when not corrected for multiple testing (all GSA-methods applied; p ≈ 0.02).

Reproductive factors have been linked to both breast cancer and DNA methylation, suggesting methylation as an important mechanism by which reproductive factors impact on disease risk. However, few studies have investigated the link between reproductive factors and DNA methylation in humans. Genome-wide methylation in peripheral blood lymphocytes of 376 healthy women from the prospective EPIC study was investigated using LUminometric Methylation Assay (LUMA). Also, methylation of 458877 CpG sites was additionally investigated in an independent group of 332 participants of the EPIC-Italy sub-cohort, using the Infinium HumanMethylation 450 BeadChip. Multivariate logistic regression and linear models were used to investigate the association between reproductive risk factors and genome wide and CpG-specific DNA methylation, respectively. Menarcheal age was inversely associated with global DNA methylation as measured with LUMA. For each yearly increase in age at menarche, the risk of having genome wide methylation below median level was increased by 32% (OR:1.32, 95%CI:1.14-1.53). When age at menarche was treated as a categorical variable, there was an inverse dose-response relationship with LUMA methylation levels (OR(12-14 vs. ≤11 yrs):1.78, 95%CI:1.01-3.17 and OR(≥15 vs. ≤11 yrs):4.59, 95%CI:2.04-10.33; P for trend

Genetic case-control association studies often include data on clinical covariates, such as body mass index (BMI), smoking status, or age, that may modify the underlying genetic risk of case or control samples. For example, in type 2 diabetes, odds ratios for established variants estimated from low-BMI cases are larger than those estimated from high-BMI cases. An unanswered question is how to use this information to maximize statistical power in case-control studies that ascertain individuals on the basis of phenotype (case-control ascertainment) or phenotype and clinical covariates (case-control-covariate ascertainment). While current approaches improve power in studies with random ascertainment, they often lose power under case-control ascertainment and fail to capture available power increases under case-control-covariate ascertainment. We show that an informed conditioning approach, based on the liability threshold model with parameters informed by external epidemiological information, fully accounts for disease prevalence and non-random ascertainment of phenotype as well as covariates and provides a substantial increase in power while maintaining a properly controlled false-positive rate. Our method outperforms standard case-control association tests with or without covariates, tests of gene x covariate interaction, and previously proposed tests for dealing with covariates in ascertained data, with especially large improvements in the case of case-control-covariate ascertainment. We investigate empirical case-control studies of type 2 diabetes, prostate cancer, lung cancer, breast cancer, rheumatoid arthritis, age-related macular degeneration, and end-stage kidney disease over a total of 89,726 samples. In these datasets, informed conditioning outperforms logistic regression for 115 of the 157 known associated variants investigated (P-value = 1 × 10(-9)). The improvement varied across diseases with a 16% median increase in χ(2) test statistics and a commensurate increase in power. This suggests that applying our method to existing and future association studies of these diseases may identify novel disease loci.

Pathway analysis has been proposed as a complement to single SNP analyses in GWAS. This study compared pathway analysis methods using two lung cancer GWAS data sets based on four studies: one a combined data set from Central Europe and Toronto (CETO); the other a combined data set from Germany and MD Anderson (GRMD). We searched the literature for pathway analysis methods that were widely used, representative of other methods, and had available software for performing analysis. We selected the programs EASE, which uses a modified Fishers Exact calculation to test for pathway associations, GenGen (a version of Gene Set Enrichment Analysis (GSEA)), which uses a Kolmogorov-Smirnov-like running sum statistic as the test statistic, and SLAT, which uses a p-value combination approach. We also included a modified version of the SUMSTAT method (mSUMSTAT), which tests for association by averaging χ(2) statistics from genotype association tests. There were nearly 18000 genes available for analysis, following mapping of more than 300,000 SNPs from each data set. These were mapped to 421 GO level 4 gene sets for pathway analysis. Among the methods designed to be robust to biases related to gene size and pathway SNP correlation (GenGen, mSUMSTAT and SLAT), the mSUMSTAT approach identified the most significant pathways (8 in CETO and 1 in GRMD). This included a highly plausible association for the acetylcholine receptor activity pathway in both CETO (FDR≤0.001) and GRMD (FDR = 0.009), although two strong association signals at a single gene cluster (CHRNA3-CHRNA5-CHRNB4) drive this result, complicating its interpretation. Few other replicated associations were found using any of these methods. Difficulty in replicating associations hindered our comparison, but results suggest mSUMSTAT has advantages over the other approaches, and may be a useful pathway analysis tool to use alongside other methods such as the commonly used GSEA (GenGen) approach.

Radiotherapy exerts part of its antineoplastic effect by generating oxidative stress, therefore genetic variation in oxidative stress-related enzymes may influence survival of rectal cancer patients. We hypothesized that genetic polymorphisms associated with higher amounts of reactive oxygen species (ROS) that exaggerate cytotoxic activity could improve survival after radiotherapy. We followed 114 rectal cancer patients who received radiotherapy for an average of 42.5 months. Associations between genotypes (GSTP1, GSTM1, GSTT1, CAT, MnSOD, MPO and eNOS) and overall survival were assessed using Kaplan-Meier curves and Cox proportional hazards regression. As hypothesized, patients carrying low ROS producing eNOS Glu298Asp asparagine allele showed an increased hazard of death compared to homozygous carriers of the glutamine allele (hazard ratio (HR): 2.10, 95% confidence interval (CI): 1.01–4.38). However, carriers of low ROS producing MPO G463A A allele had a decreased hazard of death compared to patients homozygous for the G allele (HR: 0.44, 95% CI: 0.21–0.93) although patients homozygous for the A allele had a slightly increased hazard (HR: 1.12, 95% CI: 0.25–5.08). This explorative study provides first results and highlights the need for further, larger studies to investigate association between genetic variation in oxidative stress genes and survival of rectal cancer patients who received radiotherapy.