A molecular phylogeny for the Tribe Dacini (Diptera: Tephritidae) : systematic and biogeographic implications

With well over 700 species, the Tribe Dacini is one of the most species-rich clades within the dipteran family Tephritidae, the true fruit flies. Nearly all Dacini belong to one of two very large genera, Dacus Fabricius and Bactrocera Macquart. The distribution of the genera overlap in or around the Indian subcontinent, but the greatest diversity of Dacus is in Africa and the greatest diversity of Bactrocera is in south-east Asia and the Pacific. The monophyly of these two genera has not been rigorously established, with previous phylogenies only including a small number of species and always heavily biased to one genus over the other. Moreover, the subgeneric taxonomy within both genera is complex and the monophyly of many subgenera has not been explicitly tested. Previous hypotheses about the biogeography of the Dacini based on morphological reviews and current distributions of taxa have invoked an out-of-India hypothesis; however this has not been tested in a phylogenetic framework. We attempted to resolve these issues with a dated, molecular phylogeny of 125 Dacini species generated using 16S, COI, COII and white eye genes. The phylogeny shows that Bactrocera is not monophyletic, but rather consists of two major clades: Bactrocera s.s. and the ‘Zeugodacus group of subgenera’ (a recognised, but informal taxonomic grouping of 15 Bactrocera subgenera). This ‘Zeugodacus’ clade is the sister group to Dacus, not Bactrocera and, based on current distributions, split from Dacus before that genus moved into Africa. We recommend that taxonomic consideration be given to raising Zeugodacus to genus level. Supportive of predictions following from the out-of-India hypothesis, the first common ancestor of the Dacini arose in the mid-Cretaceous approximately 80 mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50–35 mya). In contrast, diversification in Dacus, at approximately 65 mya, apparently began much earlier than predicted by the out-of-India hypothesis, suggesting that, if the Dacini arose on the Indian plate, then ancestral Dacus may have left the plate in the mid to late Cretaceous via the well documented India–Madagascar–Africa migration route. We conclude that the phylogeny does not disprove the predictions of an out-of-India hypothesis for the Dacini, although modification of the original hypothesis is required.

Population structure of Bactrocera dorsalis s.s., B. papayae and B. philippinensis (Diptera: Tephritidae) in southeast Asia: evidence for a single species hypothesis using mitochondrial DNA and wingshape data

Background Bactrocera dorsalis s.s. is a pestiferous tephritid fruit fly distributed from Pakistan to the Pacific, with the Thai/Malay peninsula its southern limit. Sister pest taxa, B. papayae and B. philippinensis, occur in the southeast Asian archipelago and the Philippines, respectively. The relationship among these species is unclear due to their high molecular and morphological similarity. This study analysed population structure of these three species within a southeast Asian biogeographical context to assess potential dispersal patterns and the validity of their current taxonomic status. Results Geometric morphometric results generated from 15 landmarks for wings of 169 flies revealed significant differences in wing shape between almost all sites following canonical variate analysis. For the combined data set there was a greater isolation-by-distance (IBD) effect under a ‘non-Euclidean’ scenario which used geographical distances within a biogeographical ‘Sundaland context’ (r2 = 0.772, P < 0.0001) as compared to a ‘Euclidean’ scenario for which direct geographic distances between sample sites was used (r2 = 0.217, P < 0.01). COI sequence data were obtained for 156 individuals and yielded 83 unique haplotypes with no correlation to current taxonomic designations via a minimum spanning network. BEAST analysis provided a root age and location of 540kya in northern Thailand, with migration of B. dorsalis s.l. into Malaysia 470kya and Sumatra 270kya. Two migration events into the Philippines are inferred. Sequence data revealed a weak but significant IBD effect under the ‘non-Euclidean’ scenario (r2 = 0.110, P < 0.05), with no historical migration evident between Taiwan and the Philippines. Results are consistent with those expected at the intra-specific level. Conclusions Bactrocera dorsalis s.s., B. papayae and B. philippinensis likely represent one species structured around the South China Sea, having migrated from northern Thailand into the southeast Asian archipelago and across into the Philippines. No migration is apparent between the Philippines and Taiwan. This information has implications for quarantine, trade and pest management.

Piecing together an integrative taxonomic puzzle: microsatellite, wing shape and aedeagus length analyses of Bactrocera dorsalis s.l. (Diptera: Tephritidae) find no evidence of multiple lineages in a proposed contact zone along the Thai/Malay Peninsula

Bactrocera dorsalis (Hendel) and B. papayae Drew & Hancock represent a closely related sibling species pair for which the biological species limits are unclear; i.e., it is uncertain if they are truely two biological species, or one biological species which has been incorrectly taxonomically split. The geographic ranges of the two taxa are thought to abut or overlap on or around the Isthmus of Kra, a recognised biogeographic barrier located on the narrowest portion of the Thai Peninsula. We collected fresh material of B. dorsalis sensu lato (i.e., B. dorsalis sensu stricto + B. papayae) in a north-south transect down the Thai Peninsula, from areas regarded as being exclusively B. dorsalis s.s., across the Kra Isthmus, and into regions regarded as exclusively B. papayae. We carried out microsatellite analyses and took measurements of male genitalia and wing shape. Both the latter morphological tests have been used previously to separate these two taxa. No significant population structuring was found in the microsatellite analysis and results were consistent with an interpretation of one, predominantly panmictic population. Both morphological datasets showed consistent, clinal variation along the transect, with no evidence for disjunction. No evidence in any tests supported historical vicariance driven by the Isthmus of Kra, and none of the three datasets supported the current taxonomy of two species. Rather, within and across the area of range overlap or abutment between the two species, only continuous morphological and genetic variation was recorded. Recognition that morphological traits previously used to separate these taxa are continuous, and that there is no genetic evidence for population segregation in the region of suspected species overlap, is consistent with a growing body of literature that reports no evidence of biological differentiation between these taxa.

Are insect frugivores always plant pests? The impact of fruit fly (Diptera: Tephritidae) larvae on host plant fitness

Herbivory is generally regarded as negatively impacting on host plant fitness. Frugivorous insects, which feed directly on plant reproductive tissues, are predicted to be particularly damaging to hosts. We tested this prediction with the fruit fly, Bactrocera tryoni, by recording the impact of larval feeding on two direct (seed number and germination) and two indirect (fruit decay rate and attraction/deterrence of vertebrate frugivores) measures of host plant fitness. Experiments were done in the laboratory, glasshouse and tropical rainforest. We found no negative impact of larval feeding on seed number or germination for three test plants: tomato, capsicum and eggplant. Further, larval feeding accelerated the initiation of decay and increased the final level of fruit decay in tomatoes, apples, pawpaw and pear, a result considered to be beneficial to the fruit. In rainforest studies, native rodents preferred infested apple and pears compared to uninfested control fruit; however, there were no differences observed between treatments for tomato and pawpaw. For our study fruits, these results demonstrate that fruit fly larval infestation has neutral or beneficial impacts on the host plant, an outcome which may be largely influenced by the physical properties of the host. These results may contribute to explaining why fruit flies have not evolved the same level of host specialization generally observed for other herbivore groups.

Beyond barcoding : a mitochondrial genomics approach to molecular phylogenetics and diagnostics of blowflies (Diptera: Calliphoridae)

Members of the Calliphoridae (blowflies) are significant for medical and veterinary management, due to the ability of some species to consume living flesh as larvae, and for forensic investigations due to the ability of others to develop in corpses. Due to the difficulty of accurately identifying larval blowflies to species there is a need for DNA-based diagnostics for this family, however the widely used DNA-barcoding marker, cox1, has been shown to fail for several groups within this family. Additionally, many phylogenetic relationships within the Calliphoridae are still unresolved, particularly deeper level relationships. Sequencing whole mt genomes has been demonstrated both as an effective method for identifying the most informative diagnostic markers and for resolving phylogenetic relationships. Twenty-seven complete, or nearly so, mt genomes were sequenced representing 13 species, seven genera and four calliphorid subfamilies and a member of the related family Tachinidae. PCR and sequencing primers developed for sequencing one calliphorid species could be reused to sequence related species within the same superfamily with success rates ranging from 61% to 100%, demonstrating the speed and efficiency with which an mt genome dataset can be assembled. Comparison of molecular divergences for each of the 13 protein-coding genes and 2 ribosomal RNA genes, at a range of taxonomic scales identified novel targets for developing as diagnostic markers which were 117–200% more variable than the markers which have been used previously in calliphorids. Phylogenetic analysis of whole mt genome sequences resulted in much stronger support for family and subfamily-level relationships. The Calliphoridae are polyphyletic, with the Polleninae more closely related to the Tachinidae, and the Sarcophagidae are the sister group of the remaining calliphorids. Within the Calliphoridae, there was strong support for the monophyly of the Chrysomyinae and Luciliinae and for the sister-grouping of Luciliinae with Calliphorinae. Relationships within Chrysomya were not well resolved. Whole mt genome data, supported the previously demonstrated paraphyly of Lucilia cuprina with respect to L. sericata and allowed us to conclude that it is due to hybrid introgression prior to the last common ancestor of modern sericata populations, rather than due to recent hybridisation, nuclear pseudogenes or incomplete lineage sorting.

Comprehensive evaluation of DNA barcoding for the molecular species identification of forensically important Australian Sarcophagidae (Diptera)

Carrion-breeding Sarcophagidae (Diptera) can be used to estimate the post-mortem interval (PMI) in forensic cases. Difficulties with accurate morphological identifications at any life stage and a lack of documented thermobiological profiles have limited their current usefulness of these flies. The molecular-based approach of DNA barcoding, which utilises a 648-bp fragment of the mitochondrial cytochrome oxidase subunit I gene, was previously evaluated in a pilot study for the discrimination between 16 Australian sarcophagids. The current study comprehensively evaluated DNA barcoding on a larger taxon set of 588 adult Australian sarcophagids. A total of 39 of the 84 known Australian species were represented by 580 specimens, which includes 92% of potentially forensically important species. A further eight specimens could not be reliably identified, but included as six unidentifable taxa. A neighbour-joining phylogenetic tree was generated and nucleotide sequence divergences were calculated using the Kimura-two-parameter distance model. All species except Sarcophaga (Fergusonimyia) bancroftorum, known for high morphological variability, were resolved as reciprocally monophyletic (99.2% of cases), with most having bootstrap support of 100. Excluding S. bancroftorum, the mean intraspecific and interspecific variation ranged from 0.00-1.12% and 2.81-11.23%, respectively, allowing for species discrimination. DNA barcoding was therefore validated as a suitable method for the molecular identification of the Australian Sarcophagidae, which will aid in the implementation of this fauna in forensic entomology.

The complete mitochondrial genome of the flesh fly, Sarcophaga impatiens Walker (Diptera: Sarcophagidae)

Approximately 2500 fly species comprise the Sarcophagidae family worldwide. The complete mitochondrial genome of the carrion-breeding, forensically important Sarcophaga impatiens Walker (Diptera: Sarcophagidae) from Australia was sequenced. The 15,169 bp circular genome contains the 37 genes found in a typical Metazoan genome: 13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes. It also contains one non-coding A+T-rich region. The arrangement of the genes was the same as that found in the ancestral insect. All the protein initiation codons are ATN, except for cox1 that begins with TCG (encoding S). The 22 tRNA anticodons of S. impatiens are consistent with those observed in Drosophila yakuba, and all form the typical cloverleaf structure, except for tRNA-Ser(AGN) that lacks the DHU arm. The mitochondrial genome of Sarcophaga presented will be valuable for resolving phylogenetic relationships within the family Sarcophagidae and the order Diptera, and could be used to identify favourable genetic markers for species identifications for forensic purposes.

Spatial and temporal foraging patterns of Queensland fruit fly, Bactrocera tryoni(Froggatt) (Diptera: Tephritidae), for protein and implications for management

Fruit flies require protein for reproductive development and actively feed upon protein sources in the field. Liquid protein baits mixed with insecticide are used routinely to manage pest fruit flies, such as Bactrocera tryoni (Froggatt). However, there are still some gaps in the underpinning science required to improve the efficacy of bait spray technology. The spatial and temporal foraging behaviour of B. tryoni in response to protein was investigated in the field. A series of linked trials using either wild flies in the open field or laboratory-reared flies in field cages and a netted orchard were undertaken using nectarines and guavas. Key questions investigated were the fly's response to protein relative to: height of protein within the canopy, fruiting status of the tree, time of day, season and size of the experimental arena. Canopy height had a significant response on B. tryoni foraging, with more flies foraging on protein in the mid to upper canopy. Fruiting status also had a significant effect on foraging, with most flies responding to protein when applied to fruiting hosts. B. tryoni demonstrated a repeatable diurnal response pattern to protein, with the peak response being between 12:00–16:00 h. Season showed significant but unpredictable effects on fruit fly response to protein in the subtropical environment where the work was undertaken. Relative humidity, but not temperature or rainfall, was positively correlated with protein response. The number of B. tryoni responding to protein decreased dramatically as the spatial scale increased from field cage through to the open field. Based on these results, it is recommend that, to be most effective, protein bait sprays should be applied to the mid to upper canopies of fruiting hosts. Overall, the results show that the protein used, an industry standard, has very low attractancy to B. tryoni and that further work is urgently needed to develop more volatile protein baits.

Five commercial citrus rate poorly as hosts of the polyphagous fruit flyBactrocera tryoni(Froggatt) (Diptera: Tephritidae) in laboratory studies

Queensland fruit fly, Bactrocera tryoni (Froggatt), is a polyphagous pest, and many citrus types are included among its hosts. While quantification of citrus host use by B. tryoni is lacking, citrus is generally considered a ‘low pressure’ crop. This paper investigates B. tryoni female oviposition preference and offspring performance in five citrus types; Murcott mandarin (Citrus reticulata), Navel orange and Valencia orange (Citrus sinensis), Eureka lemon (Citrus limon) and yellow grapefruit (Citrus paradisi). Oviposition preference was investigated in laboratory-based choice and no-choice experiments, while immature survival and offspring performance were investigated by infesting fruits in the laboratory and evaluating pupal recovery, pupal emergence and F1 fecundity. Fruit size, Brix level and peel toughness were also measured for correlation with host use. Bactrocera tryoni demonstrated an oviposition preference hierarchy among the citrus fruits tested; Murcott and grapefruit were most preferred for oviposition and lemon the least, while preference for Navel and Valencia was intermediate. Peel toughness was negatively correlated with B. tryoni oviposition preference, while no significant correlations were detected between oviposition and Brix level or fruit size. Immature survival in the tested fruit was very low. Murcott was the best host (21% pupal recovery), while all other citrus types that showed pupal recovery of 6% or lower and no pupae were recovered from Valencia orange. In pupae recovered from Navel orange and lemon, adult eclosion was greatly reduced, while in grapefruit and lemon, no eggs were recovered from F1 adults. Based on these laboratory results, many commercial citrus varieties appear to be poor hosts for B. tryoni and may pose a low post-harvest and quarantine risk. These findings need to be confirmed in the field, as they impact on both pre-harvest and post-harvest countermeasures.

Utility of COI, CAD and morphological data for resolving relationships within the genus Sarcophaga (sensu lato) (Diptera: Sarcophagidae) : a preliminary study

Currently there are ~3000 known species of Sarcophagidae (Diptera), which are classified into 173 genera in three subfamilies. Almost 25% of sarcophagids belong to the genus Sarcophaga (sensu lato) however little is known about the validity of, and relationships between the ~150 (or more) subgenera of Sarcophaga s.l. In this preliminary study, we evaluated the usefulness of three sources of data for resolving relationships between 35 species from 14 Sarcophaga s.l. subgenera: the mitochondrial COI barcode region, ~800. bp of the nuclear gene CAD, and 110 morphological characters. Bayesian, maximum likelihood (ML) and maximum parsimony (MP) analyses were performed on the combined dataset. Much of the tree was only supported by the Bayesian and ML analyses, with the MP tree poorly resolved. The genus Sarcophaga s.l. was resolved as monophyletic in both the Bayesian and ML analyses and strong support was obtained at the species-level. Notably, the only subgenus consistently resolved as monophyletic was Liopygia. The monophyly of and relationships between the remaining Sarcophaga s.l. subgenera sampled remain questionable. We suggest that future phylogenetic studies on the genus Sarcophaga s.l. use combined datasets for analyses. We also advocate the use of additional data and a range of inference strategies to assist with resolving relationships within Sarcophaga s.l.

The evolution and biogeography of the austral horse fly tribe Scionini (Diptera: Tabanidae: Pangoniinae) inferred from multiple mitochondrial and nuclear genes

Phylogenetic relationships within the Tabanidae are largely unknown, despite their considerable medical and ecological importance. The first robust phylogenetic hypothesis for the horse fly tribe Scionini is provided, completing the systematic placement of all tribes in the subfamily Pangoniinae. The Scionini consists of seven mostly southern hemisphere genera distributed in Australia, New Guinea, New Zealand and South America. A 5757. bp alignment of 6 genes, including mitochondrial (COI and COII), ribosomal (28S) and nuclear (AATS and CAD regions 1, 3 and 4) genes, was analysed for 176 taxa using both Bayesian and maximum likelihood approaches. Results indicate the Scionini are strongly monophyletic, with the exclusion of the only northern hemisphere genus Goniops. The South American genera Fidena, Pityocera and Scione were strongly monophyletic, corresponding to current morphology-based classification schemes. The most widespread genus Scaptia was paraphyletic and formed nine strongly supported monophyletic clades, each corresponding to either the current subgenera or several previously synonymised genera that should be formally resurrected. Molecular results also reveal a newly recognised genus endemic to New Zealand, formerly placed within Scaptia. Divergence time estimation was employed to assess the global biogeographical patterns in the Pangoniinae. These analyses demonstrated that the Scionini are a typical Gondwanan group whose diversification was influenced by the fragmentation of that ancient land mass. Furthermore, results indicate that the Scionini most likely originated in Australia and subsequently radiated to New Zealand and South American by both long distance dispersal and vicariance. The phylogenetic framework of the Scionini provided herein will be valuable for taxonomic revisions of the Tabanidae.

Mating compatibility among four pest members of the Bactrocera dorsalis fruit fly species complex (Diptera : Tephritidae)

Bactrocera dorsalis (Hendel), Bactrocera papayae Drew & Hancock, Bactrocera philippinensis Drew & Hancock, and Bactrocera carambolae Drew & Hancock are pest members within the B. dorsalis species complex of tropical fruit flies. The species status of these taxa is unclear and this confounds quarantine, pest management, and general research. Mating studies carried out under uniform experimental conditions are required as part of resolving their species limits. These four taxa were collected from the wild and established as laboratory cultures for which we subsequently determined levels of prezygotic compatibility, assessed by field cage mating trials for all pair-wise combinations. We demonstrate random mating among all pair-wise combinations involving B. dorsalis, B. papayae, and B. philippinensis. B. carambolae was relatively incompatible with each of these species as evidenced by nonrandom mating for all crosses. Reasons for incompatibility involving B. carambolae remain unclear; however, we observed differences in the location of couples in the field cage for some comparisons. Alongside other factors such as pheromone composition or other courtship signals, this may lead to reduced interspecific mating compatibility with B. carambolae. These data add to evidence that B. dorsalis, B. papayae, and B. philippinensis represent the same biological species, while B. carambolae remains sufficiently different to maintain its current taxonomic identity. This poses significant implications for this group's systematics, impacting on pest management, and international trade.

Multi-gene phylogenetic analysis of south-east Asian pest members of the Bactrocera dorsalis species complex (Diptera: Tephritidae) does not support current taxonomy

Bactrocera dorsalis sensu stricto, B. papayae, B. philippinensis and B. carambolae are serious pest fruit fly species of the B. dorsalis complex that predominantly occur in south-east Asia and the Pacific. Identifying molecular diagnostics has proven problematic for these four taxa, a situation that cofounds biosecurity and quarantine efforts and which may be the result of at least some of these taxa representing the same biological species. We therefore conducted a phylogenetic study of these four species (and closely related outgroup taxa) based on the individuals collected from a wide geographic range; sequencing six loci (cox1, nad4-3′, CAD, period, ITS1, ITS2) for approximately 20 individuals from each of 16 sample sites. Data were analysed within maximum likelihood and Bayesian phylogenetic frameworks for individual loci and concatenated data sets for which we applied multiple monophyly and species delimitation tests. Species monophyly was measured by clade support, posterior probability or bootstrap resampling for Bayesian and likelihood analyses respectively, Rosenberg's reciprocal monophyly measure, P(AB), Rodrigo's (P(RD)) and the genealogical sorting index, gsi. We specifically tested whether there was phylogenetic support for the four 'ingroup' pest species using a data set of multiple individuals sampled from a number of populations. Based on our combined data set, Bactrocera carambolae emerges as a distinct monophyletic clade, whereas B. dorsalis s.s., B. papayae and B. philippinensis are unresolved. These data add to the growing body of evidence that B. dorsalis s.s., B. papayae and B. philippinensis are the same biological species, which poses consequences for quarantine, trade and pest management.