Genome of Pestera Muierii skull shows high diversity and low mutational load in pre-glacial Europe.

Few complete human genomes from the European Early Upper Palaeolithic (EUP) have been sequenced. Using novel sampling and DNA extraction approaches, we sequenced the genome of a woman from "Pestera Muierii," Romania who lived ∼34,000 years ago to 13.5× coverage. The genome shows similarities to modern-day Europeans, but she is not a direct ancestor. Although her cranium exhibits both modern human and Neanderthal features, the genome shows similar levels of Neanderthal admixture (∼3.1%) to most EUP humans but only half compared to the ∼40,000-year-old Pestera Oase 1. All EUP European hunter-gatherers display high genetic diversity, demonstrating that the severe loss of diversity occurred during and after the Last Glacial Maximum (LGM) rather than just during the out-of-Africa migration. The prevalence of genetic diseases is expected to increase with low diversity; however, pathogenic variant load was relatively constant from EUP to modern times, despite post-LGM hunter-gatherers having the lowest diversity ever observed among Europeans.

Contribution of ancient human remains analysis to the understanding of the variability in HLA-B gene variants in relation to the diagnosis of spondyloarthropathies.

Genetic investigations on ancient human remains affected by rheumatological pathologies are a research field of particular interest for identifying the pathogenesis of diseases, especially those having an autoimmune background such as spondyloarthopaties (SpA). Reliable studies concerning this topic require collaboration between multiple disciplines, usually starting from paleopathologic observations up to molecular genetic screening. Here, we focused our investigation in a medieval necropolis in the Basque Country (13th-15th century, N = 163), which presents a high frequency of joint pathologies through two approaches: on the one hand, the analysis of joint manifestations for the differential diagnosis of the SpA and, on the other hand, the determination of the alleles of the HLA-B gene. The morphological analysis allowed determining that 30% of the individuals had rheumatic bone manifestations, with SpA being the most frequent (45%). The genetic analysis of individuals with and without pathologies, based on the study of the HLA-B gene, allowed finding 17 alleles for this gene, with HLA-B40, HLA-B27 and HLA-B35 being the most frequent. Although these alleles have been traditionally described as genetic markers associated to the development of SpA, in this study they were also found in individuals with other rheumatic diseases (osteoarthritis and rheumatoid arthritis) and even in individuals without pathologies. These data confirm the complexity of the relationship of the HLA-B gene variants with SpA, since it is not possible to establish a diagnosis of SpA with these variants alone. However, we suggest that allele HLA-B40, in combination with some specific rheumatic bone manifestations, facilitates the diagnosis of SpA.

A recording form for differential diagnosis of arthropathies.

The present study is focused on a group of arthropathies that may have very similar bone manifestations (rheumatoid arthritis, ankylosing spondylitis, reactive arthritis, psoriatic arthritis, osteoarthritis and diffuse idiopathic skeletal hyperostosis), which makes it more difficult to diagnose them in human remains from archaeological contexts. A stepwise recording form was designed in order to improve the identification and differential diagnosis of these pathological conditions in bone remains, particularly in joint manifestations of the spine, pelvis, hands, feet and other limb joints. This recording form was applied in the analysis of two medieval individuals from the Basque Country (Spain) who presented very severe arthropathic manifestations. The use of this recording form allowed the researchers the diagnosis of ankylosing spondylitis in one of them and diffuse idiopathic skeletal hyperostosis in the other.

Direct U-series analysis of the Lezetxiki humerus reveals a Middle Pleistocene age for human remains in the Basque Country (northern Iberia).

In 1964, a human humerus was found in a sedimentary deposit in Lezetxiki Cave (Basque Country, northern Iberia). The first studies on the stratigraphy, associated mammal faunal remains and lithic implements placed the deposits containing the humerus into the Riss glacial stage. Direct chronometric evidence has so far been missing, and the previous chronostratigraphic framework and faunal dating gave inconsistent results. Here we report laser ablation U-series analyses on the humerus yielding a minimum age of 164 ± 9 ka, corresponding to MIS 6. This is the only direct dating analysis of the Lezetxiki humerus and confirms a Middle Pleistocene age for this hominin fossil. Morphometric analyses suggest that the Lezetxiki humerus has close affinities to other Middle Pleistocene archaic hominins, such as those from La Sima de los Huesos at Atapuerca. This emphasizes the significance of the Lezetxiki fossil within the populations that predate the Neanderthals in south-western Europe. It is thus an important key fossil for the understanding of human evolution in Europe during the Middle Pleistocene, a time period when a great morphological diversity is observed but whose phylogenetic meaning is not yet fully understood.

Demographic differences between funerary caves and megalithic graves of northern Spanish Late Neolithic/Early Chalcolithic.

OBJECTIVES: The study focuses on the estimation of demographic parameters of Late Neolithic/Early Chalcolithic (mid 4th-early 3rd millenniums cal. BC) burial sites from the La Rioja region (Ebro valley, northern Spain) to identify demographic characteristics. MATERIALS AND METHODS: The human remains come from three caves (Las Yurdinas II, Peña Larga, and La Peña de Marañón) and three megalithic graves (Alto de la Huesera, San Martín, and Peña Guerra II). The total skeletal sample consists of a minimum of 261 individuals, 149 being buried in caves and 112 in megalithic graves. Data based on age and sex estimation are analyzed using abridged life tables, mortality rates, and sex ratios. RESULTS: A systematic bias against children under 5 years of age is detected both in caves (5 q0 = 187.92%) and megalithic graves (5 q0 = 71.43%) but also against some juveniles and adults compared with population models, though a statistically significant greater lack of infants is worth noting in the megaliths (t-test, P = 0.012). Moreover, a significant divergence in sex ratios (χ(2) , P = 0.002) is also identified between site types, clearly prioritizing women in caves (sex ratio = 0.45) and men in megalithic graves (sex ratio = 1.33). CONCLUSIONS: This evidence is interpreted as the result of different selective burial patterns. The mortuary variability could lie behind intragroup differential status relationships, though the hypothesis of two populations performing distinct funerary practices in a small region cannot be rejected at the present state of the research. Am J Phys Anthropol 160:284-297, 2016. © 2016 Wiley Periodicals, Inc.

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.