Genome-wide association study identifies novel loci predisposing to cutaneous melanoma

Genetic variants in glutamine metabolic pathway genes predict cutaneous melanoma-specific survival

Glutamine dependence is a unique metabolic defect seen in cutaneous melanoma (CM), directly influencing the treatment and prognosis. Here, we investigated the associations between 6025 common single-nucleotide polymorphisms (SNPs) in 77 glutamine metabolic pathway genes with CM-specific survival (CMSS) using genotyping datasets from two published genome-wide association studies (GWASs). In the single-locus analysis, 76 SNPs were found to be significantly associated with CMSS (P < .050, false-positive report probability < 0.2 and Bayesian false discovery probability < 0.8) in the discovery dataset, of which seven SNPs were replicated in the validation dataset and three SNPs (HAL rs17676826T > C, LGSN rs12663017T > A, and NOXRED1 rs8012548A > G) independently predicted CMSS, with an effect-allele attributed adjusted hazards ratio of 1.52 (95% confidence interval = 1.19-1.93) and P < .001, 0.68 (0.54-0.87) and P = .002 and 0.62 (0.46-0.83) and P = .002, respectively. The model including the number of unfavorable genotypes (NUGs) of these three SNPs and covariates improved the five-year CMSS prediction (P = .012) than the one with other covariates only. Further expression quantitative trait loci (eQTL) analysis found that the LGSN rs12663017 A allele was significantly associated with increased messenger RNA (mRNA) expression levels (P = 8.89 × 10 -11 ) in lymphoblastoid cell lines of the 1000 Genomes Project database. In the analysis of the genotype tissue expression (GTEx) project datasets, HAL rs17676826 C and NOXRED1 rs8012548 G alleles were significantly associated with their mRNA expression levels in sun-exposed skin of the lower leg (P = 6.62 × 10-6 and 1.37 × 10-7 , respectively) and in sun-not-exposed suprapubic skin (P < .001 and 1.43 × 10-8 , respectively). Taken together, these genetic variants of glutamine-metabolic pathway genes may be promising predictors of survival in patients with CM.

Role of Immune Response, Inflammation, and Tumor Immune Response-Related Cytokines/Chemokines in Melanoma Progression

To investigate the role of tumor cytokines/chemokines in melanoma immune response, we estimated the proportions of immune cell subsets in melanoma tumors from The Cancer Genome Atlas, followed by evaluation of the association between cytokine/chemokine expression and these subsets. We then investigated the association of immune cell subsets, chemokines, and cytokines with patient survival. Finally, we evaluated the immune cell tumor-infiltrating lymphocyte (TIL) score for correlation with melanoma patient outcome in a separate cohort. There was good agreement between RNA sequencing estimation of T-cell subset and pathologist-determined TIL score. Expression levels of cytokines IL-12A, IFNG, and IL-10, and chemokines CXCL9 and CXCL10 were positively correlated with PDCD1, CTLA-4, and CD8+ T-cell subset, but negatively correlated with tumor purity (Bonferroni-corrected P < 0.05). In multivariable analysis, higher expression levels of cytokines IFN-γ and TGFB1, but not chemokines, were associated with improved overall survival. A higher expression level of CD8+ T-cell subset was also associated with improved overall survival (hazard ratio [HR] = 0.06, 95% confidence interval [CI] = 0.01-0.35, P = 0.002). Finally, multivariable analysis showed that patients with a brisk TIL score had improved melanoma-specific survival than those with a nonbrisk score (HR = 0.51, 95% CI = 0.27-0.98, P = 0.0423). These results suggest that the expression of specific tumor cytokines represents important biomarkers of melanoma immune response.

Association of rs8444 polymorphism in the LASS2 3'-UTR and bladder cancer risk in Chinese population

The aim of the present study was to explore the correlations between single nucleotide polymorphisms in LASS2 gene 3'-untranslated regions and bladder cancer risk in Chinese population. We first performed PCR and sequence for LASS2-3'-UTR in 105 bladder cancer patients and 100 control subjects. Next, multivariate logistic regression analysis was used to determine the relationship between single nucleotide polymorphisms frequency and susceptibility of bladder cancer, and clinical features in 105 cases. In addition, survival curves and Cox Regression analysis were used to investigate the effect of single nucleotide polymorphisms on clinical outcome in 58 cases. Finally, quantitative reverse-transcription PCR and immunohistochemical were performed to explore the influence of single nucleotide polymorphisms on LASS2 expression. We found that a single nucleotide polymorphism (rs8444 C>T) located in the 3'-UTR of LASS2 was significantly associated with the risk of bladder cancer. We also showed the frequency of rs8444 T genotype was higher in bladder cancer group and correlated with the risk of clinical prognosis. Yet, there were no significant correlations between T/C allele frequencies and the distributions of rs8444 genotype and tumor-node-metastasis stage, histological grade and distant metastasis in bladder cancer. Furthermore, we demonstrated that rs8444 C>T could affect LASS2 expression by single nucleotide polymorphism-related mRNA stability. Our results showed that LASS2-3'-UTR rs8444 C>T polymorphism was significantly associated with the individual risk and the poor overall survival of bladder cancer, suggesting that rs8444 TT genotype maybe act as an independent risk factor of susceptibility and clinical prognosis for bladder cancer in Chinese population.

HERCing: Structural and Functional Relevance of the Large HERC Ubiquitin Ligases

Homologous to the E6AP carboxyl terminus (HECT) and regulator of chromosome condensation 1 (RCC1)-like domain-containing proteins (HERCs) belong to the superfamily of ubiquitin ligases. HERC proteins are divided into two subfamilies, Large and Small HERCs. Despite their similarities in terms of both structure and domains, these subfamilies are evolutionarily very distant and result from a convergence phenomenon rather than from a common origin. Large HERC genes, HERC1 and HERC2, are present in most metazoan taxa. They encode very large proteins (approximately 5,000 amino acid residues in a single polypeptide chain) that contain more than one RCC1-like domain as a structural characteristic. Accumulating evidences show that these unusually large proteins play key roles in a wide range of cellular functions which include neurodevelopment, DNA damage repair, and cell proliferation. To better understand the origin, evolution, and function of the Large HERC family, this minireview provides with an integrated overview of their structure and function and details their physiological implications. This study also highlights and discusses how dysregulation of these proteins is associated with severe human diseases such as neurological disorders and cancer.

Genetic variants in Barrett's esophagus and esophageal adenocarcinoma: a literature review

Surveillance of Barrett's esophagus (BE) is a clinical challenge; metaplasia of the distal esophagus increases a patient's risk of esophageal adenocarcinoma (EAC) significantly but the actual percentage of patients who progress is low. The current screening recommendations require frequent endoscopy and biopsy, which has inherent risk, high cost, and operator variation. Identifying BE patients genetically who are at high risk of progressing could deemphasize the role of endoscopic screening and create an opportunity for early therapeutic intervention. Genetic alterations in germline DNA have been identified in other disease processes and allow for early intervention or surveillance well before disease develops. The genetic component of BE remains mostly unknown and only a few genome-wide association studies exist on this topic. This review summarizes the current literature available that examines genetic alterations in BE and EAC with a particular emphasis on clinical implications.

Immunomodulatory germline variation associated with the development of multiple primary melanoma (MPM)

Multiple primary melanoma (MPM) has been associated with a higher 10-year mortality risk compared to patients with single primary melanoma (SPM). Given that 3–8% of patients with SPM develop additional primary melanomas, new markers predictive of MPM risk are needed. Based on the evidence that the immune system may regulate melanoma progression, we explored whether germline genetic variants controlling the expression of 41 immunomodulatory genes modulate the risk of MPM compared to patients with SPM or healthy controls. By genotyping these 41 variants in 977 melanoma patients, we found that rs2071304, linked to the expression of SPI1, was strongly associated with MPM risk reduction (OR = 0.60; 95% CI = 0.45–0.81; p = 0.0007) when compared to patients with SPM. Furthermore, we showed that rs6695772, a variant affecting expression of BATF3, is also associated with MPM-specific survival (HR = 3.42; 95% CI = 1.57–7.42; p = 0.0019). These findings provide evidence that the genetic variation in immunomodulatory pathways may contribute to the development of secondary primary melanomas and also associates with MPM survival. The study suggests that inherited host immunity may play an important role in MPM development.

Genome-wide association studies and polygenic risk scores for skin cancer: clinically useful yet?

BACKGROUND

Genome-wide association studies (GWAS) have identified thousands of susceptibility variants, although most have been associated with small individual risk estimates that offer little predictive value. However, combining multiple variants into polygenic risk scores (PRS) may be more informative. Multiple studies have developed PRS composed of GWAS-identified variants for cutaneous cancers. This review highlights data from these studies.

OBJECTIVES

To review published GWAS and PRS studies for melanoma, cutaneous squamous cell carcinoma (cSCC) and basal cell carcinoma (BCC), and discuss their potential clinical utility.

METHODS

We searched PubMed and the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue to identify relevant studies.

RESULTS

Results from 21 GWAS (11 melanoma, 3 cSCC, 7 BCC) and 11 PRS studies are summarized. Six loci in pigmentation genes overlap between these three cancers (ASIP/RALY, IRF4, MC1R, OCA2, SLC45A2 and TYR). Additional loci overlap for cSCC/BCC and BCC/melanoma, but no other loci are shared between cSCC and melanoma. PRS for melanoma show roughly two-to-threefold increases in risk and modest improvements in risk prediction (2-7% increases). PRS are associated with twofold and threefold increases in risk of cSCC and BCC, respectively, with small improvements (2% increase) in predictive ability.

CONCLUSIONS

Existing data indicate that PRS may offer small, but potentially meaningful, improvements to risk prediction. Additional research is needed to clarify the potential utility of PRS in cutaneous carcinomas. Clinical translation will require well-powered validation studies incorporating known risk factors to evaluate PRS as tools for screening. What's already known about this topic? Over 50 susceptibility loci for melanoma, basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC) have been identified in genome-wide association studies (GWAS). Polygenic risk scores (PRS) using variants identified from GWAS have also been developed for melanoma, BCC and cSCC, and investigated with respect to clinical risk prediction. What does this study add? This review provides an overview of GWAS findings and the potential clinical utility of PRS for melanoma, BCC and cSCC.

Genetic variants in the calcium signaling pathway genes are associated with cutaneous melanoma-specific survival

Remodeling or deregulation of the calcium signaling pathway is a relevant hallmark of cancer including cutaneous melanoma (CM). In this study, using data from a published genome-wide association study (GWAS) from The University of Texas M.D. Anderson Cancer Center, we assessed the role of 41,377 common single-nucleotide polymorphisms (SNPs) of 167 calcium signaling pathway genes in CM survival. We used another GWAS from Harvard University as the validation dataset. In the single-locus analysis, 1830 SNPs were found to be significantly associated with CM-specific survival (CMSS; P ≤ 0.050 and false-positive report probability ≤ 0.2), of which 9 SNPs were validated in the Harvard study (P ≤ 0.050). Among these, three independent SNPs (i.e. PDE1A rs6750552 T>C, ITPR1 rs6785564 A>G and RYR3 rs2596191 C>A) had a predictive role in CMSS, with a meta-analysis-derived hazards ratio of 1.52 (95% confidence interval = 1.19-1.94, P = 7.21 × 10-4), 0.49 (0.33-0.73, 3.94 × 10-4) and 0.67 (0.53-0.86, 0.0017), respectively. Patients with an increasing number of protective genotypes had remarkably improved CMSS. Additional expression quantitative trait loci analysis showed that these genotypes were also significantly associated with mRNA expression levels of the genes. Taken together, these results may help us to identify prospective biomarkers in the calcium signaling pathway for CM prognosis.

Frequent Alteration of Annexin A9 and A10 in HPV-Negative Head and Neck Squamous Cell Carcinomas: Correlation with the Histopathological Differentiation Grade

The annexin protein superfamily has been implicated in multiple physiological and pathological processes, including carcinogenesis. Altered expression of various annexins has frequently been observed and linked to the development and progression of various human malignancies. However, information is lacking on the expression and clinical significance of annexin A9 (ANXA9) and A10 (ANXA10) in head and neck squamous cell carcinomas (HNSCC). ANXA9 and ANXA10 expression was evaluated in a large cohort of 372 surgically treated HPV-negative HNSCC patients and correlated with the clinicopathologic parameters and disease outcomes. Down-regulation of ANXA9 expression was found in 42% of HNSCC tissue samples, compared to normal epithelia. ANXA9 expression in tumors was significantly associated with oropharyngeal location and histological differentiation grade (P < 0.001). In marked contrast, ANXA10 expression was absent in normal epithelium, but variably detected in the cytoplasm of cancer cells. Positive ANXA10 expression was found in 64% of tumors, and was significantly associated with differentiation grade (P < 0.001), being also more frequent in oropharyngeal tumors (P = 0.019). These results reveal that the expression of both ANXA9 and ANXA10 is frequently altered in HNSCC and associated to the tumor differentiation grade, suggesting that they could be implicated in the pathogenesis of these cancers.

A functional variant rs1537373 in 9p21.3 region is associated with pancreatic cancer risk

9p21.3 has been identified as an unexpected hot point in multiple diseases GWAS including cancers, and we performed a two-stage case-control studies integrating functional assay strategy to find the potential functional variants modified susceptibility to pancreatic cancer (PC). An expanded Illumina HumanExome Beadchip of PC including 943 cases and 3908 controls was used to examine 39 tagSNPs in 9p21.3 and the promising single nucleotide polymorphism (SNP) was validated in stage 2 comprising 624 cases and 1048 controls. The strongest signal was rs6475609 (Odds ratio, OR = 0.81, 95% confidence interval, CI = 0.72-0.91) maps to the long non-coding RNA ANRIL. Bioinformatics analysis revealed rs1537373 lies in the linkage disequilibrium (LD) block which the rs6475609 tagged might have potential function and was also associated with a decreased risk of PC in both stages (OR = 0.82, 95% CI = 0.75-0.90 in combined analysis). Dual luciferase reporter assay and the electrophoretic mobility shift assay (EMSA) verified rs1537373 as the best candidate causative variant for influencing the activity of enhancer through differential binding to certain transcription factor. The expression quantitative trait loci (e-QTL) analysis indicated the genotypes of rs1537373 were associated with expression of CDKN2B gene (P dominant = 6.00 × 10^-4 ). In conclusion, our study provided evidence that rs1537373 in ANRIL may influence transcription factor binding and regulate CDKN2B expression, thus confer the susceptibility to PC.

A Nonsynonymous Variant in the GOLM1 Gene in Cutaneous Malignant Melanoma

Background

Statistically significant linkage of melanoma to chromosome 9q21 was previously reported in a Danish pedigree resource and independently confirmed in Utah high-risk pedigrees, indicating strong evidence that this region contains a melanoma predisposition gene.

Methods

Whole-exome sequencing of pairs of related melanoma case subjects from two pedigrees with evidence of 9q21 linkage was performed to identify the responsible predisposition gene. Candidate variants were tested for association with melanoma in an independent set of 454 unrelated familial melanoma case subjects and 396 unrelated cancer-free control subjects from Utah, and 1534 melanoma case subjects and 1146 noncancer control subjects from Texas (MD Anderson) via a two-sided Fisher exact test.

Results

A rare nonsynonymous variant in Golgi Membrane Protein 1 (GOLM1), rs149739829, shared in two hypothesized predisposition carriers in one linked pedigree was observed. Segregation of this variant in additional affected relatives of the index carriers was confirmed. A statistically significant excess of carriers of the variant was observed among Utah case subjects and control subjects (odds ratio [OR] = 9.81, 95% confidence interval [CI] = 8.35 to 11.26, P < .001) and statistically significantly confirmed in Texas case subjects and control subjects (OR = 2.45, 95% CI = 1.65 to 3.25, P = .02).

Conclusion

These findings support GOLM1 as a candidate melanoma predisposition gene.

Inherited Genetic Variants Associated with Melanoma BRAF/NRAS Subtypes

BRAF and NRAS mutations arise early in melanoma development, but their associations with low-penetrance melanoma susceptibility loci remain unknown. In the Genes, Environment and Melanoma Study, 1,223 European-origin participants had their incident invasive primary melanomas screened for BRAF/NRAS mutations and germline DNA genotyped for 47 single-nucleotide polymorphisms identified as low-penetrant melanoma-risk variants. We used multinomial logistic regression to simultaneously examine each single-nucleotide polymorphism's relationship to BRAF V600E, BRAF V600K, BRAF other, and NRAS+ relative to BRAF-/NRAS- melanoma adjusted for study features. IRF4 rs12203592*T was associated with BRAF V600E (odds ratio [OR] = 0.59, 95% confidence interval [CI] = 0.43-0.79) and V600K (OR = 0.65, 95% CI = 0.41-1.03), but not BRAF other or NRAS+ melanoma. A global test of etiologic heterogeneity (Pglobal = 0.001) passed false discovery (Pglobal = 0.0026). PLA2G6 rs132985*T was associated with BRAF V600E (OR = 1.32, 95% CI = 1.05-1.67) and BRAF other (OR = 1.82, 95% CI = 1.11-2.98), but not BRAF V600K or NRAS+ melanoma. The test for etiologic heterogeneity (Pglobal) was 0.005. The IRF4 rs12203592 associations were slightly attenuated after adjustment for melanoma-risk phenotypes. The PLA2G6 rs132985 associations were independent of phenotypes. IRF4 and PLA2G6 inherited genotypes may influence melanoma BRAF/NRAS subtype development.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

Gene expression imputation identifies candidate genes and susceptibility loci associated with cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer with genetic susceptibility loci identified in recent genome-wide association studies (GWAS). Transcriptome-wide association studies (TWAS) using imputed gene expression levels can identify additional gene-level associations. Here we impute gene expression levels in 6891 cSCC cases and 54,566 controls in the Kaiser Permanente Genetic Epidemiology Research in Adult Health and Aging (GERA) cohort and 25,558 self-reported cSCC cases and 673,788 controls from 23andMe. In a discovery-validation study, we identify 19 loci containing 33 genes whose imputed expression levels are associated with cSCC at false discovery rate < 10% in the GERA cohort and validate 15 of these candidate genes at Bonferroni significance in the 23andMe dataset, including eight genes in five novel susceptibility loci and seven genes in four previously associated loci. These results suggest genetic mechanisms contributing to cSCC risk and illustrate advantages and disadvantages of TWAS as a supplement to traditional GWAS analyses.

Genetic loci linked to susceptibility for the common skin cancer cutaneous squamous cell carcinoma (cSCC) have been identified by genome wide association studies (GWAS). Here, the authors impute gene expression levels from GWAS data to perform a transcriptome wide association study (TWAS), identifying five novel genetic loci linked to cSCC susceptibility.

Characterization of Nine Cancer-Associated Variants in Human DNA Polymerase κ

Specialized DNA damage-bypass Y-family DNA polymerases contribute to cancer prevention by providing cellular tolerance to DNA damage that can lead to mutations and contribute to cancer progression by increasing genomic instability. Y-family polymerases can also bypass DNA adducts caused by chemotherapy agents. One of the four human Y-family DNA polymerases, DNA polymerase (pol) κ, has been shown to be specific for bypass of minor groove adducts and inhibited by major groove adducts. In addition, mutations in the gene encoding pol κ are associated with different types of cancers as well as with chemotherapy responses. We characterized nine variants of pol κ whose identity was inferred from cancer-associated single nucleotide polymorphisms for polymerization activity on undamaged and damaged DNA, their abilities to extend from mismatched or damaged base pairs at primer termini, and overall stability and dynamics. We find that these pol κ variants generally fall into three categories: similar activity to wild-type (WT) pol κ (L21F, I39T, P169T, F192C, and E292K), more active than WT pol κ (S423R), and less active than pol κ (R219I, R298H, and Y432S). Of these, only pol κ variants R298H and Y432S had markedly reduced thermal stability. Molecular dynamics (MD) simulations with undamaged DNA revealed that the active variant F192C and more active variant S423R with either correct or incorrect incoming nucleotide mimic WT pol κ with the correct incoming nucleotide, whereas the less active variants R219I, R298H, and Y432S with the correct incoming nucleotide mimic WT pol κ with the incorrect incoming nucleotide. Thus, the observations from MD simulations suggest a possible explanation for the observed experimental results that pol κ adopts specific active and inactive conformations that depend on both the protein variant and the identity of the DNA adduct.

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.

Genetic variants in RORA and DNMT1 associated with cutaneous melanoma survival

Cutaneous melanoma (CM) is considered as a steroid hormone-related malignancy. However, few studies have evaluated the roles of genetic variants encoding steroid hormone receptor genes and their related regulators (SHR-related genes) in CM-specific survival (CMSS). Here, we performed a pathway-based analysis to evaluate genetic variants of 191 SHR-related genes in 858 CMSS patients using a dataset from a genome-wide association study (GWAS) from The University of Texas MD Anderson Cancer Center (MDACC), and then validated the results in an additional dataset of 409 patients from the Harvard GWAS. Using multivariate Cox proportional hazards regression analysis, we identified three-independent SNPs (RORA rs782917 G > A, RORA rs17204952 C > T and DNMT1 rs7253062 G > A) as predictors of CMSS, with a variant-allele attributed hazards ratio (HR) and 95% confidence interval of 1.62 (1.25-2.09), 1.60 (1.20-2.13) and 1.52 (1.20-1.94), respectively. Combined analysis of risk genotypes of these three SNPs revealed a decreased CMSS in a dose-response manner as the number of risk genotypes increased (ptrend  < 0.001); however, no improvement in the prediction model was observed (area under the curve [AUC] = 79.6-80.8%, p = 0.656), when these risk genotypes were added to the model containing clinical variables. Our findings suggest that genetic variants of RORA and DNMT1 may be promising biomarkers for CMSS, but these results needed to be validated in future larger studies.

Founder Mutations for Early Onset Melanoma as Revealed by Whole Exome Sequencing Suggests That This is Not Associated with the Increasing Incidence of Melanoma in Poland

Purpose

Germline mutations within melanoma susceptibility genes are present only in minority of melanoma patients and it is expected that additional genes will be discovered with next generation sequence technology and whole-exome sequencing (WES).

Materials and Methods

Herein we performed WES on a cohort of 96 unrelated Polish patients with melanoma diagnosed under the age of 40 years who all screened negative for the presence of CDKN2A variants. A replication study using a set of 1,200 melanoma patient DNA samples and similarly large series of healthy controls was undertaken.

Results

We selected 21 potentially deleterious variants in 20 genes (VRK1, MYCT1, DNAH14, CASC3, MS4A12, PRC1, WWOX, CARD6, EXO5, CASC3, CASP8AP2, STK33, SAMD11, CNDP2, CPNE1, EFCAB6, CABLES1, LEKR1, NUDT17, and RRP15), which were identified by WES and confirmed by Sanger sequencing for an association study. Evaluation of the allele distribution among carriers and their relatives in available family trios revealed that these variants were unlikely to account for many familial cases of melanoma. Replication study revealed no statistically significant differences between cases and controls.

Conclusion

Although most of the changes seemed to be neutral we could not exclude an association between variants in VRK1, CREB3L3, EXO5, and STK33 with melanoma risk.

MC1R variants as melanoma risk factors independent of at-risk phenotypic characteristics: a pooled analysis from the M-SKIP project

PURPOSE

Melanoma represents an important public health problem, due to its high case-fatality rate. Identification of individuals at high risk would be of major interest to improve early diagnosis and ultimately survival. The aim of this study was to evaluate whether MC1R variants predicted melanoma risk independently of at-risk phenotypic characteristics.

MATERIALS AND METHODS

Data were collected within an international collaboration - the M-SKIP project. The present pooled analysis included data on 3,830 single, primary, sporadic, cutaneous melanoma cases and 2,619 controls from seven previously published case-control studies. All the studies had information on MC1R gene variants by sequencing analysis and on hair color, skin phototype, and freckles, ie, the phenotypic characteristics used to define the red hair phenotype.

RESULTS

The presence of any MC1R variant was associated with melanoma risk independently of phenotypic characteristics (OR 1.60; 95% CI 1.36-1.88). Inclusion of MC1R variants in a risk prediction model increased melanoma predictive accuracy (area under the receiver-operating characteristic curve) by 0.7% over a base clinical model (P=0.002), and 24% of participants were better assessed (net reclassification index 95% CI 20%-30%). Subgroup analysis suggested a possibly stronger role of MC1R in melanoma prediction for participants without the red hair phenotype (net reclassification index: 28%) compared to paler skinned participants (15%).

CONCLUSION

The authors suggest that measuring the MC1R genotype might result in a benefit for melanoma prediction. The results could be a valid starting point to guide the development of scientific protocols assessing melanoma risk prediction tools incorporating the MC1R genotype.

modSaRa: a computationally efficient R package for CNV identification

Summary

Chromosomal copy number variation (CNV) refers to a polymorphism that a DNA segment presents deletion or duplication in the population. The computational algorithms developed to identify this type of variation are usually of high computational complexity. Here we present a user-friendly R package, modSaRa, designed to perform copy number variants identification. The package is developed based on a change-point based method with optimal computational complexity and desirable accuracy. The current version of modSaRa package is a comprehensive tool with integration of preprocessing steps and main CNV calling steps.

Availability and Implementation

modSaRa is an R package written in R, C ++ and Rcpp and is now freely available for download at http://c2s2.yale.edu/software/modSaRa .

Contact

heping.zhang@yale.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Associations between smoking behavior-related alleles and the risk of melanoma

Several studies have reported that cigarette smoking is inversely associated with the risk of melanoma. This study further tested whether incorporating genetic factors will provide another level of evaluation of mechanisms underlying the association between smoking and risk of melanoma. We investigated the association between SNPs selected from genome-wide association studies (GWAS) on smoking behaviors and risk of melanoma using 2,298 melanoma cases and 6,654 controls. Among 16 SNPs, three (rs16969968 [A], rs1051730 [A] and rs2036534 [C] in the 15q25.1 region) reached significance for association with melanoma risk in men (0.01 < = P values < = 0.02; 0.85 < = Odds Ratios (ORs) <= 1.20). There was association between the genetic scores based on the number of smoking behavior-risk alleles and melanoma risk with P-trend = 0.005 among HPFS. Further association with smoking behaviors indicating those three SNPs (rs16969968 [A], rs1051730 [A] and rs2036534 [C]) significantly associated with number of cigarettes smoked per day, CPD, with P = 0.009, 0.011 and 0.001 respectively. The SNPs rs215605 in the PDE1C gene and rs6265 in the BDNF gene significantly interacted with smoking status on melanoma risk (interaction P = 0.005 and P = 0.003 respectively). Our study suggests that smoking behavior-related SNPs are likely to play a role in melanoma development and the potential public health importance of polymorphisms in the CHRNA5-A3-B4 gene cluster. Further larger studies are warranted to validate the findings.

Selection of X-chromosome Inactivation Model

To address the complexity of the X-chromosome inactivation (XCI) process, we previously developed a unified approach for the association test for X-chromosomal single-nucleotide polymorphisms (SNPs) and the disease of interest, accounting for different biological possibilities of XCI: random, skewed, and escaping XCI. In the original study, we focused on the SNP-disease association test but did not provide knowledge regarding the underlying XCI models. One can use the highest likelihood ratio (LLR) to select XCI models (max-LLR approach). However, that approach does not formally compare the LLRs corresponding to different XCI models to assess whether the models are distinguishable. Therefore, we propose an LLR comparison procedure (comp-LLR approach), inspired by the Cox test, to formally compare the LLRs of different XCI models to select the most likely XCI model that describes the underlying XCI process. We conduct simulation studies to investigate the max-LLR and comp-LLR approaches. The simulation results show that compared with the max-LLR, the comp-LLR approach has higher probability of identifying the correct underlying XCI model for the scenarios when the underlying XCI process is random XCI, escaping XCI, or skewed XCI to the deleterious allele. We applied both approaches to a head and neck cancer genetic study to investigate the underlying XCI processes for the X-chromosomal genetic variants.

Ancestry inference using principal component analysis and spatial analysis: a distance-based analysis to account for population substructure

Background

Accurate inference of genetic ancestry is of fundamental interest to many biomedical, forensic, and anthropological research areas. Genetic ancestry memberships may relate to genetic disease risks. In a genome association study, failing to account for differences in genetic ancestry between cases and controls may also lead to false-positive results. Although a number of strategies for inferring and taking into account the confounding effects of genetic ancestry are available, applying them to large studies (tens thousands samples) is challenging. The goal of this study is to develop an approach for inferring genetic ancestry of samples with unknown ancestry among closely related populations and to provide accurate estimates of ancestry for application to large-scale studies.

Methods

In this study we developed a novel distance-based approach, Ancestry Inference using Principal component analysis and Spatial analysis (AIPS) that incorporates an Inverse Distance Weighted (IDW) interpolation method from spatial analysis to assign individuals to population memberships.

Results

We demonstrate the benefits of AIPS in analyzing population substructure, specifically related to the four most commonly used tools EIGENSTRAT, STRUCTURE, fastSTRUCTURE, and ADMIXTURE using genotype data from various intra-European panels and European-Americans. While the aforementioned commonly used tools performed poorly in inferring ancestry from a large number of subpopulations, AIPS accurately distinguished variations between and within subpopulations.

Conclusions

Our results show that AIPS can be applied to large-scale data sets to discriminate the modest variability among intra-continental populations as well as for characterizing inter-continental variation. The method we developed will protect against spurious associations when mapping the genetic basis of a disease. Our approach is more accurate and computationally efficient method for inferring genetic ancestry in the large-scale genetic studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4166-8) contains supplementary material, which is available to authorized users.

Integrative genetic analysis suggests that skin color modifies the genetic architecture of melanoma

Melanoma is the deadliest form of skin cancer and presents a significant health care burden in many countries. In addition to ultraviolet radiation in sunlight, the main causal factor for melanoma, genetic factors also play an important role in melanoma susceptibility. Although genome-wide association studies have identified many single nucleotide polymorphisms associated with melanoma, little is known about the proportion of disease risk attributable to these loci and their distribution throughout the genome. Here, we investigated the genetic architecture of melanoma in 1,888 cases and 990 controls of European non-Hispanic ancestry. We estimated the overall narrow-sense heritability of melanoma to be 0.18 (P < 0.03), indicating that genetics contributes significantly to the risk of sporadically-occurring melanoma. We then demonstrated that only a small proportion of this risk is attributable to known risk variants, suggesting that much remains unknown of the role of genetics in melanoma. To investigate further the genetic architecture of melanoma, we partitioned the heritability by chromosome, minor allele frequency, and functional annotations. We showed that common genetic variation contributes significantly to melanoma risk, with a risk model defined by a handful of genomic regions rather than many risk loci distributed throughout the genome. We also demonstrated that variants affecting gene expression in skin account for a significant proportion of the heritability, and are enriched among melanoma risk loci. Finally, by incorporating skin color into our analyses, we observed both a shift in significance for melanoma-associated loci and an enrichment of expression quantitative trait loci among melanoma susceptibility variants. These findings suggest that skin color may be an important modifier of melanoma risk. We speculate that incorporating skin color and other non-genetic factors into genetic studies may allow for an improved understanding of melanoma susceptibility and guide future investigations to identify melanoma risk genes.

Germline variation in inflammation-related pathways and risk of Barrett's oesophagus and oesophageal adenocarcinoma

OBJECTIVE

Oesophageal adenocarcinoma (OA) incidence has risen sharply in Western countries over recent decades. Local and systemic inflammation is considered an important contributor to OA pathogenesis. Established risk factors for OA and its precursor, Barrett's oesophagus (BE), include symptomatic reflux, obesity and smoking. The role of inherited genetic susceptibility remains an area of active investigation. Here, we explore whether germline variation related to inflammatory processes influences susceptibility to BE/OA.

DESIGN

We used data from a genomewide association study of 2515 OA cases, 3295 BE cases and 3207 controls. Our analysis included 7863 single-nucleotide polymorphisms (SNPs) in 449 genes assigned to five pathways: cyclooxygenase (COX), cytokine signalling, oxidative stress, human leucocyte antigen and nuclear factor-κB. A principal components-based analytic framework was employed to evaluate pathway-level and gene-level associations with disease risk.

RESULTS

We identified a significant signal for the COX pathway in relation to BE risk (p=0.0059, false discovery rate q=0.03), and in gene-level analyses found an association with microsomal glutathione-S-transferase 1 (MGST1); (p=0.0005, q=0.005). Assessment of 36 MGST1 SNPs identified 14 variants associated with elevated BE risk (q<0.05). Four of these were subsequently confirmed (p<5.5×10-5) in a meta-analysis encompassing an independent set of 1851 BE cases and 3496 controls, and are known strong expression quantitative trait loci for MGST1. Three such variants were associated with similar elevations in OA risk.

CONCLUSIONS

This study provides the most comprehensive evaluation of inflammation-related germline variation in relation to risk of BE/OA and suggests that variants in MGST1 influence disease susceptibility.

Two-stage genome-wide association study identifies a novel susceptibility locus associated with melanoma

Genome-wide association studies have identified 21 susceptibility loci associated with melanoma. These loci implicate genes affecting pigmentation, nevus count, telomere maintenance, and DNA repair in melanoma risk. Here, we report the results of a two-stage genome-wide association study of melanoma. The stage 1 discovery phase consisted of 4,842 self-reported melanoma cases and 286,565 controls of European ancestry from the 23andMe research cohort and the stage 2 replication phase consisted of 1,804 melanoma cases and 1,026 controls from the University of Texas M.D. Anderson Cancer Center. We performed a combined meta-analysis totaling 6,628 melanoma cases and 287,591 controls. Our study replicates 20 of 21 previously known melanoma-loci and confirms the association of the telomerase reverse transcriptase, TERT, with melanoma susceptibility at genome-wide significance. In addition, we uncover a novel polymorphism, rs187843643 (OR = 1.96; 95% CI = [1.54, 2.48]; P = 3.53 × 10⁻⁸), associated with melanoma. The SNP rs187842643 lies within a noncoding RNA 177kb downstream of BASP1 (brain associated protein-1). We find that BASP1 expression is suppressed in melanoma as compared with benign nevi, providing additional evidence for a putative role in melanoma pathogenesis.

Variants at the OCA2/HERC2 locus affect time to first cutaneous squamous cell carcinoma in solid organ transplant recipients collected using two different study designs

BACKGROUND

Variants at the oculocutaneous albinism 2 (OCA2)/HECT and RLD domain containing E3 ubiquitin protein ligase 2 (HERC2) locus have been associated with pigmentation phenotypes and risk of developing several types of skin cancer.

OBJECTIVES

To evaluate OCA2/HERC2 locus variants for their impact on time to develop cutaneous squamous cell carcinoma (cSCC) in organ transplant recipients (OTRs) who are at elevated risk of developing cSCC.

METHODS

Participants were solid OTRs ascertained from two centres (n = 125 and 261) with an average of 13·1 years of follow-up post-transplant. DNA was available for genotyping for all participants, in addition to medical records and questionnaire data. The Ohio State University study had a case-control design with prospective follow-up, and the University of California San Francisco study was a national cross-sectional survey with retrospective chart review.

RESULTS

OCA2 variants rs12913832 and rs916977 were significantly associated with time to first cSCC post-transplant. OTRs homozygous for the brown-eye alleles of rs916977 (GG) and rs12913832 (AA) had significant delays of time to first cSCC post-transplant compared with individuals homozygous for the blue-eye alleles (hazard ratio 0·34, P < 0·001 and hazard ratio 0·54, P = 0·012, respectively). Both variants were highly associated with eye colour in the combined studies (P < 0·001).

CONCLUSIONS

This study is the first to show an association between OCA2/HERC2 variants and time to first cSCC post-transplant. This may impact dermatological screening recommendations for high-risk populations.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.

Genetic variants in the metzincin metallopeptidase family genes predict melanoma survival

Metzincins are key molecules in the degradation of the extracellular matrix and play an important role in cellular processes such as cell migration, adhesion, and cell fusion of malignant tumors, including cutaneous melanoma (CM). We hypothesized that genetic variants of the metzincin metallopeptidase family genes would be associated with CM-specific survival (CMSS). To test this hypothesis, we first performed Cox proportional hazards regression analysis to evaluate the associations between genetic variants of 75 metzincin metallopeptidase family genes and CMSS using the dataset from the genome-wide association study (GWAS) from The University of Texas MD Anderson Cancer Center (MDACC) which included 858 non-Hispanic white patients with CM, and then validated using the dataset from the Harvard GWAS study which had 409 non-Hispanic white patients with invasive CM. Four independent SNPs (MMP16 rs10090371 C>A, ADAMTS3 rs788935 T>C, TLL2 rs10882807 T>C and MMP9 rs3918251 A>G) were identified as predictors of CMSS, with a variant-allele attributed hazards ratio (HR) of 1.73 (1.32-2.29, 9.68E-05), 1.46 (1.15-1.85, 0.002), 1.68 (1.31-2.14, 3.32E-05) and 0.67 (0.51-0.87, 0.003), respectively, in the meta-analysis of these two GWAS studies. Combined analysis of risk genotypes of these four SNPs revealed a decreased CMSS in a dose-response manner as the number of risk genotypes increased (P_trend < 0.001). An improvement was observed in the prediction model (area under the curve [AUC] = 81.4% vs. 78.6%), when these risk genotypes were added to the model containing non-genotyping variables. Our findings suggest that these genetic variants may be promising prognostic biomarkers for CMSS.

Genetic variants in the genes encoding rho GTPases and related regulators predict cutaneous melanoma-specific survival

Rho GTPases control cell division, motility, adhesion, vesicular trafficking and phagocytosis, which may affect progression and/or prognosis of cancers. Here, we investigated associations between genetic variants of Rho GTPases-related genes and cutaneous melanoma-specific survival (CMSS) by re-analyzing a published melanoma genome-wide association study (GWAS) and validating the results in another melanoma GWAS. In the single-locus analysis of 36,018 SNPs in 129 Rho-related genes, 427 SNPs were significantly associated with CMSS (p < 0.050 and false-positive report probability <0.2) in the discovery dataset, and five SNPs were replicated in the validation dataset. Among these, four SNPs (i.e., RHOU rs10916352 G > C, ARHGAP22 rs3851552 T > C, ARHGAP44 rs72635537 C > T and ARHGEF10 rs7826362 A > T) were independently predictive of CMSS (a meta-analysis derived p = 9.04 × 10-4 , 9.58 × 10-4 , 1.21 × 10-4 and 8.47 × 10-4 , respectively). Additionally, patients with an increasing number of unfavorable genotypes (NUGs) of these loci had markedly reduced CMSS in both discovery dataset and validation dataset (ptrend =1.47 × 10-7 and 3.12 × 10-5 ). The model including the NUGs and clinical variables demonstrated a significant improvement in predicting the five-year CMSS. Moreover, rs10916352C and rs3851552C alleles were significantly associated with an increased mRNA expression levels of RHOU (p = 1.8 × 10-6 ) and ARHGAP22 (p = 5.0 × 10-6 ), respectively. These results may provide promising prognostic biomarkers for CM personalized management and treatment.

Genetic variants of PDGF signaling pathway genes predict cutaneous melanoma survival

To investigate whether genetic variants of platelet-derived growth factor (PDGF) signaling pathway genes are associated with survival of cutaneous melanoma (CM) patients, we assessed associations of single-nucleotide polymorphisms in PDGF pathway with melanoma-specific survival in 858 CM patients of M.D. Anderson Cancer Center (MDACC). Additional data of 409 cases from Harvard University were also included for further analysis. We identified 13 SNPs in four genes (COL6A3, NCK2, COL5A1 and PRKCD) with a nominal P < 0.05 and false discovery rate (FDR) < 0.2 in MDACC dataset. Based on linkage disequilibrium, functional prediction and minor allele frequency, a representative SNP in each gene was selected. In the meta-analysis using MDACC and Harvard datasets, there were two SNPs associated with poor survival of CM patients: rs6707820 C>T in NCK2 (HR = 1.87, 95% CI = 1.35-2.59, P_meta= 1.53E-5); and rs2306574 T>C in PRKCD (HR = 1.73, 95% CI = 1.33-2.24, P_meta= 4.56E-6). Moreover, CM patients in MDACC with combined risk genotypes of these two loci had markedly poorer survival (HR = 2.47, 95% CI = 1.58-3.84, P < 0.001). Genetic variants of rs6707820 C>T in NCK2 and rs2306574 T>C in PRKCD of the PDGF signaling pathway may be biomarkers for melanoma survival.

Genetic Variants in WNT2B and BTRC Predict Melanoma Survival

Cutaneous melanoma (CM) is the most lethal skin cancer. The Wnt pathway has an impact on development, invasion, and metastasis of CM, thus likely affecting CM prognosis. Using data from a published genome-wide association study from The University of Texas MD Anderson Cancer Center, we assessed the associations of 19,830 common single-nucleotide polymorphisms (SNPs) in 151 Wnt pathway autosomal genes with CM-specific survival and then validated significant SNPs in another genome-wide association study from Harvard University. In the single-locus analysis, 1,855 SNPs were significantly associated with CM-specific survival at P < 0.05, of which 547 SNPs were still considered noteworthy after the correction by the false-positive report probability. In the replication, two SNPs remained significantly associated with CM-specific survival after multiple comparison correction. By performing functional prediction and stepwise selection, we identified two independent SNPs (i.e., WNT2B rs1175649 G>T and BTRC rs61873997 G>A) that showed a predictive role in CM-specific survival, with an effect-allele-attributed hazards ratio (adjusted hazards ratio) of 1.99 (95% confidence interval = 1.41-2.81, P = 8.10 × 10-5) and 0.61 (0.46-0.80, 3.12×10-4), respectively. Collectively, these variants in the Wnt pathway genes may be biomarkers for outcomes of patients with CM, if validated by larger studies.

From fish bowl to bedside: The power of zebrafish to unravel melanoma pathogenesis and discover new therapeutics

Melanoma is the most aggressive and deadliest form of skin cancer. A detailed knowledge of the cellular, molecular, and genetic events underlying melanoma progression is highly relevant to diagnosis, prognosis and risk stratification, and the development of new therapies. In the last decade, zebrafish have emerged as a valuable model system for the study of melanoma. Pathway conservation, coupled with the availability of robust genetic, transgenic, and chemical tools, has made the zebrafish a powerful model for identifying novel disease genes, visualizing cancer initiation, interrogating tumor-microenvironment interactions, and discovering new therapeutics that regulate melanocyte and melanoma development. In this review, we will give an overview of these studies, and highlight recent advancements that will help unravel melanoma pathogenesis and impact human disease.

No prognostic value added by vitamin D pathway SNPs to current prognostic system for melanoma survival

The prognostic improvement attributed to genetic markers over current prognostic system has not been well studied for melanoma. The goal of this study is to evaluate the added prognostic value of Vitamin D Pathway (VitD) SNPs to currently known clinical and demographic factors such as age, sex, Breslow thickness, mitosis and ulceration (CDF). We utilized two large independent well-characterized melanoma studies: the Genes, Environment, and Melanoma (GEM) and MD Anderson studies, and performed variable selection of VitD pathway SNPs and CDF using Random Survival Forest (RSF) method in addition to Cox proportional hazards models. The Harrell's C-index was used to compare the performance of model predictability. The population-based GEM study enrolled 3,578 incident cases of cutaneous melanoma (CM), and the hospital-based MD Anderson study consisted of 1,804 CM patients. Including both VitD SNPs and CDF yielded C-index of 0.85, which provided slight but not significant improvement by CDF alone (C-index = 0.83) in the GEM study. Similar results were observed in the independent MD Anderson study (C-index = 0.84 and 0.83, respectively). The Cox model identified no significant associations after adjusting for multiplicity. Our results do not support clinically significant prognostic improvements attributable to VitD pathway SNPs over current prognostic system for melanoma survival.

MC1R signaling. Intracellular partners and pathophysiological implications

The melanocortin-1 receptor (MC1R) preferentially expressed in melanocytes is best known as a key regulator of the synthesis of epidermal melanin pigments. Its paracrine stimulation by keratinocyte-derived melanocortins also activates DNA repair pathways and antioxidant defenses to build a complex, multifaceted photoprotective response. Many MC1R actions rely on cAMP-dependent activation of two transcription factors, MITF and PGC1α, but pleiotropic MC1R signaling also involves activation of mitogen-activated kinases and AKT. MC1R partners such as β-arrestins, PTEN and the E3 ubiquitin ligase MGRN1 differentially regulate these pathways. The MC1R gene is complex and polymorphic, with frequent variants associated with skin phenotypes and increased cancer risk. We review current knowledge of signaling from canonical MC1R, its splice isoforms and natural polymorphic variants. Recently discovered intracellular targets and partners are also discussed, to highlight the diversity of mechanisms that may contribute to normal and pathological variation of pigmentation and sensitivity to solar radiation-induced damage. This article is part of a Special Issue entitled: Melanocortin Receptors - edited by Ya-Xiong Tao.

Common Variants Confer Susceptibility to Barrett's Esophagus: Insights from the First Genome-Wide Association Studies

Eight loci have been identified by the two genome-wide association studies of Barrett's esophagus that have been conducted to date. Esophageal adenocarcinoma cases were included in the second study following evidence that predisposing genetic variants for this cancer overlap with those for Barrett's esophagus. Genes with roles in embryonic development of the foregut are adjacent to 6 of the loci identified (FOXF1, BARX1, FOXP1, GDF7, TBX5, and ALDH1A2). An additional locus maps to a gene with known oncogenic potential (CREB-regulated transcription coactivator 1), but expression quantitative trait data implicates yet another gene involved in esophageal development (PBX4). These results strongly support a model whereby dysregulation of genes involved in esophageal and thoracic development increases susceptibility to Barrett's esophagus and esophageal adenocarcinoma, probably by reducing anatomical antireflux mechanisms. An additional signal at 6p21 in the major histocompatibility complex also reinforces evidence that immune and inflammatory response to reflux is involved in the development of both diseases. All of the variants identified are intronic or intergenic rather than coding and are presumed to be or to mark regulatory variants. As with genome-wide association studies of other diseases, the functional variants at each locus are yet to be identified and the genes affected need confirming. In this chapter as well as discussing the biology behind each genome-wide association signal, we review the requirements for successfully conducting genome-wide association studies and discuss how progress in understanding the genetic variants that contribute to Barrett's esophagus/esophageal adenocarcinoma susceptibility compares to other cancers.

Genetic variants in the integrin signaling pathway genes predict cutaneous melanoma survival

To identify genetic variants involved in prognosis of cutaneous melanoma (CM), we investigated associations of single nucleotide polymorphisms (SNPs) of genes in the integrin signaling pathway with CM survival by re-analyzing a published genome-wide association study (GWAS) from The University of Texas M.D. Anderson Cancer Center (MDACC) and then validated significant SNPs in another GWAS from Harvard University. In the MDACC study, 1,148 SNPs were significantly associated with CM-specific survival (CMSS) (p ≤ 0.050 and false-positive report probability ≤ 0.20), and nine SNPs were validated in the Harvard study (p ≤ 0.050). Among these, three independent SNPs (i.e., DOCK1 rs11018104 T > A, rs35748949 C > T and PAK2 rs1718404 C > T) showed a predictive role in CMSS, with an effect-allele attributed adjusted hazards ratio [adjHR of 1.50 (95% confidence interval (CI) = 1.18-1.90, p = 7.46E-04), 1.53 (1.18-1.97, 1.18E-03) and 0.58 (0.45-0.76, 5.60E-05), respectively]. Haplotype analysis revealed that a haplotype carrying two risk alleles A-T in DOCK1 was associated with the poorest survival in both MDACC (adjHR = 1.73, 95% CI = 1.19-2.50, p = 0.004) and Harvard (adjHR = 1.95, 95% CI = 1.14-3.33, p = 0.010) studies. In addition, patients with an increasing number of unfavorable genotypes (NUGs) for these three SNPs had a poorer survival. Incorporating NUGs with clinical variables showed a significantly improved ability to classify CMSS (AUC increased from 86.8% to 88.6%, p = 0.031). Genetic variants in the integrin signaling pathway may independently or jointly modulate the survival of CM patients. Further large, prospective studies are needed to validate these findings.

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.

Genetic variants in the PIWI-piRNA pathway gene DCP1A predict melanoma disease-specific survival

The Piwi-piRNA pathway is important for germ cell maintenance, genome integrity, DNA methylation and retrotransposon control and thus may be involved in cancer development. In this study, we comprehensively analyzed prognostic roles of 3,116 common SNPs in PIWI-piRNA pathway genes in melanoma disease-specific survival. A published genome-wide association study (GWAS) by The University of Texas M.D. Anderson Cancer Center was used to identify associated SNPs, which were later validated by another GWAS from the Harvard Nurses' Health Study and Health Professionals Follow-up Study. After multiple testing correction, we found that there were 27 common SNPs in two genes (PIWIL4 and DCP1A) with false discovery rate < 0.2 in the discovery dataset. Three tagSNPs (i.e., rs7933369 and rs508485 in PIWIL4; rs11551405 in DCP1A) were replicated. The rs11551405 A allele, located at the 3' UTR microRNA binding site of DCP1A, was associated with an increased risk of melanoma disease-specific death in both discovery dataset [adjusted Hazards ratio (HR) = 1.66, 95% confidence interval (CI) = 1.21-2.27, p =1.50 × 10^-3 ] and validation dataset (HR = 1.55, 95% CI = 1.03-2.34, p = 0.038), compared with the C allele, and their meta-analysis showed an HR of 1.62 (95% CI, 1.26-2.08, p =1.55 × 10^-4 ). Using RNA-seq data from the 1000 Genomes Project, we found that DCP1A mRNA expression levels increased significantly with the A allele number of rs11551405. Additional large, prospective studies are needed to validate these findings.

Genetic markers of pigmentation are novel risk loci for uveal melanoma

While the role of genetic risk factors in the etiology of uveal melanoma (UM) has been strongly suggested, the genetic susceptibility to UM is currently vastly unexplored. Due to shared epidemiological risk factors between cutaneous melanoma (CM) and UM, in this study we have selected 28 SNPs identified as risk variants in previous genome-wide association studies on CM or CM-related host phenotypes (such as pigmentation and eye color) and tested them for association with UM risk. By logistic regression analysis of 272 UM cases and 1782 controls using an additive model, we identified five variants significantly associated with UM risk, all passing adjustment for multiple testing. The three most significantly associated variants rs12913832 (OR = 0.529, 95% CI 0.415–0.673; p = 8.47E-08), rs1129038 (OR = 0.533, 95% CI 0.419–0.678; p = 1.19E-07) and rs916977 (OR = 0.465, 95% CI 0.339–0.637; p = 3.04E-07) are correlated (r² > 0.5) and map at 15q12 in the region of HERC2/OCA2, which determines eye-color in the human population. Our data provides first evidence that the genetic factors associated with pigmentation traits are risk loci of UM susceptibility.

Association of Interferon Regulatory Factor-4 Polymorphism rs12203592 With Divergent Melanoma Pathways

BACKGROUND

Solar elastosis and neval remnants are histologic markers characteristic of divergent melanoma pathways linked to differences in age at onset, host phenotype, and sun exposure. However, the association between these pathway markers and newly identified low-penetrance melanoma susceptibility loci remains unknown.

METHODS

In the Genes, Environment and Melanoma (GEM) Study, 2103 Caucasian participants had first primary melanomas that underwent centralized pathology review. For 47 single-nucleotide polymorphisms (SNPs) previously identified as low-penetrant melanoma risk variants, we used multinomial logistic regression to compare melanoma with solar elastosis and melanoma with neval remnants simultaneously to melanoma with neither of these markers, excluding melanomas with both markers. All statistical tests were two-sided.

RESULTS

IRF4 rs12203592 was the only SNP to pass the false discovery threshold in baseline models adjusted for age, sex, and study center. rs12203592*T was associated positively with melanoma with solar elastosis (odds ratio [OR] = 1.47, 95% confidence interval [CI] = 1.18 to 1.82) and inversely with melanoma with neval remnants (OR = 0.65, 95% CI = 0.48 to 0.87) compared with melanoma with neither marker (P global = 3.78 x 10(-08)). Adjusting for phenotypic characteristics and total sun exposure hours did not materially affect rs12203592's associations. Distinct early- and late-onset age distributions were observed in patients with IRF4 rs12203592 [CC] and [TT] genotypes, respectively.

CONCLUSIONS

Our findings suggest a role of IRF4 rs12203592 in pathway-specific risk for melanoma development. We hypothesize that IRF4 rs12203592 could underlie in part the bimodal age distribution reported for melanoma and linked to the divergent pathways.

The genomic landscape of cutaneous melanoma

Somatic mutation analysis of melanoma has been performed at the single gene level extensively over the past several decades. This has provided considerable insight into the critical pathways controlling melanoma initiation and progression. During the last 5 yr, next-generation sequencing (NGS) has enabled even more comprehensive mutational screening at the level of multigene panels, exomes and genomes. These studies have uncovered many new and unexpected players in melanoma development. The recent landmark study from The Cancer Genome Atlas (TCGA) consortium describing the genomic architecture of 333 cutaneous melanomas provides the largest and broadest analysis to date on the somatic aberrations underlying melanoma genesis. It thus seems timely to review the mutational landscape of melanoma and highlight the key genes and cellular pathways that appear to drive this cancer.

Biochip-Based Genotyping Assay for Detection of Polymorphisms in Pigmentation Genes Associated with Cutaneous Melanoma

AIMS

The purpose of the study was to develop a new assay for genotyping nine single nucleotide polymorphisms (SNPs) that are known to be associated with melanoma.

METHODS

Two-stage single tube polymerase chain reaction (PCR) followed by hybridization on a biochip was developed and applied in the study.

RESULTS

A total of nine SNPs were selected from five genes: MC1R (rs1805006, rs1805007, rs1805009, rs11547464), HERC2 (rs12913832), OCA2 (rs1800407), SLC45A2 (rs16891982), TYR (rs1393350), and a SNP from the intergenic locus rs12896399 were used for the synthesis of ssDNAs via a single-stage PCR process. The assays were performed on a biochip-based platform that is capable of SNP genotyping via a single reaction-tube PCR, followed by on chip hybridization. We tested 69 DNAs obtained from healthy persons and demonstrated the assays' ability to discriminate all three genotypes for almost all of the SNPs.

CONCLUSIONS

The developed approach proved robust, suggesting that it might be useful for the personalized genotyping of large cohorts of patients.

Identification of Susceptibility Loci for Cutaneous Squamous Cell Carcinoma

We report a genome-wide association study of cutaneous squamous cell carcinoma conducted among non-Hispanic white members of the Kaiser Permanente Northern California health care system. The study includes a genome-wide screen of 61,457 members (6,891 cases and 54,566 controls) genotyped on the Affymetrix Axiom European array and a replication phase involving an independent set of 6,410 additional members (810 cases and 5,600 controls). Combined analysis of screening and replication phases identified 10 loci containing single-nucleotide polymorphisms (SNPs) with P-values < 5 × 10(-8). Six loci contain genes in the pigmentation pathway; SNPs at these loci appear to modulate squamous cell carcinoma risk independently of the pigmentation phenotypes. Another locus contains HLA class II genes studied in relation to elevated squamous cell carcinoma risk following immunosuppression. SNPs at the remaining three loci include an intronic SNP in FOXP1 at locus 3p13, an intergenic SNP at 3q28 near TP63, and an intergenic SNP at 9p22 near BNC2. These findings provide insights into the genetic factors accounting for inherited squamous cell carcinoma susceptibility.