Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.     

Adaptive divergence between lake and stream populations of an East African cichlid fish

Divergent natural selection acting in different habitats may build up barriers to gene flow and initiate speciation. This speciation continuum can range from weak or no divergence to strong genetic differentiation between populations. Here, we focus on the early phases of adaptive divergence in the East African cichlid fish Astatotilapia burtoni, which occurs in both Lake Tanganyika (LT) and inflowing rivers. We first assessed the population structure and morphological differences in A. burtoni from southern LT. We then focused on four lake–stream systems and quantified body shape, ecologically relevant traits (gill raker and lower pharyngeal jaw) as well as stomach contents. Our study revealed the presence of several divergent lake–stream populations that rest at different stages of the speciation continuum, but show the same morphological and ecological trajectories along the lake–stream gradient. Lake fish have higher bodies, a more superior mouth position, longer gill rakers and more slender pharyngeal jaws, and they show a plant/algae and zooplankton‐biased diet, whereas stream fish feed more on snails, insects and plant seeds. A test for reproductive isolation between closely related lake and stream populations did not detect population‐assortative mating. Analyses of F1 offspring reared under common garden conditions indicate that the detected differences in body shape and gill raker length do not constitute pure plastic responses to different environmental conditions, but also have a genetic basis. Taken together, the A. burtoni lake–stream system constitutes a new model to study the factors that enhance and constrain progress towards speciation in cichlid fishes.     

Phylogeography of northern Dolly Varden Salvelinus malma malma based on analysis of mitochondrial DNA

The northern Dolly Varden, Salvelinus malma malma, is a typical representative of arctic fauna distributed in northeastern Asia and northwestern North America. Because its spawning habitats were affected by Pleistocene glacial advances over most of its natural range, S. m. malma is among the most interesting objects of phylogeographic and microevolutionary studies. We reconstructed the genealogy of mtDNA haplotypes from 27 Alaskan and Asian populations to study the influence of glacial and geological vicariance events on the contemporary population genetic structure, phylogeographic subdivision and distribution of the northern Dolly Varden. Analysis of restriction site states in three PCR‐amplified mtDNA regions (ND1/ND2, ND5/ND6, Cytb/D‐loop; 47% of the mitochondrial genome) resolved 75 haplotypes in 436 fish. Similar patterns of subspecific variation apparent from hierarchical diversity and nested clade analyses of mtDNA haplotypes identify weak spatial differentiation and low levels of divergence. Our results suggest that (1) demographic history has been influenced by historical range expansions and recent isolation by distance, (2) present populations from Asia and North America were colonized from one main Beringian Refugium, and (3) that this taxon's ancestral population probably experienced a bottleneck in the Beringian Refugium during the late Pleistocene (Wisconsin) glacial period.     

Contrasting geographic structure in evolutionarily divergent Lake Tanganyika catfishes

Geographic isolation is suggested to be among the most important processes in the generation of cichlid fish diversity in East Africa's Great Lakes, both through isolation by distance and fluctuating connectivity caused by changing lake levels. However, even broad scale phylogeographic patterns are currently unknown in many non‐cichlid littoral taxa from these systems. To begin to address this, we generated restriction‐site‐associated DNA sequence (RADseq) data to investigate phylogeographic structure throughout Lake Tanganyika (LT) in two broadly sympatric rocky shore catfish species from independent evolutionary radiations with differing behaviors: the mouthbrooding claroteine, Lophiobagrus cyclurus, and the brood‐parasite mochokid, Synodontis multipunctatus. Our results indicated contrasting patterns between these species, with strong lake‐wide phylogeographic signal in L. cyclurus including a deep divergence between the northern and southern lake basins. Further structuring of these populations was observed across a heterogeneous habitat over much smaller distances. Strong population growth was observed in L. cyclurus sampled from shallow shorelines, suggesting population growth associated with the colonization of new habitats following lake‐level rises. Conversely, S. multipunctatus, which occupies a broader depth range, showed little phylogeographic structure and lower rates of population growth. Our findings suggest that isolation by distance and/or habitat barriers may play a role in the divergence of non‐cichlid fishes in LT, but this effect varies by species.     

Gene flow and Andean uplift shape the diversification of Gasteracantha cancriformis (Araneae: Araneidae) in Northern South America

The Andean uplift has played a major role in shaping the current Neotropical biodiversity. However, in arthropods other than butterflies, little is known about how this geographic barrier has impacted species historical diversification. Here, we examined the phylogeography of the widespread color polymorphic spider Gasteracantha cancriformis to evaluate the effect of the northern Andean uplift on its divergence and assess whether its diversification occurred in the presence of gene flow. We inferred phylogenetic relationships and divergence times in G. cancriformis using mitochondrial and nuclear data from 105 individuals in northern South America. Genetic diversity, divergence, and population structure were quantified. We also compared multiple demographic scenarios for this species using a model‐based approach (Phrapl ) to determine divergence with or without gene flow. At last, we evaluated the association between genetic variation and color polymorphism. Both nuclear and mitochondrial data supported two well‐differentiated clades, which correspond to populations occurring on opposite sides of the Eastern cordillera of the Colombian Andes. The final uplift of this cordillera was identified as the most likely force that shaped the diversification of G. cancriformis in northern South America, resulting in a cis‐ and trans‐Andean phylogeographic structure for the species. We also found shared genetic variation between the cis‐ and trans‐Andean clades, which is better explained by a scenario of historical divergence in the face of gene flow. This has been likely facilitated by the presence of low‐elevation passes across the Eastern Colombian cordillera. Our work constitutes the first example in which the Andean uplift coupled with gene flow influenced the evolutionary history of an arachnid lineage.     

Phylogeographic signal variation in mitochondrial DNA among geographically isolated grassland butterflies

Aim We surveyed mitochondrial DNA (mtDNA) sequence variation among regionally isolated populations of 10 grassland‐associated butterfly species to determine: (1) the utility of phylogeographic comparisons among multiple species for assessing recent evolutionary patterns, and (2) the respective roles of isolation attributable to range disjunction versus isolation attributable solely to geographic distance in establishing divergence patterns. Location The Peace River grasslands of northern Alberta and British Columbia, Canada, which are isolated by 300+ km from similar communities to the south. Methods We sequenced mtDNA (1420 bp of cytochrome c oxidase subunit I) from five grassland‐restricted butterfly species that have geographically disjunct populations and from five ecologically broader species that have more continuous distributions across the same regions. Using analysis of molecular variance (AMOVA), Mantel and partial Mantel tests, and haplotype networks, we compared population structure within and between species in order to assess the validity of single‐species phylogeographic characterizations. We then contrasted variance components between disjunct and continuously distributed species to assess whether divergences were correlated more with disjunction or with geographic distance. Results Single‐species analyses varied substantially within both the disjunct and the continuous groups. One species in each of these groups had mtDNA with unusually deep intraspecific mitochondrial lineage divergences. On the whole, however, the five species with disjunct ranges exhibited greater divergence between geographically distant populations than did the five species with continuous distributions. Comparison of variance components between disjunct and continuous species indicated that isolation attributable only to geographic distance was responsible for up to half of the total sequence variation between disjunct populations of grassland butterflies. Main conclusions Our findings show that single‐species phylogeographic analyses of post‐Pleistocene butterfly distributions are inadequate for characterizing regional biogeographic divergence histories. However, comparison of mtDNA sequence divergences between groups of disjunct and continuously distributed species can allow isolation attributable to range interruption to be quantitatively distinguished from isolation attributable solely to gene flow attenuation over the same geographic area.     

Influence of invasion history on rapid morphological divergence across island populations of an exotic bird

There is increasing evidence that exotic populations may rapidly differentiate from those in their native range and that differences also arise among populations within the exotic range. Using morphological and DNA‐based analyses, we document the extent of trait divergence among native North American and exotic Hawaiian populations of northern cardinal (Cardinalis cardinalis). Furthermore, using a combination of historical records and DNA‐based analyses, we evaluate the role of founder effects in producing observed trait differences. We measured and compared key morphological traits across northern cardinal populations in the native and exotic ranges to assess whether trait divergence across the Hawaiian Islands, where this species was introduced between 1929 and 1931, reflected observed variation across native phylogeographic clades in its native North America. We used and added to prior phylogenetic analyses based on a mitochondrial locus to identify the most likely native source clade(s) for the Hawaiian cardinal populations. We then used Approximate Bayesian Computation (ABC) to evaluate the role of founder effects in producing the observed differences in body size and bill morphology across native and exotic populations. We found cardinal populations on the Hawaiian Islands had morphological traits that diverged substantially across islands and overlapped the trait space of all measured native North American clades. The phylogeographic analysis identified the eastern North American clade (C. cardinalis cardinalis) as the most likely and sole native source for all the Hawaiian cardinal populations. The ABC analyses supported written accounts of the cardinal's introduction that indicate the original 300 cardinals shipped to Hawaii were simultaneously and evenly released across Hawaii, Kauai, and Oahu. Populations on each island likely experienced bottlenecks followed by expansion, with cardinals from the island of Hawaii eventually colonizing Maui unaided. Overall, our results suggest that founder effects had limited impact on morphological trait divergence of exotic cardinal populations in the Hawaiian archipelago, which instead reflect postintroduction events.     

Genetic divergence of mitochondrial DNA in white char <Emphasis Type="Italic">Salvelinus albus</Emphasis> and northern Dolly Varden char <Emphasis Type="Italic">Salvelinus malma malma</Emphasis>

Comparative analysis of mitochondrial DNA variation was performed in white char Salvelinus albus Glubokovsky and in its putative ancestor species, northern Dolly Varden char Salvelinus malma malma Walbaum. Highly statistically significant differentiation of S. albus and S. m. malma in the areas of sympatric (Kamchatka River basin) and allopatric (Kronotskoe Lake and Kronotskaya River) residence was demonstrated. The mtDNA divergence between S. albus and S. m. malma did not exceed the range of intraspecific variation in the populations of northern Dolly Varden char. At the same time, clusterization pattern of the Salvelinus chars provides hypothesis on the common origin of two allopatric populations of white char. Genealogical analysis of haplotypes indicates that S. albus and S. m. malma currently demonstrate incomplete divergence of mitochondrial lineages. The low nucleotide divergence estimates between S. albus and S. m. malma reflect the short time period since the beginning of the divergence of ancestral lineages. These estimates are determined by ancestral polymorphism and haplotype exchange between the diverged phylogenetic groups as a result of introgressive hybridization.     

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high‐throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P. streckeri/P. illinoensis populations, but suggest the disjunct distribution of P. illinoensis occurred more recently. Our best‐supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P. illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P. illinoensis and demonstrates the utility of model‐based approaches for phylogeographic inference.     

Complete mitochondrial DNA replacement in a Lake Tanganyika cichlid fish

We used nuclear and mitochondrial DNA (mtDNA) sequences from specimens collected throughout Lake Tanganyika to clarify the evolutionary relationship between Lamprologus callipterus and Neolamprologus fasciatus . The nuclear data support the reciprocal monophyly of these two shell-breeding lamprologine cichlids. However, mtDNA sequences show that (i) L. callipterus includes two divergent and geographically disjunct (North–South) mtDNA lineages; and that (ii) N. fasciatus individuals cluster in a lineage sister group to the northern lineage of L. callipterus . The two mtDNA lineages of L. callipterus diverged c . 684 kya to 1.2 Ma, coinciding with a major water level low stand in Lake Tanganyika, which divided the lake into isolated sub-lakes. This suggests that the two mtDNA lineages originated as the result of the separation of L. callipterus populations in different sub-basins. The incongruent phylogenetic position of N. fasciatus can best be explained by an ancient unidirectional introgression from L. callipterus into N. fasciatus. Remarkably, our data indicate that this event resulted in the complete mtDNA replacement in N. fasciatus . Our data suggest that hybridization occurred soon after the divergence of the two L. callipterus mtDNA lineages, probably still during the water level low stand, and that subsequently the invading mtDNA lineage spread throughout the lake.     

Comparative multilocus phylogeography of two Palaearctic spruce bark beetles: influence of contrasting ecological strategies on genetic variation

While phylogeographic patterns of organisms are often interpreted through past environmental disturbances, mediated by climate changes, and geographic barriers, they may also be strongly influenced by species‐specific traits. To investigate the impact of such traits, we focused on two Eurasian spruce bark beetles that share a similar geographic distribution, but differ in their ecology and reproduction. Ips typographus is an aggressive tree‐killing species characterized by strong dispersal, whereas Dendroctonus micans is a discrete inbreeding species (sib mating is the rule), parasite of living trees and a poor disperser. We compared genetic variation between the two species over both beetles’ entire range in Eurasia with five independent gene fragments, to evaluate whether their intrinsic differences could have an influence over their phylogeographic patterns. We highlighted widely divergent patterns of genetic variation for the two species and argue that the difference is indeed largely compatible with their contrasting dispersal strategies and modes of reproduction. In addition, genetic structure in I. typographus divides European populations in a northern and a southern group, as was previously observed for its host plant, and suggests past allopatric divergence. A long divergence time was estimated between East Asian and other populations of both species, indicating their long‐standing presence in Eurasia, prior to the last glacial maximum. Finally, the strong population structure observed in D. micans for the mitochondrial locus provides insights into the recent colonization history of this species, from its native European range to regions where it was recently introduced.     

Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus <Emphasis Type="Italic">Tropheus</Emphasis>

Background  

Cichlid fishes are notorious for their wealth of intra- and interspecific colour pattern diversity. In Lake Tanganyika, the endemic genus Tropheus represents the most impressive example for geographic variation in the pattern and hue of integument colouration, but the taxonomy of the over 100 mostly allopatric colour morphs remains to a large degree unresolved. Previous studies of mitochondrial DNA sequence data revealed polyphyly of the six nominally described species and complex phylogeographic patterns influenced by lake level fluctuations and population admixture, and suggested the parallel evolution of similar colour patterns in divergent evolutionary lineages. A gene tree of a rapidly radiating group may be subject to incomplete and stochastic lineage sorting, and to overcome this problem we used multi-locus, nuclear AFLP data in comparison with mtDNA sequences to study diversification, migration and introgression in Tropheus colour morphs in Lake Tanganyika.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Historic speciation and recent colonization of Eurasian monkey gobies (Neogobius fluviatilis and N. pallasi) revealed by DNA sequences,microsatellites, and morphology

Aim Hidden diversity within an invasive ‘species’ can mask both invasion pathways and confound management goals. We assessed taxonomic status and population structure of the monkey goby Neogobius fluviatilis across Eurasia, comparing genetic variation across its native and invasive ranges. Location Native populations were analysed within the Black and Caspian Sea basins, including major river drainages (Dnieper, Dniester, Danube, Don and Volga rivers), along with introduced locations within the upper Danube and Vistula river systems. Methods DNA sequences and 10 nuclear microsatellite loci were analysed to test genetic diversity and divergence patterns of native and introduced populations; phylogenetic analysis of mtDNA cytochrome b and nuclear RAG‐1 sequences assessed taxonomic status of Black and Caspian Sea lineages. Multivariate analysis of morphology was used to corroborate phylogenetic patterns. Population genetic structure within each basin was evaluated with mtDNA and microsatellite data using F_ST analogues and Bayesian assignment tests. Results Phylogenetic analysis of mitochondrial and nuclear sequences discerned a pronounced genetic break between monkey gobies in the Black and Caspian Seas, indicating a long‐term species‐level separation dating to c. 3 million years. This pronounced separation further was confirmed from morphological and population genetic divergence. Bayesian inference showed congruent patterns of population structure within the Black Sea basin. Introduced populations in the Danube and Vistula River basins traced to north‐west Black Sea origins, a genetic expansion pattern matching that of other introduced Ponto‐Caspian gobiids. Main conclusions Both genetic and morphological data strongly supported two species of monkey gobies that were formerly identified as subspecies: N. fluviatilis in the Black Sea basin, Don and Volga Rivers, and the Kumo‐Manych Depression, and Neogobius pallasi in the Caspian Sea and Volga River delta. Genetic origins of introduced N. fluviatilis populations indicated a common invasion pathway shared with other introduced Ponto‐Caspian fishes and invertebrates.     

Adaptive phenotypic plasticity contributes to divergence between lake and river populations of an East African cichlid fish

Adaptive phenotypic plasticity and fixed genotypic differences have long been considered opposing strategies in adaptation. More recently, these mechanisms have been proposed to act complementarily and under certain conditions jointly facilitate evolution, speciation, and even adaptive radiations. Here, we investigate the relative contributions of adaptive phenotypic plasticity vs. local adaptation to fitness, using an emerging model system to study early phases of adaptive divergence, the generalist cichlid fish species Astatotilapia burtoni. We tested direct fitness consequences of morphological divergence between lake and river populations in nature by performing two transplant experiments in Lake Tanganyika. In the first experiment, we used wild‐caught juvenile lake and river individuals, while in the second experiment, we used F1 crosses between lake and river fish bred in a common garden setup. By tracking the survival and growth of translocated individuals in enclosures in the lake over several weeks, we revealed local adaptation evidenced by faster growth of the wild‐caught resident population in the first experiment. On the other hand, we did not find difference in growth between different types of F1 crosses in the second experiment, suggesting a substantial contribution of adaptive phenotypic plasticity to increased immigrant fitness. Our findings highlight the value of formally comparing fitness of wild‐caught and common garden‐reared individuals and emphasize the necessity of considering adaptive phenotypic plasticity in the study of adaptive divergence.     

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency‐dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color‐divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome‐wide markers from a double‐digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.     

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.     

Phylogenetic status of brown trout Salmo trutta populations in five rivers from the southern Caspian Sea and two inland lake basins,Iran: a morphogenetic approach

Interrelationships, origin and phylogenetic affinities of brown trout Salmo trutta populations from the southern Caspian Sea basin, Orumieh and Namak Lake basins in Iran were analysed from complete mtDNA control region sequences, 12 microsatellite loci and morphological characters. Among 129 specimens from six populations, seven haplotypes were observed. Based on mtDNA haplotype data, the Orumieh and southern Caspian populations did not differ significantly, but the Namak basin–Karaj population presented a unique haplotype closely related to the haplotypes of the other populations (0·1% Kimura two‐parameter, K2P divergence). All Iranian haplotypes clustered as a distinct group within the Danube phylogenetic grouping, with an average K2P distance of 0·41% relative to other Danubian haplotypes. The Karaj haplotype in the Namak basin was related to a haplotype (Da26) formerly identified in the Tigris basin in Turkey, to a Salmo trutta oxianus haplotype from the Aral Sea basin, and to haplotype Da1a with two mutational steps, as well as to other Iranian haplotypes with one to two mutational steps, which may indicate a centre of origin in the Caspian basin. In contrast to results of the mtDNA analysis, more pronounced differentiation was observed among the populations studied in the morphological and microsatellite DNA data, except for the two populations from the Orumieh basin, which were similar, possibly due to anthropogenic causes.     

Multilocus genetic diversity and historical biogeography of the endemic wall lizard from Ibiza and Formentera,Podarcis pityusensis (Squamata: Lacertidae)

Two monophyletic sister species of wall lizards inhabit the two main groups of Balearic Islands: Podarcis lilfordi from islets and small islands around Mallorca and Menorca and Podarcis pityusensis from Ibiza, Formentera and associated islets. Genetic diversity within the endangered P. lilfordi has been well characterized, but P. pityusensis has not been studied in depth. Here, 2430 bp of mtDNA and 15 microsatellite loci were analysed from P. pityusensis populations from across its natural range. Two main genetic groupings were identified, although geographical structuring differed slightly between the mtDNA and the nuclear loci. In general, individuals from islets/islands adjacent to the main island of Ibiza were genetically distinct from those from Formentera and the associated Freus islands for both mtDNA and the nuclear loci. However, most individuals from the island of Ibiza were grouped with neighbouring islets/islands for nuclear loci, but with Formentera and Freus islands for the mitochondrial locus. A time‐calibrated Bayesian tree was constructed for the principal mitochondrial lineages within the Balearics, using the multispecies coalescent model, and provided statistical support for divergence of the two main P. pityusensis lineages 0.111–0.295 Ma. This suggests a mid‐late Pleistocene intraspecific divergence, compared with an early Pleistocene divergence in P. lilfordi, and postdates some major increases in sea level between 0.4 and 0.6 Ma, which may have flooded Formentera. The program IMa2 provided a posterior divergence time of 0.089–0.221 Ma, which was similar to the multispecies coalescent tree estimate. More significantly, it indicated low but asymmetric effective gene copy migration rates, with higher migration from Formentera to Ibiza populations. Our findings suggest that much of the present‐day diversity may have originated from a late Pleistocene colonization of one island group from the other, followed by allopatric divergence of these populations. Subsequent gene flow between these insular groups seems likely to be explained by recent human introductions. Two evolutionary significant units can be defined for P. pityusensis but these units would need to exclude the populations that have been the subjects of recent admixture.     

Comparative phylogeographies of six species of hinged terrapins (Pelusios spp.) reveal discordant patterns and unexpected differentiation in the P. castaneus/P. chapini complex and P. rhodesianus

Using up to 2117 bp of mitochondrial DNA and up to 2012 bp of nuclear DNA, we analysed phylogeographic differentiation of six widely distributed species of African hinged terrapins (Pelusios spp.) representing different habitat types. Two taxa each live in savannahs or in forests and mesic savannahs, respectively, and the remaining two species occur in intermediate habitats. The species living in forests and mesic savannahs do not enter dry savannahs, whereas the savannah species may occur in dry and wet savannahs and even in semi‐arid steppe regions. We found no obvious correlation between habitat type and phylogeographic pattern: one savannah species (P. rhodesianus) shows phylogeographic structure, i.e. pronounced genetic differences among geographically distinct populations, and the other (P. nanus) not. One species inhabiting forests and mesic savannahs (P. carinatus) has phylogeographic structure, the other (P. gabonensis) not. The same pattern is true for the two ecologically intermediate species, with phylogeographic structure present in P. castaneus and absent in P. chapini. Nuclear evidence suggests that the latter two taxa with abutting and partially overlapping ranges are distinct, while mtDNA is only weakly differentiated. Pelusios castaneus shows pronounced phylogeographic structure, which could reflect Pleistocene range interruptions correlated with the fluctuating forest cover in West and Central Africa. Our results do not support the recognition of an extinct subspecies of P. castaneus for the Seychelles. Pelusios carinatus contains two well supported clades, which are separated by the Congo River. This species is closely related to P. rhodesianus, a taxon consisting of two deeply divergent mitochondrial clades. One of these clades is paraphyletic with respect to P. carinatus, but the two clades of P. rhodesianus are not differentiated in the studied nuclear markers and, again, paraphyletic with respect to P. carinatus. Using mtDNA sequences from the type material of P. rhodesianus, we were able to allocate this name to one of the two clades. However, owing to the confusing relationships of P. rhodesianus and P. carinatus, we refrain from taxonomic decisions.