Mitochondrial DNA variation in the hybridizing fire-bellied toads, Bombina bombina and B. variegata

Using five restriction enzymes, geographical variation of mitochondrial DNA (mtDNA) in Bombina bombina and B. variegata was studied in samples from 20 locations. Each restriction enzyme produced a species-specific fragment pattern. B. bombina haplotypes A and B were closely related to each other. In contrast, haplotypes A and B of B. variegata formed two distinct lineages. A very distinctive haplotype (C) was found in the Carpathian Mountains, whereas two other haplotypes, D and E (differing by a single AvaI site), were present in western Europe and the Balkans, respectively. Populations polymorphic for haplotypes D and E occurred in the central Balkans where the haplotypes could replace each other clinally. mtDNA sequence divergence between B. bombina and B. variegata was estimated as 6.0-8.1% and 4.7-5.2% between type C and types D/E of B. variegata. The latter divergence is contrary to allozyme and morphological data that place the western and Carpathian B. v. variegata together (Nei's D = 0.07) and separate them from the Balkan subspecies B. v. scabra (Nei's D = 0.18). Broad interspecific correlation among morphology, allozymes and mtDNA types in European fire-bellied toads argues that, despite continuous hybridization (interrupted perhaps during Pleistocene glacial maxima), little or no mtDNA introgression between the species has occurred outside the narrow hybrid zones that separate these parapatric species.     

Mapping Bombina mitochondrial genomes: the conundrum of Carpathian Bombina variegata (Anura: Discoglossidae)

Complete mitochondrial genomes (mtDNA) of five individuals representing two haplotypes of Bombina bombina and three of Bombina variegata were compared using restriction site maps. Phylogenetic analyses reveal three ancient mitochondrial lineages: (1) two very similar haplotypes A and B of B. bombina ; (2) almost identical haplotypes D and E of B. variegata ; and (3) haplotype C of B. variegata . Haplotype C is as different from haplotypes D/E as from A/B. These data are strikingly discordant with relationships based on morphology and allozymes. Haplotypes C and D/E represent a pre-Pleistocene mitochondrial divergence within B. variegata , nearly coincident with speciation between B. variegata and B. bombina . Geographical partitioning of the two divergent B. variegata mitochondrial lineages indicates repeated localization of the lineages in separate glacial refugia during the Pleistocene. That nuclear genes do not show a similar divergence, but rather indicate relatively free genetic exchange between populations with divergent mtDNAs, suggests that males dispersed much more widely than females during expansions from glacial refugia. Comparison of Bombina mtDNA maps with a restriction site map of Xenopus laevis mtDNA revealed 16 homologous sites; 12 of these may be nearly invariant across primitive anuran mtDNAs. Two distinct regions of heteroplasmy, representing two regions with variable numbers of sequence repeats [length variable (LV) regions], were characterized. Comparison with the Xenopus map places LV1, present in all five haplotypes, near the 5'-end of the control region, and LV2 present only in B. variegata , near the 3'-end. Although phylogenetic analyses did not group the two major B. variegata lineages together, presence of LV2 in both lineages supports placement of both within B. variegata .     

Evolution of Bombina bombina and Bombina variegata (Anura: Discoglossidae) in the Carpathian Basin: a history of repeated mt-DNA introgression across species

The structure and geographic location of hybrid zones change through time. Current patterns result from present and historical population-environment interactions that act on each of the hybridizing taxa. This is particularly evident for species involved in complex hybrid zones, such as that formed by the toad species Bombina bombina and Bombina variegata (Anura: Discoglossidae), which interact along extensive areas in Central Europe. We used data on external morphology and partial sequences of the cytochrome oxidase I (cox1) and nicotinamide adenine dinucleotid dehydrogenase subunit 4 (nad4) mitochondrial DNA (mt-DNA) genes to analyze the current patterns of genetic structure shown by both species of Bombina along their contact zone in Hungary. Phylogenetic, phylogeographic, and historical demography analyses were applied to 1.5kb mt-DNA obtained from 119 individuals representing 24 populations from Hungary and additional specimens from Slovakia, Albania, and Bosnia-Herzegovina. We use these data to infer the evolutionary history of the isolated populations of B. variegata in Hungary and to discriminate between competing biogeographic scenarios accounting for the historical interactions between species in this region. Results from the inferred phylogenetic branching pattern and sequence divergence among species and populations support the following: (i) recent population expansion has occurred in Hungarian populations of B. bombina, which are genetically very homogeneous; (ii) the Hungarian populations of B. variegata correspond to two distinct mitochondrial lineages (Carpathian and Alpine, respectively); average maximum-likelihood-corrected sequence divergence between these lineages is 8.96% for cox1 and 10.85% for nad4; (iii) mt-DNA divergence among the three isolated western populations of B. variegata from Transdanubia is low, with four closely related haplotypes, which suggests that the isolation between these populations is the result of a recent process, possibly mediated by the invasion of B. bombina; and (iv) we have detected discordances between morphology and mt-DNA data in the Transdanubia region (Bakony Mountains, Mecsek Mountains, Orség area), suggesting mt-DNA introgression across species in this regions. These results are discussed with reference to previous biogeographic hypotheses.     

Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

Limited mitochondrial DNA introgression in a Bombina hybrid zone

We analysed the variation at mitochondrial DNA (mtDNA) and six allozyme loci in two transects across a hybrid zone between the fire-bellied toads Bombina bombina and Bombina variegata in southern Poland. The mtDNA cline was narrower than allozyme clines in one transect (Przemyśl, β = 0.435, w  = 3.42 km) and shifted to the B .  bombina side in both (near Kraków by 1.32 km, α = 0.437, and near Przemyśl by 0.97 km, α = 0.319). Cytonuclear associations were weak. Narrowed mtDNA clines could be a by-product of female demography and lowered effective population size of mtDNA. Alternatively, restricted mtDNA introgression is a likely result of negative epistasis in recombinants augmented by environment-dependent selection because divergence of Bombina mitochondrial and nuclear genomes is large. The shift of the mtDNA cline, contrary to expectations from initial mating preferences and fecundity differences between the species, suggests that on the B .  bombina side of the zone hybrid females with B .  variegata mtDNA have a higher chance of leaving progeny.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 295–306.     

Reproductive biology of the fire-bellied toads, Bombina bombina and B. variegata (Anura: Discoglossidae): egg size, clutch size and larval period length differences

Clutch size, egg size, larval size and duration of larval period were estimated for samples of Bombina bombina and B. variegata from several populations. Bombina variegata lays fewer eggs per batch (mean 17.4) than B. bombina (mean 32.5). The volume of B. variegata eggs is 2.4 times greater and the dry weight 2.3 times greater than that of B. bombina. As a consequence, the volume of B. variegata hatchlings and their dry weight are respectively 1.9 times and 2.5 times greater than that of B. bombina. Analysis of variance showed that variation of egg size in B. variegata is significantly greater than in B. bombina. When reared under laboratory conditions freshly metamorphosed individuals of both species showed no difference in body size or weight. Although growth rates are equal, tadpoles of B. variegata reached metamorphosis 14 days earlier than B. bombina tadpoles because of differences in size of hatchlings. Differences in reproductive parameters found for the two European Bombina species can be interpreted as adaptations to contrasting environments in which they breed.     

Mitochondrial phylogeography of the moor frog, Rana arvalis

The moor frog Rana arvalis is a lowland species with a broad Eurasiatic distribution, from arctic tundra through forest to the steppe zone. Its present-day range suggests that glacial refugia of this frog were located outside southern European peninsulas. We studied the species-wide phylogeographical pattern using sequence variation in a 682 base pairs fragment of mtDNA cytochrome b gene; 223 individuals from 73 localities were analysed. Two main clades, A and B, differing by c. 3.6% sequence divergence were detected. The A clade is further subdivided into two subclades, AI and AII differing by 1.0%. All three lineages are present in the Carpathian Basin (CB), whereas the rest of the species range, including huge expanses of Eurasian lowlands, are inhabited solely by the AI lineage. We infer that AII and B lineages survived several glacial cycles in the CB but did not expand, at least in the present interglacial, to the north. The geographical distribution and genealogical relationships between haplotypes from the AI lineage indicate that this group had two glacial refugia, one located in the eastern part of the CB and the other probably in southern Russia. Populations from both refugia contributed to the colonization of the western part of the range, whereas the eastern part was colonized from the eastern refugium only. The effective population size as evidenced by theta(ML) is an order of magnitude higher in the AI lineage than in the AII and B lineages. Demographic expansion was detected in all three lineages.     

Quantitative immunological studies of the albumins of several species of fire bellied toads, genus Bombina.

1. Rabbit antisera against purified serum albumin of Bombina bombina were used to study relationships between B. bombina, B. variegata, and B. orientalis. 2. Quantitative micro-complement fixation tests indicated the albumins of B. bombina from central Poland and Bulgaria were indistinguishable. The albumins from several populations of B. variegata differed very slightly from that of B. bombina. The albumin of B. orientalis was quite distinct from that of B. bombina. 3. Using albumin as a molecular clock, we estimated B. bombina and B. variegata diverged within the last million years, whereas the B. orientalis lineage diverged roughly 10-12 mil yr ago.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Comparative analysis of the genetic structure of hybrid populations of European toads of the genus Bombina (Oken, 1816) from Ukraine and Croatia

The comparative analysis of hybrid populations of toads Bombina bombina and B. variegata from Ukraine and Croatia revealed the essential differences in their genetic structure. The population composition and the direction of genetic introgression are determined initially by the ecological factors (the locality type, the country relief). The clear predominance of the red-bellied toad Bombina bombina was found in two cases: in the lowland populations and in the big samples. The second parental species, B. variegata, predominates in the highlands and in the small samples, respectively.     

Post-glacial expansion into the Paleozoic Plateau: evidence of an Ozarkian refugium for the Ozark minnow Notropis nubilus (Teleostei: Cypriniformes)

Genetic variation was examined within the Ozark minnow Notropis nubilus using complete mtDNA cytochrome b gene sequences from 160 individuals representing 30 localities to test hypotheses on the origin of the distribution. Phylogenetic analyses revealed three strongly supported clades of haplotypes consistent with geographic distributions: a clade from the Western Ozarks, a clade from the Southern Ozarks and a clade from the Northern Ozarks and upper Mississippi River basin. The estimated mean ages of these clades indicated that they diverged during pre-Illinoian glacial cycles extending from the late Pliocene into the early Pleistocene. Results of demographic analyses based on coalescent approaches supported the hypothesis that the Paleozoic Plateau was not a refugium for N. nubilus during periodic glacial advances. There is evidence of a genetic signature of northern expansion into the Paleozoic Plateau from a Southern Ozarkian refugium. Populations expanded out of drainages in the Northern Ozarks into the Paleozoic Plateau during the late Pleistocene. Subsequently, the two regions were isolated due to the recent extirpation of intervening populations caused by the loss of suitable habitat.     

Post-glacial colonization of northwestern North America by the forest-associated American marten (Martes americana,Mammalia: Carnivora: Mustelidae)

Phylogeographic patterns were used to assess intraspecific diversification of American martens (Martes americana). Within martens, two morphological groups (americana and caurina) have been recognized, though the level of distinction between them has been debated. We examined mitochondrial cytochrome b gene haplotypes from 680 martens to explore the colonization history of the Pacific Northwest and found two clades that correspond to the morphological groups. The widespread americana clade extends from interior Alaska south to Montana and eastward to Newfoundland and New England (i.e. northwestern, north-central and northeastern North America). The caurina clade occurs in western North America, minimally extending from Admiralty Island (southeastern Alaska) south to Oregon and Wyoming. Our data indicated two colonization events for the Pacific Northwest (one by members of each clade) and were consistent with the persistence of populations throughout past glacial periods in eastern and western refugia. Due to vegetational and geological history following the past deglaciation, we hypothesize that martens of the caurina clade spread along the North Pacific Coast, and into southeastern Alaska, earlier than martens of the americana clade. Mismatch distributions for the americana clade were indicative of populations that recently experienced demographic expansion, while mismatch distributions for the caurina clade suggested that populations were at equilibrium. These clades are reciprocally monophyletic and distinctive (interclade divergence ranged from 2.5 to 3.0% (uncorrected p), whereas, intraclade divergence was < 0.7%), and two regions of sympatry have been identified. Genetic signatures of past admixture in hybrid zones may have been extinguished during subsequent glacial periods when ranges contracted. This recurrent pattern of relatively restricted western, or Pacific coastal, lineages and more widespread eastern, or interior continental, lineages exists across broad taxonomic groups and suggests a shared biogeographical history.     

Phylogeography of sable (Martes zibellina L. 1758) in the southeast portion of its range based on mitochondrial DNA variation: highlighting the evolutionary history of the sable

To reveal phylogeographic features of sable (Martes zibellina) in the southeast part of its range, we analyzed variability of the mitochondrial DNA (mtDNA) cytochrome b gene, tRNA (Pro), tRNA (Thr) and control region (D-loop) sequences from 78 specimens in populations of the Russian Far East, northeast China, and Mongolia. Our results revealed the presence of 49 different haplotypes split into two major phylogenetic groups—clades A and B, the latter separated into two clades, B1 and B2. Comparative analysis of D-loop haplotypes in populations originating from the southeast (Russian Far East, China and Mongolia) and the west (northern Urals) portions of sable range indicated that all three mtDNA clades were present in different regional groups. However, highest diversity of clade B1 in northeast China and its nearly complete absence from the Urals suggest that the southeast sable range, being a refuge during Pleistocene glacial periods, can be considered the center of genetic diversification and possibly origin of this species. All divergence estimates fall within the Pleistocene suggesting that Quarternary glaciations played an important role in phylogeographic differentiation of sable.     

Phylogenetic relationships among bisexual and unisexual lineages of the weevil Scepticus insularis (Coleoptera: Curculionidae) on Hokkaido, Northern Japan

We examined the geographic distributions and phylogenetic relationships of bisexual and unisexual (parthenogenetic) forms of the weevil Scepticus insularis on Hokkaido Island, northern Japan. Unisexual beetles were widespread throughout Hokkaido, whereas bisexuals were found only in three remote areas. Bisexuals (females and males) and unisexual females occurred sympatrically in two areas. We determined nucleotide sequences for part of the mitochondrial NADH dehydrogenase subunit 2 (ND2; 633 bp) gene for 104 individuals, and for part of the nuclear internal transcribed spacer region 2 (ITS2; 360 bp) for 91 individuals. In an ND2 gene tree, haplotypes of S. insularis fell into two distinct clades (A and B), which were genetically differentiated from one another by 9.1% nucleotide sequence divergence. Haplotypes of females identified as unisexual were all in clade A, whereas those of females identified as bisexual belonged to clade B. Haplotypes of males were in clade B, except for two males having a clade-A haplotype. Circumstantial evidence suggests that these two males were produced by unisexual females. The ND2 tree suggests that the current unisexual form of S. insularis on Hokkaido was of a single origin. In contrast, a gene tree for ITS2 haplotypes show no clear divergence between the two modes of reproduction, with two major haplotypes shared by unisexual females, bisexual females, and males. This incongruence between the nuclear and mitochondrial phylogenies may be attributable to occasional gene flow between the unisexual and bisexual lineages through males occasionally produced by unisexual females, but our results do not exclude the possibility that the two lineages share polymorphic ancestral ITS2 haplotypes.     

The genetic legacy of Mother Goose– phylogeographic patterns of lesser snow goose Chen caerulescens caerulescens maternal lineages

By using the polymerase chain reaction to amplify and sequence 178 bp of a rapidly evolving region of the mtDNA genome (segment I of the control region) from 81 individuals, approximately 11% of the variation present in the lesser snow goose Chen caerulescens caerulescens L. mitochondrial genome was surveyed. The 26 types of mtDNA detected formed two distinct mitochondrial clades that differ by an average of 6.7% and are distributed across the species range. Restriction analysis of amplified fragments was then used to assign the mtDNA of an additional 29 individuals to either of these clades. Within one major clade, sequence among mtDNAs was concordant with geographic location. Within the other major clade the degree of sequence divergence among haplotypes was lower and no consistent geographic structuring was evident. The two major clades presumably result from vicariant separation of lesser snow geese during the Pleistocene.     

Comparative phylogeography of the two pink salmon broodlines: an analysis based on a mitochondrial DNA genealogy

Over most of their natural northern Pacific Ocean range, pink salmon (Oncorhynchus gorbuscha) spawn in a habitat that was repeatedly and profoundly affected by Pleistocene glacial advances. A strictly two-year life cycle of pink salmon has resulted in two reproductively isolated broodlines, which spawn in alternating years and evolved as temporal replicates of the same species. To study the influence of historical events on phylogeographical and population genetic structure of the two broodlines, we first reconstructed a fine-scale mtDNA haplotype genealogy from a sample of 80 individuals and then determined the geographical distribution of the major genealogical assemblages for 718 individuals sampled from nine Alaskan and eastern Asian even- and nine odd-year pink salmon populations. Analysis of restriction site states in seven polymerase chain reaction (PCR)-amplified mtDNA regions (comprising 97% of the mitochondrial genome) using 13 endonucleases resolved 38 haplotypes, which clustered into five genealogical lineages that differed from 0.065 to 0.225% in net sequence divergence. The lineage sorting between broodlines was incomplete, which suggests a recent common ancestry. Within each lineage, haplotypes exhibited star-like genealogies indicating recent population growth. The depth of the haplotype genealogy is shallow ( approximately 0.5% of nucleotide sequence divergence) and probably reflects repeated decreases in population size due to Pleistocene glacial advances. Nested clade analysis (NCA) of geographical distances showed that the geographical distribution observed for mitochondrial DNA (mtDNA) haplotypes resulted from alternating influences of historical range expansions and episodes of restricted dispersal. Analyses of molecular variance showed weak geographical structuring of mtDNA variation, except for the strong subdivision between Asian and Alaskan populations within the even-year broodline. The genetic similarities observed among and within geographical regions probably originated from postglacial recolonizations from common sources rather than extensive gene flow. The phylogeographical and population genetic structures differ substantally between broodlines. This can be explained by stochastic lineage sorting in glacial refugia and perhaps different recolonization routes in even- and odd-year broodlines.     

Cryptic niche conservatism among evolutionary lineages of an invasive lizard

Aim There is increasing evidence that the quality and breadth of ecological niches vary among individuals, populations, evolutionary lineages and therefore also across the range of a species. Sufficient knowledge about niche divergence among clades might thus be crucial for predicting the invasion potential of species. We tested for the first time whether evolutionary lineages of an invasive species vary in their climate niches and invasive potential. Furthermore, we tested whether lineage‐specific models show a better performance than combined models. Location Europe. Methods We used species distribution models (SDMs) based on climatic information at native and invasive ranges to test for intra‐specific niche divergence among mitochondrial DNA (mtDNA) clades of the invasive wall lizard Podarcis muralis. Using DNA barcoding, we assigned 77 invasive populations in Central Europe to eight geographically distinct evolutionary lineages. Niche similarity among lineages was assessed and the predictive power of a combination of clade‐specific SDMs was compared with a combined SDM using the pooled records of all lineages. Results We recorded eight different invasive mtDNA clades in Central Europe. The analysed clades had rather similar realized niches in their native and invasive ranges, whereas inter‐clade niche differentiation was comparatively strong. However, we found only a weak correlation between geographic origin (i.e. mtDNA clade) and invasive occurrences. Clades with narrow realized niches still became successful invaders far outside their native range, most probably due to broader fundamental niches. The combined model using data for all invasive lineages achieved a much better prediction of the invasive potential. Conclusions Our results indicate that the observed niche differentiation among evolutionary lineages is mainly driven by niche realization and not by differences in the fundamental niches. Such cryptic niche conservatism might hamper the success of clade‐specific niche modelling. Cryptic niche conservatism may in general explain the invasion success of species in areas with apparently unsuitable climate.     

Phylogeography of two European newt species--discordance between mtDNA and morphology

The newts Triturus vulgaris and Triturus montandoni are sister species that exhibit contrasting levels of intraspecific morphological variation. Triturus vulgaris has a broad Eurasiatic distribution encompassing both formerly glaciated and unglaciated areas and shows substantial morphological differentiation in the southern part of its range, while T. montandoni, confined to the Carpathians, is morphologically uniform. We analysed sequence variation of two mtDNA fragments of the total length of c. 1850 bp in 285 individuals of both species collected from 103 localities. Phylogenetic analysis of 200 unique haplotypes defined 12 major clades, their age estimated at c. 4.5-1.0 million years (Myr). Most of the older clades were found in the southern part of the range, and also in central Europe, mainly in Romania. The distribution of mtDNA clades points to the existence of several glacial refugia, located in the Caucasus region, Anatolia, the Balkan Peninsula, Italy, and more to the north in central Europe. The concordance between mtDNA based phylogeny and the distribution of T. vulgaris subspecies was weak. Triturus montandoni haplotypes did not form a monophyletic group. Instead they were found in six clades, in five of them mixed with T. vulgaris haplotypes, most likely as a result of past or ongoing hybridization and multiple introgression of mtDNA from T. vulgaris to T. montandoni. Patterns of sequence variation within clades suggested long-term demographic stability in the southern groups, moderate and relatively old demographic growth in the populations inhabiting central Europe, and high growth in some of the groups that colonized northern parts of Europe after the last glacial maximum.     

Living with the genetic signature of Miocene induced change: evidence from the phylogeographic structure of the endemic angulate tortoise Chersina angulata

The phylogeographic structure of the monotypic endemic southern African angulate tortoise Chersina angulata was investigated throughout its distribution with the use of partial sequences from three mtDNA loci (COI, cyt b and ND4). Phylogeographic and phylogenetic structuring obtained for the three mtDNA markers were highly congruent and suggested the presence of two genetically distinct, reciprocally monophyletic evolutionary lineages. Group one contained two subclades with haplotypes from the north-western Cape and south-western Cape, respectively, while haplotypes from the southern Cape comprised group two. The two major clades were separated by nine and eight mutational steps for COI and ND4, respectively. Of the three mtDNA gene regions examined, the ND4 partial sequence contained the most phylogenetic signal. Haplotype diversity was generally low and we recovered 34 haplotypes for the 125 animals sequenced for the ND4 subunit. Nested clade analyses performed on the variable ND4 partial sequences suggested the presence of two major refugial areas for this species. The demographic history of the taxon was characterised by range expansion and prolonged historical fragmentation. Divergence time estimates suggest that the temporal and spatial distribution of the taxon was sculpted by changes in temperature and rainfall patterns since the late Miocene. Corroborative evidence from other reptiles is also suggestive of a late Miocene divergence, indicating that this was a major epoch for cladogenesis in southern Africa. Apart from the genetic differences between the two major clades, we also note morphometric and behavioural differences, alluding to the presence of two putative taxa nested within C. angulata.     

Exploring factors shaping population genetic structure of the freshwater fish Sinocyclocheilus grahami (Teleostei, Cyprinidae)

Phylogeographical analyses on Sinocyclocheilus grahami samples from seven localities within the Lake Dianchi Basin in China were conducted to explore the main factors shaping population structure within this species. Phylogenetic and network analyses revealed two major clades in 24 mtDNA haplotypes. One clade included three haplotypes exclusively from samples of the lower basin and another clade encompassed other haplotypes from samples of the upper basin. The estimated divergence time between the two clades predated the river capture event connecting the lower and upper lake basin and thus supported geographical isolation as the main factor shaping genetic divergence between these two clades. Furthermore, analysis of molecular variance and pair-wise Φ_ST distances revealed significant genetic differentiation within the upper basin. Mantel tests clearly supported patterns of differentiation arose purely as a result of isolation by distance. These results further highlight the importance of geographical isolation in shaping differentiation within this species.