Selection despite low genetic diversity and high gene flow in a rapid island invasion of the bumblebee, Bombus terrestris.

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.

Morphological Variation in Bumblebees (Bombus terrestris) (Hymenoptera: Apidae) After Three Decades of an Island Invasion.

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.

Latitudinal clines in sexual selection, sexual size dimorphism and sex-specific genetic dispersal during a poleward range expansion.

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.

Temporal and spatial variation in sex-specific abundance of the avian vampire fly (Philornis downsi).

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.

Female in-nest attendance predicts the number of ectoparasites in Darwin's finch species.

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.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers.

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 FST 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.

Evidence for rapid downward fecundity selection in an ectoparasite (Philornis downsi) with earlier host mortality in Darwin's finches.

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.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus.

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.

Detection of environmental and morphological adaptation despite high landscape genetic connectivity in a pest grasshopper (Phaulacridium vittatum).

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 FST  = 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, FST 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.

Signatures of local adaptation along environmental gradients in a range-expanding damselfly (Ischnura elegans).

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.

The search for loci under selection: trends, biases and progress.

Detecting genetic variants under selection using FST 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.

Baseline and stress-induced blood properties of male and female Darwin's small ground finch (Geospiza fuliginosa) of the Galapagos Islands.

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.

Morphological Variation Tracks Environmental Gradients in an Agricultural Pest, Phaulacridium vittatum (Orthoptera: Acrididae).

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.

Host-specific associations affect the microbiome of Philornis downsi, an introduced parasite to the Galápagos Islands.

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.

Applying landscape genetics to the microbial world.

Landscape genetics, which explicitly quantifies landscape effects on gene flow and adaptation, has largely focused on macroorganisms, with little attention given to microorganisms. This is despite overwhelming evidence that microorganisms exhibit spatial genetic structuring in relation to environmental variables. The increasing accessibility of genomic data has opened up the opportunity for landscape genetics to embrace the world of microorganisms, which may be thought of as 'the invisible regulators' of the macroecological world. Recent developments in bioinformatics and increased data accessibility have accelerated our ability to identify microbial taxa and characterize their genetic diversity. However, the influence of the landscape matrix and dynamic environmental factors on microorganism genetic dispersal and adaptation has been little explored. Also, because many microorganisms coinhabit or codisperse with macroorganisms, landscape genomic approaches may improve insights into how micro- and macroorganisms reciprocally interact to create spatial genetic structure. Conducting landscape genetic analyses on microorganisms requires that we accommodate shifts in spatial and temporal scales, presenting new conceptual and methodological challenges not yet explored in 'macro'-landscape genetics. We argue that there is much value to be gained for microbial ecologists from embracing landscape genetic approaches. We provide a case for integrating landscape genetic methods into microecological studies and discuss specific considerations associated with the novel challenges this brings. We anticipate that microorganism landscape genetic studies will provide new insights into both micro- and macroecological processes and expand our knowledge of species' distributions, adaptive mechanisms and species' interactions in changing environments.

Gene expression under thermal stress varies across a geographical range expansion front.

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.

Dealing with uncertainty in landscape genetic resistance models: a case of three co-occurring marsupials.

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.

Species collapse via hybridization in Darwin's tree finches.

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.

Using multilevel models to identify drivers of landscape-genetic structure among management areas.

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.

Genomic divergence and a lack of recent introgression between commercial and wild bumblebees (Bombus terrestris).

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.