Spatial–temporal analysis of species range expansion: the case of the mountain pine beetle, Dendroctonus ponderosae

Aim The spatial extent of western Canada’s current epidemic of mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae, Scolytinae), is increasing. The roles of the various dispersal processes acting as drivers of range expansion are poorly understood for most species. The aim of this paper is to characterize the movement patterns of the mountain pine beetle in areas where range expansion is occurring, in order to describe the fine‐scale spatial dynamics of processes associated with mountain pine beetle range expansion. Location Three regions of Canada’s Rocky Mountains: Kicking Horse Pass, Yellowhead Pass and Pine Pass. Methods Data on locations of mountain pine beetle‐attacked trees of predominantly lodgepole pine (Pinus contorta var. latifolia) were obtained from annual fixed‐wing aircraft surveys of forest health and helicopter‐based GPS surveys of mountain pine beetle‐damaged areas in British Columbia and Alberta. The annual (1999–2005) spatial extents of outbreak ranges were delineated from these data. Spatial analysis was conducted using the spatial–temporal analysis of moving polygons (STAMP), a recently developed pattern‐based approach. Results We found that distant dispersal patterns (spot infestations) were most often associated with marginal increases in the areal size of mountain pine beetle range polygons. When the mountain pine beetle range size increased rapidly relative to the years examined, local dispersal patterns (adjacent infestation) were more common. In Pine Pass, long‐range dispersal (> 2 km) markedly extended the north‐east border of the mountain pine beetle range. In Yellowhead Pass and Kicking Horse Pass, the extension of the range occurred incrementally via ground‐based spread. Main conclusions Dispersal of mountain pine beetle varies with geography as well as with host and beetle population dynamics. Although colonization is mediated by habitat connectivity, during periods of low overall habitat expansion, dispersal to new distant locations is common, whereas during periods of rapid invasion, locally connected spread is the dominant mode of dispersal. The propensity for long‐range transport to establish new beetle populations, and thus to be considered a driver of range expansion, is likely to be determined by regional weather patterns, and influenced by local topography. We conclude that STAMP appears to be a useful approach for examining changes in biogeograpical ranges, with the potential to reveal both fine‐ and large‐scale patterns.     

Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Environmental change has a wide range of ecological consequences, including species extinction and range expansion. Many studies have shown that insect species respond rapidly to climatic change. A mountain pine beetle epidemic of record size in North America has led to unprecedented mortality of lodgepole pine, and a significant range expansion to the northeast of its historic range. Our goal was to determine the spatial genetic variation found among outbreak population from which genetic structure, and dispersal patterns may be inferred. Beetles from 49 sampling locations throughout the outbreak area in western Canada were analysed at 13 microsatellite loci. We found significant north-south population structure as evidenced by: (i) Bayesian-based analyses, (ii) north-south genetic relationships and diversity gradients; and (iii) a lack of isolation-by-distance in the northernmost cluster. The north-south structure is proposed to have arisen from the processes of postglacial colonization as well as recent climate-driven changes in population dynamics. Our data support the hypothesis of multiple sources of origin for the outbreak and point to the need for population specific information to improve our understanding and management of outbreaks. The recent range expansion across the Rocky Mountains into the jack/lodgepole hybrid and pure jack pine zones of northern Alberta is consistent with a northern British Columbia origin. We detected no loss of genetic variability in these populations, indicating that the evolutionary potential of mountain pine beetle to adapt has not been reduced by founder events. This study illustrates a rapid range-wide response to the removal of climatic constraints, and the potential for range expansion of a regional population.     

Population structure and migration pattern of a conifer pathogen, Grosmannia clavigera, as influenced by its symbiont, the mountain pine beetle

We investigated the population structure of Grosmannia clavigera (Gc), a fungal symbiont of the mountain pine beetle (MPB) that plays a crucial role in the establishment and reproductive success of this pathogen. This insect-fungal complex has destroyed over 16 million ha of lodgepole pine forests in Canada, the largest MPB epidemic in recorded history. During this current epidemic, MPB has expanded its range beyond historically recorded boundaries, both northward and eastward, and has now reached the jack pine of Alberta, potentially threatening the Canadian boreal forest. To better understand the dynamics between the beetle and its fungal symbiont, we sampled 19 populations in western North America and genotyped individuals from these populations with eight microsatellite markers. The fungus displayed high haplotype diversity, with over 250 unique haplotypes observed in 335 single spore isolates. Linkage equilibria in 13 of the 19 populations suggested that the fungus reproduces sexually. Bayesian clustering and distance analyses identified four genetic clusters that corresponded to four major geographical regions, which suggested that the epidemic arose from multiple geographical sources. A genetic cluster north of the Rocky Mountains, where the MPB has recently become established, experienced a population bottleneck, probably as a result of the recent range expansion. The two genetic clusters located north and west of the Rocky Mountains contained many fungal isolates admixed from all populations, possibly due to the massive movement of MPB during the epidemic. The general agreement in north-south differentiation of MPB and G. clavigera populations points to the fungal pathogen's dependence on the movement of its insect vector. In addition, the patterns of diversity and the individual assignment tests of the fungal associate suggest that migration across the Rocky Mountains occurred via a northeastern corridor, in accordance with meteorological patterns and observation of MPB movement data. Our results highlight the potential of this pathogen for both expansion and sexual reproduction, and also identify some possible barriers to gene flow. Understanding the ecological and evolutionary dynamics of this fungus-beetle association is important for the modelling and prediction of MPB epidemics.     

Landscape level analysis of mountain pine beetle in British Columbia, Canada: spatiotemporal development and spatial synchrony within the present outbreak

Eruptive herbivores can exert profound landscape level influences. For example, the ongoing mountain pine beetle outbreak in British Columbia, Canada, has resulted in mortality of mature lodgepole pine over >7 million ha. Analysis of the spatio‐temporal pattern of spread can lend insights into the processes initiating and/or sustaining such phenomena. We present a landscape level analysis of the development of the current outbreak. Aerial survey assessments of tree mortality, projected onto discrete 12×12 km cells, were used as a proxy for insect population density. We examined whether the outbreak potentially originated from an epicenter and spread, or whether multiple localized populations erupted simultaneously at spatially disjunct locations. An aspatial cluster analysis of time series from 1990 to 2003 revealed four distinct time series patterns. Each time series demonstrated a general progression of increasing mountain pine beetle populations. Plotting the geographical locations of each temporal pattern revealed that the outbreak occurred first in an area of west‐central British Columbia, and then in an area to the east. The plot further revealed many localized infestations erupted in geographically disjunct areas, especially in the southern portion of the province. Autologistic regression analyses indicated a significant, positive association between areas where the outbreak first occurred and conservation lands. For example, the delineated area of west‐central British Columbia is comprised of three conservation parks and adjacent working forest. We further examined how population synchrony declines with distance at different population levels. Examination of the spatial dependence of temporal synchrony in population fluctuations during early, incipient years (i.e. 1990–1996) suggested that outbreaking mountain pine beetle populations are largely independent at scales >200 km during non‐epidemic periods. However, during epidemic years (i.e. 1999–2003), populations were clearly synchronous across the entire province, even at distances of up to 900 km. The epicentral pattern of population development can be used to identify and prioritize adjacent landscape units for both reactive and proactive management strategies intended to minimize mountain pine beetle impacts.     

Population Structure of Mountain Pine Beetle Symbiont Leptographium longiclavatum and the Implication on the Multipartite Beetle-Fungi Relationships

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.     

A spatiotemporal pattern analysis of historical mountain pine beetle outbreaks in British Columbia,Canada

This study documents the spatiotemporal patterns of mountain pine beetle infestations by applying a novel approach based on a landscape infestation dynamics conceptual model in combination with morphological spatial pattern analysis using the mountain pine beetle infested pine mortality data (1960–2010) collected by the annual British Columbia aerial overview survey. The pattern analysis at the provincial level reveals that the 1980s outbreak did not crash as originally thought. The current outbreak is most likely a result of the progressive buildup of the epidemic infestations during the transition period (1985–1995) under favourable weather conditions and substantially improved host resources. This is also true for the Northeast and Cariboo areas of the province specifically, even though the infestations in the Cariboo area remained at incipient‐epidemic levels during the transition period after the 1980s outbreak crashed in 1985. In the Southeast area, the current outbreak apparently continued from the outbreak that initiated in the late 1970s and early 1980s. The 1980s outbreak originated from multiple spatially separate locations whereas the current outbreak initiated from a single location in Tweedsmuir Provincial Park. The centralized and self‐amplifying buildup of the current outbreak implicates at least three substantial expansions that occurred in 2002, 2006, and 2008. This study suggests that at the provincial level, as well as for the Northeast and Southeast areas of the province, the current outbreak is declining but most likely will continue for many years given the ongoing and future warming climate and a large proportion of pines that remain in the habitats of mountain pine beetles. This study also suggests that dispersals, particularly long‐distance dispersal, may play a key role in driving the spread and expansion of the current outbreak although uncertainty remains due to the local dynamics of the beetle populations.     

Modeling the Expansion of an Introduced Tree Disease

Pine wilt disease is caused by the introduced pinewood nematode, Bursaphelenchus xylophilus, for which the vector is the pine sawyer beetle, Monochamus alternatus. Native Japanese pines, black pine (Pinus thunbergii) and red pine (P. densiflora), are extremely sensitive to the nematode's infection, and the parasite has been expanding nationwide in the last few decades, despite intensive control efforts. To understand the parasite's range expansion in Japan, we modeled the dynamics of the pines and the beetle that disperses the nematode, using an integro-difference equation in a one-dimensional space. Based on field data collected in Japan, we investigated the dependence of the parasite's rate of range expansion on the eradication rate of the beetle, the initial pine density, and the beetle dispersal ability. Our model predicts several results. (1) The Allee Effect operates on beetle reproduction, and consequently the parasite cannot invade a pine stand, once the beetle density decreases below a threshold. (2) The distribution of the dispersal distance of the beetles critically affects the expansion rate of the disease. As the fraction of the beetles that travel over long distance increases from zero, the range expansion accelerates sharply. (3) However, too frequent long-range dispersal results in a failure of the parasite invasion due to the Allee Effect, suggesting the importance of correctly assessing the beetle's mobility to predict the speed of range expansion of the parasite. (4) As the eradication rate is increased, the range expansion speed decreases gradually at first and suddenly drops to zero at a specific value of the eradication rate. This revised version was published online in July 2006 with corrections to the Cover Date.     

Molecular evidence for Pleistocene glacial cycles driving diversification of a North American desert spider, Agelenopsis aperta

The influence of historical climatic vs. geological changes on species diversification patterns was investigated in a widely distributed North American desert spider, Agelenopsis aperta (Araneae: Agelenidae), with particular reference to Pleistocene glacial cycles and earlier patterns of mountain building. Levels of sequence divergence obtained from the mitochondrial gene, cytochrome oxidase I, dated to the Pleistocene, eliminating Rocky Mountain orogeny as a cause of diversification, as orogeny ended 4 million years ago. The results of phylogenetic and network analyses showed the presence of three geographically defined clades, which were consistent with the presence of at least three glacial refugia: (i) east of the Rocky Mountains; (ii) between the Rocky Mountains and Sierra Nevadas; and (iii) west of the Sierra Nevadas. In addition, populations within the Rocky Mountains exhibited significantly lower genetic diversity than populations east of the Rocky Mountains and the haplotypes found within the Rockies were a subset of eastern haplotypes. These patterns suggest that a post-Pleistocene range expansion occurred out of an eastern glacial refugium into the Rocky Mountains. Examination of phylogeographical studies of other North American desert taxa indicated that mountain building explained diversification patterns more effectively for some taxa but Pleistocene climate change was more important for others, including A. aperta.     

Influences of the biophysical environment on blister rust and mountain pine beetle,and their interactions,in whitebark pine forests

Aim To understand how the biophysical environment influences patterns of infection by non‐native blister rust (caused by Cronartium ribicola) and mortality caused by native mountain pine beetles (Dendroctonus ponderosae) in whitebark pine (Pinus albicaulis) communities, to determine how these disturbances interact, and to gain insight into how climate change may influence these patterns in the future. Location High‐elevation forests in south‐west Montana, central Idaho, eastern and western Oregon, USA. Methods Stand inventory and dendroecological methods were used to assess stand structure and composition and to reconstruct forest history at sixty 0.1‐ha plots. Patterns of blister rust infection and mountain pine beetle‐caused mortality in whitebark pine trees were examined using nonparametric Kruskal–Wallis ANOVA, Mann–Whitney U‐tests, and Kolmogorov–Smirnov two‐sample tests. Stepwise regression was used to build models of blister rust infection and mountain pine beetle‐related mortality rates based on a suite of biophysical site variables. Results Occurrence of blister rust infections was significantly different among the mountain ranges, with a general gradient of decreasing blister rust occurrence from east to west. Evidence of mountain pine beetle‐caused mortality was identified on 83% of all dead whitebark pine trees and was relatively homogenous across the study area. Blister rust infected trees of all ages and sizes uniformly, while mountain pine beetles infested older, larger trees at all sites. Stepwise regressions explained 64% and 58% of the variance in blister rust infection and beetle‐caused mortality, respectively, indicating that these processes are strongly influenced by the biophysical environment. More open stand structures produced by beetle outbreaks may increase the exposure of surviving whitebark pine trees to blister rust infection. Main conclusions Variability in the patterns of blister rust infection and mountain pine beetle‐caused mortality elucidated the fundamental dynamics of these disturbance agents and suggests that the effects of climate change will be complex in whitebark pine communities and vary across the species’ range. Interactions between blister rust and beetle outbreaks may accelerate declines or facilitate the rise of rust resistance in whitebark pine depending on forest conditions at the time of the outbreak.     

Propagule pressure and colony social organization are associated with the successful invasion and rapid range expansion of fire ants in China

We characterized patterns of genetic variation in populations of the fire ant Solenopsis invicta in China using mitochondrial DNA sequences and nuclear microsatellite loci to test predictions as to how propagule pressure and subsequent dispersal following establishment jointly shape the invasion success of this ant in this recently invaded area. Fire ants in Wuchuan (Guangdong Province) are genetically differentiated from those found in other large infested areas of China. The immediate source of ants in Wuchuan appears to be somewhere near Texas, which ranks first among the southern USA infested states in the exportation of goods to China. Most colonies from spatially distant, outlying areas in China are genetically similar to one another and appear to share a common source (Wuchuan, Guangdong Province), suggesting that long‐distance jump dispersal has been a prevalent means of recent spread of fire ants in China. Furthermore, most colonies at outlier sites are of the polygyne social form (featuring multiple egg‐laying queens per nest), reinforcing the important role of this social form in the successful invasion of new areas and subsequent range expansion following invasion. Several analyses consistently revealed characteristic signatures of genetic bottlenecks for S. invicta populations in China. The results of this study highlight the invasive potential of this pest ant, suggest that the magnitude of international trade may serve as a predictor of propagule pressure and indicate that rates and patterns of subsequent range expansion are partly determined by the interplay between species traits and the trade and transportation networks.     

Evaluating future success of whitebark pine ecosystem restoration under climate change using simulation modeling

Major declines of whitebark pine forests throughout western North America from the combined effects of mountain pine beetle (Dendroctonus ponderosae) outbreaks, fire exclusion policies, and the exotic disease white pine blister rust (WPBR) have spurred many restoration actions. However, projected future warming and drying may further exacerbate the species' decline and possibly compromise long‐term success of today's restoration activities. We evaluated successes of restoration treatments under future climate using a comprehensive landscape simulation experiment. The spatially explicit, ecological process model FireBGCv2 was used to simulate whitebark pine populations on two U.S. Northern Rocky Mountain landscapes over 95 years under two climate, three restoration, and two fire management scenarios. Major findings were that (1) whitebark pine can remain on some high mountain landscapes in a future climate albeit at lower basal areas (50% decrease), (2) restoration efforts, such as thinning and prescribed burning, are vital to ensure future whitebark pine forests, and (3) climate change impacts on whitebark pine vary by local setting. Whitebark pine restoration efforts will mostly be successful in the future but only if future populations are somewhat resistant to WPBR. Results were used to develop general guidelines that address climate change impacts for planning, designing, implementing, and evaluating fine‐scale restoration activities.     

Spatial scaling of mountain pine beetle infestations

1. The relationship between occupancy and spatial contagion during the spread of eruptive and invasive species demands greater study, as it could lead to improved prediction of ecosystem damage. 2. We applied a recently developed model that links occupancy and its fractal dimension to model the spatial distribution of mountain pine beetle infestations in British Columbia, Canada. We showed that the distribution of infestation was scale-invariant in at least 24 out of 37 years (mostly in epidemic years), and presented some degree of scale-invariance in the rest. There was a general logarithmic relationship between fractal dimension and infestation occupancy. Based on the scale-invariance assumption, we further assessed the interrelationships for several landscape metrics, such as correlation length, maximum cluster size, total edge length and total number of clusters. 3. The scale-invariance assumption allows fitting the above metrics, and provides a framework to establish the scaling relationship between occupancy and spatial contagion. 4. We concluded that scale-invariance dominates the spread of mountain pine beetle. In this context, spatial aggregation can be predicted from occupancy, hence occupancy is the only variable one needs to know in order to predict the spatial distributions of populations. This supports the hypothesis that fractal dispersal kernels may be abundant among outbreaks of pests and invasive species.     

Modelling local and long-distance dispersal of invasive emerald ash borer Agrilus planipennis (Coleoptera) in North America

Limiting the damage by non-indigenous species requires rapid determination of current and potential distributions and vectors of dispersal, and development of appropriate management measures. The emerald ash borer ( Agrilus planipennis ), a wood-boring beetle native to South-East Asia, was first reported in the Great Lakes region during summer 2002. The beetle poses an enormous threat to native ash ( Fraxinus ) species of North America, as untreated trees in infested areas of Ontario, Michigan and Ohio suffer high mortality. We demonstrate that the borer has spread in North America through a combination of diffusive range extension, associated with local flights, and by long-distance 'jump' dispersal associated with human movement of infested sapling or contaminated firewood. Probability of infestation was inversely related to distance from borer epicentres but positively related to the size of human population centres. At least 9 of 39 populations that were first reported in Michigan during 2004 cannot be accounted for by local diffusion, raising the possibility that other unidentified mechanisms may be contributing to the dispersal of the beetle. In the absence of quarantine, by 2005 all of Michigan's lower peninsula was contained within the boundaries of potential diffusive range expansion. Infested ash saplings also were introduced from Michigan to Maryland during 2003, and subsequently transplanted to five sites in Maryland and Virginia. Quarantine and eradication measures have had mixed results: in the south-central USA, the species appears on the brink of eradication, whereas its distribution has continued to spread during 2005 in the Great Lakes region despite extensive containment and quarantine measures. Quarantine success in the Great Lakes region is encumbered by multiple dispersal vectors, larger borer population sizes and by the more extensive geographical distribution that was achieved prior to implementation of control measures.     

Shifting dispersal modes at an expanding species’ range margin

While it is generally recognized that noncontiguous (long‐distance) dispersal of small numbers of individuals is important for range expansion over large geographic areas, it is often assumed that colonization on more local scales proceeds by population expansion and diffusion dispersal (larger numbers of individuals colonizing adjacent sites). There are few empirical studies of dispersal modes at the front of expanding ranges, and very little information is available on dispersal dynamics at smaller geographic scales where we expect contiguous (diffusion) dispersal to be prevalent. We used highly polymorphic genetic markers to characterize dispersal modes at a local geographic scale for populations at the edge of the range of a newly invasive grass species (Brachypodium sylvaticum) that is undergoing rapid range expansion in the Pacific Northwest of North America. Comparisons of Bayesian clustering of populations, patterns of genetic diversity, and gametic disequilibrium indicate that new populations are colonized ahead of the invasion front by noncontiguous dispersal from source populations, with admixture occurring as populations age. This pattern of noncontiguous colonization was maintained even at a local scale. Absence of evidence for dispersal among adjacent pioneer sites at the edge of the expanding range of this species suggests that pioneer populations undergo an establishment phase during which they do not contribute emigrants for colonization of neighbouring sites. Our data indicate that dispersal modes change as the invasion matures: initial colonization processes appear to be dominated by noncontiguous dispersal from only a few sources, while contiguous dispersal may play a greater role once populations become established.     

Population genetic data suggest a role for mosquito‐mediated dispersal of West Nile virus across the western United States

After introduction, West Nile virus (WNV) spread rapidly across the western United States between the years 2001 and 2004. This westward movement is thought to have been mediated by random dispersive movements of resident birds. Little attention has been placed on the role of mosquito vectors in virus dispersal across North America. The mosquito vector largely responsible for WNV amplification and transmission of WNV in the western USA is Culex tarsalis. Here we present population genetic data that suggest a potential role for C. tarsalis in the dispersal of WNV across the western USA. Population genetic structure across the species range of C. tarsalis in the USA was characterized in 16 states using 12 microsatellite loci. structure and geneland analyses indicated the presence of three broad population clusters. Barriers to gene flow were resolved near the Sonoran desert in southern Arizona and between the eastern Rocky Mountains and High Plains plateau. Small genetic distances among populations within clusters indicated that gene flow was not obstructed over large portions of the West Coast and within the Great Plains region. Overall, gene flow in C. tarsalis appears to be extensive, potentially mediated by movement of mosquitoes among neighbouring populations and hindered in geographically limited parts of its range. The pattern of genetic clustering in C. tarsalis is congruent with the pattern of invasion of WNV across the western United States, raising the possibility that movement of this important vector may be involved in viral dispersal.     

Influence of ecological and geological features on rangewide patterns of genetic structure in a widespread passerine

Geological and ecological features restrict dispersal and gene flow, leading to isolated populations. Dispersal barriers can be obvious physical structures in the landscape; however microgeographic differences can also lead to genetic isolation. Our study examined dispersal barriers at both macro- and micro-geographical scales in the black-capped chickadee, a resident North American songbird. Although birds have high dispersal potential, evidence suggests dispersal is restricted by barriers. The chickadee''s range encompasses a number of physiological features which may impede movement and lead to divergence. Analyses of 913 individuals from 34 sampling sites across the entire range using 11 microsatellite loci revealed as many as 13 genetic clusters. Populations in the east were largely panmictic whereas populations in the western portion of the range showed significant genetic structure, which often coincided with large mountain ranges, such as the Cascade and Rocky Mountains, as well as areas of unsuitable habitat. Unlike populations in the central and southern Rockies, populations on either side of the northern Rockies were not genetically distinct. Furthermore, Northeast Oregon represents a forested island within the Great Basin; genetically isolated from all other populations. Substructuring at the microgeographical scale was also evident within the Fraser Plateau of central British Columbia, and in the southeast Rockies where no obvious physical barriers are present, suggesting additional factors may be impeding dispersal and gene flow. Dispersal barriers are therefore not restricted to large physical structures, although mountain ranges and large water bodies do play a large role in structuring populations in this study.     

Spread of the pinewood nematode vectored by the Japanese pine sawyer: modeling and analytical approaches

The pinewood nematode, Bursaphelenchus xylophilus, is the causative agent of pine wilt of Pinus thunbergii and P. densiflora in Japan. The nematode is vectored by cerambycid beetles of the genus Monochamus. It is inferred to have been introduced from North America early in the 1900s and then to have distributed in China, Korea, and Taiwan. Intensive and/or long-term studies of pine wilt systems have elucidated the pattern and mechanism of the nematode’s spread within a pine stand, dispersal of vector beetles, and spread pattern of pine wilt within a prefecture. The modeling of nematode spread over pine stands, which involves beetle reproduction within a pine stand, has been developing and should elucidate the factors influencing the rate at which the nematode range expands. In this review, we summarize the biologies of the nematode, beetle, and tree, and then characterize the spread of the nematode within a pine stand, locally over pine stands, and regionally over unit administrative districts. Local and regional spreading of the nematode is related primarily to long-distance dispersal by insect vectors and to the artificial transportation of pine logs infested with the nematode and its vector, respectively.     

Molecular genetic evidence for the post-Pleistocene divergence of populations of the arctic-alpine ground beetle Amara alpina (Paykull) (Coleoptera: Carabidae)

Aim This study examines the hypothesis that the biogeographic history of a species is reflected in the distribution of molecular genetic diversity and the phylogenies of extant populations. Location Populations of arctic-alpine ground beetle Amara alpina were analysed from Beringia (Alaska and northernmost British Columbia), the Hudson Bay region, the northern Appalachian Mountains, and the central Rocky Mountains of North America. Methods Mitochondrial restriction site variation of specimens from twenty-two populations were assayed by using radioactively labelled mtDNA to probe Southern membranes containing restriction enzyme digested total DNA. Restriction sites were mapped and genetic distances were calculated by pairwise comparison of presence and absence of restriction sites. Genetic distances were used in a molecular analysis of variance and to construct a minimal spanning tree. Parsimony methods were used to investigate the phylogenetic relationships between the haplotypes. These results were compared to an existing model for postglacial dispersal based on fossil and modern occurrences of arctic-alpine beetles. Results Among the twenty-two populations, fifteen haplotypes were detected. Genetic variation within each of the four regions corresponded to that expected from the palaeontologically based model. Beringian populations were the most genetically diverse. In contrast, no restriction site variation was observed in populations from the Hudson Bay region. Intermediate amounts of variation were observed in alpine populations of the Rocky and Appalachian Mountains. Maximum parsimony and cluster analysis provide evidence that at least two ancestral haplotypes existed in the Southern refugium from which the Rocky and the Appalachian Mountains populations were founded. Main conclusions The genetic results are generally consistent with the palaeontologically based model. The diversity of Beringian populations is consistent with this region having been continuously inhabited by Amara alpina throughout the Pleistocene. The Hudson Bay region was not deglaciated until about 6000 years, and its populations have no restriction site variation. The molecular genetic data support the interpretation that the Hudson Bay region was colonized from Beringia based on the occurrence of the same haplotype in both regions.     

Adaptive and neutral markers both show continent‐wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.     

The influence of a competitor on the geographic mosaic of coevolution between crossbills and lodgepole pine

The geographic mosaic theory of coevolution posits that the form of selection between interacting species varies across a landscape with coevolution important and active in some locations (i.e., coevolutionary hotspots) but not in others (i.e., coevolutionary coldspots). We tested the hypothesis that the presence of red squirrels (Tamiasciurus hudsonicus) affects the occurrence of coevolution between red crossbills (Loxia curvirostra complex) and Rocky Mountain lodgepole pine (Pinus contorta ssp. latifolia) and thereby provides a mechanism giving rise to a geographic mosaic of selection. Red squirrels are the predominant predispersal seed predator and selective agent on lodgepole pine cones. However, in four isolated mountain ranges east and west of the Rocky Mountains, red squirrels are absent and red crossbills are the main predispersal seed predator. These isolated populations of pine have apparently evolved without Tamiasciurus for about 10,000 to 12,000 years. Based on published morphological, genetic, and paleobotanical studies, we infer that cone traits in these isolated populations that show parallel differences from cones in the Rocky Mountains have changed in parallel. We used data on crossbill and conifer cone morphology and feeding preferences and efficiency to detect whether red crossbills and lodgepole pine exhibit reciprocal adaptations, which would imply coevolution. Cone traits that act to deter Tamiasciurus and result in high ratios of cone mass to seed mass were less developed in the isolated populations. Cone traits that act to deter crossbills include larger and thicker scales and perhaps increased overlap between successive scales and were enhanced in the isolated populations. In the larger, isolated mountain ranges crossbills have evolved deeper, shorter, and therefore more decurved bills to exploit these cones. This provides crossbills with higher feeding rates, and the change in bill shape has improved efficiency by reducing the concomitant increases in body mass and daily energy expenditures that would have resulted if only bill size had increased. These parallel adaptations and counter adaptations in red crossbills and lodgepole pine are interpreted as reciprocal adaptations and imply that these crossbills and pine are in coevolutionary arms races where red squirrels are absent (i.e., coevolutionary hotspots) but not where red squirrels are present (i.e., coevolutionary cold-spots).