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.

Østbye K ，Bernatchez L ，Naesje TF ，Himberg KJ ，Hindar K ... -  《MOLECULAR ECOLOGY》

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.

Babik W ，Branicki W ，Crnobrnja-Isailović J ，Cogălniceanu D ，Sas I ，Olgun K ，Poyarkov NA ，Garcia-París M ，Arntzen JW ... -  《MOLECULAR ECOLOGY》

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.

Babik W ，Branicki W ，Sandera M ，Litvinchuk S ，Borkin LJ ，Irwin JT ，Rafiński J ... -  《MOLECULAR ECOLOGY》

Phylogeography of the fire-bellied toads Bombina: independent Pleistocene histories inferred from mitochondrial genomes.

The fire-bellied toads Bombina bombina and Bombina variegata, interbreed in a long, narrow zone maintained by a balance between selection and dispersal. Hybridization takes place between local, genetically differentiated groups. To quantify divergence between these groups and reconstruct their history and demography, we analysed nucleotide variation at the mitochondrial cytochrome b gene (1096 bp) in 364 individuals from 156 sites representing the entire range of both species. Three distinct clades with high sequence divergence (K2P = 8-11%) were distinguished. One clade grouped B. bombina haplotypes; the two other clades grouped B. variegata haplotypes. One B. variegata clade included only Carpathian individuals; the other represented B. variegata from the southwestern parts of its distribution: Southern and Western Europe (Balkano-Western lineage), Apennines, and the Rhodope Mountains. Differentiation between the Carpathian and Balkano-Western lineages, K2P approximately 8%, approached interspecific divergence. Deep divergence among European Bombina lineages suggests their preglacial origin, and implies long and largely independent evolutionary histories of the species. Multiple glacial refugia were identified in the lowlands adjoining the Black Sea, in the Carpathians, in the Balkans, and in the Apennines. The results of the nested clade and demographic analyses suggest drastic reductions of population sizes during the last glacial period, and significant demographic growth related to postglacial colonization. Inferred history, supported by fossil evidence, demonstrates that Bombina ranges underwent repeated contractions and expansions. Geographical concordance between morphology, allozymes, and mtDNA shows that previous episodes of interspecific hybridization have left no detectable mtDNA introgression. Either the admixed populations went extinct, or selection against hybrids hindered mtDNA gene flow in ancient hybrid zones.

Hofman S ，Spolsky C ，Uzzell T ，Cogălniceanu D ，Babik W ，Szymura JM ... -  《MOLECULAR ECOLOGY》

DNA taxonomy, phylogeny and Pleistocene diversification of the Cicindela hybrida species group (Coleoptera: Cicindelidae).

Species delimitation is complicated where morphological variation is continuous or poorly subdivided, but for taxonomic convenience it is common practice to separate and name geographical groups to capture this variation. DNA-based approaches may be used to test if these groups in fact represent historically divided, discrete species entities. The Cicindela hybrida complex (Coleoptera: Cicindelidae) is an assemblage of up to seven morphologically recognized species and 15 subspecies with wide distribution in the Palaearctic region. We sequenced a discontinuous segment of 1899 bp of mtDNA including three regions (coxI, rrnL+trnL2+nad1, cob) for a total of 99 specimens from 36 sampling localities across Europe, revealing 48 haplotypes. Four major clades could be identified corresponding to geographical groups from central Iberia, Ukraine, central Europe, and a band from the Atlantic Iberian coast to northern Europe. Taking into account further subdivisions within these clades, four of the six named species included in the analysis were recognizable by applying various procedures for species delimitation. Age estimates from calibrated molecular clocks date the diversification of the hybrida group within the past 2 million years (Myr), and the separation of the northern clade within 0.4 Myr. Nested clade analysis revealed the rapid range expansion of the northern group consistent with postglacial dispersal, but we did not find support for specific source population(s) in the postulated southern refugia. The evolutionary framework based on mtDNA sequences is shown to identify species entities as discrete clusters of closely related sequences and provides an objective system for delineating and recognizing hierarchically structured groups. In the case of the C. hybrida complex, these groups largely coincided with those established from morphology. The study adds further support to the utility of mtDNA-based sequence profiles (the 'DNA taxonomy') as a rapid and objective synthesis of evolutionary diversity and as reference system for communication.

Cardoso A ，Vogler AP 《MOLECULAR ECOLOGY》