An Ancient Mediterranean Melting Pot: Investigating the Uniparental Genetic Structure and Population History of Sicily and Southern Italy

Due to their strategic geographic location between three different continents, Sicily and Southern Italy have long represented a major Mediterranean crossroad where different peoples and cultures came together over time. However, its multi-layered history of migration pathways and cultural exchanges, has made the reconstruction of its genetic history and population structure extremely controversial and widely debated. To address this debate, we surveyed the genetic variability of 326 accurately selected individuals from 8 different provinces of Sicily and Southern Italy, through a comprehensive evaluation of both Y-chromosome and mtDNA genomes. The main goal was to investigate the structuring of maternal and paternal genetic pools within Sicily and Southern Italy, and to examine their degrees of interaction with other Mediterranean populations. Our findings show high levels of within-population variability, coupled with the lack of significant genetic sub-structures both within Sicily, as well as between Sicily and Southern Italy. When Sicilian and Southern Italian populations were contextualized within the Euro-Mediterranean genetic space, we observed different historical dynamics for maternal and paternal inheritances. Y-chromosome results highlight a significant genetic differentiation between the North-Western and South-Eastern part of the Mediterranean, the Italian Peninsula occupying an intermediate position therein. In particular, Sicily and Southern Italy reveal a shared paternal genetic background with the Balkan Peninsula and the time estimates of main Y-chromosome lineages signal paternal genetic traces of Neolithic and post-Neolithic migration events. On the contrary, despite showing some correspondence with its paternal counterpart, mtDNA reveals a substantially homogeneous genetic landscape, which may reflect older population events or different demographic dynamics between males and females. Overall, both uniparental genetic structures and TMRCA estimates confirm the role of Sicily and Southern Italy as an ancient Mediterranean melting pot for genes and cultures.     

Genetic structure of the population of Sicily.

Genetic heterogeneity within Sicily was investigated on the basis of ACP1, ADA, ESD, GLO1, PGD, PGM1, PGM2, SODA, ABO, and MN gene frequencies, and compared to those of other regions of Italy for which these same loci have been examined. Correspondence analysis revealed no differences within the island, at least at the provincial level, but showed genetic differentiation among Italian regions, distinctly clustering northern, central, and southern populations, respectively. These data indicate a close relationship between Sicily and southern Italy. In addition, the contribution of Middle Eastern populations to the gene pool of Sicily was evident.     

Properdin factor B polymorphism in continental Italy and Sardinia

Summary A sample of healthy unrelated individuals (882 from Continental Italy and Sicily, and 139 from Sardinia) were typed for HLA and Bf polymorphisms. Bf gene frequencies in Continental Italy and Sicily show an intermediate position among European populations, with Lapps on one side and Basques on the other side; while a very peculiar pattern is observed among Sardinians, with an impressive similarity with French Basques and with the highest F1 gene frequency so far observed. The HLA: Bf gametic associations in Italy and Sardinia are also reported, which confirm the different ethnic background of these populations.This work was supported by Centro di Immunogenetica ed Istocompatibilità, C.N.R., Torino, Italy and by a grant of Ministero Pubblica Istruzione (M.P.I., Cap. 8551/1980).     

Genetic analysis of a Sicilian population using 15 short tandem repeats

The genetic structure of the population of Alia (Sicily, Italy) was analyzed using 15 short tandem repeats: TPOX, D2S1338, D3S1358, FIBRA, D5S818, CSF1PO, D7S820, D8S1179, TH01, VWA, D13S317, D16S539, D18S51, D19S433, and D21S11. Two of these markers, D2S1338 and D19S433, have never before been used in research on population genetics and only recently have they been put to use in forensic medicine. Results of the analysis underline the genetic isolation of the Alia population and show it to be a recent bottleneck as a consequence of a cholera epidemic in 1837. While comparing the Alia population with other populations from Sicily, a genetic heterogeneity within Sicily was uncovered, thus confirming previous results obtained from the analysis of classical markers. This heterogeneity underlines the existence of genetic boundaries within the island. Comparisons with other Italian, Mediterranean, and European populations highlight the differentiation of the Sicilian population, reflecting the presence of a genetic boundary that separates Sicily from northern and central Italy and from the western Mediterranean basin.     

Red cell and serum protein polymorphisms in Sardinia

Four Sardinian population samples, from the provinces of Cagliari, Oristano, Nuoro and Sassari, were studied with regard to the erythrocyte enzyme systems ACP, ESD, PGM1, ADA, AK, 6PGD and Dia, and to the serum protein systems GC and C3. The findings showed a rather high degree of genetic heterogeneity of the Sardinians compared to the other populations from the Mediterranean area (Continental Italy, Sicily, Spain, North Africa).     

The HLA system in Italy

4,902 Italians were typed for HLA-A antigens, 4,721 for HLA-B and 1,503 for HLA-C. The samples, which were composed of unrelated, healthy individuals born in Italy, were used for estimating HLA-A, HLA-B and HLA-C gene frequencies with the maximum-likelihood method. Different Italian regions showed significant differences in the HLA alleles, providing further evidence for the genetic heterogeneity of the Italian population. HLA gene frequencies place continental Italy and Sicily in a position which is similar to that of other Mediterranean populations, whereas the genetic isolation of Sardinia is quite evident. The most significant linkage disequilibrium values found in the Italian population (except for Sardinia) were in agreement with those observed in other Caucasian populations. The difference between Northern and Southern Italy and between continental Italy and Sardinia was emphasized by the linkage disequilibrium values and by the principal-component analysis as well.     

Genetic variation inOrchis papilionacea (Orchidaceae) from the Central Mediterranean region: Taxonomic inferences at the intraspecific level

The genetic differentiation ofOrchis papilionacea from Italy, Corsica, Sardinia, Sicily, and Greece was analyzed by isozyme studies. No significant differences at any of the 28 loci tested were observed between samples of the putative subspeciespapilionacea from Italy,grandiflora from Sicily andheroica from Greece, or by comparing sympatricpapilionacea-like andgrandiflora-like subsamples from Sardinia and Corsica. The genetic relationships found among the populations studied reflect more their geographic origin than their taxonomic assignment;papilionacea, grandiflora, andheroica are, therefore, at least in the study area, not to be considered as subspecies, but as simple morphs, controlled by few genes (possibly two) with Mendelian inheritance. The genetic divergence observed among the samples studied is very low and consistent with that of geographic populations of a single cohesive gene pool, linked by a high gene flow (Nm = 5.9).     

Phylogeographic inferences from the mtDNA variation of the three-toed skink, Chalcides chalcides (Reptilia: Scincidae)

Genetic diversity was analyzed in Chalcides chalcides populations from peninsular Italy, Sardinia, Sicily and Tunisia by sequencing 400 bp at the 5' end of the mitochondrial gene encoding cytochrome b (cyt b) and by restriction fragment length polymorphism (RFLP) analysis of two mitochondrial DNA segments (ND-1/2 and ND-3/4). The results of the phylogenetic analysis highlighted the presence of three main clades corresponding with three of the four main geographical areas (Tunisia, Sicily and the Italian peninsula), while Sardinia proved to be closely related to Tunisian haplotypes suggesting a colonization of this island from North Africa by human agency in historical times. On the contrary, the splitting times estimated on the basis of cyt b sequence data seem to indicate a more ancient colonization of Sicily and the Italian Peninsula, as a consequence of tectonic and climatic events that affected the Mediterranean Basin during the Pleistocene. Finally, the analysis of the genetic variability of C. chalcides populations showed a remarkable genetic homogeneity in Italian populations when compared to the Tunisian ones. This condition could be explained by a rapid post-glacial expansion from refugial populations that implied serial bottlenecking with progressive loss of haplotypes, resulting in a low genetic diversity in the populations inhabiting the more recently colonized areas.     

Late Pleistocene human evolution in Sicily: comparative morphometric analysis of Grotta di San Teodoro craniofacial remains

The paleoanthropological remains from Grotta di San Teodoro near Acquedolci (province of Messina, Italy) represent the oldest and largest skeletal collection yet found documenting human settlement of Sicily. The sample, attributed to the Late Epigravettian (between 14,000 and 10,000 years B.P.), consists of seven variously complete adult individuals (San Teodoro 1–7). We compare the cranial sample to an array of both prehistoric and recent samples using multivariate techniques including D² distance analysis, canonical variate analysis, cluster analysis, and multidimensional scaling. Overall, the San Teodoro cranial sample displays a morphometric pattern close to Western European groups of similar antiquity, in particular those from Central and Southern Italy. The morphometric affinities indicate that these people probably came from peninsular Italy by sea during the Late Epigravettian epoch. An alternative hypothesis is that they descended from immigrants that arrived by land during a low sea level episode corresponding to the maximum Würmian regression, about 18,000 years B.P, with gene flow accounting for the morphological homogeneity with the populations of peninsular Italy. The San Teodoro skeletal sample provides the first reliable evidence for human settlement of Sicily.     

Peopling of three Mediterranean islands (Corsica,Sardinia, and Sicily) inferred by Y-chromosome biallelic variability

An informative set of biallelic polymorphisms was used to study the structure of Y-chromosome variability in a sample from the Mediterranean islands of Corsica and Sicily, and compared with data on Sardinia to gain insights into the ethnogenesis of these island populations. The results were interpreted in a broader Mediterranean context by including in the analysis neighboring populations previously studied with the same methodology. All samples studied were enclosed in the comparable spectrum of European Y-chromosome variability. Pronounced differences were observed between the islands as well as in the percentages of haplotypes previously shown to have distinctive patterns of continental phylogeography. Approximately 60% of the Sicilian haplotypes are also prevalent in Southern Italy and Greece. Conversely, the Corsican sample had elevated levels of alternative haplotypes common in Northern Italy. Sardinia showed a haplotype ratio similar to that observed in Corsica, but with a remarkable difference in the presence of a lineage defined by marker M26, which approaches 35% in Sardinia but seems absent in Corsica. Although geographically adjacent, the data suggest different colonization histories and a minimal amount of recent gene flow between them. Our results identify possible ancestral continental sources of the various island populations and underscore the influence of founder effect and genetic drift. The Y-chromosome data are consistent with comparable mtDNA data at the RFLP haplogroup level of resolution, as well as linguistic and historic knowledge.     

Phylogeography of Hyalomma (Euhyalomma) lusitanicum (Acarina,Parasitiformes, Ixodidae) in Andalusia based on mitochondrial cytochrome oxidase I gene

The genetic population structure relationships of Hyalomma (Euhyalomma) lusitanicum in Andalusia (the south of the Iberian Peninsula) were examined using mtDNA sequence data from 887 bp of cytochrome oxidase subunit I (COI) gene. The sequence for the COI region was determined for 84 individuals collected in several localities of Andalusia, and 10 for other localities (i.e., five from Toledo, central Iberian Peninsula, four from Sicily (Italy) and one from Canary Island). Seventeen haplotypes were detected, including 27 polymorphic sites. The number of amino acid substitutions per site from mean diversity calculations for the entire population was 0.017. AMOVA analysis revealed a low gene flow that characterises the genetic population structure of this species in South Iberian Peninsula, with a haplotype diversity (h) value of 0.815. No geographically induced differentiation was observed, and separate evolutionary units were not detected. Our results indicate low genetic diversity across the geographical range of H. lusitanicum tick in Andalusia. Our data do not show any genetic discontinuity between the tick populations studied, including specimens from Canary Island and Sicily (Italy).

     

Emergence and Phylodynamics of Citrus tristeza virus in Sicily,Italy

Citrus tristeza virus (CTV) outbreaks were detected in Sicily island, Italy for the first time in 2002. To gain insight into the evolutionary forces driving the emergence and phylogeography of these CTV populations, we determined and analyzed the nucleotide sequences of the p20 gene from 108 CTV isolates collected from 2002 to 2009. Bayesian phylogenetic analysis revealed that mild and severe CTV isolates belonging to five different clades (lineages) were introduced in Sicily in 2002. Phylogeographic analysis showed that four lineages co-circulated in the main citrus growing area located in Eastern Sicily. However, only one lineage (composed of mild isolates) spread to distant areas of Sicily and was detected after 2007. No correlation was found between genetic variation and citrus host, indicating that citrus cultivars did not exert differential selective pressures on the virus. The genetic variation of CTV was not structured according to geographical location or sampling time, likely due to the multiple introduction events and a complex migration pattern with intense co- and re-circulation of different lineages in the same area. The phylogenetic structure, statistical tests of neutrality and comparison of synonymous and nonsynonymous substitution rates suggest that weak negative selection and genetic drift following a rapid expansion may be the main causes of the CTV variability observed today in Sicily. Nonetheless, three adjacent amino acids at the p20 N-terminal region were found to be under positive selection, likely resulting from adaptation events.     

An Overview of the Genetic Structure within the Italian Population from Genome-Wide Data

In spite of the common belief of Europe as reasonably homogeneous at genetic level, advances in high-throughput genotyping technology have resolved several gradients which define different geographical areas with good precision. When Northern and Southern European groups were considered separately, there were clear genetic distinctions. Intra-country genetic differences were also evident, especially in Finland and, to a lesser extent, within other European populations. Here, we present the first analysis using the 125,799 genome-wide Single Nucleotide Polymorphisms (SNPs) data of 1,014 Italians with wide geographical coverage. We showed by using Principal Component analysis and model-based individual ancestry analysis, that the current population of Sardinia can be clearly differentiated genetically from mainland Italy and Sicily, and that a certain degree of genetic differentiation is detectable within the current Italian peninsula population. Pair-wise F_ST statistics Northern and Southern Italy amounts approximately to 0.001 between, and around 0.002 between Northern Italy and Utah residents with Northern and Western European ancestry (CEU). The Italian population also revealed a fine genetic substructure underscoring by the genomic inflation (Sardinia vs. Northern Italy = 3.040 and Northern Italy vs. CEU = 1.427), warning against confounding effects of hidden relatedness and population substructure in association studies.     

Strong genetic cohesiveness between Italy and North Africa in four butterfly species

The sea acts as an effective dispersal barrier for most terrestrial animal species. Narrow sea straits, therefore, often represent areas where species are able to disperse from one land mass to another. In the Mediterranean Sea, the narrowest connecting points between North Africa and Europe are the Strait of Gibraltar and the Strait of Sicily. In the past, climatic oscillations caused changing sea levels and thus influenced the permeability of these sea straits. We analysed the genetic structure of four butterfly species that all occur on both sides of the Strait of Sicily. In all four species, we observed a lack of genetic differentiation between the populations of North Africa and those of Italy. Species distribution models support the strong cohesiveness in that they show a largely continuous glacial distribution over Italy and North Africa. The data obtained reveal that there was a large exchange of individuals between Italy and the eastern Maghreb during the last ice age. This might not only be the case for the species under investigation in the present study, but also might represent a more general pattern for mobile thermophilic western Palearctic species. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 818–830.     

Argas (Persicargas) persicus (Oken, 1818) (Ixodida: Argasidae) in Sicily with considerations about its Italian and West-Mediterranean distribution

Recently, in the province of Trapani (Western Sicily), some overwintering specimens of the argasid tick Argos (Persicargas) persicus (Oken, 1818) were observed and collected. Morphological and genetic analysis were utilized in order to reach a definitive identification. The species was found in two semi-natural sites where, having been found repeatedly, its presence does not appear accidental. Moreover the characteristics of the Sicilian findings seem to exclude a human-induced spread. This record, the first regarding Sicily and South Italy, is discussed together with the previous doubtful citations for Italy. These findings revalue not only all the old citations for Italy but also the hypothesis that the Mediterranean distribution of this argasid is of a natural origin.     

Distribution of ORM1, C6, C7 and APO C-II allele frequencies in populations from mainland Italy and Sardinia.

The genetic variation of the human plasma proteins ORM1, C6, C7 and APO C-II was investigated by isoelectric focusing followed by immunoblotting in populations from mainland Italy and Sardinia. The frequencies of ORM1*1 were 0.621 and 0.564, while those of C6*A were 0.657 and 0.706 on mainland Italy and in Sardinia, respectively. In the Roman sample, 8 heterozygotes with C6 variant alleles were encountered, while none were observed in Sardinians. For C7 and APO C-II a number of heterozygotes with the rare alleles C7*2 and APO C-II*2 were found, but their frequency did not reach polymorphic levels in either population. The two populations showed a significant difference in the gene frequencies distribution for ORM1.     

The response of European beech (Fagus sylvatica L.) seedlings from two Italian populations to drought and recovery

The response of gas exchange, leaf chlorophyll content, relative fluorescence to decreasing water potential and recovery was followed in European beech seedlings from two Italian populations, differing in their native precipitation amounts. A population from Sicily (southern Italy) was selected as representative of a xeric population while a population from central Italy, Abetone, represented a mesic one. Dry-matter partitioning, leaf area, hydraulic sufficiency and xylem embolism were evaluated in both well-watered control plants as well as in plants subjected to drought. With the onset of water stress, values of water potential, leaf relative water content, net photosynthesis, leaf conductance and leaf chlorophyll concentration decreased concurrently while relative fluorescence remained unchanged. The population from Sicily showed a delay in effects of the imposed drought. Within 5 days of rewatering, leaf conductance was not fully recovered while all of the other parameters recovered to control levels, in both populations. Total, shoot, stem and root dry weight tended to be higher in seedlings from Abetone, even though both populations had similar photosynthetic rates. The population from Sicily exhibited about 3% greater (even if not significant) allocation to roots than the population from Abetone. Seedlings from Abetone had higher, but not significant, leaf specific conductivity and per cent loss in hydraulic conductivity than seedlings from Sicily. Drought resulted in a reduction of hydraulic conductivity and hydraulic sufficiency in seedlings from both populations. Photosynthesis of water stressed plants from both populations appeared to be reduced primarily by carbon dioxide diffusion through stomata and perhaps secondarily by changes in chlorophyll concentration rather than by efficiency of photosystem II. The effect of hydraulic factors on gas exchange during drought and recovery was not clearly evident.     

Does genetic population structure of Ambrosina bassii L. (Araceae,Ambrosineae) attest a post-Messinian land-bridge between Sicily and Africa?

Aim of the present work is the analysis (through the study of enzyme polymorphism) of Sicilian and African (Tunisian) populations of Ambrosina bassii, a small perennial endemic to the Central-Western Mediterranean basin, in order to verify if the complex geological history of this part of the Mediterranean area left its mark in the present-day genetic structure of this taxon. Starch gel allozyme electrophoresis of seven putative loci of A. bassii was employed to estimate genetic diversity, genetic structure and gene flow. Populations from Sicily, Tunisia and Sardinia (as outgroup) were sampled. Results show that Sicily populations have 4 private alleles, Sardinia 3, Tunisia just one. One allele is private to both Sardinia and Tunisia, another one to Sardinia and Sicily. Even if there are no alleles private to Sicily and Tunisia, “cluster analysis” (based on Nei's genetic distances), “non-metric multidimensional scaling” (computed on the basis of a matrix with F_ST values between populations) and Bayesian analysis point out a clear isolation of Sardinian populations, and a greater similarity between Sicilian and Tunisian populations compared to Sardinian ones. The strong genetic affinity between populations from Sicily and Tunisia, considering the very low dispersal ability of the species, gives evidence of a recent continuity between the populations as well as between the two areas. Considering also the estimates of divergence times, a post-Messinian terrestrial connection between the two landmasses can be hypothesized.     

Parasites and the island syndrome: the colonization of the western Mediterranean islands by Heligmosomoides polygyrus (Dujardin, 1845)

Aim Populations of free‐living vertebrates on islands frequently differ from their mainland counterparts by a series of changes in morphometric, life‐history, behavioural, physiological and genetic traits, collectively referred to as the ‘island syndrome’. It is not known, however, whether the ‘island syndrome’ also affects parasitic organisms. The present study establishes the colonization pattern of the Mediterranean islands by the nematode Heligmosomoides polygyrus, a direct and specific parasite of rodent hosts of the Apodemus genus, and evaluates the effects of island colonization by this species on two components of the island syndrome: the loss of genetic diversity and the enlargement of the ecological niche. Location Heligmosomoides polygyrus was sampled on seven western Mediterranean islands ? Corsica, Crete, Elba, Majorca, Minorca, Sardinia and Sicily ? as well as in 20 continental locations covering the Mediterranean basin. Methods The mitochondrial cytochrome b gene (690 base pairs) was sequenced in 166 adult H. polygyrus individuals sampled in the 27 continental and island locations. Phylogenetic reconstructions in distance, parsimony, maximum likelihood and Bayesian posterior probabilities were carried out on the whole cytochrome b gene data set. The levels of nucleotide, haplotype and genetic divergence (Kimura two‐parameter distance estimator) diversities were estimated in each island population and in the various continental lineages. Results Phylogenetic reconstructions show that the mainland origins of H. polygyrus were continental Spain for the Balearic Islands (Majorca, Minorca), northern Italy for the Tyrrhenian Islands (Corsica, Sardinia, Elba), southern Italy for Sicily, and the Balkan region for Crete. A comparison of island H. polygyrus populations with their mainland source populations revealed two characteristic components of the island syndrome in this parasite. First, island H. polygyrus populations display a significant loss of genetic diversity, which is related (r² = 0.73) to the distance separating the island from the mainland source region. Second, H. polygyrus exhibits a niche enlargement following insularization. Indeed, H. polygyrus in Corsica is present in both A. sylvaticus and Mus musculus domesticus, while mainland H. polygyrus populations are present exclusively in Apodemus hosts. Main conclusions Our results show that H. polygyrus has undergone a loss of genetic diversity and a niche (host) enlargement following colonization of the western Mediterranean islands. To our knowledge, this study provides the first evidence for components of the ‘island syndrome’ in a parasitic nematode species.     

Genetic structure in the Mediterranean seagrass Posidonia oceanica: disentangling past vicariance events from contemporary patterns of gene flow

The Mediterranean Sea is a two‐basin system, with the boundary zone restricted to the Strait of Sicily and the narrow Strait of Messina. Two main population groups are recognized in the Mediterranean endemic seagrass Posidonia oceanica, corresponding to the Western and the Eastern basins. To address the nature of the East–West cleavage in P. oceanica, the main aims of this study were: (i) to define the genetic structure within the potential contact zone (i.e. the Strait of Sicily) and clarify the extent of gene flow between the two population groups, and (ii) to investigate the role of present water circulation patterns vs. past evolutionary events on the observed genetic pattern. To achieve these goals, we utilized SSR markers and we simulated, with respect to current regime, the possible present‐day dispersal pattern of Posidonia floating fruits using 28‐day numerical Lagrangian trajectories. The results obtained confirm the presence of the two main population groups, without any indices of reproductive isolation, with the break zone located at the level of the Southern tip of Calabria. The populations in the Strait of Sicily showed higher affinity with Western than with Eastern populations. This pattern of genetic structure probably reflects historical avenues of recolonization from relict glacial areas and past vicariance events, but seems to persist as a result of the low connectivity among populations via marine currents, as suggested by our dispersal simulation analysis.