Genome-wide association study identifies variants at 9p21 and 22q13 associated with development of cutaneous nevi

Melanocytic naevi, melanocytomas and emerging concepts

With the elucidation of the genetics of melanocytic tumours, new concepts have emerged. An important one is the identification of 'intermediate' melanocytic tumours, those with genetic progression events beyond those of melanocytic naevi but that are not fully malignant. Thus, melanocytic tumours exist on a genetic spectrum that likely corresponds to biological behaviour. There are multiple pathways to melanoma development with different initiating events and characteristic benign melanocytic neoplasms and the precise placement of tumours on these pathways remains to be established and the corresponding risks of progression quantified. In this review, I discuss the classification of melanocytic naevi based on clinical, histopathological and genetic features, as well as the concept of melanocytomas with discussion of specific recognisable subtypes.

The Complex Interplay between Nevi and Melanoma: Risk Factors and Precursors

One effort to combat the rising incidence of malignant melanoma is focused on early detection by the clinical and dermoscopic screening of melanocytic nevi. However, the interaction between nevi, which are congenital or acquired benign melanocytic proliferations, and melanoma is still enigmatic. On the one hand, the majority of melanomas are thought to form de novo, as only a third of primary melanomas are associated with a histologically identifiable nevus precursor. On the other hand, an increased number of melanocytic nevi is a strong risk factor for developing melanoma, including melanomas that do not derive from nevi. The formation of nevi is modulated by diverse factors, including pigmentation, genetic risk factors, and environmental sun exposure. While the molecular alterations that occur during the progression of a nevus to melanoma have been well characterized, many unanswered questions remain surrounding the process of nevus to melanoma evolution. In this review, we discuss clinical, histological, molecular, and genetic factors that influence nevus formation and progression to melanoma.

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

BACKGROUND

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

PLA2G6 Silencing Suppresses Melanoma Progression and Affects Ferroptosis Revealed by Quantitative Proteomics

Although phospholipase A2 group VI (PLA2G6) is involved in oncogenesis in several human tumors, its expression and role in cutaneous malignant melanoma (CMM) pathogenesis remains unclear. Here, by using the Oncomine and CCLE online database, immunohistochemistry, RT-qPCR, and western blotting analysis, we revealed that PLA2G6 was markedly up-regulated in CMM tissues compared to nevus tissues, as well as remarkably increased in vitro in SK-MEL-28 and M14 melanoma cell lines compared to human melanocytes. In vivo, PLA2G6 was also elevated in nine melanoma tissues compared to adjacent tissues. To investigate the malignant behaviors of PLA2G6 in CMM, SK-MEL-28 and M14 cell lines with PLA2G6 stable knockdown by RNAi strategy were constructed. Through CCK8 and colony formation assays in vitro and xenograft tumor experiment in vivo, we found that knockdown of PLA2G6 dramatically inhibited cell proliferation. The results of scratch-wound and transwell assays suggested that the migration and invasion of melanoma cells were prominently suppressed after silencing PLA2G6. In addition, flow cytometry showed that the knockdown of PLA2G6 promoted the apoptosis rate of melanoma cells. To further explore the potential molecular mechanism, we used liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) proteomic and bioinformatics analysis. The GO and KEGG analysis suggested that the underlying mechanism of PLA2G6 in CMM might be associated with the ferroptosis pathway, and ferroptosis-related proteins were validated to be differentially expressed in PLA2G6 knockdown SK-MEL-28 and M14 cells. Together, these results suggested that PLA2G6 knockdown significantly inhibited cell proliferation, metastasis, and promoted apoptosis in melanoma. Our findings on the biological function of PLA2G6 and the underlying association between PLA2G6 and ferroptosis in melanoma may contribute to developing a potential therapeutic strategy for melanoma.

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

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

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

The role of MC1R gene variants and phenotypical features in predicting high nevus count

Variants in the Melanocortin 1 Receptor (MC1R) gene have been associated with an increased risk of melanoma, but the role in nevus count is unclear. We investigated if specific MC1R gene variants or the number of MC1R gene variants and phenotypical features were associated with nevus count. A total of 494 participants of the 'Leiden skin cancer study' were included and the MC1R gene coding sequence was analysed by single-strand conformation polymorphism analysis followed by sequencing of unknown variants. The association between MC1R gene variants and nevus count and the association between age, gender and phenotypical features and nevus count were studied using the Chi-square test. Study of nine frequently occurring MC1R gene variants in participants without skin cancer (n = 203) showed that the 'r' Val60Leu variant was significantly associated with high nevus count (>50 nevi) (P = 0.017). This association was very strong among women (P < 0.001), but not present among men. Having one or two MC1R variants in general did not show a significant difference in the nevus count. Hair colour, skin type, eye colour and age were not significantly associated with nevus count, whereas gender showed a significant association (P = 0.008), with the highest nevus counts in female. The Val60Leu variant of the MC1R gene could be a promising candidate as an independent predictor of high nevus count, particularly in women. This information about the genetic makeup could promote personalized follow-up strategies and might help to prevent skin cancer in the future.

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.

Influence of germline genetic variants on dermoscopic features of melanoma

Nevus count is highly determined by inherited variants and has been associated with the origin of melanoma. De novo melanomas (DNMMs) are more prevalent in patients with a low nevus count and have distinctive dermoscopic features than nevus-associated melanomas. We evaluated the impact of nine single nucleotide polymorphisms (SNPs) of MTAP (rs10811629, rs2218220, rs7023329 and rs751173), PLA2G6 (rs132985 and rs2284063), IRF4 (rs12203592), and PAX3 (rs10180903 and rs7600206) genes associated with nevus count and melanoma susceptibility, and the MC1R variants on dermoscopic features of 371 melanomas from 310 patients. All MTAP variants associated with a low nevus count were associated with regression structures (peppering and mixed regression), blue-whitish veil, shiny white structures, and pigment network. SNPs of PLA2G6 (rs132985), PAX3 (rs7600206), and IRF4 (rs12203592) genes were also associated with either shiny white structures or mixed regression (all corrected p-values ≤ .06). Melanomas from red hair color MC1R variants carriers showed lower total dermoscopy score (p-value = .015) and less blotches than melanomas from non-carriers (p-value = .048). Our results provide evidence that germline variants protective for melanoma risk and/or associated with a low nevus count are associated with certain dermoscopic features, more characteristic of de novo and worse prognosis melanomas.

Lipid metabolic Reprogramming: Role in Melanoma Progression and Therapeutic Perspectives

Simple Summary

Melanoma is a devastating skin cancer characterized by an impressive metabolic plasticity. Melanoma cells are able to adapt to the tumor microenvironment by using a variety of fuels that contribute to tumor growth and progression. In this review, the authors summarize the contribution of the lipid metabolic network in melanoma plasticity and aggressiveness, with a particular attention to specific lipid classes such as glycerophospholipids, sphingolipids, sterols and eicosanoids. They also highlight the role of adipose tissue in tumor progression as well as the potential antitumor role of drugs targeting critical steps of lipid metabolic pathways in the context of melanoma.

Abstract

Metabolic reprogramming contributes to the pathogenesis and heterogeneity of melanoma. It is driven both by oncogenic events and the constraints imposed by a nutrient- and oxygen-scarce microenvironment. Among the most prominent metabolic reprogramming features is an increased rate of lipid synthesis. Lipids serve as a source of energy and form the structural foundation of all membranes, but have also emerged as mediators that not only impact classical oncogenic signaling pathways, but also contribute to melanoma progression. Various alterations in fatty acid metabolism have been reported and can contribute to melanoma cell aggressiveness. Elevated expression of the key lipogenic fatty acid synthase is associated with tumor cell invasion and poor prognosis. Fatty acid uptake from the surrounding microenvironment, fatty acid β-oxidation and storage also appear to play an essential role in tumor cell migration. The aim of this review is (i) to focus on the major alterations affecting lipid storage organelles and lipid metabolism. A particular attention has been paid to glycerophospholipids, sphingolipids, sterols and eicosanoids, (ii) to discuss how these metabolic dysregulations contribute to the phenotype plasticity of melanoma cells and/or melanoma aggressiveness, and (iii) to highlight therapeutic approaches targeting lipid metabolism that could be applicable for melanoma treatment.

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.

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.

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.

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.

Parkinson's Disease and Melanoma: Co-Occurrence and Mechanisms

Parkinson's disease (PD) is a neurodegenerative disorder that is characterized by loss of dopaminergic neurons in the substantia nigra pars compacta, depletion of dopamine in the striatum and the presence of Lewy bodies. Cancer is uncontrolled growth of cells in the body and migration of these cells from their site of origin to other parts of the body. PD and cancer are two opposite diseases, one arising from cell proliferation and the other from cell degeneration. This fundamental difference is consistent with inverse comorbidity between most cancers and neurodegenerative diseases. However, a positive association of PD and melanoma has been reported which has recently become of significant interest. A link between PD and cancer has been supported by many epidemiological studies, most of which show that PD patients have a lower risk of developing most cancers than the general population. However, the mechanisms underlying this epidemiological observation are not known. In this review we focus on epidemiological studies correlating PD and melanoma and the possible mechanisms underlying the co-occurrence of the two diseases. We explore possible explanations for the important observations that more PD patients develop melanoma that would otherwise be expected and vice-versa.

Associations of pigmentary and naevus phenotype with melanoma risk in two populations with comparable ancestry but contrasting levels of ambient sun exposure

BACKGROUND

People at high risk of developing melanoma are usually identified by pigmentary and naevus phenotypes.

OBJECTIVE

We examined whether associations of these phenotypes with melanoma risk differed by ambient sun exposure or participant characteristics in two population-based, case-control studies with comparable ancestry but different ambient sun exposure.

METHODS

Data were analysed from 616 cases and 496 controls from the Australian Melanoma Family Study and 2012 cases and 504 controls from the Leeds (UK) case-control study. Questionnaire, interview and dermatological skin examination data were collected using the same measurement protocols. Relative risks were estimated as odds ratios using unconditional logistic regression, adjusted for potential confounders.

RESULTS

Hair and skin colour were the strongest pigmentary phenotype risk factors. All associations of pigmentary phenotype with melanoma risk were similar across countries. The median number of clinically assessed naevi was approximately three times higher in Australia than Leeds, but the relative risks for melanoma associated with each additional common or dysplastic naevus were higher for Leeds than Australia, especially for naevi on the upper and lower limbs. Higher naevus counts on the head and neck were associated with a stronger relative risk for melanoma for women than men. The two countries had similar relative risks for melanoma based on self-reported naevus density categories, but personal perceptions of naevus number differed by country. There was no consistent evidence of interactions between phenotypes on risk.

CONCLUSIONS

Classifying people at high risk of melanoma based on their number of naevi should ideally take into account their country of residence, type of counts (clinical or self-reported), body site on which the naevus counts are measured and sex. The presence of naevi may be a stronger indicator of a genetic predisposition in the UK than in Australia based on less opportunity for sun exposure to influence naevus development.

Body site-specific genetic effects influence naevus count distribution in women

Body site is highly relevant for melanoma: it affects prognosis and varies according to the patient's sex. The distribution of naevi, a major risk factor for melanoma, at different body sites also varies according to sex in childhood. Using naevus counts at different body sites in 492 unrelated adults from both sexes, we observed that women have an increased number of naevi on the lower limbs compared to men (p = 8.5 × 10^-5 ), showing that a high naevus count on this site persists from childhood throughout life. Then, using data from 3,232 twins, we observed, in women, the lowest naevus count heritability on the trunk (26%), and the highest on the lower limbs (69%). Finally, we showed that, in 2,864 women, six genomic loci previously associated with both naevus count and melanoma risk (IRF4, DOCK8, MTAP, 9q31.2, KITLG and PLA2G6) have an effect on naevus count that is body site-specific, but whose effect sizes are predominantly stronger on the lower limbs. Sex-specific genetic influence on naevus count at different sites may explain differences in site-specific melanoma incidence as well as prognosis between sexes.

Phenotypic and genotypic analysis of amelanotic and hypomelanotic melanoma patients

BACKGROUND

Amelanotic/hypomelanotic melanoma is associated with poorer outcomes due to a more advanced disease stage at diagnosis.

OBJECTIVE

To determine phenotypic risks and genotypic associations with amelanotic/hypomelanotic melanoma to develop a clinical and genetic profile that could assist in identifying high-risk individuals.

METHODS

The Brisbane Naevus Morphology Study conducted from 2009 to 2016 has recruited a core of 1254 participants. Participants were drawn from a combination of volunteers from dermatology outpatient clinics, private dermatology clinics, the Brisbane Longitudinal Twin Study and QSkin study. Case participants had a personal history of melanoma and control participants no personal history of melanoma. We specifically examined seven known candidate pigmentation and melanoma genes and pigmentary phenotypic characteristics in participants with amelanotic/hypomelanotic melanoma compared to pigmented melanomas. This assayed single nucleotide polymorphisms in MC1R, TYR, HERC/OCA2, IRF4, MTAP, PLA2G6 and MITF.

RESULTS

Forty-seven participants had at least one amelanotic/hypomelanotic melanoma, and 389 had pigmented melanomas, with amelanotic/hypomelanotic melanoma patients significantly older than pigmented melanoma participants (63.3 ± 13.0 vs. 54.6 ± 15.3 years; P < 0.001). Amelanotic/hypomelanotic melanoma patients were more likely than pigmented melanoma patients to have red hair (34% vs. 15%; P = 0.01), severe hand freckling (13% vs. 5%; P = 0.01) and propensity to sunburn (63% vs. 44%; P = 0.01). MC1R R/R genotype was much more frequent in our amelanotic/hypomelanotic melanoma population (31.1% vs. 11%; P < 0.001; OR 26.4 vs. 5.9; control 1.0). Amelanotic/hypomelanotic melanoma was associated with TYR rs1126809*A/A [OR (CI 95%) 2.7 (1.1-6.8) vs. 1.2 (0.8-1.9)] and PLA2G6 rs11570734*A/A [OR (CI 95%) 3.7 (1.0-13.6) vs. 1.3 (0.9-2.0)]. The MTAP melanoma risk SNP genotype, associated with darker pigmentation, (rs4636294*A/A) was less common in amelanotic/hypomelanotic melanoma patients [OR (CI 95%) 0.8 (0.3-2.1) vs. 2.0 (1.3-3.1)].

CONCLUSIONS

Knowledge of phenotypic and genotypic associations of amelanotic/hypomelanotic melanoma can help predict risks and associations of this difficult to diagnose melanoma, which may ultimately assist clinical management and patient skin self-examination.

Systematic evaluation of cancer-specific genetic risk score for 11 types of cancer in The Cancer Genome Atlas and Electronic Medical Records and Genomics cohorts

BACKGROUND

Genetic risk score (GRS) is an odds ratio (OR)-weighted and population-standardized method for measuring cumulative effect of multiple risk-associated single nucleotide polymorphisms (SNPs). We hypothesize that GRS is a valid tool for risk assessment of most common cancers.

METHODS

Utilizing genotype and phenotype data from The Cancer Genome Atlas (TCGA) and Electronic Medical Records and Genomics (eMERGE), we tested 11 cancer-specific GRSs (bladder, breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, prostate, renal, and thyroid cancer) for association with the respective cancer type. Cancer-specific GRSs were calculated, for the first time in these cohorts, based on previously published risk-associated SNPs using the Caucasian subjects in these two cohorts.

RESULTS

Mean cancer-specific GRS in the population controls of eMERGE approximated the expected value of 1.00 (between 0.98 and 1.02) for all 11 types of cancer. Mean cancer-specific GRS was consistently higher in respective cancer patients than controls for all 11 types of cancer (P < 0.05). When subjects were categorized into low-, average-, and high-risk groups based on cancer-specific GRS (<0.5, 0.5-1.5, and >1.5, respectively), significant dose-response associations of higher cancer-specific GRS with higher OR of respective type of cancer were found for nine types of cancer (P_-trend  < 0.05). More than 64% subjects in the population controls of eMERGE can be classified as high risk for at least one type of these cancers.

CONCLUSION

Validity of GRS for predicting cancer risk is demonstrated for most types of cancer. If confirmed in larger studies, cancer-specific GRS may have the potential for developing personalized cancer screening strategy.

Melanoma susceptibility variant rs869330 in the MTAP gene is associated with melanoma outcome

The rising incidence of cutaneous melanoma (CM), an aggressive skin cancer, emphasizes the need for novel biomarkers to guide personalized care and better predict outcome. Genetic factors including germline risk variants are promising candidates for this aim. We explored the association between germline risk variants and melanoma outcome in a large genetically homogenous Belgian melanoma population, focusing on single nucleotide polymorphisms which generated the highest association with melanoma susceptibility. Between 2004 and 2014, blood samples of 1088 patients with histologically confirmed CM were collected and genotyped for nine variants. Cox proportional hazard models were used to assess the association between each single nucleotide polymorphism and relapse-free survival and overall survival, adjusted by age, sex, melanoma stage, site, and subtype. We identified significant associations for rs869330 (in the methylthioadenosine phosphorylase - MTAP gene) with overall survival (hazard ratio = 0.760, P = 0.048, 95% confidence interval: 0.580-0.998) and relapse-free survival (hazard ratio = 0.800, P = 0.020, 95% confidence interval: 0.650-0.970). This exploratory study is the first to show a significant association between the rs869330 variant (in the MTAP gene) and outcome in a large CM population.

Combining common genetic variants and non-genetic risk factors to predict risk of cutaneous melanoma

Melanoma heritability is among the highest for cancer and single nucleotide polymorphisms (SNPs) contribute to it. To date, only SNPs that reached statistical significance in genome-wide association studies or few candidate SNPs have been included in melanoma risk prediction models. We compared four approaches for building polygenic risk scores (PRS) using 12 874 melanoma cases and 23 203 controls from Melanoma Meta-Analysis Consortium as a training set, and newly genotyped 3102 cases and 2301 controls from the MelaNostrum consortium for validation. We estimated adjusted odds ratios (ORs) for melanoma risk using traditional melanoma risk factors and the PRS with the largest area under the receiver operator characteristics curve (AUC). We estimated absolute risks combining the PRS and other risk factors, with age- and sex-specific melanoma incidence and competing mortality rates from Italy as an example. The best PRS, including 204 SNPs (AUC = 64.4%; 95% confidence interval (CI) = 63-65.8%), developed using winner's curse estimate corrections, had a per-quintile OR = 1.35 (95% CI = 1.30-1.41), corresponding to a 3.33-fold increase comparing the 5th to the 1st PRS quintile. The AUC improvement by adding the PRS was up to 7%, depending on adjusted factors and country. The 20-year absolute risk estimates based on the PRS, nevus count and pigmentation characteristics for a 60-year-old Italian man ranged from 0.5 to 11.8% (relative risk  = 26.34), indicating good separation.

Novel pleiotropic risk loci for melanoma and nevus density implicate multiple biological pathways

The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.

Inherited Genetic Variants Associated with Melanoma BRAF/NRAS Subtypes

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

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

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

Construction of lncRNA-miRNA-mRNA networks reveals functional lncRNAs in abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is one of the most significant causes of morbidity and mortality in populations aged >65 years worldwide. However, the underlying mechanisms of AAA based on the competitive endogenous RNA (ceRNA) hypothesis have remained elusive. In the present study, differently expressed long non-coding RNA (lncRNA)-microRNA (miRNA)-mRNA networks in AAA were constructed by analyzing public datasets, including GSE7084, GSE24194 from rats and that of a previous study. A total of 1,219 mRNAs, 2,093 lncRNAs and 57 miRNAs were identified to differently express in AAA. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to explore the potential roles of differently expressed lncRNAs based on their regulating mRNAs. Based on the ceRNA hypothesis, lncRNA-miRNA-mRNA networks in AAA were, for the first time, constructed at a system-wide level. The present study identified 5 upregulated lncRNAs [nuclear paraspeckle assembly transcript 1, cyclin-dependent kinase inhibitor 2B antisense RNA 1, small Cajal body-specific RNA 10, AC005224.4 and SUMO1/sentrin/SMT3-specific peptidase 3-eukaryotic translation initiation factor 4A1] and the downregulated zinc ribbon domain containing 1 antisense RNA 1 as key lncRNAs in ceRNA networks. To the best of our knowledge, the present study was the first to screen ceRNA networks in AAA. In addition, key lncRNA-mRNA-biological processes analysis indicated that these key lncRNAs were involved in regulating signal transduction, protein amino acid phosphorylation, immune response, transcription, development and cell differentiation. The present study provides novel clues to explore the molecular mechanisms of AAA progression in terms of lncRNA implication.

New susceptibility loci for cutaneous melanoma risk and progression revealed using a porcine model

Despite major advances, it is estimated that a large part of melanoma predisposing genes remains to be discovered. Animal models of spontaneous diseases are valuable tools and experimental crosses can be used to identify and fine-map new susceptibility loci associated with melanoma. We performed a Genome-Wide Association Study (GWAS) of melanoma occurrence and progression (clinical ulceration and presence of metastasis) in a porcine model of spontaneous melanoma, the MeLiM pig. Five loci on chromosomes 2, 5, 7, 8 and 16 showed genome-wide significant associations (p < 5 × 10^–6) with either one of these phenotypes. Suggestive associations (p < 5 × 10^–5) were also found at 16 additional loci. Moreover, comparison of the porcine results to those reported by human melanoma GWAS indicated shared association signals notably at CDKAL1 and TERT loci but also nearby CCND1, FTO, PLA2G6 and TMEM38B-RAD23B loci. Extensive search of the literature revealed a potential key role of genes at the identified porcine loci in tumor invasion (DST, PLEKHA5, CBY1, LIMK2 and ETV5) and immune response modulation (ETV5, HERC3 and DICER1) of the progression phenotypes. These biological processes are consistent with the clinico-pathological features of MeLiM tumors and can open new routes for future melanoma research in humans.

Vitamin D in melanoma: Controversies and potential role in combination with immune check-point inhibitors

The role of vitamin D in melanoma is still controversial. Although several Authors described a correlation between vitamin D deficiency and poor survival in metastatic melanoma patients, clinical trials exploring the effects of vitamin D supplementation in this clinical setting were mostly inconclusive. However, recent evidence suggests that vitamin D exerts both anti-proliferative effects on tumor cells and immune-modulating activities, that have been widely explored in auto-immune disorders. On the one hand, vitamin D has been shown to inhibit T-helper17 lymphocytes, notoriously involved in the pathogenesis of immune-related adverse events (iAEs) which complicate immune-checkpoint inhibitor (ICI) treatment. On the other hand, vitamin D up-regulates PDL-1 expression on both epithelial and immune cells, suggesting a synergic effect in combination with ICIs, for which further investigation is needed.

Telomerase Reverse Transcriptase Polymorphism rs2736100: A Balancing Act between Cancer and Non-Cancer Disease, a Meta-Analysis

The enzyme telomerase reverse transcriptase (TERT) is essential for telomere maintenance. In replicating cells, maintenance of telomere length is important for the preservation of vital genetic information and prevention of genomic instability. A common genetic variant in TERT, rs2736100 C/A, is associated with both telomere length and multiple diseases. Carriage of the C allele is associated with longer telomere length, while carriage of the A allele is associated with shorter telomere length. Furthermore, some diseases have a positive association with the C and some with the A allele. In this study, meta-analyses were performed for two groups of diseases, cancerous diseases, e.g., lung cancer and non-cancerous diseases, e.g., pulmonary fibrosis, using data from genome-wide association studies and case-control studies. In the meta-analysis it was found that cancer positively associated with the C allele (pooled OR 1.16 [95% CI 1.09–1.23]) and non-cancerous diseases negatively associated with the C allele (pooled OR 0.81 [95% CI 0.65–0.99]). This observation illustrates that the ambiguous role of telomere maintenance in disease hinges, at least in part, on a single locus in telomerase genes. The dual role of this single nucleotide polymorphism also emphasizes that therapeutic agents aimed at influencing telomere maintenance should be used with caution.

Truncating mutations of TP53AIP1 gene predispose to cutaneous melanoma

Genetic predisposition to cutaneous malignant melanoma (CMM) involves highly penetrant predisposing genes and low and intermediate penetrant predisposing alleles. However, the missing heritability in (CMM) is still high. For such and in order to identify new genetic factors for CMM, we conducted an exome sequencing study in high-risk CMM patients. Two rounds of exome sequencing were successively performed in 33 and 27 high-risk patients. We focused on genes carrying rare nonsense, frameshift, and splice variants (allelic frequency <1%) that were present in both series of exomes. An extension study was then conducted in a large cohort (1 079 CMM patients and 1 230 Caucasian ethnically matched healthy controls), and the inactivating variants frequency was compared between groups using two-sided Fisher exact test. Two TP53AIP1 truncating mutations were identified in four patients: a frameshift c.63_64insG, p.Q22Afs*81 in two patients from the same family and in the proband of a second family; and a nonsense mutation c.95 C > A, p.Ser32Stop in a patient with multiple CMMs. In all patients, TP53AIP1 truncating variants were strongly associated with CMM risk (two-sided Fisher exact test = 0.004, OR = 3.3[1.3-8.5]). Additionally, we showed that TP53AIP1 mRNA was strongly down-regulated throughout different phases of melanoma progression. TP53AIP1 gene is a TP53 target which plays a key role by inducting apoptosis in response to UV-induced DNA damage. Constitutional mutations of TP53AIP1 had previously been involved in susceptibility to prostate cancer. Our results show that constitutional truncating TP53AIP1 mutations predispose to CMM in the French population. Replication studies in other populations should be performed.

Rare germline variants in known melanoma susceptibility genes in familial melanoma

Known high-risk cutaneous malignant melanoma (CMM) genes account for melanoma risk in <40% of melanoma-prone families, suggesting the existence of additional high-risk genes or perhaps a polygenic mechanism involving multiple genetic modifiers. The goal of this study was to systematically characterize rare germline variants in 42 established melanoma genes among 144 CMM patients in 76 American CMM families without known mutations using data from whole-exome sequencing. We identified 68 rare (<0.1% in public and in-house control datasets) nonsynonymous variants in 25 genes. We technically validated all loss-of-function, inframe insertion/deletion, and missense variants predicted as deleterious, and followed them up in 1, 559 population-based CMM cases and 1, 633 controls. Several of these variants showed disease co-segregation within families. Of particular interest, a stopgain variant in TYR was present in five of six CMM cases/obligate gene carriers in one family and a single population-based CMM case. A start gain variant in the 5'UTR region of PLA2G6 and a missense variant in ATM were each seen in all three affected people in a single family, respectively. Results from rare variant burden tests showed that familial and population-based CMM patients tended to have higher frequencies of rare germline variants in albinism genes such as TYR, TYRP1, and OCA2 (P < 0.05). Our results suggest that rare nonsynonymous variants in low- or intermediate-risk CMM genes may influence familial CMM predisposition, warranting further investigation of both common and rare variants in genes affecting functionally important pathways (such as melanogenesis) in melanoma risk assessment.

Germline Variation at CDKN2A and Associations with Nevus Phenotypes among Members of Melanoma Families

Germline mutations in CDKN2A are frequently identified among melanoma kindreds and are associated with increased atypical nevus counts. However, a clear relationship between pathogenic CDKN2A mutation carriage and other nevus phenotypes including counts of common acquired nevi has not yet been established. Using data from GenoMEL, we investigated the relationships between CDKN2A mutation carriage and 2-mm, 5-mm, and atypical nevus counts among blood-related members of melanoma families. Compared with individuals without a pathogenic mutation, those who carried one had an overall higher prevalence of atypical (odds ratio = 1.64; 95% confidence interval = 1.18-2.28) nevi but not 2-mm nevi (odds ratio = 1.06; 95% confidence interval = 0.92-1.21) or 5-mm nevi (odds ratio = 1.26; 95% confidence interval = 0.94-1.70). Stratification by case status showed more pronounced positive associations among non-case family members, who were nearly three times (odds ratio = 2.91; 95% confidence interval = 1.75-4.82) as likely to exhibit nevus counts at or above the median in all three nevus categories simultaneously when harboring a pathogenic mutation (vs. not harboring one). Our results support the hypothesis that unidentified nevogenic genes are co-inherited with CDKN2A and may influence carcinogenesis.

Melanocytic nevi and melanoma: unraveling a complex relationship

Approximately 33% of melanomas are derived directly from benign, melanocytic nevi. Despite this, the vast majority of melanocytic nevi, which typically form as a result of BRAF^V600E-activating mutations, will never progress to melanoma. Herein, we synthesize basic scientific insights and data from mouse models with common observations from clinical practice to comprehensively review melanocytic nevus biology. In particular, we focus on the mechanisms by which growth arrest is established after BRAF^V600E mutation. Means by which growth arrest can be overcome and how melanocytic nevi relate to melanoma are also considered. Finally, we present a new conceptual paradigm for understanding the growth arrest of melanocytic nevi in vivo termed stable clonal expansion. This review builds upon the canonical hypothesis of oncogene-induced senescence in growth arrest and tumor suppression in melanocytic nevi and melanoma.

Sun exposure and skin cancer, and the puzzle of cutaneous melanoma: A perspective on Fears et al. Mathematical models of age and ultraviolet effects on the incidence of skin cancer among whites in the United States. American Journal of Epidemiology 1977; 105: 420-427

Sunlight has been known as an important cause of skin cancer since around the turn of the 20th Century. A 1977 landmark paper of US scientists Fears, Scotto, and Schneiderman advanced a novel hypothesis whereby cutaneous melanoma was primarily caused by intermittent sun exposure (i.e. periodic, brief episodes of exposure to high-intensity ultraviolet radiation) while the keratinocyte cancers, squamous cell carcinoma and basal cell carcinoma, were primarily caused by progressive accumulation of sun exposure. With respect to cutaneous melanoma, this became known as the intermittent exposure hypothesis. The hypothesis stemmed from analysis of measured ambient ultraviolet radiation and age-specific incidence rates of melanoma and keratinocyte cancers collected as an extension to the US Third National Cancer Survey in several US States. In this perspective paper, we put this novel hypothesis into the context of knowledge at the time, and describe subsequent epidemiological and molecular research into melanoma that elaborated the intermittent exposure hypothesis and ultimately replaced it with a dual pathway hypothesis. Our present understanding is of two distinct biological pathways by which cutaneous melanoma might develop; a nevus prone pathway initiated by early sun exposure and promoted by intermittent sun exposure or possibly host factors; and a chronic sun exposure pathway in sun sensitive people who progressively accumulate sun exposure to the sites of future melanomas.

The Molecular Revolution in Cutaneous Biology: The Era of Genome-Wide Association Studies and Statistical, Big Data, and Computational Topics

The investigation of biological systems involving all organs of the body including the skin is in era of big data. This requires heavy-duty computational tools, and novel statistical methods. Microarrays have allowed the interrogation of thousands of common genetic markers in thousands of individuals from the same population (termed genome wide association studies or GWAS) to reveal common variation associated with disease or phenotype. These markers are usually single nucleotide polymorphisms (SNPs) that are relatively common in the population. In the case of dermatological diseases such as alopecia areata, vitiligo, psoriasis and atopic dermatitis, common variants have been identified that are associated with disease, and these provide insights into biological pathways and reveal possible novel drug targets. Other skin phenotypes such as acne, color and skin cancers are also being investigated with GWAS. Analyses of such large GWAS datasets require a consideration of a number of statistical issues including the testing of multiple markers, population substructure, and ultimately a requirement for replication. There are also issues regarding the missing heritability of disease that cannot be entirely explained with current GWAS approaches. Next generation sequencing technologies such as exome and genome sequencing of similar patient cohorts will reveal additional variants contributing to disease susceptibility. However, the data generated with these approaches will be orders of magnitude greater than that those generated with arrays, with concomitant challenges in the identification of disease causing variants.

Environmental factors, seven GWAS-identified susceptibility loci, and risk of gastric cancer and its precursors in a Chinese population

Gene–environment interactions may increase gastric cancer (GC) risk. Seven susceptibility loci identified by genome‐wide association studies (GWASs) suggest that genetic factors play a role in gastric carcinogenesis. Meanwhile, Helicobacter pylori (H. pylori) infection, smoking, and alcohol drinking are also important environmental factors for gastric cancer. However, studies to explore the role of gene–environment interactions in gastric carcinogenesis, and particularly the relationship between the seven susceptibility loci and their potential interactions with H. pylori infection, smoking, and alcohol drinking in risk of GC, and severe intestinal metaplasia (IM)/dysplasia, have been inconclusive. A total of 1273 subjects in a Chinese population were recruited, and genotyping was carried out using the competitive allele‐specific PCR (KASP) method. Unconditional logistic regression was applied to model the associations between genetic polymorphisms and the disease risk. Effect modifications by H. pylori infection, smoking and alcohol drinking were evaluated. PSCA rs2294008/rs2976392 showed a significant, multiplicative interaction with H. pylori infection in risk of GC. Meanwhile, PRKAA1 rs13361707 had an additive interaction with H. pylori infection. SLC52A3 rs13042395 showed an interaction with alcohol drinking in risk of GC. Moreover, three SNPs, MUC1 rs4072037, ZBTB20 rs9841504 and PRKAA1 rs13361707, were associated with precancerous gastric lesions (severe IM/dysplasia). Our data suggest that genetic predisposition factors identified by GWAS may interact with environmental risk factors, Particularly for H. pylori infection and alcohol consumption, to increase the risk of GC.

Lipoquality control by phospholipase A2 enzymes

The phospholipase A₂ (PLA₂) family comprises a group of lipolytic enzymes that typically hydrolyze the sn-2 position of glycerophospholipids to give rise to fatty acids and lysophospholipids. The mammalian genome encodes more than 50 PLA₂s or related enzymes, which are classified into several subfamilies on the basis of their structures and functions. From a general viewpoint, the PLA₂ family has mainly been implicated in signal transduction, producing bioactive lipid mediators derived from fatty acids and lysophospholipids. Recent evidence indicates that PLA₂s also contribute to phospholipid remodeling for membrane homeostasis or energy production for fatty acid β-oxidation. Accordingly, PLA₂ enzymes can be regarded as one of the key regulators of the quality of lipids, which I herein refer to as lipoquality. Disturbance of PLA₂-regulated lipoquality hampers tissue and cellular homeostasis and can be linked to various diseases. Here I overview the current state of understanding of the classification, enzymatic properties, and physiological functions of the PLA₂ family.

Higher Nevus Count Exhibits a Distinct DNA Methylation Signature in Healthy Human Skin: Implications for Melanoma

High nevus count is the strongest risk factor for melanoma, and although gene variants have been discovered for both traits, epigenetic variation is unexplored. We investigated 322 healthy human skin DNA methylomes associated with total body nevi count, incorporating genetic and transcriptomic variation. DNA methylation changes were identified at genes involved in melanocyte biology, such as RAF1 (P = 1.2 × 10^-6) and CTC1 (region: P = 6.3 × 10^-4), and other genes including ARRDC1 (P = 3.1 × 10^-7). A subset exhibited coordinated methylation and transcription changes within the same biopsy. The total analysis was also enriched for melanoma-associated DNA methylation variation (P = 6.33 × 10^-6). In addition, we show that skin DNA methylation is associated in cis with known genome-wide association study single nucleotide polymorphisms for nevus count, at PLA2G6 (P = 1.7 × 10^-49) and NID1 (P = 6.4 × 10^-14), as well as melanoma risk, including in or near MC1R, MX2, and TERT/CLPTM1L (P < 1 × 10^-10). Our analysis using a uniquely large dataset comprising healthy skin DNA methylomes identified known and additional regulatory loci and pathways in nevi and melanoma biology. This integrative study improves our understanding of predisposition to nevi and their potential contribution to melanoma pathogenesis.

Association of common variants in MTAP with susceptibility and overall survival of osteosarcoma: a two-stage population-based study in Han Chinese

Osteosarcoma (OS) is a common malignant tumor, which exists widely in the bone of children and adolescents, and genetic factors may influence its susceptibility. Recently, the gene MTAP has been reported to be associated with OS in a Caucasian population. To investigate the association of common variants in MTAP with OS risk in Han Chinese individuals, we designed a two-stage case-control study with 392 OS patients and 1,578 unrelated healthy controls of Han Chinese individuals. A total of 17 tagging single nucleotide polymorphisms (SNPs) were firstly genotyped in the discovery stage, and single-SNP association and haplotypic association analyses have been performed. The SNP rs7023329 was found to be strongly associated with the OS risk (adjusted P = 0.002908), and the results of odds ratios (ORs) and 95% confidence intervals (CI) revealed increased risks from A allele of the SNP on OS (OR=1.33, 95% CI=1.13-1.62). The results were confirmed with a similar pattern in the validation stage (adjusted P = 0.006737, OR=1.49, 95% CI=1.11-2.00). Moreover, haplotypic analyses indicated that one haplotype block containing rs7023329 was significantly associated with OS risk in both stages (both global P<0.0001). The statistically significant association between the rs7023329 genotype and poor survival in OS patients was also observed. To sum up, our results prove that MTAP plays an important role in the etiology of OS, suggesting this gene as a potential genetic modifier for OS development.

Nevus count associations with pigmentary phenotype, histopathological melanoma characteristics and survival from melanoma

Although nevus count is an established risk factor for melanoma, relationships between nevus number and patient and tumor characteristics have not been well studied and the influence of nevus count on melanoma-specific survival is equivocal. Using data from the Genes, Environment and Melanoma (GEM) study, a large population-based study of primary cutaneous melanoma, we evaluated associations between number of nevi and patient features, including sun-sensitivity summarized in a phenotypic index, and tumor characteristics. We also assessed the association of nevus count with melanoma-specific survival. Higher nevus counts were independently and positively associated with male gender and younger age at diagnosis, and they were inversely associated with lentigo maligna histology. We observed a borderline significant trend of poorer melanoma-specific survival with increasing quartile of nevus count, but little or no association between number of nevi and pigmentary phenotypic characteristics or prognostic tumor features.

Long noncoding RNA variations in cardiometabolic diseases

Cardiometabolic diseases are characterized as a combination of multiple risk factors for cardiovascular disease (CVD) and metabolic diseases including diabetes mellitus, dyslipidemia, hypertension and abdominal obesity. This cluster of abnormalities individually and interdependently leads to atherosclerosis and CVD morbidity and mortality. In the past decade, genome-wide association studies (GWASs) have identified a series of cardiometabolic disease-associated variants that can collectively explain a small proportion of the variability. Intriguingly, the susceptibility variants imputed from GWASs usually do not reside in the coding regions, suggesting a crucial role of the noncoding elements of the genome. In recent years, emerging evidence suggests that noncoding RNA (ncRNA) is functional for physiology and pathophysiology of human diseases. These include microRNAs and long noncoding RNAs (lncRNAs) that are now implicated in human diseases. The ncRNAs can interact with each other and with proteins, to interfere gene expressions, leading to the development of many human disorders. Although evidence suggests the functional role of lncRNAs in cardiometabolic traits, the molecular mechanisms of gene regulation underlying cardiometabolic diseases remain to be better defined. Here, we summarize the recent discoveries of lncRNA variations in the context of cardiometabolic diseases.

The genomic landscape of cutaneous melanoma

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

The role of the melanoma gene MC1R in Parkinson disease and REM sleep behavior disorder

The MC1R gene, suggested to be involved in Parkinson disease (PD) and melanoma, was sequenced in PD patients (n = 539) and controls (n = 265) from New York, and PD patients (n = 551), rapid eye movement sleep behavior disorder (RBD) patients (n = 351), and controls (n = 956) of European ancestry. Sixty-eight MC1R variants were identified, including 7 common variants with frequency > 0.01. None of the common variants was associated with PD or RBD in the different regression models. In a meta-analysis with fixed-effect model, the p.R160W variant was associated with an increased risk for PD (odds ratio = 1.22, 95% confidence interval = 1.02-1.47, p = 0.03) but with significant heterogeneity (p = 0.048). Removing one study that introduced the heterogeneity resulted in nonsignificant association (odds ratio = 1.11, 95% confidence interval, 0.92-1.35, p = 0.27, heterogeneity p = 0.57). Rare variants had similar frequencies in patients and controls (10.54% and 10.15%, respectively, p = 0.75), and no cumulative effect of carrying more than one MC1R variant was found. The present study does not support a role for the MC1R p.R160W and other variants in susceptibility for PD or RBD.

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, NM_002386:c. 456C > A (p.TYR152*), or 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.

New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk

To increase our understanding of the genetic basis of adiposity and its links to cardiometabolic disease risk, we conducted a genome-wide association meta-analysis of body fat percentage (BF%) in up to 100,716 individuals. Twelve loci reached genome-wide significance (P<5 × 10(-8)), of which eight were previously associated with increased overall adiposity (BMI, BF%) and four (in or near COBLL1/GRB14, IGF2BP1, PLA2G6, CRTC1) were novel associations with BF%. Seven loci showed a larger effect on BF% than on BMI, suggestive of a primary association with adiposity, while five loci showed larger effects on BMI than on BF%, suggesting association with both fat and lean mass. In particular, the loci more strongly associated with BF% showed distinct cross-phenotype association signatures with a range of cardiometabolic traits revealing new insights in the link between adiposity and disease risk.

The associations between Parkinson's disease and cancer: the plot thickens

Epidemiological studies support a general inverse association between the risk of cancer development and Parkinson’s disease (PD). In recent years however, increasing amount of eclectic evidence points to a positive association between PD and cancers through different temporal analyses and ethnic groups. This positive association has been supported by several common genetic mutations in SNCA, PARK2, PARK8, ATM, p53, PTEN, and MC1R resulting in cellular changes such as mitochondrial dysfunction, aberrant protein aggregation, and cell cycle dysregulation. Here, we review the epidemiological and biological advances of the past decade in the association between PD and cancers to offer insight on the recent and sometimes contradictory findings.