DNA hypomethylation in certain genes is associated with tobacco exposure but it is unknown whether these methylation changes translate into increased lung cancer risk. In an epigenome-wide study of DNA from pre-diagnostic blood samples from 132 case-control pairs in the NOWAC cohort, we observe that the most significant associations with lung cancer risk are for cg05575921 in AHRR (OR for 1 s.d.=0.37, 95% CI: 0.31-0.54, P-value=3.3 × 10(-11)) and cg03636183 in F2RL3 (OR for 1 s.d.=0.40, 95% CI: 0.31-0.56, P-value=3.9 × 10(-10)), previously shown to be strongly hypomethylated in smokers. These associations remain significant after adjustment for smoking and are confirmed in additional 664 case-control pairs tightly matched for smoking from the MCCS, NSHDS and EPIC HD cohorts. The replication and mediation analyses suggest that residual confounding is unlikely to explain the observed associations and that hypomethylation of these CpG sites may mediate the effect of tobacco on lung cancer risk.

Genome-wide association studies have identified polymorphisms linked to both smoking exposure and risk of lung cancer. The degree to which lung cancer risk is driven by increased smoking, genetics, or gene-environment interactions is not well understood.

Some single nucleotide polymorphisms (SNPs) are known to modify the risk of developing certain diseases or the reaction to drugs. Due to next generation sequencing methods the number of known human SNPs has grown. Not all SNPs lead to a modified protein, which may be the origin of a disease. Therefore, the recognition of functional SNPs is needed. Because most SNP annotation tools look for SNPs which lead to an amino acid exchange or a premature stop, we designed a new tool called AASsites which searches for SNPs which modify splicing.

Recent genome-wide association studies (GWASs) have identified common genetic variants at 5p15.33, 6p21-6p22 and 15q25.1 associated with lung cancer risk. Several other genetic regions including variants of CHEK2 (22q12), TP53BP1 (15q15) and RAD52 (12p13) have been demonstrated to influence lung cancer risk in candidate- or pathway-based analyses. To identify novel risk variants for lung cancer, we performed a meta-analysis of 16 GWASs, totaling 14 900 cases and 29 485 controls of European descent. Our data provided increased support for previously identified risk loci at 5p15 (P = 7.2 × 10(-16)), 6p21 (P = 2.3 × 10(-14)) and 15q25 (P = 2.2 × 10(-63)). Furthermore, we demonstrated histology-specific effects for 5p15, 6p21 and 12p13 loci but not for the 15q25 region. Subgroup analysis also identified a novel disease locus for squamous cell carcinoma at 9p21 (CDKN2A/p16(INK4A)/p14(ARF)/CDKN2B/p15(INK4B)/ANRIL; rs1333040, P = 3.0 × 10(-7)) which was replicated in a series of 5415 Han Chinese (P = 0.03; combined analysis, P = 2.3 × 10(-8)). This large analysis provides additional evidence for the role of inherited genetic susceptibility to lung cancer and insight into biological differences in the development of the different histological types of lung cancer.

For better lung cancer diagnosis and therapy, early detection markers of tumor dissemination are urgently needed, as most lung cancers do not show symptoms until extensive metastasis formation has already taken place. Our previous studies showed that in non-small cell lung cancer (NSCLC) early tumor dissemination is associated with a loss of chromosome 4q12-q32 and the presence of disseminated tumor cells (DTC) in the bone marrow. In order to identify the potential target gene in this region, a screen for methylation-dependent expression was performed. Lung cancer cell lines showing a loss of 4q as well as a normal bronchial epithelial cell line as control were treated with 5-aza-2'-deoxycytidine (5-aza-CdR) followed by expression profiling. Seven genes within the 4q target region, which have been associated with a positive DTC status before were found to be regulated by hypermethylation. QRT-PCR in an independent sample set identified HERC5 as a potential target gene. Quantitative methylation analysis of these lung tissue samples revealed that HERC5 promoter hypermethylation was significantly associated with positive DTC status (p = 0.020) and occurrence of brain metastases (p = 0.015). In addition, hypermethylation of the HERC5 promoter in NSCLC was identified as a predictor for poor survival for Stage I adenocarcinoma patients (p = 0.022) and also for poor overall survival in metastatic lung cancer patients (p = 0.028). In conclusion, HERC5 may function as a prognostic marker and is associated with tumor dissemination in lung cancer.

Use of the diabetes type II drug Metformin is associated with a moderately lowered risk of cancer incidence in numerous tumor entities. Studying the molecular changes associated with the tumor-suppressive action of Metformin we found that the oncogene SOX4, which is upregulated in solid tumors and associated with poor prognosis, was induced by Wnt/β-catenin signaling and blocked by Metformin. Wnt signaling inhibition by Metformin was surprisingly specific for cancer cells. Unraveling the underlying specificity, we identified Metformin and other Mitochondrial Complex I (MCI) inhibitors as inducers of intracellular acidification in cancer cells. We demonstrated that acidification triggers the unfolded protein response to induce the global transcriptional repressor DDIT3, known to block Wnt signaling. Moreover, our results suggest that intracellular acidification universally inhibits Wnt signaling. Based on these findings, we combined MCI inhibitors with H+ ionophores, to escalate cancer cells into intracellular hyper-acidification and ATP depletion. This treatment lowered intracellular pH both in vitro and in a mouse xenograft tumor model, depleted cellular ATP, blocked Wnt signaling, downregulated SOX4, and strongly decreased stemness and viability of cancer cells. Importantly, the inhibition of Wnt signaling occurred downstream of β-catenin, encouraging applications in treatment of cancers caused by APC and β-catenin mutations.

Genome-wide association studies (GWAS) identified the chromosome 15q25.1 locus as a leading susceptibility region for lung cancer. However, the pathogenic pathways, through which susceptibility SNPs within chromosome 15q25.1 affects lung cancer risk, have not been explored. We analyzed three cohorts with GWAS data consisting 42,901 individuals and lung expression quantitative trait loci (eQTL) data on 409 individuals to identify and validate the underlying pathways and to investigate the combined effect of genes from the identified susceptibility pathways. The KEGG neuroactive ligand receptor interaction pathway, two Reactome pathways, and 22 Gene Ontology terms were identified and replicated to be significantly associated with lung cancer risk, with P values less than 0.05 and FDR less than 0.1. Functional annotation of eQTL analysis results showed that the neuroactive ligand receptor interaction pathway and gated channel activity were involved in lung cancer risk. These pathways provide important insights for the etiology of lung cancer.

Persistent activation of hedgehog (HH)/GLI signaling accounts for the development of basal cell carcinoma (BCC), a very frequent nonmelanoma skin cancer with rising incidence. Targeting HH/GLI signaling by approved pathway inhibitors can provide significant therapeutic benefit to BCC patients. However, limited response rates, development of drug resistance, and severe side effects of HH pathway inhibitors call for improved treatment strategies such as rational combination therapies simultaneously inhibiting HH/GLI and cooperative signals promoting the oncogenic activity of HH/GLI. In this study, we identified the interleukin-6 (IL6) pathway as a novel synergistic signal promoting oncogenic HH/GLI via STAT3 activation. Mechanistically, we provide evidence that signal integration of IL6 and HH/GLI occurs at the level of cis-regulatory sequences by co-binding of GLI and STAT3 to common HH-IL6 target gene promoters. Genetic inactivation of Il6 signaling in a mouse model of BCC significantly reduced in vivo tumor growth by interfering with HH/GLI-driven BCC proliferation. Our genetic and pharmacologic data suggest that combinatorial HH-IL6 pathway blockade is a promising approach to efficiently arrest cancer growth in BCC patients.

Genome-wide association studies have identified hundreds of genetic variants associated with specific cancers. A few of these risk regions have been associated with more than one cancer site; however, a systematic evaluation of the associations between risk variants for other cancers and lung cancer risk has yet to be performed.

Treatment of non-small cell lung cancer (NSCLC) with radiotherapy or chemoradiotherapy is often accompanied by the development of esophagitis and pneumonitis. Identifying patients who might be at increased risk for normal tissue toxicity would help in determination of the optimal radiation dose to avoid these events. We profiled 59 single nucleotide polymorphisms (SNPs) from 37 inflammation-related genes in 173 NSCLC patients with stage IIIA/IIIB (dry) disease who were treated with definitive radiation or chemoradiation. For esophagitis risk, nine SNPs were associated with a 1.5- to 4-fold increase in risk, including three PTGS2 (COX2) variants: rs20417 (HR:1.93, 95% CI:1.10-3.39), rs5275 (HR:1.58, 95% CI:1.09-2.27), and rs689470 (HR:3.38, 95% CI:1.09-10.49). Significantly increased risk of pneumonitis was observed for patients with genetic variation in the proinflammatory genes IL1A, IL8, TNF, TNFRSF1B, and MIF. In contrast, NOS3:rs1799983 displayed a protective effect with a 45% reduction in pneumonitis risk (HR:0.55, 95% CI:0.31-0.96). Pneumonitis risk was also modulated by polymorphisms in anti-inflammatory genes, including genetic variation in IL13. rs20541 and rs180925 each resulted in increased risk (HR:2.95, 95% CI:1.14-7.63 and HR:3.23, 95% CI:1.03-10.18, respectively). The cumulative effect of these SNPs on risk was dose-dependent, as evidenced by a significantly increased risk of either toxicity with an increasing number of risk genotypes (P<0.001). These results suggest that genetic variations among inflammation pathway genes may modulate the development of radiation-induced toxicity and, ultimately, help in identifying patients who are at an increased likelihood for such events.

Apoptotic capacity (AC) in primary lymphocytes may be a marker for cancer susceptibility, and functional single nucleotide polymorphisms (SNPs) in genes involved in apoptotic pathways may modulate cellular AC in response to DNA damage. To further examine the correlation between apoptotic genotypes and phenotype, we genotyped 14 published SNPs in 11 apoptosis-related genes (i.e., p53, Bcl-2, BAX, CASP9, DR4, Fas, FasL, CASP8, CASP10, CASP3, and CASP7) and assessed the AC in response to benzo[a]pyrene-7,8-9,10-diol epoxide (BPDE) in cultured primary lymphocytes from 172 cancer-free subjects. We found that among these 14 SNPs, R72P, intron 3 16-bp del/ins, and intron 6 G>A in p53, -938C>A in Bcl-2, and I522L in CASP10 were significant predictors of the BPDE-induced lymphocytic AC in single-locus analysis. In the combined analysis of the three p53 variants, we found that the individuals with the diplotypes carrying 0-1 copy of the common p53 R-del-G haplotype had higher AC values compared to other genotypes. Although the study size may not have the statistical power to detect the role of other SNPs in AC, our findings suggest that some SNPs in genes involved in the intrinsic apoptotic pathway may modulate lymphocytic AC in response to BPDE exposure in the general population. Larger studies are needed to validate these findings for further studying individual susceptibility to cancer and other apoptosis-related diseases.

Lung cancer etiology is multifactorial, and growing evidence has indicated that long non-coding RNAs (lncRNAs) are important players in lung carcinogenesis. We performed a large-scale meta-analysis of 690,564 SNPs in 15,531 autosomal lncRNAs by using datasets from six previously published genome-wide association studies (GWASs) from the Transdisciplinary Research in Cancer of the Lung (TRICL) consortium in populations of European ancestry. Previously unreported significant SNPs (P value < 1 × 10-7) were further validated in two additional independent lung cancer GWAS datasets from Harvard University and deCODE. In the final meta-analysis of all eight GWAS datasets with 17,153 cases and 239,337 controls, a novel risk SNP rs114020893 in the lncRNA NEXN-AS1 region at 1p31.1 remained statistically significant (odds ratio = 1.17; 95% confidence interval = 1.11-1.24; P = 8.31 × 10-9). In further in silico analysis, rs114020893 was predicted to change the secondary structure of the lncRNA. Our finding indicates that SNP rs114020893 of NEXN-AS1 at 1p31.1 may contribute to lung cancer susceptibility.

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

Analyses of genome-wide association study (GWAS) data have revealed that detectable genetic mosaicism involving large (>2 Mb) structural autosomal alterations occurs in a fraction of individuals. We present results for a set of 24,849 genotyped individuals (total GWAS set II [TGSII]) in whom 341 large autosomal abnormalities were observed in 168 (0.68%) individuals. Merging data from the new TGSII set with data from two prior reports (the Gene-Environment Association Studies and the total GWAS set I) generated a large dataset of 127,179 individuals; we then conducted a meta-analysis to investigate the patterns of detectable autosomal mosaicism (n = 1,315 events in 925 [0.73%] individuals). Restricting to events >2 Mb in size, we observed an increase in event frequency as event size decreased. The combined results underscore that the rate of detectable mosaicism increases with age (p value = 5.5 × 10(-31)) and is higher in men (p value = 0.002) but lower in participants of African ancestry (p value = 0.003). In a subset of 47 individuals from whom serial samples were collected up to 6 years apart, complex changes were noted over time and showed an overall increase in the proportion of mosaic cells as age increased. Our large combined sample allowed for a unique ability to characterize detectable genetic mosaicism involving large structural events and strengthens the emerging evidence of non-random erosion of the genome in the aging population.

PARK2, a gene associated with Parkinson disease, is a tumor suppressor in human malignancies. Here, we show that c.823C>T (p.Arg275Trp), a germline mutation in PARK2, is present in a family with eight cases of lung cancer. The resulting amino acid change, p.Arg275Trp, is located in the highly conserved RING finger 1 domain of PARK2, which encodes an E3 ubiquitin ligase. Upon further analysis, the c.823C>T mutation was detected in three additional families affected by lung cancer. The effect size for PARK2 c.823C>T (odds ratio = 5.24) in white individuals was larger than those reported for variants from lung cancer genome-wide association studies. These data implicate this PARK2 germline mutation as a genetic susceptibility factor for lung cancer. Our results provide a rationale for further investigations of this specific mutation and gene for evaluation of the possibility of developing targeted therapies against lung cancer in individuals with PARK2 variants by compensating for the loss-of-function effect caused by the associated variation.

Three genetic loci for lung cancer risk have been identified by genome-wide association studies (GWAS), but inherited susceptibility to specific histologic types of lung cancer is not well established. We conducted a GWAS of lung cancer and its major histologic types, genotyping 515,922 single-nucleotide polymorphisms (SNPs) in 5739 lung cancer cases and 5848 controls from one population-based case-control study and three cohort studies. Results were combined with summary data from ten additional studies, for a total of 13,300 cases and 19,666 controls of European descent. Four studies also provided histology data for replication, resulting in 3333 adenocarcinomas (AD), 2589 squamous cell carcinomas (SQ), and 1418 small cell carcinomas (SC). In analyses by histology, rs2736100 (TERT), on chromosome 5p15.33, was associated with risk of adenocarcinoma (odds ratio [OR]=1.23, 95% confidence interval [CI]=1.13-1.33, p=3.02x10(-7)), but not with other histologic types (OR=1.01, p=0.84 and OR=1.00, p=0.93 for SQ and SC, respectively). This finding was confirmed in each replication study and overall meta-analysis (OR=1.24, 95% CI=1.17-1.31, p=3.74x10(-14) for AD; OR=0.99, p=0.69 and OR=0.97, p=0.48 for SQ and SC, respectively). Other previously reported association signals on 15q25 and 6p21 were also refined, but no additional loci reached genome-wide significance. In conclusion, a lung cancer GWAS identified a distinct hereditary contribution to adenocarcinoma.

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<0.0001). However, average levels of global methylation as measured by the Illumina technology were not significantly associated with menarcheal age. In locus by locus comparative analyses, only one CpG site had significantly different methylation depending on the menarcheal age category examined, but this finding was not replicated by pyrosequencing in an independent data set. This study suggests a link between age at menarche and genome wide DNA methylation, and the difference in results between the two arrays suggests that repetitive element methylation has a role in the association. Epigenetic changes may be modulated by menarcheal age, or the association may be a mirror of other important changes in early life that have a detectable effect on both methylation levels and menarcheal age.

This study correlated the extent of induced in vitro chromosomal damage, assessed by the mutagen sensitivity assay, with genotypes of the X-ray repair cross complementing group 1 (XRCC1) gene, which encodes for a base excision repair protein. There are two common polymorphisms that cause amino acid substitutions in XRCC1, one at codon 194 in exon 6 and another at codon 399 in exon 10. We genotyped these two polymorphisms in 524 healthy subjects and performed mutagen sensitivity assays using both bleomycin and benzo[a]pyrene-diol-epoxide (BPDE) as challenge mutagens. Our results showed that individuals with the wildtype exon 6 Arg/Arg exhibited significantly higher values of chromosomal breaks per cell (b/c) than those with one or two variant Trp alleles (P=0.005 for bleomycin and P=0.05 for BPDE). For the exon 10 polymorphism, subjects who were Gln/Gln homozygotes had higher b/c than did those with other genotypes, with evidence of a gene dosage effect. When we combined the two polymorphic sites and used the exon 6 Arg/Trp and Trp/Trp and exon 10 Arg/Arg genotypes as the reference category, these differences were enhanced for bleomycin sensitivity (P for trend = 0.032), but not for BPDE sensitivity (P for trend = 0.821). These data are biologically plausible since codon 399 is located within the BRCA1 C-terminus functional domain and codon 194 is in the linker region of the XRCC1 N-terminal functional domain. To our knowledge, this is the largest study conducted evaluating the functional relevance of these polymorphisms.

The potential for genome-wide association studies to relate phenotypes to specific genetic variation is greatly increased when data can be combined or compared across multiple studies. To facilitate replication and validation across studies, RTI International (Research Triangle Park, North Carolina) and the National Human Genome Research Institute (Bethesda, Maryland) are collaborating on the consensus measures for Phenotypes and eXposures (PhenX) project. The goal of PhenX is to identify 15 high-priority, well-established, and broadly applicable measures for each of 21 research domains. PhenX measures are selected by working groups of domain experts using a consensus process that includes input from the scientific community. The selected measures are then made freely available to the scientific community via the PhenX Toolkit. Thus, the PhenX Toolkit provides the research community with a core set of high-quality, well-established, low-burden measures intended for use in large-scale genomic studies. PhenX measures will have the most impact when included at the experimental design stage. The PhenX Toolkit also includes links to standards and resources in an effort to facilitate data harmonization to legacy data. Broad acceptance and use of PhenX measures will promote cross-study comparisons to increase statistical power for identifying and replicating variants associated with complex diseases and with gene-gene and gene-environment interactions.

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.