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

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

Massively parallel reporter assays of melanoma risk variants identify MX2 as a gene promoting melanoma.

Genome-wide association studies (GWAS) have identified ~20 melanoma susceptibility loci, most of which are not functionally characterized. Here we report an approach integrating massively-parallel reporter assays (MPRA) with cell-type-specific epigenome and expression quantitative trait loci (eQTL) to identify susceptibility genes/variants from multiple GWAS loci. From 832 high-LD variants, we identify 39 candidate functional variants from 14 loci displaying allelic transcriptional activity, a subset of which corroborates four colocalizing melanocyte cis-eQTL genes. Among these, we further characterize the locus encompassing the HIV-1 restriction gene, MX2 (Chr21q22.3), and validate a functional intronic variant, rs398206. rs398206 mediates the binding of the transcription factor, YY1, to increase MX2 levels, consistent with the cis-eQTL of MX2 in primary human melanocytes. Melanocyte-specific expression of human MX2 in a zebrafish model demonstrates accelerated melanoma formation in a BRAFV600E background. Our integrative approach streamlines GWAS follow-up studies and highlights a pleiotropic function of MX2 in melanoma susceptibility.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility.

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Interactive effect between ATPase-related genes and early-life tobacco smoke exposure on bronchial hyper-responsiveness detected in asthma-ascertained families.

BACKGROUND: A positional cloning study of bronchial hyper-responsiveness (BHR) at the 17p11 locus in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families showed significant interaction between early-life environmental tobacco smoke (ETS) exposure and genetic variants located in DNAH9. This gene encodes the heavy chain subunit of axonemal dynein, which is involved with ATP in the motile cilia function.Our goal was to identify genetic variants at other genes interacting with ETS in BHR by investigating all genes belonging to the 'ATP-binding' and 'ATPase activity' pathways which include DNAH9, are targets of cigarette smoke and play a crucial role in the airway inflammation. METHODS: Family-based interaction tests between ETS-exposed and unexposed BHR siblings were conducted in 388 EGEA families. Twenty single-nucleotide polymorphisms (SNP) showing interaction signals (p≤5.10-3) were tested in the 253 Saguenay-Lac-Saint-Jean (SLSJ) families. RESULTS: One of these SNPs was significantly replicated for interaction with ETS in SLSJ families (p=0.003). Another SNP reached the significance threshold after correction for multiple testing in the combined analysis of the two samples (p=10-5). Results were confirmed using both a robust log-linear test and a gene-based interaction test. CONCLUSION: The SNPs showing interaction with ETS belong to the ATP8A1 and ABCA1 genes, which play a role in the maintenance of asymmetry and homeostasis of lung membrane lipids.

New susceptibility loci for cutaneous melanoma risk and progression revealed using a porcine model.

Despite major advances, it is estimated that a large part of melanoma predisposing genes remains to be discovered. Animal models of spontaneous diseases are valuable tools and experimental crosses can be used to identify and fine-map new susceptibility loci associated with melanoma. We performed a Genome-Wide Association Study (GWAS) of melanoma occurrence and progression (clinical ulceration and presence of metastasis) in a porcine model of spontaneous melanoma, the MeLiM pig. Five loci on chromosomes 2, 5, 7, 8 and 16 showed genome-wide significant associations (p < 5 × 10-6) with either one of these phenotypes. Suggestive associations (p < 5 × 10-5) were also found at 16 additional loci. Moreover, comparison of the porcine results to those reported by human melanoma GWAS indicated shared association signals notably at CDKAL1 and TERT loci but also nearby CCND1, FTO, PLA2G6 and TMEM38B-RAD23B loci. Extensive search of the literature revealed a potential key role of genes at the identified porcine loci in tumor invasion (DST, PLEKHA5, CBY1, LIMK2 and ETV5) and immune response modulation (ETV5, HERC3 and DICER1) of the progression phenotypes. These biological processes are consistent with the clinico-pathological features of MeLiM tumors and can open new routes for future melanoma research in humans.

Assessing the Incremental Contribution of Common Genomic Variants to Melanoma Risk Prediction in Two Population-Based Studies.

It is unclear to what degree genomic and traditional (phenotypic and environmental) risk factors overlap in their prediction of melanoma risk. We evaluated the incremental contribution of common genomic variants (in pigmentation, nevus, and other pathways) and their overlap with traditional risk factors, using data from two population-based case-control studies from Australia (n = 1,035) and the United Kingdom (n = 1,460) that used the same questionnaires. Polygenic risk scores were derived from 21 gene regions associated with melanoma and odds ratios from published meta-analyses. Logistic regression models were adjusted for age, sex, center, and ancestry. Adding the polygenic risk score to a model with traditional risk factors increased the area under the receiver operating characteristic curve (AUC) by 2.3% (P = 0.003) for Australia and by 2.8% (P = 0.002) for Leeds. Gene variants in the pigmentation pathway, particularly MC1R, were responsible for most of the incremental improvement. In a cross-tabulation of polygenic by traditional tertile risk scores, 59% (Australia) and 49% (Leeds) of participants were categorized in the same (concordant) tertile. Of participants with low traditional risk, 9% (Australia) and 21% (Leeds) had high polygenic risk. Testing of genomic variants can identify people who are susceptible to melanoma despite not having a traditional phenotypic risk profile.

A comprehensive genome-wide analysis of melanoma Breslow thickness identifies interaction between CDC42 and SCIN genetic variants.

Breslow thickness (BT) is a major prognostic factor of cutaneous melanoma (CM), the most fatal skin cancer. The genetic component of BT has only been explored by candidate gene studies with inconsistent results. Our objective was to uncover the genetic factors underlying BT using an hypothesis-free genome-wide approach. Our analysis strategy integrated a genome-wide association study (GWAS) of single nucleotide polymorphisms (SNPs) for BT followed by pathway analysis of GWAS outcomes using the gene-set enrichment analysis (GSEA) method and epistasis analysis within BT-associated pathways. This strategy was applied to two large CM datasets with Hapmap3-imputed SNP data: the French MELARISK study for discovery (966 cases) and the MD Anderson Cancer Center study (1,546 cases) for replication. While no marginal effect of individual SNPs was revealed through GWAS, three pathways, defined by gene ontology (GO) categories were significantly enriched in genes associated with BT (false discovery rate ≤5% in both studies): hormone activity, cytokine activity and myeloid cell differentiation. Epistasis analysis, within each significant GO, identified a statistically significant interaction between CDC42 and SCIN SNPs (pmeta-int =2.2 × 10(-6) , which met the overall multiple-testing corrected threshold of 2.5 × 10(-6) ). These two SNPs (and proxies) are strongly associated with CDC42 and SCIN gene expression levels and map to regulatory elements in skin cells. This interaction has important biological relevance since CDC42 and SCIN proteins have opposite effects in actin cytoskeleton organization and dynamics, a key mechanism underlying melanoma cell migration and invasion.

Interaction between the DNAH9 gene and early smoke exposure in bronchial hyperresponsiveness.

A previous genome-wide linkage scan of bronchial hyperresponsiveness (BHR) in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families, performed in the presence of a gene×early-life environmental tobacco smoke (ETS) exposure interaction, showed the strongest interaction in the 17p11 region where linkage was detected only among unexposed siblings. Our goal was to conduct fine-scale mapping of 17p11 to identify single nucleotide polymorphisms (SNPs) interacting with ETS that influence BHR.Analyses were performed in 388 French EGEA asthmatic families, using a two-step strategy: 1) selection of SNPs displaying family-based association test (FBAT) association signals (p≤0.01) with BHR in unexposed siblings, and 2) a FBAT homogeneity test between exposed and unexposed siblings plus a robust log-linear interaction test.A single SNP reached the threshold (p≤3×10(-3)) for significant interaction with ETS using both interaction tests, after accounting for multiple testing. Results were replicated in 253 French-Canadian families, but not in 341 UK families, probably due in part to differences in phenotypic features between datasets.The SNP showing significant interaction with ETS belongs toDNAH9(dynein, axonemal, heavy chain 9), a promising candidate gene involved in respiratory cilia mobility and associated with primary ciliary dyskinesia, a disease associated with abnormalities of pulmonary function.

Genome-wide meta-analysis identifies five new susceptibility loci for cutaneous malignant melanoma.

Thirteen common susceptibility loci have been reproducibly associated with cutaneous malignant melanoma (CMM). We report the results of an international 2-stage meta-analysis of CMM genome-wide association studies (GWAS). This meta-analysis combines 11 GWAS (5 previously unpublished) and a further three stage 2 data sets, totaling 15,990 CMM cases and 26,409 controls. Five loci not previously associated with CMM risk reached genome-wide significance (P < 5 × 10(-8)), as did 2 previously reported but unreplicated loci and all 13 established loci. Newly associated SNPs fall within putative melanocyte regulatory elements, and bioinformatic and expression quantitative trait locus (eQTL) data highlight candidate genes in the associated regions, including one involved in telomere biology.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Fine mapping of genetic susceptibility loci for melanoma reveals a mixture of single variant and multiple variant regions.

At least 17 genomic regions are established as harboring melanoma susceptibility variants, in most instances with genome-wide levels of significance and replication in independent samples. Based on genome-wide single nucleotide polymorphism (SNP) data augmented by imputation to the 1,000 Genomes reference panel, we have fine mapped these regions in over 5,000 individuals with melanoma (mainly from the GenoMEL consortium) and over 7,000 ethnically matched controls. A penalized regression approach was used to discover those SNP markers that most parsimoniously explain the observed association in each genomic region. For the majority of the regions, the signal is best explained by a single SNP, which sometimes, as in the tyrosinase region, is a known functional variant. However in five regions the explanation is more complex. At the CDKN2A locus, for example, there is strong evidence that not only multiple SNPs but also multiple genes are involved. Our results illustrate the variability in the biology underlying genome-wide susceptibility loci and make steps toward accounting for some of the "missing heritability."

The effect on melanoma risk of genes previously associated with telomere length.

Telomere length has been associated with risk of many cancers, but results are inconsistent. Seven single nucleotide polymorphisms (SNPs) previously associated with mean leukocyte telomere length were either genotyped or well-imputed in 11108 case patients and 13933 control patients from Europe, Israel, the United States and Australia, four of the seven SNPs reached a P value under .05 (two-sided). A genetic score that predicts telomere length, derived from these seven SNPs, is strongly associated (P = 8.92x10(-9), two-sided) with melanoma risk. This demonstrates that the previously observed association between longer telomere length and increased melanoma risk is not attributable to confounding via shared environmental effects (such as ultraviolet exposure) or reverse causality. We provide the first proof that multiple germline genetic determinants of telomere length influence cancer risk.

Association of genetic variants in CDK6 and XRCC1 with the risk of dysplastic nevi in melanoma-prone families.

Dysplastic nevi (DN) is a strong risk factor for cutaneous malignant melanoma (CMM), and it frequently occurs in melanoma-prone families. To identify genetic variants for DN, we genotyped 677 tagSNPs in 38 melanoma candidate genes that are involved in pigmentation, DNA repair, cell cycle control, and melanocyte proliferation pathways in a total of 504 individuals (310 with DN, 194 without DN) from 53 melanoma-prone families (23 CDKN2A mutation positive and 30 negative). Conditional logistic regression, conditioning on families, was used to estimate the association between DN and each single-nucleotide polymorphism (SNP) separately, adjusted for age, sex, CMM, and CDKN2A status. P-values for SNPs in the same gene were combined to yield gene-specific P-values. Two genes, CDK6 (cyclin-dependent kinase 6) and XRCC1, were significantly associated with DN after Bonferroni correction for multiple testing (P=0.0001 and 0.00025, respectively), whereas neither gene was significantly associated with CMM. Associations for CDK6 SNPs were stronger in CDKN2A mutation-positive families (rs2079147, Pinteraction=0.0033), whereas XRCC1 SNPs had similar effects in mutation-positive and -negative families. The association for one of the associated SNPs in XRCC1 (rs25487) was replicated in two independent data sets (random-effect meta-analysis: P<0.0001). Our findings suggest that some genetic variants may contribute to DN risk independently of their association with CMM in melanoma-prone families.