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.

The importance of eco-evolutionary dynamics for predicting and managing insect range shifts.

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.

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.

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.

Avian vampire fly (Philornis downsi) mortality differs across Darwin's finch host species.

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.

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.

Conservation prioritization can resolve the flagship species conundrum.

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.

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.

An inter-island comparison of Darwin's finches reveals the impact of habitat, host phylogeny, and island on the gut microbiome.

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.

Host phylogeny, diet, and habitat differentiate the gut microbiomes of Darwin's finches on Santa Cruz Island.

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.

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.

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.

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.

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.