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.

Evolution: On the origin of Basques.

As the only speakers of a non-Indo-European language lacking connections to any other language, the origin and evolution of the Basque population have been topics of interest. A new genomic analysis of Basques shows that they are not a relict Paleolithic population as once thought.

Comparison of the Transcriptional Profiles of Melanocytes from Dark and Light Skinned Individuals under Basal Conditions and Following Ultraviolet-B Irradiation.

We analysed the whole-genome transcriptional profile of 6 cell lines of dark melanocytes (DM) and 6 of light melanocytes (LM) at basal conditions and after ultraviolet-B (UVB) radiation at different time points to investigate the mechanisms by which melanocytes protect human skin from the damaging effects of UVB. Further, we assessed the effect of different keratinocyte-conditioned media (KCM+ and KCM-) on melanocytes. Our results suggest that an interaction between ribosomal proteins and the P53 signaling pathway may occur in response to UVB in both DM and LM. We also observed that DM and LM show differentially expressed genes after irradiation, in particular at the first 6h after UVB. These are mainly associated with inflammatory reactions, cell survival or melanoma. Furthermore, the culture with KCM+ compared with KCM- had a noticeable effect on LM. This effect includes the activation of various signaling pathways such as the mTOR pathway, involved in the regulation of cell metabolism, growth, proliferation and survival. Finally, the comparison of the transcriptional profiles between LM and DM under basal conditions, and the application of natural selection tests in human populations allowed us to support the significant evolutionary role of MIF and ATP6V0B in the pigmentary phenotype.

Ancient DNA from South-East Europe Reveals Different Events during Early and Middle Neolithic Influencing the European Genetic Heritage.

The importance of the process of Neolithization for the genetic make-up of European populations has been hotly debated, with shifting hypotheses from a demic diffusion (DD) to a cultural diffusion (CD) model. In this regard, ancient DNA data from the Balkan Peninsula, which is an important source of information to assess the process of Neolithization in Europe, is however missing. In the present study we show genetic information on ancient populations of the South-East of Europe. We assessed mtDNA from ten sites from the current territory of Romania, spanning a time-period from the Early Neolithic to the Late Bronze Age. mtDNA data from Early Neolithic farmers of the Starcevo Cris culture in Romania (Cârcea, Gura Baciului and Negrilesti sites), confirm their genetic relationship with those of the LBK culture (Linienbandkeramik Kultur) in Central Europe, and they show little genetic continuity with modern European populations. On the other hand, populations of the Middle-Late Neolithic (Boian, Zau and Gumelnita cultures), supposedly a second wave of Neolithic migration from Anatolia, had a much stronger effect on the genetic heritage of the European populations. In contrast, we find a smaller contribution of Late Bronze Age migrations to the genetic composition of Europeans. Based on these findings, we propose that permeation of mtDNA lineages from a second wave of Middle-Late Neolithic migration from North-West Anatolia into the Balkan Peninsula and Central Europe represent an important contribution to the genetic shift between Early and Late Neolithic populations in Europe, and consequently to the genetic make-up of modern European populations.

Nondestructive methods for recovery of biological material from human teeth for DNA extraction.

The extraction of DNA from human skeletal remains applied to forensic, and evolutionary studies do not exclude risks, which are to be evaluated when working with unique specimens that could be damaged or even destroyed. In the present study were evaluated several nondestructive methods for recovering DNA instead of the most currently used pulverization method. Three different procedures to access inside the dental pieces (occlusal perforation, cervical perforation, and cervical cut) have been compared with the aim of recovering as many cell remains as possible to carry out a DNA extraction. Given the DNA quantitation results, a method was proposed that consists of a cervical cut to facilitate the access to the pulp cavity and a subsequent filing of the root canals down to the apex of the dental root. This methodology allows the recovery of both mitochondrial and nuclear DNA, with the minimum deterioration for the dental pieces.

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.

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.

The Carpathian range represents a weak genetic barrier in South-East Europe.

BACKGROUND: In the present study we have assessed whether the Carpathian Mountains represent a genetic barrier in East Europe. Therefore, we have analyzed the mtDNA of 128 native individuals of Romania: 62 of them from the North of Romania, and 66 from South Romania. RESULTS: We have analyzed their mtDNA variability in the context of other European and Near Eastern populations through multivariate analyses. The results show that regarding the mtDNA haplogroup and haplotype distributions the Romanian groups living outside the Carpathian range (South Romania) displayed some degree of genetic differentiation compared to those living within the Carpahian range (North Romania). CONCLUSION: The main differentiation between the mtDNA variability of the groups from North and South Romania can be attributed to the demographic movements from East to West (prehistoric or historic) that differently affected in these regions, suggesting that the Carpathian mountain range represents a weak genetic barrier in South-East Europe.

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.

The loss of functional caspase-12 in Europe is a pre-neolithic event.

BACKGROUND: Caspase-12 (CASP12) modulates the susceptibility to sepsis. In humans, the "C" allele at CASP12 rs497116 has been associated with an increased risk of sepsis. Instead, the derived "T" allele encodes for an inactive caspase-12. Interestingly, Eurasians are practically fixed for the inactive variant, whereas in Sub-Saharan Africa the active variant is still common (~24%). This marked structure has been explained as a function of the selective advantage that the inactive caspase-12 confers by increasing resistance to infection. As regards to both when positive selection started acting and as to the speed with which fixation was achieved in Eurasia, estimates depend on the method and assumptions used, and can vary substantially. Using experimental evidence, we propose that, least in Eurasia, the increase in the frequency of the T allele might be related to the selective pressure exerted by the increase in zoonotic diseases transmission caused by the interplay between increased human population densities and a closer contact with animals during the Neolithic. METHODOLOG/PRINCIPAL FINDINGS: We genotyped CASP12 rs497116 in prehistoric individuals from 6 archaeological sites from the North of the Iberian Peninsula that date from Late Upper Paleolithic to Late Neolithic. DNA extraction was done from teeth lacking cavities or breakages using standard anti-contamination procedures, including processing of the samples in a positive pressure, ancient DNA-only chamber, quantitation of DNAs by qPCR, duplication, replication, genotyping of associated animals, or cloning of PCR products. Out of 50, 24 prehistoric individuals could finally be genotyped for rs497116. Only the inactive form of CASP12 was found. CONCLUSIONS/SIGNIFICANCE: We demonstrate that the loss of caspase-12 in Europe predates animal domestication and that consequently CASP12 loss is unlikely to be related to the impact of zoonotic infections transmitted by livestock.

Ancient DNA from hunter-gatherer and farmer groups from Northern Spain supports a random dispersion model for the Neolithic expansion into Europe.

BACKGROUND/PRINCIPAL FINDINGS: The phenomenon of Neolithisation refers to the transition of prehistoric populations from a hunter-gatherer to an agro-pastoralist lifestyle. Traditionally, the spread of an agro-pastoralist economy into Europe has been framed within a dichotomy based either on an acculturation phenomenon or on a demic diffusion. However, the nature and speed of this transition is a matter of continuing scientific debate in archaeology, anthropology, and human population genetics. In the present study, we have analyzed the mitochondrial DNA diversity in hunter-gatherers and first farmers from Northern Spain, in relation to the debate surrounding the phenomenon of Neolithisation in Europe. METHODOLOGY/SIGNIFICANCE: Analysis of mitochondrial DNA was carried out on 54 individuals from Upper Paleolithic and Early Neolithic, which were recovered from nine archaeological sites from Northern Spain (Basque Country, Navarre and Cantabria). In addition, to take all necessary precautions to avoid contamination, different authentication criteria were applied in this study, including: DNA quantification, cloning, duplication (51% of the samples) and replication of the results (43% of the samples) by two independent laboratories. Statistical and multivariate analyses of the mitochondrial variability suggest that the genetic influence of Neolithisation did not spread uniformly throughout Europe, producing heterogeneous genetic consequences in different geographical regions, rejecting the traditional models that explain the Neolithisation in Europe. CONCLUSION: The differences detected in the mitochondrial DNA lineages of Neolithic groups studied so far (including these ones of this study) suggest different genetic impact of Neolithic in Central Europe, Mediterranean Europe and the Cantabrian fringe. The genetic data obtained in this study provide support for a random dispersion model for Neolithic farmers. This random dispersion had a different impact on the various geographic regions, and thus contradicts the more simplistic total acculturation and replacement models proposed so far to explain Neolithisation.

Low prevalence of lactase persistence in Neolithic South-West Europe.

The ability of humans to digest the milk component lactose after weaning requires persistent production of the lactose-converting enzyme lactase. Genetic variation in the promoter of the lactase gene (LCT) is known to be associated with lactase production and is therefore a genetic determinant for either lactase deficiency or lactase persistence during adulthood. Large differences in this genetic trait exist between populations in Africa and the Middle-East on the one hand, and European populations on the other; this is thought to be due to evolutionary pressures exerted by consumption of dairy products in Neolithic populations in Europe. In this study, we have investigated lactase persistence of 26 out of 46 individuals from Late Neolithic through analysis of ancient South-West European DNA samples, obtained from two burials in the Basque Country originating from 5000 to 4500 YBP. This investigation revealed that these populations had an average frequency of lactase persistence of 27%, much lower than in the modern Basque population, which is compatible with the concept that Neolithic and post-Neolithic evolutionary pressures by cattle domestication and consumption of dairy products led to high lactase persistence in Southern European populations. Given the heterogeneity in the frequency of the lactase persistence allele in ancient Europe, we suggest that in Southern Europe the selective advantage of lactose assimilation in adulthood most likely took place from standing population variation, after cattle domestication, at a post-Neolithic time when fresh milk consumption was already fully adopted as a consequence of a cultural influence.

The evolutionary history of TLR4 polymorphisms in Europe.

Infections exert important evolutionary pressures shaping the human genome, especially on genes involved in host defense. A crucial step for host defense is recognition of pathogens by pattern recognition receptors on innate immune cells, among which Toll-like receptor 4 (TLR4) is one of the best known. Genetic variation in TLR4 (Asp299Gly, Thr399Ile) has been recently described. Haplotype frequencies of these polymorphisms differ among African, Asian and European populations, suggesting evolutionary pressures exerted by local infections. The TLR4 299Gly/399Ile haplotype, characteristic mainly of European populations, has relatively high frequency in the Iberian peninsula. This region is also described as refuge area during the last glacial maximum 20,000 years ago, from which repopulation of Europe took place. We speculate that a genetic bottleneck in the Iberian peninsula could have promoted the increased frequency of this haplotype by genetic drift. This hypothesis is supported by three arguments: (1) the West-East gradient of prevalence in the haplotype among European populations; (2) ancient DNA from Neolithic burials in the Iberian peninsula, dated 6,600-4,500 years before present, confirmed the relatively high frequency of this haplotype in the region, and (3) no functional differences between this haplotype and wild-type TLR4 have been found. In contrast, the disappearance of the 299Gly/399Thr haplotype in Europe is most likely due to negative selection due to sepsis. In conclusion, differences in distribution of TLR4 polymorphisms Asp299Gly and Thr399Ile in European populations are most likely due to a combination of population migration events combined with selection due to sepsis.

A customized pigmentation SNP array identifies a novel SNP associated with melanoma predisposition in the SLC45A2 gene.

As the incidence of Malignant Melanoma (MM) reflects an interaction between skin colour and UV exposure, variations in genes implicated in pigmentation and tanning response to UV may be associated with susceptibility to MM. In this study, 363 SNPs in 65 gene regions belonging to the pigmentation pathway have been successfully genotyped using a SNP array. Five hundred and ninety MM cases and 507 controls were analyzed in a discovery phase I. Ten candidate SNPs based on a p-value threshold of 0.01 were identified. Two of them, rs35414 (SLC45A2) and rs2069398 (SILV/CKD2), were statistically significant after conservative Bonferroni correction. The best six SNPs were further tested in an independent Spanish series (624 MM cases and 789 controls). A novel SNP located on the SLC45A2 gene (rs35414) was found to be significantly associated with melanoma in both phase I and phase II (P<0.0001). None of the other five SNPs were replicated in this second phase of the study. However, three SNPs in TYR, SILV/CDK2 and ADAMTS20 genes (rs17793678, rs2069398 and rs1510521 respectively) had an overall p-value<0.05 when considering the whole DNA collection (1214 MM cases and 1296 controls). Both the SLC45A2 and the SILV/CDK2 variants behave as protective alleles, while the TYR and ADAMTS20 variants seem to function as risk alleles. Cumulative effects were detected when these four variants were considered together. Furthermore, individuals carrying two or more mutations in MC1R, a well-known low penetrance melanoma-predisposing gene, had a decreased MM risk if concurrently bearing the SLC45A2 protective variant. To our knowledge, this is the largest study on Spanish sporadic MM cases to date.

The diversity profile of TP53 is influenced by positive selection on the immediately upstream locus WDR79.

BACKGROUND/AIM: TP53 is an efficient central node in a signal transduction network that responds to minimize cancer. However, over 50% of tumors show some mutation in TP53. Thus, one might argue that this single central node network lacks robustness. Therefore, we wanted to investigate if natural selection has played a role in shaping the genomic region containing TP53. METHODS: We have analyzed the HapMap data for evidence of selection using F(ST) pairwise comparisons and the extended haplotype homozygosity test on a 200-kb region encompassing TP53. We have also resequenced 4 kb upstream TP53 in Europeans (including melanoma patients), Asians, Australian Aborigines and Africans. RESULTS: Genetic hitchhiking by a linked, positively selected allele at the nearby gene WDR79 may be partly responsible for the sequence diversity profile of TP53. It can help explain why the TP53 Arg72 allele is the major allele in Europeans even when the alternative allele, 72Pro, has been reported to offer an increased longevity after disease. CONCLUSIONS: Despite the important role of TP53, a complex interplay with other evolutionary forces, which are extrinsic to TP53 function, may have driven the genetic diversity pattern of this locus, and, as a consequence, its structure and function.

Functional and genetic evidence that the Mal/TIRAP allele variant 180L has been selected by providing protection against septic shock.

Adequate responses by our innate immune system toward invading pathogens were of vital importance for surviving infections, especially before the antibiotic era. Recently, a polymorphism in Mal (Ser180Leu, TIRAP rs8177374), an important adaptor protein downstream of the Toll-like receptor (TLR) 2 and 4 pathways, has been described to provide protection against a broad range of infectious pathogens. We assessed the functional effects of this polymorphism in human experimental endotoxemia, and we demonstrate that individuals bearing the TIRAP 180L allele display an increased, innate immune response to TLR4 and TLR2 ligands, but not to TLR9 stimulation. This phenotype has been related to an increased resistance to infection. However, an overshoot in the release of proinflammatory cytokines by TIRAP 180L homozygous individuals suggests a scenario of balanced evolution. We have also investigated the worldwide distribution of the Ser180Leu polymorphism in 14 populations around the globe to correlate the genetic makeup of TIRAP with the local infectious pressures. Based on the immunological, clinical, and genetic data, we propose that this mutation might have been selected in West Eurasia during the early settlement of this region after the out-of-Africa migration of modern Homo sapiens. This combination of functional and genetic data provides unique insights to our understanding of the pathogenesis of sepsis.

Planning the human variome project: the Spain report.

The remarkable progress in characterizing the human genome sequence, exemplified by the Human Genome Project and the HapMap Consortium, has led to the perception that knowledge and the tools (e.g., microarrays) are sufficient for many if not most biomedical research efforts. A large amount of data from diverse studies proves this perception inaccurate at best, and at worst, an impediment for further efforts to characterize the variation in the human genome. Because variation in genotype and environment are the fundamental basis to understand phenotypic variability and heritability at the population level, identifying the range of human genetic variation is crucial to the development of personalized nutrition and medicine. The Human Variome Project (HVP; http://www.humanvariomeproject.org/) was proposed initially to systematically collect mutations that cause human disease and create a cyber infrastructure to link locus specific databases (LSDB). We report here the discussions and recommendations from the 2008 HVP planning meeting held in San Feliu de Guixols, Spain, in May 2008.

Complex signatures of selection for the melanogenic loci TYR, TYRP1 and DCT in humans.

BACKGROUND: The observed correlation between ultraviolet light incidence and skin color, together with the geographical apportionment of skin reflectance among human populations, suggests an adaptive value for the pigmentation of the human skin. We have used Affymetrix U133a v2.0 gene expression microarrays to investigate the expression profiles of a total of 9 melanocyte cell lines (5 from lightly pigmented donors and 4 from darkly pigmented donors) plus their respective unirradiated controls. In order to reveal signatures of selection in loci with a bearing on skin pigmentation in humans, we have resequenced between 4 to 5 kb of the proximal regulatory regions of three of the most differently expressed genes, in the expectation that variation at regulatory regions might account for intraespecific morphological diversity, as suggested elsewhere. RESULTS: Contrary to our expectations, expression profiles did not cluster the cells into unirradiated versus irradiated melanocytes, or into lightly pigmented versus darkly pigmented melanocytes. Instead, expression profiles correlated with the presence of Bovine Pituitary Extract (known to contain alpha-MSH) in the media. This allowed us to differentiate between melanocytes that are synthesizing melanin and those that are not. TYR, TYRP1 and DCT were among the five most differently expressed genes between these two groups. Population genetic analyses of sequence haplotypes of the proximal regulatory flanking-regions included Tajima's D, HEW and DHEW neutrality tests analysis. These were complemented with EHH tests (among others) in which the significance was obtained by a novel approach using extensive simulations under the coalescent model with recombination. We observe strong evidence for positive selection for TYRP1 alleles in Africans and for DCT and TYRP1 in Asians. However, the overall picture reflects a complex pattern of selection, which might include overdominance for DCT in Europeans. CONCLUSION: Diversity patterns clearly evidence adaptive selection in pigmentation genes in Africans and Asians. In Europeans, the evidence is more complex, and both directional and balancing selection may be involved in light skin. As a result, different non-African populations may have acquired light skin by alternative ways, and so light skin, and perhaps dark skin too, may be the result of convergent evolution.

Overdominance in the human genome and olfactory receptor activity.

We investigate the contribution of overdominance to the maintenance of polymorphism in the human genome during the recent evolution of our species. Using the HapMap genotypic information, we have detected that the Gene Ontology term "olfactory receptor activity" is a molecular function overrepresented in genes that have SNPs (Single Nucleotide Polymorphisms) showing higher than expected number of heterozygotes in the HapMap populations. Our results suggest that the diversity of a subset of human olfactory receptors (ORs) may have been maintained by balancing selection, in the form of overdominance. This observation may suggest that the loss of OR genes during the evolution of the human lineage may have been accompanied by an increased capability to discriminate odorants with closely similar structures.

TLR4 polymorphisms, infectious diseases, and evolutionary pressure during migration of modern humans.

Infectious diseases exert a constant evolutionary pressure on the genetic makeup of our innate immune system. Polymorphisms in Toll-like receptor 4 (TLR4) have been related to susceptibility to Gram-negative infections and septic shock. Here we show that two polymorphisms of TLR4, Asp299Gly and Thr399Ile, have unique distributions in populations from Africa, Asia, and Europe. Genetic and functional studies are compatible with a model in which the nonsynonymous polymorphism Asp299Gly has evolved as a protective allele against malaria, explaining its high prevalence in subSaharan Africa. However, the same allele could have been disadvantageous after migration of modern humans into Eurasia, putatively because of increased susceptibility to severe bacterial infections. In contrast, the Asp299Gly allele, when present in cosegregation with Thr399Ile to form the Asp299Gly/Thr399Ile haplotype, shows selective neutrality. Polymorphisms in TLR4 exemplify how the interaction between our innate immune system and the infectious pressures in particular environments may have shaped the genetic variations and function of our immune system during the out-of-Africa migration of modern humans.