Phylogeography of a northeast Asian spruce, Picea jezoensis, inferred from genetic variation observed in organelle DNA markers

Range-wide genetic variation of the widespread cold-temperate spruce Picea jezoensis was studied throughout northeast Asia using maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA markers. This study assessed 33 natural populations including three varieties of the species in Japan, Russia, China, and South Korea. We depicted sharp suture zones in straits around Japan in the geographical distribution pattern of mitochondrial haplotypes (GST=0.901; NST=0.934). In contrast, we detected possible extensive pollen flow without seed flow across the straits around Japan during the past population history in the distribution pattern of chloroplast haplotypes (GST=0.233; NST=0.333). The analysis of isolation by distance of the species implied that by acting as a barrier for the movement of seeds and pollen, the sharp suture zones contributed considerably to the level of genetic differentiation between populations. Constructed networks of mitochondrial haplotypes allowed inference of the phylogeographical history of the species. We deduced that the disjunction with Kamchatka populations reflects range expansion and contraction to the north of the current distribution. Within Japan, we detected phylogeographically different types of P. jezoensis between Hokkaido and Honshu islands; P. jezoensis in Honshu Island may have colonized this region from the Asian continent via the Korean peninsula and the species in Hokkaido Island is likely to have spread from the Asian continent via Sakhalin through land bridges. Japanese endemism of mitochondrial haplotypes in Hokkaido and Honshu islands might have been promoted by separation of these islands from each other and from the Asian continent by the straits during the late Quaternary.     

Intraspecific sequence variation of chloroplast DNA inPedicularis chamissonis Steven (Scrophulariaceae) and geographic structuring of the Japanese “Alpine” plants

In order to clarify evolutionary patterns and processes of intraspecific diversification ofPedicularis chamissonis Steven, we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers betweentrnT (UGU) andtrnL (UAA) 5′exon,trnL (UAA) 3′exon andtrnF (GAA), andatpB andrbcL. In 24 populations ofP. chamissonis, 33 nucleotide substitutions and 12 insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.014. Seventeen distinct cpDNA haplotypes could be recognized and each haplotype was found to be geographically structured. Two major clades (the Northern and Southern clades) were revealed in phylogenetic analyses of cpDNA haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mts. lide of central Honshu in Japan, northward to Unalaska Island in the Aleutians. Relationships among most haplotypes were unresolved polytomies. On the other hand, the haplotypes of the Southern clade occurred from the populations of Mt. Gassan southwards to Mt. Arakawa of central Honshu. Within this clade, three subclades were clearly recognized. From these results, we concluded that the haplotypes of the Northern and Southern clades inP. chamissonis might have traveled down to Japanese Archipelago from the north in not a single glacial period.     

Further Analysis of Intraspecific Sequence Variation of Chloroplast DNA in Primula cuneifolia Ledeb. (Primulaceae): Implications for Biogeography of the Japanese Alpine Flora

Primula cuneifolia Ledeb. (Primulaceae), we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers between trnT (UGU) and the trnL (UAA) 5′exon, the trnL (UAA) 3′exon and trnF (GAA), and atpB and rbcL. In 20 populations of P. cuneifolia, 22 nucleotide substitutions and five insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.008. Eight distinct haplotypes could be recognized and each haplotype was found to be geographically structured. Three major clades (the Northern, Hokkaido and Southern clades) were revealed in phylogenetic analyses of the haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mt. Rausu and Rishiri Island of eastern and northern Hokkaido in Japan, northward to Unalaska Island in the Aleutians, and those of the Hokkaido clade were distributed in the populations of central Hokkaido and Mt. Iwaki of the northern Honshu in Japan; in addition, those of the Southern clade were observed only in the populations of the central Honshu. It was shown that the genetic diversifications of the Southern clade were higher than those of the Northern and Hokkaido clades. Furthermore, it was shown that the topology within the Southern clade was hierarchical, and the haplotypes of the Southern populations in the clade were derivative. From these results, we concluded that the cpDNA haplotypes of the three clades in P. cuneifolia arose and assumed the present distribution areas through several cycles of glacial advance and retreat in the Pleistocene. Received 24 June 1998/ Accepted in revised form 28 December 1998     

Genetic structure of Japanese populations of an ambrosia beetle,Xylosandrus germanus (Curculionidae: Scolytinae)

We examined the genetic structures of 13 Japanese populations of an ambrosia beetle, Xylosandrus germanus (Curculionidae: Scolytinae), to understand the effects of geographical barriers on the colonization dynamics of this species. The genetic structure was studied using portions of the mitochondrial cytochrome oxidase I (COI) gene. A phylogenetic analysis revealed three distinct lineages (clades A, B and C) within X. germanus. Clade A contained 21 haplotypes from all 13 populations; whereas clade B contained eight haplotypes from Hokkaido (Sapporo and Furano), Iwate and Nagano populations; and clade C contained only a single a haplotype from the Hokkaido (Furano) population. In the analysis of molecular variance (amova ), the greatest amount of genetic variation was detected between populations in Hokkaido and those in Honshu and other southern islands. Between these two groups of populations, all the values of the coefficient of gene differentiation were significantly larger than zero, except for the Hokkaido (Sapporo) versus Nagano comparison. Our results confirm that for X. germanus, gene flow has been interrupted between Hokkaido and Honshu since the last glacial maximum.     

Evolution of the Cyprinella lutrensis species group. III. Geographic variation in the mitochondrial DNA of Cyprinella lutrensis— the influence of Pleistocene glaciation on population dispersal and divergence

We employed restriction site variation in mitochondrial (mt)DNA to determine if significant phylogeographic structure occurs in the North American cyprinid fish Cyprinella lutrensis. Digestion patterns from 16 restriction endonucleases identified fifty mtDNA haplotypes among 127 individuals of Cyprinella lutrensis assayed from localities in the Gulf Coastal Plain, the Great Plains, and the Central Lowlands. Nucleotide sequence divergence among haplotypes was highly variable (mean ± SE: 2.87%± 0.08; range: 0.14–9.24%). Maximum-parsimony analysis and the neighbour joining method of tree construction revealed three major groupings (clades) of haplotypes that differed in geographic distribution. Divergence estimates between the basal clade, comprised of haplotypes primarily from the Brazos River in east Texas, and the remaining two clades, place C. lutrensis in the western Gulf Coastal Plain prior to Pleistocene glaciation. Nucleotide sequence divergence between the second clade, comprised of haplotypes from the Trinity and Calcasieu rivers in east Texas and southwestern Louisiana, respectively, and the third clade (comprised primarily of haplotypes from localities north of Texas and affected directly by Pleistocene glaciation), suggest that C. lutrensis colonized gladated regions to the north during the mid- to late Pleistocene. This hypothesis is supported by levels of intrapopulational nucleotide diversity in geographic localities outside of Texas and by geological evidence. Despite marked geographic variation in morphometries, meristics, and nuptial coloration, mtDNA variation in glaciated regions was not geographically structured, and subspecies of C. lutrensis were not identifiable by phylogenetic analysis of mtDNA.     

Complete taxon sampling of the avian genus Pica (magpies) reveals ancient relictual populations and synchronous Late‐Pleistocene demographic expansion across the Northern Hemisphere

Previous studies have suggested that bird populations in east Asia were less affected by Pleistocene climatic fluctuations than those in Europe and North America. However, this is mainly based on comparisons among species. It would be more relevant to analyse geographical populations of widespread species or species complexes. We analyzed two mitochondrial genes and two nuclear introns for all taxa of Pica to investigate 1) which Earth history factors have shaped the lineage divergence, and 2) whether different geographical populations were differently affected by the Pleistocene climatic changes. Our mitochondrial tree recovered three widespread lineages, 1) in east Asia, 2) across north Eurasia, and 3) in North America, respectively, with three isolated lineages in northwest Africa, Arabia and the Qinghai‐Tibet Plateau, respectively. Divergences among lineages took place 1.4–3.1 million yr ago. The northwest African population was sister to the others, which formed two main clades. In one of these, Arabia was sister to Qinghai‐Tibet, and these formed the sister clade to the east Asia clade. The other main clade comprised the North American and north Eurasian clades. There was no or very slight structure within these six geographical clades, including a lack of differentiation between the two North American species black‐billed magpie P. hudsonia and yellow‐billed magpie P. nutalli. Demographic expansion was recorded in the three most widespread lineages after 0.06 Ma. Asymmetric gene flow was recorded in the north Eurasian clade from southwestern Europe eastward, whereas the east Asian clade was rooted in south central China. Our results indicate that the fragmentation of the six clades of Pica was related to climatic cooling and aridification during periods of the Pliocene–Pleistocene. Populations on both sides of the Eurasian continent were similarly influenced by the Pleistocene climate changes and expanded concomitantly with the expansion of steppes. Based on results we also propose a revised taxonomy recognising seven species of Pica.     

Phylogeography of silver‐studded blue,Plebejus subsolanus (Lepidoptera: Lycaenidae), in the central mountainous regions of Japan

The silver‐studded blue, Plebejus subsolanus, is widely distributed in the Russian Altai mountains, northeastern China, the Korean Peninsula, and the Japanese archipelago. In Japan, the species is distributed across wide elevation ranges from the lowlands of Hokkaido to the subalpine zone of Honshu. Current subspecies classification in Japan is as follows: ssp. iburiensis, occurring in lowland grasslands in Hokkaido; ssp. yaginus in lower mountain grasslands in Honshu; and ssp. yarigadakeanus in higher mountain grasslands in Honshu. The habitat of this species has been markedly reduced due to recent habitat destruction and land‐use changes. Here, we undertook phylogeographic analyses of two subspecies, ssp. yaginus and yarigadakeanus in the central mountainous regions of Japan, based on two mitochondrial gene sequences, in order to collect information for establishing effective conservation strategies. From 57 samples from the four mountain ranges, we obtained a haplotype network comprised of 12 haplotypes. Because of the haplotype network topology, the geographic distribution of haplotypes and the correspondence of haplotype divergence to subspecies taxonomy, we provisionally divided the haplotypes into three haplogroups: YR1 and YR2, which comprised ssp. yarigadakeanus, and YG, which comprised ssp. yaginus. Mitochondrial DNA genetic differentiation generally agreed with morphological subspecies classification. The haplotype network suggested that ssp. yarigadakeanus populations had multiple origins, and the subspecies character of “bright blue of the male's wings” was assumed to have evolved independently in each subalpine meadow. We found that P. subsolanus was genetically differentiated depending upon the elevation at each mountain region, suggesting that each haplogroup should be a conservation unit.     

Phylogeography of Arabis serrata (Brassicaceae) in the Japanese Archipelago based on chloroplast DNA variations

Arabis serrata (Brassicaceae), a perennial plant widely distributed along the Japanese Archipelago, occurs in various habitats: for example, limestone zones, serpentine barrens, volcanic soils, and roadsides. It likely survived by adapting to its surrounding environment, resulting in great morphological and ecological variation. In this study, we performed a phylogeographic analysis to examine past changes in the distribution of A. serrata following climate oscillations during the Pleistocene. To cover most of A. serrata's range, leaves were collected from eight to ten individuals randomly selected from each of 37 populations in the Japanese Archipelago. Two chloroplast noncoding regions of the samples were amplified and sequenced: trnT(GGU)‐psbD and trnH(GUG)‐psbA spacers. Twenty‐five haplotypes were detected and distinguished by 31 substitutions. Four main haplotypes were observed in many populations distributed throughout the Japanese Archipelago. According to the genetic boundaries detected using the Monmonier algorithm, A. serrata is clustered into four groups, each including several populations: Hokkaido Island, northern mainland Honshu, central Japan, and western Japan. The boundaries, however, were not robust because all genetic parameters did not support the differentiation among groups. These results indicate the absence of an obvious geographic structure in the distribution of A. serrata, suggesting that this species has experienced a rapid range expansion in postglacial times.     

Colonization history of the Japanese water shrewChimarrogale platycephala, in the Japanese Islands

To evaluate the intraspecific evolutionary history and local differentiation of the Japanese water shrewChimarrogale platycephala (Temminck, 1842), we an a lyzed the mitochondrial cytochromeb (Cytb) sequence divergence for samples from 55 localities in the Japanese is lands of Honshu and Kyushu. According to phylogenetic trees based on theCytb data, there were fourCytb haplotype lineages, which showed rough affinities with geographic areas, namely, Eastern/Central Honshu, the Kinki District of Western Honshu, the Chugoku District of Western Honshu, and Kyushu. However, in the alpine areas of the boundary between the Kinki and Chugoku Districts, complicated distribution patterns of theCytb haplotypes were revealed. Considering the present data and geological history in the Quaternary, we hypothesized the following evolutionary scenario. First, differ entiation and division into four primary ancestral geographic colonies of the shrews occurred in hypothetical refugia in the mid — late Pleistocene. Subsequently, rapid expansion occurred and caused the complicated distribution patterns of theCytb haplotypes in the boundary areas owing to the complex topography during the late stage of the Quaternary.     

Genetic Variation and Population Structure of Hair Crab (<Emphasis Type="Italic">Erimacrus isenbeckii</Emphasis> ) in Japan Inferred from Mitochondrial DNA Sequence Analysis

Genetic variation and population structure of hair crab (Erimacrus isenbeckii) were examined using nucleotide sequence analysis of 580 base pairs (bp) in the 3′ portion of the mitochondrial cytochrome c oxidase subunit I gene (COI) of 20 samples collected from 16 locales in Japan (the Hokkaido and Honshu Islands) and one in Korea. A total of 27 haplotypes was defined by 23 variable nucleotide sites in the examined COI region. Pairwise population F _ST estimates and neighbor-joining tree inferred distinct genetic differentiation between the representative samples from the Pacific Ocean off the Eastern Hokkaido Island and the Sea of Japan, while others were intermediate between these two groups. AMOVA also showed a weak but significant differentiation among these three groups. The present results suggest a moderate population structure of hair crab, probably influenced by high gene flow between regional populations due to sea current dependent larval dispersal of this species.     

Differential geographic patterns of mitochondrial DNA variation in two sympatric species of Japanese wood mice, Apodemus speciosus and A. argenteus

We examined the gene sequences of mitochondrial cytochrome b (cyt b) in two Japanese wood mouse species, Apodemus speciosus (n = 89) and A. argenteus (n = 46), which are distributed on the four main islands of Japan (Hokkaido, Honshu, Shikoku, and Kyushu) and on the small islands surrounding them. Apodemus speciosus, the larger of the two species, showed substantial genetic variation, with a maximum of 3% sequence divergence, and remarkable phylogenetic subdivision with two major clades. One clade represents haplotypes from a central region, including Honshu, Shikoku, Kyushu, and their adjacent islands; the other clade includes haplotypes from Hokkaido and the peripheral islands, forming four subclades: a) Hokkaido, b) Sado Island, c) Satsunan Islands, and d) the Izu Islands. Sequence divergence among the four subclades was 1.0 to 1.5%, implying that A. speciosus colonized these geographic regions 0.2 to 0.3 million years ago, assuming a substitution rate of 2.4% per million years. The population on the Izu Islands has preserved haplotypes that are distinct from those in any other region, providing good evidence for the natural colonization of the volcanic islands of the Izu Islands. The cyt b sequence variation had no relation to the karyotypic dimorphism for the eastern (2n = 48) and western (2n = 46) geographic groups, between which a strict border exists at central Honshu. On the other hand, Apodemus argenteus, the smaller of the two species, showed a similar level of sequence divergence (maximum of 3%) but no substantial geographic differentiation: populations in Hokkaido, Sado, and Yakushima shared similar haplotypes with each of the central populations, suggesting that genetic exchanges occurred among the localities in the last 0.15 million years. The apparent genetic structure of the mitochondrial DNA found in the A. speciosus population might be caused solely by long-term existence in insular regions, presumably due to ecological superiority relative to A. argenteus.     

Genetic variations and phylogeography of the swallowtail butterfly Papilio machaon on the Japanese Islands

The swallowtail butterfly Papilio machaon Linnaeus, 1758 is widely distributed in the Holarctic region, including all of the main islands of Japan, as well as Sakhalin, and on other smaller islands south to Yakushima Island. The Japanese population is situated at the margin of the Eurasian distribution range of this species. It is morphologically different from other populations and has been classified as the subspecies hippocrates C. & R. Felder, 1864. The population of the Japanese Islands is considered to be genetically distinct from the continental populations in relation to the geographical history of the Japanese Islands. Therefore, we examined a part of the ND5 gene sequence of the mitochondrial DNA for P. machaon individuals of various localities in Japan and some nearby countries, and found 68 haplotypes in 400 individuals from the Japanese Islands and Sakhalin. A DNA polymorphism analysis revealed that the genetic structure of the Hokkaido population was significantly different from that of the southern populations on the main Japanese islands. These results imply that P. machaon expanded its range from the Amur region of Russia southward through Sakhalin to the Japanese Islands, and that the Tsugaru Strait between Hokkaido and Honshu may have subsequently limited their gene flow as a geographical barrier.     

Population structure and gene flow among wild populations of the black muntjac (Muntiacus crinifrons) based on mitochondrial DNA control region sequences

The black muntjac (Muntiacus crinifrons) is a rare species endemic to China and its current distribution is confined to partial mountain ranges in eastern China. To assess the population structure and gene flow among the extant populations, we sequenced a portion of the mitochondrial control region (424 bp) for 47 samples collected from its current three large populations (Huangshan, Tiammushan and Suichang). A total of 18 unique haplotypes were defined based on 22 polymorphic sites. Of these haplotypes, 15 ones were population-specific and only one haplotype was shared among the three populations. Significant genetic differentiation was detected between Suichang and Huangshan populations (?(ST)=0.1677, P<0.001) or between Suichang and Tianmushan populations (?(ST)=0.2002, P<0.001), indicating that the Suichang population may be spatially structured from other two populations along matriline. The Mantel test revealed that this significant differentiation was not driven by geographic distance (P=0.14), implying that genetic divergence of Suichang population might result from recent human disturbances. Phylogenetic analyses suggested the mitochondrial control region haplotypes were split into two well divergent clades (Clade I and Clade II). Interestingly, the two distinct haplotype clades were found to coexist in Suichang area. The nested clade analysis revealed a significant phylogeographic structure among the black muntjac populations (total cladogram: chi2=18.68; P<0.001), which was inferred to result from past fragmentation followed by range expansion. The population expansion was supported by the analysis of mismatch distribution and the tests of neutrality. Therefore, we suggest that the coexistence of distinct haplotypes in Suichang population was induced by historical population expansion after fragmentation and that the current genetic differentiation should be attributed to the reduction of female-mediated gene flow due to recent habitat fragmentation and subsequent loss.     

Comparative historical biogeography of Plateumaris leaf beetles (Coleoptera: Chrysomelidae) in Japan: interplay between fossil and molecular data

Aim In an attempt to use molecular and fossil data interactively in historical biogeography, we studied the phylogeography of five Plateumaris leaf beetles in Japan using mitochondrial cytochrome oxidase subunit I (COI) sequence data to explore interspecific differences in phylogeographical patterns and estimate the timings of colonization and geographical differentiation. Location A total of 461 beetles from five species on Hokkaido, Honshu and Kyushu islands of Japan were analysed with 117 beetles from three conspecies and two congeners from the mainland (Russia, including Sakhalin; Korea; Mongolia; Belgium; France; hereafter, the continent). Methods Using the sequence data from a 750‐bp portion of the COI gene, we studied the phylogeny of COI haplotypes, intraspecific population differentiation using analysis of molecular variance and the Mantel test, and intraspecific phylogeography using nested clade analysis. In addition, divergence times between the continental and Japanese lineages, as well as among the various Japanese lineages, were estimated using a Bayesian approach with node constraints based on fossil records of extant species. Results Three widely distributed species showed different degrees of geographical differentiation corresponding to their different colonization history in Japan. Bayesian estimates of divergence time revealed that one of two endemic species, which originated before the late Pliocene, attained intraspecific differentiation through the Pliocene and Pleistocene, whereas another endemic species has been confined in one locality, and three non‐endemic species colonized Japan after the mid‐Pleistocene. Main conclusions Molecular analyses of an insect group with relatively abundant fossil data can contribute greatly to the understanding of diverse biogeographical histories of related species in a region. Bayesian estimates of divergence time could be used to assess the variable evolutionary rates of the COI gene, and may be applied to other related insect species.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

Phylogeography of the cold‐water barnacle Chthamalus challengeri in the north‐western Pacific: effect of past population expansion and contemporary gene flow

Aim Phylogeographical patterns of marine organisms in the north‐western Pacific are shaped by the interaction of past sea‐level fluctuations during glacial maxima and present‐day gene flow. This study examines whether observed population differentiation in the barnacle Chthamalus challengeri, which is endemic to the north‐western Pacific, can be explained by the interactions between historical glacial events and patterns of contemporary gene flow. Location Eleven locations in the north‐western Pacific. Methods Partial sequences of mitochondrial cytochrome c oxidase subunit I (COI), 12S, 16S and nuclear internal transcribed spacer 1 (ITS1) were obtained from 312 individuals. Parsimony haplotype networks and analysis of molecular variance (AMOVA) were used to determine whether the observed genetic structure corresponds to marine provinces (Kuroshio Current and China Sea Coastal Provinces), zoogeographical zones (oriental and Japan warm‐temperate zones) and/or potential refugial areas (Sea of Japan and East China Sea) in the north‐western Pacific. Neutrality tests, mismatch distribution analysis and Bayesian skyline plots were used to infer the demographic history of C. challengeri. Results In total, 312, 117, 182 and 250 sequences were obtained for COI, 12S, 16S and ITS1, respectively. A panmictic population was revealed, which did not conform to the ‘isolation by distance’ model. None of the a priori population groupings based on marine provinces, zoogeographical zones or potential refugial areas was associated with observed genetic patterns. Significant negative values from neutrality tests and the unimodal mismatch distribution and expansion patterns in Bayesian skyline plots for the COI and 16S data sets indicate a population expansion in the mid‐Pleistocene (c. 200 ka). Information from the fossil record suggests that there has been a northward range shift of this species from the East China Sea or the Palaeo‐Pacific coast of Japan to the Sea of Japan since the mid‐Pleistocene. Main conclusions Chthamalus challengeri has experienced a population expansion and range shift since the mid‐Pleistocene. The observed lack of population differentiation can be explained by this past population expansion and present‐day wide‐scale larval dispersal (owing to the long planktonic larval duration) across marine provinces, which have led to the successful establishment of the species in different zoogeographical zones and habitats.     

Timing of the population dynamics of bullhead Cottus gobio (Teleostei: Cottidae) during the Pleistocene

Abstract Pleistocene genetic structure of the bullhead, Cottus gobio, was evaluated across the western Palearctic using a 771‐bp long fragment of the mitochondrial control region in 123 individuals collected at 35 sites (data set I). In total, 59 haplotypes that differed at 73 positions (9.3%) were detected. Data analysis also included sequences from Englbrecht et al. (2000 ), thus increasing the sampling to a more comprehensive data set of 529 fish and 63 control region sequences of 482 bp (data set II). A minimum spanning and phylogenetic tree identified a seventh clade (Brittany–Loire) in addition to the previously identified six clades. The geographical range of the North Sea and Lower Rhine clades was considerably larger than thought previously. Haplotype diversity was generally low, and the total fixation index high (F_ST = 0.49). Among‐group differentiation accounted for 41.7% (data set I) of the variation. Contiguous range expansions and restricted gene flow combined with isolation by distance, interspersed with past fragmentation characterize bullhead across its range. New is the knowledge that dated interglacial periods correlated with population expansions; river captures, proglacial lake systems and sea level played a significant role in the dispersal and expansion either in northern or southern direction. Hence it became possible to identify and date the colonization routes and putative palaeorefugia, most of which were located in Central and North‐west Europe. Glacial periods resulted in distinct fragmentation events and lineage sorting.     

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.     

Phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands

The phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands was analyzed on the basis of nucleotide sequences of a mitochondrial gene for cytochrome oxidase c subunit I (COI). Extremely low genetic diversity was found in populations at both the northern and southern boundaries of the geographic distribution of this species in Japan, i.e., Sendai Bay and Iriomotejima Island, respectively, which might be attributed to the population bottleneck due to historical environmental variations and/or the recent foundation of populations in the marginal part of the inhabitable range. Most populations contained unique rare haplotypes, and significant genetic differentiation on the whole was shown, while no clear geographic genetic structure was detected between the Japanese and Ryukyu Islands or over the distribution area of B. zonalis in Japan, with the exception of significant genetic divergences in Ago Bay in the central part of Honshu and the southern part of Okinawajima Island.     

Shaped by uneven Pleistocene climate: mitochondrial phylogeographic pattern and population history of white wagtail Motacilla alba (Aves: Passeriformes)

We studied the phylogeography and population history of the white wagtail Motacilla alba, which has a vast breeding range, covering areas with different Pleistocene climatic histories. The mitochondrial NADH dehydrogenase subunit II gene (ND2) and Control Region (CR) were analyzed for 273 individuals from 45 localities. Our data comprised all nine subspecies of white wagtail. Four primary clades were inferred (M, N, SW and SE), with indications of M. grandis being nested within M. alba. The oldest split was between two haplotypes from the endemic Moroccan M. a. subpersonata (clade M) and the others, at 0.63–0.96 Mya; other divergences were at 0.31–0.38 Mya. The entire differentiation falls within the part of the Pleistocene characterized by Milankovitch cycles of large amplitudes and durations. Clade N was distributed across the northern Palearctic; clade SW in southwestern Asia plus the British Isles and was predicted by Ecological niche models (ENMs) to occur also in central and south Europe; and clade SE was distributed in central and east Asia. The deep divergence within M. a. subpersonata may reflect retention of ancestral haplotypes. Regional differences in historical climates have had different impacts on different populations: clade N expanded after the last glacial maximum (LGM), whereas milder Pleistocene climate of east Asia allowed clade SE a longer expansion time (since MIS 5); clade SW expanded over a similarly long time as clade SE, which is untypical for European species. ENMs supported these conclusions in that the northern part of the Eurasian continent was unsuitable during the LGM, whereas southern parts remained suitable. The recent divergences and poor structure in the mitochondrial tree contrasts strongly with the pronounced, well defined phenotypical differentiation, indicating extremely fast plumage divergence.