Differential cycles of range contraction and expansion in European high mountain plants during the Late Quaternary: insights from Pritzelago alpina (L.) O. Kuntze (Brassicaceae)

Nuclear DNA sequence variation of the internal transcribed spacer (ITS) and amplified fragment length polymorphisms (AFLPs) were used to illuminate the evolutionary history of Pritzelago alpina, a herbaceous perennial of (sub)alpine to nival habitats of the European high mountains. Maximum likelihood analysis of ITS sequences of P. alpina, Hornungia petraea and Hymenolobus procumbens (the 'Pritzelago alliance') resolved P. alpina and H. petraea as sister taxa. ITS divergence estimates support an origin for P. alpina in the Late Tertiary, while intraspecific diversification started in the Late Quaternary (0.4-0.9 million years ago). AFLP analysis of 76 individuals of P. alpina, representing 24 localities across its entire west-east distribution, identified four mountain lineages in Cantabria, the Pyrenees, (south-) western Alps, and northeastern Alps/Tatras/Carpathians. In an analysis of molecular variance (amova), 14.3% of the total variation derived from this separation. However, relationships among these lineages remained unresolved in neighbour-joining and principal co-ordinates analyses, suggesting a population history of near simultaneous vicariance. Comparison with our previous ITS/AFLP study of Anthyllis montana (Fabaceae) indicates that the two co-distributed but altitudinally differentiated plant species exhibit temporally concordant but spatially discordant patterns of genetic variation. Moreover, levels of AFLP divergence were significantly lower in P. alpina than in the submediterranean, lower-elevation A. montana. Together, these data are consistent with a 'displacement refugia model', which predicts that European mountain plant species associated with lower- and upper-elevation habitats had a different cycle of range contraction into (long-term) glacial and (short-term) interglacial refugia, respectively.     

Chloroplast DNA diversity of the dioecious European tree Ilex aquifolium L. (English holly)

Variation in the chloroplast genome of Ilex aquifolium (English holly), a dioecious evergreen tree native to south, west and central Europe, was analysed using polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLPs) and microsatellites. Differentiation between populations was high (GST = 0.595) and evidence for phylogeographical structure was detected (NST = 0.697, significantly higher than GST). Two chloroplast lineages were inferred originating from putative glacial refugia in southern Europe (Iberia, Italy and possibly the Balkans). The GST value was higher than reported for endozoochorous hermaphrodite species and for anemochorous dioecious species investigated over a similar geographical scale. It appears that dioecy has contributed to strong differentiation between refugia and that this has been maintained following postglacial recolonization as a result of limited seed flow. Palynological records for I. aquifolium are poor, thus these results give an important insight into patterns of glacial isolation and postglacial recolonization of this species.     

Phylogeography of the intertidal copepod Tigriopus californicus reveals substantially reduced population differentiation at northern latitudes

Previous studies of the intertidal copepod Tigriopus californicus revealed one of the highest levels of mitochondrial DNA differentiation ever reported among conspecific populations. The present study extends the geographical sampling northward, adding populations from northern California to south-east Alaska. The mitochondrial phylogeny for the entire species range, based on cytochrome oxidase I sequences for a total of 49 individuals from 27 populations, again shows extreme differentiation among populations (up to 23%). However, populations from Oregon northwards appear to be derived and have interpopulation divergences five times lower than those between southern populations. Furthermore, although few individuals were sequenced from each locality, populations from Puget Sound northward had significantly reduced levels of within-population variation. These patterns are hypothesized to result from the contraction and expansion of populations driven by recent ice ages.     

Geographic analysis of nucleotide diversity and song sparrow (Aves: Emberizidae) population history

Mitochondrial DNA (mtDNA) control-region (CR) sequences were analysed to address three questions regarding the evolution of geographical variation in song sparrows. (i) Are mtDNA sequences more informative about phylogenetic relationships and population history than previously published restriction fragment (RFLP) data? (ii) Are song sparrow CR sequences evolving in a selectively neutral manner? (iii) What do the haplotype cladogram and geographical pattern of nucleotide diversity (π) suggest about the recent evolutionary history of song sparrow populations? Results from phylogenetic analyses of CR sequences corroborate RFLP results and reveal instances in which haplotypes do not group by locality. Neutrality tests ( 51 ) suggest that song sparrow mtDNA is evolving in a selectively neutral manner, although exceptions are noted. A novel geographical pattern of π suggests a model of song sparrow population history involving multiple Pleistocene refugia and colonization of some formerly glaciated regions from multiple sources. Moreover, application of coalescence theory to the haplotype cladogram suggests that two different haplotypes (48NF and 151HA) may have predominated in different parts of the song sparrow's range. This model provides insight into the current distribution of song sparrow mtDNA haplotypes and may explain the discordance between evolutionary history inferred from mtDNA and morphology in this species.     

Genetic diversity and population structure in two species of sea cucumber: differing patterns according to mode of development

The population structure of two species of sea cucumber was examined based on mitochondrial DNA sequence analysis. Cucumaria miniata , a species with pelagic nonfeeding larvae lasting less than 2 weeks, and C. pseudocurata , a brooding species lacking a pelagic phase, both occur over similar wide ranges of the northeastern Pacific between Alaska and California. No significant genetic structuring was observed among C. miniata samples with 95% of the observed nucleotide variance attributable to that within population samples. Conversely, only 3.4% of the observed variance was attributed to that within C. pseudocurata population samples, with sampling sites typically greater than 100 km apart. At a finer scale, two C. pseudocurata population samples taken 5 km apart were not statistically different. A significant genetic disjunction was observed among populations of the brooding species, but not in species with pelagic larvae, north of Vancouver Island, Canada, corresponding to the splitting of the California and Alaska currents. Given the observed high levels of genetic diversity in northern samples, this genetic break indicates survival in northern, as well as southern, Pleistocene refugia. These results clearly demonstrate the effects that changes in life-history patterns can have on dispersal, population structure, and the potential for speciation events.     

History rather than hybridization determines population structure and adaptation in Populus balsamifera

Hybridization between species is known to greatly affect their genetic diversity and, therefore, their evolution. Also, within species, there may be genetic clusters between which gene flow is limited, which may impact natural selection. However, few studies have looked simultaneously at the influence of among‐species and within‐species gene flow. Here, we study the influence of hybridization between Populus balsamifera and Populus trichocarpa on population structure and adaptation in P. balsamifera. We did this by sampling a total of 1517 individuals from across the ranges of these two species, and by genotyping them using a combination of 93 nuclear and 17 cpDNA SNPs. We found that hybridization is mostly limited to the contact zone where the species’ distributions overlap. Within P. balsamifera, we found multiple levels of population structure. Interestingly, the border between the Eastern and Central clusters is very sharp, whereas the border between the Central and Western clusters is diffuse. Outlier analysis revealed that three loci associated with the sharp border were also associated with climate. We hypothesize that the observed clusters derive from three refugia during the Pleistocene ice ages. Between the Central and Western clusters, post‐glacial long‐distance gene flow has led to the diffusion of their border. In the Eastern cluster, we hypothesize that endogenous genomic barriers have developed, leading to the sharp border and a spurious climate association. We conclude that the large‐scale genetic structure of P. balsamifera is mostly shaped by historical factors and the influence of interspecific hybridization is limited.     

Salicaceae detoxification abilities in Florida tiger swallowtail butterflies (Papilio glaucus maynardi Gauthier): Novel ability or Pleistocene holdover?

Abstract Florida populations of the eastern tiger swallowtail butterfly, Papilio glaucus L., have unique morphological features and ecological adaptations that have contributed to their subspecies status (P. g. maynardi Gauthier). We describe geographically unique abilities for detoxification of Carolina willow, Salix caroliniana Michx. (Salicaceae), for several Florida populations of P. g. maynardi. Of all the approximately 570 worldwide species of the Papilionidae, such Salicaceae detoxification abilities exist only in the allopatric North American western and northernmost species (P. rutulus Lucas, P. eurymedon Lucas and P. canadensis Rothschild & Jordan). Females of P. glaucus collected from populations in southeastern USA were examined for oviposition preference in 5‐choice assays, and displayed a low preference for Salicaceae (<5%), but larvae from Florida populations exhibited a high survival (>60%) on these plants. Detoxification abilities have previously shown to be autosomally inherited, and can be transferred via natural or hand‐paired interspecific hybrid introgression. However, these Florida populations are at least 700–1 500 km from the nearest hybrids or the hybrid species, P. appalachiensis Pavulaan & Wright, which possess these detoxification abilities. In any case, the Z ( = X)‐linked oviposition preferences for Salicaceae are lacking in these Florida populations, illustrating genetic independence of oviposition preference determination and larval survival/performance abilities. The origins of detoxification abilities are unlikely to be due to recent climate‐driven introgression, and may represent ancestral trait carry‐overs from interglacial refugium populations of the Pleistocene epoch.     

Flow- and substratum-mediated movement by a stream insect

1. In streams subject to frequent hydrologic disturbance, the ability of benthic invertebrates to disperse within the channel is key to understanding the mechanisms of flow refugium use and population persistence. This study focuses on crawling invertebrates, the effects on movement of abiotic factors (specifically, flow near the stream bed and bed micro‐topography) and the consequences for dispersal. 2. In a large flume, we observed individual cased caddisflies, Potamophylax latipennis, moving in fully turbulent flows over a precise replica of a water‐worked surface. From maps of movement paths, we quantified crawling behaviour and entrainment, and the influence of bed micro‐topography. We manipulated discharge and tested its effect on movement, linear displacement and areal dispersal. The highest discharge treatment was a disturbance to the caddis; the lowest discharge was not. Crawling behaviours were used to parameterise random walk models and estimate population dispersal, and to test the effects of abiotic factors on movement. 3. Bed micro‐topography influenced crawling in several ways. Caddis spent most of their time at the junctions between proud particles and the adjacent plane bed. The frequency distribution of turn angles was bimodal, with modal values approximating the angle required to travel around median‐sized particles. Larvae generally crawled downstream, but crawling direction relative to the flow was skewed by bed micro‐topography and was not directly downstream, unlike drift. 4. Caddis crawled for most of the time and discharge affected almost every aspect of their movement. As discharge increased, caddis crawled less often, more slowly and over shorter distances; they also became entrained more frequently and over greater distances. With increased discharge, caddis spent proportionately less time at the junctions between proud particles and the adjacent plane bed, and more time on the tops and sides of proud clasts. This is curious as most entrainment occurred from the tops and sides of clasts and entrainment is generally considered to be disadvantageous during disturbances. 5. Linear displacement (drift and entrainment combined) was downstream, but the relation between total displacement and discharge was complex. Total displacement decreased at intermediate discharge as crawling decreased, but increased at high discharge as entrainment and drift played a greater role in movement. 6. Within‐stream dispersal via crawling contained elements of both a correlated random walk (we observed directional persistence in turn angles) and a biased random walk (we observed downstream bias in move direction angles) and was best described as a biased correlated random walk. Dispersal was inversely related to discharge, suggesting that the ability of P. latipennis to crawl into flow refugia on the streambed is reduced at high flow.     

Population genetic structure of the ectomycorrhizal fungus Laccaria sp. A resembles that of its host tree Nothofagus cunninghamii

The ectomycorrhizal fungus Laccaria sp. A is restricted to temperate rainforest of southeast Australia, associated with its host tree Nothofagus cunninghamii. Eight mitochondrial microsatellite markers were used to investigate the population genetic structure of L. sp. A across its distribution in Tasmania and Victoria. The highest allelic diversity was found in Tasmania, which appeared to contain a panmictic population, whereas the more fragmented Victorian populations were characterized by low allelic diversity and differentiation between east and west. There is evidence of glacial refugia in the west and the northeast of Tasmania, and in Victoria in the Otway Ranges and Central Highlands, with postglacial migration into the Strzelecki Ranges. Narrow host-specificity may have contributed to the presence of population structure in this fungus. Allelic diversity patterns in L. sp. A are largely congruent with diversity patterns already established in populations of its host, N. cunninghamii.     

Scaling the effects of predation and disturbance in a patchy environment

The effects of hydraulic disturbances on the impact of two predatory benthic invertebrates on their prey were examined in a stream at two distinct spatial scales. At the scale of small habitat patches (0.0625 m²), hydraulic patch type was an important determinant of the microdistribution of prey and predators. Prey abundances were similar across all patch types at baseflow, but local densities were higher in patches identified as low-flow refugia after periods of high and fluctuating flow. The microdistribution pattern of predatory larvae of a caddisfly, Plectrocnemia conspersa, was similar to that of its prey, whereas predatory larvae of an alderfly, Sialis fuliginosa, did not shift their microdistribution significantly with discharge and were always most abundant in lowflow refugia. There was little evidence of an aggregative response of predators with prey, even though both predators and prey are mobile. Both predator species showed similar patch-specific patterns of per capita consumption rates: uniform consumption rates across hydraulic patch types at low and moderate flows, but highest in flow refugia during high flows. Species-specific patterns, however, were apparent in the magnitude and direction of differences between consumption rates during disturbance events, and in comparable patches at base flow: At high flow, consumption rates for P. conspersa were exaggerated (3.9 times higher) in flow refugia but at par in other patches; for S. fuliginosa they were at par in flow refugia but reduced in other patches (up to 3.3. times lower). These differences may be related to species-specific foraging behaviours (search vs ambush predators) and the influence of prey movements on feeding success. Using the patch-scale results only, it is difficult to predict the effects of physical disturbance on predation intensity at the larger scales of whole habitats, populations or communities. At the large scale (>200 m²), net predator impacts were estimated over the stream reach, using a spatially explicit model that accounts, in an additive way, for habitat heterogeneity and patch-specific responses of predators and prey. The relationship between predator impact over the whole reach and hydraulic disturbance differed for the two predators. The predator impact of S. fuliginosa decreased with increasing hydraulic disturbance, as predicted by the harsh-benign hypothesis. There was no directional trend for P. conspersa, however, and maximum predator impact may occur at intermediate disturbance levels. For the prey community in this stream, predation pressure from S. fuliginosa appears to fluctuate directly with the discharge hydrograph, whereas predation from P. conspersa may be more persistent. Flow refugia may play a dual role in the sructure of stream communities by preventing catastrophic mortality of animals (predators and prey) from physical forces during disturbances, and by maintaining (or perhaps increasing) predation pressure. Summing the effects of species interactions in small habitat patches to the larger scale of a whole stream reach indicates that the scale of approach influences the observed patterns and their implied underlying process.