Phylogeny and phylogeography of the sword-tailed newt,Cynops ensicauda (Amphibia: Caudata), as revealed by nucleotide sequences of mitochondrial DNA

We investigated the phylogenetic relationships and phylogeography among individuals of the endemic newt (Cynops ensicauda) of the Central Ryukyus, Japan based on the mitochondrial cytochrome b gene. The results of phylogenetic analyses showed the presence of remarkable differentiation between assemblages from the Amami and Okinawa Island Groups, supporting the validity of the subspecies C. e. popei described from the latter on morphological grounds. The group of individuals from Okinawa Island Group was further divided into two distinct subgroups unlike the results of previous morphological and allozyme studies. Geographic ranges of these subgroups overlapped in the northern part of Okinawajima Island. The phylogeographic pattern within the Okinawa Island Group suggests an initial division into two geographically isolated population lineages and subsequent secondary sympatry before formation of reproductive isolation. It is also likely that within the Okinawa Island Group emigration occurred from the central and northern parts of Okinawajima Island to its southern part, as well as to several small islets off its western coast. Within the Amami Island Group, recurring restricted gene flow with isolation by distance seems to have occurred around the southwestern region including three islets southwest of Amamioshima Island.     

Taxonomic study of the Burmoniscus ocellatus complex (Crustacea, Isopoda, Oniscidea) in Japan shows genetic diversification in the southern Ryukyus, southwestern Japan

To clarify the taxonomic status of the Burmoniscus ocellatus complex in Japan, we carried out morphological observations and phylogenetic analyses of specimens collected from Yonagunijima, Iriomotejima, Ishigakijima, and Miyakojima Islands of the southern Ryukyus and from Okinawajima Island of the central Ryukyus in southwestern Japan. Observations of the holotypes of Aphiloscia iriomotensis ( Nunomura, 1986 ), Ap. ishigakiensis ( Nunomura, 1986 ), and Ap. yonakuniensis ( Nunomura, 1986 ), in addition to the specimens newly collected from the five islands, indicated that these specimens belong to the genus Burmoniscus. Analyses of five morphological characters of 268 specimens collected from the five islands showed that the body size of Okinawajima specimens is distinctly smaller than those of the specimens from the southern Ryukyus. The ranges of the five morphological characters tended to overlap among the specimens from Yonagunijima, Iriomotejima, Ishigakijima and Miyakojima Islands; these morphological characters were congruent with those of 6. ocellatus (Verhoeff, 1928). The phylogenetic analyses were based on three regions of mitochondrial DNA-COI, 12S rRNA, and 16S rRNA-and indicated that the specimens collected from the southern Ryukyus constitute a monophyletic group, which is clearly distinct from the clade composed of the Okinawajima specimens. These results strongly suggest that Ap. iriomotensis, Ap. ishigakiensis, and Ap. yonakuniensis are synonymous with B. ocellatus, a species widely distributed in the southern Ryukyus. On the other hand, the species from Okinawajima Island in the central Ryukyus is considered to be an undescribed Burmoniscus species, which is closely related to but clearly distinct from S. ocellatus.     

Glacial refugia in a maritime temperate climate: Cicada (Kikihia subalpina) mtDNA phylogeography in New Zealand

Understanding the biological significance of Pleistocene glaciations requires knowledge of the nature and extent of habitat refugia during glacial maxima. An opportunity to examine evidence of glacial forest refugia in a maritime, Southern Hemisphere setting is found in New Zealand, where the extent of Pleistocene forests remains controversial. We used the mitochondrial phylogeography of a forest-edge cicada ( Kikihia subalpina ) to test the hypothesis that populations of this species survived throughout South Island during the Last Glacial Maximum. We also compared mitochondrial DNA phylogeographic patterns with male song patterns that suggest allopatric divergence across Cook Strait. Cytochrome oxidase I and II sequences were analyzed using network analysis, maximum-likelihood phylogenetic estimation, Bayesian dating and Bayesian skyline plots. K. subalpina haplotypes from North Island and South Island form monophyletic clades that are concordant with song patterns. Song divergence corresponds to approximately 2% genetic divergence, and Bayesian dating suggests that the North Island and South Island population-lineages became isolated around 761 000 years bp . Almost all South Island genetic variation is found in the north of the island, consistent with refugia in Marlborough Sounds, central Nelson and northwest Nelson. All central and southern South Island and Stewart Island haplotypes are extremely similar to northern South Island haplotypes, a 'northern richness/southern purity' pattern that mirrors genetic patterns observed in many Northern Hemisphere taxa. Proposed southern South Island forest habitat fragments may have been too small to sustain populations of K. subalpina , and/or they may have harboured ecological communities with no modern-day analogues.     

Phylogeography of the pademelons (Marsupialia: Macropodidae: Thylogale) in New Guinea reflects both geological and climatic events during the Plio‐Pleistocene

Aim Alternative hypotheses concerning genetic structuring of the widespread endemic New Guinean forest pademelons (Thylogale) based on current taxonomy and zoogeography (northern, southern and montane species groupings) and preliminary genetic findings (western and eastern regional groupings) are investigated using mitochondrial sequence data. We examine the relationship between the observed phylogeographical structure and known or inferred geological and historical environmental change during the late Tertiary and Quaternary. Location New Guinea and associated islands. Methods We used primarily museum specimen collections to sample representatives from Thylogale populations across New Guinea and three associated islands. Mitochondrial cytochrome b and control region sequence data were used to construct phylogenies and estimate the timing of population divergence. Results Phylogenetic analyses indicated subdivision of pademelons into ‘eastern’ and ‘western’ regional clades. This was largely due to the genetic distinctiveness of north‐eastern and eastern peninsula populations, as the ‘western’ clade included samples from the northern, southern and central regions of New Guinea. Two tested island groups were closely related to populations north of the Central Cordillera; low genetic differentiation of pademelon populations between north‐eastern New Guinea and islands of the Bismarck Archipelago is consistent with late Pleistocene human‐mediated translocations, while the Aru Islands population showed divergence consistent with cessation of gene flow in the mid Pleistocene. There was relatively limited genetic divergence between currently geographically isolated populations in subalpine and nearby mid‐montane or lowland regions. Main conclusions Phylogeographical structuring does not conform to zoogeographical expectations of a north/south division across the cordillera, nor to current species designations, for this generalist forest species complex. Instead, the observed genetic structuring of Thylogale populations has probably been influenced by geological changes and Pleistocene climatic changes, in particular the recent uplift of the north‐eastern Huon Peninsula and the lowering of tree lines during glacial periods. Low sea levels during glacial maxima also allowed gene flow between the continental Aru Island group and New Guinea. More work is needed, particularly multi‐taxon comparative studies, to further develop and test phylogeographical hypotheses in New Guinea.     

The biogeographical history of Asian keelback snakes of the genus Hebius (Squamata: Colubridae: Natricinae) in the Ryukyu Archipelago,Japan

The islands of the central and southern parts of the Ryukyu Archipelago and Taiwan accommodate suites of allopatric endemics for many terrestrial vertebrates. These subregions, which were isolated as three independent islands since the Late Miocene period, became connected to the continent by the Early Pleistocene land bridge (EPL). Many endemics are considered to have emerged through vicariant events accompanied with the collapse of EPL, but this is not well corroborated. We conducted divergence dating analyses for samples of four allopatric keelback snakes from 20 islands using a published timeframe for natricine snakes. Phylogenetic analyses recognised the existence of eight subclades for the four species. Four of the five subclades in the Central Ryukyus have emerged after the initial isolation of three subregions but before the EPL collapse, indicating iterative divergences within the super‐island and stagnation of diverged entities. In contrast, the Southern Ryukyus was subjected to multiple incursions from both Central Ryukyus and Taiwan, irrespective of the EPL connection. Consequently, among all seven nodes, only two were consistent with the assumption of vicariance. These findings suggest that dispersals and subsequent nonvicariant divergences have contributed to the formation of the current phylogeographical pattern in this continental island chain.     

Late Pleistocene lineage divergence among populations of Neolitsea sericea (Lauraceae) across a deep sea‐barrier in the Ryukyu Islands

Aim Our goals were: (1) to assess the population genetic structure and demographic divergence history of a bird‐dispersed tree, Neolitsea sericea, endemic to East Asian land‐bridge islands; and (2) to interpret the results in the light of controversies over the dating and configurations of land bridges through the Japanese Ryukyu Island Arc. Location Japan–Ryukyu–Taiwan Island Arc and Chinese/Korean offshore islands. Methods We applied 10 nuclear microsatellites (nSSRs) and one chloroplast (cp) DNA sequence marker (psbA–trnH intergenic spacer) to 31 populations (397 and 326 individuals, respectively) from throughout the species’ range to infer current patterns of genetic diversity and structure, and pollen‐to‐seed migration ratios (r). A coalescent‐based isolation‐with‐migration (IM) model was applied to the combined nSSR/cpDNA data set to estimate lineage divergence times and population demographic parameters. Results The geographic structure of nSSRs and the distribution of most cpDNA haplotypes revealed two distinct lineages located in areas north and south of the ‘Tokara Gap’, a narrow (c. 37 km wide) but deep (> 1000 m) sea‐strait between the northern and central Ryukyus. Based on the IM analyses, we (1) dated the divergence of these northern and southern lineages to c. 0.07 Ma (90% highest posterior density interval: 0.02–0.38 Ma); (2) estimated a slightly smaller effective population size for the northern compared to the southern lineage; and (3) recovered only trivial signals of post‐divergence gene flow between them. Main conclusions The estimated divergence time for northern and southern lineages is consistent with geological evidence for the existence of land connections in the Tokara region during cold stages of the latest Pleistocene; it is thus incompatible with an ‘ancient sea‐barrier hypothesis’ for the Ryukyu Arc, where we would have expected much older divergences related to the initial formation of the Tokara and Kerama tectonic straits during the Pliocene. Multiple factors are likely to have had a role in the divergence of N. sericea, including not only land‐bridge submergence, but also island configuration, and/or constraints on adaptation along a latitudinal temperature gradient.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Extensive genetic differentiation in Gobiomorphus breviceps from New Zealand

Partial mitochondrial DNA sequences for parts of the cytochrome b gene and control region were obtained for 89 upland bullies Gobiomorphus breviceps from 19 catchments in New Zealand. There were two highly distinctive mtDNA clades: a northern clade corresponding to the North Island, northern South Island and west coast South Island, and a south‐east clade, in the southern and eastern South Island. Within these major clades there were further distinct clades that correlated with geographic sub‐regions and catchments. The marked genetic differentiation has occurred in the absence of obvious morphological divergence. Based on cytochrome b sequence divergences and the molecular clock hypothesis, the northern and southeastern clades correspond with the uplift of the Southern Alps during the Pliocene, while populations in the North Island and northern South Island were estimated to have diverged during the Pleistocene. The widescale geographic divergences were similar to those observed in the galaxiids, Galaxias vulgaris and Galaxias divergens , but biogeographic management boundaries may not be the same, reflecting different evolutionary histories for non‐diadromous species occupying the same areas.     

Species delimitation and biogeography of the Ryukyu ground geckos,Goniurosaurus kuroiwae ssp. (Squamata: Eublepharidae), by use of mitochondrial and nuclear DNA analyses

In the Ryukyu Archipelago, Japan, Goniurosaurus geckos are currently divided into six allopatric taxa among nearby islands. Recent studies suggested the occurrence of large genetic divergence within a single island and the possible non‐monophyly of a few taxa, but their species delimitation is not well resolved. We investigated the taxonomic relationships between the possibly geographically overlapped, highly diverged entities as well as other island populations via dense sampling. Our mitochondrial and nuclear DNA analyses showed that the two genetic groups were distributed in parapatry within Okinawajima Island and that they were hybridizing in narrow area around the contact zone. Geohistorical evidence suggests that the restricted gene flow has been caused by some intrinsic isolation mechanisms. Thus, we conclude that the two lineages represent full species. Mitochondrial analysis also showed that the genetic differences among other island populations were comparable to those between these species, thereby suggesting the presence of seven full species, including one unnamed taxon. We also discuss the possible cause of this divergence and why it has occurred at such a fine geographic scale.     

Phylogeography and historical demography of the Italian treefrog, Hyla intermedia, reveals multiple refugia, population expansions and secondary contacts within peninsular Italy

We investigated the geographical patterns of genetic diversity in the Italian treefrog through sequence analysis of a mitochondrial cytochrome b gene fragment. Three main mitochondrial lineages were identified, distributed in northern, central and southern Italy, respectively. Their divergence appears indicative of a split time largely predating Late Pleistocene climatic oscillations, and syntopy between them was only observed in the geographically intermediate populations. The historical demographic reconstructions suggest that in both northern and central Italy, an expansion occurred during the last major glacial phase, when a vast widening of the lowland habitats followed the glaciation-induced fall of the sea level. Instead, in southern Italy an expansion event likely followed the end of the last glaciation, although the inference of expansion appears less reliable for the southern clade than for the others. Within this geographical area, a sharp phylogeographic discontinuity separated peninsular from Sicilian populations, and the overall pattern of diversity suggests that the latter derived from a recent colonization of the island, probably through a Late Pleistocene land bridge. Phylogenetic, phylogeographic and historical demographic analyses thus concur in delineating a scenario of multiple refugia, with four groups of populations which survived the last glacial-interglacial cycles in at least three distinct refugia arranged along peninsular Italy, and have recently come into contact following range expansions. Therefore, these results support the hypothesis that a plethora of microevolutionary processes, rather than the prolonged stability of populations, were mainly responsible for shaping the patterns of diversity within this major biodiversity hotspot.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

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.     

Geographical variation in echolocation call and body size of the Okinawan least horseshoe bat, Rhinolophus pumilus (Mammalia: Rhinolophidae), on Okinawa-jima Island, Ryukyu Archipelago, Japan

The Okinawan least horseshoe bat, Rhinolophus pumilus, is a cave-dwelling species endemic to the central and southern Ryukyus, Japan. We analyzed variation in the constant frequency (CF) of the echolocation call and in forearm length (FAL) of this species on Okinawa-jima Island on the basis of data for 479 individuals from 11 caves scattered over the island. CF values in samples from six caves, all located in the southwestern half of Okinawa-jima, were significantly higher than those in samples from five caves in the northeastern half of the island. Also, FAL was significantly greater in the latter group than in the former group, although the ranges of variation in this character substantially overlapped between the two groups. These results suggest substantial differentiation between R. pumilus populations on Okinawa-jima. The implications of our findings for the conservation of this endangered bat species are briefly discussed.     

Allozyme diversity and genetic divergence of the dolly varden Salvelinus malma Walbaum from the Kuril islands

Genetic variation was studied in the southern subspecies of the Asian Dolly Varden Salvelinus malma krascheninnikovi from the Kuril Islands. Thirty-six genetic loci controlling 19 enzyme systems were analyzed in 13 Dolly Varden populations from the Shumshu, Paramushir, Onekotan, Rasshua, Simushir, Urup, Iturup, and Kunashir islands. In the studied populations, the proportion of polymorphic loci was 35 to 85% and the mean heterozygosity was 0.104 to 0.173; populations from the Kunashir Island were characterized by maximum heterozygosity. In the island populations examined, significant inter-population heterogeneity of allele frequencies was found for all studied population pairs. For the total population of all islands, the inter-population diversity (GST = 0.188) was comparable to this parameter for the total population from the Kunashir Island (GST = 0.170). Genetic distances between populations did not correlate with the corresponding geographical distances, which indicates the lack of a pronounced gene exchange between the island populations. Cluster analysis and multidimensional scaling based on genetic distances did not reveal clear groups among the studied populations but indicated greater similarity within the Iturup-Simushir-Urup-Paramushir group and a greater genetic divergence of the Kunashir, Onekotan, Rasshua, and especially Shumshu populations. In the Shumshu population, allele frequencies indicate the admixture of genes of the northern Dolly Varden. The observed pattern of genetic differentiation was probably caused largely by genetic drift under the conditions of a limited gene flow because of homing (which is typical of the Dolly Varden) and the presence of isolated nonanadromous populations. The population-genetic analysis of the Dolly Varden from the Kuril Islands does not give grounds to distinguish any other isolated Dolly Varden species in this region than S. malma, which is represented by the southern form S. m. krascheninnikovi with an admixture of the northern form S. m. malma in the Shumshu Island.     

Genetic variation and population differentiation of Michelia formosana (Magnoliaceae) based on cpDNA variation and RAPD fingerprints: relevance to post-Pleistocene recolonization

We used sequence variation of the atpB- rbcL intergenic spacer of cpDNA and nested clade analysis to assess the phylogeographic pattern of Michelia formosana, a species restricted to Taiwan and the Ryukyus. In total, 31 haplotypes were identified and clustered into four major chlorotypes. Genetic composition of nearly all populations was heterogeneous and paraphyletic phylogenetically. Although the apportionment of cpDNA variation hardly revealed a geographic pattern due to the coancestry of dominant sequences, some chlorotypes were restrictedly distributed. According to the patterns of clade dispersion and displacement, a reconstructed minimum spanning network revealed that historical events of past fragmentation and range expansion, associated with glaciation, may have shaped the phylogeographic patterns of M. formosana. Four possible refugia were identified: the Iriomote and Ishigaki Islands (the southern Ryukyus), Wulai (northern Taiwan), and Nanjen (southern Taiwan), on the basis of the interior positions of their haplotypes in the network and the high level of nucleotide diversity. Given insufficient time for coalescence at the cpDNA locus since the late Pleistocene recolonization, lineage sorting led to low levels of genetic differentiation among populations. In contrast, hierarchical examination of the random amplified polymorphic DNA (RAPD) data scored from six populations across three geographical regions, using an analysis of molecular variance (AMOVA), indicated high genetic differentiation both among populations (Phi(ST) = 0.471) and among regions (Phi(CT) = 0.368). An unweighted pair group method with arithmetic mean (UPGMA) tree of the RAPD fingerprints revealed that populations of two offshore islands of eastern Taiwan ( M. formosana var. kotoensis) were clustered with geographically remote populations of the Ryukyus instead of those in southern Taiwan, suggesting some historical division due to geographic barriers of the central mountain range. In contrast to the paraphyly of the nearly neutral cpDNA alleles, differentiated RAPDs may have experienced diversifying selection.     

Differentiation of golden-ringed dragonfly Anotogaster sieboldii (Selys, 1854) (Cordulegastridae: Odonata) in the insular East Asia revealed by the mitochondrial gene genealogy with taxonomic implications

Molecular phylogeographical analyses of Anotogaster sieboldii (Selys, 1854) were conducted to reveal the differentiation process of insular populations. The gene genealogy based on 845 bp of the mitochondrial genes (cytochrome oxidase subunit I and subunit II) indicated that A. sieboldii includes two deeply separated lineages. These two major lineages seem to have differentiated in Miocene before the formation of the insular East Asia. One lineage includes three inner clades that correspond to the populations of northern area (the Japanese main islands, Korean Peninsula, Yakushima), Amamioshima and Okinawajima. Populations of Central Ryukyu, including Amamioshima and Okinawajima, might have been divided from the northern populations in early Pleistocene. The other major lineage includes populations of the Yaeyama Group, Taiwan and East China. The former two populations were reconstructed as a reciprocal monophyletic group. Populations of Taiwan and Yaeyama Groups would have been separated from the continental ones in Pleistocene. These two highly divergent lineages should be recognized as distinct species. Furthermore, the mitochondrial lineages revealed six genetically distinct and geographically isolated assemblages: (1) northern populations, (2) Amamioshima, (3) Okinawajima, (4) Yaeyama Group, (5) Taiwan and (6) East China.     

ALLOZYME DIFFERENTIATION AND GENETIC STRUCTURE IN ISLAND AND MAINLAND JAPANESE POPULATIONS OF CAMPANULA PUNCTATA (CAMPANULACEAE)

Allozyme variation was investigated in 17 Japanese populations of Campanula punctata, ten from the Izu Islands and seven in the mainland (Honshu). The data indicate that there are two groups, a mainland group and an island one, and that the systematically problematic Oshima Island (northernmost Izu island) populations are closely related to those of the other islands. Nei's genetic identity values among islands and among mainland populations were 0.95 and 0.97, respectively, while the value between island and mainland populations was 0.84, suggesting that the island populations are an independent species. Total genetic variation was nearly the same among island and mainland populations. However, the apportionment of variation within and among populations was considerably different; 14% of gene diversity exists among mainland populations, while 31% of the diversity exists among island populations. Mean outcrossing rates of self-incompatible mainland and Oshima populations are 0.62–0.79, supporting xenogamy; those in self-compatible island populations are 0.37–0.57 in the northern islands, indicating a mixed mating system, and 0.16–0.25 in southern ones, indicating dominant inbreeding. Total genetic diversity in each island population decreased with distance from the mainland. Genetic and geological data suggest that the ancestors of insular populations were founded on northern islands in a relatively ancient period and that they dispersed progressively to the southern ones. Chromosome number (2n = 34) and isozyme numbers indicate gene duplications in this species, which suggests it is an ancient polyploid.     

The shaping of mitochondrial DNA phylogeographic patterns of the brown hare (Lepus europaeus) under the combined influence of Late Pleistocene climatic fluctuations and anthropogenic translocations

The phylogeographic structure of the brown hare (Lepus europaeus) was studied by analysing mtDNA control region sequences of 98 individuals from continental and insular Greece, Bulgaria, Cyprus and northern Israel, together with 44 published sequences from Italy and central Europe. We found two distinct clades separated by an average nucleotide divergence of 6.6%, which may correspond to a Balkan and to an Asia Minor refugium. The estimated time of separation of the two clades was dated back to 105,000- 490,000 years ago. These two clades coexist in the area of northeastern Greece and Bulgaria, most likely as a result of a post-glacial northward expansion. Within the southern Balkan refugium, network analyses showed geographical structuring, which supports the hypothesis of several isolated Late Pleistocene populations. The central European and Italian populations appear to have originated from a non-detected northern Balkan population that was genetically closely related to some northern Greek populations, as a result of postglacial expansion, translocations or a combination of both. Moreover, several cases of ancient and recent translocations by humans were detected, especially for some island populations, while the eastern Aegean islands off the Asia Minor coast were most likely colonized naturally through Late Pleistocene land bridge connection. The genetic analysis presented here provides a framework for designing proper conservation and management guidelines for this species.     

Allozyme Diversity and Genetic Divergence of the Dolly Varden Salvelinus malma Walbaum from the Kuril Islands

Genetic variation was studied in the southern subspecies of the Asian Dolly Varden Salvelinus malma krascheninnikovi from the Kuril Islands. Thirty-six genetic loci controlling 19 enzyme systems were analyzed in 13 Dolly Varden populations from the Shumshu, Paramushir, Onekotan, Rasshua, Simushir, Urup, Iturup, and Kunashir islands. In the studied populations, the proportion of polymorphic loci was 35 to 85% and the mean heterozygosity was 0.104 to 0.173; populations from the Kunashir Island were characterized by maximum heterozygosity. In the island populations examined, significant interpopulation heterogeneity of allele frequencies was found for all studied population pairs. For the all island populations, the interpopulation diversity (G _ST = 0.188) was comparable to this parameter for the populations from the Kunashir Island only (G _ST = 0.170). Genetic distances between populations did not correlate with the corresponding geographical distances, which indicates the lack of a pronounced gene exchange between the island populations. Cluster analysis and multidimensional scaling based on genetic distances did not reveal clear groups among the studied populations but indicated greater similarity within the Iturup–Simushir–Urup–Paramushir group and a greater genetic divergence of the Kunashir, Onekotan, Rasshua, and especially Shumshu populations. In the Shumshu population, allele frequencies indicate the admixture of genes of the northern Dolly Varden. The observed pattern of genetic differentiation was probably caused largely by genetic drift under the conditions of a limited gene flow because of homing (which is typical of the Dolly Varden) and the presence of isolated nonanadromous populations. The population–genetic analysis of the Dolly Varden from the Kuril Islands does not give grounds to distinguish any other isolated char species in this region than S. malma, which is represented by the southern form S. m. krascheninnikovi with an admixture of the northern form S. m. malma in the Shumshu Island.     

Phylogeography of a Tertiary relict plant, Meconopsis cambrica (Papaveraceae), implies the existence of northern refugia for a temperate herb

The perennial herb Meconopsis cambrica, a western European endemic, is the only European species of the otherwise Himalayan genus Meconopsis and has been interpreted as a Tertiary relict species. Using rbcL and ITS sequence variation, we date the split between M. cambrica and its sister clade Papaver s.str. to the Middle to Upper Miocene (12.8 Myr, 6.4–19.2 Myr HPD). Within M. cambrica, cpDNA sequence variation reveals the existence of two groups of populations with a comparable level of genetic variation: a northern group from Great Britain, the Massif Central, the western Pyrenees and the Iberian System, and a southern group from the central and eastern Pyrenees. Populations from the Cantabrian Mountains were placed in both groups. Based on ITS sequence variation, the divergence between these two groups can be dated to 1.5 Myr (0.4–2.8 Myr HPD), and the age of the British populations is estimated as 0.37 Myr (0.0–0.9 Myr HPD). Amplified fragment length polymorphism results confirm the distinctive nature of the populations from Britain, the Massif Central and the central and eastern Pyrenees. These patterns of latitudinal variation of M. cambrica differ from patterns of longitudinal differentiation found in many other temperate species and imply glacial survival of the northern populations in northerly refugia. The primary differentiation into northern and southern cpDNA groups dates to near the onset of the Quaternary and suggests that an ancient phylogeographic pattern has survived through several glacial periods. Our data provide evidence that the species has persisted for a long period with a highly fragmented and probably very localized distribution.