POT1 and multiple primary melanomas: the dermatological phenotype.

POT1 is the second most frequently reported gene (after CDKN2A) in familial melanoma. Pathogenic variants are associated with earlier onset and/or multiple primary melanomas (MPMs). To date, POT1 phenotypical reports have been largely restricted to associated malignancies, and description of the dermatological landscape has been limited. We identified 10 variants in n=18 of 384 (4.7%) unrelated individuals (n=13 MPMs; n=5 single primary melanomas) of European ancestry. Five variants were rare (minor allele frequency <0.001) or novel (two loss-of-function (LOF), one splice acceptor and two missense) and were predicted to be functionally significant, in five unrelated probands with MPMs (≥3 melanomas). We performed three-dimensional total body photography on both individuals with confirmed pathogenic LOF variants to characterise the dermatological phenotype. Total body naevus counts (≥2 mm diameter) were significantly higher (p=7.72×10-12) in carriers compared with a control population. Majority of naevi were on the probands' back and lower limb regions, where only mild to moderate ultraviolet (UV) damage was observed. Conversely, the head/neck region, where both probands exhibited severe UV damage, had comparably fewer naevi. We hypothesise that carriage of functionally significant POT1 variants is associated with increased naevus counts generally, and naevi >5 mm in diameter specifically and the location of these are independent of UV damage.

Assessing the genetic risk of nodular melanoma using a candidate gene approach.

BACKGROUND: Nodular melanoma (NM) is a challenge to diagnose early due to its rapid growth and more atypical clinical presentation, making it the largest contributor to melanoma mortality. OBJECTIVES: Our study aim was to perform a rare-variant allele (RVA) analysis of whole-exome sequencing of patients with NM and non-NM (minor allele frequency ≤ 1% non-Finnish European) for a set of 500 candidate genes potentially implicated in melanoma. METHODS: This study recruited 131 participants with NM and 194 with non-NM from South-east Queensland and patients with NM from Victoria to perform a comparative analysis of possible genetic differences or similarities between the two melanoma cohorts. RESULTS: Phenotypic analysis revealed that a majority of patients diagnosed with NM were older males with a higher frequency of fair skin and red hair than is seen in the general population. The distribution of common melanoma polygenic risk scores was similar in patients with NM and non-NM, with over 28% in the highest quantile of scores. There was also a similar frequency of carriage of familial/high-penetrant melanoma gene and loss-of-function variants. We identified 39 genes by filtering 500 candidate genes based on the greatest frequency in NM compared with non-NM cases. The genes with RVAs of greatest frequency in NM included PTCH1, ARID2 and GHR. Rare variants in the SMO gene, which interacts with PTCH1 as ligand and receptor, were also identified, providing evidence that the Hedgehog pathway may contribute to NM risk. There was a cumulative effect in carrying multiple rare variants in the NM-associated genes. A 14.8-fold increased ratio for NM compared with non-NM was seen when two RVAs of the 39 genes were carried by a patient. CONCLUSIONS: This study highlights the importance of considering frequency of RVA to identify those at risk of NM in addition to known high penetrance genes.

The MC1R r allele does not increase melanoma risk in MITF E318K carriers.

BACKGROUND: Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored. OBJECTIVES: To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals. MATERIALS AND METHODS: Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank). RESULTS: The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals. CONCLUSIONS: RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

Amelanotic/hypopigmented melanoma in a sibship with oculocutaneous albinism.

Oculocutaneous albinism (OCA) is a rare condition characterized by hypopigmentation. A female proband and her sister, both with primary amelanotic/hypopigmented melanoma, underwent three-dimensional total-body photography and dermoscopy. Both sisters had exome sequencing along with their brother, who had OCA but no history of melanoma. Imaging analysis was consistent with OCA in terms of individual typology angle scores, degree of sun damage, and high naevus counts. Exome data filtered for variants in known OCA and melanoma/naevi susceptibility genes (n = 98) found all siblings were compound heterozygous for TYR mutations (Arg402Ter and Val275Phe), previously reported as causative OCA variants. A rare missense variant in PARP1 (p.Pro377Ser) was solely present in the melanoma-unaffected brother, which is noteworthy as this was previously reported as potentially protective in a familial melanoma pedigree positive for CDKN2A mutations. Evaluation and confirmation of functional impact in larger cohorts could personalize melanoma screening in OCA.

Genetic analysis of multiple primary melanomas arising within the boundaries of congenital nevi depigmentosa.

Here, we present a rare case of a patient who developed multiple primary melanomas within the boundaries of two nevi depigmentosa. The melanomas were excised, and as a preventive measure, the remainder of the nevi depigmentosa were removed. We performed whole-exome sequencing on excised tissue from the nevus depigmentosus, adjacent normal skin, and saliva to explain this intriguing phenomenon. We also performed a GeneTrails Comprehensive Solid Tumor Panel analysis on one of the melanoma tissues. Genetic analysis revealed germline MC1R V92M and TYR R402Q polymorphisms and a MET E168D germline mutation that may have increased the risk of melanoma development. This genetic predisposition, combined with a patient-reported history of substantial sun exposure and sunburns, which were more severe within the boundaries of the nevi depigmentosa due to the lack of photoprotective melanin, produced numerous somatic mutations in the melanocytes of the nevi depigmentosa. Fitting with this paradigm for melanoma development in chronically sun-damaged skin, the patient's melanomas harbored somatic mutations in CDKN2A (splice site), NF1, and ATRX and had a tumor mutation burden in the 90-95th percentile for melanoma.

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.

Genetic variation in IRF4 expression modulates growth characteristics, tyrosinase expression and interferon-gamma response in melanocytic cells.

A SNP within intron4 of the interferon regulatory factor4 (IRF4) gene, rs12203592*C/T, has been independently associated with pigmentation and age-specific effects on naevus count in European-derived populations. We have characterized the cis-regulatory activity of this intronic region and using human foreskin-derived melanoblast strains, we have explored the correlation between IRF4 rs12203592 homozygous C/C and T/T genotypes with TYR enzyme activity, supporting its association with pigmentation traits. Further, higher IRF4 protein levels directed by the rs12203592*C allele were associated with increased basal proliferation but decreased cell viability following UVR, an etiological factor in melanoma development. Since UVR, and accompanying IFNγ-mediated inflammatory response, is associated with melanomagenesis, we evaluated its effects in the context of IRF4 status. Manipulation of IRF4 levels followed by IFNγ treatment revealed a subset of chemokines and immuno-evasive molecules that are sensitive to IRF4 expression level and genotype including CTLA4 and PD-L1.

Molecular analysis of common polymorphisms within the human Tyrosinase locus and genetic association with pigmentation traits.

We have compared the melanogenic activities of cultured melanocytes carrying two common TYR alleles as homozygous 192S-402R wild-type, heterozygous and homozygous variant. This includes assays of TYR protein, DOPAoxidase activity, glycosylation and temperature sensitivity of protein and DOPAoxidase levels. Homozygous wild-type strains on average had higher levels of TYR protein and enzyme activity than other genotypes. Homozygous 402Q/Q melanocytes produced significantly less TYR protein, displayed altered trafficking and glycosylation, with reduced DOPAoxidase. However, near wild-type TYR activity levels could be recovered at lower growth temperature. In a sample population from Southeast Queensland, these two polymorphisms were present on four TYR haplotypes, designated as WT 192S-402R, 192Y-402R and 192S-402Q with a double-variant 192Y-402Q of low frequency at 1.9%. Based on cell culture findings and haplotype associations, we have used an additive model to assess the penetrance of the ten possible TYR genotypes derived from the combination of these haplotypes.

Melanocortin-1 receptor-mediated signalling pathways activated by NDP-MSH and HBD3 ligands.

Binding of melanocortin peptide agonists to the melanocortin-1 receptor of melanocytes results in eumelanin production, whereas binding of the agouti signalling protein inverse agonist results in pheomelanin synthesis. Recently, a novel melanocortin-1 receptor ligand was reported. A ß-defensin gene mutation was found to be responsible for black coat colour in domestic dogs. Notably, the human equivalent, ß-defensin 3, was found to bind with high affinity to the melanocortin-1 receptor; however, the action of ß-defensin as an agonist or antagonist was unknown. Here, we use in vitro assays to show that ß-defensin 3 is able to act as a weak partial agonist for cAMP signalling in human embryonic kidney (HEK) cells expressing human melanocortin-1 receptor. ß-defensin 3 is also able to activate MAPK signalling in HEK cells stably expressing either wild type or variant melanocortin-1 receptors. We suggest that ß-defensin 3 may be a novel melanocortin-1 receptor agonist involved in regulating melanocyte responses in humans.

Analysis of cultured human melanocytes based on polymorphisms within the SLC45A2/MATP, SLC24A5/NCKX5, and OCA2/P loci.

Single nucleotide polymorphisms (SNPs) within the SLC45A2/MATP, SLC24A5/NCKX5, and OCA2/P genes have been associated with natural variation of pigmentation traits in human populations. Here, we describe the characterization of human primary melanocytic cells genotyped for polymorphisms within the MATP, NCKX5, or OCA2 loci. On the basis of genotype, these cultured cells reflect the phenotypes observed by others in terms of both melanin content and tyrosinase (TYR) activity when comparing skin designated as either "White" or "Black". We found a statistically significant association of MATP-374L (darker skin) with higher TYR protein abundance that was not observed for any NCKX5-111 or OCA2 rs12913832 allele. MATP-374L/L homozygous strains displayed significantly lower MATP transcript levels compared to MATP-374F/F homozygous cells, but this did not reach statistical significance based on NCKX5 or OCA2 genotype. Similarly, we observed significantly increased levels of OCA2 mRNA in rs12913832-T (brown eye) homozygotes compared to rs12913832-C (blue eye) homozygous strains, which was not observed for MATP or NCKX5 gene transcripts. In genotype-phenotype associations performed on a collection of 226 southern European individuals using these same SNPs, we were able to show strong correlations in MATP-L374F, OCA2, and melanocortin-1 receptor with skin, eye, and hair color variation, respectively.

PPARgamma agonists attenuate proliferation and modulate Wnt/beta-catenin signalling in melanoma cells.

Peroxisome proliferator-activated receptor-gamma (PPARgamma) is a member of the nuclear hormone receptor (NHR) superfamily of ligand-activated transcriptional regulators. Accumulating evidence suggests that PPARgamma agonists such as the thiazolidinediones (TZDs) may prove to be useful anti-cancer agents exhibiting anti-proliferative and/or pro-apoptotic affects in a range of cancer cell types including melanoma, however, the mechanisms underlying this effect remain unclear. We have demonstrated the anti-proliferative effects of full and partial PPARgamma modulators in human melanoma cell lines. Ablation of PPARgamma expression in the MM96L melanoma cell line by siRNA mediated mechanisms attenuates the anti-proliferative effect of these agents suggesting this effect is directly mediated by PPARgamma. The mechanisms underlying the anti-proliferative effects of PPARgamma in melanoma cells involve the regulation of expression of a number of critical cell cycle genes and beta-catenin. Moreover, our data indicate that PPARgamma modulates Wnt/beta-catenin mediated signalling in melanoma cells in an agonist dependent manner.

Osteonectin downregulates E-cadherin, induces osteopontin and focal adhesion kinase activity stimulating an invasive melanoma phenotype.

Osteonectin is recognised as a marker of metastasis progression in melanoma and has been implicated in the transition from radial to vertical growth phase. A Tetracycline-inducible system was used to regulate Osteonectin protein levels in melanoma cell lines to examine the morphological, biochemical and invasive changes that accompany its altered expression. Assay of protein and phosphorylation changes showed a downregulation of E-cadherin, upregulation of Osteopontin and a corresponding increase in phosphorylation of Focal Adhesion Kinase on Tyr(397) and Tyr(576) concomitant with Osteonectin induction. Melanoma cells overexpressing Osteonectin displayed increased invasive potential, whereas ablation of Osteonectin gene transcription using siRNA suppressed the invasive potential of these cells and resulted in the upregulation of E-cadherin. The recently described interaction of Osteonectin with Integrin Linked Kinase leading to modulation of its activity suggests a mechanism relevant to the loss of E-cadherin and cell adhesion that occurs during melanoma progression. These results indicate a central role for Osteonectin in the regulation of gene expression changes driving the progression of melanoma toward metastasis.

Altered cell surface expression of human MC1R variant receptor alleles associated with red hair and skin cancer risk.

The human melanocortin-1 receptor gene (MC1R) encodes a G-protein coupled receptor that is primarily expressed on melanocytes, where it plays a key role in pigmentation regulation. Variant alleles are associated with red hair colour and fair skin, known as the RHC phenotype, as well as skin cancer risk. The R151C, R160W and D294H alleles, designated 'R', are strongly associated with the RHC phenotype and have been proposed to result in loss of function receptors due to impaired G-protein coupling. We recently provided evidence that the R151C and R160W variants can efficiently couple to G-proteins in response to alpha-melanocyte stimulating hormone. The possibility that altered cellular localization of the R151C and R160W variant receptors could underlie their association with RHC was therefore considered. Using immunofluorescence and ligand binding studies, we found that melanocytic cells exogenously or endogenously expressing MC1R show strong surface localization of the wild-type and D294H alleles but markedly reduced cell surface expression of the R151C and R160W receptors. In additional exogenous expression studies, the R variant D84E and the rare I155T variant, also demonstrated a significant reduction in plasma membrane receptor numbers. The V60L, V92M and R163Q weakly associated RHC alleles, designated 'r', were expressed with normal or intermediate cell surface receptor levels. These results indicate that reduced receptor coupling activity may not be the only contributing factor to the genetic association between the MC1R variants and the RHC phenotype, with MC1R polymorphisms now linked to a change in receptor localization.

Co-expression of SOX9 and SOX10 during melanocytic differentiation in vitro.

Investigations into pigment cell biology have relied on the ability to culture both murine and human melanocytes, numerous melanoma cell lines and more recently, murine and human melanoblasts. Melanoblast culture requires medium supplemented with a range of growth factors including Stem Cell Factor, Endothelin-3 and Fibroblast Growth Factor-2, withdrawal of which causes the cells to differentiate into melanocytes. Using the human melanoblast culture system, we have now examined the expression and/or DNA binding activity of several transcription factors implicated in melanocytic development and differentiation. Of these, the POU domain factor BRN2 and the SOX family member SOX10 are both highly expressed in unpigmented melanocyte precursors but are down-regulated upon differentiation. In contrast, the expression levels of the previously described MITF and PAX3 transcription factors remain relatively constant during the melanoblast-melanocyte transition. Moreover, BRN2 ablated melanoma cells lack expression of SOX10 and MITF but retain PAX3. A novel finding implicates a second SOX protein, SOX9, as a potential melanogenic transcriptional regulator, as its expression level is increased following the down-regulation of BRN2 and SOX10 in differentiated melanoblasts. Our results suggest that a complex network of transcription factor interactions requiring proper temporal coordination is necessary for acquisition and maintenance of the melanocytic phenotype.

Screening of human primary melanocytes of defined melanocortin-1 receptor genotype: pigmentation marker, ultrastructural and UV-survival studies.

Recent population studies have demonstrated an association with the red-hair and fair-skin phenotype with variant alleles of the melanocortin-1 receptor (MC1R) which result in amino acid substitutions within the coding region leading to an altered receptor activity. In particular, Arg151Cys, Arg160Trp and Asp294His were the most commonly associated variants seen in the south-east Queensland population with at least one of these alleles found in 93% of those with red hair. In order to study the individual effects of these variants on melanocyte biology and melanocytic pigmentation, we established a series of human melanocyte strains genotyped for the MC1R receptor which included wild-type consensus, variant heterozygotes, compound heterozygotes and homozygotes for Arg151Cys, Arg160Trp, Val60Leu and Val92Met alleles. These strains ranged from darkly pigmented to amelanotic, with all strains of consensus sequence having dark pigmentation. UV sensitivity was found not to be associated with either MC1R genotype or the level of pigmentation with a range of sensitivities seen across all genotypes. Ultrastructural analysis demonstrated that while consensus strains contained stage IV melanosomes in their terminal dendrites, Arg151Cys and Arg160Trp homozygote strains contained only stage II melanosomes. This was despite being able to show expression of tyrosinase and tyrosinase-related protein-1 markers, although at reduced levels and an ability to convert exogenous 3,4-dihydroxyphenyl-alanine (DOPA) to melanin in these strains.

The role of melanocortin-1 receptor polymorphism in skin cancer risk phenotypes.

We have examined melanocortin-1 receptor (MC1R) variant allele frequencies in the general population and in a collection of adolescent dizygotic and monozygotic twins to determine statistical associations of pigmentation phenotypes with increased skin cancer risk. This included hair and skin color, freckling, mole count and sun exposed skin reflectance. Nine variants were studied and designated as either strong R (OR = 63; 95% CI 32-140) or weak r (OR = 5; 95% CI 3-11) red hair alleles. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. To assess the interaction of the brown eye color gene BEY2/OCA2 on the phenotypic effects of variant MC1R alleles we imputed OCA2 genotype in the twin collection. A modifying effect of OCA2 on MC1R variant alleles was seen on constitutive skin color, freckling and mole count. In order to study the individual effects of these variants on pigmentation phenotype we have established a series of human primary melanocyte strains genotyped for the MC1R receptor. These include strains which are MC1R wild-type consensus, variant heterozygotes, and homozygotes for strong R alleles Arg151Cys and Arg160Trp. Ultrastructural analysis demonstrated that only consensus strains contained stage III and IV melanosomes in their terminal dendrites whereas Arg151Cys and Arg160Trp homozygous strains contained only immature stage I and II melanosomes. Such genetic association studies combined with the functional analysis of MC1R variant alleles in melanocytic cells should provide a link in understanding the association between pigmentary phototypes and skin cancer risk.

Osteonectin/SPARC induction by ectopic beta(3) integrin in human radial growth phase primary melanoma cells.

Expression of the beta(3) integrin subunit in melanoma in situ has been found to correlate with tumor thickness, the ability to invade and metastasize, and poor prognosis. Transition from the radial growth phase (RGP) to the vertical growth phase (VGP) is a critical step in melanoma progression and survival and is distinguished by the expression of beta(3) integrin. The molecular pathways that operate in melanoma cells associated with invasion and metastasis were examined by ectopic induction of the beta(3) integrin subunit in RGP SBcl2 and WM1552C melanoma cells, which converts these cells to a VGP phenotype. We used cDNA representational difference analysis subtractive hybridization between beta(3)-positive and -negative melanoma cells to assess gene expression profile changes accompanying RGP to VGP transition. Fourteen fragments from known genes including osteonectin (also known as SPARC and BM-40) were identified after three rounds of representational difference analysis. Induction of osteonectin was confirmed by Northern and Western blot analysis and immunohistochemistry and correlated in organotypic cultures with the beta(3)-induced progression from RGP to VGP melanoma. Expression of osteonectin was also associated with reduced adhesion to vitronectin, but not to fibronectin. Osteonectin expression was not blocked when melanoma cells were cultured with anti-alpha(v)beta(3) LM609 mAb, mitogen-activated protein kinase, or protein kinase C inhibitors, indicating that other signaling pathway(s) operate through alpha(v)beta(3) integrin during conversion from RGP to VGP.