<title>Background</title> <p>Biological invasions threaten biodiversity and ecosystem function, with pronounced negative effects on islands in particular [##REF##12374870##1##, ####UREF##0##2##, ##UREF##1##3####1##3##]. Genetic studies of invasive species can identify the adaptive potential of invaders to deal with new environmental conditions [##UREF##2##4##] or help to predict evolutionary responses to management practices (e.g. pesticides, biological control agents) [##UREF##3##5##]. Population bottlenecks affect many invasive species because they frequently experience founding effects that reduce genetic variability, but paradoxically, invasive species still manage to successfully establish and adapt to new environments [##REF##16858586##6##]. However, the effects of bottlenecks may be countered by the occurrence of multiple introductions, high reproductive rates, and subsequent migration between locally bottlenecked populations that are genetically differentiated [##REF##15602569##7##].</p> <p>For invasive arthropod parasites, these factors are inextricably linked with the distribution, genetics, and behaviour of host species [##REF##15012332##8##, ####REF##12683525##9##, ##REF##15969711##10####15969711##10##]. The recent integration of molecular ecology with parasitology has provided a path for answering a number of questions concerning the genetic structure of parasite populations, which can uncover a wealth of information regarding ecological and evolutionary processes for invasive parasites [##REF##15969711##10##]. Highly variable multilocus genotypes are particularly suited to analyses of non-equilibrium or bottlenecked populations because they provide adequate variation for assessing recent gene flow and identifying migrants [##REF##10234242##11##].</p> <p>The introduced fly, <italic>Philornis downsi</italic>, is an avian ectoparasite that is considered to be a serious threat to the persistence of endemic finch populations on the Galápagos Islands [##UREF##4##12##, ####UREF##5##13##, ##UREF##6##14####6##14##]. Recently, <italic>P. downsi </italic>was given the highest risk ranking affecting endemic fauna in the Galápagos archipelago [##UREF##1##3##]. Other pathogens affecting Galápagos birds such as avian pox virus [##UREF##7##15##] and intestinal protozoans [##UREF##8##16##] are of less concern, but may also cause high fitness impacts under certain conditions. The fly was first formally identified from Darwin finch nests in 1997 and has since been found on 11 of 13 major islands in nests of 14 endemic species [##UREF##4##12##,##UREF##5##13##]. However, <italic>P. downsi </italic>colonised the islands at least 40 years ago, as the fly was identified recently from collections made in 1964 [##UREF##5##13##]. The blood-feeding larvae of <italic>P. downsi </italic>are associated with 62–100% nestling mortality in Darwin's finches [##UREF##4##12##,##UREF##6##14##,##REF##16899139##17##], as well as physiological costs [##UREF##9##18##] and reduced growth rates in nestlings [##UREF##6##14##]. Little is known about the ecology and biology of <italic>Philornis </italic>flies and the dispersal behaviour and population genetics of the genus <italic>Philornis </italic>or of any other myiasis-causing parasite of birds [reviewed in [##UREF##5##13##]].</p> <p>One potential control method to eradicate <italic>P. downsi </italic>is the sterile insect technique (SIT), which is renowned for its effectiveness at eradicating or suppressing fruit fly and screw-worm fly populations across the globe [##UREF##10##19##,##UREF##11##20##]. SIT involves the large-scale release of laboratory-reared sterile male (and/or female) flies that eventually suppress fly populations by reducing population fecundity [reviewed in [##UREF##11##20##]]. SIT requires a thorough understanding of the reproductive ecology and population dynamics of the target species. The effectiveness of SIT is affected by the occurrence of genetically divergent 'strains' of the target species across the geographic area under control because this is detrimental to the mating success of sterile flies [##UREF##10##19##,##UREF##12##21##,##UREF##13##22##]. Specifically, high genetic divergence may reflect differences in behaviour and/or morphological characteristics that result in mating incompatibility among populations of the target species [##UREF##12##21##,##UREF##14##23##]. Thus, target populations that show low genetic divergence are not likely to show reproductive isolation and influence the success of a particular sterile strain.</p> <p>The Galápagos archipelago offers a unique system to examine the population genetics of an introduced avian parasite that causes severe fitness costs and that is still within a relatively early phase of invasion. We collected parasites in 2004, 2005 and 2006 from three islands of the Galápagos. Using mitochondrial data, we firstly determine whether the three island populations from which we sampled are of the one fly species. We then use microsatellite data to examine gene flow within and among islands to: (1) determine whether dispersal and genetic divergence are occurring among islands and between habitats within islands (wet highlands, arid lowlands), (2) determine the presence of population bottlenecks resulting from the invasion process, and (3) determine whether inter-island genetic differentiation may be of concern to the potential success of an archipelago-wide SIT program for controlling <italic>P. downsi</italic>.</p>
<title>Methods</title> <title>Study species</title> <p><italic>Philornis downsi </italic>(family Muscidae; subfamily Azeliinae; tribe Reinwardtiini) is a semi-haematophagous obligate avian parasite in its three larval stages, whereas adult flies are non-parasitic and feed on organic matter [##UREF##5##13##]. Adults lay eggs inside the nares of newly hatched nestlings (usually at one to three days old), which hatch into first instar larvae [##REF##16899139##17##,##UREF##15##24##]. Second and third instar larvae attach externally and feed on nestling blood and tissues over four to six days [##UREF##5##13##]. Most larvae of <italic>P. downsi </italic>appear to reach their third instar phase at the time of host fledging. The larvae pupariate at the base of the nesting material and remain for approximately two weeks before emerging as adult flies [##UREF##5##13##,##UREF##16##25##].</p> <title>Study area and sample collection</title> <p><italic>Philornis downsi </italic>were collected from three islands of the Galápagos: Santa Cruz (986 km<sup>2</sup>; 0° 37'S, 90° 21'W), Floreana (173 km<sup>2</sup>, 1° 28'S, 90° 48'W), and Isabela (4588 km<sup>2</sup>, 0° 58'S, 90° 58'W). Fly samples were collected from nests during the January to March finch breeding season in 2004, 2005 and 2006 from two contrasting habitats, the arid lowlands (0–100 m asl) and the humid highlands (300–600 m asl) (Table ##TAB##0##1##) [see also [##UREF##17##26##,##UREF##18##27##]]. The lowlands are characterised by low rainfall, and are dominated by the trees <italic>Acacia macracantha</italic>, <italic>Bursera graveolens</italic>, <italic>Croton scouleri</italic>, <italic>Opuntia </italic>spp., <italic>Pisonia floribunda</italic>, and <italic>Zanthoxylum fagara </italic>[##UREF##4##12##]. In contrast, the highlands have much higher rainfall [##UREF##19##28##,##UREF##20##29##], abundant moss and lichen, and are dominated by the endemic tree <italic>Scalesia pedunculata</italic>, or <italic>S. cordata </italic>(Asteraceae) on Isabela Island.</p> <p>We sampled from one site in each habitat on both Floreana (lowlands, adjacent to the town of Puerto Velasco Ibarra: 1° 16'S, 90° 29'W; highlands, base of Cerro Pajas: 1° 17'S, 090° 27'W) (Figure ##FIG##0##1##) and Isabela (lowlands: adjacent to town of Puerto Villamil: 0° 57'S, 91° 00'W; highlands: 0° 50'S, 91° 01'W), while on Santa Cruz we sampled from three sites in the lowlands: (1) Garrapatero: 0° 39'S, 90° 28'W; (2) Itabaca: 0° 29'S, 90° 17'W; (3) Punta Estrada, near Puerto Ayora: 0 ° 45'S, 90° 18'W, and one site in the highlands (Los Gemelos: 0° 37'S, 90° 22'W) (Figure ##FIG##0##1##). All sample sites were approximately 2000–4000 m<sup>2</sup>, except for the highland site on Isabela, where our sample site was only 100 m<sup>2 </sup>because habitat fragmentation has reduced the <italic>Scalesia </italic>forest to small remnant patches. The distance between highland and lowland sites was much shorter on Floreana (3–5 km) than on Santa Cruz and Isabela (both 15–25 km), while on Santa Cruz, the distance between all four sites (1 highland, 3 lowland) varied between 15 and 27 km. Data were obtained from all three islands in 2004, from just Floreana in 2005, and from Santa Cruz and Floreana in 2006 (Table ##TAB##0##1##).</p> <p>For the purpose of our study, larvae, puparia and puparia cases were sampled from 64 bird nests of five Darwin finch species (<italic>Geospiza fuliginosa</italic>, n = 25, <italic>Geospiza fortis</italic>, n = 15, <italic>Camarhynchus parvulus</italic>, n = 3, <italic>Camarhynchus pauper</italic>, n = 4; <italic>Cactospiza pallida</italic>, n = 1), while one nest was opportunistically sampled from each of the Galápagos mockingbird (<italic>Nesomimus parvulus</italic>) and the yellow warbler (<italic>Dendroica petechia aureola</italic>). Fourteen recently fledged nests were sampled for <italic>P. downsi </italic>where the finch species was unknown. GPS coordinates were recorded at each nest location. Inactive nests were collected and sealed in individual plastic bags and later dismantled for counting of <italic>P. downsi </italic>individuals. All flies were immediately preserved in 95% ethanol.</p> <title>DNA extraction and microsatellite typing</title> <p>DNA extraction was carried out using the salting out procedure described in [##UREF##21##30##] with the exception that all samples (3 mm<sup>2 </sup>tissue from each individual) were homogenised and washed three times in 10 mm TRIS prior to digestion with Proteinase K to remove traces of ethanol, excess lipids, and other potential contaminants. Across all three islands, 1012 <italic>P. downsi </italic>individuals (larvae and pupae) were genotyped (Table ##TAB##0##1##) using eight microsatellite markers [##UREF##22##31##]: Pd1 [GenBank: <ext-link ext-link-type="gen" xlink:href="EF608562">EF608562</ext-link>] Pd2 [<ext-link ext-link-type="gen" xlink:href="EF608556">EF608556</ext-link>], Pd4 [<ext-link ext-link-type="gen" xlink:href="EF608557">EF608557</ext-link>], Pd6 [<ext-link ext-link-type="gen" xlink:href="EF608564">EF608564</ext-link>], Pd7 [<ext-link ext-link-type="gen" xlink:href="EF608558">EF608558</ext-link>], Pd8 [<ext-link ext-link-type="gen" xlink:href="EF608555">EF608555</ext-link>], Pd9 [<ext-link ext-link-type="gen" xlink:href="EF608561">EF608561</ext-link>], Pd10 [<ext-link ext-link-type="gen" xlink:href="EF608563">EF608563</ext-link>]. Multiplex PCR conditions were followed as described in Dudaniec et al. [##UREF##22##31##]. Samples were genotyped on an ABI 3730 capillary electrophoresis DNA analyser (Applied Biosystems). A fluorescently labeled size standard (GS500 (-250) LIZ) was run with the samples and alleles were scored using <sc>GENEMAPPER</sc> version 3.7 (Applied Biosystems). To minimise and estimate genotyping error, each run of the DNA analyser contained eight repeated samples and a control sample run each time. In total, this resulted in 70 individual samples (14.5% of all samples genotyped) being re-amplified and genotyped at least once.</p> <title>Mitochondrial DNA sequencing</title> <p>An 822-bp region of the 3' end of the <italic>CO1 </italic>gene was amplified in five <italic>P. downsi </italic>individuals collected from Santa Cruz (1 highlands), Floreana (1 highlands, 1 lowlands), and Isabela (1 highlands, 1 lowlands). Samples were amplified using primers M202 (forwards, C1-J-1751 [##UREF##23##32##]) and M70 (reverse, UEA10 [##REF##8799733##33##]). Amplifications were performed in 10× <italic>Taq</italic>Gold buffer, 25 mM MgCl<sub>2</sub>, 10 mM total dNTP's, 200 nM each primer, 0.2 U <italic>Taq</italic>Gold polymerase, and 10–50 ng DNA. Amplification conditions were an initial denaturation at 94°C for 9 min, followed by 34 cycles of 94°C for 45 s, 55°C for 45 s, 72°C for 1 min, with a final extension of 72°C for 6 min. Sequencing was performed using the ABI Prism™ Big Dye Terminator Cycle sequencing kit (Applied Biosystems) according to the manufacturer's instructions. Products were sequenced on ABI 3700 (version 3.7) automated DNA sequencers. SeqEd (version 1.0.3) (Applied Biosystems) was used to edit chromatogram files to determine bi-directional consensus sequences and to manually align sequences across samples.</p> <title>Allele frequencies and data set construction</title> <p>We calculated allele frequencies using <sc>RELATEDNESS</sc> 5.0.8 [##UREF##24##34##] by randomly selecting one individual per sample (n = 64) to eliminate the possibility of including related individuals (a sample is defined as all <italic>P. downsi </italic>individuals collected from a single bird's nest). Exact tests were performed for each microsatellite locus to test deviation from Hardy-Weinberg equilibrium using <sc>GENEPOP</sc>[##UREF##25##35##]. All loci were in Hardy Weinberg equilibrium after sequential Bonferroni correction [##UREF##26##36##] and these allele frequencies were used for all further analyses. Genetic relatedness among <italic>P. downsi </italic>offspring within nests of Darwin's finches is low, and the individuals found within each nest are produced by up to approximately five ovipositing females that have each mated with between one and five males (as found by sib-ship reconstruction analysis by Dudaniec et al. in review). To eliminate the effect of sibs in the data, we selected unrelated individuals that were identified using the sib-ship reconstruction method implemented in the program <sc>COLONY</sc> 1.2 [##REF##15126412##37##]. Each sample of <italic>P. downsi </italic>individuals taken from an independent bird nest was run in <sc>COLONY</sc> 1.2, which uses a maximum likelihood method that partitions individuals into pure full-sib families (i.e. monogamous female parent), or full-sib families nested within half-sib families (i.e. polyandrous female parent) using progeny genotypes without known parental genotypes [##REF##15126412##37##,##REF##12969482##38##]. Three runs were performed per sample with different random seed numbers (12, 80, and 243) to ensure data convergence, and a conservative error rate of 5% was implemented based on evidence from the re-genotyping of 70 individuals, in which genotyping error ranged from 0–5% across loci.</p> <p>We selected one individual per reconstructed maternal family (i.e. one family = the offspring assigned to one putative female parent). In nested-half sib families (i.e. one mother, multiple fathers), individuals were only selected from full sib families with the largest number of members that had the highest posterior probability. Only individuals genotyped at all eight loci were included in the analysis and individuals were not sampled from families that contained Class I or Class II typing errors (identified by <sc>COLONY</sc> 1.2) [##REF##15126412##37##]. These criteria resulted in a sample size of 158 individuals sampled from 63 bird nests (with between one and six unrelated individuals per sample) (Table ##TAB##0##1##).</p> <p>To examine the probability that two randomly selected individuals from the same population will have the same multi-locus genotype, a Probability of Identity (PI) analysis was performed using <sc>GIMLET</sc>[##UREF##27##39##]. The output is a cumulative multi-locus PI value, estimated both with and without sample size correction. PI values were calculated for the dataset of 158 individuals using equations of unbiased PI, which assumes that individuals are unrelated, and PI for sibs, which assumes that all individuals are siblings [##UREF##27##39##].</p> <title>Inter-island genetic differentiation</title> <p>Heterozygosity, and pairwise <italic>F</italic>st [##UREF##28##40##] was calculated to examine genetic differentiation between islands (Santa Cruz, Floreana, Isabela) using <sc>MICROSATELLITE ANALYSER</sc> (MSA) 4.05 [##UREF##29##41##]. Genotypic differentiation was tested between islands using option 3 with 10 0000 Markov chain iterations in <sc>GENEPOP</sc>[##UREF##25##35##]. P-values for multiple tests were adjusted using sequential Bonferroni correction [##UREF##26##36##]. The AMOVA method [##REF##1644282##42##] was conducted in <sc>GENALEX</sc> version 6 [##UREF##30##43##] to partition the total genetic variation into three levels: among islands, among individuals, and within individuals using the Codom-genotypic distance calculation and 9999 permutations.</p> <title>Population bottleneck analysis</title> <p>Recently colonised species may experience a population bottleneck, resulting in a reduction in the number of alleles and expected heterozygosity at polymorphic loci. However, alleles may be lost at a faster rate than the loss of heterozygosity, so observed heterozygosity is higher than the expected heterozygosity at equilibrium [##UREF##31##44##]. The program <sc>BOTTLENECK</sc> version 1.2.02 [##REF##8978083##45##] was used to test for the presence of a recent population bottleneck for <italic>P. downsi </italic>by analysing within-population heterozygosity and allele frequency using the constructed dataset of individuals sampled from all islands. Both the stepwise mutation model (SMM) and two-phase model of mutation (TPM) were used, with the latter model being considered the most appropriate for microsatellites. The variance for the TPM was set at 5% and the proportion of SMM in TPM was set at 95% [##UREF##32##46##]. To determine differences in gene diversity across loci, the Wilcoxin sign-rank test was used as recommended for data sets with less than 20 loci (with 10 000 permutations) [##REF##8978083##45##]. We also examined the allele frequency distribution in order to see whether it is approximately L-shaped (as expected under mutation-drift equilibrium) or not (indicating that a recent bottleneck has provoked a mode shift) as described in [##UREF##33##47##].</p> <title>Genetic structure among islands</title> <p>For inferring genetic structure among the three sampled islands, we conducted two complementary individual-based Bayesian clustering analyses using <sc>STRUCTURE</sc> 2.1 [##REF##10835412##48##] and the landscape genetics program <sc>GENELAND</sc>[##UREF##34##49##] without <italic>a priori </italic>knowledge of population units and limits. Both software packages implement a Bayesian clustering method that uses a MCMC technique to define the number of populations in a sample that are at Hardy Weinberg Equilibrium. The methods implemented in these two programs differ in that <sc>GENELAND</sc> determines the optimal number of populations or 'clusters' and then allocates individuals (probabilistically) to these clusters using geographic coordinates, whereas <sc>STRUCTURE</sc> carries out the allocation sequentially for different numbers of clusters, and then flags the number of clusters with the highest likelihood [##REF##16924258##50##]. In <sc>STRUCTURE</sc>, the following run parameters were used: admixture without population information used, correlated allele frequency model, a burn-in period of 100 000 simulations followed by a run length of 1 million Markov Chain Monte Carlo (MCMC) simulations and three iterations for each number of potential clusters (defined as <italic>k </italic>= 1–5) to check for consistency of results. Estimation of <italic>k </italic>was taken to be the values of <italic>k </italic>with the highest Pr (X|<italic>k</italic>).</p> <p>In contrast to <sc>STRUCTURE</sc>, the algorithm implemented in <sc>GENELAND</sc>[##UREF##34##49##] is considered to be a powerful clustering method under conditions of low genetic differentiation among populations [##REF##14871360##51##,##REF##16629819##52##]. The model infers genetic discontinuities between populations in space from multilocus genotypes obtained from geo-referenced individuals [##UREF##34##49##,##REF##17651495##53##]. All individuals from the same sample (i.e. same bird nest) were allocated the same GPS coordinates. GPS coordinates were available for 57/63 nests (Santa Cruz: n = 18; Floreana: n = 36; Isabela: n = 3) (138 individuals in total). Samples for which GPS coordinates were missing were excluded from the analysis. To firstly infer the number of genetic clusters (<italic>k</italic>) in our data set, we used the Dirichlet model, which assumes independent allele frequencies with the following parameters: 1000 000 MCMC iterations, uncertainty attached to spatial coordinates = 0, variable number of populations = TRUE, minimum <italic>k </italic>= 1, maximum <italic>k </italic>= 5, and spatial information included in the model = TRUE. This procedure was performed three times to establish consistency of <italic>k </italic>across runs. The established <italic>k </italic>was then run five times to check the consistency of individual assignment to the inferred populations across runs. The same parameters were used but <italic>k </italic>was fixed at the modal number found in the first analyses. These five runs were post-processed (with a burn-in of 1000 × 100 iterations) to obtain posterior probabilities of population membership for each individual. Consistency of results across the five runs was checked visually.</p> <p>Inferred populations were further examined for heterozygosity, allelic richness (corrected for sample size), observed (H<sub>o</sub>) and expected (H<sub>e</sub>) heterozygosity, inbreeding coefficients (<italic>F</italic>is), and genetic differentiation (estimated using <italic>F</italic>st) using <sc>FSTAT</sc> v. 2.9.3 [##UREF##35##54##].</p>
<title>Results</title> <title>Mitochondrial sequencing</title> <p>Sequences of the <italic>CO1 </italic>mitochondrial gene fragment in five individuals across islands showed almost no variation, with two individuals (one Santa Cruz highland and one Floreana highland) having an identical single nucleotide substitution (T-G). This supports the existence of one sampled species across the three islands.</p> <title>Genetic diversity and differentiation</title> <p>Probability of identity (PI) analyses showed that the microsatellite loci had sufficient power and resolution for the analyses. The unbiased PI value was 1.333<sup>-06</sup>, and the PI for sibs was 2.610<sup>-3</sup>. This equates to one individual in approximately 751 880 having a non-unique genotype where individuals are unrelated (unbiased), and one individual in approximately 383 individuals having a non-unique genotype if all individuals are siblings.</p> <p>The total number of alleles observed at each locus was as follows; Pd1 = 4; Pd2 = 3; Pd4 = 4; Pd6 = 5; Pd7 = 3; Pd8 = 4; Pd9 = 3; Pd10 = 3 (Table ##TAB##1##2##). There was significant genotypic differentiation across the three islands (Fisher's Exact method: <italic>X</italic><sup>2 </sup>= 72.75; df = 16; P &lt; 0.001). Mean allelic richness across loci was almost identical on each island (Santa Cruz: 3.50; Floreana 3.63; Isabela: 3.5) and the range of observed heterozygosity across loci was also similar (Santa Cruz: 0.45–0.70; Floreana: 0.45–0.73; Isabela: 0.44–0.89). The number and size of alleles from each island population were the same at each locus with two exceptions: there was a unique allele at locus Pd6 on Isabela (allele frequency = 0.055), and at locus Pd7 on Floreana (allele frequency = 0.012), which were each detected only in a single individual. Pairwise <italic>F</italic>st analysis showed low, but significant levels of genetic differentiation between Santa Cruz and Floreana (<italic>F</italic>st = 0.02, P &lt; 0.02) Isabela and Floreana (<italic>F</italic>st = 0.04, P &lt; 0.02), but not between Santa Cruz and Isabela (<italic>F</italic>st = 0.01, P &gt; 0.1). The low genetic differentiation found between islands was reflected in an AMOVA, which showed that just 2% of the molecular variance was attributable to variation among islands (sum of squares (SS) = 14.23; df = 2; variance components (V) = 0.052), 4% among individuals (SS = 413.06; df = 155; V = 0.133), and 94% within individuals (SS = 385.5; df = 158; V = 2.44).</p> <title>Bottleneck analysis</title> <p>Combining individual from all islands (n = 158), a clear excess of heterozygosity (He) relative to the equilibrium heterozygosity (H<sub>eq</sub>) was observed, indicative of a population bottleneck under the TPM model (Wilcoxon sign-rank test; P &lt; 0.01) and under the SMM model (P &lt; 0.01). A mode-shift distortion in the distribution of allele frequencies was evident (Figure ##FIG##1##2##).</p> <title>Bayesian clustering analysis</title> <p>Individual-based cluster analysis using <sc>STRUCTURE</sc>[##REF##10835412##48##] did not detect any genetic structuring in <italic>P. downsi </italic>collected across the three islands (Figure ##FIG##2##3a##), with individual assignment being evenly proportioned across variable numbers of <italic>k</italic>. This implies high levels of inter-island ancestry brought about by frequent dispersal and subsequent gene flow across the three islands sampled. However, when incorporating geographic coordinates of sampling locations into Bayesian analyses using <sc>GENELAND</sc>[##UREF##34##49##], two distinct genetic clusters were consistently found across runs (Figures ##FIG##2##3b## and ##FIG##3##4##). The first cluster includes all individuals sampled from Santa Cruz and Isabela Islands (n = 62), while the second cluster includes all individuals sampled from Floreana Island (n = 76). Assignment probabilities were between 0.98 and 1.0 across all individuals.</p> <title>Genetic diversity and differentiation among inferred clusters</title> <p>The two clusters identified by <sc>GENELAND</sc> displayed comparable genetic diversity with regard to allelic richness and differed slightly in heterozygosity across loci (Tables ##TAB##1##2## and ##TAB##2##3##). Although two clusters were detected, measures of genetic differentiation (<italic>F</italic>st) between them demonstrated the low divergence between individuals on Floreana Island and those on Santa Cruz and Isabela (<italic>F</italic>st = 0.024; 95% Confidence Interval (CI) = 0.014 – 0.034; P &lt; 0.05). Tests of departure from HW equilibrium showed no significant deviation in either of the two clusters across all loci.</p>
<title>Discussion</title> <p>In combination with the microsatellite data, our mitochondrial findings are consistent with there being one species of <italic>Philornis </italic>on the islands from which we sampled. A population bottleneck was detected in the entire sample of individuals from the three islands, which is consistent with the pattern expected from an invasive, recently colonised species [##UREF##2##4##, ####UREF##3##5##, ##REF##16858586##6####16858586##6##]. We report low genetic differentiation between island populations of the invasive avian parasite <italic>P. downsi </italic>on the Galápagos archipelago. Fly populations on Santa Cruz, Floreana, and Isabela showed strong evidence for high inter-island gene flow. However, low levels of divergence were detected between individuals from Floreana Island and those from Santa Cruz and Isabela when incorporating geographic sampling information. The molecular variance was mainly explained at the level of individuals, and not by island, which further demonstrates the low genetic differentiation between islands. Bayesian clustering analysis with geographic data assigned individuals to two genetic clusters, one comprising individuals from Santa Cruz and Isabela, and the second comprising all individuals from Floreana Island (Table ##TAB##2##3##, Figure ##FIG##3##4##). This might indicate that gene flow in <italic>P. downsi </italic>between Floreana and the other islands is restricted to some extent, or that this island underwent a distinct founding process. Pairwise <italic>F</italic>st between the three islands further indicated that flies on Floreana may be slightly genetically divergent from flies on the other two islands.</p> <p>The Bayesian clustering method implemented in <sc>STRUCTURE</sc> is considered to be best able to infer correct individual assignments when genetic differentiation between populations is well defined [##REF##10835412##48##]. Furthermore, the ability to distinguish the source of an individual decreases under conditions of high dispersal and associated low genetic differentiation [##REF##14871360##51##,##REF##10581301##55##]. The level of genetic differentiation (<italic>F</italic>st) between populations is found to be a useful predictor of the performance of assignment methods [##REF##10581301##55##]. In the current study, the inability of <sc>STRUCTURE</sc> to confidently assign individuals to any cluster with certainty may reflect the lack of power to do so due to the low genetic differentiation (i.e. <italic>F</italic>st) between sampling locations. Thus, we conclude there was an insufficient signal in the data to confidently assign individuals under the model of Pritchard et al. [##REF##10835412##48##], despite reasonably high PI values across loci. Our results are therefore testament that taking the spatial context of individuals into account improved the efficiency of our analysis, as found by Fontaine et al. [##REF##17651495##53##]. Verifying the usefulness of <sc>STRUCTURE</sc> to assign individuals correctly where genetic differentiation is low and dispersal is common requires further study using empirical field data [##REF##14871360##51##,##REF##10581301##55##].</p> <p>The current study lacks genetic data from mainland <italic>P. downsi </italic>populations and data from all islands of the Galápagos where <italic>P. downsi </italic>occurs, which will be necessary for a detailed examination of founder effects, bottlenecks, introduction events and colonisation pathways. Thus, without knowing where <italic>P. downsi </italic>populations originally came from, or where they most recently arrived on the Galápagos archipelago, a comprehensive invasion history can not be constructed on a demographic or evolutionary scale [##REF##10234242##11##,##REF##18028178##56##]. However, our findings lay the foundation for a more thorough understanding of the process of <italic>P. downsi </italic>invasion on the Galápagos archipelago. It is possible that <italic>P. downsi </italic>arrived on Ecuadorian cargo ships that were transporting fruit to the islands for human consumption [##UREF##36##57##,##UREF##37##58##], while it is also suggested that the fly came with imported pigeons (discussed in [##UREF##38##59##]). Strong winds and air currents present during El Niño events on the Galápagos are believed to contribute to insect dispersal between islands [##UREF##39##60##], while transport of humans and materials is also suspected to aid inter-island insect dispersal. In four other invasive insect species, the dates of colonisation on each island suggest a wind-mediated southeast to northwest direction of colonisation across the islands [##UREF##36##57##]. Such patterns remain unexplored for <italic>P. downsi</italic>.</p> <title>Invasion processes</title> <p>Recently colonised invaders are often subject to a reduction in genetic variation and population bottlenecks because populations are not in genetic equilibrium [##UREF##2##4##, ####UREF##3##5##, ##REF##16858586##6####16858586##6##]. We provide evidence for a population bottleneck in <italic>P. downsi </italic>across the three islands examined, which could be due to a small founding population, low immigration rates, or few introduction events [##REF##16858586##6##,##UREF##40##61##]. The low allelic diversity across loci and population bottleneck in <italic>P. downsi </italic>is further evidence for a small effective population size upon initial colonisation. However, the occurrence of multiple introductions can not be excluded, particularly in the absence of comparisons with potential source populations (e.g. from Ecuador, Trinidad, or Brazil) [##UREF##5##13##]. Despite the presence of a population bottleneck and the (most likely related) low genetic diversity in <italic>P. downsi</italic>, the fly has clearly succeeded at establishing and spreading itself across the archipelago in high numbers.</p> <p>Recently established species may persist at low and possibly undetectable numbers before becoming noticeably abundant and invasive years or decades later [##UREF##3##5##], which may reflect the lag time (i.e. the time between arrival and spread) observed in many species that become invasive [##REF##12038524##62##]. This scenario seems likely concerning the invasion of <italic>P. downsi </italic>on Galápagos because the fly was not detected in finch nests and identified until 1997 [##UREF##41##63##], despite the recent discovery of specimens found in collections made in 1964 [##UREF##5##13##,##REF##16899139##17##]. The parasite has since spread successfully and in high numbers across the archipelago (11 of 13 major islands) [##UREF##38##59##], indicating that any lag period that took place has passed. Yet it is unknown how recently each island was colonised and thus, whether particular island populations are undergoing a lag period that would favour the success of an immediate eradication effort (discussed in [##UREF##42##64##]).</p> <p>Ecological [##UREF##4##12##,##UREF##19##28##,##UREF##38##59##] findings do not support the current existence of a lag period and indicate that <italic>P. downsi </italic>has spread successfully in at least 12 avian host species on the Galápagos Islands [##UREF##4##12##,##UREF##5##13##]. In the current study, we provide evidence that the <italic>P. downsi </italic>population on Floreana Island has detectable levels of genetic differentiation when compared with two other island populations, which might be the result of a separate introduction event(s) or colonisation pattern. A wider geographic sample of locations across habitats and islands is needed to examine this more definitively in combination with a larger number of highly polymorphic genetic markers. However, it is clear that <italic>P. downsi </italic>populations generally have high connectivity between islands or high shared ancestry, although variation in population processes (e.g. rates of dispersal, colonisation histories) between particular islands may allow for low levels of inter-island genetic differentiation.</p> <title>Absence of local genetic divergence</title> <p>Local populations are expected to evolve adaptive differences in response to differing environmental conditions [##UREF##42##64##]. The lack of genetic structure in <italic>P. downsi </italic>on the Galápagos archipelago may reflect the estimated short time period since the flies' introduction (~40 years ago) [##UREF##1##3##] such that populations have not yet diverged since colonisation. We document no genetic structure according to habitat type across islands, which implies high levels of fly dispersal between the two habitats. Across islands however, differences in host diversity and distribution, ecological variables, or colonisation history may result in genetic divergence due to genetic drift, as was evident from the low genetic differentiation we document on Floreana Island.</p> <p>Fly populations may show rapid evolution with geographic cline, as shown by Huey et al. [##REF##10634786##65##] who found increased wing length with latitude in <italic>Drosophila subobscura</italic>, just two decades after its introduction into North America. The evidence we present for high gene flow between habitats implies that morphological variation in <italic>P. downsi </italic>is unlikely, though other insect species on Galápagos show morphological variation and genetic differentiation between habitats and islands of the archipelago [##REF##18028178##56##,##UREF##43##66##,##REF##17304798##67##]. Clinal variation in morphology (and evidence for low dispersal) was also found for <italic>Bulimilus </italic>land snails on Galápagos [##REF##17236423##68##] and Darwin's small ground finch [##UREF##17##26##].</p> <title>Implications for control: the sterile insect technique (SIT)</title> <p>The use of SIT to control <italic>P. downsi </italic>on the Galápagos Islands is perhaps the most appropriate method for eradicating an invasive fly within this ecologically fragile island ecosystem. SIT is a non-disruptive method as it does not introduce toxic or foreign chemicals into the environment, it is species specific, and does not introduce new genetic material into populations because the released organisms are not self-replicating [##UREF##10##19##,##UREF##44##69##].</p> <p>The effectiveness of SIT is affected by population genetic differentiation within the target species because the occurrence of undetected sub-species or strain differentiation across geographic populations can be detrimental to widespread sterile male release [##UREF##10##19##]. Reinfestation of parasitic flies in SIT treated regions have been explained by genetic differentiation in the target species among allopatrically separated populations that may be experiencing reproductive isolation [e.g. [##REF##9394464##70##]]. It is therefore of great advantage to use molecular genetic techniques for species characterisation and to examine population genetic structure prior to establishing large-scale sterile male release programs. We show that gene flow in <italic>P. downsi </italic>within and between three islands of the Galápagos is high, and unlikely to result in reproductive isolation. Thus, release of a single sterile strain of <italic>P. downsi </italic>could effectively suppress and eradicate the fly across the archipelago. Captive breeding experiments of adult <italic>P. downsi </italic>from multiple island populations are necessary to determine this with high confidence.</p>
<title>Conclusion</title> <p>The wide habitat range and high dispersal capacity of <italic>P. downsi </italic>highlights the significant threat that this parasite poses to the Galápagos endemic avifauna. Our findings are concordant with the prediction that parasites with low host specificity [##UREF##19##28##], good dispersal ability and horizontal transmission will show low population genetic structure and differentiation [##REF##18028178##56##]. Ideally, it is best to eradicate invasive species before they become adapted to the local environment in which they have colonised and prior to repeated invasions with the aid of strict quarantine practices [##REF##15602569##7##]. For <italic>P. downsi</italic>, this window of opportunity appears to have passed, prompting the need for a long-term eradication program combined with sustained quarantine and monitoring practices.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>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, <italic>Philornis downsi</italic>, 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 <italic>P. downsi </italic>between three islands and across habitats (highlands, lowlands) and examine for the presence of population bottlenecks. We also examine variation at the mitochondrial gene <italic>CO1 </italic>across islands to establish if cryptic species were present.</p> <title>Results</title> <p>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 <italic>P. downsi </italic>population was found to have undergone a population bottleneck.</p> <title>Conclusion</title> <p><italic>Philornis downsi </italic>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 <italic>P. downsi</italic>. The high dispersal capacity and wide habitat use of <italic>P. downsi </italic>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 <italic>P. downsi </italic>on Galápagos, such as the sterile insect technique.</p>
<title>Authors' contributions</title> <p>RYD and MGG performed all molecular laboratory work. RYD performed all genotyping, data analysis and drafted the manuscript. MGG and SD participated in data analysis and interpretation. SK participated in the design and coordination of the study and provided funding. SK and RYD conducted field work and collection of samples. All authors read and approved the final manuscript.</p>
<title>Acknowledgements</title> <p>We thank the Galápagos National Park Service and the Charles Darwin Research Station for the opportunity to work on the Galápagos Archipelago. We thank the following people: Birgit Fessl, Carlos Vinueza, David Wiedenfeld, Gustavo Jiménez, Carlos Santos, Rebekah Christensen, Jody O'Connor, and Jeremy Robertson for their dedicated field assistance, Bradley Sinclair for fly identification and Kathy Saint for molecular technical advice. We thank the community of Floreana Island for their invaluable support throughout the project. TAME airlines provided reduced airfares to the Galápagos. This study was funded by a Flinders University Establishment Grant to SK, the Galápagos Conservation Trust, the American Bird Conservancy and Conservation International.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Map of Santa Cruz, Floreana, and Isabela Islands with sampling locations</bold>. Sampling sites on Santa Cruz: S1 = highlands; S2 = lowlands, Punta Estrada; S3 = lowlands, Garrapatero; S4 = lowlands, Itabaca. On Floreana and Isabela, one site each in the lowlands (L) and highlands (H) are indicated.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Distribution of allele frequencies indicating a mode-shift</bold>. Bars represent the proportion of alleles found in each allele frequency class. Deviation from an L-shaped distribution is indicative of a mode-shift in allele frequency due to a recent genetic bottleneck.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Estimated number of populations from <sc>STRUCTURE</sc> (a) and <sc>GENELAND</sc> (b) analyses</bold>. (a) Mean ( ± SD) probabilities of the data (LnPr [X|<italic>k</italic>]) over three replicate <sc>STRUCTURE</sc> runs plotted as a function of the putative number of clusters (<italic>k</italic>); (b) Posterior density distribution of the number of clusters estimated from <sc>GENELAND</sc> analysis in three replicates.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Genetic assignment of <italic>P. downsi </italic>individuals across three islands using Bayesian clustering analysis</bold>. Two genetic clusters are identified: (a) including all individuals from Santa Cruz (n = 62) (bottom left) and Isabela (n = 9) (centre top), and (b) all individuals from Floreana Island (n = 76) (bottom right). Black dots represent independent geographic sampling points (i.e. location of bird nests). Note that two geographic sampling points on Isabela Island were within 5 m of each other and are not distinguishable.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Sample sizes of bird nests and <italic>P. downsi </italic>individuals.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="center" colspan="2"><bold># Bird nests sampled  </bold><break/><bold>for <italic>P. downsi</italic></bold></td><td align="center" colspan="2"><bold># Individuals </bold><break/><bold>analysed</bold></td></tr><tr><td colspan="2"/><td colspan="2"><hr/></td><td colspan="2"><hr/></td></tr><tr><td align="left"><bold>Year</bold></td><td align="left"><bold>Island</bold></td><td align="center">Highland</td><td align="center">Lowland</td><td align="center">Highland</td><td align="center">Lowland</td></tr></thead><tbody><tr><td align="left">2004</td><td align="left">Santa Cruz</td><td align="center">3</td><td align="center">18</td><td align="center">11</td><td align="center">51</td></tr><tr><td/><td align="left">Floreana</td><td align="center">1</td><td align="center">4</td><td align="center">1</td><td align="center">9</td></tr><tr><td/><td align="left">Isabela</td><td align="center">2</td><td align="center">1</td><td align="center">7</td><td align="center">2</td></tr><tr><td align="left">2005</td><td align="left">Floreana</td><td align="center">11</td><td align="center">-</td><td align="center">28</td><td align="center">-</td></tr><tr><td align="left">2006</td><td align="left">Santa Cruz</td><td align="center">-</td><td align="center">2</td><td align="center">-</td><td align="center">5</td></tr><tr><td/><td align="left">Floreana</td><td align="center">11</td><td align="center">10</td><td align="center">30</td><td align="center">14</td></tr><tr><td colspan="6"><hr/></td></tr><tr><td/><td align="left"><bold>Total</bold></td><td align="center"><bold>28</bold></td><td align="center"><bold>35</bold></td><td align="center"><bold>77</bold></td><td align="center"><bold>81</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Allele frequencies for eight microsatellite loci in <italic>P. downsi </italic>within two genetic clusters.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="left"><bold>Santa Cruz </bold><break/><bold>and Isabela</bold></td><td align="left"><bold>Floreana</bold></td></tr></thead><tbody><tr><td align="left"><italic>Pd1</italic></td><td align="left">315</td><td align="left">0.319588</td><td align="left">0.219512</td></tr><tr><td/><td align="left">323</td><td align="left">0.134021</td><td align="left">0.256098</td></tr><tr><td/><td align="left">325</td><td align="left">0.154639</td><td align="left">0.109756</td></tr><tr><td/><td align="left">327</td><td align="left">0.391753</td><td align="left">0.414634</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd2</italic></td><td align="left">230</td><td align="left">0.273196</td><td align="left">0.268293</td></tr><tr><td/><td align="left">236</td><td align="left">0.582474</td><td align="left">0.439024</td></tr><tr><td/><td align="left">240</td><td align="left">0.14433</td><td align="left">0.292683</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd4</italic></td><td align="left">252</td><td align="left">0.35567</td><td align="left">0.304878</td></tr><tr><td/><td align="left">254</td><td align="left">0.298969</td><td align="left">0.353659</td></tr><tr><td/><td align="left">256</td><td align="left">0.036082</td><td align="left">0.02439</td></tr><tr><td/><td align="left">258</td><td align="left">0.309278</td><td align="left">0.317073</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd6</italic></td><td align="left">241</td><td align="left">0.391753</td><td align="left">0.390244</td></tr><tr><td/><td align="left">251</td><td align="left">0.221649</td><td align="left">0.280488</td></tr><tr><td/><td align="left">259</td><td align="left">na</td><td align="left">0.012195</td></tr><tr><td/><td align="left">261</td><td align="left">0.139175</td><td align="left">0.109756</td></tr><tr><td/><td align="left">263</td><td align="left">0.247423</td><td align="left">0.207317</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd7</italic></td><td align="left">201</td><td align="left">0.371134</td><td align="left">0.439024</td></tr><tr><td/><td align="left">207</td><td align="left">0.27835</td><td align="left">0.268293</td></tr><tr><td/><td align="left">210</td><td align="left">0.340206</td><td align="left">0.292683</td></tr><tr><td/><td align="left">213</td><td align="left">0.010309</td><td align="left">na</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd8</italic></td><td align="left">349</td><td align="left">0.139175</td><td align="left">0.195122</td></tr><tr><td/><td align="left">353</td><td align="left">0.087629</td><td align="left">0.04878</td></tr><tr><td/><td align="left">357</td><td align="left">0.134021</td><td align="left">0.097561</td></tr><tr><td/><td align="left">365</td><td align="left">0.639175</td><td align="left">0.658537</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd9</italic></td><td align="left">194</td><td align="left">0.324742</td><td align="left">0.414634</td></tr><tr><td/><td align="left">200</td><td align="left">0.469072</td><td align="left">0.45122</td></tr><tr><td/><td align="left">217</td><td align="left">0.206186</td><td align="left">0.134146</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left"><italic>Pd10</italic></td><td align="left">189</td><td align="left">0.216495</td><td align="left">0.329268</td></tr><tr><td/><td align="left">191</td><td align="left">0.231959</td><td align="left">0.134146</td></tr><tr><td/><td align="left">193</td><td align="left">0.551546</td><td align="left">0.536585</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Genetic variation at the eight microsatellite loci for the two <italic>P. downsi </italic>populations inferred from cluster analysis in <sc>GENELAND</sc>.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="3"><bold>Santa Cruz and Isabela</bold></td><td align="center" colspan="3"><bold>Floreana</bold></td></tr><tr><td/><td colspan="3"><hr/></td><td colspan="3"><hr/></td></tr><tr><td align="left"><bold>Locus</bold></td><td align="left"><bold>A</bold></td><td align="left"><bold>H</bold><sub>o</sub><bold>/H</bold><sub>e</sub></td><td align="left"><bold><italic>F</italic></bold><bold>is</bold></td><td align="left"><bold>A</bold></td><td align="left"><bold>Ho/He</bold></td><td align="left"><bold>Fis</bold></td></tr></thead><tbody><tr><td align="left"><bold>Pd1</bold></td><td align="left">4.0</td><td align="left">0.63/0.71</td><td align="left">0.038</td><td align="left">4.0</td><td align="left">0.71/0.71</td><td align="left">0.091</td></tr><tr><td align="left"><bold>Pd2</bold></td><td align="left">3.0</td><td align="left">0.46/0.57</td><td align="left">0.114</td><td align="left">3.0</td><td align="left">0.68/0.66</td><td align="left">0.091</td></tr><tr><td align="left"><bold>Pd4</bold></td><td align="left">4.0</td><td align="left">0.67/0.69</td><td align="left">0.105</td><td align="left">3.0</td><td align="left">0.64/0.67</td><td align="left">0.037</td></tr><tr><td align="left"><bold>Pd6</bold></td><td align="left">5.0</td><td align="left">0.68/0.72</td><td align="left">0.090</td><td align="left">4.0</td><td align="left">0.73/0.72</td><td align="left">-0.012</td></tr><tr><td align="left"><bold>Pd7</bold></td><td align="left">3.0</td><td align="left">0.56/0.67</td><td align="left">0.131</td><td align="left">3.0</td><td align="left">0.73/0.66</td><td align="left">0.037</td></tr><tr><td align="left"><bold>Pd8</bold></td><td align="left">4.0</td><td align="left">0.48/0.55</td><td align="left">0.023</td><td align="left">4.0</td><td align="left">0.61/0.52</td><td align="left">0.034</td></tr><tr><td align="left"><bold>Pd9</bold></td><td align="left">3.0</td><td align="left">0.67/0.64</td><td align="left">-0.083</td><td align="left">3.0</td><td align="left">0.68/0.61</td><td align="left">-0.094</td></tr><tr><td align="left"><bold>Pd10</bold></td><td align="left">3.0</td><td align="left">0.57/0.61</td><td align="left">0.186</td><td align="left">3.0</td><td align="left">0.68/0.59</td><td align="left">-0.154</td></tr><tr><td align="left"><bold>Mulitlocus</bold></td><td align="left">3.6</td><td align="left">0.59/0.65</td><td align="left">0.074</td><td align="left">3.6</td><td align="left">0.67/0.65</td><td align="left">0.005</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>The number of nests and the number of individuals analysed (following construction of dataset comprising unrelated individuals) for each island and habitat (highland/lowland) across three islands, Santa Cruz, Floreana, and Isabela in 2004, 2005 and 2006.</p></table-wrap-foot> <table-wrap-foot><p>Clusters were inferred using landscape genetic analysis: Cluster 1 = Santa Cruz and Isabela Island (n = 62); Cluster 2 = Floreana Island (n = 76); na = not applicable because allele was absent.</p></table-wrap-foot> <table-wrap-foot><p>Sample size is 62 for the Santa Cruz and Isabela cluster and 76 for the Floreana cluster for all loci. A: allelic richness (estimated for a sample size of 75 individuals); Ho:observed heterozygosity; He: expected heterozygosity; Fis (inbreeding coefficients) were calculated after Weir and Cockerham [##UREF##28##40##]. None of the loci had a significant heterozygote deficiency or excess after sequential Bonferroni correction.</p></table-wrap-foot>
<graphic xlink:href="1472-6785-8-13-1"/> <graphic xlink:href="1472-6785-8-13-2"/> <graphic xlink:href="1472-6785-8-13-3"/> <graphic xlink:href="1472-6785-8-13-4"/>
[{"surname": ["O'Dowd", "Green", "Lake"], "given-names": ["DJ", "PT", "PS"], "article-title": ["Invasional 'meltdown' on an oceanic island"], "source": ["Ecol Letters"], "year": ["2003"], "volume": ["6"], "fpage": ["812"], "lpage": ["817"], "pub-id": ["10.1046/j.1461-0248.2003.00512.x"]}, {"surname": ["Causton", "Peck", "Sinclair", "Roque-Albelo", "Hodgson", "Landry"], "given-names": ["C", "SB", "BJ", "L", "CJ", "B"], "article-title": ["Alien insects: Threats and implications for conservation of Gal\u00e1pagos Islands"], "source": ["Ann Entomol Soc Am"], "year": ["2006"], "volume": ["99"], "fpage": ["121"], "lpage": ["143"], "pub-id": ["10.1603/0013-8746(2006)099[0121:AITAIF]2.0.CO;2"]}, {"surname": ["Azzurro", "Golani", "Bucciarelli", "Bernardi"], "given-names": ["E", "D", "G", "G"], "article-title": ["Genetics of the early stages of invasion of the Lessepsian rabbitfish "], "italic": ["Siganus luridus"], "source": ["J Exp Mar Biol Ecol"], "year": ["2006"], "volume": ["333"], "fpage": ["190"], "lpage": ["201"], "pub-id": ["10.1016/j.jembe.2005.12.002"]}, {"surname": ["Sakai", "Allendorf", "Holt", "Lodge", "Molofsky", "With", "Baughman", "Cabin", "Cohen", "Ellstrand", "McCauley", "O'Neil", "Parker", "Thompson", "Weller"], "given-names": ["AK", "FW", "JS", "DM", "J", "KA", "S", "RJ", "JE", "NC", "DE", "P", "IM", "JN", "SG"], "article-title": ["The population biology of invasive species"], "source": ["Ann Rev Ecol Systemat"], "year": ["2001"], "volume": ["32"], "fpage": ["305"], "lpage": ["332"], "pub-id": ["10.1146/annurev.ecolsys.32.081501.114037"]}, {"surname": ["Fessl", "Tebbich"], "given-names": ["B", "S"], "italic": ["Philornis downsi "], "source": ["Ibis"], "year": ["2002"], "volume": ["144"], "fpage": ["445"], "lpage": ["451"], "pub-id": ["10.1046/j.1474-919X.2002.00076.x"]}, {"surname": ["Dudaniec", "Kleindorfer"], "given-names": ["RY", "S"], "article-title": ["Effects of the parasitic flies of the  genus "], "italic": ["Philornis"], "source": ["Emu"], "year": ["2006"], "volume": ["106"], "fpage": ["13"], "lpage": ["20"], "pub-id": ["10.1071/MU04040"]}, {"surname": ["Fessl", "Kleindorfer", "Tebbich"], "given-names": ["B", "S", "S"], "article-title": ["An experimental study on the effects of an introduced parasite in Darwin's finches"], "source": ["Biol Conserv"], "year": ["2006"], "volume": ["127"], "fpage": ["55"], "lpage": ["61"], "pub-id": ["10.1016/j.biocon.2005.07.013"]}, {"surname": ["Kleindorfer", "Dudaniec"], "given-names": ["S", "RY"], "article-title": ["Increasing prevalence of avian poxvirus in Darwin's finches and its effect on male pairing success"], "source": ["J Avian Biol"], "year": ["2006"], "volume": ["37"], "fpage": ["69"], "lpage": ["76"], "pub-id": ["10.1111/j.2006.0908-8857.03503.x"]}, {"surname": ["Dudaniec", "Hallas", "Kleindorfer"], "given-names": ["RY", "G", "S"], "article-title": ["Blood and intestinal parasitism in Darwin's finches: positive and negative findings"], "source": ["Acta Zool Sinica"], "year": ["2005"], "volume": ["51"], "fpage": ["507"], "lpage": ["512"]}, {"surname": ["Dudaniec", "Kleindorfer", "Fessl"], "given-names": ["RY", "S", "B"], "article-title": ["Effects of the introduced ectoparasite "], "italic": ["Philornis downsi ", "Geospiza fuliginosa"], "source": ["Austral Ecol"], "year": ["2006"], "volume": ["31"], "fpage": ["88"], "lpage": ["94"], "pub-id": ["10.1111/j.1442-9993.2006.01553.x"]}, {"surname": ["Hendrichs", "Robinson", "Cayol", "Enkerlin"], "given-names": ["J", "AS", "JP", "W"], "article-title": ["Medfly areawide sterile insect technique programs for prevention, suppression or eradication: the importance of mating behaviour studies"], "source": ["Florida Entomol"], "year": ["2002"], "volume": ["85"], "fpage": ["1"], "lpage": ["13"], "pub-id": ["10.1653/0015-4040(2002)085[0001:MASITP]2.0.CO;2"]}, {"surname": ["Vreysen", "Hendrichs", "Enkerlin"], "given-names": ["MJ", "J", "WR"], "article-title": ["The sterile insect technique as a component of sustainable area-wide integrated pest management of selected horticultural insect pests"], "source": ["J Fruit Ornamental Plant Res"], "year": ["2006"], "volume": ["14"], "fpage": ["107"], "lpage": ["131"]}, {"surname": ["Cayol", "Coronado", "Taher"], "given-names": ["JP", "P", "M"], "article-title": ["Sexual compatability in medfly (Diptera: Tephritidae) from different origins"], "source": ["Florida Entomol"], "year": ["2002"], "volume": ["85"], "fpage": ["51"], "lpage": ["57"], "pub-id": ["10.1653/0015-4040(2002)085[0051:SCIMDT]2.0.CO;2"]}, {"surname": ["Pereira", "Silva", "Quintal", "Abreu", "Andrade", "Dantas"], "given-names": ["R", "N", "C", "R", "J", "L"], "article-title": ["Sexual performance of mass reared and wild Mediterranean fruit flies (Diptera: Tephritidae) from various origins of the Madeira Islands"], "source": ["Florida Entomol"], "year": ["2007"], "volume": ["90"], "fpage": ["10"], "lpage": ["14"], "pub-id": ["10.1653/0015-4040(2007)90[10:SPOMRA]2.0.CO;2"]}, {"surname": ["McInnes", "Lance", "Jackson"], "given-names": ["DO", "DR", "CG"], "article-title": ["Behavioural resistance to the sterile insect technique by the Mediterranean fruit fly (Diptera: Tephritidae) in Hawaii"], "source": ["Ann Entomol Soc Am"], "year": ["1996"], "volume": ["89"], "fpage": ["739"], "lpage": ["744"]}, {"surname": ["Muth"], "given-names": ["A"], "article-title": ["Control of "], "italic": ["Philornis downsi"], "source": ["Report for Department of Terrestrial Invertebrates"], "year": ["2007"], "publisher-name": ["Charles Darwin Research Station, Gal\u00e1pagos Islands"]}, {"surname": ["Teixeira", "Guimar\u00e3es JH, Papavero N"], "given-names": ["DM"], "article-title": ["Myiasis caused by obligatory parasites Ib. General observations on the biology of species of the genus "], "italic": ["Philornis "], "source": ["Myiasis in man and animals in the Neotropical region; Bibliographic Database"], "year": ["1999"], "publisher-name": ["S\u00e3o Paulo: Editora Pl\u00eaiade/FAPESP"], "fpage": ["71"], "lpage": ["96"]}, {"surname": ["Kleindorfer", "Chapman", "Winkler", "Sulloway"], "given-names": ["S", "TW", "H", "FJ"], "article-title": ["Adaptive divergence in contiguous populations of Darwin's Small Ground Finch ("], "italic": ["Geospiza fuliginosa"], "source": ["Evol Ecol Res"], "year": ["2006"], "volume": ["8"], "fpage": ["357"], "lpage": ["372"]}, {"surname": ["Kleindorfer"], "given-names": ["S"], "article-title": ["The ecology of clutch size variation in Darwin's small ground finch: comparison between lowland and highland habitats"], "source": ["Ibis"], "year": ["2007"], "volume": ["149"], "fpage": ["730"], "lpage": ["741"], "pub-id": ["10.1111/j.1474-919X.2007.00694.x"]}, {"surname": ["Dudaniec", "Kleindorfer", "Fessl"], "given-names": ["RY", "S", "B"], "article-title": ["Interannual and interspecific variation in intensity of the parasitic fly, "], "italic": ["Philornis downsi"], "source": ["Biol Conserv"], "year": ["2007"], "volume": ["139"], "fpage": ["325"], "lpage": ["332"], "pub-id": ["10.1016/j.biocon.2007.07.006"]}, {"surname": ["Kleindorfer"], "given-names": ["S"], "article-title": ["Nesting success in Darwin's small tree finch, "], "italic": ["Camarhynchus parvulus "], "source": ["Anim Behav"], "year": ["2007"], "volume": ["74"], "fpage": ["795"], "lpage": ["804"], "pub-id": ["10.1016/j.anbehav.2007.01.020"]}, {"surname": ["Miller", "Dykes", "Polesky"], "given-names": ["SA", "DD", "HF"], "article-title": ["A simple salting out procedure for extracting DNA from human nucleated cells"], "source": ["Nuc Acids Res"], "year": ["1988"], "volume": ["16"], "fpage": ["1215"], "pub-id": ["10.1093/nar/16.3.1215"]}, {"surname": ["Dudaniec", "Gardner", "Kleindorfer"], "given-names": ["RY", "MG", "S"], "article-title": ["Isolation, characterisation and multiplex polymerase chain reaction of novel microsatellite loci for the avian parasite "], "italic": ["Philornis downsi "], "source": ["Mol Ecol Resources"], "year": ["2008"], "volume": ["8"], "fpage": ["142"], "lpage": ["144"], "pub-id": ["10.1111/j.1471-8286.2007.01900.x"]}, {"surname": ["Simon", "Frati", "Beckenbach", "Crespi", "Liu", "Flook"], "given-names": ["C", "F", "A", "B", "H", "P"], "article-title": ["Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers"], "source": ["Ann Entomol Soc Am"], "year": ["1994"], "volume": ["87"], "fpage": ["651"], "lpage": ["701"]}, {"surname": ["Queller", "Goodnight"], "given-names": ["DC", "KF"], "article-title": ["Estimating relatedness using genetic markers"], "source": ["Evolution"], "year": ["1989"], "volume": ["43"], "fpage": ["258"], "lpage": ["275"], "pub-id": ["10.2307/2409206"]}, {"surname": ["Raymond", "Rousset"], "given-names": ["M", "F"], "article-title": ["GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism"], "source": ["J Heredity"], "year": ["1995"], "volume": ["86"], "fpage": ["248"], "lpage": ["249"]}, {"surname": ["Hochberg"], "given-names": ["Y"], "article-title": ["A sharper Bonferroni procedure for multiple tests of significance"], "source": ["Biometrika"], "year": ["1988"], "volume": ["75"], "fpage": ["800"], "lpage": ["802"], "pub-id": ["10.1093/biomet/75.4.800"]}, {"surname": ["Vali\u00e8re"], "given-names": ["N"], "article-title": ["GIMLET: a computer program for analysing genetic individual identification data"], "source": ["Mol Ecol Notes"], "year": ["2002"], "volume": ["2"], "fpage": ["377"], "lpage": ["379"], "pub-id": ["10.1046/j.1471-8286.2002.00228.x"]}, {"surname": ["Weir", "Cockerham"], "given-names": ["BS", "CC"], "article-title": ["Estimating "], "italic": ["F"], "source": ["Evolution"], "year": ["1984"], "volume": ["38"], "fpage": ["1358"], "lpage": ["1370"], "pub-id": ["10.2307/2408641"]}, {"surname": ["Dieringer", "Schl\u00f6tterer"], "given-names": ["D", "C"], "article-title": ["Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets"], "source": ["Mol Ecol Notes"], "year": ["2003"], "volume": ["3"], "fpage": ["167"], "lpage": ["169"], "pub-id": ["10.1046/j.1471-8286.2003.00351.x"]}, {"surname": ["Peakall", "Smouse"], "given-names": ["R", "PE"], "article-title": ["Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research"], "source": ["Mol Ecol Notes"], "year": ["2006"], "volume": ["6"], "fpage": ["288"], "lpage": ["295"], "pub-id": ["10.1111/j.1471-8286.2005.01155.x"]}, {"surname": ["Luikart", "Cornuet"], "given-names": ["G", "JM"], "article-title": ["Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data"], "source": ["Conserv Biol"], "year": ["1998"], "volume": ["12"], "fpage": ["228"], "lpage": ["237"], "pub-id": ["10.1046/j.1523-1739.1998.96388.x"]}, {"surname": ["Piry", "Luikart", "Cornuet"], "given-names": ["S", "G", "JM"], "article-title": ["Bottleneck: a computer program for detecting recent reductions in the effective population size using allele frequency data"], "source": ["J Hered"], "year": ["1999"], "volume": ["90"], "fpage": ["502"], "lpage": ["503"], "pub-id": ["10.1093/jhered/90.4.502"]}, {"surname": ["Luikart", "Allendorf", "Cornuet", "Sherwin"], "given-names": ["G", "FW", "JM", "WB"], "article-title": ["Distortion of allele frequency distributions provides a test for recent population bottlenecks"], "source": ["J Heredity"], "year": ["1998"], "volume": ["89"], "fpage": ["238"], "lpage": ["247"], "pub-id": ["10.1093/jhered/89.3.238"]}, {"surname": ["Guillot", "Mortier", "Estoup"], "given-names": ["G", "F", "A"], "article-title": ["GENELAND: a compute package for landscape genetics"], "source": ["Mol Ecol Notes"], "year": ["2005"], "volume": ["5"], "fpage": ["712"], "lpage": ["715"], "pub-id": ["10.1111/j.1471-8286.2005.01031.x"]}, {"surname": ["Goudet"], "given-names": ["J"], "article-title": ["FSTAT, a program to estimate and test gene diversities and fixation indexes (version 2.3.9)"], "year": ["2001"]}, {"surname": ["Roque-Albelo", "Causton"], "given-names": ["L", "C"], "article-title": ["El Ni\u00f1o and introduced insects in the Gal\u00e1pagos Islands: different dispersal strategies, similar effects"], "source": ["Noticias de Gal\u00e1pagos"], "year": ["1999"], "volume": ["60"], "fpage": ["30"], "lpage": ["36"]}, {"surname": ["Silberglied"], "given-names": ["RE"], "article-title": ["Inter-island transport of insects aboard ships in the Gal\u00e1pagos Islands"], "source": ["Biol Conserv"], "year": ["1978"], "volume": ["13"], "fpage": ["273"], "lpage": ["275"], "pub-id": ["10.1016/0006-3207(78)90037-X"]}, {"surname": ["Wiedenfeld", "Jim\u00e9nez", "Fessl", "Kleindorfer", "Valarezo"], "given-names": ["DA", "UGA", "B", "S", "JC"], "article-title": ["Distribution of the introduced parasitic fly "], "italic": ["Philornis downsi "], "source": ["Pac Conserv Biol"], "year": ["2007"], "volume": ["13"], "fpage": ["14"], "lpage": ["19"]}, {"surname": ["Peck"], "given-names": ["SB"], "article-title": ["Aerial dispersal of insects between and to islands in the Gal\u00e1pagos Archipelago, Ecuador"], "source": ["Ann J Entomol Soc Am"], "year": ["1994"], "volume": ["87"], "fpage": ["218"], "lpage": ["224"]}, {"surname": ["Lambrinos"], "given-names": ["JG"], "article-title": ["How interactions between ecology and evolution influence contemporary invasion dynamics"], "source": ["Ecology"], "year": ["2004"], "volume": ["85"], "fpage": ["2061"], "lpage": ["2070"], "pub-id": ["10.1890/03-8013"]}, {"surname": ["Fessl", "Couri", "Tebbich"], "given-names": ["B", "MS", "S"], "italic": ["Philornis downsi "], "source": ["Stud Dipterol"], "year": ["2001"], "volume": ["8"], "fpage": ["317"], "lpage": ["322"]}, {"surname": ["Allendorf", "Lundquist"], "given-names": ["FW", "LL"], "article-title": ["Introduction: Population biology, evolution, and control of invasive species"], "source": ["Conserv Biol"], "year": ["2003"], "volume": ["17"], "fpage": ["24"], "lpage": ["30"], "pub-id": ["10.1046/j.1523-1739.2003.02365.x"]}, {"surname": ["Finston", "Peck"], "given-names": ["TL", "SB"], "article-title": ["Population structure and gene flow in Stomion: a species swarm of flightless beetles of the Gal\u00e1pagos Islands"], "source": ["Heredity"], "year": ["1995"], "volume": ["75"], "fpage": ["390"], "lpage": ["397"], "pub-id": ["10.1038/hdy.1995.151"]}, {"surname": ["Krafsur"], "given-names": ["ES"], "article-title": ["Sterile insect technique for suppressing and eradicating insect population: 55 years and counting"], "source": ["J Agric Entomol"], "year": ["1998"], "volume": ["15"], "fpage": ["303"], "lpage": ["317"]}]
BMC Ecol. 2008 Jul 31; 8:13