Identification of genetically predicted DNA methylation markers associated with non-small cell lung cancer risk among 34,964 cases and 448,579 controls

BACKGROUND

Although the associations between genetic variations and lung cancer risk have been explored, the epigenetic consequences of DNA methylation in lung cancer development are largely unknown. Here, the genetically predicted DNA methylation markers associated with non-small cell lung cancer (NSCLC) risk by a two-stage case-control design were investigated.

METHODS

The genetic prediction models for methylation levels based on genetic and methylation data of 1595 subjects from the Framingham Heart Study were established. The prediction models were applied to a fixed-effect meta-analysis of screening data sets with 27,120 NSCLC cases and 27,355 controls to identify the methylation markers, which were then replicated in independent data sets with 7844 lung cancer cases and 421,224 controls. Also performed was a multi-omics functional annotation for the identified CpGs by integrating genomics, epigenomics, and transcriptomics and investigation of the potential regulation pathways.

RESULTS

Of the 29,894 CpG sites passing the quality control, 39 CpGs associated with NSCLC risk (Bonferroni-corrected p ≤ 1.67 × 10-6 ) were originally identified. Of these, 16 CpGs remained significant in the validation stage (Bonferroni-corrected p ≤ 1.28 × 10-3 ), including four novel CpGs. Multi-omics functional annotation showed nine of 16 CpGs were potentially functional biomarkers for NSCLC risk. Thirty-five genes within a 1-Mb window of 12 CpGs that might be involved in regulatory pathways of NSCLC risk were identified.

CONCLUSIONS

Sixteen promising DNA methylation markers associated with NSCLC were identified. Changes of the methylation level at these CpGs might influence the development of NSCLC by regulating the expression of genes nearby.

PLAIN LANGUAGE SUMMARY

The epigenetic consequences of DNA methylation in lung cancer development are still largely unknown. This study used summary data of large-scale genome-wide association studies to investigate the associations between genetically predicted levels of methylation biomarkers and non-small cell lung cancer risk at the first time. This study looked at how well larotrectinib worked in adult patients with sarcomas caused by TRK fusion proteins. These findings will provide a unique insight into the epigenetic susceptibility mechanisms of lung cancer.

Lung Cancer in Ever- and Never-Smokers: Findings from Multi-Population GWAS Studies

BACKGROUND

Clinical, molecular, and genetic epidemiology studies displayed remarkable differences between ever- and never-smoking lung cancer.

METHODS

We conducted a stratified multi-population (European, East Asian, and African descent) association study on 44,823 ever-smokers and 20,074 never-smokers to identify novel variants that were missed in the non-stratified analysis. Functional analysis including expression quantitative trait loci (eQTL) colocalization and DNA damage assays, and annotation studies were conducted to evaluate the functional roles of the variants. We further evaluated the impact of smoking quantity on lung cancer risk for the variants associated with ever-smoking lung cancer.

RESULTS

Five novel independent loci, GABRA4, intergenic region 12q24.33, LRRC4C, LINC01088, and LCNL1 were identified with the association at two or three populations (P < 5 × 10-8). Further functional analysis provided multiple lines of evidence suggesting the variants affect lung cancer risk through excessive DNA damage (GABRA4) or cis-regulation of gene expression (LCNL1). The risk of variants from 12 independent regions, including the well-known CHRNA5, associated with ever-smoking lung cancer was evaluated for never-smokers, light-smokers (packyear ≤ 20), and moderate-to-heavy-smokers (packyear > 20). Different risk patterns were observed for the variants among the different groups by smoking behavior.

CONCLUSIONS

We identified novel variants associated with lung cancer in only ever- or never-smoking groups that were missed by prior main-effect association studies.

IMPACT

Our study highlights the genetic heterogeneity between ever- and never-smoking lung cancer and provides etiologic insights into the complicated genetic architecture of this deadly cancer.

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

An inter-laboratory effort to harmonize the cell-delivered in vitro dose of aerosolized materials

Air-liquid interface (ALI) lung cell models cultured on permeable transwell inserts are increasingly used for respiratory hazard assessment requiring controlled aerosolization and deposition of any material on ALI cells. The approach presented herein aimed to assess the transwell insert-delivered dose of aerosolized materials using the VITROCELL® Cloud12 system, a commercially available aerosol-cell exposure system. An inter-laboratory comparison study was conducted with seven European partners having different levels of experience with the VITROCELL® Cloud12. A standard operating procedure (SOP) was developed and applied by all partners for aerosolized delivery of materials, i.e., a water-soluble molecular substance (fluorescence-spiked salt) and two poorly soluble particles, crystalline silica quartz (DQ₁₂) and titanium dioxide nanoparticles (TiO₂ NM-105). The material dose delivered to transwell inserts was quantified with spectrofluorometry (fluorescein) and with the quartz crystal microbalance (QCM) integrated in the VITROCELL® Cloud12 system. The shape and agglomeration state of the deposited particles were confirmed with transmission electron microscopy (TEM). Inter-laboratory comparison of the device-specific performance was conducted in two steps, first for molecular substances (fluorescein-spiked salt), and then for particles. Device- and/or handling-specific differences in aerosol deposition of VITROCELL® Cloud12 systems were characterized in terms of the so-called deposition factor (DF), which allows for prediction of the transwell insert-deposited particle dose from the particle concentration in the aerosolized suspension. Albeit DF varied between the different labs from 0.39 to 0.87 (mean (coefficient of variation (CV)): 0.64 (28%)), the QCM of each VITROCELL® Cloud 12 system accurately measured the respective transwell insert-deposited dose. Aerosolized delivery of DQ₁₂ and TiO₂ NM-105 particles showed good linearity (R² > 0.95) between particle concentration of the aerosolized suspension and QCM-determined insert-delivered particle dose. The VITROCELL® Cloud 12 performance for DQ₁₂ particles was identical to that for fluorescein-spiked salt, i.e., the ratio of measured and salt-predicted dose was 1.0 (29%). On the other hand, a ca. 2-fold reduced dose was observed for TiO₂ NM-105 (0.54 (41%)), which was likely due to partial retention of TiO₂ NM-105 agglomerates in the vibrating mesh nebulizer of the VITROCELL® Cloud12. This inter-laboratory comparison demonstrates that the QCM integrated in the VITROCELL® Cloud 12 is a reliable tool for dosimetry, which accounts for potential variations of the transwell insert-delivered dose due to device-, handling- and/or material-specific effects. With the detailed protocol presented herein, all seven partner laboratories were able to demonstrate dose-controlled aerosolization of material suspensions using the VITROCELL® Cloud12 exposure system at dose levels relevant for observing in vitro hazard responses. This is an important step towards regulatory approved implementation of ALI lung cell cultures for in vitro hazard assessment of aerosolized materials.

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.

Gene-gene interaction of AhRwith and within the Wntcascade affects susceptibility to lung cancer

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.

Iam hiQ-a novel pair of accuracy indices for imputed genotypes

BACKGROUND

Imputation of untyped markers is a standard tool in genome-wide association studies to close the gap between directly genotyped and other known DNA variants. However, high accuracy with which genotypes are imputed is fundamental. Several accuracy measures have been proposed and some are implemented in imputation software, unfortunately diversely across platforms. In the present paper, we introduce Iam hiQ, an independent pair of accuracy measures that can be applied to dosage files, the output of all imputation software. Iam (imputation accuracy measure) quantifies the average amount of individual-specific versus population-specific genotype information in a linear manner. hiQ (heterogeneity in quantities of dosages) addresses the inter-individual heterogeneity between dosages of a marker across the sample at hand.

RESULTS

Applying both measures to a large case-control sample of the International Lung Cancer Consortium (ILCCO), comprising 27,065 individuals, we found meaningful thresholds for Iam and hiQ suitable to classify markers of poor accuracy. We demonstrate how Manhattan-like plots and moving averages of Iam and hiQ can be useful to identify regions enriched with less accurate imputed markers, whereas these regions would by missed when applying the accuracy measure info (implemented in IMPUTE2).

CONCLUSION

We recommend using Iam hiQ additional to other accuracy scores for variant filtering before stepping into the analysis of imputed GWAS data.

The leukocyte telomere length, single nucleotide polymorphisms near TERC gene and risk of COPD

Chronic obstructive pulmonary disease (COPD) is characterized by irreversible airflow obstruction and is associated with chronic local and systemic inflammation and oxidative stress. The enhanced oxidative stress and inflammation have been reported to affect telomere length (TL). Furthermore, a number of SNPs at loci encoding the main components of the telomerase genes, TERT and TERC have been shown to correlate with TL. We aimed to explore the leukocyte TL and genotypes for single nucleotide polymorphisms, rs12696304 (C > G) and rs10936599 (C > T) near TERC in COPD cases and matched healthy controls using q-PCR technologies. Successful assessment of TL was performed for 91 patients and 88 controls. The patients had shorter TL (17919.36 ± 1203.01 bp) compared to controls (21 271.48 ± 1891.36 bp) although not significant (p = 0.137). The TL did not associate with the gender, age, spirometric indexes, smoking habits but tended to correlate negatively with BMI (Rho = − 0.215, p = 0.076) in the controls, but not in COPD patients. The genotype frequencies of the SNPs rs12696304 and rs10936599 were compared between patients and controls and the odds ratios (OR) for developing COPD were calculated. The carriers of the common homozygous (CC) genotypes of the SNPs had higher risk for COPD, compared to carriers of the variants alleles (rs12696304 CG+GG vs. CC; OR: 0.615, 95% CI [0.424–0.894], p = 0.011 and for rs10936599 CT+TT vs. CC OR = 0.668, 95% CI [0.457–0.976], p = 0.044). Analysis on the combined effects of the TERCrs12696304 (C > G) and rs10936599 (C > T) genotypes, CC/CC genotype combination was associated with higher risk for COPD (p < 0.0001) and marginally lower FEV1% pr. in patients with GOLD II (p = 0.052). There was no association between the SNP genotypes and TL. In summary, our results suggest that COPD patients may have shorter TL, and rs12696304 and rs10936599 near TERC may affect the risk of COPD independently of TL.

Long-Term Exposure to Nanosized TiO2 Triggers Stress Responses and Cell Death Pathways in Pulmonary Epithelial Cells

There is little in vitro data available on long-term effects of TiO₂ exposure. Such data are important for improving the understanding of underlying mechanisms of adverse health effects of TiO₂. Here, we exposed pulmonary epithelial cells to two doses (0.96 and 1.92 µg/cm²) of TiO₂ for 13 weeks and effects on cell cycle and cell death mechanisms, i.e., apoptosis and autophagy were determined after 4, 8 and 13 weeks of exposure. Changes in telomere length, cellular protein levels and lipid classes were also analyzed at 13 weeks of exposure. We observed that the TiO₂ exposure increased the fraction of cells in G1-phase and reduced the fraction of cells in G2-phase, which was accompanied by an increase in the fraction of late apoptotic/necrotic cells. This corresponded with an induced expression of key apoptotic proteins i.e., BAD and BAX, and an accumulation of several lipid classes involved in cellular stress and apoptosis. These findings were further supported by quantitative proteome profiling data showing an increase in proteins involved in cell stress and genomic maintenance pathways following TiO₂ exposure. Altogether, we suggest that cell stress response and cell death pathways may be important molecular events in long-term health effects of TiO₂.

Characterization and toxicity evaluation of air-borne particles released by grinding from two dental resin composites in vitro

OBJECTIVE

The project aims to evaluate whether inhalation of particles released upon grinding of dental composites may pose a health hazard to dentists. The main objective of the study was to characterize the dust from polymer-based dental composites ground with different grain sized burs and investigate particle uptake and the potential cytotoxic effects in human bronchial cells.

METHODS

Polymerized blocks of two dental composites, Filtek™ Z250 and Filtek™ Z500 from 3M™ ESPE, were ground with super coarse (black) and fine (red) burs inside a glass chamber. Ultrafine airborne dust concentration and particle size distribution was measured real-time during grinding with a scanning mobility particle sizer (SMPS). Filter-collected airborne particles were characterized with dynamic light scattering (DLS) and scanning electron microscopy (SEM). Human bronchial epithelial cells (HBEC-3KT) were exposed to the dusts in dose-effect experiments. Toxicity was measured with lactate dehydrogenase (LDH) assay and cell counting kit-8 (CCK8). Cellular uptake was observed with transmission electron microscopy (TEM).

RESULTS

Airborne ultrafine particles showed that most particles were in the size range 15-35 nm (SMPS). SEM analysis proved that more than 80% of the particles have a minimum Feret diameter less than 1 μm. In solution (DLS), the particles have larger diameters and tend to agglomerate. Cell toxicity (LDH, CCK8) is shown after 48 h and 72 h exposure times and at the highest doses. TEM showed presence of the particles within the cell cytoplasm.

SIGNIFICANCE

Prolonged and frequent exposure through inhalation may have negative health implications for dentists.

Causal relationships between body mass index, smoking and lung cancer: Univariable and multivariable Mendelian randomization

At the time of cancer diagnosis, body mass index (BMI) is inversely correlated with lung cancer risk, which may reflect reverse causality and confounding due to smoking behavior. We used two-sample univariable and multivariable Mendelian randomization (MR) to estimate causal relationships of BMI and smoking behaviors on lung cancer and histological subtypes based on an aggregated genome-wide association studies (GWASs) analysis of lung cancer in 29 266 cases and 56 450 controls. We observed a positive causal effect for high BMI on occurrence of small-cell lung cancer (odds ratio (OR) = 1.60, 95% confidence interval (CI) = 1.24-2.06, P = 2.70 × 10-4 ). After adjustment of smoking behaviors using multivariable Mendelian randomization (MVMR), a direct causal effect on small cell lung cancer (ORMVMR = 1.28, 95% CI = 1.06-1.55, PMVMR = .011), and an inverse effect on lung adenocarcinoma (ORMVMR = 0.86, 95% CI = 0.77-0.96, PMVMR = .008) were observed. A weak increased risk of lung squamous cell carcinoma was observed for higher BMI in univariable Mendelian randomization (UVMR) analysis (ORUVMR = 1.19, 95% CI = 1.01-1.40, PUVMR = .036), but this effect disappeared after adjustment of smoking (ORMVMR = 1.02, 95% CI = 0.90-1.16, PMVMR = .746). These results highlight the histology-specific impact of BMI on lung carcinogenesis and imply mediator role of smoking behaviors in the association between BMI and lung cancer.

Genome-wide association meta-analysis identifies pleiotropic risk loci for aerodigestive squamous cell cancers

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<5x10-8) aerodigestive SqCC susceptibility loci in the 2q33.1 region (rs56321285, TMEM273). Additionally, three previously unknown loci reached suggestive significance (Pmeta<5x10-7): 1q32.1 (rs12133735, near MDM4), 5q31.2 (rs13181561, TMEM173) and 19p13.11 (rs61494113, ABHD8). Multiple previously identified loci for aerodigestive SqCC also showed evidence of pleiotropy in at least another SqCC site, these include: 4q23 (ADH1B), 6p21.33 (STK19), 6p21.32 (HLA-DQB1), 9p21.33 (CDKN2B-AS1) and 13q13.1(BRCA2). Gene-based association and gene set enrichment identified a set of 48 SqCC-related genes rel to DNA damage and epigenetic regulation pathways. Our study highlights the importance of cross-cancer analyses to identify pleiotropic risk loci of histology-related cancers arising at distinct anatomical sites.

A new efficient method to detect genetic interactions for lung cancer GWAS

BACKGROUND

Genome-wide association studies (GWAS) have proven successful in predicting genetic risk of disease using single-locus models; however, identifying single nucleotide polymorphism (SNP) interactions at the genome-wide scale is limited due to computational and statistical challenges. We addressed the computational burden encountered when detecting SNP interactions for survival analysis, such as age of disease-onset. To confront this problem, we developed a novel algorithm, called the Efficient Survival Multifactor Dimensionality Reduction (ES-MDR) method, which used Martingale Residuals as the outcome parameter to estimate survival outcomes, and implemented the Quantitative Multifactor Dimensionality Reduction method to identify significant interactions associated with age of disease-onset.

METHODS

To demonstrate efficacy, we evaluated this method on two simulation data sets to estimate the type I error rate and power. Simulations showed that ES-MDR identified interactions using less computational workload and allowed for adjustment of covariates. We applied ES-MDR on the OncoArray-TRICL Consortium data with 14,935 cases and 12,787 controls for lung cancer (SNPs = 108,254) to search over all two-way interactions to identify genetic interactions associated with lung cancer age-of-onset. We tested the best model in an independent data set from the OncoArray-TRICL data.

RESULTS

Our experiment on the OncoArray-TRICL data identified many one-way and two-way models with a single-base deletion in the noncoding region of BRCA1 (HR 1.24, P = 3.15 × 10-15), as the top marker to predict age of lung cancer onset.

CONCLUSIONS

From the results of our extensive simulations and analysis of a large GWAS study, we demonstrated that our method is an efficient algorithm that identified genetic interactions to include in our models to predict survival outcomes.

Effects of mild steel welding fume particles on pulmonary epithelial inflammation and endothelial activation

Welders have an increased risk for cardiovascular disease (CVD) following exposure to welding fumes. The underlying mechanisms are largely unknown; however, oxidative stress, systemic inflammation, and endothelial dysfunction have been suggested as contributing factors to particle-induced CVD. We investigated effects of mild steel welding fume (MSWF) on three target cell types: macrophages, pulmonary epithelial, and vascular endothelial cells. Cells were exposed to MSWF at nontoxic doses for 6 h/day, for five consecutive days. The expression of 40 genes involved in inflammation, fibrosis, and endothelial activation was analyzed. Moreover, changes in the reactive oxygen species production and migration capacity of cells were assessed. The expression of matrix metallopeptidase 1 (MMP1) was induced in both epithelial and endothelial cells following repeated exposure to MSWF. Although MMP1 is important in inflammatory responses in vivo, this effect was not concurrent with changes in the inflammatory status, cell proliferation, and migration capacities, nor did it induce oxidative stress in the cells. Thus, repeated exposure with low doses of MSWF was sufficient neither for inducing inflammatory stress in epithelial cells and macrophages nor for endothelial activation, and higher concentrations of MSWF or the nonparticle fraction of MSWF may be critical in causing the increased risk of CVD observed among welders.

Association Analysis of Driver Gene-Related Genetic Variants Identified Novel Lung Cancer Susceptibility Loci with 20,871 Lung Cancer Cases and 15,971 Controls

BACKGROUND

A substantial proportion of cancer driver genes (CDG) are also cancer predisposition genes. However, the associations between genetic variants in lung CDGs and the susceptibility to lung cancer have rarely been investigated.

METHODS

We selected expression-related single-nucleotide polymorphisms (eSNP) and nonsynonymous variants of lung CDGs, and tested their associations with lung cancer risk in two large-scale genome-wide association studies (20,871 cases and 15,971 controls of European descent). Conditional and joint association analysis was performed to identify independent risk variants. The associations of independent risk variants with somatic alterations in lung CDGs or recurrently altered pathways were investigated using data from The Cancer Genome Atlas (TCGA) project.

RESULTS

We identified seven independent SNPs in five lung CDGs that were consistently associated with lung cancer risk in discovery (P < 0.001) and validation (P < 0.05) stages. Among these loci, rs78062588 in TPM3 (1q21.3) was a new lung cancer susceptibility locus (OR = 0.86, P = 1.65 × 10-6). Subgroup analysis by histologic types further identified nine lung CDGs. Analysis of somatic alterations found that in lung adenocarcinomas, rs78062588[C] allele (TPM3 in 1q21.3) was associated with elevated somatic copy number of TPM3 (OR = 1.16, P = 0.02). In lung adenocarcinomas, rs1611182 (HLA-A in 6p22.1) was associated with truncation mutations of the transcriptional misregulation in cancer pathway (OR = 0.66, P = 1.76 × 10-3).

CONCLUSIONS

Genetic variants can regulate functions of lung CDGs and influence lung cancer susceptibility.

IMPACT

Our findings might help unravel biological mechanisms underlying lung cancer susceptibility.

Systematic analyses of regulatory variants in DNase I hypersensitive sites identified two novel lung cancer susceptibility loci

DNase I hypersensitive sites (DHS) are abundant in regulatory elements, such as promoter, enhancer and transcription factor binding sites. Many studies have revealed that disease-associated variants were concentrated in DHS-related regions. However, limited studies are available on the roles of DHS-related variants in lung cancer. In this study, we performed a large-scale case-control study with 20 871 lung cancer cases and 15 971 controls to evaluate the associations between regulatory genetic variants in DHS and lung cancer susceptibility. The expression quantitative trait loci (eQTL) analysis and pathway-enrichment analysis were performed to identify the possible target genes and pathways. In addition, we performed motif-based analysis to explore the lung-cancer-related motifs using sequence kernel association test. Two novel variants, rs186332 in 20q13.3 (C>T, odds ratio [OR] = 1.17, 95% confidence interval [95% CI]: 1.10-1.24, P = 8.45 × 10-7) and rs4839323 in 1p13.2 (T>C, OR = 0.92, 95% CI: 0.89-0.95, P = 1.02 × 10-6) showed significant association with lung cancer risk. The eQTL analysis suggested that these two SNPs might regulate the expression of MRGBP and SLC16A1, respectively. What's more, the expression of both MRGBP and SLC16A1 was aberrantly elevated in lung tumor tissues. The motif-based analysis identified 10 motifs related to the risk of lung cancer (P < 1.71 × 10-4). Our findings suggested that variants in DHS might modify lung cancer susceptibility through regulating the expression of surrounding genes. This study provided us a deeper insight into the roles of DHS-related genetic variants for lung cancer.

Immune-mediated genetic pathways resulting in pulmonary function impairment increase lung cancer susceptibility

Impaired lung function is often caused by cigarette smoking, making it challenging to disentangle its role in lung cancer susceptibility. Investigation of the shared genetic basis of these phenotypes in the UK Biobank and International Lung Cancer Consortium (29,266 cases, 56,450 controls) shows that lung cancer is genetically correlated with reduced forced expiratory volume in one second (FEV1: rg = 0.098, p = 2.3 × 10-8) and the ratio of FEV1 to forced vital capacity (FEV1/FVC: rg = 0.137, p = 2.0 × 10-12). Mendelian randomization analyses demonstrate that reduced FEV1 increases squamous cell carcinoma risk (odds ratio (OR) = 1.51, 95% confidence intervals: 1.21-1.88), while reduced FEV1/FVC increases the risk of adenocarcinoma (OR = 1.17, 1.01-1.35) and lung cancer in never smokers (OR = 1.56, 1.05-2.30). These findings support a causal role of pulmonary impairment in lung cancer etiology. Integrative analyses reveal that pulmonary function instruments, including 73 novel variants, influence lung tissue gene expression and implicate immune-related pathways in mediating the observed effects on lung carcinogenesis.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

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

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.

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.

Cellulose nanocrystals modulate alveolar macrophage phenotype and phagocytic function

Nanocellulose is a promising bio-nanomaterial with attractive properties suitable for multiple industrial applications. The increased use of nanocellulose may lead to occupational exposure and negative health outcomes. However, knowledge on its health effects is limited, and while nanocellulose exposure may induce acute inflammatory responses in the lung, the underlying mechanisms are unknown. Alveolar macrophages are key cells in alveolar particle clearance. Their activation and function may be affected by various particles. Here, we investigated the uptake of pristine cellulose nanocrystals (CNC), and their effects on alveolar macrophage polarization and biological function. CNC uptake enhanced the secretion of several cytokines but did not on its own induce a complete macrophage polarization. In presence of macrophage activators, such as LPS/IFNG and IL4/IL13, CNC exposure enhanced the expression of M1 phenotype markers and the secretion of pro-inflammatory cytokines and chemokines, while decreasing M2 markers. CNC exposure also affected the function of activated alveolar macrophages resulting in a prominent cytokine burst and altered phagocytic activity. In conclusion, CNC exposure may result in dysregulation of macrophage activation and function that are critical in inflammatory responses in the lung.

Elevated Platelet Count Appears to Be Causally Associated with Increased Risk of Lung Cancer: A Mendelian Randomization Analysis

BACKGROUND

Platelets are a critical element in coagulation and inflammation, and activated platelets are linked to cancer risk through diverse mechanisms. However, a causal relationship between platelets and risk of lung cancer remains unclear.

METHODS

We performed single and combined multiple instrumental variable Mendelian randomization analysis by an inverse-weighted method, in addition to a series of sensitivity analyses. Summary data for associations between SNPs and platelet count are from a recent publication that included 48,666 Caucasian Europeans, and the International Lung Cancer Consortium and Transdisciplinary Research in Cancer of the Lung data consisting of 29,266 cases and 56,450 controls to analyze associations between candidate SNPs and lung cancer risk.

RESULTS

Multiple instrumental variable analysis incorporating six SNPs showed a 62% increased risk of overall non-small cell lung cancer [NSCLC; OR, 1.62; 95% confidence interval (CI), 1.15-2.27; P = 0.005] and a 200% increased risk for small-cell lung cancer (OR, 3.00; 95% CI, 1.27-7.06; P = 0.01). Results showed only a trending association with NSCLC histologic subtypes, which may be due to insufficient sample size and/or weak effect size. A series of sensitivity analysis retained these findings.

CONCLUSIONS

Our findings suggest a causal relationship between elevated platelet count and increased risk of lung cancer and provide evidence of possible antiplatelet interventions for lung cancer prevention.

IMPACT

These findings provide a better understanding of lung cancer etiology and potential evidence for antiplatelet interventions for lung cancer prevention.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Genetic modifiers of radon-induced lung cancer risk: a genome-wide interaction study in former uranium miners

PURPOSE

Radon is a risk factor for lung cancer and uranium miners are more exposed than the general population. A genome-wide interaction analysis was carried out to identify genomic loci, genes or gene sets that modify the susceptibility to lung cancer given occupational exposure to the radioactive gas radon.

METHODS

Samples from 28 studies provided by the International Lung Cancer Consortium were pooled with samples of former uranium miners collected by the German Federal Office of Radiation Protection. In total, 15,077 cases and 13,522 controls, all of European ancestries, comprising 463 uranium miners were compared. The DNA of all participants was genotyped with the OncoArray. We fitted single-marker and in multi-marker models and performed an exploratory gene-set analysis to detect cumulative enrichment of significance in sets of genes.

RESULTS

We discovered a genome-wide significant interaction of the marker rs12440014 within the gene CHRNB4 (OR = 0.26, 95% CI 0.11-0.60, p = 0.0386 corrected for multiple testing). At least suggestive significant interaction of linkage disequilibrium blocks was observed at the chromosomal regions 18q21.23 (p = 1.2 × 10-6), 5q23.2 (p = 2.5 × 10-6), 1q21.3 (p = 3.2 × 10-6), 10p13 (p = 1.3 × 10-5) and 12p12.1 (p = 7.1 × 10-5). Genes belonging to the Gene Ontology term "DNA dealkylation involved in DNA repair" (GO:0006307; p = 0.0139) or the gene family HGNC:476 "microRNAs" (p = 0.0159) were enriched with LD-blockwise significance.

CONCLUSION

The well-established association of the genomic region 15q25 to lung cancer might be influenced by exposure to radon among uranium miners. Furthermore, lung cancer susceptibility is related to the functional capability of DNA damage signaling via ubiquitination processes and repair of radiation-induced double-strand breaks by the single-strand annealing mechanism.

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.

Identification of susceptibility pathways for the role of chromosome 15q25.1 in modifying lung cancer risk

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.

Inflammation in the pleural cavity following injection of multi-walled carbon nanotubes is dependent on their characteristics and the presence of IL-1 genes

Upon inhalation, multi-walled carbon nanotubes (MWCNTs) may reach the subpleura and pleural spaces, and induce pleural inflammation and/or mesothelioma in humans. However, the mechanisms of MWCNT-induced pathology after direct intrapleural injections are still only partly elucidated. In particular, a role of the proinflammatory interleukin-1 (IL-1) cytokines in pleural inflammation has so far not been published. We examined the MWCNT-induced pleural inflammation, gene expression abnormalities, and the modifying role of IL-1α and β cytokines following intrapleural injection of two types of MWCNTs (CNT-1 and CNT-2) compared with crocidolite asbestos in IL-1 wild-type (WT) and IL-1α/β KO (IL1-KO) mice. Histopathological examination of the pleura 28 days post-exposure revealed mesothelial cell hyperplasia, leukocyte infiltration, and fibrosis occurring in the CNT-1 (Mitsui-7)-exposed group. The pleura of these mice also showed the greatest changes in mRNA and miRNA expression levels, closely followed by CNT-2. In addition, the CNT-1-exposed group also presented the greatest infiltrations of leukocytes and proliferation of fibrous tissue. WT mice were more prone to development of sustained inflammation and fibrosis than IL1-KO mice. Prominent differences in genetic and epigenetic changes were also observed between the two genotypes. In conclusion, the fibrotic response to MWCNTs in the pleura depends on the particles' physico-chemical properties and on the presence or absence of the IL-1 genes. Furthermore, we found that CNT-1 was the most potent inducer of inflammatory responses, followed by CNT-2 and crocidolite asbestos.

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.

Characterization of the proteome and lipidome profiles of human lung cells after low dose and chronic exposure to multiwalled carbon nanotubes

The effects of long-term chronic exposure of human lung cells to multi-walled carbon nanotubes (MWCNT) and their impact upon cellular proteins and lipids were investigated. Since the lung is the major target organ, an in vitro normal bronchial epithelial cell line model was used. Additionally, to better mimic exposure to manufactured nanomaterials at occupational settings, cells were continuously exposed to two non-toxic and low doses of a MWCNT for 13-weeks. MWCNT-treatment increased ROS levels in cells without increasing oxidative DNA damage and resulted in differential expression of multiple anti- and pro-apoptotic proteins. The proteomic analysis of the MWCNT-exposed cells showed that among more than 5000 identified proteins; more than 200 were differentially expressed in the treated cells. Functional analyses revealed association of these differentially regulated proteins to cellular processes such as cell death and survival, cellular assembly, and organization. Similarly, shotgun lipidomic profiling revealed accumulation of multiple lipid classes. Our results indicate that long-term MWCNT-exposure of human normal lung cells at occupationally relevant low-doses may alter both the proteome and the lipidome profiles of the target epithelial cells in the lung.

Correction to: Epigenetic Basis of Circadian Rhythm Disruption in Cancer

The Chapter was inadvertently published with incorrect figure (Fig. 1). The same has been corrected in the chapter.

Circadian gene methylation in rotating-shift nurses: a cross-sectional study

Investigating the methylation status of the circadian genes may contribute to a better understanding of the shift work-related circadian disruption in individuals exposed to artificial light at night. In the present study, we determined the methylation status of the circadian genes associated with a shift work pattern among nurses and midwives participating in a cross-sectional study in Lodz, Poland. Quantitative methylation polymerase chain reaction assays were used to assess promoter CpG methylation in PER1, PER2, PER3, CRY1, CRY2, BMAL1, CLOCK, and NPAS2 in genomic DNA from whole blood of 347 women having a rotating-shift work schedule and 363 women working days only. The percentage of methylated reference (PMR) was assessed using fluorescent probes for PER1, PER2, PER3, CRY1, and NPAS2, and the percentage of gene methylation, as the methylation index (MI), using two sets of primers for BMAL1, CLOCK, and CRY2. We tested the possible association between current and lifetime rotating night-shift work characteristics and circadian gene methylation by using proportional odds regression model with blood DNA methylation, categorized into tertiles, and adjusted for age, current smoking status, folate intake and blood collection time. The findings indicated that CpG methylation in PER2 promoter was significantly decreased (P < 0.004) among nurses and midwives currently working rotating shifts, as compared with day-working nurses and midwives. The lower percentage of PER2 methylation was associated with a higher monthly frequency of current night duties (2-7 night shifts, and eight or more night shifts per month) (P = 0.012) and was associated at borderline significance (P = 0.092) with the lifetime duration of shift work (>10 ≤ 20 years and >20 ≤ 43 years of rotating-shift work) among nurses and midwives (N = 710). Moreover, women with a longer lifetime duration of shift work presented a lower status of PER1 methylation (P = 0.040) than did the women with up to 10 years of rotating-shift work. Long lifetime duration of shift work (> 10 years) among current rotating night-shift workers (N =  347) was associated with BMAL1 hypomethylation (P = 0.013). Among eight of the investigated circadian genes, only PER1, PER2, and BMAL1 showed differential methylation attributable to the rotating-shift work of nurses and midwives. The findings on blood-based DNA methylation in the circadian genes may provide a better insight into the mechanistic principles underlying the possible health effects of night-shift work but these should be verified in further studies recruiting larger populations of shift workers.

Sleep quality and methylation status of selected tumor suppressor genes among nurses and midwives

Chronic sleep restriction may affect metabolism, hormone secretion patterns and inflammatory responses. Limited reports suggest also epigenetic effects, such as changes in DNA methylation profiles. The study aims to assess the potential association between poor sleep quality or sleep duration and the levels of 5-methylcytosine in the promoter regions of selected tumor suppressor genes. A cross-sectional study was conducted on 710 nurses and midwives aged 40-60 years. Data from interviews regarding sleep habits and potential confounders were used. The methylation status of tumor suppressor genes was determined via qMSP reactions using DNA samples derived from leucocytes. No significant findings were observed in the total study population or in the two subgroups of women stratified by the current system of work. A borderline significance association was observed between a shorter duration of sleep and an increased methylation level in CDKN2A among day working nurses and midwives. Further studies are warranted to explore this under-investigated topic.

Sleep quality and methylation status of core circadian rhythm genes among nurses and midwives

ABSTARCT Poor sleep quality or sleep restriction is associated with sleepiness and concentration problems. Moreover, chronic sleep restriction may affect metabolism, hormone secretion patterns and inflammatory responses. Limited recent reports suggest a potential link between sleep deprivation and epigenetic effects such as changes in DNA methylation profiles. The aim of the present study was to assess the potential association between poor sleep quality or sleep duration and the levels of 5-methylcytosine in the promoter regions of PER1, PER2, PER3, BMAL1, CLOCK, CRY1 CRY2 and NPAS2 genes, taking into account rotating night work and chronotype as potential confounders or modifiers. A cross-sectional study was conducted on 710 nurses and midwives (347 working on rotating nights and 363 working only during the day) aged 40-60 years. Data from in-person interviews about sleep quality, chronotype and potential confounders were used. Sleep quality and chronotype were assessed using Pittsburgh Sleep Quality Questionnaire (PSQI) and Morningness-Eveningness Questionnaire (MEQ), respectively. Morning blood samples were collected. The methylation status of the circadian rhythm genes was determined via quantitative methylation-specific real-time PCR assays (qMSP) reactions using DNA samples derived from leucocytes. The proportional odds regression model was fitted to quantify the relationship between methylation index (MI) as the dependent variable and sleep quality or sleep duration as the explanatory variable. Analyses were carried out for the total population as well as for subgroups of women stratified by the current system of work (rotating night shift/day work) and chronotype (morning type/intermediate type/evening type). A potential modifying effect of the system of work or the chronotype was examined using the likelihood ratio test. No significant findings were observed in the total study population. Subgroup analyses revealed two statistically significant associations between a shorter sleep duration and 1) methylation level in PER2 among day workers, especially those with the morning chronotype (OR = 2.31, 95%CI:1.24-4.33), and 2) methylation level in CRY2 among subjects with the intermediate chronotype, particularly among day workers (OR = 0.52, 95%CI:0.28-0.96). The study results demonstrated a positive association between average sleep duration of less than 6 hours and the methylation level of PER2 among morning chronotype subjects, and an inverse association for CRY2 among intermediate chronotype subjects, but only among day workers. Both the system of work and the chronotype turned out to be important confounders and modifiers in a number of analyses, making it necessary to consider them as potential covariates in future research on sleep deficiency outcomes. Further studies are warranted to explore this under-investigated topic.

Cellular responses of human astrocytoma cells to dust from the Acheson process: An in vitro study

Silicon carbide (SiC) is largely used in various products such as diesel particulate filters and solar panels. It is produced through the Acheson process where aerosolized fractions of SiC and other by-products are generated in the work environment and may potentially affect the workers' health. In this study, dust was collected directly on a filter in a furnace hall over a time period of 24h. The collected dust was characterized by scanning electron microscopy and found to contain a high content of graphite particles, and carbon and silicon containing particles. Only 6% was classified as SiC, whereof only 10% had a fibrous structure. To study effects of exposure beyond the respiratory system, neurotoxic effects on human astrocytic cells, were investigated. Both low, occupationally relevant, and high doses from 9E-6μg/cm² up to 4.5μg/cm² were used, respectively. Cytotoxicity assay indicated no effects of low doses but an effect of the higher doses after 24h. Furthermore, investigation of intracellular reactive oxygen species (ROS) indicated no effects with low doses, whereas a higher dose of 0.9μg/cm² induced a significant increase in ROS and DNA damage. In summary, low doses of dust from the Acheson process may exert no or little toxic effects, at least experimentally in the laboratory on human astrocytes. However, higher doses have implications and are likely a result of the complex composition of the dust.

Mechanisms of Breast Cancer in Shift Workers: DNA Methylation in Five Core Circadian Genes in Nurses Working Night Shifts

Shift work has been suggested to be associated with breast cancer risk, and circadian disruption in shift workers is hypothesized as one of the mechanisms of increased cancer risk. There is, however, insufficient molecular evidence supporting this hypothesis. Using the quantitative methodology of pyrosequencing, epigenetic changes in 5-methyl cytosine (5mC) in five circadian genes CLOCK, BMAL1, CRY1, PER1 and PER2 in female nurses working night shift work (278 breast cancer cases, 280 controls) were analyzed. In breast cancer cases, a medium exposure to night work was associated with increased methylation levels of the CLOCK (p=0.050), BMAL1 (p=0.001) and CRY1 (p=0.040) genes, compared with controls. Within the cases, analysis of the effects of shift work on the methylation patterns showed that methylation of CRY1 was lower in those who had worked night shift and had a high exposure (p=0.006) compared with cases that had worked only days. For cases with a medium exposure to night work, an increase in BMAL1 (p=0.003) and PER1 (p=0.035) methylation was observed compared with day working (unexposed) cases. The methylation levels of the five core circadian genes were also analyzed in relation to the estrogen and progesterone receptors status of the tumors in the cases, and no correlations were observed. Furthermore, nineteen polymorphisms in the five circadian genes were assessed for their effects on the methylation levels of the respective genes, but no associations were found. In summary, our data suggest that epigenetic regulation of CLOCK, BMAL1, CRY1 and PER1 may contribute to breast cancer in shift workers.

Mechanisms of breast cancer risk in shift workers: association of telomere shortening with the duration and intensity of night work

Occupational factors such as shiftwork and especially night work that involves disruption of the circadian rhythm may contribute to increased breast cancer risk. Circadian disruption may also affect telomere length (TL). While short TL generally is associated with increased cancer risk, its association with breast cancer risk is inconclusive. We suggest that working schedules might be an important factor in assessment of effects of TL on breast cancer risk. Moreover, telomere shortening might be a potential mechanism for night work‐related breast cancer. In this study, effects of shift work on TL and its association with breast cancer risk were investigated in a nested breast cancer case–control study of Norwegian nurses. TL was assessed by qPCR in DNA from 563 breast cancer patients and 619 controls. Here, we demonstrate that TL is affected by intensive night work schedules, as work with six consecutive night for a period of more than 5 years was associated with decreased telomere lengths (−3.18, 95% CI: −6.46 to −0.58, P = 0.016). Furthermore, telomere shortening is associated with increased breast cancer risk in workers with long periods of consecutive night shifts. Thus, nurses with longer telomere lengths had a lower risk for breast cancer if they had worked more than four (OR: 0.37, 95% CI: 0.16–0.79, P = 0.014) or five (OR: 0.31, 95% CI: 0.10–0.83, P = 0.029) consecutive night shifts for a period of 5 years or more. These data suggest that telomere shortening is associated with the duration and intensity of night work and may be a contributing factor for breast cancer risk among female shift workers.

Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes

Although several lung cancer susceptibility loci have been identified, much of the heritability for lung cancer remains unexplained. Here 14,803 cases and 12,262 controls of European descent were genotyped on the OncoArray and combined with existing data for an aggregated genome-wide association study (GWAS) analysis of lung cancer in 29,266 cases and 56,450 controls. We identified 18 susceptibility loci achieving genome-wide significance, including 10 new loci. The new loci highlight the striking heterogeneity in genetic susceptibility across the histological subtypes of lung cancer, with four loci associated with lung cancer overall and six loci associated with lung adenocarcinoma. Gene expression quantitative trait locus (eQTL) analysis in 1,425 normal lung tissue samples highlights RNASET2, SECISBP2L and NRG1 as candidate genes. Other loci include genes such as a cholinergic nicotinic receptor, CHRNA2, and the telomere-related genes OFBC1 and RTEL1. Further exploration of the target genes will continue to provide new insights into the etiology of lung cancer.

Copy number variation, increased gene expression, and molecular mechanisms of neurofascin in lung cancer

Metastasis and cell adhesion are key aspects of cancer progression. Neurofascin (NFASC) is a member of the immunoglobulin superfamily of adhesion molecules and, while studies on NFASC are inadequate, other members have been indicated pivotal roles in cancer progression and metastasis. This study aimed at increasing the knowledge on the involvement of adhesion molecules in lung cancer progression by studying the regulation and role of NFASC in non‐small cell lung cancer (NSCLC). Here, copy number variations in the NFASC gene were analyzed in tumor and non‐tumorous lung tissues of 204 NSCLC patients. Frequent gene amplifications (OR = 4.50, 95%CI: 2.27‐8.92, P ≤ 0.001) and increased expression of NFASC (P = 0.034) were identified in tumors of NSCLC patients. Furthermore, molecular mechanisms of NFASC in lung cancer progression were evaluated by investigating the effects of NFASC silencing on cell proliferation, viability, migration, and invasion using siRNA technology in four NSCLC cell lines. Silencing of NFASC did not affect cell proliferation or viability but rather decreased NSCLC cell migration (P ≤ 0.001) and led to morphological changes, rearrangements in the actin cytoskeleton and changes in F‐actin networks in migrating NSCLC cell lines. This study is the first to report frequent copy number gain and increased expression of NFASC in NSCLC. Moreover, these data suggest that NFASC is a novel regulator of NSCLC cell motility and support a role of NFASC in the regulation of NSCLC progression.

In vitro genotoxicity testing of four reference metal nanomaterials, titanium dioxide, zinc oxide, cerium oxide and silver: towards reliable hazard assessment

There is serious concern about the potential harmful effects of certain nanomaterials (NMs), on account of their ability to penetrate cell membranes and the increased reactivity that results from their increased surface area compared with bulk chemicals. To assess the safety of NMs, reliable tests are needed. We have investigated the possible genotoxicity of four representative NMs, derived from titanium dioxide, zinc oxide, cerium oxide and silver, in two human cell lines, A549 alveolar epithelial cells and lymphoblastoid TK6 cells. A high-throughput version of the comet assay was used to measure DNA strand beaks (SBs) as well as oxidised purines (converted to breaks with the enzyme formamidopyrimidine DNA glycosylase). In parallel, cytotoxicity was measured with the alamarBlue® assay, and the ability of NM-treated cells to survive was assessed by their colony-forming efficiency. TiO₂ and CeO₂ NMs were only slightly cytotoxic by the alamarBlue® test, and had no long-term effect on colony-forming efficiency. However, both induced DNA damage at non-cytotoxic concentrations; the damage decreased from 3 to 24-h exposure, except in the case of CeO₂-treated A549 cells. ZnO and Ag NMs affected cell survival, and induced high levels of DNA damage at cytotoxic concentrations. At lower concentrations, there was significant damage, which tended to persist over 24 h. The implication is that all four reference metal NMs tested-whether cytotoxic or not-are genotoxic. A full assessment of NM toxicity should include tests on different cell types, different times of incubation and a wide range of (especially non-cytotoxic) concentrations; a test for cell viability should be performed in parallel. Inclusion of Fpg in the comet assay allows detection of indirect genotoxic effects via oxidative stress.

Night shift work and other determinants of estradiol, testosterone, and dehydroepiandrosterone sulfate among middle-aged nurses and midwives

OBJECTIVES

The aims of our study were to (i) investigate the association between rotating night shift work and blood concentrations of estradiol, testosterone and dehydroepiandrosterone sulfate (DHEAS) and (2) evaluate the role of their non-occupational determinants.

METHODS

A cross-sectional study was conducted on 345 premenopausal and 187 postmenopausal nurses and midwives (263 women working rotating night shifts and 269 women working during days). Data from in-person interviews were used, anthropometric measurements were performed, and body mass index (BMI) and waist- to-hip ratio were calculated. Morning blood and spot urine samples were collected. Multiple linear regression models were fitted with hormone concentrations as dependent variables, and night shift work characteristics and demographic, reproductive, lifestyle and anthropometric determinants as independent variables. Modification of the effect by chronotype was examined.

RESULTS

Among postmenopausal women, we observed a statistically significant positive association between the total duration of night shift work >15 years and estradiol level (P<0.05 when compared to night work duration <5 years). Night shift work characteristics were significantly associated with estradiol among morning-type postmenopausal women. The well-established associations between hormones and their major determinants, such as age and BMI, were confirmed.

CONCLUSIONS

The findings of our study imply that prolonged night shift work may be associated with increased estradiol levels among postmenopausal women, especially among the morning-type postmenopausal women.

A rare missense mutation in CHRNA4 associates with smoking behavior and its consequences

Using Icelandic whole-genome sequence data and an imputation approach we searched for rare sequence variants in CHRNA4 and tested them for association with nicotine dependence. We show that carriers of a rare missense variant (allele frequency=0.24%) within CHRNA4, encoding an R336C substitution, have greater risk of nicotine addiction than non-carriers as assessed by the Fagerstrom Test for Nicotine Dependence (P=1.2 × 10(-4)). The variant also confers risk of several serious smoking-related diseases previously shown to be associated with the D398N substitution in CHRNA5. We observed odds ratios (ORs) of 1.7-2.3 for lung cancer (LC; P=4.0 × 10(-4)), chronic obstructive pulmonary disease (COPD; P=9.3 × 10(-4)), peripheral artery disease (PAD; P=0.090) and abdominal aortic aneurysms (AAAs; P=0.12), and the variant associates strongly with the early-onset forms of LC (OR=4.49, P=2.2 × 10(-4)), COPD (OR=3.22, P=2.9 × 10(-4)), PAD (OR=3.47, P=9.2 × 10(-3)) and AAA (OR=6.44, P=6.3 × 10(-3)). Joint analysis of the four smoking-related diseases reveals significant association (P=6.8 × 10(-5)), particularly for early-onset cases (P=2.1 × 10(-7)). Our results are in agreement with functional studies showing that the human α4β2 isoform of the channel containing R336C has less sensitivity for its agonists than the wild-type form following nicotine incubation.

Effects of carbon nanotubes on intercellular communication and involvement of IL-1 genes

An increasing amount of products containing engineered nanoparticles is emerging. Among these particles are carbon nanotubes (CNTs) which are of interest for a wide range of industrial and biomedical applications. There have been raised concerns over the effects of CNTs on human health. Some types of CNTs are classified as group 2B carcinogens by the International Agency for Research on Cancer. CNTs may also induce pulmonary inflammatory and fibrotic effects. By utilizing CNTs of different lengths, we investigated the role of the proinflammatory cytokine, interleukin-1 (IL-1) on gap junctional intercellular communication (GJIC) by using IL-1 wild-type (IL1-WT) and IL-1 knock-out (IL1-KO) cells. GJIC decreased equally in both cell types after CNT exposure. Immunofluorescence staining showed Gja1 and Gjb2 in gap junctions and hemichannels for both cell types. Gjb1 and Gjb2 expression was low in IL1-KO cells, which was confirmed by protein analysis. Gja1 was upregulated with both CNTs, whereas Gjb1 was down-regulated by CNT-2 in IL1-WT cells. Connexin mRNA expression was regulated differently by the CNTs. CNT-1 affected Gja1 and Gjb2, whereas CNT-2 had an effect on Gjb1. CNTs negatively affect GJIC through gap junctions independently of the length of CNT and IL-1 status. Furthermore, connexin gene expression was affected by IL-1 at transcriptional and translational levels. As both CNTs used in this study are cytotoxic to the cells and reduce cell survival, we suggest that CNT-induced reduction in GJIC may be important for inhibiting transfer of cell survival signals between cells.

Loss of MKK3 and MK2 Copy Numbers in Non-Small Cell Lung Cancer

Identification of genetic alterations in members of the p38 mitogen-activated protein kinase (MAPK) pathway is important as these proteins have dynamic roles in tumor progression and may serve as potential therapeutic targets in cancer. We analyzed tumor and non-tumorous lung tissue of 233 non-small cell lung cancer (NSCLC) patients for the presence of copy number alterations (CNAs) in the MAPK kinase 3 (MKK3) and MAPK-activated kinase 2 (MK2) genes. We report frequent CNAs in MKK3 and MK2 genes in NSCLC. Copy number losses were detected in 31% of NSCLC tumors (odds ratio: 7.08, 95% confidence interval: 3.2-15.6, P<0.001) for the MKK3 gene and in 28% of tumors for the MK2 gene (odds ratio: 3.68, 95% confidence interval: 1.9-7.2, P<0.001). Several of the non-tumorous tissues showed an elevated MKK3 copy number, with a concurrent loss of this in 89% of the paired tumors. MKK3 gene deletions were significantly more frequent in squamous and large cell carcinoma than in adenocarcinoma. These data demonstrate a novel loss of MKK3 and MK2 genomic copy numbers in NSCLC tumors, and suggest these genes as interesting therapeutic candidates in NSCLC.

Abstract 1214: TP53 affects SOX2 copy number alterations and expression in non-small cell lung cancer

Background: Amplifications of the transcription factor, SRY (sex determining region Y)-box 2 (SOX2), are common in non-small cell lung cancer (NSCLC). SOX2-signaling is important in maintaining the stem cell-like phenotype of cancer cells and contributes to the pathogenesis of lung cancer. TP53 is recognized as a critical regulator of stem cell pluripotency, and it is known that TP53 represses many stem cell associated genes following DNA damage. We hypothesized that SOX2 copy number alterations in lung tumors could be correlated to mutational status of TP53 gene in tumors and that TP53 played a role in regulation of SOX2 gene expression. Material and Methods: Early-stage lung cancer cases (n = 258) of Norwegian Caucasians origin were recruited when admitted to Haukeland University Hospital in Bergen, Norway between 1986 and 1994. Samples of adjacent histologically non-tumorous lung tissue were collected in the lobectomi specimens at the time of surgery. Tumor histology was confirmed by an experienced pathologist and samples containing &gt;80% of tumor cells were analyzed in this study. After resection tumor and non-tumorous tissues were snap-frozen in liquid nitrogen and kept at -80oC until further processing. All subjects gave written consent, and the study was approved by the regional ethical committee in accordance with the Helsinki Declaration. SOX2 copy number alterations were evaluated by quantitative RT- PCR using SYBR Green I Technology. The TP53 and SOX2 expression levels were assessed by q-PCR in the tumors where mRNA was available. Furthermore, to understand mechanisms, TP53 and SOX2 genes were knocked down by siRNA in the A427 human lung adenocarcinoma cell line. Additionally, some miRNAs were also analyzed by q-PCR. TP53 mutations in tumors were analyzed by sequencing and PAH-DNA adduct levels were determined in non-tumorous lung tissue by 382P-postlabelling.

Results: SOX2 copy number alterations were more frequent in tumor tissues (34%) than in the adjacent non-tumorous tissues (3%) and TP53 mutations were associated with an increased risk of acquiring a SOX2 copy number alteration (OR = 2.08, 95% CI: 1.14-3.79, p = 0.017). SOX2 and TP53 expressions were strongly correlated in lung tumors and reduction in TP53 decreased SOX2 expression in lung adenocarcinoma cells. TP53 knockdown also reduced the miRNA hsa-miR-145, which has previously been shown to regulate SOX2 expression. We propose that SOX2 may be an important target of TP53 signaling in lung cancer and that hsa-miR-145-5p may be one potential driver miRNA involved in this regulation requiring further studies on the mechanisms behind the TP53-induced regulation of SOX2 expression and the possible importance of hsa-miR-145 in lung cancer.

Citation Format: Johanna Samulin-Erdem, Vidar Skaug, Per Bakke, Amund Gulsvik, Aage Haugen, Shanbeh Zienolddiny. TP53 affects SOX2 copy number alterations and expression in non-small cell lung cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1214. doi:10.1158/1538-7445.AM2015-1214

CHRNA5 risk variant predicts delayed smoking cessation and earlier lung cancer diagnosis--a meta-analysis

BACKGROUND

Recent meta-analyses show strong evidence of associations among genetic variants in CHRNA5 on chromosome 15q25, smoking quantity, and lung cancer. This meta-analysis tests whether the CHRNA5 variant rs16969968 predicts age of smoking cessation and age of lung cancer diagnosis.

METHODS

Meta-analyses examined associations between rs16969968, age of quitting smoking, and age of lung cancer diagnosis in 24 studies of European ancestry (n = 29 072). In each dataset, we used Cox regression models to evaluate the association between rs16969968 and the two primary phenotypes (age of smoking cessation among ever smokers and age of lung cancer diagnosis among lung cancer case patients) and the secondary phenotype of smoking duration. Heterogeneity across studies was assessed with the Cochran Q test. All statistical tests were two-sided.

RESULTS

The rs16969968 allele (A) was associated with a lower likelihood of smoking cessation (hazard ratio [HR] = 0.95, 95% confidence interval [CI] = 0.91 to 0.98, P = .0042), and the AA genotype was associated with a four-year delay in median age of quitting compared with the GG genotype. Among smokers with lung cancer diagnoses, the rs16969968 genotype (AA) was associated with a four-year earlier median age of diagnosis compared with the low-risk genotype (GG) (HR = 1.08, 95% CI = 1.04 to 1.12, P = 1.1*10(-5)).

CONCLUSION

These data support the clinical significance of the CHRNA5 variant rs16969968. It predicts delayed smoking cessation and an earlier age of lung cancer diagnosis in this meta-analysis. Given the existing evidence that this CHRNA5 variant predicts favorable response to cessation pharmacotherapy, these findings underscore the potential clinical and public health importance of rs16969968 in CHRNA5 in relation to smoking cessation success and lung cancer risk.

Association of the FAM46A gene VNTRs and BAG6 rs3117582 SNP with non small cell lung cancer (NSCLC) in Croatian and Norwegian populations

We analyzed for associations between a variable number of tandem repeat (VNTR) polymorphism in the Family with sequence similarity 46, member A (FAM46A) gene and a single nucleotide polymorphism (rs3117582) in the BCL2-Associated Athanogene 6 (BAG6) with non small cell lung cancer in Croatian and Norwegian subjects. A total of 503 (262 Croatian and 241Norwegian) non small cell lung cancer patients and 897 controls (568 Croatian and 329 Norwegian) were analyzed. We found that the frequency of allele b (three VNTR repeats) of FAM46A gene was significantly increased in the patients compared to the healthy controls in the Croatian and the combined Croatian and Norwegian subjects. Genotype frequencies of cd (four and five VNTR repeats) and cc (four VNTR repeats homozygote) of the FAM46A gene were significantly decreased in the patients compared to the healthy controls in the Croatian and Norwegian subjects, respectively. Logistic regression analyses revealed FAM46A genotype cc to be an independent predictive factor for non small cell lung cancer risk in the Norwegian subjects after adjustment for age, gender and smoking status. This is the first study to suggest an association between the FAM46A gene VNTR polymorphisms and non small cell lung cancer. We found also that BAG6 rs3117582 SNP was associated with non small cell lung cancer in the Norwegian subjects and the combined Croatian-Norwegian subjects corroborating the earlier finding that BAG6 rs3117582 SNP was associated with lung cancer in Europeans. Logistic regression analyses revealed that genotypes and alleles of BAG6 were independent predictive factor for non small cell lung cancer risk in the Norwegian and combined Croatian-Norwegian subjects, after adjustment for age and gender.

Analysis of polymorphisms in the circadian-related genes and breast cancer risk in Norwegian nurses working night shifts

Introduction

Some studies have suggested that night work may be associated with an increased risk of breast cancer in nurses. We aimed to explore the role of circadian gene polymorphisms in the susceptibility to night work-related breast cancer risk.

Methods

We conducted a nested case-control study of Norwegian nurses comprising 563 breast cancer cases and 619 controls within a cohort of 49,402 Norwegian nurses ages 35 to 74 years. We studied 60 single-nucleotide polymorphisms (SNPs) in 17 genes involved in the regulation of the circadian rhythm in cases and controls. The data were analyzed in relation to the two exposure variables "maximum number of consecutive night shifts ever worked" and "maximum number of consecutive night shifts worked for at least 5 years." The odds of breast cancer associated with each SNP was calculated in the main effects analysis and in relation to night shift work. The statistically significant odds ratios were tested for noteworthiness using two Bayesian tests: false positive report probability (FPRP) and Bayesian false discovery probability (BFDP).

Results

In the main effects analysis, CC carriers of rs4238989 and GG carriers of rs3760138 in the AANAT gene had increased risk of breast cancer, whereas TT carriers of BMAL1 rs2278749 and TT carriers of CLOCK rs3749474 had reduced risk. The associations were found to be noteworthy using both the FPRP and BFDP tests. With regard to the effect of polymorphisms and night work, several significant associations were observed. After applying FPRP and BFDP in women with at least four night shifts, an increased risk of breast cancer was associated with variant alleles of SNPs in the genes AANAT (rs3760138, rs4238989), BMAL1 (rs2290035, rs2278749, rs969485) and ROR-b (rs3750420). In women with three consecutive night shifts, a reduced risk of breast cancer was associated with carriage of variant alleles of SNPs in CLOCK (rs3749474), BMAL1 (rs2278749), BMAL2 (rs2306074), CSNK1E (rs5757037), NPAS2 (rs17024926), ROR-b (rs3903529, rs3750420), MTNR1A (rs131113549) and PER3 (rs1012477).

Conclusions

Significant and noteworthy associations between several polymorphisms in circadian genes, night work and breast cancer risk were found among nurses who had worked at least three consecutive night shifts.