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

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

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

OB-Folds and Genome Maintenance: Targeting Protein-DNA Interactions for Cancer Therapy

Genome stability and maintenance pathways along with their requisite proteins are critical for the accurate duplication of genetic material, mutation avoidance, and suppression of human diseases including cancer. Many of these proteins participate in these pathways by binding directly to DNA, and a subset employ oligonucleotide/oligosaccharide binding folds (OB-fold) to facilitate the protein-DNA interactions. OB-fold motifs allow for sequence independent binding to single-stranded DNA (ssDNA) and can serve to position specific proteins at specific DNA structures and then, via protein-protein interaction motifs, assemble the machinery to catalyze the replication, repair, or recombination of DNA. This review provides an overview of the OB-fold structural organization of some of the most relevant OB-fold containing proteins for oncology and drug discovery. We discuss their individual roles in DNA metabolism, progress toward drugging these motifs and their utility as potential cancer therapeutics. While protein-DNA interactions were initially thought to be undruggable, recent reports of success with molecules targeting OB-fold containing proteins suggest otherwise. The potential for the development of agents targeting OB-folds is in its infancy, but if successful, would expand the opportunities to impinge on genome stability and maintenance pathways for more effective cancer treatment.

Concomitance of a novel RMDN2-ALK fusion and an EML4-ALK fusion in a lung adenocarcinoma

The anaplastic lymphoma kinase (ALK) fusions/rearrangements in non-small cell lung cancer (NSCLC) act as oncogenic driver mutations. ALK tyrosine kinase inhibitors have anti-tumor activities in ALK-positive NSCLC. Although the EML4-ALK fusion is common in NSCLC, concomitance of an additional ALK fusion together with an EML4-ALK fusion is not common. Here, we present a lung adenocarcinoma with two ALK fusions, a novel RMDN2-ALK fusion accompanied by an EML4-ALK fusion, detected by a targeted next generation sequencing assay. The genomic translocation breakpoints of the RMDN2-ALK fusion were mapped to intron 2 for RMDN2 and exon 15 for ALK, and EML4-ALK breakpoints were mapped to intron 13 for EML4 and intron 19 for ALK. ALK break-apart FISH detected multiple ALK rearrangements, a gene fusion panel (NanoString) test confirmed the EML4-ALK fusion, and RNA-sequencing revealed two ALK fusions. The RMDN2 gene locates at the short arm of chromosome 2 between ALK and EML4 genes. The intact ALK kinase domain fused to RMDN2. Genome-wide copy number variants were found in multiple chromosome arms and the short arm of chromosome 2, suggestive of complex rearrangements. Further detailed analyses of breakpoints and copy number variants may shed light on mechanisms of their formation and pathogenesis in lung malignancies.

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.

Bladder cancer risk associated with family history of cancer

Twin studies suggest a familial aggregation of bladder cancer, but elements of this increased familial risk of bladder cancer are not well understood. To characterize familial risk of bladder cancer, we examined the relationship between family history of bladder and other types of cancer among first-degree relatives and risk of bladder cancer in 1193 bladder cancer cases and 1418 controls in a large population-based case-control study. Multivariate logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between family history of bladder cancer (defined as at least one first-degree family member with bladder cancer or a cancer of any other site). We also evaluated cancer aggregation of specific sites in family members. Participants with a first-degree relative with bladder cancer had nearly double the risk of bladder cancer (OR = 1.8, 95% CI 1.2-2.9) as those without a family history of bladder cancer. Risk was increased for having a sibling with bladder cancer (OR = 2.6, 95% CI 1.3-5.3) compared to no siblings with cancer. Bladder cancer risk was elevated when participants reported a first-degree relative with a history of female genital cancer (OR = 1.5, 95% CI 1.1-2.1), melanoma (OR = 1.9, 95% CI 1.02-3.6), and tobacco-associated cancer (OR = 1.3, 95% CI 1.06-1.6). These findings add to evidence of a familial predisposition to bladder cancer. Clarification of the aggregation of bladder cancer in families and with other cancer sites will be of interest as many loci and common polymorphisms related to bladder cancer have yet to be identified in large genomic studies.

Risk factors for melanoma by anatomical site: an evaluation of aetiological heterogeneity

BACKGROUND

Melanoma aetiology has been proposed to have two pathways, which are determined by naevi and type of sun exposure and related to the anatomical site where melanoma develops.

OBJECTIVES

We examined associations with melanoma by anatomical site for a comprehensive set of risk factors including pigmentary and naevus phenotypes, ultraviolet radiation exposure and polygenic risk.

METHODS

We analysed harmonized data from 2617 people with incident first invasive melanoma and 975 healthy controls recruited through two population-based case-control studies in Australia and the UK. Questionnaire data were collected by interview using a single protocol, and pathway-specific polygenic risk scores were derived from DNA samples. We estimated adjusted odds ratios using unconditional logistic regression that compared melanoma cases at each anatomical site with all controls.

RESULTS

When cases were compared with control participants, there were stronger associations for many naevi vs. no naevi for melanomas on the trunk, and upper and lower limbs than on the head and neck (P-heterogeneity < 0·001). Very fair skin (vs. olive/brown skin) was more weakly related to melanoma on the trunk than to melanomas at other sites (P-heterogeneity = 0·04). There was no significant difference by anatomical site for polygenic risk. Increased weekday sun exposure was positively associated with melanoma on the head and neck but not on other sites.

CONCLUSIONS

We found evidence of aetiological heterogeneity for melanoma, supporting the dual pathway hypothesis. These findings enhance understanding of risk factors for melanoma and can guide prevention and skin examination education and practices.

Genetic markers for characterization and prediction of prognosis of melanoma subtypes: a 2021 update

In this article we examined the most important genetic markers involved in melanoma susceptibility, initiation and progression, and their impact on the prognosis of the disease. Current knowledge in melanoma genetics identifies distinct pathways to the development of different melanoma subtypes characterized by specific clinico-pathological features and partially known genetic markers, modulated by high, low or absence of cumulative sun damage. The most prevalent somatic mutations are related to the activation of the MAPK pathway, which are classified into four major subtypes: BRAF mutant, NRAS mutant, NF1 mutant and triple wild type. Moreover, germinal mutations are also involved in the characterization and predictions of prognosis in melanoma. Currently, CDKN2A is seen as the main high-risk gene involved in melanoma susceptibility being mutated in around 20% of melanoma-prone families. Other high-risk susceptibility genes described include CDK4, POT1, BAP1, TERT promoter, ACD, and TERF2IP. Melanoma is one of the most genetically predisposed among all cancers in humans, and ultraviolet light from the sun is the main environmental factor. This genetic predisposition is starting to be understood, impacting not only on the risk of developing melanoma but also on the risk of developing other types of cancer, as well as on the prognosis of the disease.

Gene Discovery Using Twins

One of Nick's key early achievements at QIMR was to establish a twin study on melanoma risk factors. The Brisbane Twin Nevus Study (BTNS) had an initial focus on nevus (mole) count in adolescents but, reflecting Nick's broad interests, expanded in scope enormously over the decades. In the skin cancer arena, BTNS was essential to genetic discoveries in melanoma, eye color and pigmentation. Later studies amassed data on thousands of phenotypes, ranging from molecular phenotypes such as gene expression to studies where gene mapping findings in adolescents turned out to have translational potential in late-onset diseases. Nick's twin data have formed the basis for an enormous range of discoveries, with Nick and his colleagues continuing to capitalize on these data.

Body mass index and height and risk of cutaneous melanoma: Mendelian randomization analyses

BACKGROUND

Height and body mass index (BMI) have both been positively associated with melanoma risk, although findings for BMI have been less consistent than height. It remains unclear, however, whether these associations reflect causality or are due to residual confounding by environmental and lifestyle risk factors. We re-evaluated these associations using a two-sample Mendelian randomization (MR) approach.

METHODS

We identified single nucleotide polymorphisms (SNPs) for BMI and height from separate genome-wide association study (GWAS) meta-analyses. We obtained melanoma SNPs from the most recent melanoma GWAS meta-analysis comprising 12 874 cases and 23 203 controls. We used the inverse variance-weighted estimator to derive separate causal risk estimates across all SNP instruments for BMI and height.

RESULTS

Based on the combined estimate derived from 730 SNPs for BMI, we found no evidence of an association between genetically predicted BMI and melanoma [odds ratio (OR) per one standard deviation (1 SD) (4.6 kg/m2) increase in BMI 1.00, 95% confidence interval (CI): 0.91-1.11]. In contrast, we observed a positive association between genetically-predicted height (derived from a pooled estimate of 3290 SNPs) and melanoma risk [OR 1.08, 95% CI: 1.02-1.13, per 1 SD (9.27 cm) increase in height]. Sensitivity analyses using two alternative MR methods yielded similar results.

CONCLUSIONS

These findings provide no evidence for a causal association between higher BMI and melanoma, but support the notion that height is causally associated with melanoma risk. Mechanisms through which height influences melanoma risk remain unclear, and it remains possible that the effect could be mediated through diverse pathways including growth factors and even socioeconomic status.

Genetic and expression variations of cell cycle pathway genes in brain tumor patients

The present study was designed to determine the association between the genetic polymorphisms/expression variations of RB1 and CCND1 genes and brain tumor risk. For this purpose, 250 blood samples of brain tumor patients along with 250 controls (cohort I) and 96 brain tumor tissues (cohort II) with adjacent control section were collected. Mutation analysis of RB1 (rs137853294, rs121913300) and CCND1 (rs614367, rs498136) genes was performed using ARMS-PCR followed by sequencing, and expression analysis was performed using real-time PCR and immunohistochemistry. The results showed homozygous mutant genotype of RB1 gene polymorphism, rs121913300 (P=0.003) and CCND1 gene polymorphism rs614367 (P=0.01) were associated significantly with brain tumor risk. Moreover, significant down-regulation of RB1 (P=0.005) and up-regulation of CCND1 (P=0.0001) gene was observed in brain tumor sections vs controls. Spearman correlation showed significant negative correlation between RB1 vs proliferation marker, Ki-67 (r = -0.291*, P<0.05) in brain tumors. Expression levels of selected genes were also assessed at protein level using immunohistochemical analysis (IHC) and signification down-regulation of RB1 (P=0.0001) and up-regulation of CCND1 (P=0.0001) was observed in brain tumor compared with control sections. In conclusion, it is suggested that polymorphisms/expression variations of RB1 and CCND1 genes may be associated with increased risk of brain tumor.

On Naevi and Melanomas: Two Sides of the Same Coin?

Benign naevi are closely linked to melanoma, as risk factors, simulators, or sites of melanoma formation. There is a heavy genetic overlap between the two lesions, a shared environmental influence of ultraviolet radiation, and many similar cellular features, yet naevi remain locally situated while melanomas spread from their primary site and may progress systemically to distal organs. Untangling the overlapping contributors and predictors of naevi and melanoma is an ongoing area of research and should eventually lead to more personalized prevention and treatment strategies, through the development of melanoma risk stratification tools and early detection of evolving melanomas. This will be achieved through a range of complementary strategies: risk-adjusted primary prevention counseling; the use of lesion imaging technologies such as sequential 3D total body photography and consumer-performed lesion imaging; artificial intelligence deep phenotyping and clinical assistance; a better understanding of genetic drivers of malignancy, risk variants, clinical genetics, and polygenic effects; and the interplay between genetics, phenotype and the environment.

Molecular Cytogenetic Characterization of the Murine Melanoma Cell Lines S91 Clone M3 and B16-F1 with Variant B16-4A5

Melanoma is considered to be one of the most aggressive human tumors. Thus, early molecular diagnosis with risk factor stratification could be an efficacious strategy to increase the survival rates in affected patients. Murine cell lines B16-F1, B16-4A5, and S91 clone M3 are the ones most commonly applied in melanoma research. However, genetic peculiarities of these 3 cell lines have not been studied in detail before. Here, we closed this gap by molecular cytogenetic and array-comparative genomic hybridization studies and the translation of the characterized imbalances into the human genome. This study revealed severely rearranged karyotypes with in parts similar imbalances for all 3 cell lines. Interestingly, they involve genes known to play major roles in human melanoma. These are specifically the oncogenes and tumor suppressor genes, being associated with aggressive forms of melanoma. B16-F1, B16-4A5, and S91 clone M3 revealed aberrations which were similarly observed in human eye and skin but not in human uveal melanoma. Thus, they can be considered as model systems for advanced eye and skin melanoma.

A comprehensive re-assessment of the association between vitamin D and cancer susceptibility using Mendelian randomization

Previous Mendelian randomization (MR) studies on 25-hydroxyvitamin D (25(OH)D) and cancer have typically adopted a handful of variants and found no relationship between 25(OH)D and cancer; however, issues of horizontal pleiotropy cannot be reliably addressed. Using a larger set of variants associated with 25(OH)D (74 SNPs, up from 6 previously), we perform a unified MR analysis to re-evaluate the relationship between 25(OH)D and ten cancers. Our findings are broadly consistent with previous MR studies indicating no relationship, apart from ovarian cancers (OR 0.89; 95% C.I: 0.82 to 0.96 per 1 SD change in 25(OH)D concentration) and basal cell carcinoma (OR 1.16; 95% C.I.: 1.04 to 1.28). However, after adjustment for pigmentation related variables in a multivariable MR framework, the BCC findings were attenuated. Here we report that lower 25(OH)D is unlikely to be a causal risk factor for most cancers, with our study providing more precise confidence intervals than previously possible.

Human CTC1 promotes TopBP1 stability and CHK1 phosphorylation in response to telomere dysfunction and global replication stress

CST (CTC1-STN1-TEN1) is a heterotrimeric, RPA-like complex that binds to single-stranded DNA (ssDNA) and functions in the replication of telomeric and non-telomeric DNA. Previous studies demonstrated that deletion of CTC1 results in decreased cell proliferation and telomere DNA damage signaling. However, a detailed analysis of the consequences of conditional CTC1 knockout (KO) has not been fully elucidated. Here, we investigated the effects of CTC1 KO on cell cycle progression, genome-wide replication and activation of the DNA damage response. Consistent with previous findings, we demonstrate that CTC1 KO results in decreased cell proliferation, G2 arrest and RPA-bound telomeric ssDNA. However, despite the increased levels of telomeric RPA-ssDNA, global ATR-dependent CHK1 and p53 phosphorylation was not detected in CTC1 KO cells. Nevertheless, we show that RPA-ssDNA does activate ATR, leading to the phosphorylation of RPA and autophosphorylation of ATR. Further analysis determined that inactivation of ATR, but not CHK1 or ATM, suppressed the accumulation of G2 arrested cells and phosphorylated RPA following CTC1 removal. These results suggest that ATR is localized and active at telomeres but is unable to elicit a global checkpoint response through CHK1. Furthermore, CTC1 KO inhibited CHK1 phosphorylation following hydroxyurea-induced replication stress. Additional studies revealed that this suppression of CHK1 phosphorylation, following replication stress, is caused by decreased levels of the ATR activator TopBP1. Overall, our results identify CST as a novel regulator of the ATR-CHK1 pathway.

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.

A penalized regression framework for building polygenic risk models based on summary statistics from genome-wide association studies and incorporating external information

Large-scale genome-wide association (GWAS) studies provide opportunities for developing genetic risk prediction models that have the potential to improve disease prevention, intervention or treatment. The key step is to develop polygenic risk score (PRS) models with high predictive performance for a given disease, which typically requires a large training data set for selecting truly associated single nucleotide polymorphisms (SNPs) and estimating effect sizes accurately. Here, we develop a comprehensive penalized regression for fitting l ₁ regularized regression models to GWAS summary statistics. We propose incorporating Pleiotropy and ANnotation information into PRS (PANPRS) development through suitable formulation of penalty functions and associated tuning parameters. Extensive simulations show that PANPRS performs equally well or better than existing PRS methods when no functional annotation or pleiotropy is incorporated. When functional annotation data and pleiotropy are informative, PANPRS substantially outperforms existing PRS methods in simulations. Finally, we applied our methods to build PRS for type 2 diabetes and melanoma and found that incorporating relevant functional annotations and GWAS of genetically related traits improved prediction of these two complex diseases.

Multiplex melanoma families are enriched for polygenic risk

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

Identification of Germline Mutations in Melanoma Patients with Early Onset, Double Primary Tumors, or Family Cancer History by NGS Analysis of 217 Genes

Cutaneous melanoma is the deadliest skin malignity with a rising prevalence worldwide. Patients carrying germline mutations in melanoma-susceptibility genes face an increased risk of melanoma and other cancers. To assess the spectrum of germline variants, we analyzed 264 Czech melanoma patients indicated for testing due to early melanoma (at <25 years) or the presence of multiple primary melanoma/melanoma and other cancer in their personal and/or family history. All patients were analyzed by panel next-generation sequencing targeting 217 genes in four groups: high-to-moderate melanoma risk genes, low melanoma risk genes, cancer syndrome genes, and other genes with an uncertain melanoma risk. Population frequencies were assessed in 1479 population-matched controls. Selected POT1 and CHEK2 variants were characterized by functional assays. Mutations in clinically relevant genes were significantly more frequent in melanoma patients than in controls (31/264; 11.7% vs. 58/1479; 3.9%; p = 2.0 × 10⁻⁶). A total of 9 patients (3.4%) carried mutations in high-to-moderate melanoma risk genes (CDKN2A, POT1, ACD) and 22 (8.3%) patients in other cancer syndrome genes (NBN, BRCA1/2, CHEK2, ATM, WRN, RB1). Mutations in high-to-moderate melanoma risk genes (OR = 52.2; 95%CI 6.6–413.1; p = 3.2 × 10⁻⁷) and in other cancer syndrome genes (OR = 2.3; 95%CI 1.4–3.8; p = 0.003) were significantly associated with melanoma risk. We found an increased potential to carry these mutations (OR = 2.9; 95%CI 1.2–6.8) in patients with double primary melanoma, melanoma and other primary cancer, but not in patients with early age at onset. The analysis revealed affected genes in Czech melanoma patients and identified individuals who may benefit from genetic testing and future surveillance management of mutation carriers.

Current landscape of Oculocutaneous Albinism in Japan

Oculocutaneous albinism (OCA), which is roughly divided into non-syndromic and syndromic OCA, is a group of autosomal recessive disorders caused by mutations in genes associated with pigmentation. Patients with OCA have hypopigmentation and ocular manifestations such as photophobia, amblyopia, and nystagmus. Hermansky-Pudlak syndrome (HPS), the most common syndromic OCA, is characterized by the additional features of a bleeding tendency and other critical systemic comorbidities such as pulmonary fibrosis and immunodeficiency. NGS-based gene analyses have identified several new causative genes for OCA and have detected rare subtypes of OCA with high accuracy including Japanese patients. In our survey of 190 Japanese OCA patients/families, OCA4 is the most common subtype (25.3%) followed by OCA1 (20.0%), HPS1 (14.7%), and OCA2 (8.4%). Similar to the A481T variant in OCA2, which is associated with a mild form of OCA2 and skin color variation, the c.-492_489delAATG variant located in the promoter region of SLC45A2 has been uniquely identified in Japanese patients with a mild form of OCA4. Further, rare OCA subtypes, including OCA3, HPS2, HPS3, HPS4, HPS5, HPS6, and HPS9, have also been identified in Japanese patients. The clinical characteristics and underlying molecular mechanisms of each subtype of OCA are concisely summarized in this review.

The aryl hydrocarbon receptor promotes inflammation-induced dedifferentiation and systemic metastatic spread of melanoma cells

The aryl hydrocarbon receptor (AHR) is a ligand binding-transcription factor of the basic helix-loop-helix family regulating multiple cellular functions such as differentiation, cell cycle, apoptosis, and inflammatory reactions. In neoplastic diseases, the AHR has been described to modulate proliferation and differentiation in dichotomous ways, either inhibiting or augmenting the growth of tumors. The precise role of AHR in melanoma is mostly unknown. Here, we report a functional effect of AHR activation on inflammation-induced melanoma cell dedifferentiation and the development of lung metastases in a mouse model. Via in silico analyses of "The Cancer Genome Atlas" human melanoma cohort, we detected a correlation between AHR expression levels and a dedifferentiated melanoma cell phenotype with an invasive gene signature, which we were able to functionally recapitulate in a panel of human melanoma cell lines. Both human and mouse melanoma cell lines upregulated AHR expression after inflammatory stimulation with tumor necrosis factor-α (TNF-α). Activation of AHR in human and mouse melanoma cell lines with the endogenous ligand formylindolo(3,2-b)carbazole (FICZ) promoted inflammation-induced dedifferentiation in vitro. Importantly, mouse melanoma cells with CRISPR/Cas9-mediated disruption of the AHR gene showed impaired in vivo tumor growth after transplantation in the skin as well as decreased numbers of spontaneous lung metastases. Taken together, our results demonstrate a functional role for AHR expression in melanoma development and metastatic progression. This provides a scientific basis for future experiments that further dissect the underlying molecular mechanisms and assess the potential for AHR inhibition as part of multimodal melanoma treatment strategies.

Inherited variations in human pigmentation-related genes modulate cutaneous melanoma risk and clinicopathological features in Brazilian population

Ultraviolet light exposure and cutaneous pigmentation are important host risk factors for cutaneous melanoma (CM), and it is well known that inherited ability to produce melanin varies in humans. The study aimed to identify single-nucleotide variants (SNVs) on pigmentation-related genes with importance in risk and clinicopathological aspects of CM. The study was conducted in two stages. In stage 1, 103 CM patients and 103 controls were analyzed using Genome-Wide Human SNV Arrays in order to identify SNVs in pigmentation-related genes, and the most important SNVs were selected for data validation in stage 2 by real-time polymerase-chain reaction in 247 CM patients and 280 controls. ADCY3 c.675+9196T>G, CREB1 c.303+373G>A, and MITF c.938-325G>A were selected for data validation among 74 SNVs. Individuals with CREB1 GA or AA genotype and allele "A" were under 1.79 and 1.47-fold increased risks of CM than others, respectively. Excesses of CREB1 AA and MITF AA genotype were seen in patients with tumors at Clark levels III to V (27.8% versus 13.7%) and at III or IV stages (46.1% versus 24.9%) compared to others, respectively. When compared to others, patients with ADCY3 TT had 1.89 more chances of presenting CM progression, and those with MITF GA or AA had 2.20 more chances of evolving to death by CM. Our data provide, for the first time, preliminary evidence that inherited abnormalities in ADCY3, CREB1, and MITF pigmentation-related genes, not only can increase the risk to CM, but also influence CM patients' clinicopathological features.

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.

Assessment of polygenic architecture and risk prediction based on common variants across fourteen cancers

Genome-wide association studies (GWAS) have led to the identification of hundreds of susceptibility loci across cancers, but the impact of further studies remains uncertain. Here we analyse summary-level data from GWAS of European ancestry across fourteen cancer sites to estimate the number of common susceptibility variants (polygenicity) and underlying effect-size distribution. All cancers show a high degree of polygenicity, involving at a minimum of thousands of loci. We project that sample sizes required to explain 80% of GWAS heritability vary from 60,000 cases for testicular to over 1,000,000 cases for lung cancer. The maximum relative risk achievable for subjects at the 99th risk percentile of underlying polygenic risk scores (PRS), compared to average risk, ranges from 12 for testicular to 2.5 for ovarian cancer. We show that PRS have potential for risk stratification for cancers of breast, colon and prostate, but less so for others because of modest heritability and lower incidence.

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

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

There are more than 20 known melanoma susceptibility genes. Here, using a massively parallel reporter assay, the authors identify risk-associated variants that alter gene transcription, and demonstrate that expression of one such gene, MX2, leads to the promotion of melanoma in a zebrafish model.

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.

Emerging functions and clinical prospects of connexins and pannexins in melanoma

Cellular communication through gap junctions and hemichannels formed by connexins and through channels made by pannexins allows for metabolic cooperation and control of cellular activity and signalling. These channel proteins have been described to be tumour suppressors that regulate features such as cell death, proliferation and differentiation. However, they display cancer type-dependent and stage-dependent functions and may facilitate tumour progression through junctional and non-junctional pathways. The accumulated knowledge and emerging strategies to target connexins and pannexins are providing novel clinical opportunities for the treatment of cancer. Here, we provide an updated overview of the role of connexins and pannexins in malignant melanoma. We discuss how targeting of these channel proteins may be used to potentiate antitumour effects in therapeutic settings, including through improved immune-mediated tumour elimination.

Neural crest-derived tumor neuroblastoma and melanoma share 1p13.2 as susceptibility locus that shows a long-range interaction with the SLC16A1 gene

Neuroblastoma (NB) and malignant cutaneous melanoma (CMM) are neural crest cells (NCC)-derived tumors and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association studies (GWAS). We took a three-staged approach to conduct cross-disease meta-analysis of GWAS for NB and CMM (2101 NB cases and 4202 controls; 12 874 CMM cases and 23 203 controls) to identify shared loci. Findings were replicated in 1403 NB cases and 1403 controls of European ancestry and in 636 NB, 508 CMM cases and 2066 controls of Italian origin. We found a cross-association at locus 1p13.2 (rs2153977, odds ratio = 0.91, P = 5.36 × 10-8). We also detected a suggestive (P < 10-7) NB-CMM cross-association at 2q37.1 with opposite effect on cancer risk. Pathway analysis of 110 NB-CMM risk loci with P < 10-4 demonstrated enrichment of biological processes such as cell migration, cell cycle, metabolism and immune response, which are essential of human NCC development, underlying both tumors. In vitro and in silico analyses indicated that the rs2153977-T protective allele, located in an NB and CMM enhancer, decreased expression of SLC16A1 via long-range loop formation and altered a T-box protein binding site. Upon depletion of SLC16A1, we observed a decrease of cellular proliferation and invasion in both NB and CMM cell lines, suggesting its role as oncogene. This is the largest study to date examining pleiotropy across two NC cell-derived tumors identifying 1p13.2 as common susceptibility locus for NB and CMM risk. We demonstrate that combining genome-wide association studies results across cancers with same origins can identify new loci common to neuroblastoma and melanoma arising from tissues which originate from neural crest cells. Our results also show 1p13.2 confer risk to neuroblastoma and melanoma by regulating SLC16A1.

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.

MC1R variants and associations with pigmentation characteristics and genetic ancestry in a Hispanic, predominately Puerto Rican, population

Skin cancer risk information based on melanocortin-1 receptor (MC1R) variants could inform prevention and screening recommendations for Hispanics, but limited evidence exists on the impact of MC1R variants in Hispanic populations. We studied Hispanic subjects, predominately of Puerto Rican heritage, from Tampa, Florida, US, and Ponce, PR. Blood or saliva samples were collected by prospective recruitment or retrieved from biobanks for genotyping of MC1R variants and ancestry informative markers. Participant demographic and self-reported phenotypic information was collected via biobank records or questionnaires. We determined associations of MC1R genetic risk categories and phenotypic variables and genetic ancestry. Over half of participants carried MC1R variants known to increase risk of skin cancer, and there was diversity in the observed variants across sample populations. Associations between MC1R genetic risk groups and some pigmentation characteristics were identified. Among Puerto Ricans, the proportion of participants carrying MC1R variants imparting elevated skin cancer risk was consistent across quartiles of European, African, and Native American genetic ancestry. These findings demonstrate that MC1R variants are important for pigmentation characteristics in Hispanics and that carriage of high risk MC1R alleles occurs even among Hispanics with stronger African or Native American genetic ancestry.

Genetic Variation in Steroid and Xenobiotic Metabolizing Pathways and Enterolactone Excretion Before and After Flaxseed Intervention in African American and European American Women

BACKGROUND

Metabolism and excretion of the phytoestrogen enterolactone (ENL), which has been associated with breast cancer risk, may be affected by variation in steroid hormone and xenobiotic-metabolizing genes.

METHODS

We conducted a randomized, crossover flaxseed intervention study in 252 healthy, postmenopausal women [137 European ancestry (EA) and 115 African ancestry (AA)] from western New York. Participants were randomly assigned to maintain usual diet or consume 10 g/day ground flaxseed for 6 weeks. After a 2-month washout period, participants crossed over to the other diet condition for an additional 6 weeks. Urinary ENL excretion was measured by gas chromatography-mass spectrometry and 70 polymorphisms in 29 genes related to steroid hormone and xenobiotic metabolism were genotyped. Mixed additive genetic models were constructed to examine association of genetic variation with urinary ENL excretion at baseline and after the flaxseed intervention.

RESULTS

SNPs in several genes were nominally (P < 0.05) associated with ENL excretion at baseline and/or after intervention: ESR1, CYP1B1, COMT, CYP3A5, ARPC1A, BCL2L11, SHBG, SLCO1B1, and ZKSCAN5. A greater number of SNPs were associated among AA women than among EA women, and no SNPs were associated in both races. No SNP-ENL associations were statistically significant after correction for multiple comparisons.

CONCLUSIONS

Variation in several genes related to steroid hormone metabolism was associated with lignan excretion at baseline and/or after flaxseed intervention among postmenopausal women.

IMPACT

These findings may contribute to our understanding of the differences observed in urinary ENL excretion among AA and EA women and thus hormone-related breast cancer risk.

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.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci

The keratinocyte cancers (KC), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common cancers in fair-skinned people. KC treatment represents the second highest cancer healthcare expenditure in Australia. Increasing our understanding of the genetic architecture of KC may provide new avenues for prevention and treatment. We first conducted a series of genome-wide association studies (GWAS) of KC across three European ancestry datasets from Australia, Europe and USA, and used linkage disequilibrium (LD) Score regression (LDSC) to estimate their pairwise genetic correlations. We employed a multiple-trait approach to map genes across the combined set of KC GWAS (total N = 47 742 cases, 634 413 controls). We also performed meta-analyses of BCC and SCC separately to identify trait specific loci. We found substantial genetic correlations (generally 0.5–1) between BCC and SCC suggesting overlapping genetic risk variants. The multiple trait combined KC GWAS identified 63 independent genome-wide significant loci, 29 of which were novel. Individual separate meta-analyses of BCC and SCC identified an additional 13 novel loci not found in the combined KC analysis. Three new loci were implicated using gene-based tests. New loci included common variants in BRCA2 (distinct to known rare high penetrance cancer risk variants), and in CTLA4, a target of immunotherapy in melanoma. We found shared and trait specific genetic contributions to BCC and SCC. Considering both, we identified a total of 79 independent risk loci, 45 of which are novel.

Large miRNA survival analysis reveals a prognostic four-biomarker signature for triple negative breast cancer

Triple negative breast cancer (TNBC) is currently the only major breast tumor subtype without effective targeted therapy and, as a consequence, usually presents a poor outcome. Due to its more aggressive phenotype, there is an urgent clinical need to identify novel biomarkers that discriminate individuals with poor prognosis. We hypothesize that miRNAs can be used to this end because they are involved in the initiation and progression of tumors by altering the expression of their target genes. To identify a prognostic biomarker in TNBC, we analyzed the miRNA expression of a cohort composed of 185 patients diagnosed with TNBC using penalized Cox regression models. We identified a four-biomarker signature based on miR-221, miR-1305, miR-4708, and RMDN2 expression levels that allowed for the subdivision of TNBC into high- or low-risk groups (Hazard Ratio – HR = 0.32; 95% Confidence Interval - CI = 0.11–0.91; p = 0.03) and are also statistically associated with survival outcome in subgroups of postmenopausal status (HR = 0.19; 95% CI = 0.04–0.90; p= 0.016), node negative status (HR = 0.12; 95% CI = 0.01–1.04; p = 0.026), and tumors larger than 2cm (HR = 0.21; 95% CI = 0.05–0.81; p = 0.021). This four-biomarker signature was significantly associated with TNBC as an independent prognostic factor for survival.

Cutaneous Photoprotection: A Review of the Current Status and Evolving Strategies

Ultraviolet radiation (UVR) exposure is well established as the major environmental risk factor for the development of melanoma, cutaneous squamous cell carcinoma (cSCC), and basal cell carcinoma (BCC). Additional risk factors including genetic mutations, other environmental agents, and immune status are important in modulating the effects of UVR. Dermatologists advocate a multi-pronged approach to minimizing UVR exposure including lifestyle modifications, UVR protective clothing, and topically applied sun-protective products, i.e. sunscreen. New Federal Drug Administration (FDA) regulations on sunscreen have brought certain long-standing ingredients in sunscreen products under scrutiny. The FDA's proposed rule for over the counter (OTC) monograph states that the inorganic sunscreens, zinc oxide and titanium dioxide, were found to be "generally recognized as safe and effective," but cite insufficient evidence to grant organic sunscreens the same designation. This proposed rule by the FDA and our increasing understanding of multifactorial mechanisms of UVR damage are an impetus for innovation and advances in sun protective technology. A complete set of strategies designed to limit the risk of UV-induced skin cell malignant transformation and tumor development must address the fuller consideration of genetic, environmental, and immune factors that cooperatively drive cutaneous carcinogenesis. Recent advances in our understanding of the biochemical processes underpinning UVR associated cutaneous cellular damage, genotoxicity, and clonal expansion provide investigators with a spectrum of opportunities for technologic innovation in the prevention of skin cancer. Strategies to improve upon current topical sunscreen formulations have strived for broader UVR spectral coverage, more favorable aesthetics, increased adherence, and minimal penetration into the living epidermis. In addition to improved sunscreens, future topical therapies may target processes within the epidermis that contribute to carcinogenesis. These include reactive species quenching, delivery of DNA repair enzymes, and targeting of cytokines essential to the proliferation of mutant keratinocytes.

The interplay of sun damage and genetic risk in Australian multiple and single primary melanoma cases and controls

BACKGROUND

Skin phenotype, host genotype and ultraviolet (UV) damage play a role in the development of melanoma.

OBJECTIVES

To ascertain whether the level of UV damage at the site of melanomas was associated with genetic polymorphisms.

METHODS

Deep phenotyping was performed on 1244 individuals; 281 with multiple primary melanomas (MPMs), 304 with single primary melanoma (SPM) and 659 convenience controls. Genotype data was generated using the Illumina CoreExome microarray platform, assaying over 500 000 single-nucleotide polymorphisms. A subset of variants were combined to assess a polygenic risk score (PRS) for melanoma.

RESULTS

Most MPM cases were diagnosed in patients aged > 40 years, in sites with visible chronic UV damage. Women and those diagnosed at age ≤ 40 years were less likely to have perilesional UV damage. Patients with MPM had higher frequencies of MITF E318K, MC1R R-alleles and the ASIP risk haplotype. Individuals who had melanoma in a visibly UV-damaged site were more likely to carry MC1R rs75570604 [odds ratio (OR) 2·5], 9q31.2 rs10816595 (OR 1·4) and MTAP rs869329 (OR 1·4). These same alleles were more common in patients with MPM who were diagnosed at age ≤ 40 years. The mean PRS was significantly higher in MPM than in SPM and controls. Naevus count was comparable in early-onset MPM cases and those diagnosed at age > 40 years.

CONCLUSIONS

Our cohort demonstrated higher frequencies of previously reported alleles associated with melanoma. MPM melanomas more commonly occur in UV-damaged areas, and these individuals are more likely to carry MC1R red hair colour alleles. Awareness of the interplay of genetic vulnerability with UV damage can stratify risk and guide recommendations for melanoma screening. What's already known about this topic? Skin phenotype, host genotype and ultraviolet (UV) damage all play a role in melanoma development. One of the main risk factors is a personal history of melanoma; second and subsequent primary melanomas account for over 20% of all melanomas registered in Queensland. Multiple loci are associated with melanoma risk, including many low-penetrance loci, which may have a cumulatively significant risk. Population-wide screening programmes for melanoma are not yet economically viable. What does this study add? Patients diagnosed with melanoma at age ≤ 40 years were more likely than older patients to have melanomas in non-UV-damaged sites. Patients with multiple melanomas had higher frequencies of MITF E318K, MC1R R-alleles, and the ASIP extended risk haplotype than patients with single melanoma. CDKN2A, MC1R and MTAP variants were more frequent in patients who developed melanomas at a younger age, but also in those whose melanomas were all on visibly UV-damaged sites. What is the translational message? Incorporating these genetic findings into the known risk factors of skin phenotype and visible UV damage may allow for a more customized and economically feasible approach to early detection of melanoma, particularly in younger patients. Plain language summary available online.

Overlapping genetic architecture between Parkinson disease and melanoma

Epidemiologic studies have reported inconsistent results regarding an association between Parkinson disease (PD) and cutaneous melanoma (melanoma). Identifying shared genetic architecture between these diseases can support epidemiologic findings and identify common risk genes and biological pathways. Here, we apply polygenic, linkage disequilibrium-informed methods to the largest available case-control, genome-wide association study summary statistic data for melanoma and PD. We identify positive and significant genetic correlation (correlation: 0.17, 95% CI 0.10-0.24; P = 4.09 × 10-06) between melanoma and PD. We further demonstrate melanoma and PD-inferred gene expression to overlap across tissues (correlation: 0.14, 95% CI 0.06 to 0.22; P = 7.87 × 10-04) and highlight seven genes including PIEZO1, TRAPPC2L, and SOX6 as potential mediators of the genetic correlation between melanoma and PD. These findings demonstrate specific, shared genetic architecture between PD and melanoma that manifests at the level of gene expression.

The Medical Genome Reference Bank contains whole genome and phenotype data of 2570 healthy elderly

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing.

Healthspan and healthy aging are areas of research with potential socioeconomic impact. Here, the authors present the Medical Genome Reference Bank (MGRB) which consist of over 4,000 individuals aged 70 years and older without a history of the major age-related diseases and report on results from whole-genome sequencing and association analyses.

Role of the Aryl Hydrocarbon Receptor in Environmentally Induced Skin Aging and Skin Carcinogenesis

The skin is constantly exposed to a variety of environmental threats, including solar electromagnetic radiation, microbes, airborne particulate matter, and chemicals. Acute exposure to these environmental factors results in the activation of different signaling pathways that orchestrate adaptive stress responses to maintain cell and tissue homeostasis. Chronic exposure of skin to these factors, however, may lead to the accumulation of damaged macromolecules and loss of cell and tissue integrity, which, over time, may facilitate aging processes and the development of aging-related malignancies. One transcription factor that is expressed in all cutaneous cells and activated by various environmental stressors, including dioxins, polycyclic aromatic hydrocarbons, and ultraviolet radiation, is the aryl hydrocarbon receptor (AHR). By regulating keratinocyte proliferation and differentiation, epidermal barrier function, melanogenesis, and immunity, a certain degree of AHR activity is critical to maintain skin integrity and to adapt to acute stress situations. In contrast, a chronic activation of cutaneous AHR signaling critically contributes to premature aging and the development of neoplasms by affecting metabolism, extracellular matrix remodeling, inflammation, pigmentation, DNA repair, and apoptosis. This article provides an overview of the detrimental effects associated with sustained AHR activity in chronically stressed skin and pinpoints AHR as a promising target for chemoprevention.

Does polygenic risk influence associations between sun exposure and melanoma? A prospective cohort analysis

BACKGROUND

Melanoma develops as the result of complex interactions between sun exposure and genetic factors. However, data on these interactions from prospective studies are scant.

OBJECTIVES

To quantify the association between ambient and personal ultraviolet exposure and incident melanoma in a large population-based prospective study of men and women residing in a setting of high ambient ultraviolet radiation, and to examine potential gene-environment interactions.

METHODS

Data were obtained from the QSkin Sun and Health Study, a prospective cohort study of men and women aged 40-69 years, randomly sampled from the Queensland population in 2011. Participants were genotyped and assessed for ultraviolet exposure.

RESULTS

Among participants with genetic data (n = 15 373), 420 (2·7%) developed cutaneous melanoma (173 invasive, 247 in situ) during a median follow-up time of 4·4 years. Country of birth, age at migration, having > 50 sunburns in childhood or adolescence, and a history of keratinocyte cancer or actinic lesions were significantly associated with melanoma risk.

CONCLUSIONS

An interaction with polygenic risk was suggested: among people at low polygenic risk, markers of cumulative sun exposure (as measured by actinic damage) were associated with melanoma. In contrast, among people at high polygenic risk, markers of high-level early-life ambient exposure (as measured by place of birth) were associated with melanoma (hazard ratio for born in Australia vs. overseas 3·16, 95% confidence interval 1·39-7·22). These findings suggest interactions between genotype and environment that are consistent with divergent pathways for melanoma development. What's already known about this topic? The relationship between sun exposure and melanoma is complex, and exposure effects are highly modified by host factors and behaviours. The role of genotype on the relationship between ultraviolet radiation exposure and melanoma risk is poorly understood. What does this study add? We found that country of birth, age at migration, sunburns in childhood or adolescence, and history of keratinocyte cancer or actinic lesions were significantly associated with melanoma risk, while other measures of continuous or more intermittent patterns of sun exposure were not. We found evidence for gene-environment interactions that are consistent with divergent pathways for melanoma development. Linked Comment: Cust. Br J Dermatol 2020; 183:205-206. Plain language summary available online.

Germline variants in oculocutaneous albinism genes and predisposition to familial cutaneous melanoma

Approximately 1%-2% of cutaneous melanoma (CM) is classified as strongly familial. We sought to investigate unexplained CM predisposition in families negative for the known susceptibility genes using next-generation sequencing of affected individuals. Segregation of germline variants of interest within families was assessed by Sanger sequencing. Several heterozygous variants in oculocutaneous albinism (OCA) genes: TYR, OCA2, TYRP1 and SLC45A2, were present in our CM cohort. OCA is a group of autosomal recessive genetic disorders, resulting in pigmentation defects of the eyes, hair and skin. Missense variants classified as pathogenic for OCA were present in multiple families and some fully segregated with CM. The functionally compromised TYR p.T373K variant was present in three unrelated families. In OCA2, known pathogenic variants: p.V443I and p.N489D, were present in three families and one family, respectively. We identified a likely pathogenic SLC45A2 frameshift variant that fully segregated with CM in a family of four cases. Another four-case family harboured cosegregating variants (p.A24T and p.R153C) of uncertain functional significance in TYRP1. We conclude that rare, heterozygous variants in OCA genes confer moderate risk for CM.

Is there a causal relationship between vitamin D and melanoma risk? A Mendelian randomization study

BACKGROUND

Several preclinical studies have identified the antiproliferative effects of 25-hydroxyvitamin D [25(OH)D; vitamin D]. Ultraviolet radiation (UVR) is essential for vitamin D synthesis yet increases the risk of melanoma. Observational studies on the association of vitamin D levels with melanoma risk have reported inconclusive results, and are difficult to interpret owing to the potential confounding from the dual role of UVR.

OBJECTIVES

To determine whether there is a causal association between genetically predicted 25(OH)D concentrations and melanoma using a Mendelian randomization (MR) approach.

METHODS

We performed MR using summary data from a large genome-wide association study (GWAS) meta-analysis of melanoma risk, consisting of 12 874 cases and 23 203 controls. Five single nucleotide polymorphisms associated with 25(OH)D concentration - rs12785878, rs10741657, rs2282679, rs6013897 and rs116970203 - were selected as instrumental variables. An inverse variance weighted method was used to access the evidence for causality. MR results from the melanoma meta-analysis were combined with results from an MR study based on a melanoma risk GWAS using UK Biobank data.

RESULTS

A 20 nmol L^-1 decrease in 25(OH)D was not associated with melanoma risk [odds ratio (OR) 1·06, 95% confidence interval (CI) 0·95-1·19]. Results from the UK Biobank were concordant with this, with meta-analysis of our and UK Biobank-derived MR causal estimates showing no association (OR 1·02, 95% CI 0·92-1·13 for a 20 nmol L^-1 decrease).

CONCLUSIONS

The results suggest that vitamin D levels may not be causally associated with the risk of melanoma. What's already known about this topic? Antitumour activity of vitamin D has been identified in preclinical studies. Observational studies link vitamin D deficiency with an increased risk of a range of cancers. There is a growing public interest for vitamin D supplementation. Observational studies of melanoma are fraught with difficulties because while higher ultraviolet radiation levels increase vitamin D levels, such exposure is also associated with increased melanoma risk. Results from observational studies are inconclusive regarding the effect of vitamin D on melanoma risk. What does this study add? Using Mendelian randomization, an approach to causal inference, which is analogous to a natural randomized controlled trial, we found no causal association between vitamin D levels and melanoma.