Phylogeography and ecological niche modelling implicate coastal refugia and trans‐alpine dispersal of a New Zealand fungus beetle

The Last Glacial Maximum (LGM) severely restricted forest ecosystems on New Zealand’s South Island, but the extent of LGM distribution for forest species is still poorly understood. We used mitochondrial DNA phylogeography (COI) and ecological niche modelling (ENM) to identify LGM refugia for the mycophagous beetle Agyrtodes labralis (Leiodidae), a forest edge species widely distributed in the South Island. Both the phylogenetic analyses and the ENM indicate that A. labralis refuged in Kaikoura, Nelson, and along much of the South Island’s west coast. Phylogeography of this species indicates that recolonization of the largely deforested east and southeast South Island occurred in a west–east direction, with populations moving through the Southern Alps, and that the northern refugia participated little in interglacial population expansion. This contradicts published studies of other New Zealand species, in which recolonization occurs in a north–south fashion from many of the same refugia.     

Phylogeographic history of the New Zealand stick insect Niveaphasma annulata (Phasmatodea) estimated from mitochondrial and nuclear loci

We have assessed the utility of a single-copy nuclear locus and mitochondrial DNA (mtDNA) in a phylogeographic study of the New Zealand stick insect Niveaphasma annulata (Hutton). We amplified sequences from the mitochondrial cytochrome oxidase subunit I (COI) gene and the single-copy nuclear gene elongation factor-1α (EF1α) from 97 individuals. Allelic phase at the EF1α locus was determined using Denaturing Gradient Gel Electrophoresis. Phylogenetic analyses showed broad congruence between the geographic distribution of three major COI clades and EF1α alleles, which suggested that the phylogenetic patterns reflect population history rather than lineage sorting. However, the geographic boundaries of these clades were not always in exact agreement between the two loci. Our data indicate that Niveaphasma annulata was most likely separated into a number of refugia during Pleistocene glacial advances. Subsequent to glacial retreat these refugial populations have expanded and now form a number of zones of secondary contact. We contrast these patterns with those observed from other New Zealand taxa. Our study offers compelling evidence for the use of nuclear genes alongside mtDNA for future phylogeographic studies.     

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.     

Determining the origin and age of the Westland beech (Nothofagus) gap, New Zealand, using fungus beetle genetics

The formation and maintenance of the Nothofagus beech gap in the South Island, New Zealand, has been the focus of biogeographical debate since the 1920s. We examine the historical process of gap formation by investigating the population genetics of fungus beetles: Brachynopus scutellaris (Staphylinidae) inhabits logs and is absent from the beech gap, and Hisparonia hystrix (Nitidulidae) is contiguous through the gap and is found commonly on sooty mould growing on several plant species. Both species show distinctive northern and southern haplotype distributions while H. hystrix recolonized the gap as shown by definitive mixing. B. scutellaris shows two major haplotype clades with strong geographical concordance, and unlike H. hystrix, has clearly defined lineages that can be partitioned for molecular dating. Based on coalescence dating methods, disjunct lineages of B. scutellaris indicate that the gap was formed less than 200 000 years ago. Phylogenetic imprints from both species reveal similar patterns of population divergence corresponding to recent glacial cycles, favouring a glacial explanation for the origin of the gap. Post-gap colonization by H. hystrix may have been facilitated by the spread of Leptospermum scoparium host trees to the area, and they may be better at dispersing than B. scutellaris which may be constrained by fungal host and/or microhabitat. The gap-excluded species B. scutellaris is found in both beech and podocarp-broadleaf forests flanking the Westland gap and its absence in the gap may be related to incomplete recolonization following glacial retreat. We also discuss species status and an ancient polymorphism within B. scutellaris .     

Investigating hybridization in the parthenogenetic New Zealand stick insect Acanthoxyla (Phasmatodea) using single-copy nuclear loci

The New Zealand stick insect genus Acanthoxyla Uvarov is extremely unusual among higher taxa of animals in that all known species are obligate parthenogens. We have used a combination of the mitochondrial DNA genes cytochrome oxidase subunits I and II, 28S nuclear ribosomal RNA, and the two single-copy nuclear genes elongation factor 1alpha and phosphoglucose isomerase to test hypotheses on the role of hybridization in the evolution of this genus. Alleles at the single-copy nuclear loci in three sampled species of Acanthoxyla were resolved by cloning the PCR products. Analysis of multilocus genotypes shows that most sampled individuals of Acanthoxyla possess three alleles at the single-copy nuclear loci, which we have interpreted to indicate triploidy. Because most of the alleles from Acanthoxyla form a monophyletic group, including sets of alleles possessed by the putative triploids, we have inferred that the extant parthenogenetic lineages formed via hybridization between species of Acanthoxyla, at least one of which must have been sexual. More recently, there have been multiple introgression events from the related species Clitarchus hookeri White, although C. hookeri does not appear to be involved with the origin of parthenogenesis in Acanthoxyla. Our study demonstrates the utility of cloning alleles from multiple single-copy nuclear genes for resolving the origins of parthenogenetic lineages.     

Inferring the past to predict the future: climate modelling predictions and phylogeography for the freshwater gastropod Radix balthica (Pulmonata, Basommatophora)

Understanding the impact of past climatic events on species may facilitate predictions of how species will respond to future climate change. To this end, we sampled populations of the common pond snail Radix balthica over the entire species range (northwestern Europe). Using a recently developed analytical framework that employs ecological niche modelling to obtain hypotheses that are subsequently tested with statistical phylogeography, we inferred the range dynamics of R. balthica over time. A Maxent modelling for present-day conditions was performed to infer the climate envelope for the species, and the modelled niche was used to hindcast climatically suitable range at the last glacial maximum (LGM) c . 21 000 years ago. Ecological niche modelling predicted two suitable areas at the LGM within the present species range. Phylogeographic model selection on a COI mitochondrial DNA data set confirmed that R. balthica most likely spread from these two disjunct refuges after the LGM. The match observed between the potential range of the species at the LGM given its present climatic requirements and the phylogeographically inferred refugial areas was a clear argument in favour of niche conservatism in R. balthica , thus allowing to predict the future range. The subsequent projection of the potential range under a global change scenario predicts a moderate pole-ward shift of the northern range limits, but a dramatic loss of areas currently occupied in France, western Great Britain and southern Germany.     

Postglacial recolonization at a snail's pace (Trochulus villosus): confronting competing refugia hypotheses using model selection

The localization of Last Glacial Maximum (LGM) refugia is crucial information to understand a species' history and predict its reaction to future climate changes. However, many phylogeographical studies often lack sampling designs intensive enough to precisely localize these refugia. The hairy land snail Trochulus villosus has a small range centred on Switzerland, which could be intensively covered by sampling 455 individuals from 52 populations. Based on mitochondrial DNA sequences (COI and 16S), we identified two divergent lineages with distinct geographical distributions. Bayesian skyline plots suggested that both lineages expanded at the end of the LGM. To find where the origin populations were located, we applied the principles of ancestral character reconstruction and identified a candidate refugium for each mtDNA lineage: the French Jura and Central Switzerland, both ice-free during the LGM. Additional refugia, however, could not be excluded, as suggested by the microsatellite analysis of a population subset. Modelling the LGM niche of T. villosus , we showed that suitable climatic conditions were expected in the inferred refugia, but potentially also in the nunataks of the alpine ice shield. In a model selection approach, we compared several alternative recolonization scenarios by estimating the Akaike information criterion for their respective maximum-likelihood migration rates. The 'two refugia' scenario received by far the best support given the distribution of genetic diversity in T. villosus populations. Provided that fine-scale sampling designs and various analytical approaches are combined, it is possible to refine our necessary understanding of species responses to environmental changes.     

Intraspecific allozymatic differentiation reveals the glacial refugia and the postglacial expansions of European Erebia medusa (Lepidoptera: Nymphalidae)

Allozyme analysis of Erebia medusa over large regions of Europe revealed a significant population differentiation (F_ST: 0.149 ± 0.016). A UPGMA-analysis showed a division into four major lineages with mean inter-group genetic distances ranging from 0.051 (±0.010) to 0.117 (±0.024). An AMOVA revealed that rather more than two-thirds of the variance between samples was being between these lineages and less than one-third within lineages. An eastern group included the samples from the Czech Republic, Slovakia and north-eastern Hungary. This genetic lineage expressed significantly higher genetic diversity than the other three. A second lineage was formed by the samples from France and Germany. The two samples from western Hungary represent a third delimited lineage and the sample from northern Italy a fourth. We suppose that this genetic differentiation took place during the last ice-age in four disjunct refugia. The genetically more diverse eastern genetic lineage might have evolved in a relatively large refugium in south-eastern Europe. We assume that the other three lineages developed in relatively small relict areas around the Alps. It is likely for the western lineage that its ice-age distribution showed at least one disjunction in late Würm with the consequence of further genetic differentiation. Most probably, the eastern lineage colonized postglacial Central Europe using two alternative routes: one north and one south of the Carpathians. Up to now, neither similar glacial refugia, nor comparable secondary disjunctions in late Würm, are reported for any other animal or plant species.     

Early stages of divergence: phylogeography, climate modeling, and morphological differentiation in the South American lizard Liolaemus petrophilus (Squamata: Liolaemidae)

This study examines the phylogeographic structure within the Patagonian lizard Liolaemus petrophilus and tests for patterns of between-clade morphological divergence and sexual dimorphism, as well as demographic and niche changes associated with Pleistocene climate changes. We inferred intraspecific relationships, tested hypotheses for historical patterns of population expansion, and incorporated ecological niche modeling (ENM) with standard morphological and geometric morphometric analyses to examine between-clade divergence as indirect evidence for adaptation to different niches. The two inferred haploclades diverged during the early Pleistocene with the Southern clade depicting the genetic signature of a recent population increase associated with expanding niche envelope, whereas the Northern clade shows stable populations in a shrinking niche envelope. The combination of molecular evidence for postisolation demographic change and ENM, suggest that the two haploclades have responded differently to Pleistocene climatic events.     

External macro- and micromorphology of the male of the stick insect Hermarchus leytensis (Insecta: Phasmatodea) with phylogenetic considerations

External morphological features of adult males are described in the stick insect Hermarchus leytensis from the Philippines, a species belonging to the little-known euphasmatodean lineage Stephanacridini. Mouthparts are characterized by some likely specialized features: (1) a dentate dorsal cutting edge on the mandibles; (2) distinct differences in size and shape between the galeae; (3) absence of an apical field of trichomes on the galeae; and (4) lacinial setae not protruding from the mesal margin, which features three bearing-like protuberances. The latter character state represents a very unusual condition, not known in other phasmatodeans or even in polyneopteran insects. A distinctive characteristic of attachment devices is that each euplantula is divided into two separated pads with a smooth microstructure, as it also occurs in some members of the clade Schizodecema. Male terminalia exhibit character states previously unknown in Stephanacridini, including (1) a symmetrical type of vomer and (2) claspers equipped with ∼70 very minute ventral teeth on each thorn pad. Potential implications for the systematic placement of H. leytensis are discussed. The results also underline the importance of microanatomical investigations as source of substantial characters for future analyses on phasmatodean systematics.     

A SINE of restricted gene flow across the Alpine Fault: phylogeography of the New Zealand common skink (Oligosoma nigriplantare polychroma)

New Zealand has experienced a complex climatic and geological history since the Pliocene. Thus, identifying the processes most important in having driven the evolution of New Zealand's biota has proven difficult. Here we examine the phylogeography of the New Zealand common skink ( Oligosoma nigriplantare polychroma ) which is distributed throughout much of New Zealand and crosses many putative biogeographical boundaries. Using mitochondrial DNA sequence data, we revealed five geographically distinct lineages that are highly differentiated (pairwise Φ_ST 0.54–0.80). The phylogeographical pattern and inferred age of the lineages suggests Pliocene mountain building along active fault lines promoted their divergence 3.98–5.45 million years ago. A short interspersed nuclear element (SINE) polymorphism in the myosin gene intron ( MYH-2 ) confirmed a pattern of restricted gene flow between lineages on either side of the mountain ranges associated with the Alpine Fault that runs southwest to northeast across the South Island of New Zealand. An analysis of molecular variance confirmed that ~40% of the genetic differentiation in O. n. polychroma is distributed across this major fault line. The straits between the main islands of New Zealand accounted for much less of the variation found within O. n. polychroma , most likely due to the repeated existence of landbridges between islands during periods of the Pleistocene that allowed migration. Overall, our findings reveal the relative roles of different climatic and geological processes, and in particular, demonstrate the importance of the Alpine Fault in the evolution of New Zealand's biota.     

Phylogeography and genetic structure of the red-legged partridge (Alectoris rufa): more evidence for refugia within the Iberian glacial refugium

The Pleistocene climatic oscillations promoted the diversification in avian species during the last glacial period. The red‐legged partridge (Alectoris rufa, Family Phasianidae) has a large natural distribution extending from the Mediterranean to humid temperate zones. However, the genetic structure for this species is unknown. The present study investigates the phylogeography, genetic structure and demographic history of A. rufa across its distribution, employing both mitochondrial DNA control region sequences and nuclear microsatellite loci. Our results propose that this species was greatly affected by Pleistocene glaciations. The mismatch analyses suggest that the current populations resulted from post‐glacial expansion and subsequent differentiation resulting in five diagnosable genetic clusters: Southwestern, Central‐eastern, Northwestern, Balearic and French and Italian. Further, we found evidence of three glacial refugia within the currently recognized Iberian glacial refugium. The intraspecific structure revealed by both maternal and biparental phylogeographic analyses was not resolved in the phylogenetic analyses. Based on all considerations, we recommended that five management units be recognized.     

Combined ecological niche modelling and molecular phylogeography revealed the evolutionary history of Hordeum marinum (Poaceae)--niche differentiation, loss of genetic diversity, and speciation in Mediterranean Quaternary refugia

The Hordeum marinum species group consists of two annual grasses of western Eurasian saline meadows or marshes. The two grasses split in the Quaternary about two million years ago. Hordeum marinum and the diploid of Hordeum gussoneanum (2x) co-occur throughout the Mediterranean basin, while the autotetraploid cytotype of H. gussoneanum (4x) overlaps with its diploid progenitor geographically only in the utmost Eastern Mediterranean, extending from there eastwards into Asia. Using chloroplast sequences of the trnL-F region, six newly developed chloroplast microsatellite loci, ecological predictive models based on climate data, and the present geographical distribution of the two species we analysed differentiation processes in the H. marinum group. The chloroplast data indicated clear differences in the history of both species. For H. marinum we found a subdivision between genetically variable populations from the Iberian Peninsula and the more uniform populations from the remaining Mediterranean. As an explanation, we assume Pleistocene fragmentation of an earlier widespread population and survival in an Iberian and a Central Mediterranean glacial refuge. Chloroplast variation was completely absent within the cytotypes of H. gussoneanum, indicating a severe and recent genetic bottleneck. Due to this lack of chloroplast variation only the combination of ecological habitat modelling with molecular data analyses allowed conclusions about the history of this taxon. The distribution areas of the two cytotypes of H. gussoneanum overlap today in parts of Turkey, indicating an area with similar climate conditions during polyploid formation. However, after its origin the polyploid cytotype underwent a pronounced ecological shift, compared to its diploid progenitor, allowing it to colonize mountainous inland habitats between the Mediterranean basin and Afghanistan. The extant sympatric occurrence of H. marinum and H. gussoneanum 2x in the Mediterranean region is interpreted as a result of secondary contact after fast Holocene range expansion out of different ice age refugia.     

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae)

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.     

Persistence and diversification of the Holarctic shrew,Sorex tundrensis (Family Soricidae), in response to climate change

Environmental processes govern demography, species movements, community turnover and diversification and yet in many respects these dynamics are still poorly understood at high latitudes. We investigate the combined effects of climate change and geography through time for a widespread Holarctic shrew, Sorex tundrensis. We include a comprehensive suite of closely related outgroup taxa and three independent loci to explore phylogeographic structure and historical demography. We then explore the implications of these findings for other members of boreal communities. The tundra shrew and its sister species, the Tien Shan shrew (Sorex asper), exhibit strong geographic population structure across Siberia and into Beringia illustrating local centres of endemism that correspond to Late Pleistocene refugia. Ecological niche predictions for both current and historical distributions indicate a model of persistence through time despite dramatic climate change. Species tree estimation under a coalescent process suggests that isolation between populations has been maintained across timeframes deeper than the periodicity of Pleistocene glacial cycling. That some species such as the tundra shrew have a history of persistence largely independent of changing climate, whereas other boreal species shifted their ranges in response to climate change, highlights the dynamic processes of community assembly at high latitudes.     

Diversification of New Zealand weta (Orthoptera: Ensifera: Anostostomatidae) and their relationships in Australasia

New Zealand taxa from the Orthopteran family Anostostomatidae have been shown to consist of three broad groups, Hemiandrus (ground weta), Anisoura/Motuweta (tusked weta) and Hemideina-Deinacrida (tree-giant weta). The family is also present in Australia and New Caledonia, the nearest large land masses to New Zealand. All genera are endemic to their respective countries except Hemiandrus that occurs in New Zealand and Australia. We used nuclear and mitochondrial DNA sequence data to study within genera and among species-level genetic diversity within New Zealand and to examine phylogenetic relationships of taxa in Australasia. We found the Anostostomatidae to be monophyletic within Ensifera, and justifiably distinguished from the Stenopelmatidae among which they were formerly placed. However, the New Zealand Anostostomatidae are not monophyletic with respect to Australian and New Caledonian species in our analyses. Two of the New Zealand groups have closer allies in Australia and one in New Caledonia. We carried out maximum-likelihood and Bayesian analyses to reveal several well supported subgroupings. Our analysis included the most extensive sampling to date of Hemiandrus species and indicate that Australian and New Zealand Hemiandrus are not monophyletic. We used molecular dating approaches to test the plausibility of alternative biogeographic hypotheses for the origin of the New Zealand anostostomatid fauna and found support for divergence of the main clades at, or shortly after, Gondwanan break-up, and dispersal across the Tasman much more recently.     

Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate

Genetic data are currently providing a large amount of new information on past distribution of species and are contributing to a new vision of Pleistocene ice ages. Nonetheless, an increasing number of studies on the 'time dependency' of mutation rates suggest that date assessments for evolutionary events of the Pleistocene might be overestimated. We analysed mitochondrial (mt) DNA (COI) sequence variation in 225 Parnassius mnemosyne individuals sampled across central and eastern Europe in order to assess (i) the existence of genetic signatures of Pleistocene climate shifts; and (ii) the timescale of demographic and evolutionary events. Our analyses reveal a phylogeographical pattern markedly influenced by the Pleistocene/Holocene climate shifts. Eastern Alpine and Balkan populations display comparatively high mtDNA diversity, suggesting multiple glacial refugia. On the other hand, three widely distributed and spatially segregated lineages occupy most of northern and eastern Europe, indicating postglacial recolonization from different refugial areas. We show that a conventional 'phylogenetic' substitution rate cannot account for the present distribution of genetic variation in this species, and we combine phylogeographical pattern and palaeoecological information in order to determine a suitable intraspecific rate through a Bayesian coalescent approach. We argue that our calibrated 'time-dependent' rate (0.096 substitutions/ million years), offers the most convincing time frame for the evolutionary events inferred from sequence data. When scaled by the new rate, estimates of divergence between Balkan and Alpine lineages point to c. 19 000 years before present (last glacial maximum), and parameters of demographic expansion for northern lineages are consistent with postglacial warming (5-11 000 years before present).     

Biodiversity in the New Zealand divaricating tree daisies: Olearia sect. nov. (Compositae)

The trees and shrubs in Olearia sect. Divaricaster, sect. nov., are found in North, South, and Stewart Islands of New Zealand. They are distinguished by their small, opposite leaves borne on brachyblasts (short shoots), long shoots which abort apically, solitary or fascicled capitula, flowers with purple style arms, and a very distinctive insect fauna. The ten species and three subspecies, including the new O. fimbriata, O. gardneri and O. virgata subsp. centralis are revised and a key and distribution maps provided. The section includes four threatened species: O. hectorii, O. polita and the two new species. Species density is mapped in grid cells 1° latitude by 1° longitude. Areas west of the geological Median Tectonic Zone (Nelson, Westland and Fiordland) show relatively low diversity. The highest diversity, six species in one grid cell, is recorded in north-west Otago, around Mt Aspiring and the Humboldt Mountains. This region is also geologically significant as several allochthonous terranes are juxtaposed there. In the central South Island O. laxiflora, O. fimbriata and O. odorata have concentric distributions around what were Tertiary basins and inland seas. Members of the section favour fertile, lowland sites on alluvium and have suffered reduction in population size over the last century due to agricultural development. Several species, such as O. hectorii, the 'swamp gum', occur in sites which are periodically flooded. Olearia solandri is a coastal species, often growing as a mangrove associate, but occasionally found inland. The anomalous inland populations are discussed and their location attributed to geological uplift.     

Phylogeography of a widespread terrestrial vertebrate in a barely-studied Palearctic region: green toads (<Emphasis Type=Italic>Bufo viridis</Emphasis> subgroup) indicate glacial refugia in Eastern Central Asia

The phylogeography of western Palearctic species is relatively well studied, but data on Eastern Central Asia are scarce. We present one of the first data sets from a widespread terrestrial vertebrate (Bufo pewzowi) inhabiting Eastern Central Asian mountains and deserts to gain knowledge on its phylogeography in this region. We applied combined phylogenetic and demographic analyses to understand the evolutionary history using mitochondrial DNA D-loop variation of toads from 37 locations. Genetic structure of Bufo pewzowi is strongly affected by landscape: we found three haplotype groups in eastern Kazakhstan, Dzungaria and Tarim Basin, divided by the Tian Shan and Dzungarian Alatau ranges. A vicariant hypothesis may explain divergence among groups. The divergence time of the three major clades was estimated about 0.9 million years ago (confidence interval 0.5–1.4), and is discussed with respect to Quaternary uplifting and glaciation in the Tian Shan. Demographic analyses provided evidence for both historical bottlenecks and population expansions and suggested Pleistocene signatures. Glacial refugia were inferred in the Tarim Basin (around the Turpan depression), in southern Dzungaria (Urumqui region), at the northern foot of the Tian Shan (Gongnaisi) and perhaps at the Altai range (Terekti). Regional Post-Last Glacial Maximum dispersal patterns are proposed. A taxonomic hypothesis is presented. This study provides a detailed history of how a widespread terrestrial vertebrate responded to geological change and Quaternary glacial events in Eastern Central Asia and may have significance for future phylogeographic research in this understudied region.     

Late Pleistocene freshwater fish (Cottidae) trackways from New England (USA) glacial lakes and a reinterpretation of the ichnogenus Broomichnium Kuhn

New trace fossil material from Late Pleistocene glaciolacustrine varves of the Connecticut River Valley, New England, USA represent the first evidence of freshwater sculpin in glacial Lake Hitchcock. Paleobiogeographic data constrain the timing of the last reinhabitation of freshwater cottids from a Wisconsinan glacial refugium. Freshwater sculpin were present in the area of study by ~ 13.7 ky BP, moving approximately 400 km in 5000 yr, and following the ice margin at distances as close as 35 km. The trace fossils warrant erection of a new ichnospecies, Broomichnium flirii isp. nov. Comparison of this new ichnospecies to Broomichnium permianum reveals distinct similarities, and it is possible that the Permian examples of B. permianum also were made by fish, which would reconcile a long-running controversy. Many groups of fish are demersal and make ventral body contact with the substrate and could potentially leave similar traces. Identifying new forms of trace fossils made by fish that use alternative modes of locomotion will prove useful in paleoenvironmental interpretations.