SLC45A2: a novel malignant melanoma-associated gene

Mapping chromatin interactions at melanoma susceptibility loci and cell-type specific dataset integration uncovers distant gene targets of cis-regulation

Genome-wide association studies (GWAS) of melanoma risk have identified 68 independent signals at 54 loci. For most loci, specific functional variants and their respective target genes remain to be established. Capture-HiC is an assay that links fine-mapped risk variants to candidate target genes by comprehensively mapping cell-type specific chromatin interactions. We performed a melanoma GWAS region-focused capture-HiC assay in human primary melanocytes to identify physical interactions between fine-mapped risk variants and potential causal melanoma susceptibility genes. Overall, chromatin interaction data alone nominated potential causal genes for 61 of the 68 melanoma risk signals, identifying many candidates beyond those reported by previous studies. We further integrated these data with cell-type specific epigenomic (chromatin state, accessibility), gene expression (eQTL/TWAS), DNA methylation (meQTL/MWAS), and massively parallel reporter assay (MPRA) data to prioritize potentially cis-regulatory variants and their respective candidate gene targets. From the set of fine-mapped variants across these loci, we identified 140 prioritized candidate causal variants linked to 195 candidate genes at 42 risk signals. In addition, we developed an integrative scoring system to facilitate candidate gene prioritization, integrating melanocyte and melanoma datasets. Notably, at several GWAS risk signals we observed long-range chromatin connections (500 kb to >1 Mb) with distant candidate target genes. We validated several such cis-regulatory interactions using CRISPR inhibition, providing evidence for known cancer driver genes MDM4 and CBL, as well as the SRY-box transcription factor SOX4, as likely melanoma risk genes.

Modeling recent positive selection using identity-by-descent segments

Recent positive selection can result in an excess of long identity-by-descent (IBD) haplotype segments overlapping a locus. The statistical methods that we propose here address three major objectives in studying selective sweeps: scanning for regions of interest, identifying possible sweeping alleles, and estimating a selection coefficient s. First, we implement a selection scan to locate regions with excess IBD rates. Second, we estimate the allele frequency and location of an unknown sweeping allele by aggregating over variants that are more abundant in an inferred outgroup with excess IBD rate versus the rest of the sample. Third, we propose an estimator for the selection coefficient and quantify uncertainty using the parametric bootstrap. Comparing against state-of-the-art methods in extensive simulations, we show that our methods are more precise at estimating s when s≥0.015. We also show that our 95% confidence intervals contain s in nearly 95% of our simulations. We apply these methods to study positive selection in European ancestry samples from the Trans-Omics for Precision Medicine project. We analyze eight loci where IBD rates are more than four standard deviations above the genome-wide median, including LCT where the maximum IBD rate is 35 standard deviations above the genome-wide median. Overall, we present robust and accurate approaches to study recent adaptive evolution without knowing the identity of the causal allele or using time series data.

Unveiling the genetic landscape of hereditary melanoma: From susceptibility to surveillance

The multifactorial etiology underlying melanoma development involves an array of genetic, phenotypic, and environmental factors. Genetic predisposition for melanoma is further influenced by the complex interplay between high-, medium-, and low-penetrance genes, each contributing to varying degrees of susceptibility. Within this network, high-penetrance genes, including CDKN2A, CDK4, BAP1, and POT1, are linked to a pronounced risk for disease, whereas medium- and low-penetrance genes, such as MC1R, MITF, and others, contribute only moderately to melanoma risk. Notably, these genetic factors not only heighten the risk of melanoma but may also increase susceptibility towards internal malignancies, such as pancreatic cancer, renal cell cancer, or neural tumors. Genetic testing and counseling hold paramount importance in the clinical context of suspected hereditary melanoma, facilitating risk assessment, personalized surveillance strategies, and informed decision-making. As our understanding of the genomic landscape deepens, this review paper aims to comprehensively summarize the genetic underpinnings of hereditary melanoma, as well as current screening and management strategies for the disease.

What Makes a Mimic? Orange, Red, and Black Color Production in the Mimic Poison Frog (Ranitomeya imitator)

Aposematic organisms rely on their conspicuous appearance to signal that they are defended and unpalatable. Such phenotypes are strongly tied to survival and reproduction. Aposematic colors and patterns are highly variable; however, the genetic, biochemical, and physiological mechanisms producing this conspicuous coloration remain largely unidentified. Here, we identify genes potentially affecting color variation in two color morphs of Ranitomeya imitator: the orange-banded Sauce and the redheaded Varadero morphs. We examine gene expression in black and orange skin patches from the Sauce morph and black and red skin patches from the Varadero morph. We identified genes differentially expressed between skin patches, including those that are involved in melanin synthesis (e.g. mlana, pmel, tyrp1), iridophore development (e.g. paics, ppat, ak1), pteridine synthesis (e.g. gch1, pax3-a, xdh), and carotenoid metabolism (e.g. dgat2, rbp1, scarb2). In addition, using weighted correlation network analysis, we identified the top 50 genes with high connectivity from the most significant network associated with gene expression differences between color morphs. Of these 50 genes, 13 were known to be related to color production (gch1, gmps, gpr143, impdh1, mc1r, pax3-a, pax7, ppat, rab27a, rlbp1, tfec, trpm1, xdh).

Identifying in vivo genetic dependencies of melanocyte and melanoma development

The advent of large-scale sequencing in both development and disease has identified large numbers of candidate genes that may be linked to important phenotypes. Validating the function of these candidates in vivo is challenging, due to low efficiency and low throughput of most model systems. We have developed a rapid, scalable system for assessing the role of candidate genes using zebrafish. We generated transgenic zebrafish in which Cas9 was knocked-in to the endogenous mitfa locus, a master transcription factor of the melanocyte lineage. We used this system to identify both cell-autonomous and non-cell autonomous regulators of normal melanocyte development. We then applied this to the melanoma setting to demonstrate that loss of genes required for melanocyte survival can paradoxically promote more aggressive phenotypes, highlighting that in vitro screens can mask in vivo phenotypes. Our high-efficiency genetic approach offers a versatile tool for exploring developmental processes and disease mechanisms that can readily be applied to other cell lineages.

Forensic DNA Phenotyping: Genes and Genetic Variants for Eye Color Prediction

In recent decades, the use of genetic polymorphisms related to specific phenotypes, such as eye color, has greatly contributed to the development of the research field called forensic DNA phenotyping (FDP), enabling the investigators of crime cases to reduce the number of suspects, making their work faster and more precise. Eye color is a polygenic phenotype, and many genetic variants have been highlighted, with the major contributor being the HERC2-OCA2 locus, where many single nucleotide variations (SNPs) were identified. Interestingly, the HERC2-OCA2 locus, containing the intronic SNP rs12913832, the major eye color determinant, shows a high level of evolutionary conservation across many species of vertebrates. Currently, there are some genetic panels to predict eye color by genomic DNA analysis, even if the exact role of the SNP variants in the formation of eye color is still poorly understood, with a low level of predictivity in the so-called intermediate eye color. Many variants in OCA2, HERC2, and other genes lie in introns or correspond to synonymous variants, highlighting greater complexity in the mechanism of action of such genes than a simple missense variation. Here, we show the main genes involved in oculocutaneous pigmentation and their structural and functional features, as well as which genetic variants show the highest level of eye color predictivity in currently used FDP assays. Despite the great recent advances and impact of FDP in criminal cases, it is necessary to enhance scientific research to better understand the mechanism of action behind each genetic variant involved in eye color, with the goal of obtaining higher levels of prediction.

The relationship between single nucleotide polymorphisms and skin cancer susceptibility: A systematic review and network meta-analysis

Background

Single nucleotide polymorphisms (SNPs) interfere with the function of certain genes and thus may influence the probability of skin cancer. The correlation between SNPs and skin cancer (SC) lacks statistical power, however. Therefore, the purpose of this study was to identify the gene polymorphisms involved in skin cancer susceptibility using network meta-analysis and to determine the relationship between SNPs and SC risk.

Methods

PubMed, Embase, and Web of Science were searched for articles including "SNP" and different types of SC as keywords between January 2005 and May 2022. The Newcastle-Ottawa Scale was used to assess bias judgments. The odds ratio (ORs) and their 95% confidence intervals (CIs) were determined to estimate heterogeneity within and between studies. Meta-analysis and network meta-analysis were carried out to identify the SNPs associated with SC. The P-score of each SNP was compared to obtain the rank of probability. Subgroup analyses were performed by cancer type.

Results

A total of 275 SNPs from 59 studies were included in the study. Two subgroup SNP networks using the allele model and dominant model were analyzed. The alternative alleles of rs2228570 (FokI) and rs13181 (ERCC2) were the first-ranking SNPs in both subgroups one and two of the allele model, respectively. The homozygous dominant genotype and heterozygous genotype of rs475007 in subgroup one and the homozygous recessive genotype of rs238406 in subgroup two were most likely to be associated with skin cancer based on the dominant model.

Conclusions

According to the allele model, SNPs FokI rs2228570 and ERCC2 rs13181 and, according to the dominant model, SNPs MMP1 rs475007 and ERCC2 rs238406 are closely linked to SC risk.

The potential impact of melanosomal pH and metabolism on melanoma

Melanin is synthesized in melanocytes and is transferred into keratinocytes to block the effects of ultraviolet (UV) radiation and is important for preventing skin cancers including melanoma. However, it is known that after melanomagenesis and melanoma invasion or metastases, melanin synthesis still occurs. Since melanoma cells are no longer involved in the sun tanning process, it is unclear why melanocytes would maintain melanin synthesis after melanomagenesis has occurred. Aside from blocking UV-induced DNA mutation, melanin may provide other metabolic functions that could benefit melanoma. In addition, studies have suggested that there may be a selective advantage to melanin synthesis in melanoma; however, mechanisms regulating melanin synthesis outside the epidermis or hair follicle is unknown. We will discuss how melanosomal pH controls melanin synthesis in melanocytes and how melanosomal pH control of melanin synthesis might function in melanoma. We will also discuss potential reasons why melanin synthesis might be beneficial for melanoma cellular metabolism and provide a rationale for why melanin synthesis is not limited to benign melanocytes.

A systematic review of skin ageing genes: gene pleiotropy and genes on the chromosomal band 16q24.3 may drive skin ageing

Skin ageing is the result of intrinsic genetic and extrinsic lifestyle factors. However, there is no consensus on skin ageing phenotypes and ways to quantify them. In this systematic review, we first carefully identified 56 skin ageing phenotypes from multiple literature sources and sought the best photo-numeric grading scales to evaluate them. Next, we conducted a systematic review on all 44 Genome-wide Association Studies (GWAS) on skin ageing published to date and identified genetic risk factors (2349 SNPs and 366 genes) associated with skin ageing. We identified 19 promising SNPs found to be significantly (p-Value < 1E-05) associated with skin ageing phenotypes in two or more independent studies. Here we show, using enrichment analyses strategies and gene expression data, that (1) pleiotropy is a recurring theme among skin ageing genes, (2) SNPs associated with skin ageing phenotypes are mostly located in a small handful of 44 pleiotropic and hub genes (mostly on the chromosome band 16q24.3) and 32 skin colour genes. Since numerous genes on the chromosome band 16q24.3 and skin colour genes show pleiotropy, we propose that (1) genes traditionally identified to contribute to skin colour have more than just skin pigmentation roles, and (2) further progress towards understand the development of skin pigmentation requires understanding the contributions of genes on the chromosomal band 16q24.3. We anticipate our systematic review to serve as a hub to locate primary literature sources pertaining to the genetics of skin ageing and to be a starting point for more sophisticated work examining pleiotropic genes, hub genes, and skin ageing phenotypes.

Oncogenic Activity of Solute Carrier Family 45 Member 2 and Alpha-Methylacyl-Coenzyme A Racemase Gene Fusion Is Mediated by Mitogen-Activated Protein Kinase

Chromosome rearrangement is one of the hallmarks of human malignancies. Gene fusion is one of the consequences of chromosome rearrangements. In this report, we show that gene fusion between solute carrier family 45 member 2 (SLC45A2) and alpha-methylacyl-coenzyme A racemase (AMACR) occurs in eight different types of human malignancies, with frequencies ranging from 45% to 97%. The chimeric protein is translocated to the lysosomal membrane and activates the extracellular signal-regulated kinase signaling cascade. The fusion protein promotes cell growth, accelerates migration, resists serum starvation-induced cell death, and is essential for cancer growth in mouse xenograft cancer models. Introduction of SLC45A2-AMACR into the mouse liver using a sleeping beauty transposon system and somatic knockout of phosphatase and TENsin homolog (Pten) generated spontaneous liver cancers within a short period. Conclusion: The gene fusion between SLC45A2 and AMACR may be a driving event for human liver cancer development.

Canine coat pigmentation genetics: a review

Our understanding of canine coat colour genetics and the associated health implications is developing rapidly. To date, there are 15 genes with known roles in canine coat colour phenotypes. Many coat phenotypes result from complex and/or epistatic genetic interactions among variants within and between loci, some of which remain unidentified. Some genes involved in canine pigmentation have been linked to aural, visual and neurological impairments. Consequently, coat pigmentation in the domestic dog retains considerable ethical and economic interest. In this paper we discuss coat colour phenotypes in the domestic dog, the genes and variants responsible for these phenotypes and any proven coat colour-associated health effects.

caAtlas: An immunopeptidome atlas of human cancer

Comprehensive characterization of tumor antigens is essential for the design of cancer immunotherapies, and mass spectrometry (MS)-based immunopeptidomics enables high-throughput identification of major histocompatibility complex (MHC)-bound peptide antigens in vivo. Here we construct an immunopeptidome atlas of human cancer through an extensive collection of 43 published immunopeptidomic datasets and standardized analysis of 81.6 million MS/MS spectra using an open search engine. Our analysis greatly expands the current knowledge of MHC-bound antigens, including an unprecedented characterization of post-translationally modified antigens and their cancer-association. We also perform systematic analysis of cancer-testis antigens, cancer-associated antigens, and neoantigens. We make all these data together with annotated MS/MS spectra supporting identification of each antigen in an easily browsable web portal named cancer antigen atlas (caAtlas). caAtlas provides a central resource for the selection and prioritization of MHC-bound peptides for in vitro HLA binding assay and immunogenicity testing, which will pave the way to eventual development of cancer immunotherapies.

•

Extensive collection of 43 immunopeptidomic datasets with 1018 samples

•

Standardized and rigorous identification of HLA-bound peptides, including PTM peptides

•

Comprehensive annotation of CT antigens and cancer-associated antigens

•

User-friendly data dissemination through the caAtlas web portal

Computational Investigation of the pH Dependence of Stability of Melanosome Proteins: Implication for Melanosome formation and Disease

Intravesicular pH plays a crucial role in melanosome maturation and function. Melanosomal pH changes during maturation from very acidic in the early stages to neutral in late stages. Neutral pH is critical for providing optimal conditions for the rate-limiting, pH-sensitive melanin-synthesizing enzyme tyrosinase (TYR). This dramatic change in pH is thought to result from the activity of several proteins that control melanosomal pH. Here, we computationally investigated the pH-dependent stability of several melanosomal membrane proteins and compared them to the pH dependence of the stability of TYR. We confirmed that the pH optimum of TYR is neutral, and we also found that proteins that are negative regulators of melanosomal pH are predicted to function optimally at neutral pH. In contrast, positive pH regulators were predicted to have an acidic pH optimum. We propose a competitive mechanism among positive and negative regulators that results in pH equilibrium. Our findings are consistent with previous work that demonstrated a correlation between the pH optima of stability and activity, and they are consistent with the expected activity of positive and negative regulators of melanosomal pH. Furthermore, our data suggest that disease-causing variants impact the pH dependence of melanosomal proteins; this is particularly prominent for the OCA2 protein. In conclusion, melanosomal pH appears to affect the activity of multiple melanosomal proteins.

Measurement of Melanin Metabolism in Live Cells by [U-13C]-L-Tyrosine Fate Tracing Using Liquid Chromatography-Mass Spectrometry

Melanin synthesis occurs within a specialized organelle called the melanosome. Traditional methods for measuring melanin levels rely on the detection of chemical degradation products of melanin by high-performance liquid chromatography. Although these methods are robust, they are unable to distinguish between melanin synthesis and degradation and are best suited to measure melanin changes over long periods of time. We developed a method that actively measures both eumelanin and pheomelanin synthesis by fate tracing [U-13C] L-tyrosine using liquid chromatography-mass spectrometry. Using this method, we confirmed the previous reports of the differences in melanin synthesis between melanocytes derived from individuals with different skin colors and MC1R genotype and uncovered new information regarding the differential de novo synthesis of eumelanin and pheomelanin, also called mixed melanogenesis. We also revealed that distinct mechanisms that alter melanosomal pH differentially induce new eumelanin and pheomelanin synthesis. Finally, we revealed that the synthesis of L-3,4-dihydroxyphenylalanine, an important metabolite of L-tyrosine, is differentially controlled by multiple factors. Because L-tyrosine fate tracing is compatible with untargeted liquid chromatography-mass spectrometry‒based metabolomics, this approach enables the broad measurement of cellular metabolism in combination with melanin metabolism, and we anticipate that this approach will shed new light on multiple mechanisms of melanogenesis.

Membrane Transporters and Channels in Melanoma

Recent research has revealed that ion channels and transporters can be important players in tumor development, progression, and therapy resistance in melanoma. For example, members of the ABC family were shown to support cancer stemness-like features in melanoma cells, while several members of the TRP channel family were reported to act as tumor suppressors.Also, many transporter proteins support tumor cell viability and thus suppress apoptosis induction by anticancer therapy. Due to the high number of ion channels and transporters and the resulting high complexity of the field, progress in understanding is often focused on single molecules and is in total rather slow. In this review, we aim at giving an overview about a broad subset of ion transporters, also illustrating some aspects of the field, which have not been addressed in detail in melanoma. In context with the other chapters in this special issue on "Transportome Malfunctions in the Cancer Spectrum," a comparison between melanoma and these tumors will be possible.

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

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

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.

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.

Genetic variants associated with skin photosensitivity in a southern European population from Spain

BACKGROUND/PURPOSE

Recent GWAS studies, mostly performed in populations of North European origin, have identified the genetic loci associated with pigmentation, sun sensitivity, freckling and skin cancer susceptibility. Here, we aimed at addressing the genetic determinants of sunlight sensitivity in Spain, a southern European population.

METHODS

Nine SNPs located in 8 pigmentation-related genes (IRF4, TYR, ASP, HERC2, OCA2, BNC2, SLC24A4 and SLC45A2) were genotyped in 456 Spaniards. Additionally, the complete sequence of the MC1R gene was obtained, testing each nonsynonymous mutation supported by the classification as R or r alleles. A standardised questionnaire was used to collect demographic characteristics, pigmentation and sun sensitivity traits, as well as sun exposure habits.

RESULTS

MC1R R alleles and IRF4 rs12203592 were significantly associated with sunlight sensitivity at the Bonferroni-corrected level (P-value < 4.54 × 10^-3 ). Genetic variants in SLC45A2 (rs16891982) and HERC2 (rs12913832) were also found to be significantly associated with skin photosensitivity in our Spanish sample. Interaction analysis using the MDR method revealed epistatic effects when these four variants were considered together.

CONCLUSION

MC1R, IRF4, HERC2 and SLC45A2 play a significant role in skin sensitivity to sunlight in the Spanish population. Moreover, interaction among these four loci seems to modulate the ability of the skin to respond to UV radiation.

Associations between sun sensitive pigmentary genes and serum prostate specific antigen levels

Background

Melanoma and prostate cancer may share risk factors. This study examined the association between serum PSA levels, which is a risk factor for prostate cancer, and variants in some melanoma-associated pigmentary genes.

Methods

We studied participants, all aged 70+ years, in the Concord Health and Ageing in Men Project who had no history of prostatitis or received treatment for prostate disease (n = 1033). We genotyped variants in MC1R (rs1805007, rs1805008), ASIP (rs4911414, rs1015362), SLC45A2 (rs28777, rs16891982), IRF4 (rs12203592), TYRP1 (rs1408799), TYR (rs1126809, rs1042602), SLC24A2 (rs12896399), and OCA2 (rs7495174). Generalised linear dominant models with Poisson distribution, log link functions and robust variance estimators estimated adjusted percentage differences (%PSA) in mean serum PSA levels (ng/mL) between variant and wildtype (0%PSA = reference) genotypes, adjusting for age, body mass index, serum 25OHD levels and birth regions (Australia or New Zealand (ANZ), Europe or elsewhere).

Results

Serum PSA levels were strongly associated with advancing age and birth regions: mean PSA levels were lower in Europe-born (-29.7%) and elsewhere-born (-11.7%) men than ANZ-born men (reference). Lower %PSA was observed in men with variants in SLC45A2: rs28777 (-19.6;95%CI: -33.5, -2.7), rs16891982 (-17.3;95%CI:-30.4,-1.7) than in wildtype men (reference). There were significant interactions between birth regions and PSA levels in men with variants in MC1R (rs1805007; p-interaction = 0.0001) and ASIP (rs4911414; p-interaction = 0.007). For these genes %PSA was greater in ANZ-born men and lower in Europe- and elsewhere-born men with the variant than it was in wildtype men. In a post hoc analysis, serum testosterone levels were increased in men with MC1R rs1805007 and serum dihydrotestosterone in men with ASIP rs1015362.

Conclusion

Men with SNPs in SLC45A2, who have less sun sensitive skin, have lower PSA levels. Men with SNPs in MC1R and ASIP, who have more sun sensitive skin, and were born in ANZ, have higher PSA levels. Androgens may modify these apparent associations of pigmentary genes and sun exposure with PSA levels.

Impact

PSA levels and possibly prostate cancer risk may vary with sun sensitivity and sun exposure, the effects of which might be modified by androgen levels.

Two-stage genome-wide association study identifies a novel susceptibility locus associated with melanoma

Genome-wide association studies have identified 21 susceptibility loci associated with melanoma. These loci implicate genes affecting pigmentation, nevus count, telomere maintenance, and DNA repair in melanoma risk. Here, we report the results of a two-stage genome-wide association study of melanoma. The stage 1 discovery phase consisted of 4,842 self-reported melanoma cases and 286,565 controls of European ancestry from the 23andMe research cohort and the stage 2 replication phase consisted of 1,804 melanoma cases and 1,026 controls from the University of Texas M.D. Anderson Cancer Center. We performed a combined meta-analysis totaling 6,628 melanoma cases and 287,591 controls. Our study replicates 20 of 21 previously known melanoma-loci and confirms the association of the telomerase reverse transcriptase, TERT, with melanoma susceptibility at genome-wide significance. In addition, we uncover a novel polymorphism, rs187843643 (OR = 1.96; 95% CI = [1.54, 2.48]; P = 3.53 × 10⁻⁸), associated with melanoma. The SNP rs187842643 lies within a noncoding RNA 177kb downstream of BASP1 (brain associated protein-1). We find that BASP1 expression is suppressed in melanoma as compared with benign nevi, providing additional evidence for a putative role in melanoma pathogenesis.

Global skin colour prediction from DNA

Human skin colour is highly heritable and externally visible with relevance in medical, forensic, and anthropological genetics. Although eye and hair colour can already be predicted with high accuracies from small sets of carefully selected DNA markers, knowledge about the genetic predictability of skin colour is limited. Here, we investigate the skin colour predictive value of 77 single-nucleotide polymorphisms (SNPs) from 37 genetic loci previously associated with human pigmentation using 2025 individuals from 31 global populations. We identified a minimal set of 36 highly informative skin colour predictive SNPs and developed a statistical prediction model capable of skin colour prediction on a global scale. Average cross-validated prediction accuracies expressed as area under the receiver-operating characteristic curve (AUC) ± standard deviation were 0.97 ± 0.02 for Light, 0.83 ± 0.11 for Dark, and 0.96 ± 0.03 for Dark-Black. When using a 5-category, this resulted in 0.74 ± 0.05 for Very Pale, 0.72 ± 0.03 for Pale, 0.73 ± 0.03 for Intermediate, 0.87±0.1 for Dark, and 0.97 ± 0.03 for Dark-Black. A comparative analysis in 194 independent samples from 17 populations demonstrated that our model outperformed a previously proposed 10-SNP-classifier approach with AUCs rising from 0.79 to 0.82 for White, comparable at the intermediate level of 0.63 and 0.62, respectively, and a large increase from 0.64 to 0.92 for Black. Overall, this study demonstrates that the chosen DNA markers and prediction model, particularly the 5-category level; allow skin colour predictions within and between continental regions for the first time, which will serve as a valuable resource for future applications in forensic and anthropologic genetics.

Genetic 3'UTR variation is associated with human pigmentation characteristics and sensitivity to sunlight

Sunlight exposure induces signalling pathways leading to the activation of melanin synthesis and tanning response. MicroRNAs (miRNAs) can regulate the expression of genes involved in pigmentation pathways by binding to the complementary sequence in their 3'untranslated regions (3'UTRs). Therefore, 3'UTR SNPs are predicted to modify the ability of miRNAs to target genes, resulting in differential gene expression. In this study, we investigated the role in pigmentation and sun-sensitivity traits, as well as in melanoma susceptibility, of 38 different 3'UTR SNPs from 38 pigmentation-related genes. A total of 869 individuals of Spanish origin (526 melanoma cases and 343 controls) were analysed. The association of genotypic data with pigmentation traits was analysed via logistic regression. Web-based tools for predicting the effect of genetic variants in microRNA-binding sites in 3'UTR gene regions were also used. Seven 3'UTR SNPs showed a potential implication in melanoma risk phenotypes. This association is especially noticeable for two of them, rs2325813 in the MLPH gene and rs752107 in the WNT3A gene. These two SNPs were predicted to disrupt a miRNA-binding site and to impact on miRNA-mRNA interaction. To our knowledge, this is the first time that these two 3'UTR SNPs have been associated with sun-sensitivity traits. We state the potential implication of these SNPs in human pigmentation and sensitivity to sunlight, possibly as a result of changes in the level of gene expression through the disruption of putative miRNA-binding sites.

Genetic predisposition to melanoma

Malignant melanoma is a rare, often fatal form of skin cancer with a complex multigenic etiology. The incidence of melanoma is increasing at an alarming rate. A number of heritable factors contribute to a patient's overall melanoma risk, including response to ultraviolet light, nevus number, and pigmentation characteristics, such as eye and hair color. Approximately 5%-10% of melanoma cases are familial, yet the majority of familial cases lack identifiable germ-line mutations in known susceptibility genes. Additionally, most familial melanomas lack germ-line mutations in genes that are commonly mutated in sporadic melanoma. Candidate and systematic genome-wide association studies have led to an improved understanding of the risk factors for melanoma and the identification of susceptibility genes. In this review, we provide an overview of the major risk factors and known genes implicated in familial melanoma susceptibility.

Genome-wide association study identifies novel susceptibility loci for cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma represents the second most common cutaneous malignancy, affecting 7–11% of Caucasians in the United States. The genetic determinants of susceptibility to cutaneous squamous cell carcinoma remain largely unknown. Here we report the results of a two-stage genome-wide association study of cutaneous squamous cell carcinoma, totalling 7,404 cases and 292,076 controls. Eleven loci reached genome-wide significance (P<5 × 10⁻⁸) including seven previously confirmed pigmentation-related loci: MC1R, ASIP, TYR, SLC45A2, OCA2, IRF4 and BNC2. We identify an additional four susceptibility loci: 11q23.3 CADM1, a metastasis suppressor gene involved in modifying tumour interaction with cell-mediated immunity; 2p22.3; 7p21.1 AHR, the dioxin receptor involved in anti-apoptotic pathways and melanoma progression; and 9q34.3 SEC16A, a putative oncogene with roles in secretion and cellular proliferation. These susceptibility loci provide deeper insight into the pathogenesis of squamous cell carcinoma.

Cutaneous squamous cell carcinoma is the second most common type of skin cancer. In this genome-wide association study, which includes over 7,000 cases, the authors identify 4 new susceptibility loci for this cancer and also provide independent replication of 9 previously reported susceptibility loci.

Biochip-Based Genotyping Assay for Detection of Polymorphisms in Pigmentation Genes Associated with Cutaneous Melanoma

AIMS

The purpose of the study was to develop a new assay for genotyping nine single nucleotide polymorphisms (SNPs) that are known to be associated with melanoma.

METHODS

Two-stage single tube polymerase chain reaction (PCR) followed by hybridization on a biochip was developed and applied in the study.

RESULTS

A total of nine SNPs were selected from five genes: MC1R (rs1805006, rs1805007, rs1805009, rs11547464), HERC2 (rs12913832), OCA2 (rs1800407), SLC45A2 (rs16891982), TYR (rs1393350), and a SNP from the intergenic locus rs12896399 were used for the synthesis of ssDNAs via a single-stage PCR process. The assays were performed on a biochip-based platform that is capable of SNP genotyping via a single reaction-tube PCR, followed by on chip hybridization. We tested 69 DNAs obtained from healthy persons and demonstrated the assays' ability to discriminate all three genotypes for almost all of the SNPs.

CONCLUSIONS

The developed approach proved robust, suggesting that it might be useful for the personalized genotyping of large cohorts of patients.

Sex-specific genetic effects associated with pigmentation, sensitivity to sunlight, and melanoma in a population of Spanish origin

BACKGROUND

Human pigmentation is a polygenic quantitative trait with high heritability. In addition to genetic factors, it has been shown that pigmentation can be modulated by oestrogens and androgens via up- or down-regulation of melanin synthesis. Our aim was to identify possible sex differences in pigmentation phenotype as well as in melanoma association in a melanoma case-control population of Spanish origin.

METHODS

Five hundred and ninety-nine females (316 melanoma cases and 283 controls) and 458 males (234 melanoma cases and 224 controls) were analysed. We genotyped 363 polymorphisms (single nucleotide polymorphisms (SNPs)) from 65 pigmentation gene regions.

RESULTS

When samples were stratified by sex, we observed more SNPs associated with dark pigmentation and good sun tolerance in females than in males (107 versus 75; P = 2.32 × 10(-6)), who were instead associated with light pigmentation and poor sun tolerance. Furthermore, six SNPs in TYR, SILV/CDK2, GPR143, and F2RL1 showed strong differences in melanoma risk by sex (P < 0.01).

CONCLUSIONS

We demonstrate that these genetic variants are important for pigmentation as well as for melanoma risk, and also provide suggestive evidence for potential differences in genetic effects by sex.

Human adaptation and population differentiation in the light of ancient genomes

The influence of positive selection sweeps in human evolution is increasingly debated, although our ability to detect them is hampered by inherent uncertainties in the timing of past events. Ancient genomes provide snapshots of allele frequencies in the past and can help address this question. We combine modern and ancient genomic data in a simple statistic (DAnc) to time allele frequency changes, and investigate the role of drift and adaptation in population differentiation. Only 30% of the most strongly differentiated alleles between Africans and Eurasians changed in frequency during the colonization of Eurasia, but in Europe these alleles are enriched in genic and putatively functional alleles to an extent only compatible with local adaptation. Adaptive alleles—especially those associated with pigmentation—are mostly of hunter-gatherer origin, although lactose persistence arose in a haplotype present in farmers. These results provide evidence for a role of local adaptation in human population differentiation.

Detecting the targets of positive selection in the human genome is challenging. Here, the authors combine modern and ancient genomes to show that alleles strongly differentiated between Africans and Europeans mediated local adaptation in European populations, and were mostly contributed by ancient hunter-gatherers.

Association of Interferon Regulatory Factor-4 Polymorphism rs12203592 With Divergent Melanoma Pathways

BACKGROUND

Solar elastosis and neval remnants are histologic markers characteristic of divergent melanoma pathways linked to differences in age at onset, host phenotype, and sun exposure. However, the association between these pathway markers and newly identified low-penetrance melanoma susceptibility loci remains unknown.

METHODS

In the Genes, Environment and Melanoma (GEM) Study, 2103 Caucasian participants had first primary melanomas that underwent centralized pathology review. For 47 single-nucleotide polymorphisms (SNPs) previously identified as low-penetrant melanoma risk variants, we used multinomial logistic regression to compare melanoma with solar elastosis and melanoma with neval remnants simultaneously to melanoma with neither of these markers, excluding melanomas with both markers. All statistical tests were two-sided.

RESULTS

IRF4 rs12203592 was the only SNP to pass the false discovery threshold in baseline models adjusted for age, sex, and study center. rs12203592*T was associated positively with melanoma with solar elastosis (odds ratio [OR] = 1.47, 95% confidence interval [CI] = 1.18 to 1.82) and inversely with melanoma with neval remnants (OR = 0.65, 95% CI = 0.48 to 0.87) compared with melanoma with neither marker (P global = 3.78 x 10(-08)). Adjusting for phenotypic characteristics and total sun exposure hours did not materially affect rs12203592's associations. Distinct early- and late-onset age distributions were observed in patients with IRF4 rs12203592 [CC] and [TT] genotypes, respectively.

CONCLUSIONS

Our findings suggest a role of IRF4 rs12203592 in pathway-specific risk for melanoma development. We hypothesize that IRF4 rs12203592 could underlie in part the bimodal age distribution reported for melanoma and linked to the divergent pathways.

Identification of Susceptibility Loci for Cutaneous Squamous Cell Carcinoma

We report a genome-wide association study of cutaneous squamous cell carcinoma conducted among non-Hispanic white members of the Kaiser Permanente Northern California health care system. The study includes a genome-wide screen of 61,457 members (6,891 cases and 54,566 controls) genotyped on the Affymetrix Axiom European array and a replication phase involving an independent set of 6,410 additional members (810 cases and 5,600 controls). Combined analysis of screening and replication phases identified 10 loci containing single-nucleotide polymorphisms (SNPs) with P-values < 5 × 10(-8). Six loci contain genes in the pigmentation pathway; SNPs at these loci appear to modulate squamous cell carcinoma risk independently of the pigmentation phenotypes. Another locus contains HLA class II genes studied in relation to elevated squamous cell carcinoma risk following immunosuppression. SNPs at the remaining three loci include an intronic SNP in FOXP1 at locus 3p13, an intergenic SNP at 3q28 near TP63, and an intergenic SNP at 9p22 near BNC2. These findings provide insights into the genetic factors accounting for inherited squamous cell carcinoma susceptibility.

A probabilistic method for testing and estimating selection differences between populations

Human populations around the world encounter various environmental challenges and, consequently, develop genetic adaptations to different selection forces. Identifying the differences in natural selection between populations is critical for understanding the roles of specific genetic variants in evolutionary adaptation. Although numerous methods have been developed to detect genetic loci under recent directional selection, a probabilistic solution for testing and quantifying selection differences between populations is lacking. Here we report the development of a probabilistic method for testing and estimating selection differences between populations. By use of a probabilistic model of genetic drift and selection, we showed that logarithm odds ratios of allele frequencies provide estimates of the differences in selection coefficients between populations. The estimates approximate a normal distribution, and variance can be estimated using genome-wide variants. This allows us to quantify differences in selection coefficients and to determine the confidence intervals of the estimate. Our work also revealed the link between genetic association testing and hypothesis testing of selection differences. It therefore supplies a solution for hypothesis testing of selection differences. This method was applied to a genome-wide data analysis of Han and Tibetan populations. The results confirmed that both the EPAS1 and EGLN1 genes are under statistically different selection in Han and Tibetan populations. We further estimated differences in the selection coefficients for genetic variants involved in melanin formation and determined their confidence intervals between continental population groups. Application of the method to empirical data demonstrated the outstanding capability of this novel approach for testing and quantifying differences in natural selection.

Germline gene polymorphisms predisposing domestic mammals to carcinogenesis

Cancer is a complex disease caused in part by predisposing germline gene polymorphisms. Knowledge of carcinogenesis in companion mammals (dog and cat) and some livestock species (pig and horse) is quite advanced. The prevalence of certain cancers varies by breed in these species, suggesting the presence of predisposing genetic variants in susceptible breeds. This review summarizes the present understanding of germline gene polymorphisms, including BRCA1, BRCA2, MC1R, KIT, NRAS and RAD51, associated with predisposition to melanoma, mammary cancer, osteosarcoma and histiocytic sarcoma in dogs, cats, pigs and horses. The predisposing variants in these species are discussed in the context of human germline gene polymorphisms associated with the same types of cancer.

Melanoma genetics

Approximately 10% of melanoma cases report a relative affected with melanoma, and a positive family history is associated with an increased risk of developing melanoma. Although the majority of genetic alterations associated with melanoma development are somatic, the underlying presence of heritable melanoma risk genes is an important component of disease occurrence. Susceptibility for some families is due to mutation in one of the known high penetrance melanoma predisposition genes: CDKN2A, CDK4, BAP1, POT1, ACD, TERF2IP and TERT. However, despite such mutations being implicated in a combined total of approximately 50% of familial melanoma cases, the underlying genetic basis is unexplained for the remainder of high-density melanoma families. Aside from the possibility of extremely rare mutations in a few additional high penetrance genes yet to be discovered, this suggests a likely polygenic component to susceptibility, and a unique level of personal melanoma risk influenced by multiple low-risk alleles and genetic modifiers. In addition to conferring a risk of cutaneous melanoma, some 'melanoma' predisposition genes have been linked to other cancers, with cancer clustering observed in melanoma families at rates greater than expected by chance. The most extensively documented association is between CDKN2A germ line mutations and pancreatic cancer, and a cancer syndrome including cutaneous melanoma, uveal melanoma and mesothelioma has been proposed for BAP1 germ line mutations. Other medium to high penetrance melanoma predisposition genes have been associated with renal cell carcinoma (MITF, BAP1) and glioma (POT1). These associations between melanoma and other cancers hint at the possibility of common pathways for oncogenesis, and better knowledge of these pathways may improve understanding of the genetic basis underpinning familial melanoma. It is likely that 'melanoma' risk genes will impact on mutation screening and genetic counselling not only for melanoma but also a range of other cancers.

Comorbidity in the Tunisian population

Genetic diseases in the Tunisian population represent a real problem of public health as their spectrum encompasses more than 400 disorders. Their frequency and distribution in the country have been influenced by demographic, economic and social features especially consanguinity. In this article, we report on genetic disease association referred to as comorbidity and discuss factors influencing their expressivity. Seventy-five disease associations have been reported among Tunisian families. This comorbidity could be individual or familial. In 39 comorbid associations, consanguinity was noted. Twenty-one founder and 11 private mutations are the cause of 34 primary diseases and 13 of associated diseases. As the information dealing with this phenomenon is fragmented, we proposed to centralize it in this report in order to draw both clinicians' and researcher's attention on the occurrence of such disease associations in inbred populations as it makes genetic counseling and prenatal diagnosis challenging even when mutations are known.

Inherited genetic variants associated with occurrence of multiple primary melanoma

Recent studies, including genome-wide association studies, have identified several putative low-penetrance susceptibility loci for melanoma. We sought to determine their generalizability to genetic predisposition for multiple primary melanoma in the international population-based Genes, Environment, and Melanoma (GEM) Study. GEM is a case-control study of 1,206 incident cases of multiple primary melanoma and 2,469 incident first primary melanoma participants as the control group. We investigated the odds of developing multiple primary melanoma for 47 SNPs from 21 distinct genetic regions previously reported to be associated with melanoma. ORs and 95% confidence intervals were determined using logistic regression models adjusted for baseline features (age, sex, age by sex interaction, and study center). We investigated univariable models and built multivariable models to assess independent effects of SNPs. Eleven SNPs in 6 gene neighborhoods (TERT/CLPTM1L, TYRP1, MTAP, TYR, NCOA6, and MX2) and a PARP1 haplotype were associated with multiple primary melanoma. In a multivariable model that included only the most statistically significant findings from univariable modeling and adjusted for pigmentary phenotype, back nevi, and baseline features, we found TERT/CLPTM1L rs401681 (P = 0.004), TYRP1 rs2733832 (P = 0.006), MTAP rs1335510 (P = 0.0005), TYR rs10830253 (P = 0.003), and MX2 rs45430 (P = 0.008) to be significantly associated with multiple primary melanoma, while NCOA6 rs4911442 approached significance (P = 0.06). The GEM Study provides additional evidence for the relevance of these genetic regions to melanoma risk and estimates the magnitude of the observed genetic effect on development of subsequent primary melanoma.

Proton-associated sucrose transport of mammalian solute carrier family 45: an analysis in Saccharomyces cerevisiae

The members of the solute carrier 45 (SLC45) family have been implicated in the regulation of glucose homoeostasis in the brain (SLC45A1), with skin and hair pigmentation (SLC45A2), and with prostate cancer and myelination (SLC45A3). However, apart from SLC45A1, a proton-associated glucose transporter, the function of these proteins is still largely unknown, although sequence similarities to plant sucrose transporters mark them as a putative sucrose transporter family. Heterologous expression of the three members SLC45A2, SLC45A3 and SLC45A4 in Saccharomyces cerevisiae confirmed that they are indeed sucrose transporters. [(14)C]Sucrose-uptake measurements revealed intermediate transport affinities with Km values of approximately 5 mM. Transport activities were best under slightly acidic conditions and were inhibited by the protonophore carbonyl cyanide m-chlorophenylhydrazone, demonstrating an H(+)-coupled transport mechanism. Na(+), on the other hand, had no effect on sucrose transport. Competitive inhibition assays indicated a possible transport also of glucose and fructose. Real-time PCR of mouse tissues confirmed mRNA expression of SLC45A2 in eyes and skin and of SLC45A3 primarily in the prostate, but also in other tissues, whereas SLC45A4 showed a predominantly ubiquitous expression. Altogether the results provide new insights into the physiological significance of SLC45 family members and challenge existing concepts of mammalian sugar transport, as they (i) transport a disaccharide, and (ii) perform secondary active transport in a proton-dependent manner.

Tyrosinase gene mutations in the Chinese Han population with OCA1

Oculocutaneous albinism (OCA) is a heterogeneous autosomal recessive genetic disorder that affects melanin synthesis. OCA results in reduced or absent pigmentation in the hair, skin and eyes. Type 1 OCA (OCA1) is the result of tyrosinase (TYR) gene mutations and is a severe disease type. This study investigated TYR mutations in a Chinese cohort with OCA1. This study included two parts: patient genetic study and prenatal genetic diagnosis. A total of 30 OCA1 patients were subjected to TYR gene mutation analysis. Ten pedigrees were included for prenatal genetic diagnosis. A total of 100 unrelated healthy Chinese individuals were genotyped for controls. The coding sequence and the intron/exon junctions of TYR were analysed by bidirectional DNA sequencing. In this study, 20 mutations were identified, four of which were novel. Of these 30 OCA1 patients, 25 patients were TYR compound heterozygous; two patients carried homozygous TYR mutations; and three were heterozygous. Among the ten prenatally genotyped fetuses, three fetuses carried compound heterozygous mutations and seven carried no mutation or only one mutant allele of TYR and appeared normal at birth. In conclusion, we identified four novel TYR mutations and showed that molecular-based prenatal screening to detect TYR mutations in a fetus at risk for OCA1 provided essential information for genetic counselling of couples at risk.

Association of cancer susceptibility variants with risk of multiple primary cancers: The population architecture using genomics and epidemiology study

BACKGROUND

Multiple primary cancers account for approximately 16% of all incident cancers in the United States. Although genome-wide association studies (GWAS) have identified many common genetic variants associated with various cancer sites, no study has examined the association of these genetic variants with risk of multiple primary cancers (MPC).

METHODS

As part of the National Human Genome Research Institute (NHGRI) Population Architecture using Genomics and Epidemiology (PAGE) study, we used data from the Multiethnic Cohort (MEC) and Women's Health Initiative (WHI). Incident MPC (IMPC) cases (n = 1,385) were defined as participants diagnosed with more than one incident cancer after cohort entry. Participants diagnosed with only one incident cancer after cohort entry with follow-up equal to or longer than IMPC cases served as controls (single-index cancer controls; n = 9,626). Fixed-effects meta-analyses of unconditional logistic regression analyses were used to evaluate the associations between 188 cancer risk variants and IMPC risk. To account for multiple comparisons, we used the false-positive report probability (FPRP) to determine statistical significance.

RESULTS

A nicotine dependence-associated and lung cancer variant, CHRNA3 rs578776 [OR, 1.16; 95% confidence interval (CI), 1.05-1.26; P = 0.004], and two breast cancer variants, EMBP1 rs11249433 and TOX3 rs3803662 (OR, 1.16; 95% CI, 1.04-1.28; P = 0.005 and OR, 1.13; 95% CI, 1.03-1.23; P = 0.006), were significantly associated with risk of IMPC. The associations for rs578776 and rs11249433 remained (P < 0.05) after removing subjects who had lung or breast cancers, respectively (P ≤ 0.046). These associations did not show significant heterogeneity by smoking status (Pheterogeneity ≥ 0.53).

CONCLUSIONS

Our study has identified rs578776 and rs11249433 as risk variants for IMPC.

IMPACT

These findings may help to identify genetic regions associated with IMPC risk.

Exploratory investigation of genetic associations with basal cell carcinoma risk: genome-wide association study in Jeju Island, Korea

AIM

Little is known about the genetic associations with Basal cell carcinoma (BCC) risk in non-Caucasian populations, in which BCC is rare, as in Korea. We here conducted a pilot genome-wide association study (GWAS) in 12 patients and 48 standard controls.

METHOD

A total of 263,511 SNPs were analyzed with the Illumina HumanOmni1 Quad v1.0 DNA Analysis BeadChip for cases and Korean HapMap 570K for controls.

RESULTS

SNP-based analyses, based on the allele genetic model with adjustment for sex and age showed suggestive associations with BCC risk for 6 SNPs with a P-value (P < 0.0005). However, these associations were not statistically significant after Bonferroni correction: rs1040503, rs2216491, rs13407683, rs4751072, rs9891263, and rs1368474. In addition, results from gene-based analyses showed suggestive associations with BCC risk for 33 candidate genes with a P-value (P <0.0005). Consistent with previous GWAS and replication studies in Caucasian populations, PADI6, RHOU and SLC45A2 were identified as having null associations with BCC (P > 0.05), likely due to the smaller sample size.

CONCLUSIONS

Although this was a small-scale negative study, to our knowledge, we have conducted the first GWAS for BCC risk in an Asian population. Further large studies in non-Caucasian populations are required to achieve statistical significance and confirm these findings.

The interplay between natural selection and susceptibility to melanoma on allele 374F of SLC45A2 gene in a South European population

We aimed to study the selective pressures interacting on SLC45A2 to investigate the interplay between selection and susceptibility to disease. Thus, we enrolled 500 volunteers from a geographically limited population (Basques from the North of Spain) and by resequencing the whole coding region and intron 5 of the 34 most and the 34 least pigmented individuals according to the reflectance distribution, we observed that the polymorphism Leu374Phe (L374F, rs16891982) was statistically associated with skin color variability within this sample. In particular, allele 374F was significantly more frequent among the individuals with lighter skin. Further genotyping an independent set of 558 individuals of a geographically wider population with known ancestry in the Spanish population also revealed that the frequency of L374F was significantly correlated with the incident UV radiation intensity. Selection tests suggest that allele 374F is being positively selected in South Europeans, thus indicating that depigmentation is an adaptive process. Interestingly, by genotyping 119 melanoma samples, we show that this variant is also associated with an increased susceptibility to melanoma in our populations. The ultimate driving force for this adaptation is unknown, but it is compatible with the vitamin D hypothesis. This shows that molecular evolution analysis can be used as a useful technology to predict phenotypic and biomedical consequences in humans.

Increased risk of developing cutaneous malignant melanoma is associated with variation in pigmentation genes and VDR, and may involve epistatic effects

Cutaneous malignant melanoma (CMM) is a malicious human skin cancer that primarily affects individuals with light pigmentation and heavy sun exposure, but also has a known familial association. Multiple genes and polymorphisms have been reported as low-penetrance susceptibility loci for CMM. Here, we examined 33 candidate polymorphisms located in 11 pigmentation genes and the vitamin D receptor gene (VDR) in a population of 130 cutaneous melanoma patients and 707 healthy controls. The genotypes obtained were evaluated for main association effects and potential gene-gene interactions. MC1R, TYR, VDR and SLC45A2 genes were found to be associated with CMM in our population. The results obtained for major function MC1R mutations were the most significant [with odds ratio (OR)=1.787, confidence interval (CI)=1.320-2.419 and P=1.715(-4)], followed by TYR (rs1393350) (with OR=1.569, CI=1.162-2.118, P=0.003), VDR (GCCC haplotype in rs2238136-rs4516035-rs7139166-rs11568820 block) (with OR=5.653, CI=1.794-17.811, P=0.003) and SLC45A2 (rs16891982) (with OR=0.238, CI=0.057-0.987, P=0.048). The study also detected significant intermolecular epistatic effects between MC1R and TYR, SLC45A2 and VDR, HERC2 and VDR, OCA2 and TPCN2, as well as intramolecular interactions between variants within the genes MC1R and VDR. In the final multivariate logistic regression model for CMM development, only the gene-gene interactions discovered remained significant, showing that epistasis may be an important factor in the risk of melanoma.

Putative role of HIF transcriptional activity in melanocytes and melanoma biology

Hypoxia-inducible factor-1α (HIF-1α) is a highly oxygen sensitive bHLH protein that is part of the heterodimeric HIF-1 transcription factor. Under hypoxic stress, HIF-1 activity is induced to control expression of multiple downstream target genes, including vascular endothelial growth factor (VEGF). The normal epidermis exists in a constant mild hypoxic microenvironment and constitutively expresses HIF-1α and HIF-2α. Expression of HIF-1α and/or HIF-2α has been suggested to correlate with the increased malignant potential of melanocytes, therefore, failures of melanoma therapies may be partially linked to high HIF activity. Notably, melanomas that have the V600E BRAF mutation exhibit increased HIF-1α expression. We have utilized a bioinformatics approach to identify putative hypoxia response elements (HREs) in a set of genes known to participate in the process of melanogenesis (includingTRPM1, SLC45A2, HRAS, C-KIT, PMEL and CRH). While some of the mechanistic links between these genes and the HIF pathway have been previously explored, others await further investigation. Although agents targeting HIF activity have been proposed as novel treatment modalities for melanoma, there are currently no clinical trials in progress to test their efficacy in melanoma.