General method to unravel ancient population structures through surnames, final validation on Italian data.

We analyze the geographic location of 77,451 different Italian surnames (17,579,891 individuals) obtained from the lists of telephone subscribers of the year 1993. By using a specific neural network analysis (Self-Organizing Maps, SOMs), we automatically identify the geographic origin of 49,117 different surnames. To validate the methodology, we compare the results to a study, previously conducted, on the same database, with accurate supervised methods. By comparing the results, we find an overlap of 97%, meaning that the SOMs methodology is highly reliable and well traces back the geographic origin of surnames at the time of their introduction (Late Middle Ages/Renaissance in Italy). SOMs results enables one to distinguish monophyletic surnames from polyphyletic ones, that is surnames having had a single geographic and historic origin from those that started to be in use, with an identical spelling, in different locations (respectively, 76.06% and 21.05% of the total). As we are interested in geographic origins, polyphyletic surnames are excluded from further analyses. By comparing the present location of each monophyletic surname to its inferred geographic origin in late Middle Ages/Renaissance, we measure the extent of the migrations having occurred in Italy since that time. We find that the percentage of individuals presently living in the very area where their surname started to be in use centuries ago is extremely variable (ranging from 22.77% to 77.86% according to the province), thus meaning that self-assessed regional identities seldom correspond to the "autochthony" they imply. For example the upper part of the Thyrennian coast (Northern Latium, Tuscany) has a strong identity but few "autochthonous" inhabitants (∼28%) having been a passageway from the North to the South of Italy.

The family name as socio-cultural feature and genetic metaphor: from concepts to methods.

A recent workshop entitled "The Family Name as Socio-Cultural Feature and Genetic Metaphor: From Concepts to Methods" was held in Paris in December 2010, sponsored by the French National Centre for Scientific Research (CNRS) and by the journal Human Biology. This workshop was intended to foster a debate on questions related to the family names and to compare different multidisciplinary approaches involving geneticists, historians, geographers, sociologists and social anthropologists. This collective paper presents a collection of selected communications.

Gestation length, mode of delivery, and neonatal line-thickness variation.

The transition from an intra- to extra-uterine environment leaves its mark in deciduous teeth (and first permanent molars) as an accentuated enamel incremental ring called the neonatal line (NL). This prominent microfeature separates the enamel formed during intrauterine life from that formed after leaving the womb. However, while the physical structure of this scar is well known, the bases of its formation are still a matter of investigation. In particular, besides the influence of the birth-related abrupt environmental and dietary changes and the role played by physiological factors such as hypocalcaemia, a direct relationship between NL thickness variation and the physical was trauma implied by the birth dynamics, the Caesarean, and the operative modes are apparently associated with the thinnest and the thickest lines, respectively. By using the histological record from a deciduous dental sample (exfoliated crowns) of 100 modern healthy school-aged children (47 males and 53 females) of reported birth histories (normal delivery mode: 55 cases; Caesarean: 40; operative: 5), we investigated the relationships between birth dynamics and NL thickness variation. The Tukey Honest Significant Difference method was used to test the differences between the means of the grouping levels. The results of our histo-morphological investigation do not support the suggestion that Caesarean-born children display, on average, a thinner enamel scar compared to children associated to a normal delivery mode. Rather, our study points to the influence exerted by factors intimately related to gestational length variation on the degree of expression of the line.

Autochthony and HLA frequencies in the Netherlands: when surnames are useless markers.

To study the genetic variability of the HLA loci A, B, DR, and DQ in the Netherlands, we analyzed more than 13,000 typings provided by the Dutch National Reference Laboratory for Histocompatibility. To investigate any possibly existing population structure, we subdivided the typings by the geographic location of residency of donors and by the historical belonging of their surnames to given provinces. Concerning possible geographic patterns, we found no significant differences between the four provinces examined (North Holland, South Holland, Utrecht, Zeeland). To assess whether such a negative result was related to recent immigration to the area (the richest of the country) that erased possible preexisting patterns of HLA diversity, we reprocessed the database according to the surnames of HLA donors. We obtained two groups: (1) those having a surname typical of the four provinces they inhabit and (2) those with surnames coming from elsewhere. Such an analysis was made possible because of the availability of a database concerning the geographic origin of most Dutch surnames. Even with this surname-based approach, no major differences were found. We conclude that either the western part of the Netherlands was genetically homogeneous before the official introduction of Dutch surnames two centuries ago by Napoleon or surnames have no power in dissecting HLA variability; that is, such variability is the result of recombination phenomena that surnames cannot mirror because they are transmitted virtually unchanged generation after generation. A comparable study by other investigators recommended the use of family names to identify rare HLA haplotypes in France, but now, concerning the Netherlands, we find opposite results. We suggest that a few typing centers may be sufficient to type bone marrow donors, because HLA genetic differences between the different provinces of the Netherlands are extremely low. To maximize the number of donors, such centers should be located in areas providing the easiest access to the largest population of possible donors, thus disregarding the search for a local variability that we did not find.

From surnames to linguistic and genetic diversity: five centuries of internal migrations in Spain.

In a previous study concerning 33,753 single Spanish surnames (considered as tokens) occurring 51,419,788 times we have shown that the present-day geography of contemporary surname variability in Spain still corresponds to the political geography of the country at the end of the Middle Ages. Here we reprocess the same database, by clustering surnames with Self-Organizing Maps (SOMs) according to their geographic distribution, to identify the monophyletic surnames showing a geo-historical origin in one of the 47 provinces of continental Spain. They are 25,714, and they occur 12,348,109 times, meaning that about 75% of the Spanish population bears a surname that had a polyphyletic origin. From monophyletic surnames we compute migration matrices accounting for the internal migrations that took place since five centuries ago, when Spanish surnames started to be patrilineally inherited. The mono/ polyphyletic classification we obtain fits ancient census data and is compatible with published molecular diversity of the Y-chromosomes associated to selected Spanish surnames. Monophyletic surnames indicate that i) the provinces exhibiting a higher percentage of autochthonous surnames are also ii) those from which emigration corresponds to a local isolation-by-distance model of diffusion and iii) those that attracted a lower number of immigrants. These are also the provinces where languages other than Castilian are spoken. We suggest that demographic stability explains linguistic resilience, as people prefer to move to areas in which the linguistic variety is more similar to their own. So far the reciprocal influence of migration and language has been investigated at local scales, here we outline how to investigate it at national scales and for time-depths of centuries.

Footprints of middle ages kingdoms are still visible in the contemporary surname structure of Spain.

To assess whether the present-day geographical variability of Spanish surnames mirrors historical phenomena occurred at the times of their introduction (13th-16th century), and to infer the possible effect of foreign immigration (about 11% of present-day) on the observed patterns of diversity, we have analyzed the frequency distribution of 33,753 unique surnames (tokens) occurring 51,419,788 times, according to the list of Spanish residents of the year 2008. Isonymy measures and surname distances have been computed for, and between, the 47 mainland Spanish provinces and compared to a numerical classification of corresponding language varieties spoken in Spain. The comparison of the two bootstrap consensus trees, representing surname and linguistic variability, suggests a similar picture; major clusters are located in the east (Aragón, Cataluña, Valencia), and in the north of the country (Asturias, Galicia, León). Remaining regions appear to be considerably homogeneous. We interpret this pattern as the long-lasting effect of the surname and linguistic normalization actively led by the Christian kingdoms of the north (Reigns of Castilla y León and Aragón) during and after the southwards reconquest (Reconquista) of the territories ruled by the Arabs from the 8th century to the late 15th century, that is when surnames became transmitted in a fixed way and when Castilian linguistic varieties became increasingly prestigious and spread out. The geography of contemporary surname and linguistic variability in Spain corresponds to the political geography at the end of the Middle-Ages. The synchronicity between surname adoption and the political and cultural effects of the Reconquista have permanently forged a Spanish identity that subsequent migrations, internal or external, did not deface.

Molecular anthropology in the genomic era.

Molecular Anthropology is a relatively young field of research. In fact, less than 50 years have passed since the symposium "Classification and Human Evolution" (1962, Burg Wartenstein, Austria), where the term was formally introduced by Emil Zuckerkandl. In this time, Molecular Anthropology has developed both methodologically and theoretically and extended its applications, so covering key aspects of human evolution such as the reconstruction of the history of human populations and peopling processes, the characterization of DNA in extinct humans and the role of adaptive processes in shaping the genetic diversity of our species. In the current scientific panorama, molecular anthropologists have to face a double challenge. As members of the anthropological community, we are strongly committed to the integration of biological findings and other lines of evidence (e.g. linguistic and archaeological), while keeping in line with methodological innovations which are moving the approach from the genetic to the genomic level. In this framework, the meeting "DNA Polymorphisms in Human Populations: Molecular Anthropology in the Genomic Era" (Rome, December 3-5, 2009) offered an opportunity for discussion among scholars from different disciplines, while paying attention to the impact of recent methodological innovations. Here we present an overview of the meeting and discuss perspectives and prospects of Molecular Anthropology in the genomic era.

Do surname differences mirror dialect variation?

Our focus in this paper is the analysis of surnames, which have been proven to be reliable genetic markers because in patrilineal systems they are transmitted along generations virtually unchanged, similarly to a genetic locus on the Y chromosome. We compare the distribution of surnames to the distribution of dialect pronunciations, which are clearly culturally transmitted. Because surnames, at the time of their introduction, were words subject to the same linguistic processes that otherwise result in dialect differences, one might expect their geographic distribution to be correlated with dialect pronunciation differences. In this paper we concentrate on the Netherlands, an area of only 40,000 km2, where two official languages are spoken, Dutch and Frisian. We analyze 19,910 different surnames, sampled in 226 locations, and 125 different words, whose pronunciation was recorded in 252 sites. We find that, once the collinear effects of geography on both surname and cultural transmission are taken into account, there is no statistically significant association between the two, suggesting that surnames cannot be taken as a proxy for dialect variation, even though they can be safely used as a proxy for Y-chromosome genetic variation. We find the results historically and geographically insightful, hopefully leading to a deeper understanding of the role that local migrations and cultural diffusion play in surname and dialect diversity.

Evidence for new Neanderthal teeth in Tabun Cave (Israel) by the application of self-organizing maps (SOMs).

Morphological and metrical study suggested that seven human teeth from Tabun Cave, Israel were part of the upper dentition of a single, probably Neanderthal, individual renumbered as Tabun BC7. An enamel fragment gave ESR age estimates of 82+/-14 ka (early U-uptake) and 92+/-18 ka (linear uptake) and an age estimate of 90(+30)(-16) ka using U-series disequilibrium. Although metrical analyses suggested Neanderthal affinities, definitive assessment was difficult as the values often fell into the ranges of both Neanderthal and Levantine early modern human samples. Therefore, two further classification analyses were conducted (neural networks using self-organizing maps and homogeneity analysis). Both identify Tabun BC7 as a Neanderthal. Neural networks are a promising tool for paleoanthropological studies as they can provide reliable classifications even with incomplete data.

New method for surname studies of ancient patrilineal population structures, and possible application to improvement of Y-chromosome sampling.

Several studies showed that surnames are good markers to infer patrilineal genetic structures of populations, both on regional and microregional scales. As a case study, the spatial patterns of the 9,929 most common surnames of the Netherlands were analyzed by a clustering method called self-organizing maps (SOMs). The resulting clusters grouped surnames with a similar geographic distribution and origin. The analysis was shown to be in agreement with already known features of Dutch surnames, such as 1) the geographic distribution of some well-known locative suffixes, 2) historical census data, 3) the distribution of foreign surnames, and 4) polyphyletic surnames. Thus, these results validate the SOM clustering of surnames, and allow for the generalization of the technique. This method can be applied as a new strategy for a better Y-chromosome sampling design in retrospective population genetics studies, since the idenfication of surnames with a defined geographic origin enables the selection of the living descendants of those families settled, centuries ago, in a given area. In other words, it becomes possible to virtually sample the population as it was when surnames started to be in use. We show that, in a given location, the descendants of those individuals who inhabited the area at the time of origin of surnames can be as low as approximately 20%. This finding suggests 1) the major role played by recent migrations that are likely to have distorted or even defaced ancient genetic patterns, and 2) that standard-designed samplings can hardly portray a reliable picture of the ancient Y-chromosome variability of European populations.

Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by using Monmonier's algorithm.

When sampling locations are known, the association between genetic and geographic distances can be tested by spatial autocorrelation or regression methods. These tests give some clues to the possible shape of the genetic landscape. Nevertheless, correlation analyses fail when attempting to identify where genetic barriers exist, namely, the areas where a given variable shows an abrupt rate of change. To this end, a computational geometry approach is more suitable because it provides the locations and the directions of barriers and because it can show where geographic patterns of two or more variables are similar. In this frame we have implemented Monmonier's (1973) maximum difference algorithm in a new software package to identify genetic barriers. To provide a more realistic representation of the barriers in a genetic landscape, we implemented in the software a significance test by means of bootstrap matrices analysis. As a result, the noise associated with genetic markers can be visualized on a geographic map and the areas where genetic barriers are more robust can be identified. Moreover, this multiple matrices approach can visualize the patterns of variation associated with different markers in the same overall picture. This improved Monmonier's method is highly reliable and can be applied to nongenetic data whenever sampling locations and a distance matrix between corresponding data are available.

Detecting recombination in TT virus: a phylogenetic approach.

TT virus (TTV) has a remarkable genetic heterogeneity. To study TTV evolution, phylogenetic analyses were performed on 739 DNA sequences mapping in the N22 region of ORF1. Analysis of neighbor-joining consensus trees shows significant differences between DNA and protein phylogeny. Median joining networks phylogenetic clustering indicates that DNA sequence analysis is biased by homoplasy (i.e., genetic variability not originated by descent), indicative of either hypermutation or recombination. Statistical analysis shows that the significant excess of homoplasy is due to frequent recombination among closely related strains. Recombination events imply that the transmission of TTV is not clonal and provide the necessary basis to explain (i) the high degree of genetic divergence between TTV isolates, (ii) the lack of population structure on a world scale, and (iii) the number of highly divergent strains that seems typical of this virus. We show that recombination phenomena can be detected by phylogenetic analyses in very short sequences when a sufficiently large data set is available.

Y-chromosome analysis in Egypt suggests a genetic regional continuity in Northeastern Africa.

The geographic location of Egypt, at the interface between North Africa, the Middle East, and southern Europe, prompted us to investigate the genetic diversity of this population and its relationship with neighboring populations. To assess the extent to which the modern Egyptian population reflects this intermediate geographic position, ten Unique Event Polymorphisms (UEPs), mapping to the nonrecombining portion of the Y chromosome, have been typed in 164 Y chromosomes from three North African populations. The analysis of these binary markers, which define 11 Y-chromosome lineages, were used to determine the haplogroup frequencies in Egyptians, Moroccan Arabs, and Moroccan Berbers and thereby define the Y-chromosome background in these regions. Pairwise comparisons with a set of 15 different populations from neighboring European, North African, and Middle Eastern populations and geographic analysis showed the absence of any significant genetic barrier in the eastern part of the Mediterranean area, suggesting that genetic variation and gene flow in this area follow the "isolation-by-distance" model. These results are in sharp contrast with the observation of a strong north-south genetic barrier in the western Mediterranean basin, defined by the Gibraltar Strait. Thus, the Y-chromosome gene pool in the modern Egyptian population reflects a mixture of European, Middle Eastern, and African characteristics, highlighting the importance of ancient and recent migration waves, followed by gene flow, in the region.