Discovery of copy number variants by multiplex amplifiable probe hybridization (MAPH) in candidate pigmentation genes.

BACKGROUND: Copy Number Variants (CNVs) contribute to a large fraction of genetic diversity and some of them have been reported to offer an evolutionary advantage. AIM: To identify CNVs in pigmentary loci that could contribute to human skin pigmentation diversity. SUBJECTS AND METHODS: This study assessed the existence of CNVs in every exon of candidate genes: TYR, TYRP1, DCT, MC1R and SLC24A5, using the Multiplex Amplifiable Probe Hybridization technique (MAPH). This study analysed a total of 99 DNA samples of unrelated individuals from different populations. Validation and further analysis in a larger Spanish sample were performed by RT-qPCR. RESULTS: Five CNVs were identified by MAPH: DCT exons 4 and 8, TYR exon 1 and SLC24A5 exons 1 and 4. Real-time quantitative PCR (RT-qPCR) confirmed the CNV in exon 1 of SLC24A5. This study further analysed the 5' promoter region of SLC24A5 and found another CNV in this region. However, no association was found between the CNV and the degree of pigmentation. CONCLUSION: Although the functional role of these structural variants in pigmentation should be the subject of future work, the results emphasize the need to consider all classes of variation (both SNPs and CNVs) when exploring the genetics of skin pigmentation.

Simultaneous purifying selection on the ancestral MC1R allele and positive selection on the melanoma-risk allele V60L in south Europeans.

In humans, the geographical apportionment of the coding diversity of the pigmentary locus melanocortin-1 receptor (MC1R) is, unusually, higher in Eurasians than in Africans. This atypical observation has been interpreted as the result of purifying selection due to functional constraint on MC1R in high UV-B radiation environments. By analyzing 3,142 human MC1R alleles from different regions of Spain in the context of additional haplotypic information from the 1000 Genomes (1000G) Project data, we show that purifying selection is also strong in southern Europe, but not so in northern Europe. Furthermore, we show that purifying and positive selection act simultaneously on MC1R. Thus, at least in Spain, regions at opposite ends of the incident UV-B radiation distribution show significantly different frequencies for the melanoma-risk allele V60L (a mutation also associated to red hair and fair skin and even blonde hair), with higher frequency of V60L at those regions of lower incident UV-B radiation. Besides, using the 1000G south European data, we show that the V60L haplogroup is also characterized by an extended haplotype homozygosity (EHH) pattern indicative of positive selection. We, thus, provide evidence for an adaptive value of human skin depigmentation in Europe and illustrate how an adaptive process can simultaneously help to maintain a disease-risk allele. In addition, our data support the hypothesis proposed by Jablonski and Chaplin (Human skin pigmentation as an adaptation to UVB radiation. Proc Natl Acad Sci U S A. 2010;107:8962-8968), which posits that habitation of middle latitudes involved the evolution of partially depigmented phenotypes that are still capable of suitable tanning.

IL-12/15/18-induced cell death and mitochondrial dynamics of human NK cells.

Natural killer (NK) cells are lymphocytes with potent antitumor functions and, consequently, several NK cell-based strategies have been developed for cancer immunotherapy. A remarkable therapeutic approach is the adoptive transfer of NK cells stimulated with IL-12, IL-15 and IL-18. This cytokine stimulation endows NK cells with properties that resemble immunological memory and, for this reason, they are known as cytokine-induced memory-like (CIML) NK cells. Very promising results have been reported in clinical trials and yet, there are still unknown aspects of CIML NK cells. Here, we have conducted a preliminary study of their mitochondrial dynamics. Our results show that upon IL-12/15/18 stimulation the viability of NK cells decreased and an increment in mitochondrial superoxide levels was observed. In addition, we found that mitochondria appeared slightly elongated and their cristae density decreased following IL-12/15/18 stimulation, possibly in a process mediated by the low levels of optic atrophy type 1 (OPA1) protein. Interestingly, although mitophagy was slightly impaired, an increase in autophagic flux was observed, which might explain the reduced viability and the accumulation of unfit mitochondria. Our findings could be of relevance in order to design new strategies intended to improve the mitochondrial fitness of IL-12/15/18-stimulated NK cells with the aim of improving their therapeutic efficacy.

Melanoma Clinical Decision Support System: An Artificial Intelligence-Based Tool to Diagnose and Predict Disease Outcome in Early-Stage Melanoma Patients.

This study set out to assess the performance of an artificial intelligence (AI) algorithm based on clinical data and dermatoscopic imaging for the early diagnosis of melanoma, and its capacity to define the metastatic progression of melanoma through serological and histopathological biomarkers, enabling dermatologists to make more informed decisions about patient management. Integrated analysis of demographic data, images of the skin lesions, and serum and histopathological markers were analyzed in a group of 196 patients with melanoma. The interleukins (ILs) IL-4, IL-6, IL-10, and IL-17A as well as IFNγ (interferon), GM-CSF (granulocyte and macrophage colony-stimulating factor), TGFß (transforming growth factor), and the protein DCD (dermcidin) were quantified in the serum of melanoma patients at the time of diagnosis, and the expression of the RKIP, PIRIN, BCL2, BCL3, MITF, and ANXA5 proteins was detected by immunohistochemistry (IHC) in melanoma biopsies. An AI algorithm was used to improve the early diagnosis of melanoma and to predict the risk of metastasis and of disease-free survival. Two models were obtained to predict metastasis (including "all patients" or only patients "at early stages of melanoma"), and a series of attributes were seen to predict the progression of metastasis: Breslow thickness, infiltrating BCL-2 expressing lymphocytes, and IL-4 and IL-6 serum levels. Importantly, a decrease in serum GM-CSF seems to be a marker of poor prognosis in patients with early-stage melanomas.

MGRN1 as a Phenotypic Determinant of Human Melanoma Cells and a Potential Biomarker.

Mahogunin Ring Finger 1 (MGRN1), a ubiquitin ligase expressed in melanocytes, interacts with the α melanocyte-stimulating hormone receptor, a well-known melanoma susceptibility gene. Previous studies showed that MGRN1 modulates the phenotype of mouse melanocytes and melanoma cells, with effects on pigmentation, shape, and motility. Moreover, MGRN1 knockdown augmented the burden of DNA breaks in mouse cells, indicating that loss of MGRN1 promoted genomic instability. However, data concerning the roles of MGRN1 in human melanoma cells remain scarce. We analyzed MGRN1 knockdown in human melanoma cells. Transient MGRN1 depletion with siRNA or permanent knockdown in human melanoma cells by CRISPR/Cas9 caused an apparently MITF-independent switch to a more dendritic phenotype. Lack of MGRN1 also increased the fraction of human cells in the S phase of the cell cycle and the burden of DNA breaks but did not significantly impair proliferation. Moreover, in silico analysis of publicly available melanoma datasets and estimation of MGRN1 in a cohort of clinical specimens provided preliminary evidence that MGRN1 expression is higher in human melanomas than in normal skin or nevi and pointed to an inverse correlation of MGRN1 expression in human melanoma with patient survival, thus suggesting potential use of MGRN1 as a melanoma biomarker.

Association of TYR SNP rs1042602 with Melanoma Risk and Prognosis.

Cutaneous melanoma is the most aggressive of skin tumors. In order to discover new biomarkers that could help us improve prognostic prediction in melanoma patients, we have searched for germline DNA variants associated with melanoma progression. Thus, after exome sequencing of a set of melanoma patients and healthy control individuals, we identified rs1042602, an SNP within TYR, as a good candidate. After genotyping rs1042602 in 1025 patients and 773 healthy donors, we found that the rs1042602-A allele was significantly associated with susceptibility to melanoma (CATT test: p = 0.0035). Interestingly, we also observed significant differences between patients with good and bad prognosis (5 years of follow-up) (n = 664) (CATT test for all samples p = 0.0384 and for men alone p = 0.0054). Disease-free-survival (DFS) analyses also showed that patients with the A allele had shorter DFS periods. In men, the association remained significant even in a multivariate Cox Proportional-hazards model, which was adjusted for age at diagnosis, Breslow thickness, ulceration and melanoma subtype (HR 0.4; 95% confidence interval (CI) 0.20-0.83; p = 0.0139). Based on our results, we propose that rs1042602-A is a risk allele for melanoma, which also seems to be responsible for a poorer prognosis of the disease, particularly in men.

BRAF V600E mutational load as a prognosis biomarker in malignant melanoma.

Analyzing the mutational load of driver mutations in melanoma could provide valuable information regarding its progression. We aimed at analyzing the heterogeneity of mutational load of BRAF V600E in biopsies of melanoma patients of different stages, and investigating its potential as a prognosis factor. Mutational load of BRAF V600E was analyzed by digital PCR in 78 biopsies of melanoma patients of different stages and 10 nevi. The BRAF V600E load was compared among biopsies of different stages. Results showed a great variability in the load of V600E (0%-81%). Interestingly, we observed a significant difference in the load of V600E between the early and late melanoma stages, in the sense of an inverse correlation between BRAF V600E mutational load and melanoma progression. In addition, a machine learning approach showed that the mutational load of BRAF V600E could be a good predictor of metastasis in stage II patients. Our results suggest that BRAF V600E is a promising biomarker of prognosis in stage II patients.

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.