Genome-wide association study identifies three new melanoma susceptibility loci

Germline mutations predisposing to melanoma and associated malignancies and syndromes: a narrative review

The pathogenesis of melanoma is influenced by a complex combination of environmental factors and individual genetic susceptibility. Familial melanoma refers to cases where there are two first-degree relatives with a melanoma diagnosis. Less strict definitions include second-degree relatives or even three or more of any degree from the same family, although this is not clearly defined in the literature. The term hereditary melanoma is reserved for sporadic or familial melanomas linked to high-risk genes with high penetrance. The first genes related to melanoma were CDKN2A and CDK4, but recently, other genes, mostly tumor suppressor genes, have been described. Internal malignancies, particularly pancreatic cancer, have also been associated with melanoma. Recent studies suggest that there could be a link between melanoma and other neoplasms and tumor predisposition syndromes. This review presents an updated overview of familial melanoma criteria and genes involved in melanoma pathogenesis, emphasizing their clinicopathological aspects and other associated malignancies.

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

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

Malignant Melanoma: An Overview, New Perspectives, and Vitamin D Signaling

Melanoma, originating through malignant transformation of melanin-producing melanocytes, is a formidable malignancy, characterized by local invasiveness, recurrence, early metastasis, resistance to therapy, and a high mortality rate. This review discusses etiologic and risk factors for melanoma, diagnostic and prognostic tools, including recent advances in molecular biology, omics, and bioinformatics, and provides an overview of its therapy. Since the incidence of melanoma is rising and mortality remains unacceptably high, we discuss its inherent properties, including melanogenesis, that make this disease resilient to treatment and propose to use AI to solve the above complex and multidimensional problems. We provide an overview on vitamin D and its anticancerogenic properties, and report recent advances in this field that can provide solutions for the prevention and/or therapy of melanoma. Experimental papers and clinicopathological studies on the role of vitamin D status and signaling pathways initiated by its active metabolites in melanoma prognosis and therapy are reviewed. We conclude that vitamin D signaling, defined by specific nuclear receptors and selective activation by specific vitamin D hydroxyderivatives, can provide a benefit for new or existing therapeutic approaches. We propose to target vitamin D signaling with the use of computational biology and AI tools to provide a solution to the melanoma problem.

Genetic markers of cardiac autonomic neuropathy in the Kazakh population

BACKGROUND

Cardiac autonomic neuropathy (CAN) is a complication of diabetes mellitus (DM) that increases the risk of morbidity and mortality by disrupting cardiac innervation. Recent evidence suggests that CAN may manifest even before the onset of DM, with prediabetes and metabolic syndrome potentially serving as precursors. This study aims to identify genetic markers associated with CAN development in the Kazakh population by investigating the SNPs of specific genes.

MATERIALS AND METHODS

A case-control study involved 82 patients with CAN (cases) and 100 patients without CAN (controls). A total of 182 individuals of Kazakh nationality were enrolled from a hospital affiliated with the RSE "Medical Center Hospital of the President's Affairs Administration of the Republic of Kazakhstan". 7 SNPs of genes FTO, PPARG, SNCA, XRCC1, FLACC1/CASP8 were studied. Statistical analysis was performed using Chi-square methods, calculation of odds ratios (OR) with 95% confidence intervals (CI), and logistic regression in SPSS 26.0.

RESULTS

Among the SNCA gene polymorphisms, rs2737029 was significantly associated with CAN, almost doubling the risk of CAN (OR 2.03(1.09-3.77), p = 0.03). However, no statistically significant association with CAN was detected with the rs2736990 of the SNCA gene (OR 1.00 CI (0.63-1.59), p = 0.99). rs12149832 of the FTO gene increased the risk of CAN threefold (OR 3.22(1.04-9.95), p = 0.04), while rs1801282 of the PPARG gene and rs13016963 of the FLACC1 gene increased the risk twofold (OR 2.56(1.19-5.49), p = 0.02) and (OR 2.34(1.00-5.46), p = 0.05) respectively. rs1108775 and rs1799782 of the XRCC1 gene were associated with reduced chances of developing CAN both before and after adjustment (OR 0.24, CI (0.09-0.68), p = 0.007, and OR 0.43, CI (0.22-0.84), p = 0.02, respectively).

CONCLUSION

The study suggests that rs2737029 (SNCA gene), rs12149832 (FTO gene), rs1801282 (PPARG gene), and rs13016963 (FLACC1 gene) may be predisposing factors for CAN development. Additionally, SNPs rs1108775 and rs1799782 (XRCC1 gene) may confer resistance to CAN. Only one polymorphism rs2736990 of the SNCA gene was not associated with CAN.

The ATM Ser49Cys Variant Effects ATM Function as a Regulator of Oncogene-Induced Senescence

An apical component of the cell cycle checkpoint and DNA damage repair response is the ataxia-telangiectasia mutated (ATM) Ser/Thr protein kinase. A variant of ATM, Ser49Cys (rs1800054; minor allele frequency = 0.011), has been associated with an elevated risk of melanoma development; however, the functional consequence of this variant is not defined. ATM-dependent signalling in response to DNA damage has been assessed in a panel of patient-derived lymphoblastoid lines and primary human melanocytic cell strains heterozygous for the ATM Ser49Cys variant allele. The ATM Ser49Cys allele appears functional for acute p53-dependent signalling in response to DNA damage. Expression of the variant allele did reduce the efficacy of oncogene expression in inducing senescence. These findings demonstrate that the ATM 146C>G Ser49Cys allele has little discernible effect on the acute response to DNA damage but has reduced function observed in the chronic response to oncogene over-expression. Analysis of melanoma, naevus and skin colour genomics and GWAS analyses have demonstrated no association of this variant with any of these outcomes. The modest loss of function detected suggest that the variant may act as a modifier of other variants of ATM/p53-dependent signalling.

ATM Variant as a Cause of Hereditary Cutaneous Melanoma in a Spanish Family: Case Report

Introduction

Ataxia-Telangiectasia Mutated (ATM) is a cancer predisposition gene; carriers of germline pathogenic variants have an increased risk of developing malignancies, including breast, prostate, pancreatic, and ovarian cancer. Most ATM variants are of uncertain significance. Findings from genome-wide association studies (GWAS) suggest that ATM may be a low-risk melanoma susceptibility locus.

Case Report

We report the case of a Hispanic family whose members who have presented cutaneous melanoma have been found to be carriers for the ATM pathogenic variant c.3747-1G>C (rs730881364), one of whom was diagnosed at 24 years old.

Discussion

We describe for the first time the possible clinical association between ATM (c.3747-1G>C) and familial melanoma. In silico splice site analysis predicts that this alteration will weaken the native splice acceptor site and will result in the creation or strengthening of a novel splice acceptor site, assuming a variant that entails loss of functionality that is probably pathogenic and related to oncogenesis. However, we cannot exclude that cutaneous melanoma in both members and at an early age is the result of chance, environmental interaction, other uncontrolled external factors, or the interaction of other genetic alterations other than the ATM variant described in this study.

Susceptibility Genes Associated with Multiple Primary Cancers

With advancements in treatment and screening techniques, we have been witnessing an era where more cancer survivors harbor multiple primary cancers (MPCs), affecting approximately one in six patients. Identifying MPCs is crucial for tumor staging and subsequent treatment choices. However, the current clinicopathological criteria for clinical application are limited and insufficient, making it challenging to differentiate them from recurrences or metastases. The emergence of next-generation sequencing (NGS) technology has provided a genetic perspective for defining multiple primary cancers. Researchers have found that, when considering multiple tumor pairs, it is crucial not only to examine well-known essential mutations like MLH1/MSH2, EGFR, PTEN, BRCA1/2, CHEK2, and TP53 mutations but also to explore certain pleiotropic loci. Moreover, specific deleterious mutations may serve as regulatory factors in second cancer development following treatment. This review aims to discuss these susceptibility genes and provide an explanation of their functions based on the signaling pathway background. Additionally, the association network between genetic signatures and different tumor pairs will be summarized.

Pan-cancer and cross-population genome-wide association studies dissect shared genetic backgrounds underlying carcinogenesis

Integrating genomic data of multiple cancers allows de novo cancer grouping and elucidating the shared genetic basis across cancers. Here, we conduct the pan-cancer and cross-population genome-wide association study (GWAS) meta-analysis and replication studies on 13 cancers including 250,015 East Asians (Biobank Japan) and 377,441 Europeans (UK Biobank). We identify ten cancer risk variants including five pleiotropic associations (e.g., rs2076295 at DSP on 6p24 associated with lung cancer and rs2525548 at TRIM4 on 7q22 nominally associated with six cancers). Quantifying shared heritability among the cancers detects positive genetic correlations between breast and prostate cancer across populations. Common genetic components increase the statistical power, and the large-scale meta-analysis of 277,896 breast/prostate cancer cases and 901,858 controls identifies 91 newly genome-wide significant loci. Enrichment analysis of pathways and cell types reveals shared genetic backgrounds across said cancers. Focusing on genetically correlated cancers can contribute to enhancing our insights into carcinogenesis.

A Comprehensive Analysis of Cutaneous Melanoma Patients in Greece Based on Multi-Omic Data

Cutaneous melanoma (CM) is the most aggressive type of skin cancer, and it is characterised by high mutational load and heterogeneity. In this study, we aimed to analyse the genomic and transcriptomic profile of primary melanomas from forty-six Formalin-Fixed, Paraffin-Embedded (FFPE) tissues from Greek patients. Molecular analysis for both germline and somatic variations was performed in genomic DNA from peripheral blood and melanoma samples, respectively, exploiting whole exome and targeted sequencing, and transcriptomic analysis. Detailed clinicopathological data were also included in our analyses and previously reported associations with specific mutations were recognised. Most analysed samples (43/46) were found to harbour at least one clinically actionable somatic variant. A subset of samples was profiled at the transcriptomic level, and it was shown that specific melanoma phenotypic states could be inferred from bulk RNA isolated from FFPE primary melanoma tissue. Integrative bioinformatics analyses, including variant prioritisation, differential gene expression analysis, and functional and gene set enrichment analysis by group and per sample, were conducted and molecular circuits that are implicated in melanoma cell programmes were highlighted. Integration of mutational and transcriptomic data in CM characterisation could shed light on genes and pathways that support the maintenance of phenotypic states encrypted into heterogeneous primary tumours.

Ataxia-Telangiectasia Mutated Loss of Heterozygosity in Melanoma

ATM germline pathogenic variants were recently found enriched in high-risk melanoma patients. However, ATM loss of heterozygosity (LOH) has never been investigated in melanoma and, therefore, a causal association with melanoma development has not been established yet. The purpose of this study was to functionally characterize 13 germline ATM variants found in high-risk melanoma patients-and classified by in silico tools as pathogenic, uncertain significance, or benign-using multiple assays evaluating ATM/pATM expression and/or LOH in melanoma tissues and cell lines. We assessed ATM status by Immunohistochemistry (IHC), Western Blot, Whole-Exome Sequencing/Copy Number Variation analysis, and RNA sequencing, supported by Sanger sequencing and microsatellite analyses. For most variants, IHC results matched those obtained with in silico classification and LOH analysis. Two pathogenic variants (p.Ser1135_Lys1192del and p.Ser1993ArgfsTer23) showed LOH and complete loss of ATM activation in melanoma. Two variants of unknown significance (p.Asn358Ile and p.Asn796His) showed reduced expression and LOH, suggestive of a deleterious effect. This study, showing a classic two-hit scenario in a well-known tumor suppressor gene, supports the inclusion of melanoma in the ATM-related cancer spectrum.

A multi-phenotype analysis reveals 19 susceptibility loci for basal cell carcinoma and 15 for squamous cell carcinoma

Basal cell carcinoma and squamous cell carcinoma are the most common skin cancers, and have genetic overlap with melanoma, pigmentation traits, autoimmune diseases, and blood biochemistry biomarkers. In this multi-trait genetic analysis of over 300,000 participants from Europe, Australia and the United States, we reveal 78 risk loci for basal cell carcinoma (19 previously unknown and replicated) and 69 for squamous cell carcinoma (15 previously unknown and replicated). The previously unknown risk loci are implicated in cancer development and progression (e.g. CDKL1), pigmentation (e.g. TPCN2), cardiometabolic (e.g. FADS2), and immune-regulatory pathways for innate immunity (e.g. IFIH1), and HIV-1 viral load modulation (e.g. CCR5). We also report an optimised polygenic risk score for effective risk stratification for keratinocyte cancer in the Canadian Longitudinal Study of Aging (794 cases and 18139 controls), which could facilitate skin cancer surveillance e.g. in high risk subpopulations such as transplantees.

Implementation of individualised polygenic risk score analysis: a test case of a family of four

BACKGROUND

Polygenic risk scores (PRS) have been widely applied in research studies, showing how population groups can be stratified into risk categories for many common conditions. As healthcare systems consider applying PRS to keep their populations healthy, little work has been carried out demonstrating their implementation at an individual level.

CASE PRESENTATION

We performed a systematic curation of PRS sources from established data repositories, selecting 15 phenotypes, comprising an excess of 37 million SNPs related to cancer, cardiovascular, metabolic and autoimmune diseases. We tested selected phenotypes using whole genome sequencing data for a family of four related individuals. Individual risk scores were given percentile values based upon reference distributions among 1000 Genomes Iberians, Europeans, or all samples. Over 96 billion allele effects were calculated in order to obtain the PRS for each of the individuals analysed here.

CONCLUSIONS

Our results highlight the need for further standardisation in the way PRS are developed and shared, the importance of individual risk assessment rather than the assumption of inherited averages, and the challenges currently posed when translating PRS into risk metrics.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Defining novel causal SNPs and linked phenotypes at melanoma-associated loci

A number of genomic regions have been associated with melanoma risk through genome-wide association studies; however, the causal variants underlying the majority of these associations remain unknown. Here, we sequenced either the full locus or the functional regions including exons of 19 melanoma-associated loci in 1959 British melanoma cases and 737 controls. Variant filtering followed by Fisher's exact test analyses identified 66 variants associated with melanoma risk. Sequential conditional logistic regression identified the distinct haplotypes on which variants reside, and massively parallel reporter assays provided biological insights into how these variants influence gene function. We performed further analyses to link variants to melanoma risk phenotypes and assessed their association with melanoma-specific survival. Our analyses replicate previously known associations in the melanocortin 1 receptor (MC1R) and tyrosinase (TYR) loci, while identifying novel potentially causal variants at the MTAP/CDKN2A and CASP8 loci. These results improve our understanding of the architecture of melanoma risk and outcome.

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

Simple Summary

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

Abstract

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

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

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

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

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

Tumor Microenvironment Profiling Identifies Prognostic Signatures and Suggests Immunotherapeutic Benefits in Neuroblastoma

The tumor microenvironment (TME) influences disease initiation and progression. Cross-talks of cells within TME can affect the efficacy of immunotherapies. However, a precise, concise, and comprehensive TME landscape in neuroblastoma (NB) has not been established. Here, we profiled the TME landscape of 498 NB-related patients on a self-curated gene list and identified three prognostic TMEsubgroups. The differentially expressed genes in these three TMEsubgroups were used to construct a genetic signature of the TME landscape and characterize three GeneSubgroups. The subgroup with the worst overall survival prognosis, the TMEsubgroup/GeneSubgroup3, lacked immune cell infiltration and received the highest scores of MYCN- and ALK-related signatures and lowest scores of immune pathways. Additionally, we found that the GeneSubgroup3 might be benefited from anti-GD2 instead of anti-PD-1 therapy. We further created a 48-gene signature, the TMEscore, to infer prognosis and validated it in three independent NB cohorts and a pan-cancer cohort of 9,460 patients. We did RNA-seq on 16 samples and verified that TMEscore was higher in patients with stage 3/4 than stage 1/2 diseases. The TMEscore could also predict responses for several immunotherapies. After adding clinical features, we found that the nomogram-based score system, the TMEIndex, surpassed the current risk system at predicting survivals. Our analysis explained TME at the transcriptome level and paved the way for immunotherapies in NB.

First international workshop of the ATM and cancer risk group (4-5 December 2019)

The first International Workshop of the ATM and Cancer Risk group focusing on the role of Ataxia-Telangiectasia Mutated (ATM) gene in cancer was held on December 4 and 5, 2019 at Institut Curie in Paris, France. It was motivated by the fact that germline ATM pathogenic variants have been found to be associated with different cancer types. However, due to the lack of precise age-, sex-, and site-specific risk estimates, no consensus on management guidelines for variant carriers exists, and the clinical utility of ATM variant testing is uncertain. The meeting brought together epidemiologists, geneticists, biologists and clinicians to review current knowledge and on-going challenges related to ATM and cancer risk. This report summarizes the meeting sessions content that covered the latest results in family-based and population-based studies, the importance of accurate variant classification, the effect of radiation exposures for ATM variant carriers, and the characteristics of ATM-deficient tumors. The report concludes that ATM variant carriers outside of the context of Ataxia-Telangiectasia may benefit from effective cancer risk management and therapeutic strategies and that efforts to set up large-scale studies in the international framework to achieve this goal are necessary.

Citrus-Gene interaction and melanoma risk in the UK Biobank

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

Integration of DNA Microarray with Clinical and Genomic Data

DNA microarrays have been widely employed to understand cancer development. This technology is able to measure expression levels of a large numbers of genes or to genotype multiple regions of a genome in a massively parallel experiment. In addition, the detection of methylation patterns and gene copy number variations are also performed. Clinicians began to apply these findings in personalized medicine for the selection of cancer therapy according to the individual's cancer genomic profile. Because cancer is a complex disease it is of great value to integrate microarray data with genomic and clinical data. Here, we presented an overview of DNA microarray technology and discuss about benefits and challenging of microarray data integration.

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

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

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

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

Association between a 46-SNP Polygenic Risk Score and melanoma risk in Dutch patients with familial melanoma

BACKGROUND

Familial clustering of melanoma suggests a shared genetic predisposition among family members, but only 10%-40% of familial cases carry a pathogenic variant in a known high-risk melanoma susceptibility gene. We investigated whether a melanoma-specific Polygenic Risk Score (PRS) is associated with melanoma risk in patients with genetically unexplained familial melanoma.

METHODS

Dutch familial melanoma cases (n=418) were genotyped for 46 SNPs previously identified as independently associated with melanoma risk. The 46-SNP PRS was calculated and standardised to 3423 healthy controls (sPRS) and the association between PRS and melanoma risk was modelled using logistic regression. Within the case series, possible differences were further explored by investigating the PRS in relation to (1) the number of primary melanomas in a patient and (2) the extent of familial clustering of melanoma.

RESULTS

The PRS was significantly associated with melanoma risk, with a per-SD OR of 2.12 (95% CI 1.90 to 2.35, p<0.001), corresponding to a 5.70-fold increased risk (95% CI 3.93 to 8.28) when comparing the top 90th to the middle 40-60th PRS percentiles. The mean PRS was significantly higher in cases with multiple primary melanomas than in cases with a single melanoma (sPRS 1.17 vs 0.71, p=0.001). Conversely, cases from high-density melanoma families had a lower (but non-significant) mean PRS than cases from low-density families (sPRS 0.60 vs 0.94, p=0.204).

CONCLUSION

Our work underlines the significance of a PRS in determining melanoma susceptibility and encourages further exploration of the diagnostic value of a PRS in genetically unexplained melanoma families.

A pleiotropic ATM variant (rs1800057 C>G) is associated with risk of multiple cancers

ATM (ataxia-telangiectasia mutated) is an important cell-cycle checkpoint kinase required for cellular response to DNA damage. Activated by DNA double strand breaks, ATM regulates the activities of many downstream proteins involved in various carcinogenic events. Therefore, ATM or its genetic variants may have a pleiotropic effect in cancer development. We conducted a pleiotropic analysis to evaluate associations between genetic variants of ATM and risk of multiple cancers. With genotyping data extracted from previously published genome-wide association studies of various cancers, we performed multivariate logistic regression analysis, followed by a meta-analysis for each cancer site, to identify cancer risk-associated single-nucleotide polymorphisms (SNPs). In the ASSET two-sided analysis, we found that two ATM SNPs were significantly associated with risk of multiple cancers. One tagging SNP (rs1800057 C>G) was associated with risk of multiple cancers (two-sided P=5.27×10 -7). Because ATM rs1800057 is a missense variant, we also explored the intermediate phenotypes through which this variant may confer risk of multiple cancers and identified a possible immune-mediated effect of this variant. Our findings indicate that genetic variants of ATM may have a pleiotropic effect on cancer risk and thus provide an important insight into common mechanisms of carcinogenesis.

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

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

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

A genome-wide association study identifies novel gene associations with facial skin wrinkling and mole count in Latin Americans

BACKGROUND

Genome-wide association studies (GWASs) have identified genes influencing skin ageing and mole count in Europeans, but little is known about the relevance of these (or other genes) in non-Europeans.

OBJECTIVES

To conduct a GWAS for facial skin ageing and mole count in adults < 40 years old, of mixed European, Native American and African ancestry, recruited in Latin America.

METHODS

Skin ageing and mole count scores were obtained from facial photographs of over 6000 individuals. After quality control checks, three wrinkling traits and mole count were retained for genetic analyses. DNA samples were genotyped with Illumina's HumanOmniExpress chip. Association testing was performed on around 8 703 729 single-nucleotide polymorphisms (SNPs) across the autosomal genome.

RESULTS

Genome-wide significant association was observed at four genome regions: two were associated with wrinkling (in 1p13·3 and 21q21·2), one with mole count (in 1q32·3) and one with both wrinkling and mole count (in 5p13·2). Associated SNPs in 5p13·2 and in 1p13·3 are intronic within SLC45A2 and VAV3, respectively, while SNPs in 1q32·3 are near the SLC30A1 gene, and those in 21q21·2 occur in a gene desert. Analyses of SNPs in IRF4 and MC1R are consistent with a role of these genes in skin ageing.

CONCLUSIONS

We replicate the association of wrinkling with variants in SLC45A2, IRF4 and MC1R reported in Europeans. We identify VAV3 and SLC30A1 as two novel candidate genes impacting on wrinkling and mole count, respectively. We provide the first evidence that SLC45A2 influences mole count, in addition to variants in this gene affecting melanoma risk in Europeans.

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

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

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.

Pan-cancer study detects genetic risk variants and shared genetic basis in two large cohorts

Deciphering the shared genetic basis of distinct cancers has the potential to elucidate carcinogenic mechanisms and inform broadly applicable risk assessment efforts. Here, we undertake genome-wide association studies (GWAS) and comprehensive evaluations of heritability and pleiotropy across 18 cancer types in two large, population-based cohorts: the UK Biobank (408,786 European ancestry individuals; 48,961 cancer cases) and the Kaiser Permanente Genetic Epidemiology Research on Adult Health and Aging cohorts (66,526 European ancestry individuals; 16,001 cancer cases). The GWAS detect 21 genome-wide significant associations independent of previously reported results. Investigations of pleiotropy identify 12 cancer pairs exhibiting either positive or negative genetic correlations; 25 pleiotropic loci; and 100 independent pleiotropic variants, many of which are regulatory elements and/or influence cross-tissue gene expression. Our findings demonstrate widespread pleiotropy and offer further insight into the complex genetic architecture of cross-cancer susceptibility.

Polymorphism rs1801516 (G > A) in the ATM gene is not associated with overall cancer risk: an updated meta-analysis

OBJECTIVE

The ataxia telangiectasia mutated (ATM) gene contains a functional single nucleotide polymorphism (SNP) rs1801516 (G > A) that may be associated with cancer risk. This meta-analysis aimed to interrogate the relationship between rs1801516 and cancer occurrence and disease etiology.

METHODS

We retrieved and identified the available case-control studies that met the inclusion criteria from the PubMed, Web of Science, and Embase databases. Odds ratio (OR) and 95% confidence intervals (CIs) were used to measure the association between rs1801516 and cancer risk. Additionally, we performed sensitivity, subgroup, and publication bias analyses.

RESULTS

After inclusion criteria were met, the meta-analysis included 29 studies, with 9,453 cancer patients (cases) and 14,646 controls. No association was found between rs1801516 and cancer risk (pooled OR = 0.911; 95% CI, 0.740-1.123). Concordantly, no association was found between rs1801516 and cancer risk after subgroup analysis by source of controls, cancer type, or ethnicity, which confirmed the finding of the dominant model that this SNP is not involved in the occurrence of cancer.

CONCLUSIONS

Through this meta-analysis, we found no association between rs1801516 and cancer occurrence as a risk factor. These data provide useful information for future case-control studies on cancer etiology.

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.

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

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

Evolution of approaches to identify melanoma missing heritability

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

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10^–7, minor allele frequency of 0.2–0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10⁻⁶), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10^–4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Insights into Genetic Susceptibility to Melanoma by Gene Panel Testing: Potential Pathogenic Variants in ACD, ATM, BAP1, and POT1

The contribution of recently established or candidate susceptibility genes to melanoma missing heritability has yet to be determined. Multigene panel testing could increase diagnostic yield and better define the role of candidate genes. We characterized 273 CDKN2A/ARF and CDK4-negative probands through a custom-designed targeted gene panel that included CDKN2A/ARF, CDK4, ACD, BAP1, MITF, POT1, TERF2IP, ATM, and PALB2. Co-segregation, loss of heterozygosity (LOH)/protein expression analysis, and splicing characterization were performed to improve variant classification. We identified 16 (5.9%) pathogenic and likely pathogenic variants in established high/medium penetrance cutaneous melanoma susceptibility genes (BAP1, POT1, ACD, MITF, and TERF2IP), including two novel variants in BAP1 and 4 in POT1. We also found four deleterious and five likely deleterious variants in ATM (3.3%). Thus, including potentially deleterious variants in ATM increased the diagnostic yield to about 9%. Inclusion of rare variants of uncertain significance would increase the overall detection yield to 14%. At least 10% of melanoma missing heritability may be explained through panel testing in our population. To our knowledge, this is the highest frequency of putative ATM deleterious variants reported in melanoma families, suggesting a possible role in melanoma susceptibility, which needs further investigation.

The genetics of human ageing

The past two centuries have witnessed an unprecedented rise in human life expectancy. Sustaining longer lives with reduced periods of disability will require an understanding of the underlying mechanisms of ageing, and genetics is a powerful tool for identifying these mechanisms. Large-scale genome-wide association studies have recently identified many loci that influence key human ageing traits, including lifespan. Multi-trait loci have been linked with several age-related diseases, suggesting shared ageing influences. Mutations that drive accelerated ageing in prototypical progeria syndromes in humans point to an important role for genome maintenance and stability. Together, these different strands of genetic research are highlighting pathways for the discovery of anti-ageing interventions that may be applicable in humans.

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

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

The effects of DNA repair polymorphisms on chromosome aberrations in the population of Kazakhstan

Purpose: To analyze the effects of DNA repair polymorphism and other factors on the frequency chromosome aberrations in an irradiated cohort of subjects living around the Semipalatinsk nuclear test site and non-exposed group of subjects from ecologically favorable zones of Kazakhstan.Materials and methods: Blood samples were collected in the rural areas of the East Kazakhstan district around the Semipalatinsk nuclear test site and ecologically favorable zones of Almaty region of Kazakhstan. Chromosome aberrations in the fresh and cryopreserved peripheral blood lymphocyte cultures were analyzed by Giemsa staining. Single nucleotide polymorphisms at eight DNA repair genes (XRCC1 rs1799782, XRCC1 rs25487, XRCC3 rs861539, ATM rs1801516, XPD rs1799793, XPD rs13181, APEX1 rs1130409, and hOGG1 rs1052133) were determined by PCR-RFLP method.Results: The age of donors and smoking significantly affected the frequency of chromosome aberrations among the irradiated and control subjects. In the irradiated and control cohorts, the frequency of chromosome aberrations was significantly increased in the heterozygous ATM rs1801516 (1853 Asp/Asn) individuals; for the rest of the loci no significant associations between polymorphism and the frequency of chromosome aberrations were detected.Conclusions: The age of donors, smoking, and the ATM rs1801516 polymorphism significantly affect the frequency of chromosome aberrations among individuals inhabiting contaminated area around the Semipalatinsk nuclear weapon test site, as well as among those inhabiting ecologically favorable zones of Kazakhstan.

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.

Combined analysis of keratinocyte cancers identifies novel genome-wide loci

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

Long telomeres and cancer risk: the price of cellular immortality

The distribution of telomere length in humans is broad, but it has finite upper and lower boundaries. Growing evidence shows that there are disease processes that are caused by both short and long telomere length extremes. The genetic basis of these short and long telomere syndromes may be linked to mutations in the same genes, such as the telomerase reverse transcriptase (TERT), but through differential effects on telomere length. Short telomere syndromes have a predominant degenerative phenotype marked by organ failure that most commonly manifests as pulmonary fibrosis and are associated with a relatively low cancer incidence. In contrast, insights from studies of cancer-prone families as well as genome-wide association studies (GWAS) have identified both rare and common variants that lengthen telomeres as being strongly associated with cancer risk. We have hypothesized that these cancers represent a long telomere syndrome that is associated with a high penetrance of cutaneous melanoma and chronic lymphocytic leukemia. In this Review, we will synthesize the clinical and human genetic observations with data from mouse models to define the role of telomeres in cancer etiology and biology.

Immunomodulatory germline variation associated with the development of multiple primary melanoma (MPM)

Multiple primary melanoma (MPM) has been associated with a higher 10-year mortality risk compared to patients with single primary melanoma (SPM). Given that 3–8% of patients with SPM develop additional primary melanomas, new markers predictive of MPM risk are needed. Based on the evidence that the immune system may regulate melanoma progression, we explored whether germline genetic variants controlling the expression of 41 immunomodulatory genes modulate the risk of MPM compared to patients with SPM or healthy controls. By genotyping these 41 variants in 977 melanoma patients, we found that rs2071304, linked to the expression of SPI1, was strongly associated with MPM risk reduction (OR = 0.60; 95% CI = 0.45–0.81; p = 0.0007) when compared to patients with SPM. Furthermore, we showed that rs6695772, a variant affecting expression of BATF3, is also associated with MPM-specific survival (HR = 3.42; 95% CI = 1.57–7.42; p = 0.0019). These findings provide evidence that the genetic variation in immunomodulatory pathways may contribute to the development of secondary primary melanomas and also associates with MPM survival. The study suggests that inherited host immunity may play an important role in MPM development.

Genome-wide association studies and polygenic risk scores for skin cancer: clinically useful yet?

BACKGROUND

Genome-wide association studies (GWAS) have identified thousands of susceptibility variants, although most have been associated with small individual risk estimates that offer little predictive value. However, combining multiple variants into polygenic risk scores (PRS) may be more informative. Multiple studies have developed PRS composed of GWAS-identified variants for cutaneous cancers. This review highlights data from these studies.

OBJECTIVES

To review published GWAS and PRS studies for melanoma, cutaneous squamous cell carcinoma (cSCC) and basal cell carcinoma (BCC), and discuss their potential clinical utility.

METHODS

We searched PubMed and the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue to identify relevant studies.

RESULTS

Results from 21 GWAS (11 melanoma, 3 cSCC, 7 BCC) and 11 PRS studies are summarized. Six loci in pigmentation genes overlap between these three cancers (ASIP/RALY, IRF4, MC1R, OCA2, SLC45A2 and TYR). Additional loci overlap for cSCC/BCC and BCC/melanoma, but no other loci are shared between cSCC and melanoma. PRS for melanoma show roughly two-to-threefold increases in risk and modest improvements in risk prediction (2-7% increases). PRS are associated with twofold and threefold increases in risk of cSCC and BCC, respectively, with small improvements (2% increase) in predictive ability.

CONCLUSIONS

Existing data indicate that PRS may offer small, but potentially meaningful, improvements to risk prediction. Additional research is needed to clarify the potential utility of PRS in cutaneous carcinomas. Clinical translation will require well-powered validation studies incorporating known risk factors to evaluate PRS as tools for screening. What's already known about this topic? Over 50 susceptibility loci for melanoma, basal cell carcinoma (BCC) and cutaneous squamous cell carcinoma (cSCC) have been identified in genome-wide association studies (GWAS). Polygenic risk scores (PRS) using variants identified from GWAS have also been developed for melanoma, BCC and cSCC, and investigated with respect to clinical risk prediction. What does this study add? This review provides an overview of GWAS findings and the potential clinical utility of PRS for melanoma, BCC and cSCC.