Search

Invasive species are predicted to adjust their morphological, physiological and life-history traits to adapt to their non-native environments. Although a loss of genetic variation during invasion may restrict local adaptation, introduced species often thrive in novel environments. Despite being founded by just a few individuals, Bombus terrestris (Hymenoptera: Apidae) has in less than 30 years successfully spread across the island of Tasmania (Australia), becoming abundant and competitive with native pollinators. We use RADseq to investigate what neutral and adaptive genetic processes associated with environmental and morphological variation allow B. terrestris to thrive as an invasive species in Tasmania. Given the widespread abundance of B. terrestris, we expected little genetic structure across Tasmania and weak signatures of environmental and morphological selection. We found high gene flow with low genetic diversity, although with significant isolation-by-distance and spatial variation in effective migration rates. Restricted migration was evident across the mid-central region of Tasmania, corresponding to higher elevations, pastural land, low wind speeds and low precipitation seasonality. Tajima's D indicated a recent population expansion extending from the south to the north of the island. Selection signatures were found for loci in relation to precipitation, wind speed and wing loading. Candidate loci were annotated to genes with functions related to cuticle water retention and insect flight muscle stability. Understanding how a genetically impoverished invasive bumblebee has rapidly adapted to a novel island environment provides further understanding about the evolutionary processes that determine successful insect invasions, and the potential for invasive hymenopteran pollinators to spread globally., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Introduced social insects can be highly invasive outside of their native range. Around the world, the introduction and establishment of the eusocial bumblebee Bombus terrestris (L. 1758) (Hymenoptera: Apidae) has negatively impacted native pollinators and ecosystems. Understanding how morphological variation is linked to environmental variation across invasive ranges can indicate how rapidly species may be diverging or adapting across novel ranges and may assist with predicting future establishment and spread. Here we investigate whether B. terrestris shows morphological variation related to environmental variation across the island of Tasmania (Australia) where it was introduced three decades ago. We collected 169 workers from 16 sites across Tasmania and related relative abundance and morphology to landscape-wide climate, land use, and vegetation structure. We found weak morphological divergence related to environmental conditions across Tasmania. Body size of B. terrestris was positively associated with the percentage of urban land cover, a relationship largely driven by a single site, possibly reflecting high resource availability in urban areas. Proboscis length showed a significant negative relationship with the percentage of pasture. Wing loading and local abundance were not related to the environmental conditions within sites. Our results reflect the highly adaptable nature of B. terrestris and its ability to thrive in different environments, which may have facilitated the bumblebee's successful invasion across Tasmania., (© The Author(s) 2023. Published by Oxford University Press on behalf of Entomological Society of America.)

Evolutionary change impacts the rate at which insect pests, pollinators, or disease vectors expand or contract their geographic ranges. Although evolutionary changes, and their ecological feedbacks, strongly affect these risks and associated ecological and economic consequences, they are often underappreciated in management efforts. Greater rigor and scope in study design, coupled with innovative technologies and approaches, facilitates our understanding of the causes and consequences of eco-evolutionary dynamics in insect range shifts. Future efforts need to ensure that forecasts allow for demographic and evolutionary change and that management strategies will maximize (or minimize) the adaptive potential of range-shifting insects, with benefits for biodiversity and ecosystem services., Competing Interests: Conflict of interest statement The authors declare no conflict of interest., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Range expansions can be shaped by sex differences in behaviours and other phenotypic traits affecting dispersal and reproduction. Here, we investigate sex differences in morphology, behaviour and genomic population differentiation along a climate-mediated range expansion in the common bluetail damselfly (Ischnura elegans) in northern Europe. We sampled 65 sites along a 583-km gradient spanning the I. elegans range in Sweden and quantified latitudinal gradients in site relative abundance, sex ratio and sex-specific shifts in body size and mating status (a measure of sexual selection). Using single nucleotide polymorphism (SNP) data for 426 individuals from 25 sites, we further investigated sex-specific landscape and climatic effects on neutral genetic connectivity and migration patterns. We found evidence for sex differences associated with the I. elegans range expansion, namely (a) increased male body size with latitude, but no latitudinal effect on female body size, resulting in reduced sexual dimorphism towards the range limit, (b) a steeper decline in male genetic similarity with increasing geographic distance than in females, (c) male-biased genetic migration propensity and (d) a latitudinal cline in migration distance (increasing migratory distances towards the range margin), which was stronger in males. Cooler mean annual temperatures towards the range limit were associated with increased resistance to gene flow in both sexes. Sex ratios became increasingly male biased towards the range limit, and there was evidence for a changed sexual selection regime shifting from favouring larger males in the south to favouring smaller males in the north. Our findings suggest sex-specific spatial phenotype sorting at the range limit, where larger males disperse more under higher landscape resistance associated with cooler climates. The combination of latitudinal gradients in sex-biased dispersal, increasing male body size and (reduced) sexual size dimorphism should have emergent consequences for sexual selection dynamics and the mating system at the expanding range front. Our study illustrates the importance of considering sex differences in the study of range expansions driven by ongoing climate change., (© 2021 British Ecological Society.)

The global movement of bees for agricultural pollination services can affect local pollinator populations via hybridization. When commercial bumblebees are of the same species but of different geographic origin, intraspecific hybridization may result in beneficial integration of new genetic variation, or alternatively may disrupt locally adapted gene complexes. However, neither the existence nor the extent of genomic introgression and evolutionary divergence between wild and commercial bumblebees is fully understood. We obtained whole-genome sequencing data from wild and commercial Bombus terrestris collected from sites in Southern Sweden with and without long-term use of commercially imported B .  terrestris . We search for evidence of introgression, dispersal and genome-wide differentiation in a comparative genomic analysis of wild and commercial bumblebees. Commercial B .  terrestris were found in natural environments near sites where commercial bumblebees were used, as well as drifting wild B .  terrestris in commercial bumblebee colonies. However, we found no evidence for widespread, recent genomic introgression of commercial B .  terrestris into local wild conspecific populations. We found that wild B .  terrestris had significantly higher nucleotide diversity (Nei's pi, π) , while the number of segregating sites (Watterson's theta, θw ) was higher in commercial B .  terrestris . A highly divergent region on chromosome 11 was identified in commercial B .  terrestris and found to be enriched with structural variants. The genes present in this region are involved in flight muscle contraction and structure and pathogen immune response, providing evidence for differing evolutionary processes operating in wild and commercial B .  terrestris . We did not find evidence for recent introgression, suggesting that co-occurring commercial B .  terrestris have not disrupted evolutionary processes in wild B .  terrestris populations., Competing Interests: None declared., (© 2022 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Understanding the range and behaviour of an invasive species is critical to identify key habitat areas to focus control efforts. Patterns of range use in parasites can differ temporally, across life stages and between sexes. The invasive avian vampire fly, Philornis downsi, spends the larval stage of its life within bird nests, feeding on developing nestlings and causing high levels of mortality and deformation. However, little is known of the ecology and behaviour of the non-parasitic adult fly life stage. Here, we document sex-specific temporal and spatial patterns of abundance of adult avian vampire flies during a single Darwin's finch breeding season. We analyse fly trapping data collected across 7 weeks in the highlands (N = 405 flies) and lowlands (N = 12 flies) of Floreana Island (Galápagos). Lowland catches occurred later in the season, which supports the hypothesis that flies may migrate from the food-rich highlands to the food-poor lowlands once host breeding has commenced. Fly abundance was not correlated with host nesting density (oviposition site) but was correlated with distance to the agricultural zone (feeding site). We consistently caught more males closer to the agricultural zone and more females further away from the agricultural zone. These sex differences suggest that males may be defending or lekking at feeding sites in the agricultural zone for mating. This temporal and sex-specific habitat use of the avian vampire fly is relevant for developing targeted control methods and provides insight into the behavioural ecology of this introduced parasite on the Galápagos Archipelago., (© 2021. The Author(s).)

Selection should act on parental care and favour parental investment decisions that optimize the number of offspring produced. Such predictions have been robustly tested in predation risk contexts, but less is known about alternative functions of parental care under conditions of parasitism. The avian vampire fly ( Philornis downsi ) is a myasis-causing ectoparasite accidentally introduced to the Galápagos Islands, and one of the major mortality causes in Darwin's finch nests. With an 11-year dataset spanning 21 years, we examine the relationship between parental care behaviours and number of fly larvae and pupae in Darwin's finch nests. We do so across three host species ( Camarhynchus parvulus , C. pauper , Geospiza fuliginosa ) and one hybrid Camarhynchus group. Nests with longer female brooding duration (minutes per hour spent sitting on hatchlings to provide warmth) had fewer parasites, and this effect depended on male food delivery to chicks. Neither male age nor number of nest provisioning visits were directly associated with number of parasites. While the causal mechanisms remain unknown, we provide the first empirical study showing that female brooding duration is negatively related to the number of ectoparasites in nests. We predict selection for coordinated host male and female behaviour to reduce gaps in nest attendance, especially under conditions of novel and introduced ectoparasites.

In invasive parasites, generalism is considered advantageous during the initial phase of introduction. Thereafter, fitness costs to parasites, such as host-specific mortality, can drive parasites towards specialism to avoid costly hosts. It is important to determine changes in host specificity of invasive populations to understand host-parasite dynamics and their effects on vulnerable host populations. We examined changes in mortality in the introduced avian vampire fly (Philornis downsi) (Diptera: Muscidae), a generalist myasis-causing ectoparasite, between 2004 and 2020 on Floreana Island (Galápagos). Mortality was measured as the proportion of immature larvae found upon host nest termination. Over the time period, the avian vampire fly was most abundant and had low mortality in nests of the critically endangered medium tree finch (Camarhynchus pauper) and had the highest mortality in nests of hybrid tree finches (Camarhynchus spp.). Low larval mortality was also found in small tree (Camarhynchus parvulus) and small ground finch (Geospiza fuliginosa) nests. Selection could favour avian vampire flies that select medium tree finch nests and/or avoid hybrid nests. Overall, the finding of differences in avian vampire fly survival across host species is parsimonious with the idea that the introduced fly may be evolving towards host specialisation., (© 2021. The Author(s).)

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation., (© 2020 John Wiley & Sons Ltd.)

Fecundity selection is a critical component of fitness and a major driver of adaptive evolution. Trade-offs between parasite mortality and host resources are likely to impose a selection pressure on parasite fecundity, but this is little studied in natural systems. The 'fecundity advantage hypothesis' predicts female-biased sexual size dimorphism whereby larger females produce more offspring. Parasitic insects are useful for exploring the interplay between host resource availability and parasite fecundity, because female body size is a reliable proxy for fecundity in insects. Here we explore temporal changes in body size in the myiasis-causing parasite Philornis downsi (Diptera: Muscidae) on the Galápagos Islands under conditions of earlier in-nest host mortality. We aim to investigate the effects of decreasing host resources on parasite body size and fecundity. Across a 12-year period, we observed a mean of c. 17% P. downsi mortality in host nests with 55 ± 6.2% host mortality and a trend of c. 66% higher host mortality throughout the study period. Using specimens from 116 Darwin's finch nests (Passeriformes: Thraupidae) and 114 traps, we found that over time, P. downsi pupae mass decreased by c. 32%, and male (c. 6%) and female adult size (c. 11%) decreased. Notably, females had c. 26% smaller abdomens in later years, and female abdomen size was correlated with number of eggs. Our findings imply natural selection for faster P. downsi pupation and consequently smaller body size and lower parasite fecundity in this newly evolving host-parasite system., (© 2020 The Authors. Journal of Evolutionary Biology published by John Wiley & Sons Ltd on behalf of European Society for Evolutionary Biology.)

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

Conservation strategies based on charismatic flagship species, such as tigers, lions, and elephants, successfully attract funding from individuals and corporate donors. However, critics of this species-focused approach argue it wastes resources and often does not benefit broader biodiversity. If true, then the best way of raising conservation funds excludes the best way of spending it. Here we show that this conundrum can be resolved, and that the flagship species approach does not impede cost-effective conservation. Through a tailored prioritization approach, we identify places containing flagship species while also maximizing global biodiversity representation (based on 19,616 terrestrial and freshwater species). We then compare these results to scenarios that only maximized biodiversity representation, and demonstrate that our flagship-based approach achieves 79-89% of our objective. This provides strong evidence that prudently selected flagships can both raise funds for conservation and help target where these resources are best spent to conserve biodiversity.

Darwin's finch species in the Galapagos Archipelago are an iconic adaptive radiation that offer a natural experiment to test for the various factors that influence gut microbiome composition. The island of Floreana has the longest history of human settlement within the archipelago and offers an opportunity to compare island and habitat effects on Darwin's finch microbiomes. In this study, we compare gut microbiomes in Darwin's finch species on Floreana Island to test for effects of host phylogeny, habitat (lowlands, highlands), and island (Floreana, Santa Cruz). We used 16S rRNA Illumina sequencing of fecal samples to assess the gut microbiome composition of Darwin's finches, complemented by analyses of stable isotope values and foraging data to provide ecological context to the patterns observed. Overall bacterial composition of the gut microbiome demonstrated co-phylogeny with Floreana hosts, recapitulated the effect of habitat and diet, and showed differences across islands. The finch phylogeny uniquely explained more variation in the microbiome than did foraging data. Finally, there were interaction effects for island × habitat, whereby the same Darwin's finch species sampled on two islands differed in microbiome for highland samples (highland finches also had different diets across islands) but not lowland samples (lowland finches across islands had comparable diet). Together, these results corroborate the influence of phylogeny, age, diet, and sampling location on microbiome composition and emphasize the necessity for comprehensive sampling given the multiple factors that influence the gut microbiome in Darwin's finches, and by extension, in animals broadly., Competing Interests: The authors have declared that no competing interests exist.

Darwin's finches are an iconic example of an adaptive radiation with well-characterized evolutionary history, dietary preferences, and biogeography, offering an unparalleled opportunity to disentangle effects of evolutionary history on host microbiome from other factors like diet and habitat. Here, we characterize the gut microbiome in Darwin's finches, comparing nine species that occupy diverse ecological niches on Santa Cruz island. The finch phylogeny showed moderate congruence with the microbiome, which was comprised mostly of the bacterial phyla Firmicutes, Actinobacteria, and Proteobacteria. Diet, as measured with stable isotope values and foraging observations, also correlated with microbiome differentiation. Additionally, each gut microbial community could easily be classified by the habitat of origin independent of host species. Altogether, these findings are consistent with a model of microbiome assembly in which environmental filtering via diet and habitat are primary determinants of the bacterial taxa present with lesser influence from the evolutionary history between finch species.

Widespread species that exhibit both high gene flow and the capacity to occupy heterogeneous environments make excellent models for examining local selection processes along environmental gradients. Here we evaluate the influence of temperature and landscape variables on genetic connectivity and signatures of local adaptation in Phaulacridium vittatum, a widespread agricultural pest grasshopper, endemic to Australia. With sampling across a 900-km latitudinal gradient, we genotyped 185 P. vittatum from 19 sites at 11,408 single nucleotide polymorphisms (SNPs) using ddRAD sequencing. Despite high gene flow across sites (pairwise F ST  = 0.0003-0.08), landscape genetic resistance modelling identified a positive nonlinear effect of mean annual temperature on genetic connectivity. Urban areas and water bodies had a greater influence on genetic distance among sites than pasture, agricultural areas and forest. Together, F ST outlier tests and environmental association analysis (EAA) detected 242 unique SNPs under putative selection, with the highest numbers associated with latitude, mean annual temperature and body size. A combination of landscape genetic connectivity analysis together with EAA identified mean annual temperature as a key driver of both neutral gene flow and environmental selection processes. Gene annotation of putatively adaptive SNPs matched with gene functions for olfaction, metabolic detoxification and ultraviolet light shielding. Our results imply that this widespread agricultural pest has the potential to spread and adapt under shifting temperature regimes and land cover change., (© 2019 John Wiley & Sons Ltd.)

Invertebrate pests often show high morphological variation and wide environmental tolerances. Knowledge of how phenotypic variation is associated with environmental heterogeneity can elucidate the processes underpinning these patterns. Here we examine morphological variation and relative abundance along environmental gradients in a widespread agricultural pest, native to Australia, the wingless grasshopper Phaulacridium vittatum (Sjöstedt). We test for correlations between body size, wing presence, and stripe polymorphism with environmental variables. Using multiple regression and mixed-effects modeling, body size and stripe polymorphism were positively associated with solar radiation, and wing presence was positively associated with foliage projective cover (FPC). There were no associations between body size or morphological traits with relative abundance. However, relative abundance was positively associated with latitude, soil moisture, and wind speed, but was negatively associated with FPC. Therefore, sites with low relative abundance and high forest cover were more likely to contain winged individuals. Overall, our results suggest that environmental and climatic conditions strongly influence the relative abundance and the distribution of morphotypes in P. vittatum, which is likely to affect dispersal and fitness in different landscapes. This knowledge is useful for informing how environmental change might influence the future spread and impact of this agricultural pest.

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range-expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site-associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V-ATPase) and visual processes (long-wavelength-sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate-mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients., (© 2018 John Wiley & Sons Ltd.)

Birds are renowned for exhibiting marked sex-specific differences in activity levels and reproductive investment during the breeding season, potentially impacting circulating blood parameters associated with stress and energetics. Males of many passerines often do not incubate, but they experience direct exposure to intruder threat and exhibit aggressive behaviour during the nesting phase in order to defend territories against competing males and predators. Nesting females often have long bouts of inactivity during incubation, but they must remain vigilant of the risks posed by predators and conspecific intruders approaching the nest. Here, we use 33 free-living male (n = 16) and female (n = 17) Darwin's small ground finches (Geospiza fuliginosa) on Floreana Island (Galapagos Archipelago) to better understand how sex-specific roles during the reproductive period impact baseline and stress-induced levels of plasma corticosterone (CORT), blood glucose and haematocrit. Specifically, we hypothesise that males are characterised by higher baseline values given their direct and relatively frequent exposure to intruder threat, but that a standardised stress event (capture and holding) overrides any sex-specific differences. In contrast with expectations, baseline levels of all blood parameters were similar between sexes (13.4 ± 1.9 ng ml -1 for CORT, 13.7 ± 0.4 mmol l -1 for glucose, 58.3 ± 0.8% for haematocrit). Interestingly, females with higher body condition had lower baseline haematocrit. All blood parameters changed with time since capture (range 1.2-41.3 min) in both sexes, whereby CORT increased linearly, haematocrit decreased linearly, and glucose increased to a peak at ∼20 min post-capture and declined to baseline levels thereafter. Our results do not support the hypothesis that sex-specific roles during the reproductive period translate to differences in blood parameters associated with stress and energetics, but we found some evidence that blood oxygen transport capacity may decline as finches increase in body condition., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Detecting genetic variants under selection using F ST outlier analysis (OA) and environmental association analyses (EAAs) are popular approaches that provide insight into the genetic basis of local adaptation. Despite the frequent use of OA and EAA approaches and their increasing attractiveness for detecting signatures of selection, their application to field-based empirical data have not been synthesized. Here, we review 66 empirical studies that use Single Nucleotide Polymorphisms (SNPs) in OA and EAA. We report trends and biases across biological systems, sequencing methods, approaches, parameters, environmental variables and their influence on detecting signatures of selection. We found striking variability in both the use and reporting of environmental data and statistical parameters. For example, linkage disequilibrium among SNPs and numbers of unique SNP associations identified with EAA were rarely reported. The proportion of putatively adaptive SNPs detected varied widely among studies, and decreased with the number of SNPs analysed. We found that genomic sampling effort had a greater impact than biological sampling effort on the proportion of identified SNPs under selection. OA identified a higher proportion of outliers when more individuals were sampled, but this was not the case for EAA. To facilitate repeatability, interpretation and synthesis of studies detecting selection, we recommend that future studies consistently report geographical coordinates, environmental data, model parameters, linkage disequilibrium, and measures of genetic structure. Identifying standards for how OA and EAA studies are designed and reported will aid future transparency and comparability of SNP-based selection studies and help to progress landscape and evolutionary genomics., (© 2018 John Wiley & Sons Ltd.)

The composition and diversity of bacteria forming the microbiome of parasitic organisms have implications for differential host pathogenicity and host-parasite co-evolutionary interactions. The microbiome of pathogens can therefore have consequences that are relevant for managing disease prevalence and impact on affected hosts. Here, we investigate the microbiome of an invasive parasitic fly Philornis downsi, recently introduced to the Galápagos Islands, where it poses extinction threat to Darwin's finches and other land birds. Larvae infest nests of Darwin's finches and consume blood and tissue of developing nestlings, and have severe mortality impacts. Using 16s rRNA sequencing data, we characterize the bacterial microbiota associated with P. downsi adults and larvae sourced from four finch host species, inhabiting two islands and representing two ecologically distinct groups. We show that larval and adult microbiomes are dominated by the phyla Proteobacteria and Firmicutes, which significantly differ between life stages in their distributions. Additionally, bacterial community structure significantly differed between larvae retrieved from strictly insectivorous warbler finches (Certhidea olivacea) and those parasitizing hosts with broader dietary preferences (ground and tree finches, Geospiza and Camarhynchus spp., respectively). Finally, we found no spatial effects on the larval microbiome, as larvae feeding on the same host (ground finches) harboured similar microbiomes across islands. Our results suggest that the microbiome of P. downsi changes during its development, according to dietary composition or nutritional needs, and is significantly affected by host-related factors during the larval stage. Unravelling the ecological significance of bacteria for this parasite will contribute to the development of novel, effective control strategies., (© 2017 John Wiley & Sons Ltd.)

Many ectothermic species are currently expanding their distributions polewards due to anthropogenic global warming. Molecular genetic mechanisms facilitating range expansion under these conditions are largely unknown, but understanding these could help mitigate expanding pests and disease vectors, or help explain why some species fail to track changing climates. Here, using RNA-seq data, we examine genomewide changes in gene expression under heat and cold stress in the range-expanding damselfly Ischnura elegans in northern Europe. We find that both the number of genes involved and levels of gene expression under heat stress have become attenuated during the expansion, consistent with a previously reported release from selection on heat tolerances as species move polewards. Genes upregulated under cold stress differed between core and edge populations, corroborating previously reported rapid adaptation to cooler climates at the expansion front. Expression of sixty-nine genes exhibited a region x treatment effect; these were primarily upregulated in response to heat stress in core populations but in response to cold stress at the range edge, suggesting that some cellular responses originally adapted to heat stress may switch to cold-stress functionality upon encountering novel thermal selection regimes during range expansion. Transcriptional responses to thermal stress involving heat-shock and neural function genes were largely geographically conserved, while retrotransposon, regulatory, muscle function and defence gene expression patterns were more variable. Flexible mechanisms of cold-stress response and the ability of some genes to shift their function between heat and cold stress might be key mechanisms facilitating rapid poleward expansion in insects., (© 2016 John Wiley & Sons Ltd.)

Landscape genetics lacks explicit methods for dealing with the uncertainty in landscape resistance estimation, which is particularly problematic when sample sizes of individuals are small. Unless uncertainty can be quantified, valuable but small data sets may be rendered unusable for conservation purposes. We offer a method to quantify uncertainty in landscape resistance estimates using multimodel inference as an improvement over single model-based inference. We illustrate the approach empirically using co-occurring, woodland-preferring Australian marsupials within a common study area: two arboreal gliders (Petaurus breviceps, and Petaurus norfolcensis) and one ground-dwelling antechinus (Antechinus flavipes). First, we use maximum-likelihood and a bootstrap procedure to identify the best-supported isolation-by-resistance model out of 56 models defined by linear and non-linear resistance functions. We then quantify uncertainty in resistance estimates by examining parameter selection probabilities from the bootstrapped data. The selection probabilities provide estimates of uncertainty in the parameters that drive the relationships between landscape features and resistance. We then validate our method for quantifying uncertainty using simulated genetic and landscape data showing that for most parameter combinations it provides sensible estimates of uncertainty. We conclude that small data sets can be informative in landscape genetic analyses provided uncertainty can be explicitly quantified. Being explicit about uncertainty in landscape genetic models will make results more interpretable and useful for conservation decision-making, where dealing with uncertainty is critical., (© 2015 John Wiley & Sons Ltd.)

Species hybridization can lead to fitness costs, species collapse, and novel evolutionary trajectories in changing environments. Hybridization is predicted to be more common when environmental conditions change rapidly. Here, we test patterns of hybridization in three sympatric tree finch species (small tree finch Camarhynchus parvulus, medium tree finch Camarhynchus pauper, and large tree finch: Camarhynchus psittacula) that are currently recognized on Floreana Island, Galápagos Archipelago. Genetic analysis of microsatellite data from contemporary samples showed two genetic populations and one hybrid cluster in both 2005 and 2010; hybrid individuals were derived from genetic population 1 (small morph) and genetic population 2 (large morph). Females of the large and rare species were more likely to pair with males of the small common species. Finch populations differed in morphology in 1852-1906 compared with 2005/2010. An unsupervised clustering method showed (a) support for three morphological clusters in the historical tree finch sample (1852-1906), which is consistent with current species recognition; (b) support for two or three morphological clusters in 2005 with some (19%) hybridization; and (c) support for just two morphological clusters in 2010 with frequent (41%) hybridization. We discuss these findings in relation to species demarcations of Camarhynchus tree finches on Floreana Island.

Landscape genetics offers a powerful approach to understanding species' dispersal patterns. However, a central obstacle is to account for ecological processes operating at multiple spatial scales, while keeping research outcomes applicable to conservation management. We address this challenge by applying a novel multilevel regression approach to model landscape drivers of genetic structure at both the resolution of individuals and at a spatial resolution relevant to management (i.e. local government management areas: LGAs) for the koala (Phascolartos cinereus) in Australia. Our approach allows for the simultaneous incorporation of drivers of landscape-genetic relationships operating at multiple spatial resolutions. Using microsatellite data for 1106 koalas, we show that, at the individual resolution, foliage projective cover (FPC) facilitates high gene flow (i.e. low resistance) until it falls below approximately 30%. Out of six additional land-cover variables, only highways and freeways further explained genetic distance after accounting for the effect of FPC. At the LGA resolution, there was significant variation in isolation-by-resistance (IBR) relationships in terms of their slopes and intercepts. This was predominantly explained by the average resistance distance among LGAs, with a weaker effect of historical forest cover. Rates of recent landscape change did not further explain variation in IBR relationships among LGAs. By using a novel multilevel model, we disentangle the effect of landscape resistance on gene flow at the fine resolution (i.e. among individuals) from effects occurring at coarser resolutions (i.e. among LGAs). This has important implications for our ability to identify appropriate scale-dependent management actions., (© 2013 John Wiley & Sons Ltd.)

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a 'flat' landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

The Coastal Giant Salamander (Dicamptodon tenebrosus) is classified as threatened at the northern periphery of its range in British Columbia (BC), Canada, primarily due to forestry practices and habitat fragmentation. Characterising dispersal behaviour and population connectivity is therefore a priority for this region, while genetic differentiation in core versus peripheral locations remains unstudied in this wide-ranging species. We present seven new polymorphic microsatellite markers for use in population genetic analyses of D. tenebrosus. We examine locus characteristics and genetic variation in 12 streams at the species' northern range limit in BC, and within two regions representing sub-peripheral (North Cascades) and core localities (South Cascades) in Washington State, United States. In BC, the number of alleles per locus ranged from 2-5 and observed heterozygosity ranged from 0.044-0.825. Genetic differentiation was highest between BC and the South Cascades, and intermediate between BC and the North Cascades. Across loci, mean allelic richness was similar across regions, while private allelic richness was highest in the core locality (corrected for sample size). These new microsatellite loci will be a valuable addition to existing markers for detailed landscape and population genetic analyses of D. tenebrosus across its range.

Background: Understanding the dispersal and genetic structure of invasive insects across islands is important for designing management plans that are appropriate at spatial and temporal scales. For invasive parasites, population dynamics are largely determined by the distribution and density of their host species. The introduced parasitic fly, Philornis downsi, parasitises nestlings of endemic birds on all major islands of the Galápagos archipelago. The fly's high mortality and fitness impacts are of conservation concern for vulnerable and declining species of Darwin's finches. Using microsatellite data in Bayesian clustering and landscape genetic analyses, we examine gene flow and dispersal in P. downsi between three islands and across habitats (highlands, lowlands) and examine for the presence of population bottlenecks. We also examine variation at the mitochondrial gene CO1 across islands to establish if cryptic species were present., Results: Both the mitochondrial and microsatellite data were consistent with there being a single species across islands. We found low genetic differentiation between islands and strong evidence for inter-island gene flow, or shared recent ancestry among individuals. Landscape genetic analysis identified two genetic clusters: one encompassing Santa Cruz and Isabela, and one on Floreana Island. There was no evidence of genetic differentiation between habitats and molecular variance was mainly attributable to within individuals. The combined P. downsi population was found to have undergone a population bottleneck., Conclusion: Philornis downsi populations have high connectivity within and between islands, with low levels of genetic differentiation between Floreana and the other two islands examined. The genetic bottleneck found across islands suggests there was a small founding population or few introduction events of P. downsi. The high dispersal capacity and wide habitat use of P. downsi highlights the significant threat that this parasite poses to the Galápagos avifauna. Our findings are relevant for assessing the viability of methods to control P. downsi on Galápagos, such as the sterile insect technique.

An enrichment technique was used to isolate 11 di-, tri-, and tetra microsatellites for the parasitic fly Philornis downsi (Diptera: Muscidae). These loci were polymerase chain reaction amplified in singleplexes or two-plexes for P. downsi. The loci showed low to moderate polymorphism, exhibited between three and four alleles, and observed heterozygosity ranged from 0.05 to 0.86. These new markers will be useful for population-level and paternity analyses and will provide valuable information about the ecology of this high-impact parasite of vulnerable bird species., (© 2007 The Authors.)