Molecular Characterization of Melanoma Cases in Denmark Suspected of Genetic Predisposition

Both environmental and host factors influence risk of cutaneous melanoma (CM), and worldwide, the incidence varies depending on constitutional determinants of skin type and pigmentation, latitude, and patterns of sun exposure. We performed genetic analysis of CDKN2A, CDK4, BAP1, MC1R, and MITFp.E318K in Danish high-risk melanoma cases and found CDKN2A germline mutations in 11.3% of CM families with three or more affected individuals, including four previously undescribed mutations. Rare mutations were also seen in CDK4 and BAP1, while MC1R variants were common, occurring at more than twice the frequency compared to Danish controls. The MITF p.E318K variant similarly occurred at an approximately three-fold higher frequency in melanoma cases than controls. To conclude, we propose that mutation screening of CDKN2A and CDK4 in Denmark should predominantly be performed in families with at least 3 cases of CM. In addition, we recommend that testing of BAP1 should not be conducted routinely in CM families but should be reserved for families with CM and uveal melanoma, or mesothelioma.     

Inherited Variation at MC1R and Histological Characteristics of Primary Melanoma

Variation in the melanocortin-1receptor (MC1R) gene is associated with pigmentary phenotypes and risk of malignant melanoma. Few studies have reported on MC1R variation with respect to tumor characteristics, especially clinically important prognostic features. We examined associations between MC1R variants and histopathological melanoma characteristics. Study participants were enrolled from nine geographic regions in Australia, Canada, Italy and the United States and were genotyped for MC1R variants classified as high-risk [R] (D84E, R142H, R151C, R160W, and D294H, all nonsense and insertion/deletion) or low-risk [r] (all other nonsynonymous) variants. Tissue was available for 2,160 white participants of the Genes, Environment and Melanoma (GEM) Study with a first incident primary melanoma diagnosis, and underwent centralized pathologic review. No statistically significant associations were observed between MC1R variants and AJCC established prognostic tumor characteristics: Breslow thickness, presence of mitoses or presence of ulceration. However, MC1R was significantly associated with anatomic site of melanoma (p = 0.002) and a positive association was observed between carriage of more than one [R] variant and melanomas arising on the arms (OR = 2.39; 95% CI: 1.40, 4.09). We also observed statistically significant differences between sun-sensitive and sun-resistant individuals with respect to associations between MC1R genotype and AJCC prognostic tumor characteristics. Our results suggest inherited variation in MC1R may play an influential role in anatomic site presentation of melanomas and may differ with respect to skin pigmentation phenotype.     

MC1R variation and melanoma risk in the Swedish population in relation to clinical and pathological parameters

The genetic background of cutaneous malignant melanoma (CMM) includes both germ line aberrations in high‐penetrance genes, like CDKN2A, and allelic variation in low‐penetrance genes like the melanocortin‐1 receptor gene, MC1R. Red‐hair colour associated MC1R alleles (RHC) have been associated with red hair, fair skin and risk of CMM. We investigated MC1R and CDKN2A variation in relation to phenotype, clinical factors and CMM risk in the Swedish population. The study cohort consisted of sporadic primary melanoma patients, familial melanoma patients and a control group. An allele‐dose dependent increase in melanoma risk for carriers of variant MC1R alleles (after adjusting for phenotype), with an elevated risk among familial CMM patients, was observed. This elevated risk was found to be significantly associated with an increased frequency of dysplastic nevi (DN) among familial patients compared to sporadic patients. MC1R variation was found to be less frequent among acral lentiginous melanomas (ALM) and dependent on tumour localisation. No association was found between CDKN2A gene variants and general melanoma risk. Two new variants in the POMC gene were identified in red haired individuals without RHC alleles.     

Replication and Predictive Value of SNPs Associated with Melanoma and Pigmentation Traits in a Southern European Case-Control Study

Background

Genetic association studies have revealed numerous polymorphisms conferring susceptibility to melanoma. We aimed to replicate previously discovered melanoma-associated single-nucleotide polymorphisms (SNPs) in a Greek case-control population, and examine their predictive value.

Methods

Based on a field synopsis of genetic variants of melanoma (MelGene), we genotyped 284 patients and 284 controls at 34 melanoma-associated SNPs of which 19 derived from GWAS. We tested each one of the 33 SNPs passing quality control for association with melanoma both with and without accounting for the presence of well-established phenotypic risk factors. We compared the risk allele frequencies between the Greek population and the HapMap CEU sample. Finally, we evaluated the predictive ability of the replicated SNPs.

Results

Risk allele frequencies were significantly lower compared to the HapMap CEU for eight SNPs (rs16891982 – SLC45A2, rs12203592 – IRF4, rs258322 – CDK10, rs1805007 – MC1R, rs1805008 - MC1R, rs910873 - PIGU, rs17305573- PIGU, and rs1885120 - MTAP) and higher for one SNP (rs6001027 – PLA2G6) indicating a different profile of genetic susceptibility in the studied population. Previously identified effect estimates modestly correlated with those found in our population (r = 0.72, P<0.0001). The strongest associations were observed for rs401681-T in CLPTM1L (odds ratio [OR] 1.60, 95% CI 1.22–2.10; P = 0.001), rs16891982-C in SCL45A2 (OR 0.51, 95% CI 0.34–0.76; P = 0.001), and rs1805007-T in MC1R (OR 4.38, 95% CI 2.03–9.43; P = 2×10⁻⁵). Nominally statistically significant associations were seen also for another 5 variants (rs258322-T in CDK10, rs1805005-T in MC1R, rs1885120-C in MYH7B, rs2218220-T in MTAP and rs4911442-G in the ASIP region). The addition of all SNPs with nominal significance to a clinical non-genetic model did not substantially improve melanoma risk prediction (AUC for clinical model 83.3% versus 83.9%, p = 0.66).

Conclusion

Overall, our study has validated genetic variants that are likely to contribute to melanoma susceptibility in the Greek population.     

Pigmentary Markers in Danes – Associations with Quantitative Skin Colour,Nevi Count,Familial Atypical Multiple-Mole,and Melanoma Syndrome

To investigate whether pigmentation genes involved in the melanogenic pathway (melanogenesis) contributed to melanoma predisposition, we compared pigmentary genetics with quantitative skin pigmentation measurements, the number of atypical nevi, the total nevus count, and the familial atypical multiple mole and melanoma (FAMMM) syndrome. We typed 32 pigmentary SNP markers and sequenced MC1R in 246 healthy individuals and 116 individuals attending periodic control for malignant melanoma development, 50 of which were diagnosed with FAMMM. It was observed that individuals with any two grouped MC1R variants (missense, {"type":"entrez-nucleotide","attrs":{"text":"NM_002386","term_id":"193083133"}}NM_002386:c. 456C > A (p.TYR152*), or {"type":"entrez-nucleotide","attrs":{"text":"NM_002386","term_id":"193083133"}}NM_002386:c.83_84insA (p.Asn29Glnfs*14) had significantly (p<0.001) lighter skin pigmentation of the upper-inner arm than those with none or one MC1R variant. We did not observe any significant association of the MC1R variants with constitutive pigmentation measured on the buttock area. We hypothesize that the effect of MC1R variants on arm pigmentation is primarily reflecting the inability to tan when subjected to UVR. A gender specific effect on skin pigmentation was also observed, and it was found that the skin pigmentation of females on average were darker than that of males (p<0.01). We conclude that MC1R variants are associated with quantitative skin colour in a lightly pigmented Danish population. We did not observe any association between any pigmentary marker and the FAMMM syndrome. We suggest that the genetics of FAMMM is not related to the genetics of the pigmentary pathway.     

<Emphasis Type="Italic">CDKN2A</Emphasis> and <Emphasis Type="Italic">MC1R</Emphasis> variants found in Cypriot patients diagnosed with cutaneous melanoma

The prevalence of genetic variants associated to cutaneous melanoma (CM) has never been determined within Cypriot melanomas. This study evaluates the frequency of variants in cyclin-dependent kinase inhibitor 2A (CDKN2A) and melanocortin-1 receptor (MC1R) in 32 patients diagnosed with CM. Other characteristics and risk factors were also assessed. CDKN2A p.Ala148Thr was detected in three of 32 patients, while the control group revealed no variations within CDKN2A. MC1R screening in 32 patients revealed the following variations: p.Val60Leu in 11 patients, p.Arg142His in four patients, p.Thr314Thr in one patient, p.Arg160Trp in one patient, p.Val92Met/p.Thr314Thr in one patient and p.Val92Met/p.Arg142His/p.Thr314Thr in one patient. The control group revealed only p.Val60Leu (in 10 of 45 individuals), which is frequently found in general populations. Two unrelated patients carried CDKN2A p.Ala148Thr in combination with MC1R p.Arg142His, suggesting digenic inheritance that may provide evidence of different gene variants acting synergistically to contribute for CM development. This study confirms the presence of CDKN2A and MC1R variants among Cypriot melanomas and supports existing evidence of a role for these variants in susceptibility to melanoma.     

Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.     

Vitamin D receptor variants and the malignant melanoma risk: A population-based study

Background: There is continuing interest in identifying low-penetrance genes which are associated with an increased susceptibility to common types of cancer, including malignant melanoma. Methods: We sought to examine the association between four VDR common variants (rs1544410, rs731236, rs10735810, rs4516035) and the risk of melanoma in the Polish population. We also determined the prevalence of compound carriers of VDR and known MM genetic risk factors MC1R and CDKN2A (A148T) variants. We examined 763 unselected melanoma cases, 763 healthy adults matched for sex and age with the melanoma cases and 777 newborns. Results: None of the VDR variants alone or as compound carriers of two or more of the VDR genotypes were associated with MM risk. There were no major differences between the prevalences of the examined variants among patients with MM on UV-exposed and UV-non exposed skin areas, as well as among early-onset and late-onset cases. We found no association between VDR and MC1R or between VDR and CDKN2A common variants. A statistically significant over-representation of one VDR haplotype: rs731236_A + rs1544410_T (OR = 3.2, p = 0.02) was detected. Linkage disequilibrium of rs1544410 and rs731236 was confirmed. Conclusion: To answer the question, whether VDR can be regarded as melanoma susceptibility gene, additional, large multi-center association studies have to be performed.     

Natural and orthogonal model for estimating gene–gene interactions applied to cutaneous melanoma

Epistasis, or gene–gene interaction, results from joint effects of genes on a trait; thus, the same alleles of one gene may display different genetic effects in different genetic backgrounds. In this study, we generalized the coding technique of a natural and orthogonal interaction (NOIA) model for association studies along with gene–gene interactions for dichotomous traits and human complex diseases. The NOIA model which has non-correlated estimators for genetic effects is important for estimating influence from multiple loci. We conducted simulations and data analyses to evaluate the performance of the NOIA model. Both simulation and real data analyses revealed that the NOIA statistical model had higher power for detecting main genetic effects and usually had higher power for some interaction effects than the usual model. Although associated genes have been identified for predisposing people to melanoma risk: HERC2 at 15q13.1, MC1R at 16q24.3 and CDKN2A at 9p21.3, no gene–gene interaction study has been fully explored for melanoma. By applying the NOIA statistical model to a genome-wide melanoma dataset, we confirmed the previously identified significantly associated genes and found potential regions at chromosomes 5 and 4 that may interact with the HERC2 and MC1R genes, respectively. Our study not only generalized the orthogonal NOIA model but also provided useful insights for understanding the influence of interactions on melanoma risk.     

Family of melanocortin receptor (MCR) genes in mammals—mutations, polymorphisms and phenotypic effects

The melanocortin receptor gene family consists of five single-exon members, which are located on autosomes. Three genes (MC2R, MC4R and MC5R) are syntenic in the human, mouse, cattle and dog genomes, while in the pig, the syntenic group comprises MC1R, MC2R and MC5R. Two genes (MC1R and MC4R) have been extensively studied due to their function in melanogenesis (MC1R) and energy control (MC4R). Conservative organisation of these genes in five mammalian species (human, mouse, cattle, pig and dog), in terms of the encoded amino acid sequence, is higher in the case of MC4R compared to MC1R. Polymorphisms of these two genes are responsible or associated with variation of pigmentation (MC1R) and adipose tissue deposition (MC4R). Polymorphic variants in MC1R, causing coat colour variation, were described in humans and domestic mammals (cattle, horse, pig, sheep, dog), as well as farm red and arctic foxes. The MC4R gene is very polymorphic in humans and it is well known that some variants cause monogenic obesity or significantly contribute to the development of polygenic obesity. Such relationships are not so evident in domestic mammals; however, at least one missense substitution (298Asp?>?Asn) in the porcine MC4R significantly contributes, at least in some breeds, to fat tissue accumulation, feed conversion ratio and daily weight gain. Knowledge on the phenotypic effects of polymorphisms of MC2R, MC3R and MC5R in domestic mammals is scarce, probably due to the small number of reports addressing these genes. Thus, further studies focused on these genes should be undertaken.     

Genetic Association Study of Adiposity and Melanocortin-4 Receptor (MC4R) Common Variants: Replication and Functional Characterization of Non-Coding Regions

Common genetic variants 3′ of MC4R within two large linkage disequilibrium (LD) blocks spanning 288 kb have been associated with common and rare forms of obesity. This large association region has not been refined and the relevant DNA segments within the association region have not been identified. In this study, we investigated whether common variants in the MC4R gene region were associated with adiposity-related traits in a biracial population-based study. Single nucleotide polymorphisms (SNPs) in the MC4R region were genotyped with a custom array and a genome-wide array and associations between SNPs and five adiposity-related traits were determined using race-stratified linear regression. Previously reported associations between lower BMI and the minor alleles of rs2229616/Val103Ile and rs52820871/Ile251Leu were replicated in white female participants. Among white participants, rs11152221 in a proximal 3′ LD block (closer to MC4R) was significantly associated with multiple adiposity traits, but SNPs in a distal 3′ LD block (farther from MC4R) were not. In a case-control study of severe obesity, rs11152221 was significantly associated. The association results directed our follow-up studies to the proximal LD block downstream of MC4R. By considering nucleotide conservation, the significance of association, and proximity to the MC4R gene, we identified a candidate MC4R regulatory region. This candidate region was sequenced in 20 individuals from a study of severe obesity in an attempt to identify additional variants, and the candidate region was tested for enhancer activity using in vivo enhancer assays in zebrafish and mice. Novel variants were not identified by sequencing and the candidate region did not drive reporter gene expression in zebrafish or mice. The identification of a putative insulator in this region could help to explain the challenges faced in this study and others to link SNPs associated with adiposity to altered MC4R expression.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

In vitro behavior and UV response of melanocytes derived from carriers of CDKN2A mutations and MC1R variants

Coinheritance of germline mutation in cyclin‐dependent kinase inhibitor 2A (CDKN2A) and loss‐of‐function (LOF) melanocortin 1 receptor (MC1R) variants is clinically associated with exaggerated risk for melanoma. To understand the combined impact of these mutations, we established and tested primary human melanocyte cultures from different CDKN2A mutation carriers, expressing either wild‐type MC1R or MC1RLOF variant(s). These cultures expressed the CDKN2A product p16 (INK4A) and functional MC1R. Except for 32ins24 mutant melanocytes, the remaining cultures showed no detectable aberrations in proliferation or capacity for replicative senescence. Additionally, the latter cultures responded normally to ultraviolet radiation (UV) by cell cycle arrest, JNK, p38, and p53 activation, hydrogen peroxide generation, and repair of DNA photoproducts. We propose that malignant transformation of melanocytes expressing CDKN2A mutation and MC1RLOF allele(s) requires acquisition of somatic mutations facilitated by MC1R genotype or aberrant microenvironment due to CDKN2A mutation in keratinocytes and fibroblasts.     

MC1R Is a Prognostic Marker and Its Expression Is Correlated with MSI in Colorectal Cancer

Melanocortin 1 receptor (MC1R) is thought to be a marker of poor prognosis and a potential target for the treatment of melanoma. Studies have found that MC1R promotes several tumor behaviors, including cell proliferation and differentiation, pigment formation, and genome damage repair. Some single-nucleotide polymorphisms (SNPs) of MC1R are involved in the occurrence and development of melanoma. A few studies have reported a relationship between MC1R and colorectal cancer (CRC). In this research, our objective was to examine MC1R expression and MC1R SNPs and investigate their correlation with the clinicopathological features of human CRC tissues. We evaluated MC1R mRNA expression by performing bioinformatic analyses on human CRC expression datasets. We used Western blotting and RT-qPCR to compare MC1R expression in CRC tissues with that in normal tissues, and MC1R SNPs in CRC tissues were detected by PCR-direct sequencing (DS). The expression of MC1R was significantly decreased in CRC tissues compared with normal tissue, and its expression was negatively associated with P53 expression, MLH1 expression, and PMS2 expression, and high MC1R expression was significantly associated with microsatellite instability (MSI). MC1R SNPs were also associated with the clinicopathological characteristics of CRC; for example, the rs2228479 locus genotype was correlated with Ki67 status, and the rs885479 locus genotype was correlated with age and T stage. In conclusion, MC1R plays a crucial role in the progression of CRC and may be a marker of poor prognosis in CRC.     

Association between Common Polymorphism near the MC4R Gene and Obesity Risk: A Systematic Review and Meta-Analysis

Background

Genome-wide association studies on Europeans have shown that two polymorphisms (rs17782313, rs12970134) near the melanocortin 4 receptor (MC4R) gene were associated with increased risk of obesity. Subsequently studies among different ethnic populations have shown mixed results with some confirming and others showing inconsistent results, especially among East Asians and Africans. We performed a comprehensive meta-analysis of various studies from different ethnic populations to assess the association of the MC4R polymorphism with obesity risk.

Methods

We retrieved all published literature that investigated association of MC4R variants with obesity from PubMed and Embase. Pooled odds ratio (OR) with 95% confidence interval (CI) was calculated using fixed- or random-effects model.

Results

A total of 61 studies (80,957 cases/220,223 controls) for rs17782313 polymorphism (or proxy) were included in the meta-analysis. The results suggested that rs17782313 polymorphism was significantly associated with obesity risk (OR = 1.18, 95%CI = 1.15–1.21, p<0.001). Similar trends were observed among subgroups of Europeans and East Asians, adults and children, studies with high quality score, and for each five MC4R polymorphisms independently.

Conclusions

The present meta-analysis confirms the significant association of MC4R polymorphism with risk of obesity. Further studies should be conducted to identify the causal variant and the underlying mechanisms of the identified association.