Phylogeography and postglacial dispersion of the chub (Leuciscus cephalus) in Europe

A phylogeographic analysis of mitochondrial DNA variation was performed in order to test the hypothesis of a postglacial recolonization of mid- and north-European rivers from a Danubian refuge. Over 345 chub specimens from European rivers covering most of the species' native range were investigated using 600 bp of the cytochrome b gene. Chub in European rivers belong to four highly divergent mitochondrial groups (lineages) differing by mean divergence estimates from 5.2% to 7.89%. These four lineages have a largely allopatric distribution, implying four geographical sets: two Mediterranean, and two north-European sets. This pattern provided strong evidence for: (i) the eradication of this species from most of Europe during maximum ice extent; (ii) its survival in four refugia (Adriatic side of the Balkans, eastern Greece (Aegean Rivers), southern tributaries of the Danube, and periphery of Black and Caspian Seas); (iii) a differential postglacial recolonization of mid- and northern Europe from the last two refugia only; (iv) the occurrence of this recolonization in two steps for the Danubian (western) lineage that entered western Europe (Rhine-Rhone-Loire drainages) during the Riss-Würm interglacial period and survived the last glaciation there before colonizing Garonne, UK and German drainages up to the Elbe during the Holocene; and (v) the occurrence of this recolonization in a single step for the Ponto-Caspian (eastern) lineage that entered the Baltic area as far as the Oder in the Holocene. Both lineages came into contact in the River Elbe without evident mixing.     

The Balkans and the colonization of Europe: the post‐glacial range expansion of the wild boar,Sus scrofa

Aim We focus on the biogeographical role of the Balkan Peninsula as a glacial refugium and source of northward post‐glacial dispersal for many European taxa. Specifically, we analysed the genetic structure and variation of wild boar (Sus scrofa) samples primarily from Greece, a region that has repeatedly served as a glacial refugium within the Balkan Peninsula. Location Continental Greece, the Aegean island of Samos and Bulgaria. Methods We analysed wild boar samples from 18 localities. Samples from common domestic breeds were also examined to take into account interactions between wild and domesticated animals. Phylogenetic analyses were carried out on a 637‐bp fragment of the mitochondrial DNA control region in 200 wild boar and 27 domestic pigs. The sequences were also compared with 791 Eurasian wild boar and domestic pig D‐loop sequences obtained from GenBank. Results Ninety‐four haplotypes were identified in the European wild boar data set, of which 68 were found in the Balkan samples and assigned to two previously described clades: the E1 European and Near Eastern clades. All of the continental samples clustered in the E1 clade and the samples from Samos fell into the Near Eastern clade, consistent with the island’s proximity to Asia Minor. Intriguingly, 62 novel haplotypes were identified and are found exclusively in the Balkans. Only six haplotypes were shared between wild boar and domestic pigs. Main conclusions Our data reveal numerous novel and geographically restricted haplotypes in wild boar populations, suggesting the presence of separate refugia in the Balkans. Our analyses support the hypothesis of a post‐glacial wild boar expansion consistent with the leading edge model, north and west from modern day Greece, and suggest little maternal introgression of Near Eastern and domestic haplotypes into wild Balkan populations.     

Phylogeography of Rhinichthys cataractae (Teleostei: Cyprinidae): pre‐glacial colonization across the Continental Divide and Pleistocene diversification within the Rio Grande drainage

The longnose dace, Rhinichthys cataractae, is a primary freshwater fish inhabiting riffle habitats in small headwater rivers and streams across the North American continent, including drainages east and west of the Continental Divide. The mitochondrially encoded cytochrome b gene (1140 bp) and 2298–2346 bp of the nuclear‐encoded genes S7 and RAG1 were obtained from 87 individuals of R. cataractae (collected from 17 sites throughout its range) and from several close relatives. Phylogenetic analyses recovered a monophyletic R. cataractae species‐group that contained Rhinichthys evermanni, Rhinichthys sp. ‘Millicoma dace’, and a non‐exclusive R. cataractae. Within the R. cataractae species‐group, two well‐supported lineages were identified, including a western lineage (containing R. evermanni, R. sp. ‘Millicoma dace’ and individuals of R. cataractae from Pacific slope drainages) and an eastern lineage (containing individuals of R. cataractae from Arctic, Atlantic, and Gulf slope drainages). Within the eastern lineage of R. cataractae, two well‐supported groups were recovered: a south‐eastern group, containing individuals from the Atlantic slope, southern tributaries to the Mississippi River, and the Rio Grande drainage; and a north‐eastern group, containing individuals from the Arctic slope and northern tributaries to the Mississippi River. Estimates of the timing of divergence within the R. cataractae species‐group, combined with ancestral area‐reconstruction methods, indicate a separation between the eastern and western lineages during the Pliocene to early‐Pleistocene, with a direction of colonization from the west of the Continental Divide eastward. Within the southern portion of its range, R. cataractae likely entered the Rio Grande drainage during the Pleistocene via stream capture events between the Arkansas River (Mississippi River drainage) and headwaters of the Rio Grande. A close relationship between populations of R. cataractae in the Rio Grande drainage and the adjacent Canadian River (Mississippi River drainage) is consistent with hypothesized stream capture events between the Pecos (Rio Grande drainage) and Canadian rivers during the late‐Pleistocene. The population of R. cataractae in the lower Rio Grande may have become separated from other populations in the Rio Grande drainage (upper Rio Grande and Pecos River) and Canadian River during the late‐Pleistocene, well before initiation of recent and significant anthropogenic disturbance within the Rio Grande drainage. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 317–333.     

Phylogeography of the olive sea snake, Aipysurus laevis (Hydrophiinae) indicates Pleistocene range expansion around northern Australia but low contemporary gene flow

Pleistocene sea-level fluctuations profoundly changed landmass configurations around northern Australia. The cyclic emergence of the Torres Strait land bridge and concomitant shifts in the distribution of shallow-water marine habitats repeatedly sundered east and west coast populations. These biogeographical perturbations invoke three possible scenarios regarding the directions of interglacial range expansion: west to east, east to west, or bidirectional. We evaluated these scenarios for the olive sea snake, Aipysurus laevis, by exploring its genetic structure around northern Australia based on 354 individuals from 14 locations in three regions (Western Australia, WA; Gulf of Carpentaria, GoC; Great Barrier Reef, GBR). A 726-bp fragment of the mitochondrial DNA ND4 region revealed 41 variable sites and 38 haplotypes, with no shared haplotypes among the three regions. Population genetic structure was strong overall, phiST=0.78, P<0.001, and coalescent analyses revealed no migration between regions. Genetic diversity was low in the GBR and GoC and the genetic signatures of these regions indicated range or population expansions consistent with their recent marine transgressions around 7000 years ago. By contrast, genetic diversity on most WA reefs was higher and there were no signals of recent expansion events on these reefs. Phylogenetic analyses indicated that GBR and GoC haplotypes were derived from WA haplotypes; however, statistical parsimony suggested that recent range expansion in the GBR-GoC probably occurred from east coast populations, possibly in the Coral Sea. Levels of contemporary female-mediated gene flow varied within regions and reflected potential connectivity among populations afforded by the different regional habitat types.     

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.     

Subantarctic flowering plants: pre‐glacial survivors or post‐glacial immigrants?

Aim The aim here was to assess whether the present‐day assemblage of subantarctic flowering plants is the result of a rapid post‐Last Glacial Maximum (LGM) colonization or whether subantarctic flowering plants survived on the islands in glacial refugia throughout the LGM. Location The circumpolar subantarctic region, comprising six remote islands and island groups between latitudes 46° and 55° S, including South Georgia in the South Atlantic Ocean, the Prince Edward Islands, Îles Crozet, Îles Kerguelen, the Heard Island group in the South Indian Ocean and Macquarie Island in the South Pacific Ocean. Methods Floristic affinities between the subantarctic islands were assessed by cluster analysis applied to an up‐to‐date dataset of the phanerogamic flora in order to test for the existence of provincialism within the subantarctic. A review of the primary literature on the palaeobotany, geology and glacial history of the subantarctic islands was carried out and supplemented with additional palaeobotanical data and new field observations from South Georgia, Île de la Possession (Îles Crozet) and Îles Kerguelen. Results First, a strong regionalism was observed, with different floras characterizing the islands in each of the ocean basins, and endemic species being present in the South Indian Ocean and South Pacific Ocean provinces. Second, the majority of the plant species were present at the onset of accumulation of post‐glacial organic sediment and there is no evidence for the natural arrival of new immigrants during the subsequent period. Third, a review of geomorphological data suggested that the ice cover was incomplete during the LGM on the majority of the islands, and ice‐free biological refugia were probably present even on the most glaciated islands. Main conclusions Several independent lines of evidence favour the survival of a native subantarctic phanerogamic flora in local refugia during the LGM rather than a post‐LGM colonization from more distant temperate landmasses in the Southern Hemisphere.     

Multiple Pleistocene refugia and post‐glacial colonization in the European chub (Squalius cephalus) revealed by combined use of nuclear and mitochondrial markers

Aim To analyse patterns of nuclear and mitochondrial genetic variation in the European chub, Squalius cephalus (Linnaeus, 1758), in order to understand the evolutionary history of this species and to test biogeographical hypotheses for the existence of co‐distributed European freshwater fish species. Location Rivers in Europe (Finland, Poland, Czech Republic, France, Bulgaria, Spain, Italy). Methods We genotyped 12 polymorphic microsatellite markers derived from 310 individuals collected from across the distribution of S. cephalus in Europe (including a total of 15 populations) and sequenced mitochondrial DNA (mtDNA) from a subset of 75 individuals. Sequences of mtDNA cytochrome b were analysed using both phylogenetic (median‐joining networks) and population genetic methods (tests for demographic history, mismatch distributions, Bayesian coalescent analysis). Geographical structure in microsatellite loci was examined using a distance method (F_ST), factorial correspondence analysis (FCA) and a Bayesian clustering method (structure ). Results The mtDNA network showed a clear split into four different haplogroup lineages: Western (separated into Atlantic and Danubian sublineages), Eastern, Aegean (occurring in two distinct sublineages in the Balkans and in Spain) and Adriatic. Our results indicate recent population expansion in the Eastern and Western Atlantic lineages and the admixture of two previously separate sublineages (Atlantic and Danubian) in the Western lineage. Bayesian structure analysis as well as FCA results roughly corresponded to the mtDNA‐based structure, separating the sampled individuals into almost non‐overlapping groups. Main conclusions Our results support hypotheses suggesting origins of extant lineages of freshwater fishes in multiple refugia and the subsequent post‐glacial colonization of Europe via different routes. We confirmed the previously proposed two‐step expansion scenario from the Danube refuge, the existence of a secondary (Atlantic) refuge during the last glaciation (probably in the Rhone River) and population expansion of this lineage. Conspicuous divergences among Mediterranean populations reflect their different origin, as well as their low contribution to the recent genetic pool of chub in central Europe.     

Pleistocene glacial cycles drove lineage diversification and fusion in the Yosemite toad (Anaxyrus canorus)

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this “glacial pulse” model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed (“fused”) lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double‐digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low‐ and high‐elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.     

Phylogeography of the tidewater goby, Eucyclogobius newberryi (Teleostei, Gobiidae), in coastal California

The tidewater goby, Eucyclogobius newberryi, inhabits discrete, seasonally closed estuaries and lagoons along approximately 1500 km of California coastline. This species is euryhaline but has no explicit marine stage, yet population extirpation and recolonization data suggest tidewater gobies disperse intermittently via the sea. Analyses of mitochondrial control region and cytochrome b sequences demonstrate a deep evolutionary bifurcation in the vicinity of Los Angeles that separates southern California populations from all more northerly populations. Shallower phylogeographic breaks, in the vicinities of Seacliff, Point Buchon, Big Sur, and Point Arena segregate the northerly populations into five groups in three geographic clusters: the Point Conception and Ventura groups between Los Angeles and Point Buchon, a lone Estero Bay group from central California, and San Francisco and Cape Mendocino groups from northern California. The phylogenetic relationships between and patterns of molecular diversity within the six groups are consistent with repeated, and sometimes rapid, northward and southward range expansions out of central California caused by Quaternary climate change. Plio-Pleistocene tectonism, Quaternary coastal geography and hydrography, and historical human activities probably also influenced the modern geographic and genetic structure of E. newberryi. The phylogeography of E. newberryi is concordant with phylogeographic patterns in several other coastal California taxa, suggesting common extrinsic factors have had similar effects on different species. However, there is no evidence of a phylogeographic break coincident with a biogeographic boundary at Point Conception.     

Phylogeography of the northern hogsucker, Hypentelium nigricans (Teleostei: Cypriniformes): genetic evidence for the existence of the ancient Teays River

Aim To assess the roles of dispersal and vicariance in shaping the present distribution and diversity within Hypentelium nigricans, the northern hogsucker (Teleostei: Cypriniformes). Location Eastern United States. Methods Parsimony analyses, Bayesian analyses, pairwise genetic divergence and mismatch plots are used to examine patterns of genetic variation across H. nigricans. Results Species relationships within the genus Hypentelium were consistent with previous hypotheses. However, relationships between haplotypes within H. nigricans revealed two deeply divergent groups, a clade containing haplotypes from the New and Roanoke rivers (Atlantic Slope) plus Interior Highlands and upper Mississippi River and a clade containing haplotypes from the Eastern Highlands, previously glaciated regions of the Ohio and Wabash rivers, and the Amite and Homochitto rivers of south‐western Mississippi. Main conclusions The phylogenetic history of Hypentelium was shaped by old vicariant events associated with erosion of the Blue Ridge and separation of the Mobile and Mississippi river basins. Within H. nigricans two clades existed prior to the Pleistocene; a widespread clade in the pre‐glacial Teays‐Mississippi River system and a clade in Cumberland and Tennessee rivers. Pleistocene events fragmented the Teays‐Mississippi fauna. Following the retreat of the glaciers H. nigricans dispersed northward into previously glaciated regions. These patterns are replicated in other clades of fishes and are consistent with some of the predictions of Mayden's (Systematic Zoology, 37, 329, 1988) pre‐Pleistocene vicariance hypothesis.     

Post-pleistocene demographic history of the North Atlantic endemic Irish moss Chondrus crispus: glacial survival, spatial expansion and gene flow

Range expansions and gene flow as micro-evolutionary processes played a leading role in the population demographic history of marine organisms. Herein, we sequenced partial mtDNA Cox1 gene from 26 assigned geographical populations to understand how Irish moss (Chondrus crispus) responded to severe climatic oscillations during the Pleistocene glaciations and contemporary forces such as gene flow. Phylogeographic patterns indicated that haplotype frequency distributions were strongly skewed, with nearly half found only in single samples and thus restricted to a single population. Analysis of molecular variance revealed that most of the variation was within populations with no significant genetic structuring on either side of the Atlantic. Demographic analyses indicated that ISI (Irish Sea and Ireland) and NS (the North Sea) areas experienced a slight trend of increase in population size over time, whereas EC (the English Channel) area experienced expansion beginning approximately 170,000-360,000 BP. The observed complex genetic pattern of C. crispus is consistent with a scenario of multiple unrelated founding events by survival of this species in at least three putative Pleistocene refugia along the European coastline, and subsequent trans-Atlantic dispersal combined with contiguous northward population expansion predating the LGM and geographically gene flow.     

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.