Search

The tumor microenvironment (TME), including immune cells, cancer‐associated fibroblasts, endothelial cells, adjacent normal cells, and others, plays a crucial role in influencing tumor behavior and progression. Here, we characterized the TME in 83 primary renal tumors and matched metastatic or recurrence tissue samples (n = 15) from papillary renal cell carcinoma (pRCC) types 1 (n = 20) and 2 (n = 49), collecting duct carcinomas (CDC; n = 14), and high‐grade urothelial carcinomas (HGUC; n = 5). We investigated 10 different markers of immune infiltration, vasculature, cell proliferation, and epithelial‐to‐mesenchymal transition by using machine learning image analysis in conjunction with immunohistochemistry. Marker expression was compared by Mann–Whitney and Kruskal–Wallis tests and correlations across markers using Spearman's rank correlation coefficient. Multivariable Poisson regression analysis was used to compare marker expression between histological types, while accounting for variation in tissue size. Several immune markers showed different rates of expression across histological types of renal carcinoma. Using pRCC1 as reference, the incidence rate ratio (IRR) of CD3+ T cells (IRR [95% confidence interval, CI] = 2.48 [1.53–4.01]) and CD20+ B cells (IRR [95% CI] = 4.38 [1.22–5.58]) was statistically significantly higher in CDC. In contrast, CD68+ macrophages predominated in pRCC1 (IRR [95% CI] = 2.35 [1.42–3.9]). Spatial analysis revealed CD3+ T‐cell and CD20+ B‐cell expressions in CDC to be higher at the proximal (p

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations, we identified a novel NRAS variant (c.170A>C, p.D57A) in a melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11 273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2+ binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A>C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted.

Although genome-wide association studies have identified more than eighty genetic variants associated with non-small cell lung cancer (NSCLC) risk, biological mechanisms of these variants remain largely unknown. By integrating a large-scale genotype data of 15 581 lung adenocarcinoma (AD) cases, 8350 squamous cell carcinoma (SqCC) cases, and 27 355 controls, as well as multiple transcriptome and epigenomic databases, we conducted histology-specific meta-analyses and functional annotations of both reported and novel susceptibility variants. We identified 3064 credible risk variants for NSCLC, which were overrepresented in enhancer-like and promoter-like histone modification peaks as well as DNase I hypersensitive sites. Transcription factor enrichment analysis revealed that USF1 was AD-specific while CREB1 was SqCC-specific. Functional annotation and gene-based analysis implicated 894 target genes, including 274 specifics for AD and 123 for SqCC, which were overrepresented in somatic driver genes (ER = 1.95, P = 0.005). Pathway enrichment analysis and Gene-Set Enrichment Analysis revealed that AD genes were primarily involved in immune-related pathways, while SqCC genes were homologous recombination deficiency related. Our results illustrate the molecular basis of both well-studied and new susceptibility loci of NSCLC, providing not only novel insights into the genetic heterogeneity between AD and SqCC but also a set of plausible gene targets for post-GWAS functional experiments.

Epidemiologic studies often rely on questionnaire data, exposure measurement tools, and/or biomarkers to identify risk factors and the underlying carcinogenic processes. An emerging and promising complementary approach to investigate cancer etiology is the study of somatic “mutational signatures” that endogenous and exogenous processes imprint on the cellular genome. These signatures can be identified from a complex web of somatic mutations thanks to advances in DNA sequencing technology and analytical algorithms. This approach is at the core of the Sherlock-Lung study (2018–ongoing), a retrospective case-only study of over 2,000 lung cancers in never-smokers (LCINS), using different patterns of mutations observed within LCINS tumors to trace back possible exposures or endogenous processes. Whole genome and transcriptome sequencing, genome-wide methylation, microbiome, and other analyses are integrated with data from histological and radiological imaging, lifestyle, demographic characteristics, environmental and occupational exposures, and medical records to classify LCINS into subtypes that could reveal distinct risk factors. To date, we have received samples and data from 1,370 LCINS cases from 17 study sites worldwide and whole-genome sequencing has been completed on 1,257 samples. Here, we present the Sherlock-Lung study design and analytical strategy, also illustrating some empirical challenges and the potential for this approach in future epidemiologic studies.

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene–sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency

BACKGROUND: CDKN2A, CDK4, and POT1 are well-established melanoma-susceptibility genes. OBJECTIVE: We evaluated melanoma histopathology for individuals with germline mutations of CDKN2A, CDK4, and POT1. METHODS: We assessed histopathology for melanomas diagnosed in melanoma-prone families (≥2 individuals with melanoma) from the United States, Italy, and Spain. Comparisons between mutation carriers and non-carriers (no mutation) were adjusted for age, sex, Breslow depth, and correlations among individuals within the same family. RESULTS: Histologic slides were evaluated for 290 melanomas (139 from 132 non-carriers, 122 from 68 CDKN2A carriers, 10 from 6 CDK4 carriers, 19 from 16 POT1 carriers). Superficial spreading was the predominant subtype for all groups. Spitzoid morphology (>25% of tumor) was observed in the majority of invasive melanomas from POT1 carriers (10 of 15 [67%], P

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.

Intratumor heterogeneity (ITH) of genomic alterations may impact prognosis of lung adenocarcinoma (LUAD). Here, we investigate ITH of somatic copy number alterations (SCNAs), DNA methylation, and point mutations in lung cancer driver genes in 292 tumor samples from 84 patients with LUAD. LUAD samples show substantial SCNA and methylation ITH, and clonal architecture analyses present congruent evolutionary trajectories for SCNAs and DNA methylation aberrations. Methylation ITH mapping to gene promoter areas or tumor suppressor genes is low. Moreover, ITH composed of genetic and epigenetic mechanisms altering the same cancer driver genes is shown in several tumors. To quantify ITH for valid statistical association analyses, we develope an average pairwise ITH index (APITH), which does not depend on the number of samples per tumor. Both APITH indexes for SCNAs and methylation aberrations show significant associations with poor prognosis. This study further establishes the important clinical implications of genetic and epigenetic ITH in LUAD., Many tumors are known to be heterogeneous. Here, the authors examined multiple samples from 84 patients with lung adenocarcinoma and demonstrate that the intratumor heterogeneity of methylation and copy number associates with poor prognosis.

Introduction: Aberrant Wnt signalling, regulating cell development and stemness, is observed in many cancer entities. Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated to 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 genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2 and other lung cancer-related genes. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer. Further, decision trees were created to investigate variant x variant interaction. All analyses were performed for overall lung cancer and for subgroups. Results: No association 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. Although predictability is poor, AhR/Wnt-variants are unexpected overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkable, 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, no assigned to any CHRN gene. 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.

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

Epidemiological studies have observed that risk of endometriosis is associated with history of cutaneous melanoma and vice versa. Evidence for shared biological mechanisms between the two traits is limited. The aim of this study was to investigate the genetic correlation and causal relationship between endometriosis and melanoma. Summary statistics from genome-wide association meta-analyses (GWAS) for endometriosis and melanoma were used to estimate the genetic correlation between the traits and Mendelian randomization was used to test for a causal association. When using summary statistics from separate female and male melanoma cohorts we identified a significant positive genetic correlation between melanoma in females and endometriosis (rg = 0.144, se = 0.065, p = 0.025). However, we find no evidence of a correlation between endometriosis and melanoma in males or a combined melanoma dataset. Endometriosis was not genetically correlated with skin color, red hair, childhood sunburn occasions, ease of skin tanning, or nevus count suggesting that the correlation between endometriosis and melanoma in females is unlikely to be influenced by pigmentary traits. Mendelian Randomization analyses also provided evidence for a relationship between the genetic risk of melanoma in females and endometriosis. Colocalization analysis identified 27 genomic loci jointly associated with the two diseases regions that contain different causal variants influencing each trait independently. This study provides evidence of a small genetic correlation and relationship between the genetic risk of melanoma in females and endometriosis. Genetic risk does not equate to disease occurrence and differences in the pathogenesis and age of onset of both diseases means it is unlikely that occurrence of melanoma causes endometriosis. This study instead provides evidence that having an increased genetic risk for melanoma in females is related to increased risk of endometriosis. Larger GWAS studies with increased power will be required to further investigate these associations.

ObjectivesAsbestos is the main risk factor for peritoneal mesothelioma (PeM). However, due to its rarity, PeM has rarely been investigated in community-based studies. We examined the association between asbestos exposure and PeM risk in a general population in Lombardy, Italy.MethodsFrom the regional mesothelioma registry, we selected PeM cases diagnosed in 2000–2015. Population controls (matched by area, gender and age) came from two case–control studies in Lombardy on lung cancer (2002–2004) and pleural mesothelioma (2014). Assessment of exposure to asbestos was performed through a quantitative job-exposure matrix (SYN-JEM) and expert evaluation based on a standardised questionnaire. We calculated period-specific and gender-specific OR and 90% CI using conditional logistic regression adjusted for age, province of residence and education.ResultsWe selected 68 cases and 2116 controls (2000–2007) and 159 cases and 205 controls (2008–2015). The ORs for ever asbestos exposure (expert-based, 2008–2015 only) were 5.78 (90% CI 3.03 to 11.0) in men and 8.00 (2.56 to 25.0) in women; the ORs for definite occupational exposure were 12.3 (5.62 to 26.7) in men and 14.3 (3.16 to 65.0) in women. The ORs for ever versus never occupational asbestos exposure based on SYN-JEM (both periods) were 2.05 (90% CI 1.39 to 3.01) in men and 1.62 (0.79 to 3.27) in women. In men, clear positive associations were found for duration, cumulative exposure (OR 1.33 (1.19 to 1.48) per fibres/mL-years) and latency.ConclusionsUsing two different methods of exposure assessment we provided evidence of a clear association between asbestos exposure and PeM risk in the general population.

L.K. is a fellow in the Canadian Institutes of Health Research (CIHR) Strategic Training in Advanced Genetic Epidemiology (STAGE) programme and is supported by the CIHR Doctoral Research Award from the Frederick Banting and Charles Best Canada Graduate Scholarships (GSD-137441). Transdisciplinary Research for Cancer in Lung (TRICL) of the International Lung Cancer Consortium (ILCCO) was supported by the National Institutes of Health (U19-CA148127, CA148127S1). Genotyping for the TRICL-ILCCO OncoArray was supported by in-kind genotyping at Centre for Inherited Disease Research (CIDR) (26820120008i-0–6800068-1). Genotyping for the Head and Neck Cancer OncoArray performed at CIDR was funded by the US National Institute of Dental and Craniofacial Research (NIDCR) grant 1X01HG007780–0. CAPUA study was supported by FIS-FEDER/Spain grant numbers FIS-01/310, FIS-PI03–0365 and FIS-07-BI060604, FICYT/Asturias grant numbers FICYT PB02–67 and FICYT IB09–133, and the University Institute of Oncology (IUOPA), of the University of Oviedo and the Ciber de Epidemiologia y Salud Publica. CIBERESP, SPAIN. The work performed in the CARET study was supported by the National Institute of Health (NIH)/National Cancer Institute (NCI): UM1 CA167462 (PI: Goodman), National Institute of Health UO1-CA6367307 (PIs Omen, Goodman); National Institute of Health R01 CA111703 (PI Chen), National Institute of Health 5R01 CA151989 (PI Doherty). The Liverpool Lung Project is supported by the Roy Castle Lung Cancer Foundation. The Harvard Lung Cancer Study was supported by the NIH (National Cancer Institute) grants CA092824, CA090578 and CA074386. The Multiethnic Cohort Study was partially supported by NIH Grants CA164973, CA033619, CA63464 and CA148127. The work performed in MSH-PMH study was supported by the Canadian Cancer Society Research Institute (020214), Ontario Institute of Cancer and Cancer Care Ontario Chair Award to R.J.H. and G.L. and the Alan Brown Chair and Lusi Wong Programs at the Princess Margaret Hospital Foundation. The Norway study was supported by Norwegian Cancer Society, Norwegian Research Council. The work in TLC study has been supported in part the James & Esther King Biomedical Research Program (09KN-15), National Institutes of Health Specialized Programs of Research Excellence (SPORE) Grant (P50 CA119997) and by a Cancer Center Support Grant (CCSG) at the H. Lee Moffitt Cancer Center and Research Institute, an NCI designated Comprehensive Cancer Center (grant number P30-CA76292). The dataset(s) used for the analyses described were obtained from Vanderbilt University Medical Center’s BioVU, which is supported by institutional funding and by the Vanderbilt CTSA grant UL1 TR000445 from NCATS/NIH. Dr Melinda Aldrich is supported by the by NIH/National Cancer Institute 5K07CA172294. The Copenhagen General Population Study (CGPS) was supported by the Chief Physician Johan Boserup and Lise Boserup Fund, the Danish Medical Research Council and Herlev Hospital. The NELCS study: Grant Number P20RR018787 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH). Kentucky Lung Cancer Research Initiative (KLCRI) was supported by the Department of Defense (Congressionally Directed Medical Research Program, U.S. Army Medical Research and Materiel Command Program) under award number: 10153006 (W81XWH-11–1-0781). Views and opinions of, and endorsements by the author(s) do not reflect those of the US Army or the Department of Defense. This research was also supported by unrestricted infrastructure funds from the UK Center for Clinical and Translational Science, NIH grant UL1TR000117 and Markey Cancer Center NCI Cancer Center Support Grant (P30 CA177558) Shared Resource Facilities: Cancer Research Informatics, Biospecimen and Tissue Procurement, and Biostatistics and Bioinformatics. The research undertaken by M.D.T., L.V.W. and M.S.A. was partly funded by the National Institute for Health Research (NIHR). The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. M.D.T. holds a Medical Research Council Senior Clinical Fellowship (G0902313). The Tampa study was funded by Public Health Service grants P01-CA68384 and R01-DE13158 from the National Institutes of Health. The University of Pittsburgh head and neck cancer case–control study is supported by US National Institutes of Health grants P50 CA097190 and P30 CA047904. The Carolina Head and Neck Cancer Study (CHANCE) was supported by the National Cancer Institute (R01CA90731). The Head and Neck Genome Project (GENCAPO) was supported by the Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP; grants 04/12054–9 and 10/51168–0). The authors thank all the members of the GENCAPO team. This publication presents data from the Head and Neck 5000 study. The study was a component of independent research funded by the National Institute for Health Research (NIHR) under its Programme Grants for Applied Research scheme (RP-PG-0707–10034). The views expressed in this publication are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. Human papillomavirus (HPV) serology was supported by a Cancer Research UK Programme Grant, the Integrative Cancer Epidemiology Programme (grant number: C18281/A19169). The Alcohol-Related Cancers and Genetic Susceptibility Study in Europe (ARCAGE) was funded by the European Commission’s fifth framework programme (QLK1– 2001-00182), the Italian Association for Cancer Research, Compagnia di San Paolo/FIRMS, Region Piemonte and Padova University (CPDA057222). The Rome Study was supported by the Associazione Italiana per la Ricerca sul Cancro (AIRC) awards IG 2011 10491 and IG 2013 14220 to S.B. and by Fondazione Veronesi to S.B. The IARC Latin American study was funded by the European Commission INCO-DC programme (IC18-CT97–0222), with additional funding from Fondo para la Investigacion Cientifica y Tecnologica (Argentina) and the Fundacao de Amparo a Pesquisa do Estado de Sao Paulo (01/01768–2). The IARC Central Europe study was supported by the European Commission’s INCO-COPERNICUS Program (IC15-CT98–0332), US NIH/National Cancer Institute grant CA92039 and World Cancer Research Foundation grant WCRF 99A28. The IARC Oral Cancer Multicenter study was funded by grant S06 96 202489 05F02 from Europe against Cancer; grants FIS 97/0024, FIS 97/0662 and BAE 01/5013 from Fondo de Investigaciones Sanitarias, Spain; the UICC Yamagiwa-Yoshida Memorial International Cancer Study; the National Cancer Institute of Canada; Associazione Italiana per la Ricerca sul Cancro; and the Pan-American Health Organization. Coordination of the EPIC study is financially supported by the European Commission (DG SANCO) and the International Agency for Research on Cancer.

Background: Germline mutations in telomere-related genes such as POT1 and TERT predispose individuals to familial melanoma. Objectives: To evaluate the prevalence of germline mutations in POT1 and TERT in a large cohort of Spanish melanoma-prone families (at least two affected first- or second-degree relatives). Methods: Overall, 228 CDKN2A wild-type melanoma-prone families were included in the study. Screening of POT1 was performed in one affected person from each family and TERT was sequenced in one affected patient from 202 families (26 families were excluded owing to DNA exhaustion/degradation). TERT promoter sequencing was extended to an additional 30 families with CDKN2A mutation and 70 patients with sporadic multiple primary melanoma (MPM) with a family history of other cancers. Results: We identified four families with potentially pathogenic POT1 germline mutations: a missense variant c.233T>C (p.Ile78Thr); a nonsense variant c.1030G>T (p.Glu344*); and two other variants, c.255G>A (r.125_255del) and c.1792G>A (r.1791_1792insAGTA, p.Asp598Serfs*22), which we confirmed disrupted POT1 mRNA splicing. A TERT promoter variant of unknown significance (c.-125C>A) was detected in a patient with MPM, but no germline mutations were detected in TERT promoter in cases of familial melanoma. Conclusions: Overall, 1·7% of our CDKN2A/CDK4-wild type Spanish melanoma-prone families carry probably damaging mutations in POT1. The frequency of TERT promoter germline mutations in families with melanoma in our population is extremely rare.

Genome-wide association studies have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP, and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). As ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWAS identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor gene (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. As AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus, and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression, and altering melanocyte growth phenotypes upon exposure.

While expression quantitative trait loci (eQTL) have been powerful in identifying susceptibility genes from genome-wide association studies (GWAS) findings, most trait-associated loci are not explained by eQTL alone. Alternative QTLs including DNA methylation QTL (meQTL) are emerging, but cell-type-specific meQTL using cells of disease origin has been lacking. Here we established an meQTL dataset using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization using meQTL, eQTL, and mRNA splice-junction QTL from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified susceptibility genes at 63% of melanoma GWAS loci. Among three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes are preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTL and melanoma meQTL identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTL identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIPseq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTL.

Although expression quantitative trait loci (eQTLs) have been powerful in identifying susceptibility genes from genome-wide association study (GWAS) findings, most trait-associated loci are not explained by eQTLs alone. Alternative QTLs, including DNA methylation QTLs (meQTLs), are emerging, but cell-type-specific meQTLs using cells of disease origin have been lacking. Here, we established an meQTL dataset by using primary melanocytes from 106 individuals and identified 1,497,502 significant cis-meQTLs. Multi-QTL colocalization with meQTLs, eQTLs, and mRNA splice-junction QTLs from the same individuals together with imputed methylome-wide and transcriptome-wide association studies identified candidate susceptibility genes at 63% of melanoma GWAS loci. Among the three molecular QTLs, meQTLs were the single largest contributor. To compare melanocyte meQTLs with those from malignant melanomas, we performed meQTL analysis on skin cutaneous melanomas from The Cancer Genome Atlas (n = 444). A substantial proportion of meQTL probes (45.9%) in primary melanocytes is preserved in melanomas, while a smaller fraction of eQTL genes is preserved (12.7%). Integration of melanocyte multi-QTLs and melanoma meQTLs identified candidate susceptibility genes at 72% of melanoma GWAS loci. Beyond GWAS annotation, meQTL-eQTL colocalization in melanocytes suggested that 841 unique genes potentially share a causal variant with a nearby methylation probe in melanocytes. Finally, melanocyte trans-meQTLs identified a hotspot for rs12203592, a cis-eQTL of a transcription factor, IRF4, with 131 candidate target CpGs. Motif enrichment and IRF4 ChIP-seq analysis demonstrated that these target CpGs are enriched in IRF4 binding sites, suggesting an IRF4-mediated regulatory network. Our study highlights the utility of cell-type-specific meQTLs.

Introduction The optimal length for clinical follow-up of renal cell carcinoma (RCC) patients is unclear. We evaluated the impact of ISUP/WHO tumor grade and histological subtype on short- and long-term survival and risk of recurrence/metastasis in a large cohort of RCC patients. Patients and Methods We studied 1679 RCC patients from a single referral center in Italy. Adjusted hazard ratios for overall survival were estimated using Cox regression models. Adjusted absolute risk of developing recurrence or metastasis was computed considering competing risks of mortality. Results During up to 13 years of follow-up, 175 (10.4%) RCC patients died, of whom 92% beyond 5 years. Hazard ratio of grade IV clear cell carcinomas (ccRCC) was 3.82 compared to grade II. Notably, 33% of recurrences and 56% of distant metastases occurred beyond 5 years of follow-up. The estimated probabilities of recurrence/metastasis were 15% and 45% within and beyond 5 years of follow-up, respectively. After 5 years, the absolute risk of recurrences increased also for papillary renal cell carcinoma type I (35.2%) and grade I ccRCC (17%). Conclusion After 5 years of follow-up, both risk of mortality and recurrences or metastases were high and were modified by histological types and tumor grade. These data strongly support histology- and grade-tailored surveillance strategies and long-term follow-up for RCC patients.

We previously identified 10 lung adenocarcinoma susceptibility loci in a genome-wide association study (GWAS) conducted in the Female Lung Cancer Consortium in Asia (FLCCA), the largest genomic study of lung cancer among never-smoking women to date. Furthermore, household coal use for cooking and heating has been linked to lung cancer in Asia, especially in Xuanwei, China. We investigated the potential interaction between genetic susceptibility and coal use in FLCCA. We analyzed GWAS-data from Taiwan, Shanghai, and Shenyang (1472 cases; 1497 controls), as well as a separate study conducted in Xuanwei (152 cases; 522 controls) for additional analyses. We summarized genetic susceptibility using a polygenic risk score (PRS), which was the weighted sum of the risk-alleles from the 10 previously identified loci. We estimated associations between a PRS, coal use (ever/never), and lung adenocarcinoma with multivariable logistic regression models, and evaluated potential gene-environment interactions using likelihood ratio tests. There was a strong association between continuous PRS and lung adenocarcinoma among never coal users (Odds Ratio (OR) = 1.69 (95% Confidence Interval (CI) = 1.53, 1.87), p=1 × 10−26). This effect was attenuated among ever coal users (OR = 1.24 (95% CI: 1.03, 1.50), p = 0.02, p-interaction = 6 × 10−3). We observed similar attenuation among coal users from Xuanwei. Our study provides evidence that genetic susceptibility to lung adenocarcinoma among never-smoking Asian women is weaker among coal users. These results suggest that lung cancer pathogenesis may differ, at least partially, depending on exposure to coal combustion products. Notably, these novel findings are among the few instances of sub-multiplicative gene-environment interactions in the cancer literature.

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1β pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

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, Author summary Squamous cell carcinomas are specific type of cancer that can arise in multiple organs of the aerodigestive tract including the lung, oral cavity, oropharynx, larynx and esophagus. Previous studies have shown that aerodigestive squamous cell carcinomas share common environmental risk factors (tobacco smoking and alcohol intake). Here, we investigate genetic factors involved in the risk of aerodigestive squamous cell carcinomas as a group in a large genetic association study involving 13,887 cancer cases and 61,961 controls. We identified one genome-wide significant region within 2q33.1 and 3 other suggestive regions at 1q32.1, 5q31.2 and 19p13.11. Gene-based analyses also identify a list of SqCC-related genes that are involved in DNA damage response and epigenetic regulation. Our results suggest some overlap in the genetic factors influencing the risk of aerodigestive squamous cell carcinomas in European populations and highlights the importance of cross-cancer studies.

The German Social Accident Insurance supported this study (grant number FP 271)., Hovanec J, Siemiatycki J, Conway DI, Olsson A, Guenel P, Luce D, Jöckel KH, Pohlabeln H, Ahrens W, Karrasch S, Wichmann HE, Gustavsson P, Consonni D, Merletti F, Richiardi L, Lorenzo S, Fortes C, Parent MÉ, McLaughlin JR, Demers P, Landi MT, Caporaso N, Fernández-Tardón G, Zaridze D, Świątkowska B, Pándics T, Lissowska J, Fabianova E, Field JK, Mates D, Bencko V, Foretova L, Janout V, Kromhout H, Vermeulen R, Boffetta P, Straif K, Schüz J, Casjens S, Pesch B, Brüning T, Behrens T.

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

SummaryLung cancer is the leading cause of cancer death worldwide. Genome-wide association studies have revealed genetic risk factors, highlighting the role of smoking, family history, telomere regulation, and DNA damage-repair in lung cancer etiology. Many studies have focused on a single ethnic group to avoid confounding from variability in allele frequencies across populations; however, comprehensive multi-ethnic analyses may identify variants that are more likely to be causal. This large-scale, multi- ethnic meta-analyses identified 28 novel risk loci achieving genome-wide significance. Leading candidates were further studied using single-cell methods for evaluating DNA-damage. DNA-damage promoting activities were confirmed for selected genes by knockdown genes and overexpression studies.