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.

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.

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.

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.

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.

Severe coronavirus disease 2019 (COVID-19) often lead to acute respiratory distress syndrome and multi-organ dysfunction, driven by a dysregulated immune response, including a cytokine storm with elevated proinflammatory cytokine levels. Natural killer (NK) cells are part of the innate immune system with a fundamental role in the defense against viral infections. However, during COVID-19 acute infection, they exhibit an altered phenotype and impaired functionality contributing to the immunopathogenesis of the disease. In this work, we have studied a cohort of patients with COVID-19 (ranging from mild to severe) by analyzing IL-15, TGF-β, PlGF and GDF-15 plasma levels and performing multiparametric flow cytometry studies. Our results revealed that severe COVID-19 patients exhibited high levels of IL-15, PlGF and GDF-15, along with an enrichment of a NK cell subset expressing the CD151 tetraspanin, which correlated with IL-15 plasma levels and disease severity. In patients, these CD151+ NK cells displayed a more activated phenotype characterized by an increased expression of HLA-DR, CD38 and granzyme B, a distinct receptor repertoire, with lower levels of CD160 and CD31 and higher levels of CD55 and, remarkably, a higher expression of tissue-resident markers CD103 and the NK cell decidual marker CD9. Last of all, in individuals with severe disease, we identified an expansion of a CD151^brightCD9+ NK cell subset, suggesting that these cells play a specific role in COVID-19. Altogether, our findings suggest that CD151+ NK cells may have a relevant role in COVID-19 immunopathogenesis.