Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

Founder lineages in the Iberian Roma mitogenomes recapitulate the Roma diaspora and show the effects of demographic bottlenecks

The Roma are the largest ethnic minority in Europe. With a Northwestern Indian origin around ~ 1.5 kya, they travelled throughout West Asia until their arrival in Europe around the eleventh century CE. Their diaspora through Europe is characterized by population bottlenecks and founder events which have contributed to their present day genetic and cultural diversity. In our study, we focus on the effects of founder effects in the mitochondrial DNA (mtDNA) pool of Iberian Roma by producing and analyzing 144 novel whole mtDNA sequences of Iberian Roma. Over 60% of their mtDNA pool is composed by founder lineages of South Asian origin or acquired by gene flow during their diaspora in the Middle East or locally in Europe in Europe. The TMRCA of these lineages predates the historical record of the Roma arrival in Spain. The abundance of founder lineages is in contrast with ~ 0.7% of autochthonous founder lineages present in the non-Roma Iberian population. Within those founder lineages, we found a substantial amount of South Asian M5a1b1a1 haplotypes and high frequencies of West Eurasian founder lineages (U3b1c, J2b1c, J1c1b, J1b3a, H88, among others), which we characterized phylogenetically and put in phylogeographical context. Besides, we found no evidence of genetic substructure of Roma within the Iberian Peninsula. These results show the magnitude of founder effects in the Iberian Roma and further explain the Roma history and genetic diversity from a matrilineal point of view.

Eurasiaplex-2: Shifting the focus to SNPs with high population specificity increases the power of forensic ancestry marker sets

To compile a new South Asian-informative panel of forensic ancestry SNPs, we changed the strategy for selecting the most powerful markers for this purpose by targeting polymorphisms with near absolute specificity - when the South Asian-informative allele identified is absent from all other populations or present at frequencies below 0.001 (one in a thousand). More than 120 candidate SNPs were identified from 1000 Genomes datasets satisfying an allele frequency screen of ≥ 0.1 (10 % or more) allele frequency in South Asians, and ≤ 0.001 (0.1 % or less) in African, East Asian, and European populations. From the candidate pool of markers, a final panel of 36 SNPs, widely distributed across most autosomes, were selected that had allele frequencies in the five 1000 Genomes South Asian populations ranging from 0.4 to 0.15. Slightly lower average allele frequencies, but consistent patterns of informativeness were observed in gnomAD South Asian datasets used to validate the 1000 Genomes variant annotations. We named the panel of 36 South Asian-specific SNPs Eurasiaplex-2, and the informativeness of the panel was evaluated by compiling worldwide population data from 4097 samples in four genome variation databases that largely complement the global sampling of 1000 Genomes. Consistent patterns of allele frequency distribution, which were specific to South Asia, were observed in all populations in, or closely sited to, the Indian sub-continent. Pakistani populations from the HGDP-CEPH panel had markedly lower allele frequencies, highlighting the need to develop a statistical system to evaluate the ancestry inference value of counting the number of population-specific alleles present in an individual.

Paternal genetic structure of the Bosnian-Herzegovinian Roma: A Y-chromosomal STR study

OBJECTIVES

Studies indicate the complex nature of the genetic structure of the European Roma which has been shaped by different effects of their demographic history, while preserving their ancestral Indian origin. The primary aims of this study were to present for the first time the paternal profiles of the Roma from Bosnia and Herzegovina based on the data from Y-chromosome STR loci, identify the components of non-Roma paternal gene flow into the Roma, and evaluate the genetic relationships with other European Roma populations.

MATERIALS AND METHODS

In this study, 110 DNA samples of unrelated males from Roma populations residing in different regions of Bosnia and Herzegovina were genotyped using the 23 Y-STR loci included in the PowerPlex Y23 system.

RESULTS

The analysis of the genetic structure of the Bosnian-Herzegovinian Roma revealed intra-country population substructuring and indicated differing genetic affinities between the Bosnian-Herzegovinian Roma and other European Roma populations. The paternal genetic structure of the Bosnian-Herzegovinian Roma has two components: an ancestral component represented by haplogroup H1a1a-M82, and European component presented by haplogroups I1-M253, I2a1a2b-L621, J2a1a-L26, J2a1a1a2b2a3~Z7671, J2b2a-M241, G2a2b2a1a1b-L497, and E1b1b-M215.

CONCLUSION

Genetic relations between the Bosnian-Herzegovinian Roma and other European Roma are shaped by different influences on their demographic history. The data suggest that the paternal gene pool of the Roma from Bosnia and Herzegovina might be a consequence of an early separation of the proto-Roma population and the later gene flow as well as factors of the isolation that accompany the Roma populations in some Bosnian-Herzegovinian regions.

The genetic scenario of Mercheros: an under-represented group within the Iberian Peninsula

BACKGROUND

The general picture of human genetic variation has been vastly depicted in the last years, yet many populations remain broadly understudied. In this work, we analyze for the first time the Merchero population, a Spanish minority ethnic group that has been scarcely studied and historically persecuted. Mercheros have been roughly characterised by an itinerant history, common traditional occupations, and the usage of their own language.

RESULTS

Here, we examine the demographic history and genetic scenario of Mercheros, by using genome-wide array data, whole mitochondrial sequences, and Y chromosome STR markers from 25 individuals. These samples have been complemented with a wide-range of present-day populations from Western Eurasia and North Africa. Our results show that the genetic diversity of Mercheros is explained within the context of the Iberian Peninsula, evidencing a modest signal of Roma admixture. In addition, Mercheros present low genetic isolation and intrapopulation heterogeneity.

CONCLUSIONS

This study represents the first genetic characterisation of the Merchero population, depicting their fine-scale ancestry components and genetic scenario within the Iberian Peninsula. Since ethnicity is not only influenced by genetic ancestry but also cultural factors, other studies from multiple disciplines are needed to further explore the Merchero population. As with Mercheros, there is a considerable gap of underrepresented populations and ethnic groups in publicly available genetic data. Thus, we encourage the consideration of more ethnically diverse population panels in human genetic studies, as an attempt to improve the representation of human populations and better reconstruct their fine-scale history.

Representing vulnerable populations in genetic studies: The case of the Roma

Moreau () has raised concerns about the use of DNA data obtained from vulnerable populations, such as the Uighurs in China. We discuss another case, situated in Europe and with a research history dating back 100 years: genetic investigations of Roma. In our article, we focus on problems surrounding representativity in these studies. We claim that many of the circa 440 publications in our sample neglect the methodological and conceptual challenges of representativity. Moreover, authors do not account for problematic misrepresentations of Roma resulting from the conceptual frameworks and sampling schemes they use. We question the representation of Roma as a "genetic isolate" and the underlying rationales, with a strong focus on sampling strategies. We discuss our results against the optimistic prognosis that the "new genetics" could help to overcome essentialist understandings of groups.

Case Report: Two Monochorionic Twins With a Critically Different Course of Progressive Osseus Heteroplasia

Progressive osseous heteroplasia (POH; OMIM 166350) is a rare autosomal-dominant genetic disorder in which extra-skeletal bone forms within skin and muscle tissue. POH is one of the clinical manifestations of an inactivating mutation in the GNAS gene. GNAS gene alterations are difficult matter to address, as GNAS alleles show genetic imprinting and produce several transcript products, and the same mutation may lead to strikingly different phenotypes. Also, most of the publications concerning POH patients are either clinical depictions of a case (or a case series), descriptions of their genetic background, or a tentative correlation of both clinical and molecular findings. Treatment for POH is rarely addressed, and POH still lacks therapeutic options. We describe a unique case of POH in two monochorionic twins, who presented an almost asymptomatic vs. the severe clinical course, despite sharing the same mutation and genetic background. We also report the results of the therapeutic interventions currently available for heterotopic ossification in the patient with the severe course. This article not only critically supports the assumption that the POH course is strongly influenced by factors beyond genetic background but also remarks the lack of options for patients suffering an orphan disease, even after testing drugs with promising in vitro results.

Sex-biased patterns shaped the genetic history of Roma

The Roma population is a European ethnic minority characterized by recent and multiple dispersals and founder effects. After their origin in South Asia around 1,500 years ago, they migrated West. In Europe, they diverged into ethnolinguistically distinct migrant groups that spread across the continent. Previous genetic studies based on genome-wide data and uniparental markers detected Roma founder events and West-Eurasian gene flow. However, to the best of our knowledge, it has not been assessed whether these demographic processes have equally affected both sexes in the population. The present study uses the largest and most comprehensive dataset of complete mitochondrial and Y chromosome Roma sequences to unravel the sex-biased patterns that have shaped their genetic history. The results show that the Roma maternal genetic pool carries a higher lineage diversity from South Asia, as opposed to a single paternal South Asian lineage. Nonetheless, the European gene flow events mainly occurred through the maternal lineages; however, a signal of this gene flow is also traceable in the paternal lineages. We also detect a higher female migration rate among European Roma groups. Altogether, these results suggest that sociocultural factors influenced the emergence of sex-biased genetic patterns at global and local scales in the Roma population through time.

Analysis of the mitochondrial CYTB gene sequence in human populations of northeastern Bosnia

This study offers the first report on variation sequence of the mitochondrial cytochrome b (MTCYTB) gene in populations from Bosnia (northeastern Bosnia). This study was designed on the analysis of the genetic diversity of two populations of different cultural-anthropological and genetic origin, Roma population and native/non-Roma population. The main aim of our study was to estimate the usefulness of the CYTB sequence in the analysis of genetic categorization of different populations and intergroup diversity, as well as to provide some additional information on haplogroup-associated polymorphisms within the CYTB region in defining haplogroup status. Estimation of the genetic diversity was done using intra and intergroup genetic indices. The population-specific polymorphisms have been found in both categories of the populations. The results of the analysis of genetic differentiation show significant pairwise Fst differences between the Romani and native populations. Also, registered significant genetic differentiation is illustrated on the level of genetic variation between subpopulations of the Roma and non-Roma origin. The important result in our study is the confirmation of the significance of the triad of polymorphisms T14783C-G15043A-G15301A, indicating the influence of Asian component of the maternal gene pool on the genetic structure of the studied population of the Roma. Our data show that the haplogroup polymorphisms exist in the CYTB region and can provide useful information on the haplogroups that were defined only by the control region of the mtDNA. The results of this study indicate the region of CYTB gene can be a benefit in providing some additional information in the analysis of genetic structure of human populations and can be additionally applied in population studies.

Mitochondrial DNA diversity of the Roma from northeastern Bosnia, Bosnia and Herzegovina

This study is the first report on the mtDNA diversity in the Roma population from Bosnia. The main aim of this study was to analyse the mtDNA diversity in the studied population, evaluate the genetic relations with other European Roma populations, and analyse the influences of the Roma gene flow on the mitochondrial profile of the Roma from northeastern Bosnia. MtDNA variability in the analysed population has been studied by means of hypervariable segment I and II (HVSI/II) of the control region sequencing and analysis of restriction fragment-length polymorphisms of the coding region. Our results show that genetic structure of the Roma from northeastern Bosnia has a combination of lineages of three main layers: specific founder of Indian origin (M5a1 and M35b) and founder non-M lineages of Indian/European origin (H7a, X2b and X2d) and lineages of European/Middle East origin (H, H1, H11a, V, T2b, K1b and W). The distribution of the haplogroups in the gene pool of the comparative European Roma populations indicates the separate origin of the Bosnian Roma and the Bulgarian Vlax and Croatian Vlax Roma. The data suggest that mitochondrial gene pool of the Roma population from northeastern Bosnia might be a consequence of early parting and the later different migration routes that are part of their demographic history. Our data confirm the high genetic heterogeneity of the Roma populations that can be shaped by effects of genetic drift, isolation and low effective population size, and this correlates with the migratory history of the Roma.

Whole Exome Sequencing Identifies New Host Genomic Susceptibility Factors in Empyema Caused by Streptococcus pneumoniae in Children: A Pilot Study

Pneumonia is the leading cause of death amongst infectious diseases. Streptococcus pneumoniae is responsible for about 25% of pneumonia cases worldwide, and it is a major cause of childhood mortality. We carried out a whole exome sequencing (WES) study in eight patients with complicated cases of pneumococcal pneumonia (empyema). An initial assessment of statistical association of WES variation with pneumonia was carried out using data from the 1000 Genomes Project (1000G) for the Iberian Peninsula (IBS) as reference controls. Pseudo-replication statistical analyses were carried out using different European control groups. Association tests pointed to single nucleotide polymorphism (SNP) rs201967957 (gene MEIS1; chromosome 2; p-value_IBS = 3.71 × 10^-13) and rs576099063 (gene TSPAN15; chromosome 10; p-value_IBS = 2.36 × 10^-8) as the best candidate variants associated to pneumococcal pneumonia. A burden gene test of pathogenicity signaled four genes, namely, OR9G9, MUC6, MUC3A and APOB, which carry significantly increased pathogenic variation when compared to controls. By analyzing various transcriptomic data repositories, we found strong supportive evidence for the role of MEIS1, TSPAN15 and APOBR (encoding the receptor of the APOB protein) in pneumonia in mouse and human models. Furthermore, the association of the olfactory receptor gene OR9G9 has recently been related to some viral infectious diseases, while the role of mucin genes (MUC6 and MUC3A), encoding mucin glycoproteins, are well-known factors related to chronic obstructive airway disease. WES emerges as a promising technique to disentangle the genetic basis of host genome susceptibility to infectious respiratory diseases.

Whole Exome Sequencing reveals new candidate genes in host genomic susceptibility to Respiratory Syncytial Virus Disease

Respiratory syncytial virus (RSV) is an important cause of serious lower respiratory tract disease in infants. Several studies have shown evidence pointing to the genome of the host as an important factor determining susceptibility to respiratory disease caused by RSV. We sequenced the complete exomes of 54 patients infected by RSV that needed hospitalization due to development of severe bronchiolitis. The Iberian sample (IBS) from The 1000 Genomes Project (1000G) was used as control group; all the association results were pseudo-replicated using other 1000G-European controls and Spanish controls. The study points to SNP rs199665292 in the olfactory receptor (OR) gene OR13C5 as the best candidate variant (P-value = 1.16 × 10-12; OR = 5.56). Genetic variants at HLA genes (HLA-DQA1, HLA-DPB1), and in the mucin 4 gene (MUC4) also emerge as susceptibility candidates. By collapsing rare variants in genes and weighing by pathogenicity, we obtained confirmatory signals of association in the OR gene OR8U1/OR8U8, the taste receptor TAS2R19, and another mucin gene (MUC6). Overall, we identified new predisposition variants and genes related to RSV infection. Of special interest is the association of RSV to olfactory and taste receptors; this finding is in line with recent evidence pointing to their role in viral infectious diseases.

A genetic chronology for the Indian Subcontinent points to heavily sex-biased dispersals

BACKGROUND

India is a patchwork of tribal and non-tribal populations that speak many different languages from various language families. Indo-European, spoken across northern and central India, and also in Pakistan and Bangladesh, has been frequently connected to the so-called "Indo-Aryan invasions" from Central Asia ~3.5 ka and the establishment of the caste system, but the extent of immigration at this time remains extremely controversial. South India, on the other hand, is dominated by Dravidian languages. India displays a high level of endogamy due to its strict social boundaries, and high genetic drift as a result of long-term isolation which, together with a very complex history, makes the genetic study of Indian populations challenging.

RESULTS

We have combined a detailed, high-resolution mitogenome analysis with summaries of autosomal data and Y-chromosome lineages to establish a settlement chronology for the Indian Subcontinent. Maternal lineages document the earliest settlement ~55-65 ka (thousand years ago), and major population shifts in the later Pleistocene that explain previous dating discrepancies and neutrality violation. Whilst current genome-wide analyses conflate all dispersals from Southwest and Central Asia, we were able to tease out from the mitogenome data distinct dispersal episodes dating from between the Last Glacial Maximum to the Bronze Age. Moreover, we found an extremely marked sex bias by comparing the different genetic systems.

CONCLUSIONS

Maternal lineages primarily reflect earlier, pre-Holocene processes, and paternal lineages predominantly episodes within the last 10 ka. In particular, genetic influx from Central Asia in the Bronze Age was strongly male-driven, consistent with the patriarchal, patrilocal and patrilineal social structure attributed to the inferred pastoralist early Indo-European society. This was part of a much wider process of Indo-European expansion, with an ultimate source in the Pontic-Caspian region, which carried closely related Y-chromosome lineages, a smaller fraction of autosomal genome-wide variation and an even smaller fraction of mitogenomes across a vast swathe of Eurasia between 5 and 3.5 ka.

Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia

Background

From a mtDNA dominant perspective, the exit from Africa of modern humans to colonize Eurasia occurred once, around 60 kya, following a southern coastal route across Arabia and India to reach Australia short after. These pioneers carried with them the currently dominant Eurasian lineages M and N. Based also on mtDNA phylogenetic and phylogeographic grounds, some authors have proposed the coeval existence of a northern route across the Levant that brought mtDNA macrohaplogroup N to Australia. To contrast both hypothesis, here we reanalyzed the phylogeography and respective ages of mtDNA haplogroups belonging to macrohaplogroup M in different regions of Eurasia and Australasia.

Results

The macrohaplogroup M has a historical implantation in West Eurasia, including the Arabian Peninsula. Founder ages of M lineages in India are significantly younger than those in East Asia, Southeast Asia and Near Oceania. Moreover, there is a significant positive correlation between the age of the M haplogroups and its longitudinal geographical distribution. These results point to a colonization of the Indian subcontinent by modern humans carrying M lineages from the east instead the west side.

Conclusions

The existence of a northern route, previously proposed for the mtDNA macrohaplogroup N, is confirmed here for the macrohaplogroup M. Both mtDNA macrolineages seem to have differentiated in South East Asia from ancestral L3 lineages. Taking this genetic evidence and those reported by other disciplines we have constructed a new and more conciliatory model to explain the history of modern humans out of Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0816-8) contains supplementary material, which is available to authorized users.

Genomic continuity of Argentinean Mennonites

Mennonites are Anabaptist communities that originated in Central Europe about 500 years ago. They initially migrated to different European countries, and in the early 18^th century they established their first communities in North America, from where they moved to other American regions. We aimed to analyze an Argentinean Mennonite congregation from a genome-wide perspective by way of investigating >580.000 autosomal SNPs. Several analyses show that Argentinean Mennonites have European ancestry without signatures of admixture with other non-European American populations. Among the worldwide datasets used for population comparison, the CEU, which is the best-subrogated Central European population existing in The 1000 Genome Project, is the dataset showing the closest genome affinity to the Mennonites. When compared to other European population samples, the Mennonites show higher inbreeding coefficient values. Argentinean Mennonites show signatures of genetic continuity with no evidence of admixture with Americans of Native American or sub-Saharan African ancestry. Their genome indicates the existence of an increased endogamy compared to other Europeans most likely mirroring their lifestyle that involve small communities and historical consanguineous marriages.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Mapping the genomic mosaic of two 'Afro-Bolivians' from the isolated Yungas valleys

Background

Unraveling the ancestry of ‘Afro-American’ communities is hampered by the complex demographic processes that took place during the Transatlantic Slave Trade (TAST) and the (post-)colonization periods. ‘Afro-Bolivians’ from the subtropical Yungas valleys constitute small and isolated communities that live surrounded by the predominant Native American community of Bolivia. By genotyping >580,000 SNPs in two ‘Afro-Bolivians’, and comparing these genomic profiles with data compiled from more than 57 African groups and other reference ancestral populations (n = 1,161 in total), we aimed to disentangle the complex admixture processes undergone by ‘Afro-Bolivians’.

Results

The data indicate that these two genomes constitute a complex mosaic of ancestries that is approximately 80 % of recent African origin; the remaining ~20 % being European and Native American. West-Central Africa contributed most of the African ancestry to ‘Afro-Bolivians’, and this component is related to populations living along the Atlantic coast (i.e. Senegal, Ghana, Nigeria). Using tract length distribution of genomic segments attributable to distinct ancestries, we could date the time of admixture in about 400 years ago. This time coincides with the maximum importation of slaves to Bolivia to compensate the diminishing indigenous labor force needed for the development of the National Mint of Potosí.

Conclusions

Overall, the data indicate that the genome of ‘Afro-Bolivians’ was shaped by a complex process of admixture occurring in America among individuals originating in different West-Central African populations; their genomic mosaics received additional contributions of Europeans and local Native Americans (e.g. Aymaras).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2520-x) contains supplementary material, which is available to authorized users.

Origins, admixture and founder lineages in European Roma

The Roma, also known as 'Gypsies', represent the largest and the most widespread ethnic minority of Europe. There is increasing evidence, based on linguistic, anthropological and genetic data, to suggest that they originated from the Indian subcontinent, with subsequent bottlenecks and undetermined gene flow from/to hosting populations during their diaspora. Further support comes from the presence of Indian uniparentally inherited lineages, such as mitochondrial DNA M and Y-chromosome H haplogroups, in a significant number of Roma individuals. However, the limited resolution of most genetic studies so far, together with the restriction of the samples used, have prevented the detection of other non-Indian founder lineages that might have been present in the proto-Roma population. We performed a high-resolution study of the uniparental genomes of 753 Roma and 984 non-Roma hosting European individuals. Roma groups show lower genetic diversity and high heterogeneity compared with non-Roma samples as a result of lower effective population size and extensive drift, consistent with a series of bottlenecks during their diaspora. We found a set of founder lineages, present in the Roma and virtually absent in the non-Roma, for the maternal (H7, J1b3, J1c1, M18, M35b, M5a1, U3, and X2d) and paternal (I-P259, J-M92, and J-M67) genomes. This lineage classification allows us to identify extensive gene flow from non-Roma to Roma groups, whereas the opposite pattern, although not negligible, is substantially lower (up to 6.3%). Finally, the exact haplotype matching analysis of both uniparental lineages consistently points to a Northwestern origin of the proto-Roma population within the Indian subcontinent.

A genome-wide study of modern-day Tuscans: revisiting Herodotus's theory on the origin of the Etruscans

Background

The origin of the Etruscan civilization (Etruria, Central Italy) is a long-standing subject of debate among scholars from different disciplines. The bulk of the information has been reconstructed from ancient texts and archaeological findings and, in the last few years, through the analysis of uniparental genetic markers.

Methods

By meta-analyzing genome-wide data from The 1000 Genomes Project and the literature, we were able to compare the genomic patterns (>540,000 SNPs) of present day Tuscans (N = 98) with other population groups from the main hypothetical source populations, namely, Europe and the Middle East.

Results

Admixture analysis indicates the presence of 25–34% of Middle Eastern component in modern Tuscans. Different analyses have been carried out using identity-by-state (IBS) values and genetic distances point to Eastern Anatolia/Southern Caucasus as the most likely geographic origin of the main Middle Eastern genetic component observed in the genome of modern Tuscans.

Conclusions

The data indicate that the admixture event between local Tuscans and Middle Easterners could have occurred in Central Italy about 2,600–3,100 years ago (y.a.). On the whole, the results validate the theory of the ancient historian Herodotus on the origin of Etruscans.