Tephritid Integrative Taxonomy: Where We Are Now, with a Focus on the Resolution of Three Tropical Fruit Fly Species Complexes.

Accurate species delimitation underpins good taxonomy. Formalization of integrative taxonomy in the past decade has provided a framework for using multidisciplinary data to make species delimitation hypotheses more rigorous. We address the current state of integrative taxonomy by using as a case study an international project targeted at resolving three important tephritid species complexes: Bactrocera dorsalis complex, Anastrepha fraterculus complex, and Ceratitis FAR (C. fasciventris, C. anonae, C. rosa) complex. The integrative taxonomic approach has helped deliver significant advances in resolving these complexes: It has been used to identify some taxa as belonging to the same biological species as well as to confirm hidden cryptic diversity under a single taxonomic name. Nevertheless, the general application of integrative taxonomy has not been without issue, revealing challenges that must be considered when undertaking an integrative taxonomy project. Scrutiny of this international case study provides a unique opportunity to document lessons learned for the benefit of not only tephritid taxonomists, but also the wider taxonomic community.

Sexual selection in true fruit flies (Diptera: Tephritidae): transcriptome and experimental evidences for phytochemicals increasing male competitive ability.

In male tephritid fruit flies of the genus Bactrocera, feeding on secondary plant compounds (sensu lato male lures = methyl eugenol, raspberry ketone and zingerone) increases male mating success. Ingested male lures alter the male pheromonal blend, normally making it more attractive to females and this is considered the primary mechanism for the enhanced mating success. However, the male lures raspberry ketone and zingerone are known, across a diverse range of other organisms, to be involved in increasing energy metabolism. If this also occurs in Bactrocera, then this may represent an additional benefit to males as courtship is metabolically expensive and lure feeding may increase a fly's short-term energy. We tested this hypothesis by performing comparative RNA-seq analysis between zingerone-fed and unfed males of Bactrocera tryoni. We also carried out behavioural assays with zingerone- and cuelure-fed males to test whether they became more active. RNA-seq analysis revealed, in zingerone-fed flies, up-regulation of 3183 genes with homologues transcripts to those known to regulate intermale aggression, pheromone synthesis, mating and accessory gland proteins, along with significant enrichment of several energy metabolic pathways and gene ontology terms. Behavioural assays show significant increases in locomotor activity, weight reduction and successful mating after mounting; all direct/indirect measures of increased activity. These results suggest that feeding on lures leads to complex physiological changes, which result in more competitive males. These results do not negate the pheromone effect, but do strongly suggest that the phytochemical-induced sexual selection is governed by both female preference and male competitive mechanisms.

Evaluation of identification methods for cryptic Bactrocera dorsalis (Diptera: Tephritidae) specimens: combining morphological and molecular techniques.

The potential for population genomics to elucidate invasion pathways of a species is limited by taxonomic identification issues. The Oriental fruit fly pest, Bactrocera dorsalis (Hendel) belongs to a complex in which several sympatric species are attracted to the same lure used in trapping and are morphologically cryptic and/or reported to hybridize. In this study, we evaluated the taxonomic ambiguity between B. dorsalis and 2 major cryptic species, based on morphological expertise and 289 target specimens sampled across the whole distribution range. Specimens were then subjected to DNA sequence analyses of the COI mitochondrial barcode and the EIF3L nuclear marker to evaluate the potential for molecular identification, in particular for specimens for which morphological identification was inconclusive. To this aim, we produced reference datasets with DNA sequences from target specimens whose morphological identification was unambiguous, which we complemented with 56 new DNA sequences from closest relatives and 76 published and curated DNA sequences of different species in the complex. After the necessary morphological observation, about 3.5% of the target dataset and 47.6% of the specimens from Southeast Asian islands displayed ambiguous character states shared with B. carambolae and/or B. occipitalis. Critical interpretation of DNA sequence data solved morphological ambiguities only when combining both mitochondrial and nuclear markers. COI discriminated B. dorsalis from 5 species; EIF3L and ITS from another species. We recommend this procedure to ensure correct identification of B. dorsalis specimens in population genetics studies and surveillance programs.

Development of a cost-effective, morphology-preserving method for DNA isolation from bulk invertebrate trap catches: Tephritid fruit flies as an exemplar.

Insect identification and preservation of voucher specimens is integral to pest diagnostic and surveillance activities; yet bulk-trapped insects are a diagnostic challenge due to high catch numbers and the susceptibility of samples to environmental damage. Many insect trap catches rely on examination of morphological characters for species identifications, which is a time consuming and highly skilled task, hence there is a need for more efficient molecular approaches. Many bulk DNA extraction methods require destructive sampling of specimens, resulting in damaged, or fully destroyed, voucher specimens. We developed an inexpensive, rapid, bulk DNA isolation method that preserves specimens as pinned vouchers to a standard that allows for post-extraction morphological examination and inclusion in insect reference collections. Our protocol was validated using a group of insects that are time-consuming to identify when trapped in large numbers-the dacine fruit flies (Diptera: Tephritidae: Dacinae). In developing our method, we evaluated existing protocols against the following criteria: effect on morphology; suitability for large trap catches; cost; ease of handling; and application to downstream molecular diagnostic analyses such as real-time PCR and metabarcoding. We found that the optimum method for rapid isolation of DNA extraction was immersing flies in a NaOH:TE buffer at 75°C for 10 minutes, without the need for proteinase K or detergents. This HotSOAK method produced sufficient high-quality DNA whilst preserving morphological characters suitable for species-level identification with up to 20,000 flies in a sample. The lysates performed well in down-stream analyses such as loop-mediated isothermal amplification (LAMP) and real-time PCR applications, while for metabarcoding PCR the lysate required an additional column purification step. Development of this method is a key step required for upscaling our capacity to accurately detect insects captured in bulk traps, whether for biodiversity, biosecurity, or pest management objectives.

A novel diagnostic gene region for distinguishing between two pest fruit flies: Bactrocera tryoni (Froggatt) and Bactrocera neohumeralis (Hardy) (Diptera: Tephritidae).

Bactrocera tryoni and Bactrocera neohumeralis are morphologically similar sibling pest fruit fly species that possess different biological attributes, geographic distributions, and host ranges. The need to differentiate between the two species is critical for accurate pest status assessment, management, biosecurity, and maintenance of reference colonies. While morphologically similar, adults may be separated based on subtle characters; however, some characters exhibit intraspecific variability, creating overlap between the two species. Additionally, there is currently no single molecular marker or rapid diagnostic assay that can reliably distinguish between B. neohumeralis and B. tryoni; therefore, ambiguous samples remain undiagnosed. Here we report the first molecular marker that can consistently distinguish between B. tryoni and B. neohumeralis. Our diagnostic region consists of two adjacent single nucleotide polymorphisms (SNPs) within the pangolin (pan) gene region. We confirmed the genotypes of each species are consistent across their distributional range, then developed a tetra-primer amplification refractory mutation system (ARMS) PCR assay for rapid diagnosis of the species. The assay utilizes four primers in multiplex, with two outer universal primers, and two internal primers: one designed to target two adjacent SNPs (AA) present in B. tryoni and the other targeting adjacent SNPs present in B. neohumeralis (GG). The assay accurately discriminates between the two species, but their SNP genotypes are shared with other nontarget tephritid fruit fly species. Therefore, this assay is most suited to adult diagnostics where species confirmation is necessary in determining ambiguous surveillance trap catches; maintaining pure colony lines; and in Sterile Insect Technique management responses.

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 wing-shape 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' (r(2) = 0.772, P < 0.0001) as compared to a 'Euclidean' scenario for which direct geographic distances between sample sites was used (r(2) = 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 (r(2) = 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.

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 80mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50-35mya). In contrast, diversification in Dacus, at approximately 65mya, 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.

A diagnostic LAMP assay for rapid identification of an invasive plant pest, fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae).

Fall armyworm (FAW), Spodoptera frugiperda (Lepidoptera: Noctuidae), is a highly polyphagous invasive plant pest that has expanded its global geographic distribution, including recently into much of Australia. Rapid diagnostic tests are required for identification of FAW to assist subsequent management and control. We developed a new loop-mediated isothermal amplification (LAMP) assay based on the mitochondrial cytochrome c oxidase subunit I (COI) gene for accurate and timely diagnosis of FAW in the field. The specificity of the new assay was tested against a broad panel of twenty non-target noctuids, including eight other Spodoptera species. Only S. frugiperda samples produced amplification within 20 min, with an anneal derivative temperature of 78.3 ± 0.3 °C. A gBlock dsDNA fragment was developed and trialled as a synthetic positive control, with a different anneal derivative of 81 °C. The new FAW LAMP assay was able to detect FAW DNA down to 2.4 pg, similar to an existing laboratory-based real-time PCR assay. We also trialled the new FAW assay with a colorimetric master mix and found it could successfully amplify positive FAW samples in half the time compared to an existing FAW colorimetric LAMP assay. Given the high sensitivity and rapid amplification time, we recommend the use of this newly developed FAW LAMP assay in a portable real-time fluorometer for in-field diagnosis of FAW.

Loop-mediated isothermal amplification (LAMP) assays for detection of the New Guinea fruit fly Bactrocera trivialis (Drew) (Diptera: Tephritidae).

The cue-lure-responding New Guinea fruit fly, Bactrocera trivialis, poses a biosecurity risk to neighbouring countries, e.g., Australia. In trapping programs, lure caught flies are usually morphologically discriminated from non-target species; however, DNA barcoding can be used to confirm similar species where morphology is inconclusive, e.g., Bactrocera breviaculeus and B. rufofuscula. This can take days-and a laboratory-to resolve. A quicker, simpler, molecular diagnostic assay would facilitate a more rapid detection and potential incursion response. We developed LAMP assays targeting cytochrome c oxidase subunit I (COI) and Eukaryotic Translation Initiation Factor 3 Subunit L (EIF3L); both assays detected B. trivialis within 25 min. The BtrivCOI and BtrivEIF3L assay anneal derivatives were 82.7 ± 0.8 °C and 83.3 ± 1.3 °C, respectively, detecting down to 1 × 101 copies/µL and 1 × 103 copies/µL, respectively. Each assay amplified some non-targets from our test panel; however notably, BtrivCOI eliminated all morphologically similar non-targets, and combined, the assays eliminated all non-targets. Double-stranded DNA gBlocks were developed as positive controls; anneal derivatives for the COI and EIF3L gBlocks were 84.1 ± 0.7 °C and 85.8 ± 0.2 °C, respectively. We recommend the BtrivCOI assay for confirmation of suspect cue-lure-trapped B. trivialis, with BtrivEIF3L used for secondary confirmation when required.

Development of internal COI primers to improve and extend barcoding of fruit flies (Diptera: Tephritidae: Dacini).

Accurate species-level identifications underpin many aspects of basic and applied biology; however, identifications can be hampered by a lack of discriminating morphological characters, taxonomic expertise or time. Molecular approaches, such as DNA "barcoding" of the cytochrome c oxidase (COI) gene, are argued to overcome these issues. However, nuclear encoding of mitochondrial genes (numts) and poor amplification success of suboptimally preserved specimens can lead to erroneous identifications. One insect group for which these molecular and morphological problems are significant are the dacine fruit flies (Diptera: Tephritidae: Dacini). We addressed these issues associated with COI barcoding in the dacines by first assessing several "universal" COI primers against public mitochondrial genome and numt sequences for dacine taxa. We then modified a set of four primers that more closely matched true dacine COI sequence and amplified two overlapping portions of the COI barcode region. Our new primers were tested alongside universal primers on a selection of dacine species, including both fresh preserved and decades-old dry specimens. Additionally, Bactrocera tryoni mitochondrial and nuclear genomes were compared to identify putative numts. Four numt clades were identified, three of which were amplified using existing universal primers. In contrast, our new primers preferentially amplified the "true" mitochondrial COI barcode in all dacine species tested. The new primers also successfully amplified partial barcodes from dry specimens for which full length barcodes were unobtainable. Thus we recommend these new primers be incorporated into the suites of primers used by diagnosticians and quarantine labs for the accurate identification of dacine species.

Population structure of a global agricultural invasive pest, Bactrocera dorsalis (Diptera: Tephritidae).

Bactrocera dorsalis, the Oriental fruit fly, is one of the world's most destructive agricultural insect pests and a major impediment to international fresh commodity trade. The genetic structuring of the species across its entire geographic range has never been undertaken, because under a former taxonomy B. dorsalis was divided into four distinct taxonomic entities, each with their own, largely non-overlapping, distributions. Based on the extensive sampling of six a priori groups from 63 locations, genetic and geometric morphometric datasets were generated to detect macrogeographic population structure, and to determine prior and current invasion pathways of this species. Weak population structure and high genetic diversity were detected among Asian populations. Invasive populations in Africa and Hawaii are inferred to be the result of separate, single invasions from South Asia, while South Asia is also the likely source of other Asian populations. The current northward invasion of B. dorsalis into Central China is the result of multiple, repeated dispersal events, most likely related to fruit trade. Results are discussed in the context of global quarantine, trade, and management of this pest. The recent expansion of the fly into temperate China, with very few associated genetic changes, clearly demonstrates the threat posed by this pest to ecologically similar areas in Europe and North America.

Effect of Body Size, Age, and Premating Experience on Male Mating Success in Bactrocera tryoni (Diptera: Tephritidae).

Variation in male body size, age, and prior sexual experience may all influence male mating success in tephritid fruit flies. Bactrocera tryoni (Froggatt) is an Australian pest tephritid for which the sterile insect technique (SIT) is being actively pursued, and for which information on what makes males more competitive is urgently needed. Pair-wise competitive mating trials were run using laboratory-reared flies in walk-in field cages, evaluating young, large, and virgin B. tryoni males against old, small, and nonvirgin males, respectively. Analysis of male sexual competitiveness indices revealed that young and large males obtained significantly more copulations compared to old and small males; there was no significant difference between virgin and nonvirgin males in obtaining mates. While SIT programs will always release young males, the results do show that rearing programs which focus on producing larger males, rather than greater numbers of smaller males, will produce more sexually competitive males. After release, virgin SIT males will not be at a competitive disadvantage with sexually experienced males based on prior mating experience.

The mating system of the true fruit fly Bactrocera tryoni and its sister species, Bactrocera neohumeralis.

The frugivorous "true" fruit fly, Bactrocera tryoni (Queensland fruit fly), is presumed to have a nonresourced-based lek mating system. This is largely untested, and contrary data exists to suggest Bactrocera tryoni may have a resource-based mating system focused on fruiting host plants. We tested the mating system of Bactrocera tryoni, and its close sibling Bactrocera neohumeralis, in large field cages using laboratory reared flies. We used observational experiments that allowed us to determine if: (i) mating pairs were aggregated or nonaggregated; (ii) mating system was resource or nonresource based; (iii) flies utilized possible landmarks (tall trees over short) as mate-rendezvous sites; and (iv) males called females from male-dominated leks. We recorded nearly 250 Bactrocera tryoni mating pairs across all experiments, revealing that: (i) mating pairs were aggregated; (ii) mating nearly always occurred in tall trees over short; (iii) mating was nonresource based; and (iv) that males and females arrived at the mate-rendezvous site together with no evidence that males preceded females. Bactrocera neohumeralis copulations were much more infrequent (only 30 mating pairs in total), but for those pairs there was a similar preference for tall trees and no evidence of a resource-based mating system. Some aspects of Bactrocera tryoni mating behavior align with theoretical expectations of a lekking system, but others do not. Until evidence for unequivocal female choice can be provided (as predicted under a true lek), the mating system of Bactrocera tryoni is best described as a nonresource based, aggregation system for which we also have evidence that land-marking may be involved.

Effects of laboratory colonization on Bactrocera dorsalis (Diptera, Tephritidae) mating behaviour: 'what a difference a year makes'.

Laboratory-reared insects are widely known to have significantly reduced genetic diversity in comparison to wild populations; however, subtle behavioural changes between laboratory-adapted and wild or 'wildish' (i.e., within one or very few generations of field collected material) populations are less well understood. Quantifying alterations in behaviour, particularly sexual, in laboratory-adapted insects is important for mass-reared insects for use in pest management strategies, especially those that have a sterile insect technique component. We report subtle changes in sexual behaviour between 'wildish' Bactrocera dorsalis flies (F1 and F2) from central and southern Thailand and the same colonies 12 months later when at six generations from wild. Mating compatibility tests were undertaken under standardised semi-natural conditions, with number of homo/heterotypic couples and mating location in field cages analysed via compatibility indices. Central and southern populations of Bactrocera dorsalis displayed positive assortative mating in the 2010 trials but mated randomly in the 2011 trials. 'Wildish' southern Thailand males mated significantly earlier than central Thailand males in 2010; this difference was considerably reduced in 2011, yet homotypic couples from southern Thailand still formed significantly earlier than all other couple combinations. There was no significant difference in couple location in 2010; however, couple location significantly differed among pair types in 2011 with those involving southern Thailand females occurring significantly more often on the tree relative to those with central Thailand females. Relative participation also changed with time, with more southern Thailand females forming couples relative to central Thailand females in 2010; this difference was considerably decreased by 2011. These results reveal how subtle changes in sexual behaviour, as driven by laboratory rearing conditions, may significantly influence mating behaviour between laboratory-adapted and recently colonised tephritid fruit flies over a relatively short period of time.