Phylogeographic evidence for two mesic refugia in a biodiversity hotspot

Phylogeographic studies of flora in species-rich south-western Australia point to complex evolutionary histories, reflecting patterns of persistence and resilience to climatic changes during the Pleistocene. We asked whether coastal areas of the mid-west and south, as well as granite outcrops and inland ranges, have acted as major refugia within this region during Pleistocene climatic fluctuations by analysing phylogeographic patterns in the shrub Calothamnus quadrifidus R.Br. (Myrtaceae). We determined variation in chloroplast DNA data for 41 populations across the geographic range. Relationships and major clades were resolved using parsimony and Bayesian analyses. We tested for demographic and spatial expansion of the major clades and estimated clade divergence dates using an uncorrelated, lognormal relaxed clock based on two conservative chloroplast mutation rates. Two distinct phylogeographic clades were identified showing divergence during the Pleistocene, consistent with other phylogeographic studies of south-west Australian flora, emphasising the impact of climatic oscillations in driving divergence in this landscape. The southern clade was more diverse, having higher haplotype diversity and greater genetic structure, while the northern clade showed evidence of fluctuation in population size. Regions of high haplotype diversity with adjacent areas of low diversity observed in each clade indicated the locations of two coastal refugia: one on the south coast and another along the mid-west coast. This is the first evidence for major Pleistocene refugia using chloroplast genetic data in a common, widespread species from this region.     

Genetic Diversity and Geographic Distribution of Genetically Distinct Rabies Viruses in the Philippines

Background

Rabies continues to be a major public health problem in the Philippines, where 200–300 human cases were reported annually between 2001 and 2011. Understanding the phylogeography of rabies viruses is important for establishing a more effective and feasible control strategy.

Methods

We performed a molecular analysis of rabies viruses in the Philippines using rabied animal brain samples. The samples were collected from 11 of 17 regions, which covered three island groups (Luzon, Visayas, and Mindanao). Partial nucleoprotein (N) gene sequencing was performed on 57 samples and complete glycoprotein (G) gene sequencing was performed on 235 samples collected between 2004 and 2010.

Results

The Philippine strains of rabies viruses were included in a distinct phylogenetic cluster, previously named Asian 2b, which appeared to have diverged from the Chinese strain named Asian 2a. The Philippine strains were further divided into three major clades, which were found exclusively in different island groups: clades L, V, and M in Luzon, Visayas, and Mindanao, respectively. Clade L was subdivided into nine subclades (L1–L9) and clade V was subdivided into two subclades (V1 and V2). With a few exceptions, most strains in each subclade were distributed in specific geographic areas. There were also four strains that were divided into two genogroups but were not classified into any of the three major clades, and all four strains were found in the island group of Luzon.

Conclusion

We detected three major clades and two distinct genogroups of rabies viruses in the Philippines. Our data suggest that viruses of each clade and subclade evolved independently in each area without frequent introduction into other areas. An important implication of these data is that geographically targeted dog vaccination using the island group approach may effectively control rabies in the Philippines.     

Gene flow and population history in high dispersal marine invertebrates: mitochondrial DNA analysis of Holothuria nobilis (Echinodermata: Holothuroidea) populations from the Indo-Pacific

The sea cucumber, Holothuria nobilis, has a long-lived planktotrophic larvae, and previous allozyme surveys have suggested that high dispersal is realized. In contrast, recent ecological studies indicate that dispersal is low. To reconcile these data, and to investigate the evolution of this Indo-Pacific species, we screened geographical variation in 559 bp of a mitochondrial gene (COI) in 360 samples from the Australasian region and La Réunion. Sequences from La Réunion differed by > 7% from others and may constitute another species. Haplotype diversity in other samples was high (0.942, SD = 0.007), but haplotypes were closely related (mean nucleotide diversity: 0.0075, SD = 0.0041). AMOVA, pairwise FST values and exact tests did not detect significant population structure. Nested clade analysis showed that one of two main clades was over-represented in west Australia, whereas the other was more common in the northern Great Barrier Reef. Isolation-by-distance was identified as the main determinant of population structure at several clade levels. Contiguous range expansion was inferred for evolutionary older clade levels and this may correspond to a late Pleistocene (88 000-193 000 years ago) population expansion inferred from haplotype mismatch distributions. Thus, the population genetic structures detected are likely to be formed prior to the last ice age, with some indications for high dispersal on shorter time scales.     

A trend of central versus peripheral structuring in mitochondrial and nuclear gene sequences of the Japanese wood mouse, Apodemus speciosus

A phylogeographic analysis was performed on Japanese endemic wood mice (Apodemus speciosus) using nuclear interphotoreceptor retinol binding protein (IRBP) gene sequences (1,152 bp), together with previously published mitochondrial cytochrome b (cyt b) data. In the IRBP analysis, 40 haplotypes were recovered from 84 individuals by statistical and subcloning methods. Substantial sequence variation was determined from the IRBP data (pi=0.0047), and no significant evidence of recombination was detected. From the phylogenetic analysis, the 40 haplotypes fell into two major groups with geographic associations, irrespective of the karyotype groups (2n=46 and 2n=48), yielding a trend of central (Hokkaido, Honshu, Shikoku, Kyushu, and Sado) and peripheral (Izu, Oki, Tsushima, and Satsunan Is.) groupings. This geographic pattern is similar to that observed in the cyt b data, with a different insular grouping of Sado, Hokkaido, Izu, and Satsunan islands, and also to that of morphological features. In both gene data sets, nested clade analyses revealed allopatric fragmentation in the "peripheral island clades" and range expansion in the "central island clades." A mismatch analysis using cyt b data also suggested expansion of the central islands clade. Thus, the trend of central vs. peripheral structuring may be attributable to past demographic dynamics in the two distinct haplotype clades, such as range expansion of one clade in the central area of the Japanese Islands, leaving the other clade in the periphery.     

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.     

Recent volcanism and mitochondrial DNA structuring in the lizard Gallotia atlantica from the island of Lanzarote

The phylogeography of the lacertid lizard Gallotia atlantica from the small volcanic island of Lanzarote (Canary Islands) was analysed based on 1075 bp of mitochondrial DNA (mtDNA) sequence (partial cytochrome b and ND2) for 157 individuals from 27 sites (including three sites from neighbouring islets). Levels of sequence divergence were generally low, with the most distant haplotypes separated by only 14 mutational steps. MtDNA divergence appears to coincide with formation of the middle Pleistocene lowland that united formerly separate ancient islands to form the current island of Lanzarote, allowing rejection of a two-island model of phylogeographical structure. There was evidence of large-scale population expansion after island unification, consistent with the colonization of new areas. A nested clade phylogeographical analysis (NCPA) revealed significant phylogeographical structuring. Two-step and higher-level clades each had disjunct distributions, being found to the east and west of a common area with a north-south orientation that extends between coasts in the centre-east of the island (El Jable). Other clades were almost entirely restricted to the El Jable region alone. Bayesian Markov chain Monte Carlo analyses were used to separate ongoing gene flow from historical associations. These supported the NCPA by indicating recent (75,000-150,000 years ago) east-west vicariance across the El Jable region. Lava flows covered El Jable and other parts of the central lowland at this time and likely led to population extinctions and temporary dispersal barriers, although present-day evidence suggests some populations would have survived in small refugia. Expansion of the latter appears to explain the presence of a clade located between the eastern and western components of the disjunct clades. Direct relationships between mtDNA lineages and morphology were not found, although one of two morphological forms on the island has a disjunct distribution that is broadly concordant with east-west components of the phylogeographical pattern. This work demonstrates how recent volcanic activity can cause population fragmentation and thus shape genetic diversity on microgeographical scales.     

Genetic Differentiation in Insular Lowland Rainforests: Insights from Historical Demographic Patterns in Philippine Birds

Phylogeographic studies of Philippine birds support that deep genetic structure occurs across continuous lowland forests within islands, despite the lack of obvious contemporary isolation mechanisms. To examine the pattern and tempo of diversification within Philippine island forests, and test if common mechanisms are responsible for observed differentiation, we focused on three co-distributed lowland bird taxa endemic to Greater Luzon and Greater Negros-Panay: Blue-headed Fantail (Rhipidura cyaniceps), White-browed Shama (Copsychus luzoniensis), and Lemon-throated Leaf-Warbler (Phylloscopus cebuensis). Each species has two described subspecies within Greater Luzon, and a single described subspecies on Greater Negros/Panay. Each of the three focal species showed a common geographic pattern of two monophyletic groups in Greater Luzon sister to a third monophyletic group found in Greater Negros-Panay, suggesting that common or similar biogeographic processes may have produced similar distributions. However, studied species displayed variable levels of mitochondrial DNA differentiation between clades, and genetic differentiation within Luzon was not necessarily concordant with described subspecies boundaries. Population genetic parameters for the three species suggested both rapid population growth from small numbers and geographic expansion across Luzon Island. Estimates of the timing of population expansion further supported that these events occurred asynchronously throughout the Pleistocene in the focal species, demanding particular explanations for differentiation, and support that co-distribution may be secondarily congruent.     

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.     

Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae)

We sequenced 1077 bp of the mitochondrial cytochrome b gene and 511 bp of the nuclear Apolipoprotein B gene in bicoloured shrew (Crocidura leucodon, Soricidae) populations ranging from France to Georgia. The aims of the study were to identify the main genetic clades within this species and the influence of Pleistocene climatic variations on the respective clades. The mitochondrial analyses revealed a European clade distributed from France eastwards to north-western Turkey and a Near East clade distributed from Georgia to Romania; the two clades separated during the Middle Pleistocene. We clearly identified a population expansion after a bottleneck for the European clade based on mitochondrial and nuclear sequencing data; this expansion was not observed for the eastern clade. We hypothesize that the western population was confined to a small Italo-Balkanic refugium, whereas the eastern population subsisted in several refugia along the southern coast of the Black Sea.     

Trans-Arctic dispersals and the evolution of a circumpolar marine fish species complex, the capelin (Mallotus villosus)

Trans-Arctic dispersals and population and range expansions during the Pleistocene enhanced opportunities for evolutionary diversification and contributed to the process of speciation within the capelin, a northern marine-fish complex exhibiting a circumpolar distribution. Capelin is composed of four highly divergent and geographically discrete mitochondrial DNA (mtDNA) clades (609 bp; cytochrome b). Two clades occur in the North Atlantic, one associated with Canadian Atlantic waters, including Hudson Bay, and the second distributed from West Greenland to the Barents Sea. Two additional clades occur in the Arctic and northeast Pacific Oceans, representing the most recent divergence within the capelin phylogenetic tree. Judged from mtDNA diversity, capelin populations comprising all clades experienced at least one demographic and spatial reduction-expansion episode during recent Pleistocene glaciations that imprinted their molecular architecture. The large contemporary populations in the northeast Pacific and Arctic Oceans exhibited significant genetic structure whereas no such structure was detected in the equally extensive North Atlantic clades. All clades are characterized by one or two prevalent mtDNA haplotypes distributed over the entire range of the clade. Assuming a Pacific ancestor for capelin, we infer that capelin dispersed on two separate occasions to the North Atlantic. A more recent event resulted in the isolation of eastern Pacific and Arctic clades, with the Arctic clade positioned for a potential third Atlantic invasion, as revealed by the presence of this clade in the Labrador Sea. The Labrador Sea is a potential contact zone for three of the four capelin clades.     

Phylogeographic heterogeneity of the brown macroalga Sargassum horneri (Fucaceae) in the northwestern Pacific in relation to late Pleistocene glaciation and tectonic configurations

Pleistocene glacial oscillations and associated tectonic processes are believed to have influenced the historical abundances and distribution of organisms in the Asia Northwest Pacific (ANP). Accumulating evidence indicates that factors shaping tempospatial population dynamics and distribution patterns of marine taxa vary with biogeographical latitude, pelagic behaviour and oceanographic regimes. To detect what kinds of historical and contemporary factors affected genetic connectivity, phylogeographic profiles of littoral macroalga Sargassum horneri in the ANP were analysed based on mitochondrial (Cox3) and chloroplast (rbcL) data sets. Five distinct clades were recovered. A strong signature of biogeographical structure was revealed (Φ(CT) = 0.487, P < 0.0001) derived from remarkable differentiation in clade distribution, as clade I is restricted to Chinese marginal seas (Yellow-Bohai Sea, East China Sea and South China Sea), whereas clades II-V are discontinuously scattered around the main Islands of Japan. Furthermore, two secondary contact regions were identified along the south Japan-Pacific coastline. This significant differentiation between the two basins may reflect historical glacial isolation in the northwestern Pacific, which is congruent with the estimates of clade divergence and demographic expansion during the late Quaternary low sea levels. Analysis of molecular variance and the population-pair statistic F(ST) also revealed significant genetic structural differences between Chinese marginal seas and the Japanese basin. This exceptional phylogeographic architecture in S. horneri, initially shaped by historical geographic isolation during the late Pleistocene ice age and physical biogeographical barriers, can be complicated by oceanographic regimes (ocean surface currents) and relocating behaviour such as oceanic drifting.     

Phylogeography of Pseudacris regilla (Anura: Hylidae) in western North America, with a proposal for a new taxonomic rearrangement

The Baja California populations of Pseudacris regilla, a widespread species in Western North America ranging from British Columbia to southern Baja California, are characterized by extensive geographic fragmentation. We performed phylogeographic and historical demographic analyses on 609 bp of the cytochrome b mitochondrial gene of 110 individuals representing 28 populations to determine the relative influences of current and historical processes in shaping the present distribution of genetic diversity on the Baja California Peninsula. Haplotypes from this area were nested in a clade with three well-differentiated groups. Two of these groups are from Baja California Sur and another is from California and Baja California. The estimated date for the split of these groups, between 0.9-1 Ma, fits with previously proposed hypotheses of vicariance due to different transpeninsular seaways, although successive population fragmentation and expansion due to climatic oscillations during Pleistocene glaciations cannot be discarded. Historical demographic analyses detected signs of past population expansions, especially in the southernmost group. With respect to populations north of this region, two older clades were identified, one with haplotypes mainly distributed in central California, and the other corresponding to the northern half of the species range, in what apparently is a recurrent pattern in the Pacific coast of North America. Based on the concordance between mt-DNA and available allozyme data indicating that these species have a long independent evolutionary history, we propose to consider the three major clades as distinct species: P. regilla, P. pacifica, and P. hypochondriaca.     

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.     

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.     

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.     

Historical biogeography and post-glacial recolonization of South American temperate rain forest by the relictual marsupial Dromiciops gliroides

Aim Long‐term climatic variation has generated historical expansions and contractions of species ranges, with accompanying fragmentation and population bottlenecks, which are evidenced by spatial variation in genetic structure of populations. We examine here hypotheses concerning dispersal and vicariance in response to historical geoclimatic change and potential isolation produced by mountains and water barriers. Location The temperate rain forest of southern South America, which is distributed from coastal Chile, including the large continental island of Chiloé, across the Andes into Argentina. Methods We investigated our hypotheses in the phylogenetically and biogeographically relictual marsupial Dromiciops gliroides. We examined 56 specimens, which resulted from field samples and museum study skins from 21 localities. We evaluated the influence of two major barriers, the Andean cordillera and the waterway between the mainland and the large island of Chiloé, by performing Bayesian and maximum‐likelihood phylogenetic analyses on sequences of 877 base pairs of mitochondrial DNA. We further tested the contribution of the proposed geographical barriers using analysis of molecular variance (amova ). We also evaluated the responses of populations to historical north–south shifts of habitat associated with glacial history and sea‐level change. Results Our analyses revealed a phylogeny with three clades, two of which are widespread and contain nearly all the haplotypes: a northern clade (36–39° S) and a southern clade (40–43° S). These two clades contain forms from both sides of the Andes. Within the southern clade, island and mainland forms were not significantly differentiated. Tests of recent demographic change revealed that southern populations have experienced recent expansion, whereas northern populations exhibit long‐term stability. The direction of recent gene flow and range expansion is predominantly from Chile to Argentina, with a modest reciprocal exchange across the Andes. Recent gene flow from the island of Chiloé to the mainland is also supported. Main conclusions The genetic structure of contemporary D. gliroides populations suggests recent gene flow across the Andes and between the mainland and the island of Chiloé. Differences in demographic history that we detected between northern and southern populations have resulted from historical southward shifts of habitat associated with glacial recession in South America. Our results add to a growing literature that demonstrates the value of genetic data to illuminate how environmental history shapes species range and population structure.     

Phylogeographical analysis of an estuarine fish, Salanx ariakensis (Osmeridae: Salanginae) in the north-western Pacific

This study extended the geographic coverage of a previous study to explore population genetic structure and demographic history in the Ariake icefish Salanx ariakensis from three populations of continental coastlines and one island population in the north-western Pacific based on a partial sequence of the mitochondrial cytochrome b gene. The S. ariakensis showed high genetic diversity and strong genetic structure. Phylogenetic analysis showed a shallow gene tree with no clear phylogeographical structure. Contiguous range expansion and restricted gene flow were inferred to be main population events by nested-clade analysis. Significant genetic differentiations between populations could be attributable to negligible gene flow by coalescent analysis. High nucleotide diversity of each population was due to geographic mixing of heterogenous haplotypes during lowering sea levels of the Pleistocene. These findings indicate that cycles of geographic isolation and secondary contact happened in the Pleistocene glacial–interglacial cycles shaping genetic structure and population demography of S. ariakensis .     

Population genetic structure and phylogeography of the oak gall wasp Andricus chodjaii (Hymenoptera: Cynipidae) in Turkey as inferred from mitochondrial and nuclear DNA sequences

Sequence data of mitochondrial cytochrome b gene and nuclear ITS2 region were used to assess genetic diversity, intraspecific phylogeography and population genetic structure of the oak gall wasp Andricus chodjaii from Turkey. We examined 293 individuals from 21 localities which generated 57 cyt b haplotypes and 8 ITS2 alleles. The average genetic diversity was 0.575 for cyt b and 0.202 ITS2, and the average nucleotide diversity 0.015 for cyt b gene and 0.001 for the ITS2 region. Phylogenetic analyses of cyt b haplotypes produced mostly similar topologies with geographically significant groupings. The ITS2 data provided less resolution without robust and apparent geographic structure. Population demographic analysis indicated that some eastern populations expanded, however, some others underwent either expansion or decline resulting in genetically structured populations. Molecular clock applied to the mtDNA data indicated that ingroup haplotypes diversified from the outgroup haplotypes around Early Pliocene. Further diversification events throughout Pleistocene resulted in major clade formations. It appears that geographic formations and glacial and interglacial cycles of Pleistocene were crucial for shaping the phylogeographic structure of A. chodjaii in Turkey.     

Mitochondrial DNA reveals a mosaic pattern of phylogeographical structure in Atlantic and Mediterranean populations of anchovy (Engraulis encrasicolus)

This study extends the geographic coverage of a previous study of mitochondrial DNA restriction fragment length polymorphism in European anchovy. Both studies together include 24 samples representing 17 localities extending from the Black Sea, through the Mediterranean Sea to the eastern Atlantic as far south as Dakar, Senegal. Eighty-eight haplotypes define two clades (A and B) separated by 3.2% sequence divergence. Clade A has a star-like genealogy indicative of a recent population expansion. Clade B has a more complex genealogy, consisting of several haplotypes at intermediate frequencies. The distributions of these clades consist of a mosaic with abrupt changes between some areas and gradients between other areas. Clade A predominates the Black and Aegean seas, but is present throughout the Mediterranean. Unexpectedly, new data show that clade A is also at a high frequency in the Atlantic, from Portugal to at least Senegal. Overall, the level of genetic differentiation among populations is high (F(ST)=0.148, p<0.0001), with the greatest differences between basins. AMOVA reveals four main geographical groups: Atlantic, central Mediterranean, Aegean Sea, and Black Sea. Mismatch distribution clearly indicates historical bottleneck and population expansion for clade A, while for clade B such evidence is equivocal. This difference may reflect a range expansion for both clades, but with higher gene flow (Nm values) between demes for clade A. Both contemporary and historical processes are important in shaping the complex genetic population structure of European anchovy.     

The roles of rivers and Pleistocene refugia in shaping genetic diversity in Praomys misonnei in tropical Africa

Aim This study aims to elucidate the phylogeography of the murid rodent Praomys misonnei and to document whether or not rain forest refugia and rivers structure patterns of diversity within this species. Location Tropical Africa, from Ghana to Kenya. Methods Patterns of genetic structure and signatures of population history (cytochrome b gene) were assessed in a survey of 229 individuals from 54 localities. Using maximum likelihood, Bayesian, network and genetic structure analyses, we inferred intra‐specific relationships and tested hypotheses for historical patterns of gene flow within P. misonnei. Results Our phylogenetic analyses reveal a strong phylogeographical structure. We identified four major geographical clades within P. misonnei: one clade in Ghana and Benin, a Nigerian clade, a West Central African clade and a Central and East African clade. Several subclades were identified within these four major clades. A signal of population expansion was detected in most clades or subclades. Coalescence within all of the major clades of P. misonnei occurred during the Middle Pleistocene and/or the beginning of Late Pleistocene. Main conclusions Our results suggest a role for both Pleistocene refugia and rivers in structuring genetic diversity in P. misonnei. This forest‐dwelling rodent may have been isolated in a number of forest fragments during arid periods and expanded its range during wetter periods. Potential forest refugia may have been localized in Benin–Ghana, south‐western Cameroon, southern Gabon, northern Gabon and eastern Democratic Republic of Congo–western Uganda. The Niger and/or the Cross Rivers, the Oubangui‐Congo, Sanaga, Ogooue and/or Ivindo Rivers probably stopped the re‐expansion of the species from relict areas.