A new statistical framework for the quantification of covariate associations with species distributions

Identifying processes that shape species geographical ranges is a prerequisite for understanding environmental change. Currently, species distribution modelling methods do not offer credible statistical tests of the relative influence of climate factors and typically ignore other processes (e.g. biotic interactions and dispersal limitation). We use a hierarchical model fitted with Markov Chain Monte Carlo to combine ecologically plausible niche structures using regression splines to describe unimodal but potentially skewed response terms. We apply spatially explicit error terms that account for (and may help identify) missing variables. Using three example distributions of European bird species, we map model results to show sensitivity to change in each covariate. We show that the overall strength of climatic association differs between species and that each species has considerable spatial variation in both the strength of the climatic association and the sensitivity to climate change. Our methods are widely applicable to many species distribution modelling problems and enable accurate assessment of the statistical importance of biotic and abiotic influences on distributions.

Bird diversity and environmental gradients in Britain:a test of the species-energy hypothesis

1. We tested the species diversity-energy hypothesis using the British bird fauna. This predicts that temperature patterns should match diversity patterns. We also tested the hypothesis that the mechanism operates directly through effects of temperature on thermoregulatory loads; this further predicts that seasonal changes in temperature cause matching changes in patterns of diversity, and that species' body mass is influential.

2. We defined four assemblages using migration status (residents or visitors) and season (summer or winter distribution). Records of species' presence/absence in a total of 2362, 10 x 10-km, quadrats covering most of Britain were used, together with a wide selection of habitat, topographic and seasonal climatic data.

3. We fitted a logistic regression model to each species' distribution using the environmental data. We then combined these individual species models mathematically to form a diversity model. Analysis of this composite model revealed that summer temperature was the factor most strongly associated with diversity.

4. Although the species-energy hypothesis was supported, the direct mechanism, predicting an important role for body mass and matching seasonal patterns of change between diversity and temperature, was not supported.

5. However, summer temperature is the best overall explanation for bird diversity patterns in Britain. It is a better predictor of winter diversity than winter temperature. Winter diversity is predicted more precisely from environmental factors than summer diversity.

6. Climate change is likely to influence the diversity of different areas to different extents; for resident species, low diversity areas may respond more strongly as climate change progresses. For winter visitors, higher diversity areas may respond more strongly, while summer visitors are approximately neutral.

Phylogeographic insights into cryptic glacial refugia

The glacial episodes of the Quaternary (2.6 million years ago–present) were a major factor in shaping the present-day distributions of extant flora and fauna, with expansions and contractions of the ice sheets rendering large areas uninhabitable for most species. Fossil records suggest that many species survived glacial maxima by retreating to refugia, usually at lower latitudes. Recently, phylogeographic studies have given support to the existence of previously unknown, or cryptic, refugia. Here we summarise many of these insights into the glacial histories of species in cryptic refugia gained through phylogeographic approaches. Understanding such refugia might be important as the Earth heads into another period of climate change, in terms of predicting the effects on species distribution and survival.

Range-edge effects promote clonal growth in peripheral populations of the one-sided wintergreen Orthilia secunda

Existing in suboptimal conditions is a frequent occurrence for species inhabiting the cusp of their ecological range. In range-edge populations of plants, the scarcity of suitable habitat may be reflected in small population sizes which may result in increased self-pollination and/or inbreeding and an increase in the incidence of clonal reproduction. These factors may result in a decrease in levels of genetic diversity and a loss of potential adaptive variation that may compromise species' ability to cope with changes in their environment, an issue that is particularly relevant today with the current concern surrounding global climate change and its effect on species' distributional ranges. In the present study, we have compared the levels of clonal reproduction in the one-sided wintergreen Orthilia secunda (L.) House in (1) populations from its main continuous distribution range, (2) populations occurring on the limits of the continuous range, and (3) peripheral populations outwith the species' continuous distribution range. Range-edge populations in Scotland and Sweden displayed significantly lower genotypic richness and diversity than those from the main area of the species' distribution in these countries. Populations from Ireland, which occur in the temperate zone rather than the boreal conditions that are the preferred habitat for the species, and which represent relict populations left over from cooler periods in the Earth's history, displayed no within-population genetic diversity, suggesting a complete lack of sexual reproduction. Furthermore, the genetic distinctiveness of the Irish populations, which contained alleles not found in either the Scottish or the Swedish populations, highlights the value of 'trailing edge' populations and supports the concept of 'parochial conservation', namely the conservation of species that are locally rare but globally common.

Genetic differentiation of Euterpe edulis Mart. populations estimated by AFLP analysis

Heart-of-palm (Euterpe edulis Mart.) is a wild palm with a wide distribution throughout the Atlantic Rainforest. Populations of E. edulis represent important renewable natural resources but are currently under threat from predatory exploitation. Furthermore, because the species is indigenous to the Atlantic Rainforest, which is located in the most economically developed and populated region of Brazil, social and economic pressures have devastated heart-of-palm forests. In order to estimate the partitioning of genetic variation of endangered E. edulis populations, 429 AFLP markers were used to analyse 150 plants representing 11 populations of the species distribution range. Analysis of the genetic structure of populations carried out using analysis of molecular variance (AMOVA) revealed moderate genetic variation within populations (57.4%). Genetic differentiation between populations (F-ST = 0.426) was positively correlated with geographical distance. These results could be explained by the historical fragmentation of the Atlantic coastal region, together with the life cycle and mating system The data obtained in this work should have important implications for conservation and future breeding programmes of E. edulis.

Phylogeographic analysis of North American populations of the parasitic herbaceous plant Monotropa hypopitys reveals a complex history of range expansion from multiple late glacial refugia

Aim We carried out a phylogeographic study across the range of the herbaceous plant species Monotropa hypopitys L. in North America to determine whether its current disjunct distribution is due to recolonization from separate eastern and western refugia after the Last Glacial Maximum (LGM). Location North America: Pacific Northwest and north-eastern USA/south-eastern Canada. Methods Palaeodistribution modelling was carried out to determine suitable climatic regions for M. hypopitys at the LGM. We analysed between 155 and 176 individuals from 39 locations spanning the species' entire range in North America. Sequence data were obtained for the chloroplast rps2 gene (n=168) and for the nuclear ITS region (n=158). Individuals were also genotyped for eight microsatellite loci (n=176). Interpolation of diversity values was used to visualize the range-wide distribution of genetic diversity for each of the three marker classes. Minimum spanning networks were constructed showing the relationships between the rps2 and ITS haplotypes, and the geographical distributions of these haplotypes were plotted. The numbers of genetic clusters based on the microsatellite data were estimated using Bayesian clustering approaches. Results The palaeodistribution modelling indicated suitable climate envelopes for M. hypopitys at the LGM in both the Pacific Northwest and south-eastern USA. High levels of genetic diversity and endemic haplotypes were found in Oregon, the Alexander Archipelago, Wisconsin, and in the south-eastern part of the species' distribution range. Main conclusions Our results suggest a complex recolonization history for M. hypopitys in North America, involving persistence in separate eastern and western refugia. A generally high degree of congruence between the different marker classes analysed indicated the presence of multiple refugia, with at least two refugia in each area. In the west, putative refugia were identified in Oregon and the Alexander Archipelago, whereas eastern refugia may have been located in the southern part of the species' current distribution, as well as in the 'Driftless Area'. These findings are in contrast to a previous study on the related species Orthilia secunda, which has a similar disjunct distribution to M. hypopitys, but which appears to have recolonized solely from western refugia. © 2011 Blackwell Publishing Ltd.

Trends in candidemia and antifungal susceptibility in a university hospital in Northern Ireland 2001-2006

Objectives: To describe the species distribution and antifungal susceptibility trends for documented episodes of candidemia at the Royal Hospitals, Belfast, 2001-2006. Methods: Laboratory-based retrospective observational study of all episodes of candidemia. Results: There were 151 episodes of candidemia. The species recovered were: 96 C. albicans; 26 C. glabrata; 18 C. parapsilosis; five C. tropicalis; four C. guilliermondii; one C. famata and one C. dubliniensis. We separated the data into two periods 2001-2003 and 2004-2006; contrary to the findings of other investigators, there was a notable trends toward increasing frequency of C. albicans and decreasing frequency of non-albicans species over time. Although the proportion of C. albicans, C. parapsilosis and C. tropicalis isolates susceptible to fluconazole was unchanged over time, a trend of decreased susceptibility of C. glabrata to fluconazole was noted over the six-year period. Overall, 73% and 7.7% of C. glabrata isolates had susceptible-dose-dependent and resistant phenotypes, respectively. The percentage of C. glabrata isolates susceptible to fluconazole (MIC

Refining Risk Estimates Using Models (Chapter 4)

In 2004 nineteen scientists from fourteen institutions in seven countries
collaborated in the landmark study described in chapter 2 (Thomas et al., 2004a). This chapter provides an overview of results of studies published subsequently and assesses how much, and why, new results differ from those of Thomas et al.
Some species distribution modeling (SDM) studies are directly comparable to the Thomas et al. estimates. Others using somewhat different methods nonetheless illuminate whether the original estimates were of the right order of magnitude. Climate similarity models (Williams et al., 2007; Williams and Jackson, 2007), biome, and vegetation dynamic models (Perry and Enright, 2006) have also been
applied in the context of climate change, providing interesting opportunities
for comparison and cross-validation with results from SDMs.
This chapter concludes with an assessment of whether the range of extinction risk estimates presented in 2004 can be narrowed, and whether the mean estimate should be revised upward or downward. To set the stage for these analyses, the chapter begins with brief reviews of advances in climate modeling and species modeling since 2004.

Spatial autocorrelation in the response of soft-bottom marine benthos to gas extraction activities: The case of amphipods in the Ionian Sea

The spatial distributions of marine fauna and of pollution are both highly structured, and thus the resulting high levels of autocorrelation may invalidate conclusions based on classical statistical approaches. Here we analyse the close correlation observed between proxies for the disturbance associated with gas extraction activities and amphipod distribution patterns around four hydrocarbon platforms. We quantified the amount of variation independently accounted for by natural environmental variables, proxies for the disturbance caused by platforms, and spatial autocorrelation. This allowed us to demonstrate how each of these three factors significantly affects the community structure of amphipods. Sophisticated statistical techniques are required when taking into account spatial autocorrelation: nevertheless our data demonstrate that this approach not only enables the formulation of robust statistical inferences but also provides a much deeper understanding of the subtle interactions between human disturbance and natural factors affecting the structure of marine invertebrates communities. (C) 2012 Elsevier Ltd. All rights reserved.

Protected area networks and savannah bird biodiversity in the face of climate change and land degradation

The extent to which climate change might diminish the efficacy of protected areas is one of the most pressing conservation questions. Many projections suggest that climate-driven species distribution shifts will leave protected areas impoverished and species inadequately protected while other evidence suggests that intact ecosystems within protected areas will be resilient to change. Here, we tackle this problem empirically. We show how recent changes in distribution of 139 Tanzanian savannah bird species are linked to climate change, protected area status and land degradation. We provide the first evidence of climate-driven range shifts for an African bird community. Our results suggest that the continued maintenance of existing protected areas is an appropriate conservation response to the challenge of climate and environmental change.

Identifying optimal feeding habitat and proposed Marine Protected Areas (pMPAs) for the black-legged kittiwake (Rissa tridactyla) suggests a need for complementary management approaches

Marine Protected Areas (MPAs) are an important conservation tool. For marine predators, recent research has focused on the use of Species Distribution Models (SDMs) to identify proposed sites. We used a maximum entropy modelling approach based on static and dynamic oceanographic parameters to determine optimal feeding habitat for black-legged kittiwakes (Rissa tridactyla) at two colonies during two consecutive breeding seasons (2009 and 2010). A combination of Geographic Positioning System (GPS) loggers and Time-Depth Recorders (TDRs) attributed feeding activity to specific locations. Feeding areas were <30 km from the colony, <40 km from land, in productive waters, 25–175m deep. The predicted extent of optimal habitat declined at both colonies between 2009 and 2010 coincident with declines in reproductive success. Whilst the area of predicted optimal habitat changed, its location was spatially stable between years. There was a close match between observed feeding locations and habitat predicted as optimal at one colony (Lambay Island, Republic of Ireland), but a notable mismatch at the other (Rathlin Island, Northern Ireland). Designation of an MPA at Rathlin may, therefore, be less effective than a similar designation at Lambay perhaps due to the inherent variability in currents and sea state in the North Channel compared to the comparatively stable conditions in the central Irish Sea. Current strategies for designating MPAs do not accommodate likely future redistribution of resources due to climate change. We advocate the development of new approaches including dynamic MPAs that track changes in optimal habitat and non-colony specific ecosystem management.

Restrospective genetic monitoring of the threatened Yellow marsh saxifrage (Saxifraga hirculus) reveals genetic erosion but provides valuable information for conservation strategies

Aim: Retrospective genetic monitoring, comparing genetic diversity of extant populations with historical samples, can provide valuable and often unique insights into evolutionary processes informing conservation strategies. The Yellow marsh saxifrage (Saxifraga hirculus) is listed as ‘critically endangered’ in Ireland with only two extant populations. We quantified genetic changes over time and identified genotypes in extant populations that could be used as founders for reintroductions to sites where the species is extinct.

Location: Ireland.

Methods: Samples were obtained from both locations where the species is currently found, including the most threatened site at the Garron Plateau, Co. Antrim, which held only 13 individuals during 2011. Herbarium samples covering the period from 1886 to 1957 were obtained including plants from the same area as the most threatened population, as well as three extinct populations. In total, 422 individuals (319 present-day and 103 historical) were genotyped at six microsatellite loci. Species distribution modelling was used to identify areas of potentially suitable habitat for reintroductions.

Results: Level of phenotypic diversity within the most threatened population was significantly lower in the present-day compared with historical samples but levels of observed heterozygosity and number of alleles, whilst reduced, did not differ significantly. However, Bayesian clustering analysis suggested gradual lineage replacement over time. All three measures of genetic diversity were generally lower at the most threatened population compared with the more substantial extant populations in Co. Mayo. Species distribution modelling suggested that habitat at one site where the species is extinct may be suitable for reintroduction.

Main conclusions: The dominant genetic lineage in the most threatened population is rare elsewhere; thus, care needs to be taken when formulating any potential reintroduction programme. Our findings highlight both the need for genetic monitoring of threatened populations, but also for its swift implementation before levels of diversity become critically low.

Cryptic introgression into the Kidney saxifrage (Saxifraga hirsuta) from its more abundant sympatric congener Saxifraga spathularis, and the potential risk of genetic assimilation

Background and Aims: Although hybridization can play a positive role in plant evolution, it has been shown that excessive unidirectional hybridization can result in replacement of a species’ gene pool, and even the extinction ofrare species via genetic assimilation. This study examines levels of introgression between the common Saxifraga spathularis and its rarer congener S. hirsuta, which have been observed to hybridize in the wild where they occursympatrically. 
Methods: Seven species-specific single nucleotide polymorphisms (SNPs) were analysed in 1025 plants representing both species and their hybrid, S. polita, from 29 sites across their ranges in Ireland. In addition, species distributionmodelling was carried out to determine whether the relative abundance of the two parental species is likely to change under future climate scenarios. 
Key Results: Saxifraga spathularis individuals tended to be genetically pure, exhibiting little or no introgression from S. hirsuta, but significant levels of introgression of S. spathularis alleles into S. hirsuta were observed, indicatingthat populations exhibiting S. hirsuta morphology are more like a hybrid swarm, consisting of backcrosses and F2s. Populations of the hybrid, S. polita, were generally comprised of F1s or F2s, with some evidence of backcrossing. Species distribution modelling under projected future climate scenarios indicated an increase in suitable habitats for both parental species.
Conclusions: Levels of introgression observed in this study in both S. spathularis and S. hirsuta would appear to be correlated with the relative abundance of the species. Significant introgression of S. spathularis alleles was detectedin the majority of the S. hirsuta populations analysed and, consequently, ongoing introgression would appear to represent a threat to the genetic integrity of S. hirsuta, particularly in areas where the species exists sympatricallywith its congener and where it is greatly outnumbered.

Burkholderia cepacia complex epidemiology in persons with cystic fibrosis from Australia and New Zealand

The Burkholderia cepacia complex (Bcc) is a group of significant opportunistic respiratory pathogens which affect people with cystic fibrosis. In this study, we sought to ascertain the epidemiology and geographic species distribution of 116 Bcc isolates collected from people with CF in Australia and New Zealand. We performed a combination of recA-based PCR, amplified rDNA restriction analysis (ARDRA), pulsed-field gel electrophoresis and repetitive extragenic palindromic PCR on each isolate. Each Burkholderia cenocepacia isolate was also screened by PCR for the presence of the B. cepacia epidemic strain marker. One hundred and fourteen isolates were assigned to a species using recA-based PCR and ARDRA. B. cenocepacia, B. multivorans and B. cepacia accounted for 45.7%, 29.3% and 11.2% of the isolates, respectively. Strain analysis of B. cenocepacia revealed that 85.3% of the isolates were unrelated. One related B. cenocepacia strain was identified amongst 15 people. Whilst full details of person-to-person contact was not available, all patients attended CF centres in Queensland (Qld) and New South Wales (NSW). Although person-to-person transmission of B. cenocepacia strains has occurred in Australia, the majority of CF-related Bcc infections in Australia and New Zealand are most likely acquired from the environment.

MICROBIAL COMMUNITY COMPOSITION IN THE LUNGS OF PATIENTS WITH CYSTIC FIBROSIS AND NON-CF BRONCHIECTASIS

Introduction and Aims: Persistent bacterial infection is a major cause of morbidity and mortality in patients with both Cystic Fibrosis (CF) and non-CF Bronchiectasis (non-CFBX). Numerous studies have shown that CF and non-CFBX airways are colonised by a complex microbiota. However, many bacteria are difficult, if not impossible, to culture by conventional laboratory techniques. Therefore, molecular detection techniques offer a more comprehensive view of bacterial diversity within clinical specimens. The objective of this study was to characterise and compare bacterial diversity and relative abundance in patients with CF and non-CFBX during exacerbation and when clinically stable.

Methods: Sputum samples were collected from CF (n=50 samples) and non-CFBX (n=52 samples) patients at the start and end of treatment for an infective exacerbation and when clinically stable. Pyrosequencing was used to assess the microbial diversity and relative genera (or the closest possibly taxonomic order) abundance within the samples. Each sequence read was defined based on 3% difference.

Results: High-throughput pyrosequencing allowed a sensitive and detailed examination of microbial community composition. Rich microbial communities were apparent within both CF (171 species-level phylotypes per genus) and non-CFBX airways (144 species-level phylotypes per genus). Relative species distribution within those two environments was considerably different; however, relatively few genera formed a core of microorganisms, representing approximately 90% of all sequences, which dominated both environments. Relative abundance based on observed operational taxonomic units demonstrated that the most abundant bacteria in CF were Pseudomonas (28%), Burkholderia (22%), Streptococcus (13%), family Pseudomonadaceae (8%) and Prevotella (6%). In contrast, the most commonly detected operational taxonomic units in non-CFBX were Haemophilus (22%), Streptococcus (14%), other (unassigned taxa) (11%), Pseudomonas (10%), Veillonella (7%) and Prevotella (6%).

Conclusions: These results suggest that distinctive microbial communities are associated with infection and/or colonisation in patients with both CF and non-CFBX. Although relatively high species richness was observed within the two environments, each was dominated by different core taxa. This suggests that differences in the lung environment of these two diseases may affect adaptability of the relevant bacterial taxa.