Phylogenetic relationships among fruit flies, Bactrocera (Diptera, Tephritidae), based on the mitochondrial rDNA sequences

A transposon-based genetic marker for conspecific identity within the Bactrocera dorsalis species complex

Here we describe a molecular approach to assess conspecific identity that relies on the comparison of an evolved mutated transposable element sequence and its genomic insertion site in individuals from closely related species. This was explored with the IFP2 piggyBac transposon, originally discovered in Trichoplusia ni as a 2472 bp functional element, that was subsequently found as mutated elements in seven species within the Bactrocera dorsalis species complex. In a B. dorsalis [Hendel] strain collected in Kahuku, Hawaii, a degenerate 2420 bp piggyBac sequence (pBacBd-Kah) having ~ 94.5% sequence identity to IFP2 was isolated, and it was reasoned that common species, or strains within species, should share the same evolved element and its precise genomic insertion site. To test this assumption, PCR using primers to pBacBd-Kah and adjacent genomic sequences was used to isolate and compare homologous sequences in strains of four sibling species within the complex. Three of these taxa, B. papayae, B. philippinensis, and B. invadens, were previously synonymized with B. dorsalis, and found to share nearly identical pBacBd-Kah homologous elements (> 99% nucleotide identity) within the identical insertion site consistent with conspecific species. The fourth species tested, B. carambolae, considered to be a closely related yet independent species sympatric with B. dorsalis, also shared the pBacBd-Kah sequence and insertion site in one strain from Suriname, while another divergent pBacBd-Kah derivative, closer in identity to IFP2, was found in individuals from French Guiana, Bangladesh and Malaysia. This data, along with the absence of pBacBd-Kah in distantly related Bactrocera, indicates that mutated descendants of piggyBac, as well as other invasive mobile elements, could be reliable genomic markers for common species identity.

The Impact of Fast Radiation on the Phylogeny of Bactrocera Fruit Flies as Revealed by Multiple Evolutionary Models and Mutation Rate-Calibrated Clock

Several true fruit flies (Tephritidae) cause major damage to agriculture worldwide. Among them, species of the genus Bactrocera are extensively studied to understand the traits associated with their invasiveness and ecology. Comparative approaches based on a reliable phylogenetic framework are particularly effective, but several nodes of the Bactrocera phylogeny are still controversial, especially concerning the reciprocal affinities of the two major pests B. dorsalis and B. tryoni. Here, we analyzed a newly assembled genomic-scaled dataset using different models of evolution to infer a phylogenomic backbone of ten representative Bactrocera species and two outgroups. We further provide the first genome-scaled inference of their divergence by calibrating the clock using fossil records and the spontaneous mutation rate. The results reveal a closer relationship of B. dorsalis with B. latifrons than to B. tryoni, contrary to what was previously supported by mitochondrial-based phylogenies. By employing coalescent-aware and heterogeneous evolutionary models, we show that this incongruence likely derives from a hitherto undetected systematic error, exacerbated by incomplete lineage sorting and possibly hybridization. This agrees with our clock analysis, which supports a rapid and recent radiation of the clade to which B. dorsalis, B. latifrons and B. tryoni belong. These results provide a new picture of Bactrocera phylogeny that can serve as the basis for future comparative analyses.

A comprehensive phylogeny helps clarify the evolutionary history of host breadth and lure response in the Australian Dacini fruit flies (Diptera: Tephritidae)

The tribe Dacini (Diptera: Tephritidae) contains over 930 recognised species and has been widely studied due to the economic importance of some taxa, such as the Oriental fruit fly Bactrocera dorsalis. Despite the attention this group has received, very few phylogenetic reconstructions have comprehensively sampled taxa from a single biogeographic region, thereby limiting our capacity to address more targeted evolutionary questions. To study the evolution of diet breadth and male lure response, two key traits fundamental to understanding dacine diversity and the biology of pest taxa, we analysed 273 individuals representing 144 described species from Australia (80% continental coverage), the Pacific, and select close relatives from South-east Asia to estimate a dated molecular phylogenetic reconstruction of the Dacini. We utilised seven loci with a combined total of 4,332 nucleotides, to estimate both Bayesian and Maximum Likelihood phylogenies of the tribe. Consistent with other molecular phylogenies of the tribe, there was a high level of disagreement between the placement of species in the phylogeny and their current subgeneric and species-complex level taxonomies. The Australian fauna exhibit high levels of endemism, with radiations of both exclusively Australian clades, and clades that originate elsewhere (e.g. the Bactrocera dorsalis species group). Bidirectional movement of species has occurred between Papua New Guinea and Australia, with evidence for multiple incursions over evolutionary time. The Bactrocera aglaiae species group emerged sister to all other Bactrocera species examined. Divergence time estimates were ∼ 30 my younger than previously reported for this group, with the tribe diverging from its most recent common ancestor ∼ 43 mya. Ancestral trait reconstruction and tests for trait phylogenetic signal revealed a strong signal for the evolution of male lure response across the tree, with cue-lure/raspberry ketone lure response the ancestral trait. Methyl eugenol response has arisen on multiple, independent occasions. The evolution of host breadth exhibited a weaker signal; yet, basal groups were more likely to be host specialists. Both the evolution of lure response and host fruit use provide predictive information for the outbreak management of understudied pest fruit flies for which direct inference of these features may be lacking. Our results, which parallel those of earlier research into the closely-related African Dacus spp., demonstrate how geographically focussed taxon coverage allows Dacini phylogenetics to more explicitly test evolutionary hypotheses, thereby progressing our understanding of the evolution of this highly diverse and recently-radiated group of flies.

Host shifts in economically significant fruit flies (Diptera: Tephritidae) with high degree of polyphagy

Insects tend to feed on related hosts. Coevolution tends to be dominated by interactions resulting from plant chemistry in defense strategies, and evolution of secondary metabolisms being in response to insect herbivory remains a classic explanation of coevolution. The present study examines whether evolutionary constraints existing in host associations of economically important fruit flies in the species-rich tribe Dacini (Diptera: Tephritidae) and to what extent these species have evolved specialized dietary patterns. We found a strong effect of host phylogeny on associations on the 37 fruit flies tested, although the fruit fly species feeding on ripe commercially grown fruits that lost the toxic compounds after long-term domestication are mostly polyphagous. We assessed the phylogenetic signal of host breadth across the fruit fly species, showing that the results were substantially different depending on partition levels. Further, we mapped main host family associations onto the fruit fly phylogeny and Cucurbitaceae has been inferred as the most likely ancestral host family for Dacini based on ancestral state reconstruction.

Τhe Complete Mitochondrial Genome of Bactrocera carambolae (Diptera: Tephritidae): Genome Description and Phylogenetic Implications

Bactrocera carambolae is one of the approximately 100 sibling species of the Bactrocera dorsalis complex and considered to be very closely related to B. dorsalis. Due to their high morphological similarity and overlapping distribution, as well as to their economic impact and quarantine status, the development of reliable markers for species delimitation between the two taxa is of great importance. Here we present the complete mitochondrial genome of B. carambolae sourced from its native range in Malaysia and its invaded territory in Suriname. The mitogenome of B. carambolae presents the typical organization of an insect mitochondrion. Comparisons of the analyzed B. carambolae sequences to all available complete mitochondrial sequences of B. dorsalis revealed several species-specific polymorphic sites. Phylogenetic analysis based on Bactrocera mitogenomes supports that B. carambolae is a differentiated taxon though closely related to B. dorsalis. The present complete mitochondrial sequences of B. carambolae could be used, in the frame of Integrative Taxonomy, for species discrimination and resolution of the phylogenetic relationships within this taxonomically challenging complex, which would facilitate the application of species-specific population suppression strategies, such as the sterile insect technique.

Genetic characterization of Bactrocera fruit flies (Diptera: Tephritidae) from Northeastern India based on DNA barcodes

The Northeastern region of India, one of the mega biodiversity hot spots has enormous potential for the production of fruits and vegetables. Fruit flies of the genus Bactrocera Macquart are important pests of fruits and vegetables, and one of the limiting factors in successful production of these commodities. The relationship among some of the species is unclear due to their high molecular and morphological similarities. Moreover, due to the significant morphological resemblance between fruit fly species, reliable identification is very difficult task. We genetically characterized 10 fruit fly species of the genus Bactrocera by using standard DNA barcoding region of COI gene. The characterization and identification of eight species were straight forward. This study was unable to establish the molecular identity of Bactrocera sp. 2. Within the 547 bp region of partial COI gene, there were 157 variable sites of which 110 sites were parsimony informative, 153 were synonymous substitutions and 4 were non-synonymous substitutions. The estimate of genetic divergence among the ten species was in the range of 0-21.9% and the pairwise genetic distance of Bactrocera. (Bactrocera) dorsalis (Hendel) with B. (B.) carambolae was only 0.7%. Phylogenetic analysis formed separate clades for fruit and vegetable infesting fruit flies. B. (B.) aethriobasis Hardy, B. (B.) thailandica and B. (B.) tuberculata (Bezzi) have been reported for the first time from the Northeastern India. The information generated from this study would certainly have implications for pest management, taxonomy, quarantine and trade.

Genetic analysis of Bactrocera zonata (Diptera: Tephritidae) populations from India based on cox1 and nad1 gene sequences

The peach fruit fly, Bactrocera zonata, is among the most serious and polyphagous insect pest of fruit crops in many parts of the world under genus Bactrocera. In the present study, the genetic structure, diversity and demographic history of B. zonata in India were inferred from mitochondrial cytochrome oxidase 1 (cox1) and NADH dehydrogenase 1 (nad1) sequences. The efficiency of DNA barcodes for identification of B. zonata was also tested. Genetic diversity indices [number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π) and average number of nucleotide differences (k)] of B. zonata populations across India maintain high level of genetic diversity without isolation by distance among the geographic regions. Non-significant negative correlation between pairwise Fst and geographic distance suggests a high level of gene flow among studied populations of B. zonata. The possibility of sudden expansion of B. zonata revealed through mismatch distribution analysis as well as negative Tajima's D and Fu's Fs values further supported by star-like network of haplotypes. DNA barcoding analysis suggests that B. zonata specimens can be clearly differentiated from other species with 100% accuracy of identification. Therefore, cytochrome oxidase 1 (cox1) barcode sequences generated in the present study could be a valuable source for the rapid identification and global population genetic study of B. zonata.

The chromosomes and the mitogenome of Ceratitis fasciventris (Diptera: Tephritidae): two genetic approaches towards the Ceratitis FAR species complex resolution

Ceratitis fasciventris is a serious agricultural pest of the Tephritidae family that belongs to the African Ceratitis FAR species complex. Species limits within the FAR complex are obscure and multidisciplinary approaches have attempted to resolve phylogenetic relationships among its members. These studies support the existence of at least three additional species in the complex, C. anonnae, C. rosa and C. quilicii, while they indicate the presence of two structured populations (F1 and F2) within the C. fasciventris species. In the present study we present the mitotic karyotype, polytene chromosome maps, in situ hybridization data and the complete mitochondrial genome sequence of an F2 population of C. fasciventris. This is the first polytene chromosome map and complete mitogenome of a member of the FAR complex and only the second reported for the Ceratitis genus. Both polytene chromosomes and mitochondrial sequence could provide valuable information and be used as reference for comparative analysis among the members of the complex towards the clarification of their phylogenetic relationships.

Violet LED light enhances the recruitment of a thrip predator in open fields

The predatory bug Orius sauteri is an indigenous natural enemy of thrips and whiteflies in Asian countries. To put these bugs to practical use in pest management, methods to attract and retain the bugs in agricultural fields are needed. We previously showed that violet light (405 nm) attracts O. sauteri selectively. Many thrips and whiteflies are attracted to UV or green light. In this study, we examined the effect of violet-LED illumination on O. sauteri in pesticide-free eggplant (Solanum melongena L.) cultivation. In three cultivation trials, the density of O. sauteri on eggplant leaves was consistently higher in the illuminated plots; at least twice that of the non-illuminated plots. Simultaneously, the density of thrips declined markedly to less than half that of the non-illuminated plots. We identified three positive effects of violet light including an “immediate-effect” on predator attraction, a “persistent-effect” on predator reproduction, and a “secondary-effect” on the food web structure. Our results showed that illumination with violet light provides a powerful tool for integrated pest management. This is the first report on the use of illumination to manipulate the behavior of natural enemies.

Historical perspective on the synonymization of the four major pest species belonging to the Bactrocera dorsalis species complex (Diptera, Tephritidae)

An FAO/IAEA-sponsored coordinated research project on integrative taxonomy, involving close to 50 researchers from at least 20 countries, culminated in a significant breakthrough in the recognition that four major pest species, Bactroceradorsalis, Bactroceraphilippinensis, Bactrocerapapayae and Bactrocerainvadens, belong to the same biological species, Bactroceradorsalis. The successful conclusion of this initiative is expected to significantly facilitate global agricultural trade, primarily through the lifting of quarantine restrictions that have long affected many countries, especially those in regions such as Asia and Africa that have large potential for fresh fruit and vegetable commodity exports. This work stems from two taxonomic studies: a revision in 1994 that significantly increased the number of described species in the Bactroceradorsalis species complex; and the description in 2005 of Bactrocerainvadens, then newly incursive in Africa. While taxonomically valid species, many biologists considered that these were different names for one biological species. Many disagreements confounded attempts to develop a solution for resolving this taxonomic issue, before the FAO/IAEA project commenced. Crucial to understanding the success of that initiative is an accounting of the historical events and perspectives leading up to the international, multidisciplinary collaborative efforts that successfully achieved the final synonymization. This review highlights the 21 year journey taken to achieve this outcome.

The mitochondrial genome of the peach fruit fly, Bactrocera zonata (Saunders) (Diptera: Tephritidae): Complete DNA sequence, genome organization, and phylogenetic analysis with other tephritids using next generation DNA sequencing

Mitochondrial genome can provide information for genomic structure as well as for phylogenetic analysis and evolutionary biology. The complete 15,935 bp mitochondrial genome of Bactrocera zonata (Diptera: Tephritidae), is assembled from Illumina MiSeq read data. The mitogenome information for B. zonata was compared to the homologous sequences of other tephritids. Annotation indicated that the structure and orientation of 13 protein coding genes (PCGs), 22 tRNA and 2 rRNA sequences were typical of, and similar to, the ten closely related tephritid species. The nucleotide composition shows heavily biased toward As and Ts accounting 73.34% and exhibits a slightly positive AT skew, which is similar to other known tephritid species. All PCGs are initiated by ATN codons, except for cox1 with TCG and atp8 with GTG. Nine PCGs use a common stop codon of TAA or TAG, whereas the remaining four use an incomplete termination codon T or TA likely to be completed by adenylation. All tRNAs have the typical clover-leaf structure, with an exception for trnS((AGN)). Four short intergenic spacers showed high degree of conservation among B. zonata and other ten tephritids. A poly(T) stretch at the 5' end followed by [TA(A)]n-like stretch and a tandem repeats of 39 bp has been observed in CR. The analysis of gene evolutionary rate revealed that the cox1 and atp6 exhibits lowest and highest gene substitution rates, respectively than other genes. The phylogenetic relationships based on Maximum Likelihood method using all protein-coding genes and two ribosomal RNA genes confirmed that B. zonata is closely related to Bactrocera correcta, Bactrocera carambolae, Bactrocera papayae, and Bactrocera philippinensis and Bactrocera dorsalis belonging to B. dorsalis species complex forms a monophyletic clade, which is in accordance with the traditional morphological classification and recent molecular works.

The Bactrocera dorsalis species complex: comparative cytogenetic analysis in support of Sterile Insect Technique applications

Background

The Bactrocera dorsalis species complex currently harbors approximately 90 different members. The species complex has undergone many revisions in the past decades, and there is still an ongoing debate about the species limits. The availability of a variety of tools and approaches, such as molecular-genomic and cytogenetic analyses, are expected to shed light on the rather complicated issues of species complexes and incipient speciation. The clarification of genetic relationships among the different members of this complex is a prerequisite for the rational application of sterile insect technique (SIT) approaches for population control.

Results

Colonies established in the Insect Pest Control Laboratory (IPCL) (Seibersdorf, Vienna), representing five of the main economic important members of the Bactrocera dorsalis complex were cytologically characterized. The taxa under study were B. dorsalis s.s., B. philippinensis, B. papayae, B. invadens and B. carambolae. Mitotic and polytene chromosome analyses did not reveal any chromosomal characteristics that could be used to distinguish between the investigated members of the B. dorsalis complex. Therefore, their polytene chromosomes can be regarded as homosequential with the reference maps of B. dorsalis s.s.. In situ hybridization of six genes further supported the proposed homosequentiallity of the chromosomes of these specific members of the complex.

Conclusions

The present analysis supports that the polytene chromosomes of the five taxa under study are homosequential. Therefore, the use of the available polytene chromosome maps for B. dorsalis s.s. as reference maps for all these five biological entities is proposed. Present data provide important insight in the genetic relationships among the different members of the B. dorsalis complex, and, along with other studies in the field, can facilitate SIT applications targeting this complex. Moreover, the availability of 'universal' reference polytene chromosome maps for members of the complex, along with the documented application of in situ hybridization, can facilitate ongoing and future genome projects in this complex.

Rapid diagnosis of the economically important fruit fly, Bactrocera correcta (Diptera: Tephritidae) based on a species-specific barcoding cytochrome oxidase I marker

The guava fruit fly, Bactrocera correcta (Bezzi) (Diptera: Tephritidae), is an invasive pest of fruit and vegetable crops that primarily inhabits Southeast Asia and which has the potential to become a major threat within both the Oriental and Australian oceanic regions as well as California and Florida. In light of the threat posed, it is important to develop a rapid, accurate and reliable method to identify B. correcta in quarantine work in order to provide an early warning to prevent its widespread invasion. In the present study, we describe a species-specific polymerase chain reaction assay for the diagnosis of B. correcta using mitochondrial DNA cytochrome oxidase I (mtDNA COI) barcoding genes. A B. correcta-specific primer pair was designed according to variations in the mtDNA COI barcode sequences among 14 fruit fly species. The specificity and sensitivity of the B. correcta-specific primer pair was tested based on the presence or absence of a band in the gel profile. A pair of species-specific B. correcta primers was successfully designed and named BCOR-F/BCOR-R. An ∼280 bp fragment was amplified from specimens belonging to 17 geographical populations and four life stages of B. correcta, while no such diagnostic bands were present in any of the 14 other related fruit fly species examined. Sensitivity test results demonstrated that successful amplification can be obtained with as little as 1 ng μl⁻¹ of template DNA. The species-specific PCR analysis was able to successfully diagnose B. correcta, even in immature life stages, and from adult body parts. This method proved to be a robust single-step molecular technique for the diagnosis of B. correcta with respect to potential plant quarantine.

Taxonomic identity of the invasive fruit fly pest, Bactrocera invadens: concordance in morphometry and DNA barcoding

In 2003, a new fruit fly pest species was recorded for the first time in Kenya and has subsequently been found in 28 countries across tropical Africa. The insect was described as Bactrocera invadens, due to its rapid invasion of the African continent. In this study, the morphometry and DNA Barcoding of different populations of B. invadens distributed across the species range of tropical Africa and a sample from the pest's putative aboriginal home of Sri Lanka was investigated. Morphometry using wing veins and tibia length was used to separate B. invadens populations from other closely related Bactrocera species. The Principal component analysis yielded 15 components which correspond to the 15 morphometric measurements. The first two principal axes contributed to 90.7% of the total variance and showed partial separation of these populations. Canonical discriminant analysis indicated that only the first five canonical variates were statistically significant. The first two canonical variates contributed a total of 80.9% of the total variance clustering B. invadens with other members of the B. dorsalis complex while distinctly separating B. correcta, B. cucurbitae, B. oleae and B. zonata. The largest Mahalanobis squared distance (D² = 122.9) was found to be between B. cucurbitae and B. zonata, while the lowest was observed between B. invadens populations against B. kandiensis (8.1) and against B. dorsalis s.s (11.4). Evolutionary history inferred by the Neighbor-Joining method clustered the Bactrocera species populations into four clusters. First cluster consisted of the B. dorsalis complex (B. invadens, B. kandiensis and B. dorsalis s. s.), branching from the same node while the second group was paraphyletic clades of B. correcta and B. zonata. The last two are monophyletic clades, consisting of B. cucurbitae and B. oleae, respectively. Principal component analysis using the genetic distances confirmed the clustering inferred by the NJ tree.

Discovery of genes related to insecticide resistance in Bactrocera dorsalis by functional genomic analysis of a de novo assembled transcriptome

Insecticide resistance has recently become a critical concern for control of many insect pest species. Genome sequencing and global quantization of gene expression through analysis of the transcriptome can provide useful information relevant to this challenging problem. The oriental fruit fly, Bactrocera dorsalis, is one of the world's most destructive agricultural pests, and recently it has been used as a target for studies of genetic mechanisms related to insecticide resistance. However, prior to this study, the molecular data available for this species was largely limited to genes identified through homology. To provide a broader pool of gene sequences of potential interest with regard to insecticide resistance, this study uses whole transcriptome analysis developed through de novo assembly of short reads generated by next-generation sequencing (NGS). The transcriptome of B. dorsalis was initially constructed using Illumina's Solexa sequencing technology. Qualified reads were assembled into contigs and potential splicing variants (isotigs). A total of 29,067 isotigs have putative homologues in the non-redundant (nr) protein database from NCBI, and 11,073 of these correspond to distinct D. melanogaster proteins in the RefSeq database. Approximately 5,546 isotigs contain coding sequences that are at least 80% complete and appear to represent B. dorsalis genes. We observed a strong correlation between the completeness of the assembled sequences and the expression intensity of the transcripts. The assembled sequences were also used to identify large numbers of genes potentially belonging to families related to insecticide resistance. A total of 90 P450-, 42 GST-and 37 COE-related genes, representing three major enzyme families involved in insecticide metabolism and resistance, were identified. In addition, 36 isotigs were discovered to contain target site sequences related to four classes of resistance genes. Identified sequence motifs were also analyzed to characterize putative polypeptide translational products and associate them with specific genes and protein functions.

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

Distinct genetic lineages of Bactrocera caudata (Insecta: Tephritidae) revealed by COI and 16S DNA sequences

The fruit fly Bactrocera caudata is a pest species of economic importance in Asia. Its larvae feed on the flowers of Cucurbitaceae such as Cucurbita moschata. To-date it is distinguished from related species based on morphological characters. Specimens of B. caudata from Peninsular Malaysia and Indonesia (Bali and Lombok) were analysed using the partial DNA sequences of cytochrome c oxidase subunit I (COI) and 16S rRNA genes. Both gene sequences revealed that B. caudata from Peninsular Malaysia was distinctly different from B. caudata of Bali and Lombok, without common haplotype between them. Phylogenetic analysis revealed two distinct clades, indicating distinct genetic lineage. The uncorrected ‘p’ distance for COI sequences between B. caudata of Malaysia-Thailand-China and B. caudata of Bali-Lombok was 5.65%, for 16S sequences from 2.76 to 2.99%, and for combined COI and 16S sequences 4.45 to 4.46%. The ‘p’ values are distinctly different from intraspecific ‘p’ distance (0–0.23%). Both the B. caudata lineages are distinctly separated from related species in the subgenus Zeugodacus – B. ascita, B. scutellata, B. ishigakiensis, B. diaphora, B. tau, B. cucurbitae, and B. depressa. Molecular phylogenetic analysis indicates that the B. caudata lineages are closely related to B. ascita sp. B, and form a clade with B. scutellata, B. ishigakiensis, B. diaphora and B. ascita sp. A. This study provides additional baseline for the phylogenetic relationships of Bactrocera fruit flies of the subgenus Zeugodacus. Both the COI and 16S genes could be useful markers for the molecular differentiation and phylogenetic analysis of tephritid fruit flies.

Molecular phylogeny and identification of the peach fruit fly, Bactrocera zonata, established in Egypt

The genetic structure of the Egyptian peach fruit fly (Bactrocera zonata (Saunders) (Diptera: Tephritidae)) population was analyzed using total RNA from adult females. A portion of mitochondrial cytochrome oxidase I (COI), 369 bp was amplified using RT-PCR, and was sequenced and analyzed to clarify the phylogenetic relationship of B. zonata established in Egypt. The data suggested that the gene shared a similarity in sequence compared to Bactrocera COI gene found in GenBank. Molecular phylogenetic analyses were performed based on nucleotide sequences in order to examine the position of the Egyptian population among many other species of fruit flies. The results indicate that four accession numbers of B. zonata (three from New Zealand and one from India) are closely related, while the Egyptian B. zonata are close to the 71 accession numbers of Bactrocera include one B. zonata from New Zealand. These two B. zonata from Egypt and New Zealand showed a close relationship in neighbor-joining analysis using the seven accession numbers of B. zonata. In addition, a theoretical restriction map of the homology portion of the COI gene was constructed using 212 restriction enzymes obtained from the restriction enzyme database to identify the Egyptian and New Zealand B. zonata.

Development of allele-specific single-nucleotide polymorphism-based polymerase chain reaction markers in cytochrome oxidase I for the differentiation of Bactrocera papayae and Bactrocera carambolae (Diptera: Tephritidae)

Differentiation of Bactrocera papayae Drew & Hancock and Bactrocera carambolae Drew & Hancock (Diptera: Tephritidae) based on morphological characters has often been problematical. We describe here a single-nucleotide polymorphism (SNP)-based polymerase chain reaction (PCR) assay to differentiate between these two species. For detection of SNPs, fragments derived from each species were amplified using two primer pairs, COIF/COIR and UEA7/UEA10, sequenced, and aligned to obtain a contiguous 1,517-bp segment. Two new sets of primers were designed based on the 11 SNPs identified in the region. Results of the SNP-PCR test using any one of these species-specific primer sets indicate that these two species could be differentiated on basis of presence or absence of a band in the gel profile. We also tested the SNP-PCR primers on Bactrocera umbrosa F., Bactrocera cucurbitae Coquillett, Bactrocera latifrons Hendel, and Bactrocera tau (Walker) but did not detect any band in the gel, indicating the likelihood of a false positive for B. papayae is nil. This SNP-PCR method is efficient and useful, especially for immature life stages or when only adult body parts of the two species are available for identification, as encountered often in quarantine work.

Mitotic and polytene chromosomes analysis of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

The Oriental fruit fly, Batrocera dorsalis s.s. (Hendel) is one of the most destructive agricultural pests, belonging to a large group of difficult to distinguish morphologically species, referred as the B. dorsalis complex. We report here a cytogenetic analysis of two laboratory strains of the species and provide a photographic polytene chromosome map from larval salivary glands. The mitotic complement consists of six chromosome pairs including a heteromorphic sex (XX/XY) chromosome pair. Analysis of the polytene complement has shown a total of five polytene chromosomes (10 polytene arms) that correspond to the five autosomes. The most important landmarks of each polytene chromosome and characteristic asynapsis at a specific chromosomal region are presented and discussed. Chromosomal homology between B. dorsalis and Ceratitis capitata has been determined by comparing chromosome banding patterns. The detection of chromosome inversions in both B. dorsalis strains is shown and discussed. Our results show that the polytene maps presented here are suitable for cytogenetic analysis of this species and can be used for comparative studies among species of the Tephritidae family. They also provide a diagnostic tool that could accelerate species identification within the B. dorsalis complex and could shed light on the ongoing speciation in this complex. Polytene chromosome maps can facilitate the development of biological control methods and support the genome mapping project of the species that is currently in progress.

Species determination of Malaysian Bactrocera pests using PCR-RFLP analyses (Diptera: Tephritidae)

BACKGROUND

Identification of Bactrocera carambolae Drew and Hancock, B. papayae Drew and Hancock, B. tau Walker, B. latifrons Hendel, B. cucurbitae Coquillett, B. umbrosa Fabricius and B. caudata Fabricius would pose a problem if only a body part or an immature stage were available. Analysis of polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) of cytochrome oxidase I (COI) gene using primers COIR, COIF, UEA7 and UEA10 and restriction enzymes (MseI, RsaI and Alu1) was carried out. The banding profiles in the electrophoresis gel were analysed.

RESULTS

The COI gene in six Bactrocera spp. was successfully amplified by COIR and COIF, as well as UEA7 and UEA10, while B. caudata was amplified successfully only by UEA primers. Using COI amplified PCR products and restriction enzymes, distinct banding profiles for B. tau, B. latifrons, B. cucurbitae and B. umbrosa were observed, but not for B. carambolae and B. papayae. However, using UEA7, UEA10 and RsaI, B. caudata could be identified, while B. carambolae and B. papayae might possibly be separated from one another. It was also shown that adult body parts or immature life stages of B. carambolae, B. papayae, B. latifrons and B. cucurbitae produced the same banding profiles as the adults.

CONCLUSION

PCR-RFLP analyses are able to identify positively five Bactrocera species, while B. papayae and B. carambolae might possibly be separated from one another, even if immature life stages or adult body parts are used.

Highly similar piggyBac elements in Bactrocera that share a common lineage with elements in noctuid moths

The piggyBac IFP2 transposable element, originally discovered in a Trichoplusia ni cell line, also exists as nearly identical elements in other noctuid lepidopterans, and in several species of the tephritid genus Bactrocera. To further define the distribution of piggyBacs in Bactrocera, and compare their relationship to sequences found in Lepidoptera, a survey by PCR amplification was performed in a range of Bactrocera species. Highly similar piggyBac sequences were found in all B. dorsalis complex species tested, as well as in species in the B. zonata and B. frauenfeldi complexes. All nucleotide sequences had > 94% identity to corresponding sequences in the T. ni IFP2 element, and > 88% identity among the sequences. Conserved primers did not amplify any distantly related sequences that have been found by computational searches in a wider range of insect and non-insect species. Notably, 55 nucleotide substitutions relative to IFP2 were common to all the Bactrocera sequences, 44 of which exist in piggyBacs previously sequenced from moths, with 17 resulting in amino acid substitutions. These piggyBac elements, that apparently traversed orders by horizontal transfer, probably arose from a lineage separate from IFP2 and the other known elements in T. ni. Implications for the presence of nearly identical piggyBacs, in widely distributed insects, to the applied use of piggyBac vectors are discussed.

New contributions towards the understanding of the phylogenetic relationships among economically important fruit flies (Diptera: Tephritidae)

Fruit flies (Diptera: Tephritidae) are a species-rich and economically important group. The phylogenetic relationships among the many taxa are still to be fully resolved and the monophyly of several groups is still to be confirmed. This paper reports a study of the phylogenetic relationships among 23 economically important tephritid species (representing several major lineages of the family) which examines the sequence of a region of mitochondrial DNA encompassing the cytb, tRNA(Ser) and ND1 genes. Substitutions characteristic of particular taxa were found that could help classify members of the family at any developmental stage. The trees obtained by the maximum parsimony, neighbour joining and maximum likelihood methods were generally compatible with present morphological classification patterns. However, the data reveal some characteristics of the phylogenetic relationships of this family that do not agree with present classifications. The results support the probable non-monophyletic nature of the subfamily Trypetinae and suggest that Bactrocera cucurbitae (Coquillet) is more closely related to the genus Dacus than to other species of Bactrocera.

Molecular phylogenetics of the genus Ceratitis (Diptera: Tephritidae)

The Afrotropical fruit fly genus Ceratitis MacLeay is an economically important group that comprises over 89 species, subdivided into six subgenera. Cladistic analyses of morphological and host use characters have produced several phylogenetic hypotheses for the genus. Only monophyly of the subgenera Pardalaspis and Ceratitis (sensu stricto) and polyphyly of the subgenus Ceratalaspis are common to all of these phylogenies. In this study, the hypotheses developed from morphological and host use characters are tested using gene trees produced from DNA sequence data of two mitochondrial genes (cytochrome oxidase I and NADH-dehydrogenase subunit 6) and a nuclear gene (period). Comparison of gene trees indicates the following relationships: the subgenus Pardalaspis is monophyletic, subsection A of the subgenus Pterandrus is monophyletic, the subgenus Pterandrus may be either paraphyletic or polyphyletic, the subgenus Ceratalaspis is polyphyletic, and the subgenus Ceratitis s. s. might not be monophyletic. In addition, the genera Ceratitis and Trirhithrum do not form reciprocally monophyletic clades in the gene trees. Although the data statistically reject monophyly for Trirhithrum under the Shimodaira-Hasegawa test, they do not reject monophyly of Ceratitis.

Invasive phytophagous pests arising through a recent tropical evolutionary radiation: the Bactrocera dorsalis complex of fruit flies

The Bactrocera dorsalis complex of tropical fruit flies (Diptera: Tephritidae: Dacinae) contains 75 described species, largely endemic to Southeast Asia. Within the complex are a small number of polyphagous pests of international significance, including B. dorsalis sensu stricto, B. papayae, B. carambolae, and B. philippinensis. Most species within the complex were described in 1994 and since then substantial research has been undertaken in developing morphological and molecular diagnostic techniques for their recognition. Such techniques can now resolve most taxa adequately. Genetic evidence suggests that the complex has evolved in only the last few million years, and development of a phylogeny of the group is considered a high priority to provide a framework for future evolutionary and ecological studies. As model systems, mating studies on B. dorsalis s.s. and B. cacuminata have substantially advanced our understanding of insect use of plant-derived chemicals for mating, but such studies have not been applied to help resolve the limits of biological species within the complex. Although they are commonly regarded as major pests, there is little published evidence documenting economic losses caused by flies of the B. dorsalis complex. Quantification of economic losses caused by B. dorsalis complex species is urgently needed to prioritize research for quarantine and management. Although they have been documented as invaders, relatively little work has been done on the invasion biology of the complex and this is an area warranting further work.

Functional significance of phytochemical lures to dacine fruit flies (Diptera: Tephritidae): an ecological and evolutionary synthesis

While lures of plant origin are vital tools in dacine (Diptera: Tephritidae) pest management, the ecological and evolutionary significance of this lure response remains enigmatic. Two hypotheses (the ancestral host hypothesis and sexual selection by female choice) have been invoked to explain the functional significance of these chemicals to dacine fruit flies. These hypotheses are often treated as alternatives to one another and evidence favouring one is used to reject the other. This review highlights that these two hypotheses are not logical alternatives to each other as the ancestral host hypothesis attempts to explain the ultimate function of the response of Dacinae to these plant-derived parapheromones while the sexual selection hypothesis provides a proximate explanation for lure response. Research on lure response, dacine mating behaviour, functional significance of lures, plant phylogeny and biochemistry and dacine pheromone chemistry are used to evaluate the evidence in relation to both these hypotheses. Some of the key findings are that there is evidence both in support of and against these two hypotheses. Response of fruit flies to related phenyl propanoids to those commonly used as lures in pest management and distribution of phenyl propanoids attractive to dacines among plant orders strongly support the ancestral host hypothesis. Evidence from pheromone chemistry, dacine mating behaviour and the functional significance of lures both support and contradict the sexual selection hypothesis. Lures appear to have different proximate functions in different dacine species. Considerably greater research is needed to clarify the functional role of phytochemical lures to dacine fruit flies. The two prevalent hypotheses should be investigated independently. Specific research on dacine phylogeny and distribution of lures in plants in relation to ecological roles played by adult dacines is required to elucidate the ultimate roles of the chemicals. Exploration of female response to lures and the behavioural consequences of dacine response to these chemicals to both the insect and plant may shed light on the proximate functions of these chemicals.

Phylogenetic relationships among Bactrocera species (Diptera: Tephritidae) inferred from mitochondrial DNA sequences

Several members of the dipteran family Tephritdae are serious pests because females lay eggs in ripening fruit. The genus Bactrocera is one of the largest within the family with over 500 described species arranged in 28 subgenera. The phylogenetic relationships among the various species and subgenera, and the monophyly of specific groups have not been examined using a rigorous phylogenetic analysis. Therefore, phylogenetic relationships among 24 Bactrocera species belonging to 9 subgenera were inferred from DNA sequence of portions of the mitochondrial 16S rRNA, cytochrome oxidase II, tRNA(Lys), and tRNA(Asp) genes. Two morphological characters that traditionally have been used to define the four groups within the subgenus Bactrocera were evaluated in a phylogenetic context by mapping the character states onto the parsimony tree. In addition, the evolutionary trend in male-lure response was evaluated in a phylogenetic context. Maximum parsimony analyses suggested the following relationships: (1) the genus Bactrocera is monophyletic, (2) the subgenus B. (Zeugodacus) is paraphyletic, (3) the subgenus B. (Daculus) is a sister group to subgenus B. (Bactrocera), and (4) the subgenus B. (Bactrocera) is monophyletic. The mapping analyses suggested that the morphological characters exhibit a simple evolutionary transition from one character state to another. Male-lure response was identified as being a labile behavior that has been lost on multiple occasions. Cue-lure response was plesiomorphic to methyl-eugenol response, and the latter has evolved independently within the Bactrocera and Zeugodacus groups of subgenera. The implications of our results for devising a coherent, consolidated classification for Bactrocera is discussed.