Comparing patterns of nuclear and mitochondrial divergence in a cryptic species complex: the case of Iberian and North African wall lizards (Podarcis, Lacertidae)

Combining different sources of information is essential for a complete understanding of the process of genetic differentiation between species. The Iberian and North African wall lizard ( Podarcis ) species complex has been the object of several studies regarding morphological and mitochondrial DNA variation but, so far, no large-scale survey of nuclear variation within this group has been accomplished. In this study, ten polymorphic allozyme loci were studied in 569 individuals collected across the Iberian Peninsula and North Africa. The obtained data were analysed using both conventional population genetic tools and recent Bayesian model-based clustering methods. Our results show that there are several well-differentiated entities corroborating the major splits observed in mtDNA analyses. These groups correspond not only to the fully recognized species Podarcis bocagei , Podarcis carbonelli , and Podarcis vaucheri but also to multiple forms within the polytypic Podarcis hispanica , all of which have a similar level of differentiation to that observed between the acknowledged species. However, relationships between forms are weakly supported both by population and individual clustering methods, suggesting a scenario of a rapid diversification that contrasts to the clear bifurcating model assumed from previous mtDNA analyses. Individual multilocus analyses report few individuals misassigned or apparently admixed, some of which are most likely explained by the persistence of high levels of ancestral polymorphism. Other admixed individuals, however, are probably the result of limited levels of gene flow between forms.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 121–133.     

Genetic admixture between the Iberian endemic lizards Podarcis bocagei and Podarcis carbonelli: evidence for limited natural hybridization and a bimodal hybrid zone

When recently diverged taxa come into contact, the extent of introgression between them is related to the degree of differentiation that they have achieved. Studying contact zones is therefore essential to understand if differentiated taxa are reproductively isolated and, ultimately, if they are likely to remain distinct. Recent work on Iberian and North African wall lizards ( Podarcis ) has documented the existence of multiple evolutionary units, diagnosable both by genetic markers and morphology, but suggests that gene flow between distinct forms has occurred. Therefore, we were interested in evaluating how species boundaries are maintained in the areas where they meet. In this work, we study the contact zone between Podarcis bocagei and Podarcis carbonelli . We sampled a transect including the only locality where these two lizards are known to occur in syntopy and analysed a battery of 15 unlinked nuclear genetic markers and mitochondrial DNA. We also conducted a preliminary analysis of morphology and fertility. Using model-based clustering approaches, we show that the two species hybridize in the population where they have direct contact, but evidences of introgression are low for nearby populations. Although a significant number of individuals show evidence of admixture, this hybrid zone is clearly bimodal, suggesting strong barriers to gene flow, of which the putative nature are discussed. Interestingly, morphological analyses do not support the existence of intermediate forms among individuals that are admixed genetically. Taken together, these results constitute further evidence validating P. bocagei and P. carbonelli as distinct species.     

Mitochondrial DNA sequence variation in Portuguese native dog breeds: diversity and phylogenetic affinities

In an extensive survey of the genetic diversity in Portuguese dogs, we have examined an 887-bp fragment of the mitochondrial DNA (mtDNA) from 8 Portuguese, 1 Spanish, and 2 North African native dog breeds, including village dogs from Portugal and Tunisia. Forty-nine haplotypes were found in the 164 individuals analyzed, with private haplotypes being found in several breeds. For example, the Castro Laboreiro Watchdog, a rare breed from a small and isolated region in Portugal, was monomorphic for mtDNA and possessed a new haplotype, which may be provisionally considered a breed-specific marker. Phylogenetic analyses recapitulated 4 major clades identified in other studies, but new haplotypes, grouping within a clade that was previously thought as geographically restricted, were detected in Portugal and Morocco. Portuguese village dogs showed no genetic differentiation from nonnative dogs or from local breeds of the areas in which the village dogs were sampled. Although Iberian and North African dog breeds possessed breed-specific mtDNA haplotypes, no significant geographic structure could be detected among them. There is no evidence for introgression of North African haplotypes in Iberian dogs, contrary to previous results for other domestic animals.     

African female heritage in Iberia: a reassessment of mtDNA lineage distribution in present times

The Iberian peninsula is a peripheral region of Europe in close proximity to Africa. Its inhabitants have an overall mtDNA genetic landscape typical of European background, although with signs of some African influence, whose features we deemed to disclose by analyzing available mtDNA HVRI distributions and new data. We analyzed 1,045 sequences. The most relevant results are the following: (1) North African sequences (haplogroup U6) present an overall frequency of 2.39%, and sub-Saharan sequences reach 3.83%, values that are, in both cases, much higher than those generally observed in Europe; and (2) there is a substantial geographic heterogeneity in the distribution of these lineages (haplogroup L being the most frequent in the south, whereas haplogroup U6 is generally more common in the north). The analysis of the observed diversity within each haplogroup strongly suggests that both were recently introduced (in historical times). Although for haplogroup U6 the documented event that is demographically compatible is the Islamic period (beginning of the 8th century to the end of the 15th century), for haplogroup L the most probable origin is the modern slave trade (mid 15th century to the end of the 18th century). However, the observed geographic structuring for one of the haplogroups does not fit the expected distribution provided by simplistic historical considerations. In fact, although for haplogroup L the north-south increasing frequency is corroborated by historical data, the opposite trend, observed for haplogroup U6, is more difficult to reconcile with the magnitude and time span of the Islamic political and cultural influence, which lasted longer and was more intense in the south. To clarify this conundrum, we need not only a substantial increase in the amount of mtDNA data (particularly for North Africa) but also new historical data and interpretations.     

Early Holocenic and Historic mtDNA African Signatures in the Iberian Peninsula: The Andalusian Region as a Paradigm

Determining the timing, identity and direction of migrations in the Mediterranean Basin, the role of “migratory routes” in and among regions of Africa, Europe and Asia, and the effects of sex-specific behaviors of population movements have important implications for our understanding of the present human genetic diversity. A crucial component of the Mediterranean world is its westernmost region. Clear features of transcontinental ancient contacts between North African and Iberian populations surrounding the maritime region of Gibraltar Strait have been identified from archeological data. The attempt to discern origin and dates of migration between close geographically related regions has been a challenge in the field of uniparental-based population genetics. Mitochondrial DNA (mtDNA) studies have been focused on surveying the H1, H3 and V lineages when trying to ascertain north-south migrations, and U6 and L in the opposite direction, assuming that those lineages are good proxies for the ancestry of each side of the Mediterranean. To this end, in the present work we have screened entire mtDNA sequences belonging to U6, M1 and L haplogroups in Andalusians—from Huelva and Granada provinces—and Moroccan Berbers. We present here pioneer data and interpretations on the role of NW Africa and the Iberian Peninsula regarding the time of origin, number of founders and expansion directions of these specific markers. The estimated entrance of the North African U6 lineages into Iberia at 10 ky correlates well with other L African clades, indicating that U6 and some L lineages moved together from Africa to Iberia in the Early Holocene. Still, founder analysis highlights that the high sharing of lineages between North Africa and Iberia results from a complex process continued through time, impairing simplistic interpretations. In particular, our work supports the existence of an ancient, frequently denied, bridge connecting the Maghreb and Andalusia.     

Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East‐Caucasian and Central Italian populations

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East‐Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.     

Counting the founders: the matrilineal genetic ancestry of the Jewish Diaspora

The history of the Jewish Diaspora dates back to the Assyrian and Babylonian conquests in the Levant, followed by complex demographic and migratory trajectories over the ensuing millennia which pose a serious challenge to unraveling population genetic patterns. Here we ask whether phylogenetic analysis, based on highly resolved mitochondrial DNA (mtDNA) phylogenies can discern among maternal ancestries of the Diaspora. Accordingly, 1,142 samples from 14 different non-Ashkenazi Jewish communities were analyzed. A list of complete mtDNA sequences was established for all variants present at high frequency in the communities studied, along with high-resolution genotyping of all samples. Unlike the previously reported pattern observed among Ashkenazi Jews, the numerically major portion of the non-Ashkenazi Jews, currently estimated at 5 million people and comprised of the Moroccan, Iraqi, Iranian and Iberian Exile Jewish communities showed no evidence for a narrow founder effect, which did however characterize the smaller and more remote Belmonte, Indian and the two Caucasus communities. The Indian and Ethiopian Jewish sample sets suggested local female introgression, while mtDNAs in all other communities studied belong to a well-characterized West Eurasian pool of maternal lineages. Absence of sub-Saharan African mtDNA lineages among the North African Jewish communities suggests negligible or low level of admixture with females of the host populations among whom the African haplogroup (Hg) L0-L3 sub-clades variants are common. In contrast, the North African and Iberian Exile Jewish communities show influence of putative Iberian admixture as documented by mtDNA Hg HV0 variants. These findings highlight striking differences in the demographic history of the widespread Jewish Diaspora.     

A European Melting Pot of Harbour Porpoise in the French Atlantic Coasts Inferred from Mitochondrial and Nuclear Data

Field surveys have reported a global shift in harbour porpoise distribution in European waters during the last 15 years, including a return to the Atlantic coasts of France. In this study, we analyzed genetic polymorphisms at a fragment of the mitochondrial control region (mtDNA CR) and 7 nuclear microsatellite loci, for 52 animals stranded and by-caught between 2000 and 2010 along the Atlantic coasts of France. The analysis of nuclear and mitochondrial loci provided contrasting results. The mtDNA revealed two genetically distinct groups, one closely related to the Iberian and African harbour porpoises, and the second related to individuals from the more northern waters of Europe. In contrast, nuclear polymorphisms did not display such a distinction. Nuclear markers suggested that harbour porpoises behaved as a randomly mating population along the Atlantic coasts of France. The difference between the two kinds of markers can be explained by differences in their mode of inheritance, the mtDNA being maternally inherited in contrast to nuclear loci that are bi-parentally inherited. Our results provide evidence that a major proportion of the animals we sampled are admixed individuals from the two genetically distinct populations previously identified along the Iberian coasts and in the North East Atlantic. The French Atlantic coasts are clearly the place where these two previously separated populations of harbour porpoises are now admixing. The present shifts in distribution of harbour porpoises along this coast is likely caused by habitat changes that will need to be further studied.     

Complex biogeographical distribution of genetic variation within Podarcis wall lizards across the Strait of Gibraltar

Abstract Aim To examine the effect of a known geological barrier on genetic variation within a wall lizard species complex. Location The Iberian Peninsula and North Africa. Methods Sequencing of partial 12S rRNA and cytochrome b mtDNA. Results The current distribution of genetic variability is not related to the opening of the Strait of Gibraltar. Conclusions Podarcis hispanica in North Africa is probably a species complex. The Strait of Gibraltar should not be used as a known barrier to gene flow in other land based organisms without careful sampling to test for multiple crossings since its formation.     

Evolutionary patterns of the mitochondrial genome in the Moorish gecko, Tarentola mauritanica

A previous study on the evolutionary patterns of Tarentola mauritanica demonstrated that low levels of mitochondrial diversity observed in the European populations relative to nuclear markers were consistent with a selective sweep hypothesis. In order to unravel the mitochondrial evolutionary history in this European population and two other lineages of T. mauritanica (Iberian and North African clades), variation within 22 nearly complete mitogenomes was analyzed. Surprisingly, each clade seems to have a distinct evolutionary history; with both the European and Iberian clades presenting a decrease of polymorphism, which in the former is consistent with departure from neutrality of the mtDNA (positive or background selection), but in the latter seems to be the result of a bottleneck after a population expansion. The pattern exhibited by the North African clade seems to be a consequence of adaptation to certain mtDNA variants by positive selection.     

THE ORIGIN OF WEST EUROPEAN SUBSPECIES OF HONEYBEES (APIS MELLIFERA): NEW INSIGHTS FROM MICROSATELLITE AND MITOCHONDRIAL DATA

Apis mellifera is composed of three evolutionary branches including mainly African (branch A), western and northern European (branch M), and southeastern European (branch C) populations. The existence of morphological clines extending from the equator to the Polar Circle through Morocco and Spain raised the hypothesis that the branch M originated in Africa. Mitochondrial DNA analysis revealed that branches A and M were characterized by highly diverged lineages implying very remote links between both branches. It also revealed that mtDNA haplotypes from lineages A coexisted with haplotypes M in the Iberian Peninsula and formed a south-north frequency cline, suggesting that this area could be a secondary contact zone between the two branches. By analyzing 11 populations sampled along a France-Spain/Portugal-Morocco-Guinea transect at 8 microsatellite loci and the DraI RFLP of the COI-COII mtDNA marker, we show that Iberian populations do not present any trace of “africanization” and are very similar to French populations when considering microsatellite markers. Therefore, the Iberian Peninsula is not a transition area. The higher haplotype A variability observed in Spanish and Portuguese samples compared to that found in Africa is explained by a higher mutation rate and multiple and recent introductions. Selection appears to be the best explanation to the morphological and allozymic clines and to the diffusion and maintenance of African haplotypes in Spain and Portugal.     

Phylogeography of the Ibero-Maghrebian red-eyed grass snake (<Emphasis Type="Italic">Natrix astreptophora</Emphasis>)

We examined phylogeographic differentiation of the red-eyed grass snake (Natrix astreptophora) using 1984 bp of mtDNA and 13 microsatellite loci from specimens collected across its distribution range in southwestern Europe and northwestern Africa. Based on phylogenetic analyses of mtDNA, European N. astreptophora constituted the sister clade to a weakly supported North African clade comprised of two deeply divergent and well-supported clades, one corresponding to Moroccan snakes and the other to snakes from Algeria and Tunisia. This tripartite differentiation was confirmed by analyses of microsatellite loci. According to a fossil-calibrated molecular clock, European and North African N. astreptophora diverged 5.44 million years ago (mya), and the two Maghrebian clades split 4.64 mya. These dates suggest that the radiation of the three clades was initiated by the environmental changes related to the Messinian Salinity Crisis and the reflooding of the Mediterranean Basin. The differentiation of N. astreptophora, with distinct clades in the Iberian Peninsula and in the western and eastern Maghreb, corresponds to a general phylogeographic paradigm and resembles the differentiation found in another co-distributed Natrix species, the viperine snake (N. maura). Despite both species being good swimmers, the Strait of Gibraltar constitutes a significant biogeographic barrier for them. The discovery that North Africa harbours two endemic lineages of N. astreptophora necessitates more conservation efforts for these imperilled snakes.     

Genetic characterization of southwestern European bovine breeds: a historical and biogeographical reassessment with a set of 16 microsatellites

The origin of Iberian cattle has been suggested by some authors to be the product of European and north African cattle entrances during the last few thousands of years. However, these hypotheses were mainly based on morphological similarities. This study analyzed 889 unrelated individuals from 15 representative Iberian breeds and 3 French breeds for 16 microsatellite loci. Statistical tests were used to calculate interpopulation genetic distances (D(A)) and principal components analysis (PCA). To visualize the geographical distribution of the genetic differentiation between Iberian cattle breeds, data from the PCA analysis were used to construct synthetic maps. Genetic similarity among neighboring Iberian breeds is mainly caused by gene flow. However, recent demographic fluctuations and reproductive isolation in Alistana, Mirandesa, and Tudanca has increased genetic drift, which may be the main cause for the relatively high differentiation of these populations. The synthetic maps constructed with the first and second PCs revealed (1) a large differentiation between Northern Iberian breeds rather than between more geographically distant breeds, and (2) a clear east-west gradient that may be related with the model of demic diffusion of agriculture. Finally, we detected no strong evidence for an African genetic influence in the Iberian cattle breeds analyzed in this study.     

Evidence of a cenozoic Betic-Kabilian connection based on freshwater fish phylogeography (Luciobarbus, Cyprinidae)

The phylogenetic relationships among 26 species of the subgenus Luciobarbus were examined through comparison of the complete sequence of the mitochondrial genes ATPase 6 and 8 and cytochrome b. The monophyletic condition of the Luciobarbus subgenus was confirmed by several data treatment methods. Findings indicate a closer relationship among species inhabiting Caucasian, Greek, and North African areas than between the latter and those of the Iberian region. The isolation of these main clades is postulated to have occurred after the Messinian salinity crisis (5.5 MY ago) when the Iberian Peninsula broke away from the African continent. Based on a calibrated molecular clock proposed herein, the subsequent splitting of this subgenus corresponds to the isolation of the other three areas (Caucasian, Greek, and North African) about 4.5 MY ago, which would have interrupted any possible gene flow between lineages. However, the present data indicate that the North African Kabilies Mountain area maintained its contact with the south Iberian Peninsula for a longer period.     

Karyosystematic study of Potentilla L. subgen. Potentilla (Rosaceae) in the Iberian Peninsula

Chromosome numbers in 80 populations belonging to 18 species of Potentilla L. subgen. Potentilla from the Iberian Peninsula and two of P. maura, a North African endemic taxon, have been counted. The basic number of chromosomes is always x = 7 and these chromosomes are small (between 1 and 2 μm). For three species, the number of chromosomes is reported for the first time and, for another six, this number has been established in Iberian representatives. Moreover, new ploidy levels have been obtained for P. hispanica and P. crantzii with regard to their entire distribution area, and in P. cinerea and P. neumanniana for the Iberian Peninsula. Some taxonomic, phylogenetic and phytogeographic comments are made for several species or groups of species from the West Mediterranean region. In 13 species only one ploidy level has been found, but six species have several ploidy levels. Seven ploidy levels occur in the investigated taxa. The frequency of each ploidy level represented within Iberian Potentilla is analysed and the data are compared with those available for taxa from the rest of the distribution area of the genus.     

Genetic and morphological relationships of the Berlenga wall lizard (Podarcis bocagei berlengensis: Lacertidae)

Six biometric characters and 15 electrophoretic loci of three known subspecies of Podarcis bocagei were studied. Contrary to biometric indications, Podarcis bocagei berlengensis showed a closer genetic relationship with Podarcis bocagei carbonelli than with the nominal subspecies. The biometric results were confounded by the relative large size of P. b. berlengensis . Post-glacial colonization from two different areas might explain the greater genetic differentiation between the P. b. bocagei and P. b. carbonelli populations (Nei's D - 0.158 - near the species threshold). It might also explain the known colonization of these two lizards through the islands off the west coast of the Iberian penninsula and their present distributions to the north of (Galicia; Northern Portugal) and to the south of the valley of river Douro (central Portugal, western and Sistema Central), respectively.     

Mitochondrial haplogroup H1 in north Africa: an early holocene arrival from Iberia

The Tuareg of the Fezzan region (Libya) are characterized by an extremely high frequency (61%) of haplogroup H1, a mitochondrial DNA (mtDNA) haplogroup that is common in all Western European populations. To define how and when H1 spread from Europe to North Africa up to the Central Sahara, in Fezzan, we investigated the complete mitochondrial genomes of eleven Libyan Tuareg belonging to H1. Coalescence time estimates suggest an arrival of the European H1 mtDNAs at about 8,000-9,000 years ago, while phylogenetic analyses reveal three novel H1 branches, termed H1v, H1w and H1x, which appear to be specific for North African populations, but whose frequencies can be extremely different even in relatively close Tuareg villages. Overall, these findings support the scenario of an arrival of haplogroup H1 in North Africa from Iberia at the beginning of the Holocene, as a consequence of the improvement in climate conditions after the Younger Dryas cold snap, followed by in situ formation of local H1 sub-haplogroups. This process of autochthonous differentiation continues in the Libyan Tuareg who, probably due to isolation and recent founder events, are characterized by village-specific maternal mtDNA lineages.     

Ancestral matrilineages and mitochondrial DNA diversity of the Lidia cattle breed

To clarify the genetic ancestry and the mitochondrial DNA (mtDNA) diversity of the Lidia cattle breed, a 521-bp D-loop fragment was sequenced in 527 animals belonging to 70 herds distributed across 29 lineages. The mtDNA diversity recorded was similar to that seen for Middle Eastern breeds and greater than that recorded for the majority of European breeds. Haplotype T3 was the most common (81%), followed by the African T1 haplotype (17%); very low frequencies were recorded for haplotypes T and T2. The results agree with there being two major ancestral lines for the Lidia breed, European and African, similar to that seen for other Mediterranean breeds. A wide range of variation in haplotype frequencies was seen between the examined lineages. Haplotype T3 was present in all those analysed; in five it was the only one present, and in only one lineage (Miura) was its frequency lower than that of T1. T1*, a haplotype reported in Criollo breeds and to date in only a single European breed (the Retinta breed from Spain), was found in a single animal belonging to the Concha y Sierra lineage. Network analysis of the Lidia breed revealed the presence of two major haplotypes: T3 and T1. The Lidia breed appears to be more closely related to prehistoric Iberian and Italian than to British aurochs.     

Alu insertion polymorphisms in NW Africa and the Iberian Peninsula: evidence for a strong genetic boundary through the Gibraltar Straits

An analysis of 11 I Alu insertion polymorphisms (ACE, TPA25, PV92, APO, FXIIIB, D1, A25, B65, HS2.43, HS3.23, and HS4.65) has been performed in several NW African (Northern, Western, and Southeastern Moroccans, Saharawi; Algerians; Tunisians) and Iberian (Basques, Catalans, and Andalusians) populations. Genetic distances and principal component analyses show a clear differentiation of NW African and Iberian groups of samples, suggesting a strong genetic barrier matching the geographical Mediterranean Sea barrier. The restriction to gene flow may be attributed to the navigational hazards across the Straits, but cultural factors must also have played a role. Some degree of gene flow from sub-Saharan Africa can be detected in the southern part of North Africa and in Saharawi and Southeastern Moroccans, as a result of a continuous gene flow across the Sahara desert that has created a south-north cline of sub-Saharan Africa influence in North Africa. Iberian samples show a substantial degree of homogeneity and fall within the cluster of European-based genetic diversity.     

Multilocus and morphological analysis of south-eastern Iberian Wall lizards (Squamata,Podarcis)

The phylogenetic relationships among the wall lizards of the Podarcis hispanicus complex that inhabit the south-east (SE) of the Iberian Peninsula and other lineages of the complex remain unclear. In this study, four mitochondrial and two nuclear markers were used to study genetic relationships within this complex. The phylogenetic analyses based on mtDNA gene trees constructed with ML and BI, and a species tree using *BEAST support three divergent clades in this region: the Valencia, Galera and Albacete/Murcia lineages. These three lineages were also corroborated in species delimitation analyses based on mtDNA using bPTP, mPTP, GMYC, ABGD and BAPS. Bayesian inference species delimitation method (BPP) based on both nuclear data and a combined data set (mtDNA + nuclear) showed high posterior probabilities for these three SE lineages (≥0.94) and another Bayesian analysis (STACEY) based on combined data set recovered the same three groups in this region. Divergence time dating of the species tree provided an estimated divergence of the Galera lineage from the other SE group (Podarcis vaucheri, (Albacete/Murcia, Valencia)) at 12.48 Ma. During this period, the Betic–Rifian arc was isolated, which could have caused the isolation of the Galera form distributed to the south of the Betic Corridor. Although lizards from the Albacete/Murcia and Galera lineage are morphologically similar, they clearly represent distinct genetic lineages. The noteworthy separation of the Galera lineage enables us to conclude that this lineage must be considered as a new full species.