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

Mapping chromatin interactions at melanoma susceptibility loci uncovers distant cis-regulatory gene targets

Genome-wide association studies (GWASs) of melanoma risk have identified 68 independent signals at 54 loci. For most loci, specific functional variants and their respective target genes remain to be established. Capture-HiC is an assay that links fine-mapped risk variants to candidate target genes by comprehensively mapping chromatin interactions. We performed a melanoma GWAS region-focused capture-HiC assay in human primary melanocytes to identify physical interactions between fine-mapped risk variants and potential causal melanoma-susceptibility genes. Overall, chromatin-interaction data alone nominated potential causal genes for 61 of the 68 melanoma risk signals, identifying many candidates beyond those reported by previous studies. We further integrated these data with epigenomic (chromatin state, accessibility), gene expression (expression quantitative trait locus [eQTL]/transcriptome-wide association study [TWAS]), DNA methylation (methylation QTL [meQTL]/methylome-wide association study [MWAS]), and massively parallel reporter assay (MPRA) data generated from melanoma-relevant cell types to prioritize potentially cis-regulatory variants and their respective candidate gene targets. From the set of fine-mapped variants across these loci, we identified 140 prioritized credible causal variants linked to 195 candidate genes at 42 risk signals. In addition, we developed an integrative scoring system to facilitate candidate gene prioritization, integrating melanocyte and melanoma datasets. Notably, at several GWAS risk signals, we observed long-range chromatin connections (500 kb to >1 Mb) with distant candidate target genes. We validated several such cis-regulatory interactions using CRISPR inhibition, providing evidence for known cancer driver genes MDM4 and CBL, as well as the SRY-box transcription factor SOX4, as likely melanoma risk genes.

Polygenic Risk Score Improves Melanoma Risk Assessment in a Patient Cohort from the Veneto Region of Italy

Recent genome-wide association studies (GWASs) have identified many single nucleotide polymorphisms (SNPs) that alone weakly affect melanoma risk, but their combined effect on a polygenic risk score (PRS) can have a far bigger impact on estimating risk. However, the PRS is not yet at the stage of being utilized in clinical practice, and further evidence is needed. In this study, 270 melanoma patients fulfilling the criteria for a suspected genetic predisposition but with a negative genetic test for high/medium-penetrance genes were genotyped for 57 SNPs selected in previous GWASs to construct a PRS model. We found a significantly higher mean PRS57 in all melanoma cases than in controls (0.58 vs. 0.00, p < 0.001), and the mean PRS57 in multiple primary melanoma cases was twice that in single melanoma cases (0.689 vs. 0.362, p = 0.025). Interestingly, our results confirm the association of the PRS57 not only with other melanoma risk factors but also with a younger age at diagnosis. This evidence supports the potentially powerful discriminative role of PRS in the selection of high-risk patients who should undergo stricter surveillance protocols.

Mapping chromatin interactions at melanoma susceptibility loci and cell-type specific dataset integration uncovers distant gene targets of cis-regulation

Genome-wide association studies (GWAS) of melanoma risk have identified 68 independent signals at 54 loci. For most loci, specific functional variants and their respective target genes remain to be established. Capture-HiC is an assay that links fine-mapped risk variants to candidate target genes by comprehensively mapping cell-type specific chromatin interactions. We performed a melanoma GWAS region-focused capture-HiC assay in human primary melanocytes to identify physical interactions between fine-mapped risk variants and potential causal melanoma susceptibility genes. Overall, chromatin interaction data alone nominated potential causal genes for 61 of the 68 melanoma risk signals, identifying many candidates beyond those reported by previous studies. We further integrated these data with cell-type specific epigenomic (chromatin state, accessibility), gene expression (eQTL/TWAS), DNA methylation (meQTL/MWAS), and massively parallel reporter assay (MPRA) data to prioritize potentially cis-regulatory variants and their respective candidate gene targets. From the set of fine-mapped variants across these loci, we identified 140 prioritized candidate causal variants linked to 195 candidate genes at 42 risk signals. In addition, we developed an integrative scoring system to facilitate candidate gene prioritization, integrating melanocyte and melanoma datasets. Notably, at several GWAS risk signals we observed long-range chromatin connections (500 kb to >1 Mb) with distant candidate target genes. We validated several such cis-regulatory interactions using CRISPR inhibition, providing evidence for known cancer driver genes MDM4 and CBL, as well as the SRY-box transcription factor SOX4, as likely melanoma risk genes.

Melanocortin-1 receptor (MC1R): a review for dermatologists

Melanocortin-1 receptor (MC1R) and its variants have a pivotal role in melanin synthesis. However, MC1R has been associated to non-pigmentary pathways related to DNA-repair activities and inflammation. The aim of this review is to provide an up-to-date overview about the role of MC1R in the skin. Specifically, after summarizing the current knowledge about MC1R structure and polymorphisms, we report data concerning the correlation between MC1R, phenotypic traits, skin aging, other diseases and skin cancers and their risk assessment through genetic testing.

OSCAA: A two-dimensional Gaussian mixture model for copy number variation association analysis

Copy number variants (CNVs) are prevalent in the human genome and are found to have a profound effect on genomic organization and human diseases. Discovering disease-associated CNVs is critical for understanding the pathogenesis of diseases and aiding their diagnosis and treatment. However, traditional methods for assessing the association between CNVs and disease risks adopt a two-stage strategy conducting quantitative CNV measurements first and then testing for association, which may lead to biased association estimation and low statistical power, serving as a major barrier in routine genome-wide assessment of such variation. In this article, we developed One-Stage CNV-disease Association Analysis (OSCAA), a flexible algorithm to discover disease-associated CNVs for both quantitative and qualitative traits. OSCAA employs a two-dimensional Gaussian mixture model that is built upon the PCs from copy number intensities, accounting for technical biases in CNV detection while simultaneously testing for their effect on outcome traits. In OSCAA, CNVs are identified and their associations with disease risk are evaluated simultaneously in a single step, taking into account the uncertainty of CNV identification in the statistical model. Our simulations demonstrated that OSCAA outperformed the existing one-stage method and traditional two-stage methods by yielding a more accurate estimate of the CNV-disease association, especially for short CNVs or CNVs with weak signals. In conclusion, OSCAA is a powerful and flexible approach for CNV association testing with high sensitivity and specificity, which can be easily applied to different traits and clinical risk predictions.

Impact of iris color on uveal melanoma-related outcomes in 7245 patients at a single ocular oncology center

PURPOSE

To determine the relationship between iris color and uveal melanoma (UM)-related metastasis and death in a large cohort of patients from a single ocular oncology center.

DESIGN

Retrospective case series.

SUBJECTS

Patients diagnosed with UM between February 1971 and August 2007.

METHODS

Patient information was obtained from chart documentation.

MAIN OUTCOME MEASURES

UM-related metastasis and death.

RESULTS

Out of 7245 patients, iris color was blue in 3702 (51%), green in 1458 (20%), and brown in 2085 (29%). Mean age was 58 ± 15 years and mean tumor thickness was 5.5 ± 3.3 millimeters. Some clinical features differed between iris color groups, with the blue irides group having a larger proportion of post-equatorial tumors with significantly closer proximity to the foveola and optic disc compared to the brown irides group. At a mean follow-up of 75 months, there was no statistically significant difference in metastasis between the various iris color groups. On univariate analysis, those with blue irides showed a higher incidence of UM-related death compared to the green and brown irides groups (8.3%, 5.9% and 7.5% respectively, p value = 0.02). Kaplan-Meier event free survival from UM-related death significantly differed only between the blue and green irides groups (p value = 0.007) with the green irides group showing the highest survival. However, on multivariate analysis, iris color was not predictive of UM-related death.

CONCLUSION

Iris color was not predictive of UM-related metastasis or death. However, Kaplan-Meier survival at 20 years was poorest for blue irides group compared to green.

Identification of genomic-wide genetic links between cutaneous melanoma and obesity-related physical traits via cFDR

BACKGROUND

Both epidemiological and clinical studies have suggested the comorbidity between cutaneous melanoma (CM) and obesity-related physical traits. However, it remains unclear about their shared genetic architecture.

OBJECTIVE

To determine the shared genetic architecture between CM and obesity-related physical traits through conditional false discovery rate (cFDR) analysis.

METHOD

Quantile-quantile plots were firstly built to assess the pleiotropic enrichment of shared single nucleotide polymorphisms between CM and each trait. Then, cFDR and conjunctional cFDR (ccFDR) were used to identify the shared risk loci between CM and each trait. Moreover, the functional evaluation of shared risk genes was carried out through analyses of expression quantitative trait loci (eQTL), Kyoto Encyclopedia of Genes and Genomes and gene ontology, respectively. Finally, single-cell sequence analysis was performed to locate the expression of eQTL-mapped genes in tissues.

RESULTS

Successive pleiotropic enrichment was found between CM and 5 obesity-related traits or height. 24 shared risk loci were identified between CM and 13 traits except appendicular lean mass using ccFDR analysis, with 17 novel and 4 validated loci. The functions of ccFDR-identified and eQTL-mapped genes were revealed to be mainly involved in cellular senescence, proliferation, meiotic nuclear division, cell cycle, and the metabolism of lipid, cholesterol and glucose. Single-cell sequence analysis showed that keratinocytes contribute to the occurrence and aggressiveness of CM through secreting paracrine cytokines.

CONCLUSION

Our findings demonstrate the significant genetic correlation between CM and obesity-related physical traits, which may provide a novel genetical basis for the pathogenesis and treatment of CM.

Malignancies in Prader-Willi Syndrome: Results From a Large International Cohort and Literature Review

CONTEXT

Prader-Willi syndrome (PWS) is a complex disorder combining hypothalamic dysfunction, neurodevelopmental delay, hypotonia, and hyperphagia with risk of obesity and its complications. PWS is caused by the loss of expression of the PWS critical region, a cluster of paternally expressed genes on chromosome 15q11.2-q13. As life expectancy of patients with PWS increases, age-related diseases like malignancies might pose a new threat to health.

OBJECTIVE

To investigate the prevalence and risk factors of malignancies in patients with PWS and to provide clinical recommendations for cancer screening.

METHODS

We included 706 patients with PWS (160 children, 546 adults). We retrospectively collected data from medical records on past or current malignancies, the type of malignancy, and risk factors for malignancy. Additionally, we searched the literature for information about the relationship between genes on chromosome 15q11.2-q13 and malignancies.

RESULTS

Seven adults (age range, 18-55 years) had been diagnosed with a malignancy (acute lymphoblastic leukemia, intracranial hemangiopericytoma, melanoma, stomach adenocarcinoma, biliary cancer, parotid adenocarcinoma, and colon cancer). All patients with a malignancy had a paternal 15q11-13 deletion. The literature review showed that several genes on chromosome 15q11.2-q13 are related to malignancies.

CONCLUSION

Malignancies are rare in patients with PWS. Therefore, screening for malignancies is only indicated when clinically relevant symptoms are present, such as unexplained weight loss, loss of appetite, symptoms suggestive of paraneoplastic syndrome, or localizing symptoms. Given the increased cancer risk associated with obesity, which is common in PWS, participation in national screening programs should be encouraged.

OSCAA: A Two-Dimensional Gaussian Mixture Model for Copy Number Variation Association Analysis

Copy number variants (CNVs) are prevalent in the human genome which provide profound effect on genomic organization and human diseases. Discovering disease associated CNVs is critical for understanding the pathogenesis of diseases and aiding their diagnosis and treatment. However, traditional methods for assessing the association between CNVs and disease risks adopt a two-stage strategy conducting quantitative CNV measurements first and then testing for association, which may lead to biased association estimation and low statistical power, serving as a major barrier in routine genome wide assessment of such variation. In this article, we developed OSCAA, a flexible algorithm to discover disease associated CNVs for both quantitative and qualitative traits. OSCAA employs a two-dimensional Gaussian mixture model that is built upon the principal components from copy number intensities, accounting for technical biases in CNV detection while simultaneously testing for their effect on outcome traits. In OSCAA, CNVs are identified and their associations with disease risk are evaluated simultaneously in a single step, taking into account the uncertainty of CNV identification in the statistical model. Our simulations demonstrated that OSCAA outperformed the existing one-stage method and traditional two-stage methods by yielding a more accurate estimate of the CNV-disease association, especially for short CNVs or CNVs with weak signal. In conclusion, OSCAA is a powerful and flexible approach for CNV association testing with high sensitivity and specificity, which can be easily applied to different traits and clinical risk predictions.

Molecular Classifiers in Skin Cancers: Challenges and Promises

Skin cancers are common and heterogenous malignancies affecting up to two in three Australians before age 70. Despite recent developments in diagnosis and therapeutic strategies, the mortality rate and costs associated with managing patients with skin cancers remain high. The lack of well-defined clinical and histopathological features makes their diagnosis and classification difficult in some cases and the prognostication difficult in most skin cancers. Recent advancements in large-scale "omics" studies, including genomics, transcriptomics, proteomics, metabolomics and imaging-omics, have provided invaluable information about the molecular and visual landscape of skin cancers. On many occasions, it has refined tumor classification and has improved prognostication and therapeutic stratification, leading to improved patient outcomes. Therefore, this paper reviews the recent advancements in omics approaches and appraises their limitations and potential for better classification and stratification of skin cancers.

Skin Phototype and Disease: A Comprehensive Genetic Approach to Pigmentary Traits Pleiotropy Using PRS in the GCAT Cohort

Human pigmentation has largely been associated with different disease prevalence among populations, but most of these studies are observational and inconclusive. Known to be genetically determined, pigmentary traits have largely been studied by Genome-Wide Association Study (GWAS), mostly in Caucasian ancestry cohorts from North Europe, identifying robustly, several loci involved in many of the pigmentary traits. Here, we conduct a detailed analysis by GWAS and Polygenic Risk Score (PRS) of 13 pigmentary-related traits in a South European cohort of Caucasian ancestry (n = 20,000). We observed fair phototype strongly associated with non-melanoma skin cancer and other dermatoses and confirmed by PRS-approach the shared genetic basis with skin and eye diseases, such as melanoma (OR = 0.95), non-melanoma skin cancer (OR = 0.93), basal cell carcinoma (OR = 0.97) and darker phototype with vitiligo (OR = 1.02), cataracts (OR = 1.04). Detailed genetic analyses revealed 37 risk loci associated with 10 out of 13 analyzed traits, and 16 genes significantly associated with at least two pigmentary traits. Some of them have been widely reported, such as MC1R, HERC2, OCA2, TYR, TYRP1, SLC45A2, and some novel candidate genes C1QTNF3, LINC02876, and C1QTNF3-AMACR have not been reported in the GWAS Catalog, with regulatory potential. These results highlight the importance of the assess phototype as a genetic proxy of skin functionality and disease when evaluating open mixed populations.

Association of TYR SNP rs1042602 with Melanoma Risk and Prognosis

Cutaneous melanoma is the most aggressive of skin tumors. In order to discover new biomarkers that could help us improve prognostic prediction in melanoma patients, we have searched for germline DNA variants associated with melanoma progression. Thus, after exome sequencing of a set of melanoma patients and healthy control individuals, we identified rs1042602, an SNP within TYR, as a good candidate. After genotyping rs1042602 in 1025 patients and 773 healthy donors, we found that the rs1042602-A allele was significantly associated with susceptibility to melanoma (CATT test: p = 0.0035). Interestingly, we also observed significant differences between patients with good and bad prognosis (5 years of follow-up) (n = 664) (CATT test for all samples p = 0.0384 and for men alone p = 0.0054). Disease-free-survival (DFS) analyses also showed that patients with the A allele had shorter DFS periods. In men, the association remained significant even in a multivariate Cox Proportional-hazards model, which was adjusted for age at diagnosis, Breslow thickness, ulceration and melanoma subtype (HR 0.4; 95% confidence interval (CI) 0.20-0.83; p = 0.0139). Based on our results, we propose that rs1042602-A is a risk allele for melanoma, which also seems to be responsible for a poorer prognosis of the disease, particularly in men.

Common genetic variants associated with melanoma risk or naevus count in patients with wildtype MC1R melanoma

BACKGROUND

Hypomorphic MC1R variants are the most prevalent genetic determinants of melanoma risk in the white population. However, the genetic background of patients with wildtype (WT) MC1R melanoma is poorly studied.

OBJECTIVES

To analyse the role of candidate common genetic variants on the melanoma risk and naevus count in Spanish patients with WT MC1R melanoma.

METHODS

We examined 753 individuals with WT MC1R from Spain (497 patients and 256 controls). We used OpenArray reverse-transcriptase polymerase chain reaction to genotype a panel of 221 common genetic variants involved in melanoma, naevogenesis, hormonal pathways and proinflammatory pathways. Genetic models were tested using multivariate logistic regression models. Nonparametric multifactor dimensionality reduction (MDR) was used to detect gene-gene interactions within each biological subgroup of variants.

RESULTS

We found that variant rs12913832 in the HERC2 gene, which is associated with blue eye colour, increased melanoma risk in individuals with WT MC1R [odds ratio (OR) 1·97, 95% confidence interval (CI) 1·48-2·63; adjusted P < 0·001; corrected P < 0·001]. We also observed a trend between the rs3798577 variant in the oestrogen receptor alpha gene (ESR1) and a lower naevus count, which was restricted to female patients with WT MC1R (OR 0·51, 95% CI 0·33-0·79; adjusted P = 0·002; corrected P = 0·11). This sex-dependent association was statistically significant in a larger cohort of patients with melanoma regardless of their MC1R status (n = 1497; OR 0·71, 95% CI 0·57-0·88; adjusted P = 0·002), reinforcing the hypothesis of an association between hormonal pathways and susceptibility to melanocytic proliferation. Last, the MDR analysis revealed four genetic combinations associated with melanoma risk or naevus count in patients with WT MC1R.

CONCLUSIONS

Our data suggest that epistatic interaction among common variants related to melanocyte biology or proinflammatory pathways might influence melanocytic proliferation in individuals with WT MC1R. What is already known about this topic? Genetic variants in the MC1R gene are the most prevalent melanoma genetic risk factor in the white population. Still, 20-40% of cases of melanoma occur in individuals with wildtype MC1R. Multiple genetic variants have a pleiotropic effect in melanoma and naevogenesis. Additional variants in unexplored pathways might also have a role in melanocytic proliferation in these patients. Epidemiological evidence suggests an association of melanocytic proliferation with hormonal pathways and proinflammatory pathways. What does this study add? Variant rs12913832 in the HERC2 gene, which is associated with blue eye colour, increases the melanoma risk in individuals with wildtype MC1R. Variant rs3798577 in the oestrogen receptor gene is associated with naevus count regardless of the MC1R status in female patients with melanoma. We report epistatic interactions among common genetic variants with a role in modulating the risk of melanoma or the number of naevi in individuals with wildtype MC1R. What is the translational message? We report a potential role of hormonal signalling pathways in melanocytic proliferation, providing a basis for better understanding of sex-based differences observed at the epidemiological level. We show that gene-gene interactions among common genetic variants might be responsible for an increased risk for melanoma development in individuals with a low-risk phenotype, such as darkly pigmented hair and skin.

Assessment of genetic susceptibility to multiple primary cancers through whole-exome sequencing in two large multi-ancestry studies

BACKGROUND

Up to one of every six individuals diagnosed with one cancer will be diagnosed with a second primary cancer in their lifetime. Genetic factors contributing to the development of multiple primary cancers, beyond known cancer syndromes, have been underexplored.

METHODS

To characterize genetic susceptibility to multiple cancers, we conducted a pan-cancer, whole-exome sequencing study of individuals drawn from two large multi-ancestry populations (6429 cases, 165,853 controls). We created two groupings of individuals diagnosed with multiple primary cancers: (1) an overall combined set with at least two cancers across any of 36 organ sites and (2) cancer-specific sets defined by an index cancer at one of 16 organ sites with at least 50 cases from each study population. We then investigated whether variants identified from exome sequencing were associated with these sets of multiple cancer cases in comparison to individuals with one and, separately, no cancers.

RESULTS

We identified 22 variant-phenotype associations, 10 of which have not been previously discovered and were significantly overrepresented among individuals with multiple cancers, compared to those with a single cancer.

CONCLUSIONS

Overall, we describe variants and genes that may play a fundamental role in the development of multiple primary cancers and improve our understanding of shared mechanisms underlying carcinogenesis.

Coronary artery disease and cancer: a significant resemblance

Cancer and coronary artery disease (CAD) are two of the most common causes of death, and they frequently coexist, especially as the world's population ages. CAD can develop prior to or following cancer diagnosis, as well as a side effect of cancer treatment. CAD develops as complex interactions of lifestyle and hereditary variables, just like the development of the most complex and non-communicable diseases. Cancer is caused by both external/acquired factors (tobacco, food, physical activity, alcohol consumption, epigenetic alterations) and internal/inherited factors (genetic mutations, hormones, and immunological diseases). The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated protein 9 (Cas9) system has recently emerged as a strong tool for gene therapy for both cancer as well as CAD treatment due to its great accuracy and efficiency. A deeper understanding of the complex link between CAD and cancer should lead to better prevention, faster detection, and safer treatment strategies.

Genetic variants in CYP2B6 and HSD17B12 associated with risk of squamous cell carcinoma of the head and neck

Polycyclic aromatic hydrocarbons (PAH) and tobacco-specific nitrosamines (TSNA) metabolism-related genes play an important role in the development of cancers. We assessed the associations of genetic variants in genes involved in the metabolism of PAHs and TSNA with risk of squamous cell carcinoma of the head and neck (SCCHN) in European populations using two published genome-wide association study datasets. In the single-locus analysis, we identified two SNPs (rs145533669 and rs35246205) in CYP2B6 to be associated with risk of SCCHN (P = 1.57 × 10-4 and .004, respectively), two SNPs (EPHX1 rs117522494 and CYP2B6 rs145533669) to be associated with risk of oropharyngeal cancer (P = .001 and .004, respectively), and one SNP (rs4359199 in HSD17B12) to be associated with risk of oral cancer (P = .006). A significant interaction effect was found between rs4359199 and drinking status on risks of SCCHN and oropharyngeal cancer (P < .05). eQTL and sQTL analyzes revealed that two SNPs (CYP2B6 rs35246205 and HSD17B12 rs4359199) were correlated with alternative splicing or mRNA expression levels of the corresponding genes in liver cells (P < .05 for both). In silico functional annotation suggested that these two SNPs may regulate mRNA expression by affecting the binding of transcription factors. Results from phenome-wide association studies presented significant associations between these genes and risks of other cancers, smoking behavior and alcohol dependence (P < .05). Thus, our study provided some insight into the underlying genetic mechanism of head and neck cancer, which warrants future functional validation.

Genome-Wide Association Study Suggests the Variant rs7551288*A within the DHCR24 Gene Is Associated with Poor Overall Survival in Melanoma Patients

Simple Summary

The aim of this work was to investigate prognostic genetic factors in melanoma patients. Phenotypic and disease data as well as biomaterial were collected after informed consent from patients followed up in a Skin Cancer Center of a University clinic. Genome-wide analysis (GWAS) was performed with survival data of 556 melanoma patients and genetic data including more than 300,000 common polymorphisms. The SNP rs7551288 reached suggestive genome-wide significance (p = 2 × 10⁻⁶). This intronic variant of the DHCR24 gene is involved in the cholesterol synthesis pathway. Further analyses and a literature review suggest an important role of this locus for the clinical course of disease in melanoma patients.

Abstract

Melanoma incidence rates are high among individuals with fair skin and multiple naevi. Established prognostic factors are tumour specific, and less is known about prognostic host factors. A total of 556 stage I to stage IV melanoma patients from Germany with phenotypic and disease-specific data were analysed; 64 of these patients died of melanoma after a median follow-up time of 8 years. Germline DNA was assessed by the HumanCoreExome BeadChip and data of 356,384 common polymorphisms distributed over all 23 chromosomes were used for a genome-wide analysis. A suggestive genome-wide significant association of the intronic allele rs7551288*A with diminished melanoma-specific survival was detected (p = 2 × 10⁻⁶). The frequency of rs7551288*A was 0.43 and was not associated with melanoma risk, hair and eye colour, tanning and total naevus count. Cox regression multivariate analyses revealed a 5.31-fold increased risk of melanoma-specific death for patients with the rs7551288 A/A genotype, independent of tumour thickness, ulceration and stage of disease at diagnoses. The variant rs7551288 belongs to the DHCR24 gene, which encodes Seladin-1, an enzyme involved in the biosynthesis of cholesterol. Further investigations are needed to confirm this genetic variant as a novel prognostic biomarker and to explore whether specific treatment strategies for melanoma patients might be derived from it.

Molecular Pathology of Skin Melanoma: Epidemiology, Differential Diagnostics, Prognosis and Therapy Prediction

Similar to other malignancies, TCGA network efforts identified the detailed genomic picture of skin melanoma, laying down the basis of molecular classification. On the other hand, genome-wide association studies discovered the genetic background of the hereditary melanomas and the susceptibility genes. These genetic studies helped to fine-tune the differential diagnostics of malignant melanocytic lesions, using either FISH tests or the myPath gene expression signature. Although the original genomic studies on skin melanoma were mostly based on primary tumors, data started to accumulate on the genetic diversity of the progressing disease. The prognostication of skin melanoma is still based on staging but can be completed with gene expression analysis (DecisionDx). Meanwhile, this genetic knowledge base of skin melanoma did not turn to the expected wide array of target therapies, except the BRAF inhibitors. The major breakthrough of melanoma therapy was the introduction of immune checkpoint inhibitors, which showed outstanding efficacy in skin melanoma, probably due to their high immunogenicity. Unfortunately, beyond BRAF, KIT mutations and tumor mutation burden, no clinically validated predictive markers exist in melanoma, although several promising biomarkers have been described, such as the expression of immune-related genes or mutations in the IFN-signaling pathway. After the initial success of either target or immunotherapies, sooner or later, relapses occur in the majority of patients, due to various induced genetic alterations, the diagnosis of which could be developed to novel predictive genetic markers.

Genome-wide association study of actinic keratosis identifies new susceptibility loci implicated in pigmentation and immune regulation pathways

Actinic keratosis (AK) is a common precancerous cutaneous neoplasm that arises on chronically sun-exposed skin. AK susceptibility has a moderate genetic component, and although a few susceptibility loci have been identified, including IRF4, TYR, and MC1R, additional loci have yet to be discovered. We conducted a genome-wide association study of AK in non-Hispanic white participants of the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (n = 63,110, discovery cohort), with validation in the Mass-General Brigham (MGB) Biobank cohort (n = 29,130). We identified eleven loci (P < 5 × 10-8), including seven novel loci, of which four novel loci were validated. In a meta-analysis (GERA + MGB), one additional novel locus, TRPS1, was identified. Genes within the identified loci are implicated in pigmentation (SLC45A2, IRF4, BNC2, TYR, DEF8, RALY, HERC2, and TRPS1), immune regulation (FOXP1 and HLA-DQA1), and cell signaling and tissue remodeling (MMP24) pathways. Our findings provide novel insight into the genetics and pathogenesis of AK susceptibility.

Citrus-Gene interaction and melanoma risk in the UK Biobank

High citrus consumption may increase melanoma risk; however, little is known about the biological mechanisms of this association, or whether it is modified by genetic variants. We conducted a genome-wide analysis of gene-citrus consumption interactions on melanoma risk among 1563 melanoma cases and 193 296 controls from the UK Biobank. Both the 2-degrees-of-freedom (df) joint test of genetic main effect and gene-environment (G-E) interaction and the standard 1-df G-E interaction test were performed. Three index SNPs (lowest P-value SNP among highly correlated variants [r²  > .6]) were identified from among the 365 genome-wide significant 2-df test results (rs183783391 on chromosome 3 [MITF], rs869329 on chromosome 9 [MTAP] and rs11446223 on chromosome 16 [DEF8]). Although all three were statistically significant for the 2-df test (4.25e-08, 1.98e-10 and 4.93e-13, respectively), none showed evidence of interaction according to the 1-df test (P = .73, .24 and .12, respectively). Eight nonindex, 2-df test significant SNPs on chromosome 16 were significant (P < .05) according to the 1-df test, providing evidence of citrus-gene interaction. Seven of these SNPs were mapped to AFG3L1P (rs199600347, rs111822773, rs113178244, rs3803683, rs73283867, rs78800020, rs73283871), and one SNP was mapped to GAS8 (rs74583214). We identified several genetic loci that may elucidate the association between citrus consumption and melanoma risk. Further studies are needed to confirm these findings.

Association of Melanoma-Risk Variants with Primary Melanoma Tumor Prognostic Characteristics and Melanoma-Specific Survival in the GEM Study

Genome-wide association studies (GWAS) and candidate pathway studies have identified low-penetrant genetic variants associated with cutaneous melanoma. We investigated the association of melanoma-risk variants with primary melanoma tumor prognostic characteristics and melanoma-specific survival. The Genes, Environment, and Melanoma Study enrolled 3285 European origin participants with incident invasive primary melanoma. For each of 47 melanoma-risk single nucleotide polymorphisms (SNPs), we used linear and logistic regression modeling to estimate, respectively, the per allele mean changes in log of Breslow thickness and odds ratios for presence of ulceration, mitoses, and tumor-infiltrating lymphocytes (TILs). We also used Cox proportional hazards regression modeling to estimate the per allele hazard ratios for melanoma-specific survival. Passing the false discovery threshold (p = 0.0026) were associations of IRF4 rs12203592 and CCND1 rs1485993 with log of Breslow thickness, and association of TERT rs2242652 with presence of mitoses. IRF4 rs12203592 also had nominal associations (p < 0.05) with presence of mitoses and melanoma-specific survival, as well as a borderline association (p = 0.07) with ulceration. CCND1 rs1485993 also had a borderline association with presence of mitoses (p = 0.06). MX2 rs45430 had nominal associations with log of Breslow thickness, presence of mitoses, and melanoma-specific survival. Our study indicates that further research investigating the associations of these genetic variants with underlying biologic pathways related to tumor progression is warranted.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Genetic variants of SDCCAG8 and MAGI2 in mitosis-related pathway genes are independent predictors of cutaneous melanoma-specific survival

Mitosis is a prognostic factor for cutaneous melanoma (CM), but accurate mitosis detection in CM tissues is difficult. Therefore, the 8th Edition of the American Joint Committee on Cancer staging system has removed the mitotic rate as a category criterion of the tumor T-category, based on the evidence that the mitotic rate was not an independent prognostic factor for melanoma survival. As single-nucleotide polymorphisms (SNPs) have been shown to be potential predictors for cutaneous melanoma-specific survival (CMSS), we investigated the potential prognostic value of SNPs in mitosis-related pathway genes in CMSS by analyzing their associations with outcomes of 850 CM patients from The University of Texas MD Anderson Cancer Center in a discovery dataset and validated the findings in another dataset of 409 CM patients from the Harvard University Nurses' Health Study and Health Professionals Follow-up Study. In both datasets, we identified two SNPs (SDCCAG8 rs10803138 G>A and MAGI2 rs3807694 C>T) as independent prognostic factors for CMSS, with adjusted allelic hazards ratios of 1.49 (95% confidence interval = 1.17-1.90, P = .001) and 1.45 (1.13-1.86, P = .003), respectively. Furthermore, their combined unfavorable alleles also predicted a poor survival in both discovery and validation datasets in a dose-response manner (Ptrend  = .0006 and .0001, respectively). Additional functional analysis revealed that both SDCCAG8 rs10803138 A and MAGI2 rs3807694 T alleles were associated with elevated mRNA expression levels in normal tissues. Therefore, these findings suggest that SDCCAG8 rs10803138 G>A and MAGI2 rs3807694 C>T are independent prognostic biomarkers for CMSS, possibly by regulating the mRNA expression of the corresponding genes involved in mitosis.

Genomic Risk Score for Melanoma in a Prospective Study of Older Individuals

BACKGROUND

Recent genome-wide association meta-analysis for melanoma doubled the number of previously identified variants. We assessed the performance of an updated polygenic risk score (PRS) in a population of older individuals, where melanoma incidence and cumulative ultraviolet radiation exposure is greatest.

METHODS

We assessed a PRS for cutaneous melanoma comprising 55 variants in a prospective study of 12 712 individuals in the ASPirin in Reducing Events in the Elderly Trial. We evaluated incident melanomas diagnosed during the trial and prevalent melanomas diagnosed preenrolment (self-reported). Multivariable models examined associations between PRS as a continuous variable (per SD) and categorical (low-risk [0%-20%], medium-risk [21%-80%], high-risk [81%-100%] groups) with incident melanoma. Logistic regression examined the association between PRS and prevalent melanoma.

RESULTS

At baseline, mean participant age was 75 years; 55.0% were female, and 528 (4.2%) had prevalent melanomas. During follow-up (median = 4.7 years), 120 (1.0%) incident cutaneous melanomas occurred, 98 of which were in participants with no history. PRS was associated with incident melanoma (hazard ratio = 1.46 per SD, 95% confidence interval [CI] = 1.20 to 1.77) and prevalent melanoma (odds ratio [OR] = 1.55 per SD, 95% CI = 1.42 to 1.69). Participants in the highest-risk PRS group had increased risk compared with the low-risk group for incident melanoma (OR = 2.51, 95% CI = 1.28 to 4.92) and prevalent melanoma (OR = 3.66, 95% CI = 2.69 to 5.05). When stratifying by sex, only males had an association between the PRS and incident melanoma, whereas both sexes had an association between the PRS and prevalent melanoma.

CONCLUSIONS

A genomic risk score is associated with melanoma risk in older individuals and may contribute to targeted surveillance.

Behind the Scene: Exploiting MC1R in Skin Cancer Risk and Prevention

Melanoma and non-melanoma skin cancers (NMSCs) are the most frequent cancers of the skin in white populations. An increased risk in the development of skin cancers has been associated with the combination of several environmental factors (i.e., ultraviolet exposure) and genetic background, including melanocortin-1 receptor (MC1R) status. In the last few years, advances in the diagnosis of skin cancers provided a great impact on clinical practice. Despite these advances, NMSCs are still the most common malignancy in humans and melanoma still shows a rising incidence and a poor prognosis when diagnosed at an advanced stage. Efforts are required to underlie the genetic and clinical heterogeneity of melanoma and NMSCs, leading to an optimization of the management of affected patients. The clinical implications of the impact of germline MC1R variants in melanoma and NMSCs' risk, together with the additional risk conferred by somatic mutations in other peculiar genes, as well as the role of MC1R screening in skin cancers' prevention will be addressed in the current review.

Identification of OCA2 as a novel locus for the co-morbidity of asthma-plus-eczema

BACKGROUND

Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co-morbidity of two of these diseases, have not been considered. This may partly explain missing heritability.

OBJECTIVE

To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema.

METHODS

We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co-morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status ("asthma-plus-eczema" vs. the presence of only one disease "asthma only or eczema only"). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses.

RESULTS

Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms.

CONCLUSION

Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes.

Genetic variants of EML1 and HIST1H4E in myeloid cell-related pathway genes independently predict cutaneous melanoma-specific survival

Both in vivo and in vitro evidence has supported a key role of myeloid cells in immune suppression in melanoma and in promoting melanocytic metastases. Some single-nucleotide polymorphisms (SNPs) have been shown to predict cutaneous melanoma-specific survival (CMSS), but the association between genetic variation in myeloid cell-related genes and cutaneous melanoma (CM) patient survival remains unknown.

METHODS

we investigated associations between SNPs in myeloid cell-related pathway genes and CMSS in a discovery dataset of 850 CM patients and replicated the findings in another dataset of 409 CM patients.

RESULTS

we identified two SNPs (EML1 rs10151787 A>G and HIST1H4E rs2069018 T>C) as independent prognostic factors for CMSS, with adjusted allelic hazards ratios of 1.56 (95% confidence interval =1.19-2.05, P=0.001) and 1.66 (1.22-2.26, P=0.001), respectively; so were their combined unfavorable alleles in a dose-response manner in both discovery and replication datasets (P trend<0.001 and 0.002, respectively). Additional functional analysis revealed that both EML1 rs10151787 G and HIST1H4E rs2069018 C alleles were associated with elevated mRNA expression levels in normal tissues.

CONCLUSIONS

Our findings suggest that EML1 rs10151787 A>G and HIST1H4E rs2069018 T>C are independent prognostic biomarkers for CMSS.

The α-melanocyte-stimulating hormone/melanocortin-1 receptor interaction: A driver of pleiotropic effects beyond pigmentation

Melanocortin-1 Receptor (MC1R), when stimulated by alpha-melanocyte-stimulating hormone (α-MSH), is a driver of eumelanogenesis. Brown/black eumelanin is an effective filter against ultraviolet radiation (UVR) and is a scavenger of free radicals. Several polymorphic variants of MC1R are frequent in red-head people. These polymorphisms reduce the ability of MC1R to promote eumelanogenesis after its activation and spontaneous pheomelanogenesis take place. Since pheomelanin can act as an endogenous photosensitizer, people carrying MC1R polymorphisms are more susceptible to skin cancer. Here, we summarize current knowledge on the biology of MC1R beyond its ability to drive eumelanogenesis. We analyze its capacity to cope with oxidative insult and consequent DNA damage. We describe its ability to transduce through different pathways. We start from the canonical pathway, the cAMP/protein kinase A (PKA) pathway mainly involved in promoting eumelanogenesis, and protection from oxidative damage, and we then move on to describe more recent knowledge concerning ERK pathways, phosphoinositide 3-kinase (PI3K) pathway/AKT, and α-MSH/Peroxisome proliferators activated receptor-γ (PPAR-γ) connection. We describe MC1R polymorphic variants associated with melanoma risk which represent an open window of clinical relevance.

MX2 mediates establishment of interferon response profile, regulates XAF1, and can sensitize melanoma cells to targeted therapy

MX2 is an interferon inducible gene that is mostly known for its antiviral activity. We have previously demonstrated that MX2 is also associated with the tumorigenesis process in melanoma. However, it remains unknown which molecular mechanisms are regulated by MX2 in response to interferon signaling in this disease. Here, we report that MX2 is necessary for the establishment of an interferon‐induced transcriptional profile partially through regulation of STAT1 phosphorylation and other interferon‐related downstream factors, including proapoptotic tumor suppressor XAF1. MX2 and XAF1 expression tightly correlate in both cultured melanoma cell lines and in patient‐derived primary and metastatic tumors, where they also are significantly related with survival. MX2 mediates IFN growth‐inhibitory signals in both XAF1 dependent and independent ways and in a cell type and context‐dependent manner. Higher MX2 expression renders melanoma cells more sensitive to targeted therapy drugs such as vemurafenib and trametinib; however, this effect is XAF1 independent. In summary, we uncovered a new mechanism in the complex regulation of interferon signaling in melanoma that can influence both survival and response to therapy.

Association of genetic variants of TMEM135 and PEX5 in the peroxisome pathway with cutaneous melanoma-specific survival

Background

Peroxisomes are ubiquitous and dynamic organelles that are involved in the metabolism of reactive oxygen species (ROS) and lipids. However, whether genetic variants in the peroxisome pathway genes are associated with survival in patients with melanoma has not been established. Therefore, our aim was to identify additional genetic variants in the peroxisome pathway that may provide new prognostic biomarkers for cutaneous melanoma (CM).

Methods

We assessed the associations between 8,397 common single-nucleotide polymorphisms (SNPs) in 88 peroxisome pathway genes and CM disease-specific survival (CMSS) in a two-stage analysis. For the discovery, we extracted the data from a published genome-wide association study from The University of Texas MD Anderson Cancer Center (MDACC). We then replicated the results in another dataset from the Nurse Health Study (NHS)/Health Professionals Follow-up Study (HPFS).

Results

Overall, 95 (11.1%) patients in the MDACC dataset and 48 (11.7%) patients in the NHS/HPFS dataset died of CM. We found 27 significant SNPs in the peroxisome pathway genes to be associated with CMSS in both datasets after multiple comparison correction using the Bayesian false-discovery probability method. In stepwise Cox proportional hazards regression analysis, with adjustment for other covariates and previously published SNPs in the MDACC dataset, we identified 2 independent SNPs (TMEM135 rs567403 C>G and PEX5 rs7969508 A>G) that predicted CMSS (P=0.003 and 0.031, respectively, in an additive genetic model). The expression quantitative trait loci analysis further revealed that the TMEM135 rs567403 GG and PEX5 rs7969508 GG genotypes were associated with increased and decreased levels of mRNA expression of their genes, respectively.

Conclusions

Once our findings are replicated by other investigators, these genetic variants may serve as novel biomarkers for the prediction of survival in patients with CM.

Genetically Raised Circulating Bilirubin Levels and Risk of Ten Cancers: A Mendelian Randomization Study

Bilirubin, an endogenous antioxidant, may play a protective role in cancer development. We applied two-sample Mendelian randomization to investigate whether genetically raised bilirubin levels are causally associated with the risk of ten cancers (pancreas, kidney, endometrium, ovary, breast, prostate, lung, Hodgkin's lymphoma, melanoma, and neuroblastoma). The number of cases and their matched controls of European descent ranged from 122,977 and 105,974 for breast cancer to 1200 and 6417 for Hodgkin's lymphoma, respectively. A total of 115 single-nucleotide polymorphisms (SNPs) associated (p < 5 × 10-8) with circulating total bilirubin, extracted from a genome-wide association study in the UK Biobank, were used as instrumental variables. One SNP (rs6431625) in the promoter region of the uridine-diphosphoglucuronate glucuronosyltransferase1A1 (UGT1A1) gene explained 16.9% and the remaining 114 SNPs (non-UGT1A1 SNPs) explained 3.1% of phenotypic variance in circulating bilirubin levels. A one-standarddeviation increment in circulating bilirubin (≈ 4.4 µmol/L), predicted by non-UGT1A1 SNPs, was inversely associated with risk of squamous cell lung cancer and Hodgkin's lymphoma (odds ratio (OR) 0.85, 95% confidence interval (CI) 0.73-0.99, P 0.04 and OR 0.64, 95% CI 0.42-0.99, p 0.04, respectively), which was confirmed after removing potential pleiotropic SNPs. In contrast, a positive association was observed with the risk of breast cancer after removing potential pleiotropic SNPs (OR 1.12, 95% CI 1.04-1.20, p 0.002). There was little evidence for robust associations with the other seven cancers investigated. Genetically raised bilirubin levels were inversely associated with risk of squamous cell lung cancer as well as Hodgkin's lymphoma and positively associated with risk of breast cancer. Further studies are required to investigate the utility of bilirubin as a low-cost clinical marker to improve risk prediction for certain cancers.

Functional annotation of melanoma risk loci identifies novel susceptibility genes

Genome-wide association study (GWAS)-identified single-nucleotide polymorphisms (SNPs) are tag SNPs located in both transcribed and non-coding regulatory DNA regions, rather than representing causal or functional variants for disease. To identify functional variants or genes for melanoma susceptibility, we used functional mapping and annotation (FUMA) to perform functional annotation of the summary statistics of 2541 significant melanoma risk SNPs (P < 5 × 10-8) identified by GWAS. The original GWAS melanoma study included 15 990 cases and 26 409 controls, representing the largest international meta-analysis of melanoma susceptibility. We prioritized 330 unique genes, including those in immune cytokine signaling pathways, from 19 loci through positional, expression quantitative trait locus, and chromatin interaction mapping. In comparison, only 38 melanoma-related genes were identified in the original meta-analysis. In addition to the well-known melanoma susceptibility genes confirmed in the meta-analysis (MC1R, CDKN2A, TERT, OCA2 and ARNT/SETDB1), we also identified additional novel genes using FUMA to map SNPs to genes. Through chromatin interaction mapping, we prioritized IFNA7, IFNA10, IFNA16, IFNA17, IFNA14, IFNA6, IFNA21, IFNA4, IFNE and IFNA5; these 10 most significant genes are all involved in immune system and cytokine signaling pathways. In the gene analysis, we identified 72 genes with a P < 2.5 × 10-6. The genes associated with melanoma risk were DEF8 (P = 1.09 × 10-57), DBNDD1 (P = 2.19 × 10-42), SPATA33 (P = 3.54 × 10-38) and MC1R (P = 1.04 × 10-36). In summary, this study identifies novel putative melanoma susceptibility genes and provides a guide for further experimental validation of functional variants and disease-related genes.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma.

Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness.

To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness.

The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry.

We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10⁻⁸) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

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

BACKGROUND

Folate metabolism plays an important role in DNA methylation and nucleic acid synthesis and thus may function as a regulatory factor in cancer development. Genome-wide association studies (GWASs) have identified some single-nucleotide polymorphisms (SNPs) associated with cutaneous melanoma-specific survival (CMSS), but no SNPs were found in genes involved in the folate metabolic pathway.

OBJECTIVES

To examine associations between SNPs in folate metabolic pathway genes and CMSS.

METHODS

We comprehensively evaluated 2645 (422 genotyped and 2223 imputed) common SNPs in folate metabolic pathway genes from a published GWAS of 858 patients from The University of Texas MD Anderson Cancer Center and performed the validation in another GWAS of 409 patients from the Nurses' Health Study and Health Professionals Follow-up Study, in which 95/858 (11·1%) and 48/409 (11·7%) patients died of cutaneous melanoma, respectively.

RESULTS

We identified two independent SNPs (MTHFD1 rs1950902 G>A and ALPL rs10917006 C>T) to be associated with CMSS in both datasets, and their meta-analysis yielded an allelic hazards ratio of 1·75 (95% confidence interval 1·32-2·32, P = 9·96 × 10-5 ) and 2·05 (1·39-3·01, P = 2·84 × 10-4 ), respectively. The genotype-phenotype correlation analyses provided additional support for the biological plausibility of these two variants' roles in tumour progression, suggesting that variation in SNP-related mRNA expression levels is likely to be the mechanism underlying the observed associations with CMSS.

CONCLUSIONS

Two possibly functional genetic variants, MTHFD1 rs1950902 and ALPL rs10917006, were likely to be independently or jointly associated with CMSS, which may add to personalized treatment in the future, once further validated. What is already known about this topic? Existing data show that survival rates vary among patients with melanoma with similar clinical characteristics; therefore, it is necessary to identify additional complementary biomarkers for melanoma-specific prognosis. A hypothesis-driven approach, by pooling the effects of single-nucleotide polymorphisms (SNPs) in a specific biological pathway as genetic risk scores, may provide a prognostic utility, and genetic variants of genes in folate metabolism have been reported to be associated with cancer risk. What does this study add? Two genetic variants in the folate metabolic pathway genes, MTHFD1 rs1950902 and ALPL rs10917006, are significantly associated with cutaneous melanoma-specific survival (CMSS). What is the translational message? The identification of genetic variants will make a risk-prediction model possible for CMSS. The SNPs in the folate metabolic pathway genes, once validated in larger studies, may be useful in the personalized management and treatment of patients with cutaneous melanoma.

Novel genetic variants of PIP5K1C and MVB12B of the endosome-related pathway predict cutaneous melanoma-specific survival

Endosomes regulate cell polarity, adhesion, signaling, immunity, and tumor progression, which may influence cancer outcomes. Here we evaluated associations between 36,068 genetic variants of 228 endosome-related pathway genes and cutaneous melanoma disease-specific survival (CMSS) using genotyping data from two previously published genome-wide association studies. The discovery dataset included 858 CM patients with 95 deaths from The University of Texas MD Anderson Cancer Center, and the replication dataset included 409 CM patients with 48 deaths from the Nurses' Health Study (NHS) and the Health Professionals Follow-up Study (HPFS). In multivariate Cox proportional hazards regression analysis, we found that two novel SNPs (PIP5K1C rs11666894 A>C and MVB12B rs12376285 C>T) predicted CMSS, with adjusted hazards ratios of 1.47 (95% confidence interval = 1.15-1.89 and P = 0.002) and 1.73 (1.30-2.31 and 0.0002), respectively. Combined analysis of risk genotypes of these two SNPs revealed a dose-dependent decrease in CMSS associated with an increased number of risk genotypes (P _trend = 0.0002). Subsequent expression quantitative trait loci (eQTL) analysis revealed that PIP5K1C rs11666894 was associated with mRNA expression levels in lymphoblastoid cell lines from 373 European descendants (P<0.0001) and that MVB12B rs12376285 was associated with mRNA expression levels in cultured fibroblasts from 605 European-Americans (P<0.0001). Our findings suggest that novel genetic variants of PIP5K1C and MVB12B in the endosome-related pathway genes may be promising prognostic biomarkers for CMSS, but these results need to be validated in future larger studies.

Germline and somatic albinism variants in amelanotic/hypomelanotic melanoma: Increased carriage of TYR and OCA2 variants

Amelanotic/hypomelanotic melanoma is a clinicopathologic subtype with absent or minimal melanin. This study assessed previously reported coding variants in albinism genes (TYR, OCA2, TYRP1, SLC45A2, SLC24A5, LRMDA) and common intronic, regulatory variants of OCA2 in individuals with amelanotic/hypomelanotic melanoma, pigmented melanoma cases and controls. Exome sequencing was available for 28 individuals with amelanotic/hypomelanotic melanoma and 303 individuals with pigmented melanoma, which were compared to whole exome data from 1144 Australian controls. Microarray genotyping was available for a further 17 amelanotic/hypomelanotic melanoma, 86 pigmented melanoma, 147 melanoma cases (pigmentation unknown) and 652 unaffected controls. Rare deleterious variants in TYR/OCA1 were more common in amelanotic/hypomelanotic melanoma cases than pigmented melanoma cases (set mixed model association tests P = 0.0088). The OCA2 hypomorphic allele p.V443I was more common in melanoma cases (1.8%) than controls (1.0%, X2 P = 0.02), and more so in amelanotic/hypomelanotic melanoma (4.4%, X2 P = 0.007). No amelanotic/hypomelanotic melanoma cases carried an eye and skin darkening haplotype of OCA2 (including rs7174027), present in 7.1% of pigmented melanoma cases (P = 0.0005) and 9.4% controls. Variants in TYR and OCA2 may play a role in amelanotic/hypomelanotic melanoma susceptibility. We suggest that somatic loss of function at these loci could contribute to the loss of tumor pigmentation, consistent with this we found a higher rate of somatic mutation in TYR/OCA2 in amelanotic/hypomelanotic melanoma vs pigmented melanoma samples (28.6% vs 3.0%; P = 0.021) from The Cancer Genome Atlas Skin Cutaneous Melanoma collection.

New Insights into the Role of Sphingolipid Metabolism in Melanoma

Cutaneous melanoma is a deadly skin cancer whose aggressiveness is directly linked to its metastatic potency. Despite remarkable breakthroughs in term of treatments with the emergence of targeted therapy and immunotherapy, the prognosis for metastatic patients remains uncertain mainly because of resistances. Better understanding the mechanisms responsible for melanoma progression is therefore essential to uncover new therapeutic targets. Interestingly, the sphingolipid metabolism is dysregulated in melanoma and is associated with melanoma progression and resistance to treatment. This review summarises the impact of the sphingolipid metabolism on melanoma from the initiation to metastatic dissemination with emphasis on melanoma plasticity, immune responses and resistance to treatments.

Genetic variants in TKT and DERA in the nicotinamide adenine dinucleotide phosphate pathway predict melanoma survival

BACKGROUND

Cutaneous melanoma (CM) is the most lethal type of skin cancers. Nicotinamide adenine dinucleotide phosphate (NADPH) plays an important role in anabolic reactions and tumorigenesis, but many genes are involved in the NADPH system.

METHODS

We used 10,912 single-nucleotide polymorphisms (SNPs) (2018 genotyped and 8894 imputed) in 134 NADPH-related genes from a genome-wide association study (GWAS) of 858 patients from The University of Texas MD Anderson Cancer Center (MDACC) in a single-locus analysis to predict CM survival. We then replicated the results in another GWAS data set of 409 patients from the Nurses' Health Study (NHS) and the Health Professionals Follow-up Study (HPFS).

RESULTS

There were 95 of 858 (11.1%) and 48 of 409 (11.7%) patients who died of CM, respectively. In multivariable Cox regression analyses, we identified two independent SNPs (TKT rs9864057 G > A and deoxyribose phosphate aldolase (DERA) rs12297652 A > G) to be significantly associated with CM-specific survival [hazards ratio (HR) of 1.52, 95% confidence interval (CI) = 1.18-1.96, P = 1.06 × 10^-3 and 1.51 (1.19-1.91, 5.89 × 10^-4)] in the meta-analysis, respectively. Furthermore, an increasing number of risk genotypes of these two SNPs was associated with a higher risk of death in the MDACC, the NHS/HPFS, and their combined data sets (P_trend<0.001, = 0.004 and <0.001, respectively). In the expression quantitative trait loci analysis, TKT rs9864057 G > A and DERA rs12297652 A > G were also significantly associated with higher mRNA expression levels in sun-exposed lower-leg skin (P = 0.043 and 0.006, respectively).

CONCLUSIONS

These results suggest that these two potentially functional SNPs may be valuable prognostic biomarkers for CM survival, but larger studies are needed to validate these findings.

HERC Ubiquitin Ligases in Cancer

HERC proteins are ubiquitin E3 ligases of the HECT family. The HERC subfamily is composed of six members classified by size into large (HERC1 and HERC2) and small (HERC3-HERC6). HERC family ubiquitin ligases regulate important cellular processes, such as neurodevelopment, DNA damage response, cell proliferation, cell migration, and immune responses. Accumulating evidence also shows that this family plays critical roles in cancer. In this review, we provide an integrated view of the role of these ligases in cancer, highlighting their bivalent functions as either oncogenes or tumor suppressors, depending on the tumor type. We include a discussion of both the molecular mechanisms involved and the potential therapeutic strategies.

A Genome-Wide Association Study Identifies Two Novel Susceptible Regions for Squamous Cell Carcinoma of the Head and Neck

To identify genetic variants for risk of squamous cell carcinoma of the head and neck (SCCHN), we conducted a two-phase genome-wide association study consisting of 7,858,089 SNPs in 2,171 cases and 4,493 controls of non-Hispanic white, of which, 434,839 typed and 7,423,250 imputed SNPs were used as the discovery. SNPs with P < 1 × 10-3 were further validated in the OncoArray study of oral and pharynx cancer (5,205 cases and 3,232 controls of European ancestry) from databases of Genotypes and Phenotypes. Meta-analysis of the discovery and replication studies identified one novel locus 6p22.1 (P = 2.96 × 10-9 for the leading rs259919) and two cancer susceptibility loci 6p21.32 (rs3135001, HLA-DQB1) and 6p21.33 (rs1265081, CCHCR1) associated with SCCHN risk. Further stratification by tumor site revealed four known cancer loci (5p15.33, 6p21.32, 6p21.33, and 2p23.1) associated with oral cavity cancer risk and oropharyngeal cancer risk, respectively. In addition, one novel locus 18q22.2 (P = 2.54 × 10-9 for the leading SNP rs142021700) was identified for hypopharynx and larynx cancer risk. For SNPs in those reported or novel loci, we also performed functional annotations by bioinformatics prediction and expression quantitative trait loci analysis. Collectively, our identification of four reported loci (2p23.1, 5p15.33, 6p21.32, and 6p21.33) and two novel loci (6p22.1 and 18q22.2) for SCCHN risk highlight the importance of human leukocyte antigen loci for oropharyngeal cancer risk, suggesting that immunologic mechanisms are implicated in the etiology of this subset of SCCHN. SIGNIFICANCE: Two novel risk loci for SCCHN in non-Hispanic white individuals highlight the importance of immunologic mechanism in the disease etiology.

An accurate and powerful method for copy number variation detection

MOTIVATION

Integration of multiple genetic sources for copy number variation detection (CNV) is a powerful approach to improve the identification of variants associated with complex traits. Although it has been shown that the widely used change point based methods can increase statistical power to identify variants, it remains challenging to effectively detect CNVs with weak signals due to the noisy nature of genotyping intensity data. We previously developed modSaRa, a normal mean-based model on a screening and ranking algorithm for copy number variation identification which presented desirable sensitivity with high computational efficiency. To boost statistical power for the identification of variants, here we present a novel improvement that integrates the relative allelic intensity with external information from empirical statistics with modeling, which we called modSaRa2.

RESULTS

Simulation studies illustrated that modSaRa2 markedly improved both sensitivity and specificity over existing methods for analyzing array-based data. The improvement in weak CNV signal detection is the most substantial, while it also simultaneously improves stability when CNV size varies. The application of the new method to a whole genome melanoma dataset identified novel candidate melanoma risk associated deletions on chromosome bands 1p22.2 and duplications on 6p22, 6q25 and 19p13 regions, which may facilitate the understanding of the possible roles of germline copy number variants in the etiology of melanoma.

AVAILABILITY AND IMPLEMENTATION

http://c2s2.yale.edu/software/modSaRa2 or https://github.com/FeifeiXiaoUSC/modSaRa2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic variants in PDSS1 and SLC16A6 of the ketone body metabolic pathway predict cutaneous melanoma-specific survival

A few single-nucleotide polymorphisms (SNPs) have been identified to be associated with cutaneous melanoma (CM) survival through genome-wide association studies, but stringent multiple testing corrections required for the hypothesis-free testing may have masked some true associations. Using a hypothesis-driven analysis approach, we sought to evaluate associations between SNPs in ketone body metabolic pathway genes and CM survival. We comprehensively assessed associations between 4196 (538 genotyped and 3658 imputed) common SNPs in 44 ketone body metabolic pathway genes and CM survival, using a dataset of 858 patients of a case-control study from The University of Texas M.D. Anderson Cancer Center as the discovery set and another dataset of 409 patients from the Nurses' Health Study and the Health Professionals Follow-up Study as the replication set. There were 95/858 (11.1%) and 48/409 (11.7%) patients who died of CM, respectively. We identified two independent SNPs (ie, PDSS1 rs12254548 G>C and SLC16A6 rs71387392 G>A) that were associated with CM survival, with allelic hazards ratios of 0.58 (95% confidence interval [CI] = 0.44-0.76, P = 9.00 × 10-5 ) and 1.98 (95% CI = 1.34-2.94, P = 6.30 × 10-4 ), respectively. Additionally, associations between genotypes of the SNPs and messenger RNA expression levels of their corresponding genes support the biologic plausibility of a role for these two variants in CM tumor progression and survival. Once validated by other larger studies, PDSS1 rs12254548 and SLC16A6 rs71387392 may be valuable biomarkers for CM survival.

Common germline variants of the human APOE gene modulate melanoma progression and survival

Common germline variants of the APOE gene are major risk modifiers of neurodegenerative and atherosclerotic diseases^1-3, but their effect on cancer outcome is poorly defined. Here we report that, in a reversal of their effect on Alzheimer's disease, the APOE4 and APOE2 variants confer favorable and poor outcomes in melanoma, respectively. Mice expressing the human APOE4 allele exhibited reduced melanoma progression and metastasis relative to APOE2 mice. APOE4 mice exhibited enhanced anti-tumor immune activation relative to APOE2 mice, and T cell depletion experiments showed that the effect of APOE genotype on melanoma progression was mediated by altered anti-tumor immunity. Consistently, patients with melanoma carrying the APOE4 variant experienced improved survival in comparison to carriers of APOE2. Notably, APOE4 mice also showed improved outcomes under PD1 immune checkpoint blockade relative to APOE2 mice, and patients carrying APOE4 experienced improved anti-PD1 immunotherapy survival after progression on frontline regimens. Finally, enhancing APOE expression via pharmacologic activation of liver X receptors, previously shown to boost anti-tumor immunity⁴, exhibited therapeutic efficacy in APOE4 mice but not in APOE2 mice. These findings demonstrate that pre-existing hereditary genetics can impact progression and survival outcomes of a future malignancy and warrant prospective investigation of APOE genotype as a biomarker for melanoma outcome and therapeutic response.

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.

MX 2 is a novel regulator of cell cycle in melanoma cells

MX2 protein is a dynamin-like GTPase2 that has recently been identified as an interferon-induced restriction factor of HIV-1 and other primate lentiviruses. A single nucleotide polymorphism (SNP), rs45430, in an intron of the MX2 gene, was previously reported as a novel melanoma susceptibility locus in genome-wide association studies. Functionally, however, it is still unclear whether and how MX2 contributes to melanoma susceptibility and tumorigenesis. Here, we show that MX2 is differentially expressed in melanoma tumors and cell lines, with most metastatic cell lines showing lower MX2 expression than primary melanoma cell lines and melanocytes. Furthermore, high expression of MX2 RNA in primary melanoma tumors is associated with better patient survival. Overexpression of MX2 reduces in vivo proliferation partially through inhibition of AKT activation, suggesting that it can act as a tumor suppressor in melanoma. However, we have also identified a subset of melanoma cell lines with high endogenous MX2 expression where downregulation of MX2 leads to reduced proliferation. In these cells, MX2 downregulation interfered with DNA replication and cell cycle processes. Collectively, our data for the first time show that MX2 is functionally involved in the regulation of melanoma proliferation but that its function is context-dependent.

Evolution of approaches to identify melanoma missing heritability

Introduction: Around 10% of melanoma patients have a positive family history of melanoma and/or related cancers. Although a germline pathogenic variant in a high-risk gene can be identified in up to 40% of these patients, the remaining part of melanoma heritability remains largely unexplained.Areas covered: The aim of this review is to provide an overview of the impact that new technologies and new research approaches had and are having on finding more efficient ways to unravel the missing heritability in melanoma.Expert opinion: High-throughput sequencing technologies have been crucial in increasing the number of genes/loci that might be implicated in melanoma predisposition. However, results from these approaches may have been inferior to the expectations, due to an increase in quantitative information which hasn't been followed at the same speed by an improvement of the methods to correctly interpret these data. Optimal approaches for improving our knowledge on melanoma heritability are currently based on segregation analysis coupled with functional assessment of candidate genes. An improvement of computational methods to infer genotype-phenotype correlations could help address the issue of missing heritability.

Novel Genetic Variants of ALG6 and GALNTL4 of the Glycosylation Pathway Predict Cutaneous Melanoma-Specific Survival

Because aberrant glycosylation is known to play a role in the progression of melanoma, we hypothesize that genetic variants of glycosylation pathway genes are associated with the survival of cutaneous melanoma (CM) patients. To test this hypothesis, we used a Cox proportional hazards regression model in a single-locus analysis to evaluate associations between 34,096 genetic variants of 227 glycosylation pathway genes and CM disease-specific survival (CMSS) using genotyping data from two previously published genome-wide association studies. The discovery dataset included 858 CM patients with 95 deaths from The University of Texas MD Anderson Cancer Center, and the replication dataset included 409 CM patients with 48 deaths from Harvard University nurse/physician cohorts. In the multivariable Cox regression analysis, we found that two novel single-nucleotide polymorphisms (SNPs) (ALG6 rs10889417 G>A and GALNTL4 rs12270446 G>C) predicted CMSS, with an adjusted hazards ratios of 0.60 (95% confidence interval = 0.44-0.83 and p = 0.002) and 0.66 (0.52-0.84 and 0.004), respectively. Subsequent expression quantitative trait loci (eQTL) analysis revealed that ALG6 rs10889417 was associated with mRNA expression levels in the cultured skin fibroblasts and whole blood cells and that GALNTL4 rs12270446 was associated with mRNA expression levels in the skin tissues (all p < 0.05). Our findings suggest that, once validated by other large patient cohorts, these two novel SNPs in the glycosylation pathway genes may be useful prognostic biomarkers for CMSS, likely through modulating their gene expression.

Regulation of the MDM2-p53 pathway by the ubiquitin ligase HERC2

The p53 tumor suppressor protein is a transcription factor that plays a prominent role in protecting cells from malignant transformation. Protein levels of p53 and its transcriptional activity are tightly regulated by the ubiquitin E3 ligase MDM2, the gene expression of which is transcriptionally regulated by p53 in a negative feedback loop. The p53 protein is transcriptionally active as a tetramer, and this oligomerization state is modulated by a complex formed by NEURL4 and the ubiquitin E3 ligase HERC2. Here, we report that MDM2 forms a complex with oligomeric p53, HERC2, and NEURL4. HERC2 knockdown results in a decline in MDM2 protein levels without affecting its protein stability, as it reduces its mRNA expression by inhibition of its promoter activation. DNA damage induced by bleomycin dissociates MDM2 from the p53/HERC2/NEURL4 complex and increases the phosphorylation and acetylation of oligomeric p53 bound to HERC2 and NEURL4. Moreover, the MDM2 promoter, which contains p53-response elements, competes with HERC2 for binding of oligomeric, phosphorylated and acetylated p53. We integrate these findings in a model showing the pivotal role of HERC2 in p53-MDM2 loop regulation. Altogether, these new insights in p53 pathway regulation are of great interest in cancer and may provide new therapeutic targets.

Genetic variants in ELOVL2 and HSD17B12 predict melanoma-specific survival

Fatty acids play a key role in cellular bioenergetics, membrane biosynthesis and intracellular signaling processes and thus may be involved in cancer development and progression. In the present study, we comprehensively assessed associations of 14,522 common single-nucleotide polymorphisms (SNPs) in 149 genes of the fatty-acid synthesis pathway with cutaneous melanoma disease-specific survival (CMSS). The dataset of 858 cutaneous melanoma (CM) patients from a published genome-wide association study (GWAS) by The University of Texas M.D. Anderson Cancer Center was used as the discovery dataset, and the identified significant SNPs were validated by a dataset of 409 CM patients from another GWAS from the Nurses' Health and Health Professionals Follow-up Studies. We found 40 noteworthy SNPs to be associated with CMSS in both discovery and validation datasets after multiple comparison correction by the false positive report probability method, because more than 85% of the SNPs were imputed. By performing functional prediction, linkage disequilibrium analysis, and stepwise Cox regression selection, we identified two independent SNPs of ELOVL2 rs3734398 T>C and HSD17B12 rs11037684 A>G that predicted CMSS, with an allelic hazards ratio of 0.66 (95% confidence interval = 0.51-0.84 and p = 8.34 × 10-4 ) and 2.29 (1.55-3.39 and p = 3.61 × 10-5 ), respectively. Finally, the ELOVL2 rs3734398 variant CC genotype was found to be associated with a significantly increased mRNA expression level. These SNPs may be potential markers for CM prognosis, if validated by additional larger and mechanistic studies.