BACKGROUND
wolbachia pipientis  is a common endosymbiotic bacterium found in  <dig> to 50 % of insect, arachnid and terrestrial isopod species  <cit> .
in many host species wolbachia is a reproductive parasite, while other host species may benefit from infections  <cit> .
wolbachia is mostly maternally inherited, but occasional horizontal transmission into uninfected lineages contributes to the large number of infected species .
such switches of host lineages may result from hybridisation and introgression between closely related species  <cit> ; sharing of ecological niches  <cit> ; or transfer by or from parasitoids and predators  <cit> .
the initial contact between wolbachia of an infected with an uninfected host is the first essential step for the establishment of new wolbachia-host associations.
it must be followed by the colonisation of the host ovaries as a requirement for inheritance, as well as the induction of a reproductive manipulation or, alternatively, a positive fitness effect to promote its spread through new host populations  <cit> .
from an epidemiological perspective, this first exposure may be seen as a spillover of infectious agents from an infected source species to a new host species, and it may occur with or without its propagation in its new host species  <cit> .
for wolbachia, a model system for microbial symbiosis research, this initial spillover has so far not been differentiated from the more general concept of horizontal transfer itself.
most studies dealing with this phenomenon investigated horizontal transmission as the establishment of new heritable infections in new host lineages that may then express a reproductive phenotype leading to increases in infection frequencies in populations.
however, not every spillover will result in the establishment of an inherited infection.
spillovers may still be detectable as non-heritable somatic and transient wolbachia infections, or wolbachia dna  <cit>  e.g.
in digestive systems.
it can be expected that the frequency of spillover is larger than the actual establishment of newly inherited infections, and this difference may inflate global estimates of actual infection frequencies.
another overestimation of infection frequencies  may stem from lateral transfer of wolbachia genes into host genomes that can be seen as foot-prints of historic infections  <cit>  that may occur without the presence of current infections, yet this has so far rarely been addressed in wolbachia surveys  <cit> .
wolbachia spillover and with it wolbachia horizontal transmission may be driven by ecological interactions between infected and uninfected host species, and may therefore be specific to habitat and climate.
the diversity of ecological interactions is higher in tropical than temperate regions due to the increased biodiversity in the tropics  <cit> .
within this existing conceptual framework it can be argued that wolbachia spillover and horizontal transfer from infected to uninfected species should be more frequently observed in tropical than in temperate insect communities.
however, this is not reflected in studies comparing the overall incidence of wolbachia-infected species in arthropod communities from tropical and temperate regions that indicated parity of wolbachia incidence across climatic zones .
on a different scale, wolbachia prevalence within individual host species and populations was shown to vary from very low  <cit>  to fixation  <cit>  but little information is available on how this varying infection prevalence within species is distributed across extensive latitudinal clines and different climatic regions, except for the studies on the wolbachia prevalence in drosophila melanogaster and drosophila simulans populations along the east coast of australia .
we have previously investigated wolbachia horizontal transfer between host species and potential lateral gene transfer in australian tephritid species by analysing fruit fly communities of tropical australia  <cit> .
true fruit flies of the family tephritidae encompass approximately  <dig>  species worldwide, including key pests of the genera anastrepha, bactrocera, ceratitis, dacus and rhagoletis  <cit> .
the largest proportion of australian tephritids belongs to the genus bactrocera with  <dig> endemic species  <cit> .
the majority originate from and are restricted to the tropical regions  <cit> .
however, since establishment of horticultural production in australia in the 19th century, several fruit fly species have expanded into more temperate regions, in particular due to invasive expansion and planting of host plants  <cit> .
double infections of two wolbachia strains were found in individuals of six australian tephritid fruit fly species of the genera bactrocera and dacus  <cit> .
these two wolbachia strains were characterised with two distinct and complete multi locus sequence typing  profiles, st- <dig> and st- <dig>  individuals of a seventh fruit fly species also carried st- <dig> as a single infection.
two more wolbachia strains, st- <dig> and st- <dig>  also characterised with full sets of mlst genes, were detected in one individual of these bactrocera species, bactrocera frauenfeldi, and both of these strains were also found in individuals of a parasitoid wasp species, fopius arisanus , that is commonly found parasitising australian fruit fly larvae .
all infected tephritid species were sampled from tropical northern queensland where they have overlapping plant host range  <cit> , food sources  <cit>  and parasitoid complexes  <cit> , perhaps providing more opportunities for wolbachia spillover and horizontal transfer into new lineages within that environment than further south.
this earlier mlst based study also revealed the presence of pseudogenised mlst genes in two species, bactrocera peninsularis and bactrocera perkinsi  <cit> .
here we have expanded on our previous work by surveying more fruit fly individuals  collected over a significantly greater area of eastern australia, from tropical to temperate regions and over a distance of about  <dig>  km, with a focus on the two widespread and economically-relevant species bactrocera tryoni and bactrocera neohumeralis.
the main objective of our study was to assess incidence and prevalence of the four previously characterised wolbachia strains in australian tephritid fruit flies from geographically distant populations in order to test whether these high levels of shared wolbachia strains found in the tropics can also be seen in other regions of australia.
we hypothesised that, due to increased species diversity and interactions, the detection of wolbachia, including the detection of pseudogenes, may be higher in the more tropical regions of australia.
we also hypothesised that an overall low prevalence of identical wolbachia strains in individuals of several species in the tropics without proliferation in other regions may be a manifestation of a wolbachia spillover effect while the detection of wolbachia pseudogenes and failed amplification of other marker genes in all specimens of a species may suggest lateral transfer of wolbachia genes into the fly genome and the loss of infections in such species.
methods
insect samples
this study included field collected and laboratory specimens of fruit flies belonging to the genera bactrocera, dacus, dirioxa and ceratitis.
a total of  <dig> flies from  <dig> species were field-collected in new south wales, queensland, the northern territory and the torres strait islands during two sampling periods, from  <dig> to  <dig>  and from  <dig> to  <dig> .
while the majority of tephritid species reside in the equatorial and tropical regions of australia, ten of the  <dig> species included in this study also occur in the subtropical and temperate regions  <cit> .
most flies were male and collected in the summer by trapping with male attractants cue lure  <cit> , methyl eugenol for bactrocera visenda and zingerone for bactrocera jarvisi  <cit> .
both sexes of b. tryoni, bactrocera cacuminata and dirioxa pornia were collected directly on or from infested fruit.
fly specimens were identified under a stereomicroscope using identification keys  <cit> .
samples were selected based on availability and to canvas a range of species and populations.
prior to  <dig>  bactrocera aquilonis was recognised as a distinct species that is morphologically similar to b. tryoni but found only in the northern territory and western australia  <cit> .
hybridisation between these two species is considered to have occurred in the northern territory  <cit>  and it has been proposed that this species be synonymised with b. tryoni  <cit> ; therefore, here all flies originally classified as b. aquilonis collected from the northern territory have been listed as b. tryoni.
the species bactrocera papayae  has also recently been synonymised with bactrocera dorsalis  <cit> .
a subset of the fly samples  were selected for mlst analysis and reported in morrow et al.
<cit> .
in addition to field-caught flies, we also screened eight females from each of the following laboratory lines: b. tryoni, b. neohumeralis, b. jarvisi, and b. cacuminata kept at western sydney university, richmond, new south wales; two independent b. neohumeralis lines from cairns, queensland, and ceratitis capitata  from perth, western australia   <cit> .table  <dig>
wolbachia prevalence in australian fruit flies
bactrocera allwoodi
ba
bactrocera bryoniae
bb
bactrocera cacuminata
bca
bactrocera chorista
bch
bactrocera decurtans
bde
bactrocera dorsalis
bdo
bactrocera fallacis
bfa
bactrocera frauenfeldi
bfr
bactrocera jarvisi
bj
bactrocera manskii
bma
bactrocera murrayi
bmu
bactrocera neohumeralis
bn
bactrocera quadrata
bq
bactrocera strigifinis
bs
bactrocera tryoni
bt
bactrocera visenda
bv
dacus absonifacies
dab
dacus aequalis
dae
dacus axanus
dax
dacus bellulus
dbe
dacus newmani
dn
dirioxa pornia
dip
pseudogenesa
bactrocera peninsularis
bpen
bactrocera perkinsi
bper
a
bactrocera peninsularis amplified strongly for 16s rdna and fbpa but not for wsp, and the loci are therefore possibly pseudogenes.
bactrocera perkinsi amplified both 16s rdna and wsp, however co-amplified pseudogenes for wsp, coxa, hcpa and fbpa, but lacked ftsz  <cit> .
neither species was considered to harbour a true wolbachia infection for any analyses
dna extraction, pcr, cloning and sequencing
genomic dna was extracted from individual fruit fly abdomens, containing the reproductive organs, while the remainder of the specimen was stored in ethanol at −80 °c for subsequent independent confirmation of positive results.
prior to dna extraction, specimens were treated with 4 % sodium hypochlorite  for 5 min, then triton-x  and then thoroughly rinsed with milli-q water to reduce surface contamination.
dna from d. melanogaster line w <dig>  infected with wolbachia strain wmelpop  <cit> , was used as a positive control.
insect tissue was ground in  <dig>  ml microcentrifuge tubes with microtube pestles  and cell lysis performed overnight followed by extraction according to the genelute mammalian genomic dna miniprep kit  protocol.
elution of dna from spin columns was performed with 100 μl nuclease-free water.
risk of contamination was minimised by routinely replacing stock solutions, dispensing aliquots of stock reagents and using filter tips for all dna extractions and pcrs.
although cross contamination of flies caught in the same trap or stored in the same tube of ethanol after collection has previously been shown to be unlikely  <cit> , we have further minimised this risk by surface treatment of samples with sodium hypochlorite prior to dna extraction , plus independent extraction and pcr experiments in different laboratories.
pcr-based screening of fruit fly dna  was undertaken using the wolbachia surface protein  and 16s rdna loci.
primers for wsp were 81f and 691r  <cit>  and wsp-f and wsp-r  <cit> .
16s rdna was amplified with wspecf and wspecr  <cit> .
host mitochondrial cytochrome c oxidase subunit i  fragments were amplified with dick and pat  <cit>  as dna quality control.
veracity of pcr results was tested by inclusion of no-template controls.
all positive amplicons were confirmed by replication and further screening with other primer sets  <cit> .
coi and wsp amplicons were prepared for direct sequencing by treatment with a combination of  <dig> u exonuclease i  and  <dig> u shrimp alkaline phosphatase , with incubation at 37 °c for 30 min, then 95 °c for 5 min, prior to sequencing by macrogen .
for products displaying multiple sequences through direct sequencing, wsp was pcr amplified for cloning.
amplicons were either gel-extracted using the wizard sv gel and pcr clean-up system  and eluted in 25 μl nuclease-free water; or used directly in the ligation reaction.
ligation was with  <dig>  μl pgem-t easy vector , 1x rapid ligation buffer and 3u t <dig> dna ligase .
transformation of jm <dig> competent cells  was performed according to the manufacturer’s protocol.
colonies were smeared into a pcr tube using a sterile pipette tip and subjected to pcr using standard t <dig> promoter and sp <dig> primers with reaction and cycling conditions as described for insect coi .
positive clones, recognised by appropriately sized pcr products, were prepared for direct sequencing as described above.
a minimum of three clones, but usually eight clones for each transformed ligate were selected for sequencing in both directions, using t <dig> promoter and sp <dig> primers.
dig southern hybridisation
pcr amplicons were also authenticated by southern hybridisation using dig dna labelling and detection kit  based on the higher sensitivity method outlined in arthofer et al.
<cit> .
the wsp probe was generated as described in additional file  <dig>  using d. melanogaster w <dig> dna as template.
this method was applied to all individuals, to also confirm the absence of wsp amplicons.
flies were classified as uninfected when repeated amplification attempts with wsp and 16s rdna were negative but coi was positive.
individuals were considered wolbachia infected when wsp and 16s rdna primers amplified appropriately sized fragments from species for which we had established complete mlst profiles; wsp amplicons hybridised to the wsp probe; and wsp amplicons produced sequence homologues.
the coi locus of selected infected and uninfected flies was sequenced to determine the host-mitochondrial association of wolbachia sequences.
pcr-rflp
single restriction enzyme digestion was performed on wsp amplicons to test for multiple infections and to confirm the presence of no more than the two sequence types revealed via clone sequencing.
the sequence differences within two of the wsp alleles found as multiple infections enabled taqi  and spei  to distinguish the alleles.
taqi and spei  reactions were performed according to the manufacturer’s protocols for 3 h. the samples were electrophoresed on  <dig>  % agarose gels.
uncut, taqi cut and spei uncut bands were independently excised from the gel.
samples were purified using the wizard sv gel and pcr clean-up system and sequenced.
analyses of wolbachia incidence and prevalence
wolbachia incidence was defined as the percentage of species with individuals that were infected with a fully characterised wolbachia strain .
wolbachia prevalence was defined as the percentage of infected individuals within a species  <cit> .
furthermore, species for which all individuals had wolbachia pseudogenes and incomplete marker gene sets were considered as uninfected.
while we sampled an average of  <dig> individuals per species, sample size was limited for some species.
thus, we restricted inferences about wolbachia incidence to species for which we had at least ten individuals from within the same collection period.
this was to reduce the risk of underestimating wolbachia incidences in species where wolbachia occurs at low prevalence and was comparable to other recent studies of infection frequencies  <cit> .
for analysis of incidence, sampling locations were placed into five climate groups as defined by the australian bureau of meteorology, based on their shared köppen climate classifications: equatorial, tropical, subtropical, temperate and grasslands .
the latitudinal gradient was then split in half, with the midpoint set at 24°, south of gladstone, queensland.
wolbachia incidence between the northern and southern halves of the gradient was then tested through fisher’s exact test on the numbers of species with infected individuals versus species without infected individuals.fig.
<dig> incidence of wolbachia infection in tephritid species.
the köppen map of australia represents the climate groups and latitudes as approximate climate divisions along the east coast of australia.
wolbachia incidence is presented as filled sections of pie graphs, with the number of infected species per total number of tested species within each of the six regions.
grey pie graphs represent wolbachia incidence in all studied species of individual regions.
black pie graphs represent wolbachia incidence in species for which at least ten individuals were available north or south of 24° latitude .
wolbachia incidence was higher in the northern half than the southern half of the latitudinal gradient
prevalence was determined for twelve species only; species were selected based on either wolbachia infection or availability of a minimum of  <dig> sampled individuals if uninfected.
to incorporate the large differences in species distributions, a linear model  was fitted to the relationship of wolbachia prevalence with individual species distribution , based on the midpoint latitude of each species’ geographic range as per hancock et al.
<cit>  and royer and hancock  <cit> .
prevalence of wolbachia across geographic regions was only assessed for two widespread species, b. tryoni  and b. neohumeralis , as they were sampled from equatorial to temperate regions and had evidence of wolbachia infections.
the five köppen climate groups were further divided into  <dig> latitudinal zones based on common sampling locations .
infection prevalence across latitude and species was tested using multivariate generalised linear models  available in the package mvabund specifically designed for multivariate abundance data in community ecology  <cit>  and this was performed in r  <dig>   <cit> .table  <dig>
wolbachia prevalence across climate zones in two widely-distributed australian fruit fly species
b. neohumeralis
bn
b. tryoni
bt
the latitudinal gradient of eastern australia was divided into regions and represented by a major town or city.
darwin and alice springs are included as separate areas because they are not along the east coast of australia.
regions are classified into climate groups according to the köppen classification .
bactrocera neohumeralis does not occur in darwin, richmond or alice springs  <cit> , whereas non-determined  regions were not sampled
fisher’s exact test was also applied to test wolbachia prevalence over time within five species that were polymorphic for wolbachia infections and were collected across the two sampling periods.
phylogenetic analyses
phylogenetic analyses were performed for coi and wsp sequences in order to verify that mitochondrial lineages had independently acquired wolbachia strains through horizontal transfer .
it is also noted that phylogenetic analyses were not undertaken to infer species relationships, as single gene approaches are insufficient for such questions.
dna sequences were trimmed and edited in sequencher  <dig>   and then analysed in mega  <dig>   <cit> .
coi and wsp genes were independently aligned .
pairwise distance matrices were calculated for coi using number of differences and p-distance models.
substitution models were selected using find best dna model , which calculated the lowest bayesian information criterion score .
bayesian inference phylogenies were produced by mrbayes  <dig>   <cit>  running  <dig> generations with a sample frequency of  <dig>  the first 25 % of trees were discarded, and a 50 % majority rule consensus tree returned.
RESULTS
wolbachia-specific primers for wsp and 16s rdna were used to screen  <dig> field-collected australian fruit flies representing  <dig> species of bactrocera, dacus and dirioxa.
overall, individuals of eight species were positive for both wsp and 16s rdna, while all individuals of one species, b. peninsularis, were positive for 16s rdna without amplification of wsp .
all individuals from the independently established laboratory lines were negative for wolbachia .
initial screening using wsp primers 81f and 691r  <cit>  appeared to produce false positives for some flies.
thus, primers wsp-f and wsp-r  <cit>  were chosen, as their amplicons were more consistent and only occasionally produced spurious bands.
southern hybridisation with wsp was applied to all  <dig> samples, including those that initially appeared to be negative for wolbachia; this confirmed specificity as well as improved detection sensitivity.
in this way, four additional individuals, one each from four species  with very faint wsp pcr amplification and undetectable 16s rdna fragments using standard pcr were confirmed to carry wsp dna by southern hybridisation, while other individuals in these species were positive for both loci after standard pcr.
we considered these four individuals infected with low titre wolbachia because we had previously characterised complete mlst profiles in other individuals of these host species.
furthermore, control southern hybridisation to samples with high titre infections demonstrated that the dig-labelled wsp probe did not bind to primer-dimers or spurious wsp products, but hybridised to wsp amplicons that were verified by sequencing.
sequence analysis of wsp and 16s rdna
thirty-eight flies amplified at the wsp locus and  <dig> of these also successfully amplified at the 16s rdna locus.
sequence analysis identified five different wsp alleles of supergroup a. two wsp alleles  were found to co-occur in individuals of six species ; beyond this, wsp  <dig> occurred as a single sequence in one b. frauenfeldi, four b. tryoni and three b. neohumeralis individuals; and wsp  <dig> in one b. neohumeralis individual .
wolbachia mlst characterisation linked wsp  <dig> and wsp  <dig> with strains st- <dig> and st- <dig> respectively, demonstrating a high degree of sharing of wolbachia strains among six tephritid species  <cit> .
pcr-rflp and sequencing confirmed the absence of any other detectable wsp variants in these individuals.table  <dig>
wsp allele numbers and genbank accession numbers of analysed individuals
661
661
a
661
661
661
661
661
661
661
662
amplicons were sequenced and confirmed by pcr-rflp ; †wsp allele numbers from wolbachia wsp database, new wsp alleles in bold; a
wsp allele number not assigned
fig.
<dig> bayesian inference tree of wsp sequences from wolbachia within tephritid fruit flies and other invertebrates.
highlighted sequences are from this study .
all other sequences were retrieved from genbank or the wolbachia wsp database .
supergroup classification of wsp sequences is marked at the nodes.
posterior probabilities > <dig>  are marked at the nodes; tree was rooted by supergroup d and f strains.
scale bar represents the number of nucleotide substitutions per site
cloning of wsp from one b. frauenfeldi individual  revealed two alleles unlike those found in the other b. frauenfeldi individuals: one sequence identical to wsp  <dig> of drosophila simulans strain wri  <cit> ; the other  with high similarity to two wsp sequences detected in b. dorsalis from china  <cit> .
furthermore, the wsp and mlst alleles were shared with the wolbachia infections found in the parasitoid f. arisanus  <cit> .
all five b. perkinsi individuals produced wsp and 16s rdna amplicons.
cloning and sequencing of the wsp fragments revealed a novel allele, wsp  <dig>  as well as multiple clones of a wsp sequence with a single base insertion.
this insertion disrupts the open reading frame by incorporating stop codons, and may represent a wolbachia pseudogene transferred to the host genome of this species.
amplification of the ftsz locus in all individuals of b. perkinsi have previously failed during the mlst characterisation of wolbachia in this species while all other loci were positive with a group like sequences  <cit> .
we interpreted this as presence of multiple pseudogenes and absence of true wolbachia infections in the analysed specimens of this species.
for b. peninsularis, the 16s rdna fragment was consistently amplified from all seven individuals, whereas the wsp locus failed to amplify.
direct sequencing of the 16s rdna amplicon showed multiple peaks in the sequence chromatogram indicating two distinct sequences: one full length b group 16s rdna sequence, and one sequence with a 31 bp deletion at positions 182– <dig>  indicative of a potential pseudogene.
the sequence of the b. peninsularis full-length 16s rdna fragment was deposited in genbank, accession number kc <dig>  failure to detect a wsp sequence in any of these individuals, and prior failure to detect the five mlst markers except fbpa   <cit> , suggested that this sequence may not represent a genuine wolbachia infection; consequently this species was classified as uninfected.
as a control for the 16s rdna pcr assays, one b. neohumeralis sample, id <dig>  was chosen for 16s rdna sequencing to confirm homology to other wolbachia 16s rdna sequences in genbank .
analysis of wolbachia incidence and prevalence
geographically as well as irrespective of host species and wolbachia strains, wolbachia was detected in five of  <dig> samples from equatorial queensland ,  <dig> of  <dig>  samples from tropical queensland and the northern territory, and  <dig> of  <dig>  samples from subtropical queensland.
however, none of the  <dig> samples from temperate new south wales and none of the  <dig> samples from central australia had detectable wolbachia .
the analysis of wolbachia incidence within species revealed that it was restricted to the northern half of australia ; analysis of species for which at least ten individuals per species were available for each half revealed wolbachia incidence in at least one specimen for a number of these species in the northern half of the gradient while none was infected in the southern half.
therefore, incidence of wolbachia was higher in the northern half of the gradient, with wolbachia incidence in seven out of nine species for which a minimum of ten individuals were tested; all four species from the southern half for which a minimum of ten individuals per species were available in that region, were negative .
a linear model was applied in order to appropriately represent both widespread and more tropically restricted species and this detected a subtle relationship between the prevalence of wolbachia and the midpoint latitude of twelve individual species .
for the seven species that had individuals with a full wolbachia mlst profile plus wsp, the prevalence ranged from  <dig>  %  to  <dig>  % .
multivariate analysis of infection prevalence in the only two widespread species that are endemic from equatorial to temperate climes, b. neohumeralis and b. tryoni , detected a significant interaction of wolbachia prevalence by latitude  with wolbachia being more common in the northern regions of east coast australia.
univariate analysis with adjusted p-values  to account for family-wise error across both species did not detect a significant interaction , demonstrating that the detectable effect on a single species is weak, but the pattern is strengthened when analysing both species together.table  <dig> analyses of wolbachia prevalence at different latitudes of b. neohumeralis and b. tryoni
b. neohumeralis
b. tryoni
table presents analysis of deviance for multivariate and univariate analyses of wolbachia prevalence  at different latitudes of two widespread species b. neohumeralis and b. tryoni, as per output of the mvabund package.
p-value ) significance codes:  <dig> ‘***’  <dig>  ‘**’  <dig>  ‘*’  <dig>  ‘.’  <dig>  ‘’ 1
furthermore, temporal effects on wolbachia prevalence within b. bryoniae, b. frauenfeldi, b. neohumeralis, b. strigifinis and b. tryoni were tested for equatorial and tropical samples by using a fisher’s exact test  over the sampling years.
no significant temporal change within the tested regions was detected, except for b. tryoni , but this represented neither an overall increase nor decrease in prevalence over time .
analysis of mitochondrial dna
most tephritid species produced clear coi sequences.
bactrocera jarvisi, b. murrayi, b. dorsalis and b. perkinsi individuals produced ambiguous sequences, indicative of potential nuclear mitochondrial  dna in these species.
the latter sequences were not included in phylogenetic analysis.
instead, b. jarvisi and b. dorsalis sequences were retrieved from genbank.
bayesian analysis of  <dig> sequences  over 571 bp returned a well-supported consensus tree .fig.
<dig> bayesian inference tree of coi sequences from tephritid fruit flies.
representatives of the five major tephritid fruit fly genera, anastrepha, bactrocera, ceratitis, dacus and rhagoletis as well as dirioxa are shown.
dark grey shading highlights the individuals from this study harbouring wolbachia, light grey highlights specimens that were negative for wolbachia.
sequences from species without shading were retrieved from genbank .
tree was rooted with drosophila melanogaster, and nodes are labelled with posterior probability values.
scale bar represents the number of nucleotide substitutions per site
mitochondrial haplotypes of morphologically identified individuals of different species harbouring identical wsp alleles were different, except for the two sibling species b. neohumeralis and b. tryoni; both had individuals with identical wolbachia and mitochondrial sequences but without any linkage between wolbachia infection and haplotypes.
coi distance measures of b. neohumeralis and b. tryoni  showed little difference within infected  and uninfected individuals , and between infected and uninfected individuals .
discussion
we surveyed  <dig> australian tephritid fruit fly species throughout their natural range, and detected wolbachia sequences in individuals of nine species.
in two species, all individuals had pseudogenes and incomplete sets of marker genes – indicative of potential lateral gene transfer to the host genome without the presence of an existing infection.
both species are restricted to the tropics.
the remaining seven fruit fly species had both 16s rdna and wsp sequences characteristic for four wolbachia strains for which complete mlst profiles had previously been established  <cit> .
the four wolbachia strains occurred as paired or single infections in more than one host species, indicative of horizontal wolbachia transmission events into these species.
overall the incidence of wolbachia strains  was restricted to tropical regions, and the prevalence of the four wolbachia strains was very low in infected species .
in five infected species , wolbachia prevalence appeared to not have changed over time.
in b. neohumeralis and b. tryoni, wolbachia was not linked with any particular mitochondrial haplotypes.
we thus concluded that the detected strains may not necessarily manipulate reproduction of their host species and maternal inheritance may also be unreliable.
therefore, the detection of identical wolbachia strains at low prevalence in several species of this fruit fly community may be due to a spillover from a yet unknown source species, without the propagation of wolbachia in any of the analysed host species, and this spillover appears to be restricted to the tropical regions.
wolbachia strains and pseudogenes
the previous mlst characterisation in australian tropical tephritids revealed four complete mlst profiles of supergroup a wolbachia in seven species  <cit> .
here, however, 16s rdna and wsp sequences were isolated from nine and eight species, respectively while complete marker gene sets were not found for all nine species, which is likely to be due to lateral gene transfer of wolbachia into the host genome.
more specifically, all individuals of b. perkinsi produced 16s rdna and wsp sequences .
in the previous study, no individuals of this species amplified ftsz and produced both pseudogenic and novel orf sequences for other mlst loci  <cit> .
another species, b. peninsularis, also had signs of wolbachia lateral gene transfer, as all individuals only amplified 16s rdna but not wsp in this study while previously fbpa was detected as the only positive mlst locus in all individuals of this species  <cit> .
similarly, pseudogenised 16s rdna and fbpa sequences were found in one out of eight tested glossina fly species  <cit> .
we thus conclude that the wolbachia signals in b. peninsularis and b. perkinsi are due to lateral gene transfer from wolbachia into these host genomes, as has been reported for other insects  <cit> .
it may also be noted that the wolbachia pseudogenic sequences were different enough to conclude that they did not derive from the four strains detected in this study.
two of the previously characterised wolbachia strains st- <dig> and st- <dig>  <cit>  were found with their 16s rdna and wsp sequences in individuals of six fruit fly species.
single infections by both of these strains were detected in b. neohumeralis individuals and single infections by st- <dig> in all four infected b. tryoni individuals and one b. frauenfeldi individual.
two additional strains, st- <dig> and st- <dig>  previously characterised by morrow et al.
<cit>  were detected with their 16s rdna and wsp sequence variants in another b. frauenfeldi individual.
it was demonstrated that this double infection is identical to wolbachia from the fruit fly parasitoid f. arisanus, and this could be a spillover from this parasitoid species.
this conclusion can be drawn based on the result of our earlier study in which all individuals of this parasitoid species were infected, including individuals from the sydney region where fruit fly host populations were uninfected  <cit> .
of additional interest is that both f. arisanus and b. frauenfeldi are not endemic to australia; f. arisanus was introduced to australia in  <dig> and now occurs from northern queensland to sydney  while b. frauenfeldi was first detected on cape york peninsula in  <dig> from where it has spread throughout northern queensland  <cit> .
this also means that st- <dig> and st- <dig> were only found in two host species that were originally not native to australia.
latitudinal distribution of wolbachia in australian tephritids
after excluding pseudogenes in two species as potential wolbachia infections, we detected wolbachia in 29 % of tested australian tephritid species, and this is lower than previous estimates of 40 % incidence of infections in arthropods, mostly based on a european dataset  <cit> .
our results varied between climatic regions.
wolbachia incidence was restricted to equatorial and tropical regions north of gladstone, while wolbachia was absent from individuals caught in the more southern subtropical, temperate and arid inland regions.
specifically, wolbachia was found in species whose range is limited to tropical and equatorial queensland ; in species with a broader climatic range, wolbachia occurred as far south as cairns for b. bryoniae, mackay for b. tryoni and gladstone for b. neohumeralis, but was absent from individuals in more southerly subtropical and temperate regions.
within all of these species, wolbachia prevalence was very low , and in the more widely distributed species b. neohumeralis and b. tryoni, wolbachia was also restricted to the northern ranges of each species.
wolbachia in australian tephritids: transient infections resulting from a spillover?
previously, identical wolbachia strains were found in multiple species of fruit flies that share ecological niches, but are unlikely to hybridise in nature, except possibly for b. tryoni and b. neohumeralis  <cit> .
this study found that only low frequency wolbachia signals were detected in samples from tropical environments, including for those species that are widespread from tropical to temperate climes.
furthermore, no wolbachia was detected in any laboratory colonies of these species, including those originating from tropical locations, a result that was supported by  <dig> pyrosequencing of 16s rdna amplicon libraries from these same individuals  <cit> .
this suggests that there is ample opportunity for the horizontal transfer of wolbachia, but the low prevalence, and lack of increasing prevalence over time, does not unequivocally show that these infections confer any selective advantage to infected females e.g.
through reproductive manipulation such as cytoplasmic incompatibility  nor that they are being reliably inherited.
experimental evidence of a heritable wolbachia infection is necessary to confirm an established infection.
in drosophila suzukii, a similarly low frequency wolbachia infection occurs with variable maternal transmission; however no ci or fitness benefits  were observed, and maintenance of wolbachia in this species has yet to be explained  <cit> .
here, in the absence of any evidence of vertical inheritance, the finding of two strains  in six host species at low prevalence is difficult to explain, and perhaps more likely due to a spillover of wolbachia from an as yet unknown host species without any proliferation in any of the studied species.
this interpretation still lacks direct experimental evidence, but transient infections are frequently observed in wolbachia microinjection studies that detect wolbachia in many individuals emerging from microinjected embryos.
however, only a few injected embryos result in the establishment of stably inherited matrilines while wolbachia is lost from most other matrilines  <cit> , or remains undetectable in some matrilines for many generations  <cit>  .
in our study, the potential spillover of wolbachia is restricted to the tropical regions, and is perhaps due to the higher biodiversity, resulting in increased ecological interactions in the tropics.
the finding of two strains  that are shared between fruit fly and parasitoid individuals supports the interpretation of a potential spillover with the parasitoid being the source in this instance; st- <dig> and st- <dig> were only found in one b. frauenfeldi individual while all f. arisanus individuals were infected with both strains, including from regions where host fruit flies were uninfected  <cit> .
if, however, the infections are heritable in the studied fruit fly species, then the latitudinal restriction of wolbachia incidence and prevalence in australian fruit flies to northern australia may be due to  wolbachia loss due to environmental conditions of their hosts in the southern range or  an ongoing wolbachia invasion of fruit fly species from north to south, driven by a reproductive phenotype.
an increasing number of studies have demonstrated that wolbachia responds to the climatic environment of its host insects.
for example, precipitation frequency appeared to correlate with the distribution of single and double infections in leaf beetle individuals in panama, with multiple infections restricted to wetter regions  <cit> .
in another study, multiple infections were found to be more frequent in one tropical habitat when compared with two temperate habitats  <cit> .
temperature is known to affect wolbachia titres in a number of insect hosts.
for example, the virulent wmelpop strain appears to over-replicate at elevated temperatures  <cit> .
however, high temperatures may also effectively reduce wolbachia densities in hosts through increased wolbachia-bacteriophage activity  <cit> , and thus result in diminished penetrance of ci or male-killing phenotypes  <cit> .
it is possible that temperature fluctuations and extremes, commonly found in southern regions of australia, may be selecting more strongly against wolbachia than the more constant temperature conditions in the tropics.
such climatic effects could also be contributing to the distribution of wolbachia in previously reported analyses such as d. melanogaster  <cit>  and cat fleas  <cit> .
similarly, climate and latitude were found to determine ranges of other microbial symbioses, for example in marine invertebrates  <cit> , terrestrial insects  and humans  <cit> .
the second scenario supposes that the wolbachia in the australian tephritid species, in particular infections of the widespread b. tryoni and b. neohumeralis, is relatively recently acquired, and a progressive ci driven sweep may therefore be in its infancy.
however this hypothesis is less likely for the following reasons; a ci driven invasion should result in an increase of infection prevalence over time that was not detected in these species.
furthermore, wolbachia-induced ci is expected to cause a selective sweep of infected mitochondrial haplotypes  <cit> , yet there does not appear to be support for this in the b. tryoni and b. neohumeralis species complex  <cit>  because both species share mitochondrial haplotypes across infected and uninfected individuals.
we were not able to directly assess the ci phenotype of wolbachia in australian tephritids due to the lack of infected laboratory populations.
the ci characterisation of the wolbachia strains will thus require future field collection efforts, in particular of females, in order to set up infected laboratory colonies for crossing experiments.
testing of field females will also avoid a potential underestimate of wolbachia infection rates due to our male-biased sampling approach that would not detect male-killing wolbachia strains.
however, prevalence of male-killing wolbachia may be generally low as found for another fruit fly family, drosophilidae  <cit> .
evidence for numtdna in bactrocera
blacket et al.
<cit>  have previously detected numtdna amplicons for b. tryoni and b. neohumeralis by using primers for a different section of coi; our different primer set did not yield numtdna for b. tryoni and b. neohumeralis but revealed potential numtdna sequences in other species such as b. jarvisi, b. dorsalis, b. murrayi and b. perkinsi.
interestingly, potential transfer of both mitochondrial and wolbachia genes into the nuclear host genome were detected for b. perkinsi.
CONCLUSIONS
in this study of australian tephritids we detected several wolbachia pseudogenes that may be host-genomic footprints of historic infections that have since been lost from host species.
within the species of relevance , these pseudogenes were fixed at 100 % while prevalence of wolbachia infections in other species was overly low .
under the assumption that the detected wolbachia strains in seven fruit fly species may not be reliably inherited nor cause a reproductive phenotype, these infections could be considered transient and a result of a spillover from a yet unknown source .
in our study system, this spillover appears to be restricted to northern australia and has not resulted in further proliferation within any of the tested host lineages, and this warrants future investigation.
the fact that we have not found a single case of a high prevalence infected australian tephritid fruit fly species in  <dig> tested species raises the question whether these taxa exhibit traits that reduce wolbachia infectivity while two species had wolbachia infections in their evolutionary past.
other related pest taxa such as olive fly bactrocera oleae appear uninfected in natural populations while laboratory populations were successfully infected  <cit> ; a high incidence of wsp genes was found in thai bactrocera species  <cit>  while prevalence of wsp in b. dorsalis was found to be extremely low  <cit> , although these strains were not further characterised.
perhaps some of these bactrocera species also contain wolbachia signals that are either pseudogenes and/or transient somatic infections, and similar questions may be valid for wolbachia signals found in other host species.
our study highlights that tropical insect communities have signs of lateral gene transfer and may be exposed to increased horizontal wolbachia transfer.
it also emphasises that global estimates for wolbachia infection frequencies should be interpreted with caution, and account for lateral gene transfer, variable prevalence within species and transient spillover effects.
this may be of particular relevance if wolbachia surveys are based on a single or a few marker genes for detection of infections in a host species, and without a broader sampling effort from the host’s distribution .
availability of supporting data
data sets supporting the results of this article are included within the article and its additional files.
dna sequences generated in this study were submitted to genbank: accession numbers kc581371–kc <dig>  kc <dig>  kc775785-kc <dig>  and kc668319-kc <dig>  kc668323-kc <dig>  kc <dig>  and kc775793-kc <dig> .
new wsp alleles have also been deposited in the wolbachia wsp database .
additional files
additional file 1:
collection data for field-caught fruit flies from this study.
information for collection location, climate group and year, specimen specific id number, sex, result for wsp and 16s rdna screening and genbank accession number of coi partial gene are included.
additional file 2:
laboratory strains screened for
wolbachia
using
wsp
and 16s rdna.
additional file 3:
pcr conditions and primer sequences used in this research.
additional file 4:
relationship between prevalence of
wolbachia
infection in twelve tephritid species and the midpoint of their latitudinal distribution, based on hancock
et al.
[
32
] and royer and hancock [
33
] .
at least ten individuals were tested for uninfected species .
additional file 5:
comparison of
wolbachia
infected individuals sampled from equatorial and tropical regions for
b. bryoniae
,
b. frauenfeldi
,
b. neohumeralis,
,
b. strigifinis
and
b. tryoni
over five collection years.
fisher’s exact test shows overall no significant  temporal effect at α =  <dig> , except for b. tryoni, but this did not represent an overall increase or decrease in prevalence.
additional file 6:
distance matrix of
coi
sequences covering 551 bp positions of
b. tryoni
and
b. neohumeralis
individuals, grouped by
wolbachia
-infected  and uninfected .
below the diagonal are the pairwise number of differences and above the diagonal are p-distance values.
p-distances given for within- and between- group means are listed below the main table.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
the project was conceived by mr and jlm; the research was designed by jlm and mr with input from mf and dcas; jlm performed the research and analysed the data; mr, mf, jer and dcas supplied materials and specimens; jlm and mr wrote the paper with input by mf, jer and dcas.
all authors read and approved the final manuscript.

 wolbachia pipientis is a common endosymbiotic bacterium found in to 50 % of insect, arachnid and terrestrial isopod species. in many host species wolbachia is a reproductive parasite, while other host species may benefit from infections. wolbachia is mostly maternally inherited, but occasional horizontal transmission into uninfected lineages contributes to the large number of infected species. such switches of host lineages may result from hybridisation and introgression between closely related species; sharing of ecological niches; or transfer by or from parasitoids and predators. the initial contact between wolbachia of an infected with an uninfected host is the first essential step for the establishment of new wolbachia-host associations. it must be followed by the colonisation of the host ovaries as a requirement for inheritance, as well as the induction of a reproductive manipulation or, alternatively, a positive fitness effect to promote its spread through new host populations. from an epidemiological perspective, this first exposure may be seen as a spillover of infectious agents from an infected source species to a new host species, and it may occur with or without its propagation in its new host species. for wolbachia, a model system for microbial symbiosis research, this initial spillover has so far not been differentiated from the more general concept of horizontal transfer itself. most studies dealing with this phenomenon investigated horizontal transmission as the establishment of new heritable infections in new host lineages that may then express a reproductive phenotype leading to increases in infection frequencies in populations. however, not every spillover will result in the establishment of an inherited infection. spillovers may still be detectable as non-heritable somatic and transient wolbachia infections, or wolbachia dna e.g. in digestive systems. it can be expected that the frequency of spillover is larger than the actual establishment of newly inherited infections, and this difference may inflate global estimates of actual infection frequencies. another overestimation of infection frequencies  may stem from lateral transfer of wolbachia genes into host genomes that can be seen as foot-prints of historic infections that may occur without the presence of current infections, yet this has so far rarely been addressed in wolbachia surveys. wolbachia spillover and with it wolbachia horizontal transmission may be driven by ecological interactions between infected and uninfected host species, and may therefore be specific to habitat and climate. the diversity of ecological interactions is higher in tropical than temperate regions due to the increased biodiversity in the tropics. within this existing conceptual framework it can be argued that wolbachia spillover and horizontal transfer from infected to uninfected species should be more frequently observed in tropical than in temperate insect communities. however, this is not reflected in studies comparing the overall incidence of wolbachia-infected species in arthropod communities from tropical and temperate regions that indicated parity of wolbachia incidence across climatic zones. on a different scale, wolbachia prevalence within individual host species and populations was shown to vary from very low to fixation but little information is available on how this varying infection prevalence within species is distributed across extensive latitudinal clines and different climatic regions, except for the studies on the wolbachia prevalence in drosophila melanogaster and drosophila simulans populations along the east coast of australia. we have previously investigated wolbachia horizontal transfer between host species and potential lateral gene transfer in australian tephritid species by analysing fruit fly communities of tropical australia. true fruit flies of the family tephritidae encompass approximately species worldwide, including key pests of the genera anastrepha, bactrocera, ceratitis, dacus and rhagoletis. the largest proportion of australian tephritids belongs to the genus bactrocera with endemic species. the majority originate from and are restricted to the tropical regions. however, since establishment of horticultural production in australia in the 19th century, several fruit fly species have expanded into more temperate regions, in particular due to invasive expansion and planting of host plants. double infections of two wolbachia strains were found in individuals of six australian tephritid fruit fly species of the genera bactrocera and dacus. these two wolbachia strains were characterised with two distinct and complete multi locus sequence typing profiles, st- and st- individuals of a seventh fruit fly species also carried st- as a single infection. two more wolbachia strains, st- and st- also characterised with full sets of mlst genes, were detected in one individual of these bactrocera species, bactrocera frauenfeldi, and both of these strains were also found in individuals of a parasitoid wasp species, fopius arisanus, that is commonly found parasitising australian fruit fly larvae. all infected tephritid species were sampled from tropical northern queensland where they have overlapping plant host range, food sources and parasitoid complexes, perhaps providing more opportunities for wolbachia spillover and horizontal transfer into new lineages within that environment than further south. this earlier mlst based study also revealed the presence of pseudogenised mlst genes in two species, bactrocera peninsularis and bactrocera perkinsi. here we have expanded on our previous work by surveying more fruit fly individuals collected over a significantly greater area of eastern australia, from tropical to temperate regions and over a distance of about  km, with a focus on the two widespread and economically-relevant species bactrocera tryoni and bactrocera neohumeralis. the main objective of our study was to assess incidence and prevalence of the four previously characterised wolbachia strains in australian tephritid fruit flies from geographically distant populations in order to test whether these high levels of shared wolbachia strains found in the tropics can also be seen in other regions of australia. we hypothesised that, due to increased species diversity and interactions, the detection of wolbachia, including the detection of pseudogenes, may be higher in the more tropical regions of australia. we also hypothesised that an overall low prevalence of identical wolbachia strains in individuals of several species in the tropics without proliferation in other regions may be a manifestation of a wolbachia spillover effect while the detection of wolbachia pseudogenes and failed amplification of other marker genes in all specimens of a species may suggest lateral transfer of wolbachia genes into the fly genome and the loss of infections in such species. methods insect samples this study included field collected and laboratory specimens of fruit flies belonging to the genera bactrocera, dacus, dirioxa and ceratitis. a total of flies from species were field-collected in new south wales, queensland, the northern territory and the torres strait islands during two sampling periods, from to and from to. while the majority of tephritid species reside in the equatorial and tropical regions of australia, ten of the species included in this study also occur in the subtropical and temperate regions. most flies were male and collected in the summer by trapping with male attractants cue lure, methyl eugenol for bactrocera visenda and zingerone for bactrocera jarvisi. both sexes of b. tryoni, bactrocera cacuminata and dirioxa pornia were collected directly on or from infested fruit. fly specimens were identified under a stereomicroscope using identification keys. samples were selected based on availability and to canvas a range of species and populations. prior to bactrocera aquilonis was recognised as a distinct species that is morphologically similar to b. tryoni but found only in the northern territory and western australia. hybridisation between these two species is considered to have occurred in the northern territory and it has been proposed that this species be synonymised with b. tryoni; therefore, here all flies originally classified as b. aquilonis collected from the northern territory have been listed as b. tryoni. the species bactrocera papayae has also recently been synonymised with bactrocera dorsalis. a subset of the fly samples were selected for mlst analysis and reported in morrow et al.. in addition to field-caught flies, we also screened eight females from each of the following laboratory lines: b. tryoni, b. neohumeralis, b. jarvisi, and b. cacuminata kept at western sydney university, richmond, new south wales; two independent b. neohumeralis lines from cairns, queensland, and ceratitis capitata from perth, western australia.table wolbachia prevalence in australian fruit flies bactrocera allwoodi ba bactrocera bryoniae bb bactrocera cacuminata bca bactrocera chorista bch bactrocera decurtans bde bactrocera dorsalis bdo bactrocera fallacis bfa bactrocera frauenfeldi bfr bactrocera jarvisi bj bactrocera manskii bma bactrocera murrayi bmu bactrocera neohumeralis bn bactrocera quadrata bq bactrocera strigifinis bs bactrocera tryoni bt bactrocera visenda bv dacus absonifacies dab dacus aequalis dae dacus axanus dax dacus bellulus dbe dacus newmani dn dirioxa pornia dip pseudogenesa bactrocera peninsularis bpen bactrocera perkinsi bper a bactrocera peninsularis amplified strongly for 16s rdna and fbpa but not for wsp, and the loci are therefore possibly pseudogenes. bactrocera perkinsi amplified both 16s rdna and wsp, however co-amplified pseudogenes for wsp, coxa, hcpa and fbpa, but lacked ftsz. neither species was considered to harbour a true wolbachia infection for any analyses dna extraction, pcr, cloning and sequencing genomic dna was extracted from individual fruit fly abdomens, containing the reproductive organs, while the remainder of the specimen was stored in ethanol at −80 °c for subsequent independent confirmation of positive results. prior to dna extraction, specimens were treated with 4 % sodium hypochlorite for 5 min, then triton-x and then thoroughly rinsed with milli-q water to reduce surface contamination. dna from d. melanogaster line w infected with wolbachia strain wmelpop, was used as a positive control. insect tissue was ground in  ml microcentrifuge tubes with microtube pestles and cell lysis performed overnight followed by extraction according to the genelute mammalian genomic dna miniprep kit protocol. elution of dna from spin columns was performed with 100 μl nuclease-free water. risk of contamination was minimised by routinely replacing stock solutions, dispensing aliquots of stock reagents and using filter tips for all dna extractions and pcrs. although cross contamination of flies caught in the same trap or stored in the same tube of ethanol after collection has previously been shown to be unlikely, we have further minimised this risk by surface treatment of samples with sodium hypochlorite prior to dna extraction, plus independent extraction and pcr experiments in different laboratories. pcr-based screening of fruit fly dna was undertaken using the wolbachia surface protein and 16s rdna loci. primers for wsp were 81f and 691r and wsp-f and wsp-r. 16s rdna was amplified with wspecf and wspecr. host mitochondrial cytochrome c oxidase subunit i fragments were amplified with dick and pat as dna quality control. veracity of pcr results was tested by inclusion of no-template controls. all positive amplicons were confirmed by replication and further screening with other primer sets. coi and wsp amplicons were prepared for direct sequencing by treatment with a combination of u exonuclease i and u shrimp alkaline phosphatase, with incubation at 37 °c for 30 min, then 95 °c for 5 min, prior to sequencing by macrogen. for products displaying multiple sequences through direct sequencing, wsp was pcr amplified for cloning. amplicons were either gel-extracted using the wizard sv gel and pcr clean-up system and eluted in 25 μl nuclease-free water; or used directly in the ligation reaction. ligation was with  μl pgem-t easy vector, 1x rapid ligation buffer and 3u t dna ligase. transformation of jm competent cells was performed according to the manufacturer’s protocol. colonies were smeared into a pcr tube using a sterile pipette tip and subjected to pcr using standard t promoter and sp primers with reaction and cycling conditions as described for insect coi. positive clones, recognised by appropriately sized pcr products, were prepared for direct sequencing as described above. a minimum of three clones, but usually eight clones for each transformed ligate were selected for sequencing in both directions, using t promoter and sp primers. dig southern hybridisation pcr amplicons were also authenticated by southern hybridisation using dig dna labelling and detection kit based on the higher sensitivity method outlined in arthofer et al.. the wsp probe was generated as described in additional file using d. melanogaster w dna as template. this method was applied to all individuals, to also confirm the absence of wsp amplicons. flies were classified as uninfected when repeated amplification attempts with wsp and 16s rdna were negative but coi was positive. individuals were considered wolbachia infected when wsp and 16s rdna primers amplified appropriately sized fragments from species for which we had established complete mlst profiles; wsp amplicons hybridised to the wsp probe; and wsp amplicons produced sequence homologues. the coi locus of selected infected and uninfected flies was sequenced to determine the host-mitochondrial association of wolbachia sequences. pcr-rflp single restriction enzyme digestion was performed on wsp amplicons to test for multiple infections and to confirm the presence of no more than the two sequence types revealed via clone sequencing. the sequence differences within two of the wsp alleles found as multiple infections enabled taqi and spei to distinguish the alleles. taqi and spei reactions were performed according to the manufacturer’s protocols for 3 h. the samples were electrophoresed on  % agarose gels. uncut, taqi cut and spei uncut bands were independently excised from the gel. samples were purified using the wizard sv gel and pcr clean-up system and sequenced. analyses of wolbachia incidence and prevalence wolbachia incidence was defined as the percentage of species with individuals that were infected with a fully characterised wolbachia strain. wolbachia prevalence was defined as the percentage of infected individuals within a species. furthermore, species for which all individuals had wolbachia pseudogenes and incomplete marker gene sets were considered as uninfected. while we sampled an average of individuals per species, sample size was limited for some species. thus, we restricted inferences about wolbachia incidence to species for which we had at least ten individuals from within the same collection period. this was to reduce the risk of underestimating wolbachia incidences in species where wolbachia occurs at low prevalence and was comparable to other recent studies of infection frequencies. for analysis of incidence, sampling locations were placed into five climate groups as defined by the australian bureau of meteorology, based on their shared köppen climate classifications: equatorial, tropical, subtropical, temperate and grasslands. the latitudinal gradient was then split in half, with the midpoint set at 24°, south of gladstone, queensland. wolbachia incidence between the northern and southern halves of the gradient was then tested through fisher’s exact test on the numbers of species with infected individuals versus species without infected individuals.fig. incidence of wolbachia infection in tephritid species. the köppen map of australia represents the climate groups and latitudes as approximate climate divisions along the east coast of australia. wolbachia incidence is presented as filled sections of pie graphs, with the number of infected species per total number of tested species within each of the six regions. grey pie graphs represent wolbachia incidence in all studied species of individual regions. black pie graphs represent wolbachia incidence in species for which at least ten individuals were available north or south of 24° latitude. wolbachia incidence was higher in the northern half than the southern half of the latitudinal gradient prevalence was determined for twelve species only; species were selected based on either wolbachia infection or availability of a minimum of sampled individuals if uninfected. to incorporate the large differences in species distributions, a linear model was fitted to the relationship of wolbachia prevalence with individual species distribution, based on the midpoint latitude of each species’ geographic range as per hancock et al. and royer and hancock. prevalence of wolbachia across geographic regions was only assessed for two widespread species, b. tryoni and b. neohumeralis, as they were sampled from equatorial to temperate regions and had evidence of wolbachia infections. the five köppen climate groups were further divided into latitudinal zones based on common sampling locations. infection prevalence across latitude and species was tested using multivariate generalised linear models available in the package mvabund specifically designed for multivariate abundance data in community ecology and this was performed in r.table wolbachia prevalence across climate zones in two widely-distributed australian fruit fly species   b. neohumeralis bn  b. tryoni bt the latitudinal gradient of eastern australia was divided into regions and represented by a major town or city. darwin and alice springs are included as separate areas because they are not along the east coast of australia. regions are classified into climate groups according to the köppen classification. bactrocera neohumeralis does not occur in darwin, richmond or alice springs, whereas non-determined regions were not sampled fisher’s exact test was also applied to test wolbachia prevalence over time within five species that were polymorphic for wolbachia infections and were collected across the two sampling periods. phylogenetic analyses phylogenetic analyses were performed for coi and wsp sequences in order to verify that mitochondrial lineages had independently acquired wolbachia strains through horizontal transfer. it is also noted that phylogenetic analyses were not undertaken to infer species relationships, as single gene approaches are insufficient for such questions. dna sequences were trimmed and edited in sequencher and then analysed in mega. coi and wsp genes were independently aligned. pairwise distance matrices were calculated for coi using number of differences and p-distance models. substitution models were selected using find best dna model, which calculated the lowest bayesian information criterion score. bayesian inference phylogenies were produced by mrbayes running generations with a sample frequency of the first 25 % of trees were discarded, and a 50 % majority rule consensus tree returned.  wolbachia-specific primers for wsp and 16s rdna were used to screen field-collected australian fruit flies representing species of bactrocera, dacus and dirioxa. overall, individuals of eight species were positive for both wsp and 16s rdna, while all individuals of one species, b. peninsularis, were positive for 16s rdna without amplification of wsp. all individuals from the independently established laboratory lines were negative for wolbachia. initial screening using wsp primers 81f and 691r appeared to produce false positives for some flies. thus, primers wsp-f and wsp-r were chosen, as their amplicons were more consistent and only occasionally produced spurious bands. southern hybridisation with wsp was applied to all samples, including those that initially appeared to be negative for wolbachia; this confirmed specificity as well as improved detection sensitivity. in this way, four additional individuals, one each from four species with very faint wsp pcr amplification and undetectable 16s rdna fragments using standard pcr were confirmed to carry wsp dna by southern hybridisation, while other individuals in these species were positive for both loci after standard pcr. we considered these four individuals infected with low titre wolbachia because we had previously characterised complete mlst profiles in other individuals of these host species. furthermore, control southern hybridisation to samples with high titre infections demonstrated that the dig-labelled wsp probe did not bind to primer-dimers or spurious wsp products, but hybridised to wsp amplicons that were verified by sequencing. sequence analysis of wsp and 16s rdna thirty-eight flies amplified at the wsp locus and of these also successfully amplified at the 16s rdna locus. sequence analysis identified five different wsp alleles of supergroup a. two wsp alleles were found to co-occur in individuals of six species; beyond this, wsp occurred as a single sequence in one b. frauenfeldi, four b. tryoni and three b. neohumeralis individuals; and wsp in one b. neohumeralis individual. wolbachia mlst characterisation linked wsp and wsp with strains st- and st- respectively, demonstrating a high degree of sharing of wolbachia strains among six tephritid species. pcr-rflp and sequencing confirmed the absence of any other detectable wsp variants in these individuals.table wsp allele numbers and genbank accession numbers of analysed individuals 661 661 a 661 661 661 661 661 661 661 662 amplicons were sequenced and confirmed by pcr-rflp; †wsp allele numbers from wolbachia wsp database, new wsp alleles in bold; a wsp allele number not assigned fig. bayesian inference tree of wsp sequences from wolbachia within tephritid fruit flies and other invertebrates. highlighted sequences are from this study. all other sequences were retrieved from genbank or the wolbachia wsp database. supergroup classification of wsp sequences is marked at the nodes. posterior probabilities > are marked at the nodes; tree was rooted by supergroup d and f strains. scale bar represents the number of nucleotide substitutions per site cloning of wsp from one b. frauenfeldi individual revealed two alleles unlike those found in the other b. frauenfeldi individuals: one sequence identical to wsp of drosophila simulans strain wri; the other with high similarity to two wsp sequences detected in b. dorsalis from china. furthermore, the wsp and mlst alleles were shared with the wolbachia infections found in the parasitoid f. arisanus. all five b. perkinsi individuals produced wsp and 16s rdna amplicons. cloning and sequencing of the wsp fragments revealed a novel allele, wsp as well as multiple clones of a wsp sequence with a single base insertion. this insertion disrupts the open reading frame by incorporating stop codons, and may represent a wolbachia pseudogene transferred to the host genome of this species. amplification of the ftsz locus in all individuals of b. perkinsi have previously failed during the mlst characterisation of wolbachia in this species while all other loci were positive with a group like sequences. we interpreted this as presence of multiple pseudogenes and absence of true wolbachia infections in the analysed specimens of this species. for b. peninsularis, the 16s rdna fragment was consistently amplified from all seven individuals, whereas the wsp locus failed to amplify. direct sequencing of the 16s rdna amplicon showed multiple peaks in the sequence chromatogram indicating two distinct sequences: one full length b group 16s rdna sequence, and one sequence with a 31 bp deletion at positions 182– indicative of a potential pseudogene. the sequence of the b. peninsularis full-length 16s rdna fragment was deposited in genbank, accession number kc failure to detect a wsp sequence in any of these individuals, and prior failure to detect the five mlst markers except fbpa, suggested that this sequence may not represent a genuine wolbachia infection; consequently this species was classified as uninfected. as a control for the 16s rdna pcr assays, one b. neohumeralis sample, id was chosen for 16s rdna sequencing to confirm homology to other wolbachia 16s rdna sequences in genbank. analysis of wolbachia incidence and prevalence geographically as well as irrespective of host species and wolbachia strains, wolbachia was detected in five of samples from equatorial queensland, of samples from tropical queensland and the northern territory, and of samples from subtropical queensland. however, none of the samples from temperate new south wales and none of the samples from central australia had detectable wolbachia. the analysis of wolbachia incidence within species revealed that it was restricted to the northern half of australia; analysis of species for which at least ten individuals per species were available for each half revealed wolbachia incidence in at least one specimen for a number of these species in the northern half of the gradient while none was infected in the southern half. therefore, incidence of wolbachia was higher in the northern half of the gradient, with wolbachia incidence in seven out of nine species for which a minimum of ten individuals were tested; all four species from the southern half for which a minimum of ten individuals per species were available in that region, were negative. a linear model was applied in order to appropriately represent both widespread and more tropically restricted species and this detected a subtle relationship between the prevalence of wolbachia and the midpoint latitude of twelve individual species. for the seven species that had individuals with a full wolbachia mlst profile plus wsp, the prevalence ranged from  % to  %. multivariate analysis of infection prevalence in the only two widespread species that are endemic from equatorial to temperate climes, b. neohumeralis and b. tryoni, detected a significant interaction of wolbachia prevalence by latitude with wolbachia being more common in the northern regions of east coast australia. univariate analysis with adjusted p-values to account for family-wise error across both species did not detect a significant interaction, demonstrating that the detectable effect on a single species is weak, but the pattern is strengthened when analysing both species together.table analyses of wolbachia prevalence at different latitudes of b. neohumeralis and b. tryoni b. neohumeralis b. tryoni table presents analysis of deviance for multivariate and univariate analyses of wolbachia prevalence at different latitudes of two widespread species b. neohumeralis and b. tryoni, as per output of the mvabund package. p-value ) significance codes: ‘***’ ‘**’ ‘*’ ‘.’ ‘’ 1 furthermore, temporal effects on wolbachia prevalence within b. bryoniae, b. frauenfeldi, b. neohumeralis, b. strigifinis and b. tryoni were tested for equatorial and tropical samples by using a fisher’s exact test over the sampling years. no significant temporal change within the tested regions was detected, except for b. tryoni, but this represented neither an overall increase nor decrease in prevalence over time. analysis of mitochondrial dna most tephritid species produced clear coi sequences. bactrocera jarvisi, b. murrayi, b. dorsalis and b. perkinsi individuals produced ambiguous sequences, indicative of potential nuclear mitochondrial dna in these species. the latter sequences were not included in phylogenetic analysis. instead, b. jarvisi and b. dorsalis sequences were retrieved from genbank. bayesian analysis of sequences over 571 bp returned a well-supported consensus tree.fig. bayesian inference tree of coi sequences from tephritid fruit flies. representatives of the five major tephritid fruit fly genera, anastrepha, bactrocera, ceratitis, dacus and rhagoletis as well as dirioxa are shown. dark grey shading highlights the individuals from this study harbouring wolbachia, light grey highlights specimens that were negative for wolbachia. sequences from species without shading were retrieved from genbank. tree was rooted with drosophila melanogaster, and nodes are labelled with posterior probability values. scale bar represents the number of nucleotide substitutions per site mitochondrial haplotypes of morphologically identified individuals of different species harbouring identical wsp alleles were different, except for the two sibling species b. neohumeralis and b. tryoni; both had individuals with identical wolbachia and mitochondrial sequences but without any linkage between wolbachia infection and haplotypes. coi distance measures of b. neohumeralis and b. tryoni showed little difference within infected and uninfected individuals, and between infected and uninfected individuals. discussion we surveyed australian tephritid fruit fly species throughout their natural range, and detected wolbachia sequences in individuals of nine species. in two species, all individuals had pseudogenes and incomplete sets of marker genes – indicative of potential lateral gene transfer to the host genome without the presence of an existing infection. both species are restricted to the tropics. the remaining seven fruit fly species had both 16s rdna and wsp sequences characteristic for four wolbachia strains for which complete mlst profiles had previously been established. the four wolbachia strains occurred as paired or single infections in more than one host species, indicative of horizontal wolbachia transmission events into these species. overall the incidence of wolbachia strains was restricted to tropical regions, and the prevalence of the four wolbachia strains was very low in infected species. in five infected species, wolbachia prevalence appeared to not have changed over time. in b. neohumeralis and b. tryoni, wolbachia was not linked with any particular mitochondrial haplotypes. we thus concluded that the detected strains may not necessarily manipulate reproduction of their host species and maternal inheritance may also be unreliable. therefore, the detection of identical wolbachia strains at low prevalence in several species of this fruit fly community may be due to a spillover from a yet unknown source species, without the propagation of wolbachia in any of the analysed host species, and this spillover appears to be restricted to the tropical regions. wolbachia strains and pseudogenes the previous mlst characterisation in australian tropical tephritids revealed four complete mlst profiles of supergroup a wolbachia in seven species. here, however, 16s rdna and wsp sequences were isolated from nine and eight species, respectively while complete marker gene sets were not found for all nine species, which is likely to be due to lateral gene transfer of wolbachia into the host genome. more specifically, all individuals of b. perkinsi produced 16s rdna and wsp sequences. in the previous study, no individuals of this species amplified ftsz and produced both pseudogenic and novel orf sequences for other mlst loci. another species, b. peninsularis, also had signs of wolbachia lateral gene transfer, as all individuals only amplified 16s rdna but not wsp in this study while previously fbpa was detected as the only positive mlst locus in all individuals of this species. similarly, pseudogenised 16s rdna and fbpa sequences were found in one out of eight tested glossina fly species. we thus conclude that the wolbachia signals in b. peninsularis and b. perkinsi are due to lateral gene transfer from wolbachia into these host genomes, as has been reported for other insects. it may also be noted that the wolbachia pseudogenic sequences were different enough to conclude that they did not derive from the four strains detected in this study. two of the previously characterised wolbachia strains st- and st- were found with their 16s rdna and wsp sequences in individuals of six fruit fly species. single infections by both of these strains were detected in b. neohumeralis individuals and single infections by st- in all four infected b. tryoni individuals and one b. frauenfeldi individual. two additional strains, st- and st- previously characterised by morrow et al. were detected with their 16s rdna and wsp sequence variants in another b. frauenfeldi individual. it was demonstrated that this double infection is identical to wolbachia from the fruit fly parasitoid f. arisanus, and this could be a spillover from this parasitoid species. this conclusion can be drawn based on the result of our earlier study in which all individuals of this parasitoid species were infected, including individuals from the sydney region where fruit fly host populations were uninfected. of additional interest is that both f. arisanus and b. frauenfeldi are not endemic to australia; f. arisanus was introduced to australia in and now occurs from northern queensland to sydney while b. frauenfeldi was first detected on cape york peninsula in from where it has spread throughout northern queensland. this also means that st- and st- were only found in two host species that were originally not native to australia. latitudinal distribution of wolbachia in australian tephritids after excluding pseudogenes in two species as potential wolbachia infections, we detected wolbachia in 29 % of tested australian tephritid species, and this is lower than previous estimates of 40 % incidence of infections in arthropods, mostly based on a european dataset. our results varied between climatic regions. wolbachia incidence was restricted to equatorial and tropical regions north of gladstone, while wolbachia was absent from individuals caught in the more southern subtropical, temperate and arid inland regions. specifically, wolbachia was found in species whose range is limited to tropical and equatorial queensland; in species with a broader climatic range, wolbachia occurred as far south as cairns for b. bryoniae, mackay for b. tryoni and gladstone for b. neohumeralis, but was absent from individuals in more southerly subtropical and temperate regions. within all of these species, wolbachia prevalence was very low, and in the more widely distributed species b. neohumeralis and b. tryoni, wolbachia was also restricted to the northern ranges of each species. wolbachia in australian tephritids: transient infections resulting from a spillover? previously, identical wolbachia strains were found in multiple species of fruit flies that share ecological niches, but are unlikely to hybridise in nature, except possibly for b. tryoni and b. neohumeralis. this study found that only low frequency wolbachia signals were detected in samples from tropical environments, including for those species that are widespread from tropical to temperate climes. furthermore, no wolbachia was detected in any laboratory colonies of these species, including those originating from tropical locations, a result that was supported by pyrosequencing of 16s rdna amplicon libraries from these same individuals. this suggests that there is ample opportunity for the horizontal transfer of wolbachia, but the low prevalence, and lack of increasing prevalence over time, does not unequivocally show that these infections confer any selective advantage to infected females e.g. through reproductive manipulation such as cytoplasmic incompatibility nor that they are being reliably inherited. experimental evidence of a heritable wolbachia infection is necessary to confirm an established infection. in drosophila suzukii, a similarly low frequency wolbachia infection occurs with variable maternal transmission; however no ci or fitness benefits were observed, and maintenance of wolbachia in this species has yet to be explained. here, in the absence of any evidence of vertical inheritance, the finding of two strains in six host species at low prevalence is difficult to explain, and perhaps more likely due to a spillover of wolbachia from an as yet unknown host species without any proliferation in any of the studied species. this interpretation still lacks direct experimental evidence, but transient infections are frequently observed in wolbachia microinjection studies that detect wolbachia in many individuals emerging from microinjected embryos. however, only a few injected embryos result in the establishment of stably inherited matrilines while wolbachia is lost from most other matrilines, or remains undetectable in some matrilines for many generations. in our study, the potential spillover of wolbachia is restricted to the tropical regions, and is perhaps due to the higher biodiversity, resulting in increased ecological interactions in the tropics. the finding of two strains that are shared between fruit fly and parasitoid individuals supports the interpretation of a potential spillover with the parasitoid being the source in this instance; st- and st- were only found in one b. frauenfeldi individual while all f. arisanus individuals were infected with both strains, including from regions where host fruit flies were uninfected. if, however, the infections are heritable in the studied fruit fly species, then the latitudinal restriction of wolbachia incidence and prevalence in australian fruit flies to northern australia may be due to wolbachia loss due to environmental conditions of their hosts in the southern range or an ongoing wolbachia invasion of fruit fly species from north to south, driven by a reproductive phenotype. an increasing number of studies have demonstrated that wolbachia responds to the climatic environment of its host insects. for example, precipitation frequency appeared to correlate with the distribution of single and double infections in leaf beetle individuals in panama, with multiple infections restricted to wetter regions. in another study, multiple infections were found to be more frequent in one tropical habitat when compared with two temperate habitats. temperature is known to affect wolbachia titres in a number of insect hosts. for example, the virulent wmelpop strain appears to over-replicate at elevated temperatures. however, high temperatures may also effectively reduce wolbachia densities in hosts through increased wolbachia-bacteriophage activity, and thus result in diminished penetrance of ci or male-killing phenotypes. it is possible that temperature fluctuations and extremes, commonly found in southern regions of australia, may be selecting more strongly against wolbachia than the more constant temperature conditions in the tropics. such climatic effects could also be contributing to the distribution of wolbachia in previously reported analyses such as d. melanogaster and cat fleas. similarly, climate and latitude were found to determine ranges of other microbial symbioses, for example in marine invertebrates, terrestrial insects and humans. the second scenario supposes that the wolbachia in the australian tephritid species, in particular infections of the widespread b. tryoni and b. neohumeralis, is relatively recently acquired, and a progressive ci driven sweep may therefore be in its infancy. however this hypothesis is less likely for the following reasons; a ci driven invasion should result in an increase of infection prevalence over time that was not detected in these species. furthermore, wolbachia-induced ci is expected to cause a selective sweep of infected mitochondrial haplotypes, yet there does not appear to be support for this in the b. tryoni and b. neohumeralis species complex because both species share mitochondrial haplotypes across infected and uninfected individuals. we were not able to directly assess the ci phenotype of wolbachia in australian tephritids due to the lack of infected laboratory populations. the ci characterisation of the wolbachia strains will thus require future field collection efforts, in particular of females, in order to set up infected laboratory colonies for crossing experiments. testing of field females will also avoid a potential underestimate of wolbachia infection rates due to our male-biased sampling approach that would not detect male-killing wolbachia strains. however, prevalence of male-killing wolbachia may be generally low as found for another fruit fly family, drosophilidae. evidence for numtdna in bactrocera blacket et al. have previously detected numtdna amplicons for b. tryoni and b. neohumeralis by using primers for a different section of coi; our different primer set did not yield numtdna for b. tryoni and b. neohumeralis but revealed potential numtdna sequences in other species such as b. jarvisi, b. dorsalis, b. murrayi and b. perkinsi. interestingly, potential transfer of both mitochondrial and wolbachia genes into the nuclear host genome were detected for b. perkinsi.  in this study of australian tephritids we detected several wolbachia pseudogenes that may be host-genomic footprints of historic infections that have since been lost from host species. within the species of relevance, these pseudogenes were fixed at 100 % while prevalence of wolbachia infections in other species was overly low. under the assumption that the detected wolbachia strains in seven fruit fly species may not be reliably inherited nor cause a reproductive phenotype, these infections could be considered transient and a result of a spillover from a yet unknown source. in our study system, this spillover appears to be restricted to northern australia and has not resulted in further proliferation within any of the tested host lineages, and this warrants future investigation. the fact that we have not found a single case of a high prevalence infected australian tephritid fruit fly species in tested species raises the question whether these taxa exhibit traits that reduce wolbachia infectivity while two species had wolbachia infections in their evolutionary past. other related pest taxa such as olive fly bactrocera oleae appear uninfected in natural populations while laboratory populations were successfully infected; a high incidence of wsp genes was found in thai bactrocera species while prevalence of wsp in b. dorsalis was found to be extremely low, although these strains were not further characterised. perhaps some of these bactrocera species also contain wolbachia signals that are either pseudogenes and/or transient somatic infections, and similar questions may be valid for wolbachia signals found in other host species. our study highlights that tropical insect communities have signs of lateral gene transfer and may be exposed to increased horizontal wolbachia transfer. it also emphasises that global estimates for wolbachia infection frequencies should be interpreted with caution, and account for lateral gene transfer, variable prevalence within species and transient spillover effects. this may be of particular relevance if wolbachia surveys are based on a single or a few marker genes for detection of infections in a host species, and without a broader sampling effort from the host’s distribution. availability of supporting data data sets supporting the results of this article are included within the article and its additional files. dna sequences generated in this study were submitted to genbank: accession numbers kc581371–kc kc kc775785-kc and kc668319-kc kc668323-kc kc and kc775793-kc. new wsp alleles have also been deposited in the wolbachia wsp database. additional files additional file 1: collection data for field-caught fruit flies from this study. information for collection location, climate group and year, specimen specific id number, sex, result for wsp and 16s rdna screening and genbank accession number of coi partial gene are included. additional file 2: laboratory strains screened for wolbachia using wsp and 16s rdna. additional file 3: pcr conditions and primer sequences used in this research. additional file 4: relationship between prevalence of wolbachia infection in twelve tephritid species and the midpoint of their latitudinal distribution, based on hancock et al. [ 32 ] and royer and hancock [ 33 ]. at least ten individuals were tested for uninfected species. additional file 5: comparison of wolbachia infected individuals sampled from equatorial and tropical regions for b. bryoniae, b. frauenfeldi, b. neohumeralis,, b. strigifinis and b. tryoni over five collection years. fisher’s exact test shows overall no significant temporal effect at α = , except for b. tryoni, but this did not represent an overall increase or decrease in prevalence. additional file 6: distance matrix of coi sequences covering 551 bp positions of b. tryoni and b. neohumeralis individuals, grouped by wolbachia -infected and uninfected. below the diagonal are the pairwise number of differences and above the diagonal are p-distance values. p-distances given for within- and between- group means are listed below the main table. competing interests the authors declare that they have no competing interests. authors’ contributions the project was conceived by mr and jlm; the research was designed by jlm and mr with input from mf and dcas; jlm performed the research and analysed the data; mr, mf, jer and dcas supplied materials and specimens; jlm and mr wrote the paper with input by mf, jer and dcas. all authors read and approved the final manuscript. 
for the two remaining species all individuals had incomplete mlst profiles and wolbachia pseudogenes that may be indicative of lateral gene transfer into host genomes. instead, the findings may be indicative of transient infections that result from spillover events from a yet unknown source. it also emphasises that global estimates of wolbachia frequencies may need to consider lateral gene transfer and wolbachia spillover that may be regionally restricted, transient and not inherited. in a previous study of tephritid fruit fly species of tropical australia we detected a high incidence of identical wolbachia strains in several species as well as wolbachia pseudogenes. they can also be transferred horizontally into uninfected host lineages. these spillover events appear to be restricted to northern australia, without proliferation in host lineages further south. here, we have investigated this further by analysing field specimens of species collected along a  km climate gradient of eastern australia. the presence of wolbachia pseudogenes at high prevalence in two species in absence of complete mlst profiles may represent footprints of historic infections that have been lost. our study highlights that tropical fruit fly communities contain wolbachia pseudogenes and may be exposed to frequent horizontal wolbachia transfer. a wolbachia spillover from an infected source population must occur prior to the establishment of heritable infections, but this spillover may be transient. the detection of wolbachia was restricted to northern australia, including in five species that only occur in the tropics. the detection of identical low prevalence strains in a small number of individuals of seven species may question their role as reproductive manipulator and their vertical inheritance. within the more widely distributed bactrocera tryoni and bactrocera neohumeralis, wolbachia also only occurred in the north, and was not linked to any particular mitochondrial haplotypes. wolbachia sequences were detected in individuals of nine of the species. seven species displayed four distinct wolbachia strains based on characterisation of full multi locus sequencing profiles; the strains occurred as single and double infections in a small number of individuals.
 maternally inherited wolbachia bacteria infect many insect species. they can also be transferred horizontally into uninfected host lineages. a wolbachia spillover from an infected source population must occur prior to the establishment of heritable infections, but this spillover may be transient. in a previous study of tephritid fruit fly species of tropical australia we detected a high incidence of identical wolbachia strains in several species as well as wolbachia pseudogenes. here, we have investigated this further by analysing field specimens of species collected along a  km climate gradient of eastern australia.  wolbachia sequences were detected in individuals of nine of the species. seven species displayed four distinct wolbachia strains based on characterisation of full multi locus sequencing profiles; the strains occurred as single and double infections in a small number of individuals. for the two remaining species all individuals had incomplete mlst profiles and wolbachia pseudogenes that may be indicative of lateral gene transfer into host genomes. the detection of wolbachia was restricted to northern australia, including in five species that only occur in the tropics. within the more widely distributed bactrocera tryoni and bactrocera neohumeralis, wolbachia also only occurred in the north, and was not linked to any particular mitochondrial haplotypes.  the presence of wolbachia pseudogenes at high prevalence in two species in absence of complete mlst profiles may represent footprints of historic infections that have been lost. the detection of identical low prevalence strains in a small number of individuals of seven species may question their role as reproductive manipulator and their vertical inheritance. instead, the findings may be indicative of transient infections that result from spillover events from a yet unknown source. these spillover events appear to be restricted to northern australia, without proliferation in host lineages further south. our study highlights that tropical fruit fly communities contain wolbachia pseudogenes and may be exposed to frequent horizontal wolbachia transfer. it also emphasises that global estimates of wolbachia frequencies may need to consider lateral gene transfer and wolbachia spillover that may be regionally restricted, transient and not inherited. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. issue-copyright-statement© the author 2015 