Impacts of climate change on non-native species

Anthropogenic changes to climate and extreme weather events have already led to the introduction of non-native species (NNS) to the North Atlantic. Regional climate models predict that there will be a continuation of the current trend of warming throughout the 21st century providing enhanced opportunities for NNS at each stage of the invasion process. Increasing evidence is now available to show that climate change has led to the northwards range expansion of a number of NNS in the UK and Ireland, such as the Asian club tunicate Styela clava and the Pacific oyster Crassostrea gigas. Providing definitive evidence though of the direct linkage between climate change and the spread of the majority of NNS is extremely challenging, due to other confounding factors, such as anthropogenic activity. Localised patterns of water movement and food supply may also be complicating the overall pattern of northwards range expansion, by preventing the expansion of some NNS, such as the slipper limpet Crepidula fornicata and the Chilean oyster Ostrea chilensis, from a particular region. A greater understanding of the other aspects of climate change and increased atmospheric CO2, such as increased rainfall, heat waves, frequency of storm events, and ocean acidification may aid in increasing the confidence that scientists have in predicting the long term influence of climate change on the introduction, spread and establishment of NNS.

Contrasting transcriptome response to thermal stress in two key zooplankton species, Calanus finmarchicus and C. glacialis

Climate change has already led to the range expansion of warm-water plankton assemblages in the northeast Atlantic and the corresponding range contraction of colder-water species. The temperate copepod Calanus finmarchicus is predicted to shift farther northward into polar waters traditionally dominated by the arctic copepod C. glacialis. To identify temperaturemediated changes in gene expression that may be critical for the thermal acclimation and resilience of the 2 Calanus spp., we conducted a whole transcriptome profiling using RNA-seq on an Ion Torrent platform. Transcriptome responses of C. finmarchicus and C. glacialis from Disko Bay, west Greenland, were investigated under realistic thermal stresses (at + 5, +10 and +15°C) for 4 h and 6 d. C. finmarchicus showed a strong response to temperature and duration of stress, involving up-regulation of genes related to protein folding, transcription, translation and metabolism. In sharp contrast, C. glacialis displayed only low-magnitude changes in gene expression in response to temperature and duration of stress. Differences in the thermal responses of the 2 species, particularly the lack of thermal stress response in C. glacialis, are in line with laboratory and field observations and suggest a vulnerability of C. glacialis to climate change.

The Genetics of Postglacial Colonization and Peripherality in Northwestern Populations of the Spring Peeper, Pseudacris crucifer

Glaciation over the Pleistocene induced dramatic range fluctuations for species across North America such that postglacial recolonization by southern refugial lineages has characterized the genetic structure of northern North American species. Based on the leading edge model of postglacial range expansion, dispersal and rapid population growth in these northern taxa is expected to produce vast areas of genetic homogeneity. Previous work on the widely distributed spring peeper (Pseudacris crucifer) revealed six distinct mitochondrial lineages that diverged between 3-11 mya, expanding and contracting with glacial cycles. Beginning 16,000 yBP, receding glaciers permitted Eastern lineage refugia residing in the southern Appalachians to migrate northward into the St. Lawrence Valley then westward through most of central Canada. Peripheral populations at the northwestern range limit of P. crucifer in central Manitoba are likely descended from this westward expanding Eastern lineage. According to the central-marginal hypothesis, founder effects from colonization as well as limited gene flow is expected to reveal genetic differentiation and lower genetic diversity in peripheral populations. The goal of my study is to further our understanding of peripheral range dynamics in peripheral Manitoba populations of P. crucifer by determining their genetic affinity and diversity relative to more central populations in Ontario and Minnesota. In this study I amplified and aligned cytochrome b sequences from sample sites across central Manitoba to reconstruct a Bayesian phylogeny for P. crucifer; additionally, microsatellite loci were genotyped to estimate genetic diversity. Results from this study affirmed Eastern lineage descent for peripheral Manitoba sites by aligning with Ontario. Initial colonization by the Interior lineage between glacial retreat and the appearance of arid vicariance events may explain the apparent introgression of non-Eastern lineages in Manitoba. However, genetic diversity measured in expected heterozygosity (H¬e) was not found to be significantly different in Manitoba genotypes. Greater isolation by distance and inbreeding relative to Ontario and Minnesota is likely the primary driver of genetic variation in these sites. Further sampling is necessary to generate a more complete genetic population structure for P. crucifer.

Factors influencing the distribution of native and introduced Gammarus spp. in Irish river systems

Populations of Gammarus duebeni celticus, previously the only amphipod species resident in the rivers of the Lough Neagh catchment, N. Ireland, have been subjected to invasion by G. pulex from the British mainland. Numerous previous studies have investigated the potential behavioural mechanisms, principally differential mutual predation, underlying the replacement of G. d. celticus by G. pulex in Irish waters, and the mutually exclusive distributions of these species in Britain and mainland Europe. However, the relative degree of influence of abiotic versus biotic factors in structuring these amphipod communities remains unresolved. This study used principal component analysis (PCA) to distinguish physico-chemical parameters that have significant roles in determining the current distribution of G. pulex relative to G. d. celticus in L. Neagh rivers. We show that the original domination of rivers by the native G. d, celticus has changed radically, with many sites in several rivers containing either both species or only G. pulex. G. pulex was more abundant than the G. d. celticus in sites with low dissolved oxygen levels. This was reflected in the macroinvertebrate assemblages associated with G. pulex in these sites, which tended to be those tolerant of low biological water quality. The present study thus emphasizes the importance of the habitat template, particularly water quality, for Gammarus spp. interactions. If rivers become increasingly stressed by organic pollution, it is probable the range expansion of G. pulex will continue. Because these two species are not ecological equivalents, the outcomes of G. pulex incursions into G. d. celticus sites may ultimately depend on the prevailing physico-chemical regimes in each site.

Potential roles for differential body size and microhabitat complexity in mediating biotic interactions within invasive freshwater amphipod assemblages

Cannibalism and intraguild predation (IGP) are common amongst freshwater amphipod crustacean aswsemblages, particularly between individuals of different body size, with IGP of smaller by larger species. The decline of Gammarus tigrinus Populations in mainland Europe has been accompanied by the arrival of the Ponto-Caspian invader Dikerogammarus villosus and previous studies have implicated IGP of G. tigrinus by the larger D. villosus as the principal driving force in this replacement. We examined how factors such as microhabitat and body size may mediate both cannibalism within G. tigrinus populations and IGP by D. villosus and thus contribute to field patterns of coexistence and exclusion. A field Survey of an invaded Dutch fake indicated that G. tigrinus and D. villosus differed in distribution. with D. villosus being the numerically dominant amphipod (80-96 %) on the rocky boulder Substrate of the shoreline and G. tigrinus being the dominant amphipod (100 %) in the crushed shell/sand matrix immediately adjacent to this. Laboratory microcosm experiments indicated that G. tigrinus cannibalism, particularly of smaller by larger size classes, may be common. In addition, although D. villosus predation of all G. tigrinus size classes was extreme, the smallest size classes Suffered the highest predation. Indeed, when exposed to D. villosus, predation of larger G. tigrinus was lowest when smaller G. tigrinus were also present. Increasing microhabitat complexity from a simple bare substrate littered with Dreissena polymorpha zebra mussels to a Crushed shell/sand matrix significantly reduced both cannibalism and IGP. Our Study emphasizes the need to consider both life history stages and habitat template, when considering the impacts of biotic interactions and it also emphasizes that complex, interacting factors may be mediating the range expansion of D. villosus.

Fronts, jumps and secondary introductions suggested as different invasion patterns in marine species, with an increase in spread rates over time

Not all introduced (invasive) species in a region will spread from a single point of introduction. Long-distance dispersal or further introductions can obscure the pattern of spread, but the regional importance of such processes is difficult to gauge. These difficulties are further compounded when information on the multiple scale process of invasive species range expansion is reduced to one-dimensional estimates of spread (e. g. km yr(-1)). We therefore compared the results of two different metrics of range expansion: maximum linear rate of spread and accumulation of occupied grid squares (50 x 50 km) over time. An analysis of records for 54 species of introduced marine macrophytes in the Mediterranean and northeast Atlantic revealed cases where the invasion process was probably missed (e. g. Atlantic Bonnemaisonia hamifera) and suggested cases of secondary introductions or erratic jump dispersal (Dasysiphonia sp. and Womersleyella setacea). A majority of species analysed showed evidence for an accumulation of invaded sites without a clear invasion front. Estimates of spread rate are increasing for more recent introductions. The increase is greater than can be accounted for by temporally varying search effort and implies a historical increase in vector efficiency and/or a decreased resistance of native communities to invasion.

Distribution and status of the Red Squirrel (Sciurus vulgaris) and Grey Squirrel (Sciurus carolinensis) in Ireland

The distributions of the Red and Grey Squirrel were surveyed in Northern Ireland and the Republic of Ireland between 1994 and 1996. Survey methods differed between the two studies. In the former, all suitable habitat, of at least 15 ha, was inspected for species presence or absence. In the Republic, data were gathered through questionnaires to governmental and independent wildlife bodies. The combined results indicate that the Red Squirrel remains widespread and locally abundant, and is present in all but two counties. The Grey Squirrel is now more widespread than ever before, and can be found in 22 of the 32 counties. Its range expansion has varied from 0 km/yr to an estimated 13.4 km/yr, as various geographical features, principally rivers, have hindered its progress in certain directions.

The Holocene vegetation cover of Britain and Ireland: overcoming problems of scale and discerning patterns of openness

The vegetation of Europe has undergone substantial changes during the course of the Holocene epoch, resulting from range expansion of plants following climate amelioration, competition between taxa and disturbance through anthropogenic activities. Much of the detail of this pattern is understood from
decades of pollen analytical work across Europe, and this understanding has been used to address questions relating to vegetation-climate feedback, biogeography and human impact. Recent advances in modelling the relationship between pollen and vegetation now make it possible to transform pollen
proportions into estimates of vegetation cover at both regional and local spatial scales, using the Landscape Reconstruction Algorithm (LRA), i.e. the REVEALS (Regional Estimates of VEgetation Abundance from Large Sites) and the LOVE (LOcal VEgetation) models. This paper presents the compilation and analysis of 73 pollen stratigraphies from the British Isles, to assess the application of the LRA and describe the pattern of landscape/woodland openness (i.e. the cover of low herb and bushy vegetation) through the Holocene. The results show that multiple small sites can be used as an effective replacement for a single large site for the reconstruction of regional vegetation cover. The REVEALS vegetation estimates imply that the British Isles had a greater degree of landscape/woodland openness at the regional scale than areas on the European mainland. There is considerable spatial bias in the British Isles dataset towards wetland areas and uplands, which may explain higher estimates of landscape openness compared with Europe. Where multiple estimates of regional vegetation are available from within the same region inter-regional differences are greater than intra-regional differences, supporting the use of the REVEALS model to the estimation of regional vegetation from pollen data.

A spatio-temporal contrast of the predatory impact of an invasive freshwater crustacean

Aim: Impacts of invasive species may vary across invasion gradients, owing to trait-based sorting of individuals through dispersal: those aggregating at invasion fronts may be more aggressive and voracious. We examine, in the field and laboratory, variation in the predatory impacts of an invasive Ponto-Caspian crustacean Hemimysis anomala G.O. Sars, 1907 at two sites along a spatio-temporal gradient of invasion.

Location: Republic of Ireland.

Methods: We used reciprocal transplant field-deployed mesocosms to compare predation rates of invasion front and well-established H. anomala on natural zooplankton assemblages. In the laboratory, we measured the functional response (relationship between predation rate and prey supply) of H. anomala from both sites, for a per capita mechanistic comparison of predation efficiency. We also assessed prey selectivity of H. anomala in the mesocosm experiments to further compare feeding behaviour. Finally, we used a correlative approach to assess the community impact of H. anomala across sites, including a nearby uninvaded site, by comparing zooplankton diversities and densities.

Results: Invasion front H. anomala had higher predation rates than well-established H. anomala at high in situ zooplankton densities. Invasion front H. anomala also had higher functional responses - in particular showing higher 'attack rates' - indicating a heightened ability to locate and capture prey. Prey selectivity was consistent across the spatio-temporal contrast, with positive selection for cladocerans. Zooplankton diversity and density declined with time since H. anomala invasion, both being maximal at the uninvaded site.

Main conclusions: Our study, for the first time, (1) reveals differences in predatory per capita effects and associated behavioural traits between two sites along a spatio-temporal invasion gradient and (2) shows a negative community-level impact of the invasive H. anomala in natural water bodies. Further spatio-temporal comparisons of predatory per capita effects of invaders are needed to assess the generality of these results.

Evolution of the mating system in colonizing plants.

Colonization is likely to be more successful for species with an ability to self-fertilize and thus to establish new populations as single individuals. As a result, self-compatibility should be common among colonizing species. This idea, labelled 'Baker's law', has been influential in discussions of sexual-system and mating-system evolution. However, its generality has been questioned, because models of the evolution of dispersal and the mating system predict an association between high dispersal rates and outcrossing rather than selfing, and because of many apparent counter examples to the law. The contrasting predictions made by models invoking Baker's law versus those for the evolution of the mating system and dispersal urges a reassessment of how we should view both these traits. Here, I review the literature on the evolution of mating and dispersal in colonizing species, with a focus on conceptual issues. I argue for the importance of distinguishing between the selfing or outcrossing rate and a simple ability to self-fertilize, as well as for the need for a more nuanced consideration of dispersal. Colonizing species will be characterized by different phases in their life pattern: dispersal to new habitat, implying an ecological sieve on dispersal traits; establishment and a phase of growth following colonization, implying a sieve on reproductive traits; and a phase of demographic stasis at high density, during which new trait associations can evolve through local adaptation. This dynamic means that the sorting of mating-system and dispersal traits should change over time, making simple predictions difficult.

Extreme host plant conservatism during at least 20 million years of host plant pursuit by oak gallwasps

Diversification of insect herbivores is often associated with coevolution between plant toxins and insect countermeasures, resulting in a specificity that restricts host plant shifts. Gall inducers, however, bypass plant toxins and the factors influencing host plant associations in these specialized herbivores remain unclear. We reconstructed the evolution of host plant associations in Western Palaearctic oak gallwasps (Cynipidae: Cynipini), a species-rich lineage of specialist herbivores on oak (Quercus). (1) Bayesian analyses of sequence data for three genes revealed extreme host plant conservatism, with inferred shifts between major oak lineages (sections Cerris and Quercus) closely matching the minimum required to explain observed diversity. It thus appears that the coevolutionary demands of gall induction constrain host plant shifts, both in cases of mutualism (e.g., fig wasps, yucca moths) and parasitism (oak gallwasps). (2) Shifts between oak sections occurred independently in sexual and asexual generations of the gallwasp lifecycle, implying that these can evolve independently. (3) Western Palaearctic gallwasps associated with sections Cerris and Quercus diverged at least 20 million years ago (mya), prior to the arrival of oaks in the Western Palaearctic from Asia 5-7 mya. This implies an Asian origin for Western Palaearctic gallwasps, with independent westwards range expansion by multiple lineages.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximate to 141 thousand years ago (Kya), an exit out-of-Africa approximate to 51 Kya, and a recent colonization of the Americas approximate to 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

The automation and evaluation of nested clade phylogeographic analysis

Nested clade phylogeographic analysis (NCPA) is a popular method for reconstructing the demographic history of spatially distributed populations from genetic data. Although some parts of the analysis are automated, there is no unique and widely followed algorithm for doing this in its entirety, beginning with the data, and ending with the inferences drawn from the data. This article describes a method that automates NCPA, thereby providing a framework for replicating analyses in an objective way. To do so, a number of decisions need to be made so that the automated implementation is representative of previous analyses. We review how the NCPA procedure has evolved since its inception and conclude that there is scope for some variability in the manual application of NCPA. We apply the automated software to three published datasets previously analyzed manually and replicate many details of the manual analyses, suggesting that the current algorithm is representative of how a typical user will perform NCPA. We simulate a large number of replicate datasets for geographically distributed, but entirely random-mating, populations. These are then analyzed using the automated NCPA algorithm. Results indicate that NCPA tends to give a high frequency of false positives. In our simulations we observe that 14% of the clades give a conclusive inference that a demographic event has occurred, and that 75% of the datasets have at least one clade that gives such an inference. This is mainly due to the generation of multiple statistics per clade, of which only one is required to be significant to apply the inference key. We survey the inferences that have been made in recent publications and show that the most commonly inferred processes (restricted gene flow with isolation by distance and contiguous range expansion) are those that are commonly inferred in our simulations. However, published datasets typically yield a richer set of inferences with NCPA than obtained in our random-mating simulations, and further testing of NCPA with models of structured populations is necessary to examine its accuracy.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.

Host and habitat filtering in seedling root-associated fungal communities: taxonomic and functional diversity are altered in ‘novel’ soils

Climatic and land use changes have significant consequences for the distribution of tree species, both through natural dispersal processes and following management prescriptions. Responses to these changes will be expressed most strongly in seedlings near current species range boundaries. In northern temperate forest ecosystems, where changes are already being observed, ectomycorrhizal fungi contribute significantly to successful tree establishment. We hypothesised that communities of fungal symbionts might therefore play a role in facilitating, or limiting, host seedling range expansion. To test this hypothesis, ectomycorrhizal communities of interior Douglas-fir and interior lodgepole pine seedlings were analysed in a common greenhouse environment following growth in five soils collected along an ecosystem gradient. Currently, Douglas-fir’s natural distribution encompasses three of the five soils, whereas lodgepole pine’s extends much further north. Host filtering was evident amongst the 29 fungal species encountered: 7 were shared, 9 exclusive to Douglas-fir and 13 exclusive to lodgepole pine. Seedlings of both host species formed symbioses with each soil fungal community, thus Douglas-fir did so even where those soils came from outside its current distribution. However, these latter communities displayed significant taxonomic and functional differences to those found within the host distribution, indicative of habitat filtering. In contrast, lodgepole pine fungal communities displayed high functional similarity across the soil gradient. Taxonomic and/or functional shifts in Douglas-fir fungal communities may prove ecologically significant during the predicted northward migration of this species; especially in combination with changes in climate and management operations, such as seed transfer across geographical regions for forestry purposes.