The imprint of the geographical, evolutionary and ecological context on species-area relationships

Species-area relationships (SAR) are fundamental in the understanding of biodiversity patterns and of critical importance for predicting species extinction risk worldwide. Despite the enormous attention given to SAR in the form of many individual analyses, little attempt has been made to synthesize these studies. We conducted a quantitative meta-analysis of 794 SAR, comprising a wide span of organisms, habitats and locations. We identified factors reflecting both pattern-based and dynamic approaches to SAR and tested whether these factors leave significant imprints on the slope and strength of SAR. Our analysis revealed that SAR are significantly affected by variables characterizing the sampling scheme, the spatial scale, and the types of organisms or habitats involved. We found that steeper SAR are generated at lower latitudes and by larger organisms. SAR varied significantly between nested and independent sampling schemes and between major ecosystem types, but not generally between the terrestrial and the aquatic realm. Both the fit and the slope of the SAR were scale-dependent. We conclude that factors dynamically regulating species richness at different spatial scales strongly affect the shape of SAR. We highlight important consequences of this systematic variation in SAR for ecological theory, conservation management and extinction risk predictions.

Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal

We examined a remnant host plant (Primula veris L.) habitat network that was last inhabited by the rare butterfly Hamearis lucina L. in north Wales in 1943, to assess the relative contribution of several spatial parameters to its regional extinction. We first examined relationships between P. veris characteristics and H. lucina eggs in surviving H. lucina populations, and used these to predict the suitability and potential carrying capacity of the habitat network in north Wales. This resulted in an estimate of roughly 4500 eggs (ca 227 adults). We developed a discrete space, discrete time metapopulation model to evaluate the relative contribution of dispersal distance, habitat and environmental stochasticity as possible causes of extinction. We simulated the potential persistence of the butterfly in the current network as well as in three artificial (historical and present) habitat networks that differed in the quantity (current and X3) and fragmentation of the habitat (current and aggregated). We identified that reduced habitat quantity and increased isolation would have increased the probability of regional extinction, in conjunction with environmental stochasticity and H. lucina's dispersal distance. This general trend did not change in a qualitative manner when we modified the ability of dispersing females to stay in, and find suitable habitats (by changing the size of the grid cells used in the model). Contrary to most metapopulation model predictions, system persistence declined with increasing migration rate, suggesting that the mortality of migrating individuals in fragmented landscapes may pose significant risks to system-wide persistence. Based on model predictions for the present landscape we argue that a major programme of habitat restoration would be required for a re-established metapopulation to persist for > 100 years.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.

The paradox of invasion: Reeves’ muntjac deer invade the British Isles from a limited number of founding females

High levels of genetic diversity and high propagule pressure are favoured by conservation biologists as the basis for successful reintroductions and ensuring the persistence of populations. However, invasion ecologists recognize the ‘paradox of invasion’, as successful species introductions may often be characterized by limited numbers of individuals and associated genetic bottlenecks. In the present study, we used a combination of high-resolution nuclear and mitochondrial genetic markers to investigate the invasion history of Reeves' muntjac deer in the British Isles. This invasion has caused severe economic and ecological damage, with secondary spread currently a concern throughout Europe and potentially globally. Microsatellite analysis based on eight loci grouped all 176 introduced individuals studied from across the species' range in the UK into one genetic cluster, and seven mitochondrial D-loop haplotypes were recovered, two of which were present at very low frequency and were related to more common haplotypes. Our results indicate that the entire invasion can be traced to a single founding event involving a low number of females. These findings highlight the fact that even small releases of species may, if ignored, result in irreversible and costly invasion, regardless of initial genetic diversity or continual genetic influx.

Genetic provenance and best practice woodland management: a case study in native alder (Alnus glutinosa)

In recent years, the native woodlands of Europe, including those of Britain and Ireland, have increasingly come under threat from a range of biotic and abiotic factors, and are therefore a conservation priority demanding careful management in order to realise their inherent ecological and cultural benefits. Because the distribution of genetic variation across populations and regions is increasingly considered an important component of woodland management, we carried out a population genetic analysis on black alder (Alnus glutinosa) across Northern Ireland in order to inform “best practice” strategies. Our findings suggest that populations harbour high levels of genetic diversity, with very little differentiation between populations. Significant F IS values were observed in over half of the populations analysed, however, which could reflect inbreeding as a result of the patchy occurrence of alder in Northern Ireland, with scattered, favourable damp habitats being largely isolated from each other by extensive tracts of farmland. Although there is no genetic evidence to support the broad-scale implementation of tree seed zones along the lines of those proposed for native woodlands in Great Britain, we suggest that the localised occurrence of rare chloroplast haplotypes should be taken into account on a case-by-case basis. This, coupled with the identification of populations containing high genetic diversity and that are broadly representative of the region as a whole, will provide a sound genetic basis for woodland management, both in alder and more generally for species that exhibit low levels of genetic differentiation.

Differential ecological impacts of an invader and a native under environmental change are revealed by comparative functional responses.

Predicting the ecological impacts of damaging invasive species under relevant environmental contexts is a major challenge, for which comparative functional responses (the relationship between resource availability and consumer uptake rate) have great potential. Here, the functional responses of Gammarus pulex, an ecologically damaging invader in freshwaters in Ireland and other islands, were compared with those of a native trophic equivalent Gammarus duebeni celticus. Experiments were conducted at two dissolved oxygen concentrations (80 and 50 % saturation), representative of anthropogenic water quality changes, using two larval prey, blackfly (Simuliidae spp.) and mayfly (Baetis rhodani). Overall, G. pulex had higher Type II functional responses and hence predatory impacts than G. d. celticus and the functional responses of both predators were reduced by lowered oxygen concentration. However, this reduction was of lower magnitude for the invader as compared to the native. Further, the invader functional response at low oxygen was comparable to that of the native at high oxygen. Attack rates of the two predators were similar, with low oxygen reducing these attack rates, but this effect occurred more strongly for blackfly than mayfly prey. Handling times were significantly lower for the invader compared with the native, and significantly higher at low oxygen, however, the effect of lowered oxygen on handling times was minimal for the invader and pronounced for the native. Maximum feeding rates were significantly greater for the invader compared with the native, and significantly reduced at low oxygen, with this effect again lesser for the invader as compared to the native. The greater functional responses of the invader corroborate with its impacts on recipient macroinvertebrate communities when it replaces the native. Further, our experiments predict that the impact of the invader will be less affected than the native under altered oxygen regimes driven by anthropogenic influences.

Use of an earthworm lysosomal biomarker for the ecological assessment of pollution from an industrial plastics fire

A low-cost field technique employing retention of the dye neutral-red by lysosomes in coelomocyte cells taken from earthworms (Lumbricus castaneus), was used as a means of assessing the ecological effects (if any) of an industrial accident. Earthworms and soil samples were collected at the site of a large industrial plastics fire in Thetford, UK along a 200 m transect leading from the factory perimeter fence, over a layer of molten plastic impregnated soil and into the surrounding forest. Coelomic fluid extracted from the earthworms was dye-loaded with neutral-red and lysosomal leaking observed. Metal residues in soil and earthworms were found to be highly elevated close to the factory perimeter and to rapidly drop to background levels within the first 50 m of the transect. Coelomocyte cells taken from earthworms adjacent to the factory perimeter showed the shortest period of neutral-red retention (2 min); cells taken from worms further into the surrounding forest had a longer retention time (12 min), whilst cells taken from worms from a control site showed even greater retention times (25 min). Thus, the neutral-red retention times correlated negatively with measured residues of heavy metals in the earthworms, the higher the body metal concentration the shorter the retention time. This field trial has demonstrated the validity of using an in vitro cellular biomarker technique for use in biological impact assessment along gradients of contamination.

Ecological communities are vulnerable to realistic extinction sequences

Loss of species will directly change the structure and potentially the dynamics of ecological communities, which in turn may lead to additional species loss (secondary extinctions) due to direct and/or indirect effects (e.g. loss of resources or altered population dynamics). Furthermore, the vulnerability of food webs to repeated species loss is expected to be affected by food web topology, species interactions, as well as the order in which species go extinct. Species traits such as body size, abundance and connectivity might determine a species' vulnerability to extinction and, thus, the order in which species go primarily extinct. Yet, the sequence of primary extinctions, and their effects on the vulnerability of food webs to secondary extinctions, when species abundances are allowed to respond dynamically, has only recently become the focus of attention. Here, we analyse and compare topological and dynamical robustness to secondary extinctions of model food webs, in the face of 34 extinction sequences based on species traits. Although secondary extinctions are frequent in the dynamical approach and rare in the topological approach, topological and dynamical robustness tends to be correlated for many bottom-up directed, but not for top-down directed deletion sequences. Furthermore, removing species based on traits that are strongly positively correlated to the trophic position of species (such as large body size, low abundance, high net effect) is, under the dynamical approach, found to be as destructive as removing primary producers. Such top-down oriented removal of species are often considered to correspond to realistic extinction scenarios, but earlier studies, based on topological approaches, have found such extinction sequences to have only moderate effects on the remaining community. Thus, our result suggests that the structure of ecological communities, and therefore the integrity of important ecosystem processes could be more vulnerable to realistic extinction sequences than previously believed.

A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Ecological differences and coexistence in a guild of microparasites: Bartonella in wild rodents

The study of ecological differences among coexisting microparasites has been largely neglected, but it addresses important and unusual issues because there is no clear distinction in such cases between conventional (resource) and apparent competition. Here patterns in the population dynamics are examined for four species of Bartonella (bacterial parasites) coexisting in two wild rodent hosts, bank voles (Clethrionomys glareolus) and wood mice (Apodemus sylvaticus). Using generalized linear modeling and mixed effects models, we examine, for these four species, seasonal patterns and dependencies on host density (both direct and delayed) and, having accounted for these, any differences in prevalence between the two hosts. Whereas previous studies had failed to uncover species differences, here all four were different. Two, B. doshiae and B. taylorii, were more prevalent in wood mice, and one, B. birtlesii, was more prevalent in bank voles. B. birtlesii, B. grahamii, and B. taylorii peaked in prevalence in the fall, whereas B. doshiae peaked in spring. For B. birtlesii in bank voles, density dependence was direct, but for B. taylorii in wood mice density dependence was delayed. B. birtlesii prevalence in wood mice was related to bank vole density. The implications of these differences for species coexistence are discussed.