The draft genome of the pest tephritid fruit fly Bactrocera tryoni: resources for the genomic analysis of hybridising species

Evolutionary History of the DD41D Family of Tc1/Mariner Transposons in Two Mayetiola Species

Tc1/mariner elements are ubiquitous in eukaryotic genomes including insects. They are diverse and divided into families and sub-families. The DD34D family including mauritiana and irritans subfamilies have already been identified in two closely related species of Cecidomyiids M. destructor and M. hordei. In the current study the de novo and similarity-based methods allowed the identification for the first time of seven consensuses in M. destructor and two consensuses in M. hordei belonging to DD41D family whereas the in vitro method allowed the amplification of two and three elements in these two species respectively. Most of identified elements accumulated different mutations and long deletions spanning the N-terminal region of the transposase. Phylogenetic analyses showed that the DD41D elements were clustered in two groups belonging to rosa and Long-TIR subfamilies. The age estimation of the last transposition events of the identified Tc1/mariner elements in M. destructor showed different evolutionary histories. Indeed, irritans elements have oscillated between periods of silencing and reappearance while rosa and mauritiana elements have shown regular activity with large recent bursts. The study of insertion sites showed that they are mostly intronic and that some recently transposed elements occurred in genes linked to putative DNA-binding domains and enzymes involved in metabolic chains. Thus, this study gave evidence of the existence of DD41D family in two Mayetiola species and an insight on their evolutionary history.

No single rescue recipe: genome complexities modulate insect response to climate change

Declines in insect populations have gained formidable attention. Given their crucial role in the ecosystem, the causes of declining insect populations must be investigated. However, the insect clade has been associated with low extinction and high diversification rates. It is unlikely that insects underwent mass extinctions in the past. However, the pace of current climate change could make insect populations vulnerable to extinction. We propose genome size (GS) and transposable elements (TEs) to be rough estimates to assess extinction risk. Larger GS and/or proliferating TEs have been associated with adaptation in rapid climate change scenarios. We speculate that unstable, stressful environmental conditions are strongly associated with GS and TE expansion, which could be further correlated with adaptations. Alternately, stressful conditions trigger TE bursts that are not purged in smaller populations. GS and TE loads could be indicators of small effective populations in the wild, likely experiencing bottlenecks or drastic climatic perturbations, which calls for an urgent assessment of extinction risk.

Conservation of shibire and RpII215 temperature-sensitive lethal mutations between Drosophila and Bactrocera tryoni

The sterile insect technique can suppress and eliminate population outbreaks of the Australian horticultural pest, Bactrocera tryoni, the Queensland fruit fly. Sterile males mate with wild females that produce inviable embryos, causing population suppression or elimination. Current sterile insect releases are mixed sex, as the efficient removal of unrequired factory-reared females is not yet possible. In this paper, we assessed the known Drosophila melanogaster temperature-sensitive embryonic lethal alleles shibire (G268D, shits1) and RNA polymerase II 215 (R977C, RpII215ts) for potential use in developing B. tryoni genetic sexing strains (GSS) for the conditional removal of females. Complementation tests in D. melanogaster wild-type or temperature-sensitive genetic backgrounds were performed using the GAL4-UAS transgene expression system. A B. tryoni wild-type shibire isoform partially rescued Drosophila temperature lethality at 29°C by improving survivorship to pupation, while expressing B. tryoni shits1 failed to rescue the lethality, supporting a temperature-sensitive phenotype. Expression of the B. tryoni RpII215 wild-type protein rescued the lethality of D. melanogaster RpII215ts flies at 29°C. Overexpressing the B. tryoni RpII215ts allele in the D. melanogaster wild-type background unexpectedly produced a dominant lethal phenotype at 29°C. The B. tryoni shibire and RpII215 wild-type alleles were able to compensate, to varying degrees, for the function of the D. melanogaster temperature-sensitive proteins, supporting functional conservation across species. Shibire and RpII215 hold potential for developing insect strains that can selectively kill using elevated temperatures; however, alleles with milder effects than shits1 will need to be considered.

Patterns of Variation in the Usage of Fatty Acid Chains among Classes of Ester and Ether Neutral Lipids and Phospholipids in the Queensland Fruit Fly

Modern lipidomics has the power and sensitivity to elucidate the role of insects' lipidomes in their adaptations to the environment at a mechanistic molecular level. However, few lipidomic studies have yet been conducted on insects beyond model species such as Drosophila melanogaster. Here, we present the lipidome of adult males of another higher dipteran frugivore, Bactrocera tryoni. We describe 421 lipids across 15 classes of ester neutral lipids and phospholipids and ether neutral lipids and phospholipids. Most of the lipids are specified in terms of the carbon and double bond contents of each constituent hydrocarbon chain, and more ether lipids are specified to this degree than in any previous insect lipidomic analyses. Class-specific profiles of chain length and (un)saturation are broadly similar to those reported in D. melanogaster, although we found fewer medium-length chains in ether lipids. The high level of chain specification in our dataset also revealed widespread non-random combinations of different chain types in several ester lipid classes, including deficits of combinations involving chains of the same carbon and double bond contents among four phospholipid classes and excesses of combinations of dissimilar chains in several classes. Large differences were also found in the length and double bond profiles of the acyl vs. alkyl or alkenyl chains of the ether lipids. Work on other organisms suggests some of the differences observed will be functionally consequential and mediated, at least in part, by differences in substrate specificity among enzymes in lipid synthesis and remodelling pathways. Interrogation of the B. tryoni genome showed it has comparable levels of diversity overall in these enzymes but with some gene gain/loss differences and considerable sequence divergence from D. melanogaster.

Chromosome-level genome assembly of an important wolfberry fruit fly (Neoceratitis asiatica Becker)

Tephritidae pests are significant agricultural pests with a notable impact on the economy, with a wide range of species and most having broad host ranges and strong reproductive abilities. However, the wolfberry fruit fly, Neoceratitis asiatica (Becker), is a Tephritidae fly which only harms wolfberry. Here, we assembled and annotated N. asiatica genome at the chromosome level and compared it with the genomic and transcriptomic information from other Tephritidae flies. The assembled genome of N. asiatica had a size of 563.8 Mb and achieved a completeness level of 99.1%, 18,387 genes were annotated totally. All contigs were assembled into 7 linkage groups with an N50 of 93.166 Mb assisted by the Hi-C technique. The high-quality genome developed here will provide a significant resource for exploring the genetic basis of the adaptive and reproductive differences among various Tephritidae pests, and provides an important theoretical basis for the prevention and control of Tephritidae pests.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

Genetic variation for rectal gland volatiles among recently collected isofemale lines and a domesticated strain of Queensland fruit fly, Bactrocera tryoni (Diptera: Tephritidae)

Divergence between populations in mating behaviour can function as a potent premating isolating mechanism and promote speciation. However, very few cases of inherited intraspecific variation in sexual signalling have been reported in tephritid fruit flies, despite them being a highly speciose family. We tested for such variation in one tephritid, the Queensland fruit fly, Bactrocera tryoni (Qfly). Qfly mating behaviour depends on volatiles secreted from male rectal glands but no role for the volatiles from female rectal glands has yet been reported. We previously detected over 100 volatile compounds in male rectal glands and identified over 30 of them. Similar numbers were recorded in females. However, many compounds showed presence/absence differences between the sexes and many others showed quantitative differences between them. Here we report inherited variation among 24 Qfly lines (23 isofemale lines established from recent field collections and one domesticated line) in the abundance of three esters, two alcohols, two amides, an aldehyde and 18 unidentified volatiles in male rectal glands. We did not find any compounds in female rectal glands that varied significantly among the lines, although this may at least partly reflect lower female sample numbers. Most of the 26 male compounds that differed between lines were more abundant in the domesticated line than any of the recently established isofemale lines, which concurs with other evidence for changes in mating behaviour during domestication of this species. There were also large differences in several of the 26 compounds among the isofemale lines, and some of these differences were associated with the regions from which the lines were collected. While some of the variation in different compounds was correlated across lines, much of it was not, implicating involvement of multiple genes. Our findings parallel reports of geographic variation in other Qfly traits and point to inherited differences in reproductive physiology that could provide a basis for evolution of premating isolation between ecotypes.

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.

Partial Ribosomal Nontranscribed Spacer Sequences Distinguish Rhagoletis zephyria (Diptera: Tephritidae) From the Apple Maggot, R. pomonella

The apple maggot, Rhagoletis pomonella (Walsh), was introduced into the apple-growing regions of the Pacific Northwest in the U.S.A. during the past 60-100 yr. Apple maggot (larvae, puparia, and adults) is difficult to distinguish from its morphologically similar sister species, Rhagoletis zephyria Snow, which is native and abundant in the Pacific Northwest. While morphological identifications are common practice, a simple, inexpensive assay based on genetic differences would be very useful when morphological traits are unclear. Here we report nucleotide substitution and insertion-deletion mutations in the nontranscribed spacer (NTS) of the ribosomal RNA gene cistron of R. pomonella and R. zephyria that appear to be diagnostic for these two fly species. Insertion-deletion variation is substantial and results in a 49 base-pair difference in PCR amplicon size between R. zephyria and R. pomonella that can be scored using agarose gel electrophoresis. PCR amplification and DNA sequencing of 766 bp of the NTS region from 38 R. pomonella individuals and 35 R. zephyria individuals from across their geographic ranges led to the expected PCR fragments of approx. 840 bp and 790 bp, respectively, as did amplification and sequencing of a smaller set of 26 R. pomonella and 16 R. zephyria flies from a sympatric site in Washington State. Conversely, 633 bp mitochondrial COI barcode sequences from this set of flies were polyphyletic with respect to R. pomonella and R. zephyria. Thus, differences in NTS PCR products on agarose gels potentially provide a simple way to distinguish between R. pomonella and R. zephyria.

Diversity and Evolution of pogo and Tc1/mariner Transposons in the Apoidea Genomes

Bees (Apoidea), the largest and most crucial radiation of pollinators, play a vital role in the ecosystem balance. Transposons are widely distributed in nature and are important drivers of species diversity. However, transposons are rarely reported in important pollinators such as bees. Here, we surveyed 37 bee genomesin Apoidea, annotated the pogo and Tc1/mariner transposons in the genome of each species, and performed a phylogenetic analysis and determined their overall distribution. The pogo and Tc1/mariner families showed high diversity and low abundance in the 37 species, and their proportion was significantly higher in solitary bees than in social bees. DD34D/mariner was found to be distributed in almost all species and was found in Apis mellifera, Apis mellifera carnica, Apis mellifera caucasia, and Apis mellifera mellifera, and Euglossa dilemma may still be active. Using horizontal transfer analysis, we found that DD29-30D/Tigger may have experienced horizontal transfer (HT) events. The current study displayed the evolution profiles (including diversity, activity, and abundance) of the pogo and Tc1/mariner transposons across 37 species of Apoidea. Our data revealed their contributions to the genomic variations across these species and facilitated in understanding of the genome evolution of this lineage.

Raspberry ketone diet supplement reduces attraction of sterile male Queensland fruit fly to cuelure by altering expression of chemoreceptor genes

Sterile male Queensland fruit fly, Bactrocera tryoni (Froggatt), fed as immature adults on the plant compound raspberry ketone (RK), show a reduced attraction to cuelure, a synthetic analogue of RK used as an attractant in Male Annihilation Technique. We hypothesized the reduced attraction of RK-fed adult males to cuelure may be a consequence of altered expression of chemoreception genes. A Y-tube olfactometer assay with RK-fed and RK-unfed sterile B. tryoni males tested the subsequent behavioural response to cuelure. Behavioral assays confirmed a significant decrease in attraction of RK-fed sterile males to cuelure. RK-fed, non-responders (to cue-lure) and RK-unfed, responders (to cue-lure) males were sampled and gene expression compared by de novo RNA-seq analysis. A total of 269 genes in fly heads were differentially expressed between replicated groups of RK-fed, cuelure non-responders and RK-unfed, cuelure responders. Among them, 218 genes including 4 chemoreceptor genes were up regulated and 51 genes were down regulated in RK-fed, cuelure non-responders. De novo assembly generated many genes with unknown functions and no significant BLAST hits to homologues in other species. The enriched and suppressed genes reported here, shed light on the transcriptional changes that affect the dynamics of insect responses to chemical stimuli.

Gene selection for studying frugivore-plant interactions: a review and an example using Queensland fruit fly in tomato

Fruit production is negatively affected by a wide range of frugivorous insects, among them tephritid fruit flies are one of the most important. As a replacement for pesticide-based controls, enhancing natural fruit resistance through biotechnology approaches is a poorly researched but promising alternative. The use of quantitative reverse transcription PCR (RT-qPCR) is an approach to studying gene expression which has been widely used in studying plant resistance to pathogens and non-frugivorous insect herbivores, and offers a starting point for fruit fly studies. In this paper, we develop a gene selection pipe-line for known induced-defense genes in tomato fruit, Solanum lycopersicum, and putative detoxification genes in Queensland fruit fly, Bactrocera tryoni, as a basis for future RT-qPCR research. The pipeline started with a literature review on plant/herbivore and plant/pathogen molecular interactions. With respect to the fly, this was then followed by the identification of gene families known to be associated with insect resistance to toxins, and then individual genes through reference to annotated B. tryoni transcriptomes and gene identity matching with related species. In contrast for tomato, a much better studied species, individual defense genes could be identified directly through literature research. For B. tryoni, gene selection was then further refined through gene expression studies. Ultimately 28 putative detoxification genes from cytochrome P450 (P450), carboxylesterase (CarE), glutathione S-transferases (GST), and ATP binding cassette transporters (ABC) gene families were identified for B. tryoni, and 15 induced defense genes from receptor-like kinase (RLK), D-mannose/L-galactose, mitogen-activated protein kinase (MAPK), lipoxygenase (LOX), gamma-aminobutyric acid (GABA) pathways and polyphenol oxidase (PPO), proteinase inhibitors (PI) and resistance (R) gene families were identified from tomato fruit. The developed gene selection process for B. tryoni can be applied to other herbivorous and frugivorous insect pests so long as the minimum necessary genomic information, an annotated transcriptome, is available.

Tephritid fruit flies have a large diversity of co-occurring RNA viruses

Tephritid fruit flies are amongst the most devastating pests of horticulture, and Sterile Insect Technique (SIT) programs have been developed for their control. Their interactions with viruses are still mostly unexplored, yet, viruses may negatively affect tephritid health and performance in SIT programs, and, conversely, constitute potential biological control agents. Here we analysed ten transcriptome libraries obtained from laboratory populations of nine tephritid species from Australia (six species of Bactrocera, and Zeugodacus cucumis), Asia (Bactrocera dorsalis) and Europe (Ceratitis capitata). We detected new viral diversity, including near-complete (>99%) and partially complete (>80%) genomes of 34 putative viruses belonging to eight RNA virus families. On average, transcriptome libraries included 3.7 viruses, ranging from 0 (Z. cucumis) to 9 (B. dorsalis). Most viruses belonged to the Picornavirales, represented by fourteen Dicistroviridae (DV), nine Iflaviridae (IV) and two picorna-like viruses. Others were a virus from Rhabdoviridae (RV), one from Xinmoviridae (both Mononegavirales), several unclassified Negev- and toti-like viruses, and one from Metaviridae (Ortervirales). Using diagnostic PCR primers for four viruses found in the transcriptome of the Bactrocera tryoni strain bent wings (BtDV1, BtDV2, BtIV1, and BtRV1), we tested nine Australian laboratory populations of five species (B. tryoni, Bactrocera neohumeralis, Bactrocera jarvisi, Bactrocera cacuminata, C. capitata), and one field population each of B. tryoni, B. cacuminata and Dirioxa pornia. Viruses were present in most laboratory and field populations yet their incidence differed for each virus. Prevalence and co-occurrence of viruses in B. tryoni and B. cacuminata were higher in laboratory than field populations. This raises concerns about the potential accumulation of viruses and their potential health effects in laboratory and mass-rearing environments which might affect flies used in research and control programs such as SIT.

Traveler, a New DD35E Family of Tc1/Mariner Transposons, Invaded Vertebrates Very Recently

The discovery of new members of the Tc1/mariner superfamily of transposons is expected based on the increasing availability of genome sequencing data. Here, we identified a new DD35E family termed Traveler (TR). Phylogenetic analyses of its DDE domain and full-length transposase showed that, although TR formed a monophyletic clade, it exhibited the highest sequence identity and closest phylogenetic relationship with DD34E/Tc1. This family displayed a very restricted taxonomic distribution in the animal kingdom and was only detected in ray-finned fish, anura, and squamata, including 91 vertebrate species. The structural organization of TRs was highly conserved across different classes of animals. Most intact TR transposons had a length of ∼1.5 kb (range 1,072-2,191 bp) and harbored a single open reading frame encoding a transposase of ∼340 aa (range 304-350 aa) flanked by two short-terminal inverted repeats (13-68 bp). Several conserved motifs, including two helix-turn-helix motifs, a GRPR motif, a nuclear localization sequence, and a DDE domain, were also identified in TR transposases. This study also demonstrated the presence of horizontal transfer events of TRs in vertebrates, whereas the average sequence identities and the evolutionary dynamics of TR elements across species and clusters strongly indicated that the TR family invaded the vertebrate lineage very recently and that some of these elements may be currently active, combining the intact TR copies in multiple lineages of vertebrates. These data will contribute to the understanding of the evolutionary history of Tc1/mariner transposons and that of their hosts.

Geographic distribution of sex chromosome polymorphism in Anastrepha fraterculus sp. 1 from Argentina

BACKGROUND

Anastrepha fraterculus is recognized as a quarantine pest in several American countries. This fruit fly species is native to the American continent and distributed throughout tropical and subtropical regions. It has been reported as a complex of cryptic species, and at least eight morphotypes have been described. Only one entity of this complex, formerly named Anastrepha fraterculus sp. 1, is present in Argentina. Previous cytogenetic studies on this morphotype described the presence of sex chromosome variation identified by chromosomal size and staining patterns. In this work, we expanded the cytological study of this morphotype by analyzing laboratory strains and wild populations to provide information about the frequency and geographic distribution of these sex chromosome variants. We analyzed the mitotic metaphases of individuals from four laboratory strains and five wild populations from the main fruit-producing areas of Argentina, including the northwest (Tucumán and La Rioja), northeast (Entre Ríos and Misiones), and center (Buenos Aires) of the country.

RESULTS

In wild samples, we observed a high frequency of X1X1 (0.94) and X1Y5 (0.93) karyomorphs, whereas X1X2 and X1Y6 were exclusively found at a low frequency in Buenos Aires (0.07 and 0.13, respectively), Entre Ríos (0.16 and 0.14, respectively) and Tucumán (0.03 and 0.04, respectively). X2X2 and X2Y5 karyomorphs were not found in wild populations but were detected at a low frequency in laboratory strains. In fact, karyomorph frequencies differed between wild populations and laboratory strains. No significant differences among A. fraterculus wild populations were evidenced in either karyotypic or chromosomal frequencies. However, a significant correlation was observed between Y5 chromosomal frequency and latitude.

CONCLUSIONS

We discuss the importance of cytogenetics to understand the possible route of invasion and dispersion of this pest in Argentina and the evolutionary forces acting under laboratory conditions, possibly driving changes in the chromosomal frequencies. Our findings provide deep and integral genetic knowledge of this species, which has become of relevance to the characterization and selection of valuable A. fraterculus sp. 1 strains for mass rearing production and SIT implementation.

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication

Background

The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a “common garden” approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication.

Results

Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait.

Conclusions

Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

Precise single base substitution in the shibire gene by CRISPR/Cas9-mediated homology directed repair in Bactrocera tryoni

Background

Pest eradication using the Sterile Insect Technique (SIT) involves high-density releases of sterilized males that mate with wild females and ultimately suppress the population. Sterilized females are not required for SIT and their removal or separation from males prior to release remains challenging. In order to develop genetic sexing strains (GSS), conditional traits such as temperature sensitive lethality are required.

Results

Here we introduce a known Drosophila melanogaster temperature sensitive embryonic lethal mutation into Bactrocera tryoni, a serious horticultural pest in Australia. A non-synonymous point mutation in the D. melanogaster gene shibire causes embryonic lethality at 29 °C and we successfully used CRISPR/Cas9 technology to recreate the orthologous shibire temperature sensitive-1 (shi^ts1) mutation in B. tryoni. Genotypic analyses over three generations revealed that a high fitness cost was associated with the shi^ts1 mutant allele and shi^ts1 homozygotes were not viable at 21 °C, which is a more severe phenotype than that documented in D. melanogaster.

Conclusions

We have demonstrated the first successful use of CRISPR/Cas9 to introduce precise single base substitutions in an endogenous gene via homology-directed repair in an agricultural pest insect and this technology can be used to trial other conditional mutations for the ultimate aim of generating genetic sexing strains for SIT.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00934-3.

Separating two tightly linked species-defining phenotypes in Bactrocera with hybrid recombinant analysis

Background

Bactrocera tryoni and Bactrocera neohumeralis mate asynchronously; the former mates exclusively around dusk while the latter mates during the day. The two species also differ in the colour of the post-pronotal lobe (callus), which is predominantly yellow in B. tryoni and brown in B. neohumeralis. We have examined the genetic relationship between the two characters in hybrids, backcrosses and multigeneration hybrid progeny.

Results

Our analysis of the mating time of the parental species revealed that while B. tryoni mate exclusively at dusk, B. neohumeralis females pair with B. neohumeralis males during the day and with B. tryoni males at dusk. We found considerable variance in mating time and callus colour among hybrid backcross individuals of both sexes but there was a strong although not invariant trend for callus colour to co-segregate with mating time in both sexes. To genetically separate these two phenotypes we allowed the interspecific F1 hybrids to propagate for 25 generations (F25) without selection for mating time or callus colour, finding that the advanced hybrid population had moved towards B. tryoni phenotypes for both traits. Selection for day mating in replicate lines at F25 resulted in significant phenotypic shifts in both traits towards B. neohumeralis phenotypes in F26. However, we were unable to completely recover the mating time profile of B. neohumeralis and relaxation of selection for day mating led to a shift back towards dusk mating, but not yellow callus colour, by F35.

Conclusion

We conclude that the inheritance of the two major species-defining traits is separable but tightly linked and involves more than one gene in each case. It also appears that laboratory conditions select for the B. tryoni phenotypes for mating time. We discuss our findings in relation to speciation theory and the likely effects of domestication during the generation of mass release strains for sterile insect control programmes.

Wolbachia pipientis Associated With Tephritid Fruit Fly Pests: From Basic Research to Applications

Members of the true fruit flies (family Tephritidae) are among the most serious agricultural pests worldwide, whose control and management demands large and costly international efforts. The need for cost-effective and environmentally friendly integrated pest management (IPM) has led to the development and implementation of autocidal control strategies. These approaches include the widely used sterile insect technique and the incompatible insect technique (IIT). IIT relies on maternally transmitted bacteria (namely Wolbachia) to cause a conditional sterility in crosses between released mass-reared Wolbachia-infected males and wild females, which are either uninfected or infected with a different Wolbachia strain (i.e., cytoplasmic incompatibility; CI). Herein, we review the current state of knowledge on Wolbachia-tephritid interactions including infection prevalence in wild populations, phenotypic consequences, and their impact on life history traits. Numerous pest tephritid species are reported to harbor Wolbachia infections, with a subset exhibiting high prevalence. The phenotypic effects of Wolbachia have been assessed in very few tephritid species, due in part to the difficulty of manipulating Wolbachia infection (removal or transinfection). Based on recent methodological advances (high-throughput DNA sequencing) and breakthroughs concerning the mechanistic basis of CI, we suggest research avenues that could accelerate generation of necessary knowledge for the potential use of Wolbachia-based IIT in area-wide integrated pest management (AW-IPM) strategies for the population control of tephritid pests.

Development and characterization of genomic resources for a non-model marine teleost, the red snapper (Lutjanus campechanus, Lutjanidae): Construction of a high-density linkage map, anchoring of genome contigs and comparative genomic analysis

The red snapper Lutjanus campechanus is an exploited reef fish of major economic importance in the Gulf of Mexico region. Studies of genome wide genetic variation are needed to understand the structure of wild populations and develop breeding programs for aquaculture but interpretation of these genome scans is limited by the absence of reference genome. In this work, the first draft of a reference genome was developed and characterized for the red snapper. P-454 and Illumina sequencing were conducted to produce paired-end reads that were assembled into reference contigs and scaffolds. The current assembly spans over 770 Mb, representing an estimated 69% of the red snapper genome in 67,254 scaffolds (N50 = 16,803 bp). The genome contigs were applied to map double digest Restriction-Site Associated DNA Tags and characterize Single Nucleotide Polymorphisms (SNPs) in five outbred full-sib families. The identified SNPs and 97 microsatellite loci were used to generate a high-density linkage map that includes 7,420 markers distributed across 24 linkage groups and spans 1,346.64 cM with an average inter–marker distance of 0.18 cM. Sex-specific maps revealed a 1.10:1 female to male map length ratio. A total of 4,422 genome contigs (10.5% of the assembly) were anchored to the map and used in a comparative genomic analysis of the red snapper and two model teleosts. Red snapper showed a high degree of chromosome level syntenic conservation with both medaka and spotted green puffer and a near one to one correspondence between the 24 red snapper linkage groups and corresponding medaka chromosomes was observed. This work established the first draft of a reference genome for a lutjanid fish. The obtained genomic resources will serve as a framework for the interpretation of genome scans during studies of wild populations and captive breeding programs.

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.

Identification of Y-chromosome scaffolds of the Queensland fruit fly reveals a duplicated gyf gene paralogue common to many Bactrocera pest species

Bactrocera tryoni (Queensland fruit fly) are polyphagous horticultural pests of eastern Australia. Heterogametic males contain a sex-determining Y-chromosome thought to be gene poor and repetitive. Here, we report 39 Y-chromosome scaffolds (~700 kb) from B. tryoni identified using genotype-by-sequencing data and whole-genome resequencing. Male diagnostic PCR assays validated eight Y-scaffolds, and one (Btry4096) contained a novel gene with five exons that encode a predicted 575 amino acid protein. The Y-gene, referred to as typo-gyf, is a truncated Y-chromosome paralogue of X-chromosome gene gyf (1773 aa). The Y-chromosome contained ~41 copies of typo-gyf, and expression occurred in male flies and embryos. Analysis of 13 tephritid transcriptomes confirmed typo-gyf expression in six additional Bactrocera species, including Bactrocera latifrons, Bactrocera dorsalis and Bactrocera zonata. Molecular dating estimated typo-gyf evolved within the past 8.02 million years (95% highest posterior density 10.56-5.52 million years), after the split with Bactrocera oleae. Phylogenetic analysis also highlighted complex evolutionary histories among several Bactrocera species, as discordant nuclear (116 genes) and mitochondrial (13 genes) topologies were observed. B. tryoni Y-sequences may provide useful sites for future transgene insertions, and typo-gyf could act as a Y-chromosome diagnostic marker for many Bactrocera species, although its function is unknown.

Differences in gene regulation in a tephritid model of prezygotic reproductive isolation

The two tephritid fruit fly pests, Bactrocera tryoni and Bactrocera neohumeralis, are unusually well suited to the study of the genetics of reproductive isolating mechanisms. Sequence difference between the species is no greater than between a pair of conspecific Drosophila melanogaster populations. The two species exist in close sympatry, yet do not hybridize in the field, apparently kept separate by a strong premating isolation mechanism involving the time of day at which mating occurs. This spurred us to search for key genes for which time of day expression is regulated differently between the species. Using replicated, quantitative transcriptomes from head tissues of males of the two species, sampled in the day and night, we identified 141 transcripts whose abundance showed a significant interaction between species and time of day, indicating a difference in gene regulation. The brain transcripts showing this interaction were enriched for genes with a neurone function and 90% of these were more abundant at night than day in B. tryoni. Features of the expression patterns suggest that there may be a difference in the regulation of sleep-wake cycles between the species. In particular several genes, which in D. melanogaster are expressed in circadian pacemaker cells, are promising candidates to further explore the genetic differentiation involved in this prezygotic reproductive isolation mechanism.

A genome-wide approach for uncovering evolutionary relationships of Australian Bactrocera species complexes (Diptera: Tephritidae)

Australia and Southeast Asia are hotspots of global diversity in the fruit-fly genus Bactrocera. Although a great diversity of species has been long recognised, evolutionary relationships are poorly understood, largely because previous sequencing techniques have provided insufficient historical signal for phylogenetic reconstruction. Poorly understood biogeographic history in Bactrocera has prevented a deeper understanding of migratory patterns in this economically important pest group. Using representatives from Australia and Malaysia, we tested the utility of a genome-reduction approach that generates thousands of single-nucleotide polymorphisms for phylogenetic reconstructions. This approach has high utility for species identification because of the ease of sample addition over time, and the species-level specificity able to be achieved with the markers. These data have provided a strongly supported phylogenetic tree congruent with topologies generated using more intensive sequencing approaches. In addition, our results do not support taxonomic assignments to species complex for a number of species, such as B. endiandrae in the dorsalis complex, yet find a close relationship between B. pallida and the dorsalis species. Our data have further validated non-monophyletic evolution of male response to primary attractants. We also showed at least two diversification events between Australia and Southeast Asia, indicating trans-regional dispersal in important pest species.

Near full-length 16S rRNA gene next-generation sequencing revealed Asaia as a common midgut bacterium of wild and domesticated Queensland fruit fly larvae

BACKGROUND

Gut microbiota affects tephritid (Diptera: Tephritidae) fruit fly development, physiology, behavior, and thus the quality of flies mass-reared for the sterile insect technique (SIT), a target-specific, sustainable, environmentally benign form of pest management. The Queensland fruit fly, Bactrocera tryoni (Tephritidae), is a significant horticultural pest in Australia and can be managed with SIT. Little is known about the impacts that laboratory-adaptation (domestication) and mass-rearing have on the tephritid larval gut microbiome. Read lengths of previous fruit fly next-generation sequencing (NGS) studies have limited the resolution of microbiome studies, and the diversity within populations is often overlooked. In this study, we used a new near full-length (> 1300 nt) 16S rRNA gene amplicon NGS approach to characterize gut bacterial communities of individual B. tryoni larvae from two field populations (developing in peaches) and three domesticated populations (mass- or laboratory-reared on artificial diets).

RESULTS

Near full-length 16S rRNA gene sequences were obtained for 56 B. tryoni larvae. OTU clustering at 99% similarity revealed that gut bacterial diversity was low and significantly lower in domesticated larvae. Bacteria commonly associated with fruit (Acetobacteraceae, Enterobacteriaceae, and Leuconostocaceae) were detected in wild larvae, but were largely absent from domesticated larvae. However, Asaia, an acetic acid bacterium not frequently detected within adult tephritid species, was detected in larvae of both wild and domesticated populations (55 out of 56 larval gut samples). Larvae from the same single peach shared a similar gut bacterial profile, whereas larvae from different peaches collected from the same tree had different gut bacterial profiles. Clustering of the Asaia near full-length sequences at 100% similarity showed that the wild flies from different locations had different Asaia strains.

CONCLUSIONS

Variation in the gut bacterial communities of B. tryoni larvae depends on diet, domestication, and horizontal acquisition. Bacterial variation in wild larvae suggests that more than one bacterial species can perform the same functional role; however, Asaia could be an important gut bacterium in larvae and warrants further study. A greater understanding of the functions of the bacteria detected in larvae could lead to increased fly quality and performance as part of the SIT.

Chromatin immunoprecipitation (ChIP) method for non-model fruit flies (Diptera: Tephritidae) and evidence of histone modifications

Interactions between DNA and proteins located in the cell nucleus play an important role in controlling physiological processes by specifying, augmenting and regulating context-specific transcription events. Chromatin immunoprecipitation (ChIP) is a widely used methodology to study DNA-protein interactions and has been successfully used in various cell types for over three decades. More recently, by combining ChIP with genomic screening technologies and Next Generation Sequencing (e.g. ChIP-seq), it has become possible to profile DNA-protein interactions (including covalent histone modifications) across entire genomes. However, the applicability of ChIP-chip and ChIP-seq has rarely been extended to non-model species because of a number of technical challenges. Here we report a method that can be used to identify genome wide covalent histone modifications in a group of non-model fruit fly species (Diptera: Tephritidae). The method was developed by testing and refining protocols that have been used in model organisms, including Drosophila melanogaster. We demonstrate that this method is suitable for a group of economically important pest fruit fly species, viz., Bactrocera dorsalis, Ceratitis capitata, Zeugodacus cucurbitae and Bactrocera tryoni. We also report an example ChIP-seq dataset for B. tryoni, providing evidence for histone modifications in the genome of a tephritid fruit fly for the first time. Since tephritids are major agricultural pests globally, this methodology will be a valuable resource to study taxa-specific evolutionary questions and to assist with pest management. It also provides a basis for researchers working with other non-model species to undertake genome wide DNA-protein interaction studies.

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.

Screening mitochondrial DNA sequence variation as an alternative method for tracking established and outbreak populations of Queensland fruit fly at the species southern range limit

Understanding the relationship between incursions of insect pests and established populations is critical to implementing effective control. Studies of genetic variation can provide powerful tools to examine potential invasion pathways and longevity of individual pest outbreaks. The major fruit fly pest in eastern Australia, Queensland fruit fly Bactrocera tryoni (Froggatt), has been subject to significant long‐term quarantine and population reduction control measures in the major horticulture production areas of southeastern Australia, at the species southern range limit. Previous studies have employed microsatellite markers to estimate gene flow between populations across this region. In this study, we used an independent genetic marker, mitochondrial DNA (mtDNA) sequences, to screen genetic variation in established and adjacent outbreak populations in southeastern Australia. During the study period, favorable environmental conditions resulted in multiple outbreaks, which appeared genetically distinctive and relatively geographically localized, implying minimal dispersal between simultaneous outbreaks. Populations in established regions were found to occur over much larger areas. Screening mtDNA (female) lineages proved to be an effective alternative genetic tool to assist in understanding fruit fly population dynamics and provide another possible molecular method that could now be employed for better understanding of the ecology and evolution of this and other pest species.

Helitrons shaping the genomic architecture of Drosophila: enrichment of DINE-TR1 in α- and β-heterochromatin, satellite DNA emergence, and piRNA expression

Drosophila INterspersed Elements (DINEs) constitute an abundant but poorly understood group of Helitrons present in several Drosophila species. The general structure of DINEs includes two conserved blocks that may or not contain a region with tandem repeats in between. These central tandem repeats (CTRs) are similar within species but highly divergent between species. It has been assumed that CTRs have independent origins. Herein, we identify a subset of DINEs, termed DINE-TR1, which contain homologous CTRs of approximately 150 bp. We found DINE-TR1 in the sequenced genomes of several Drosophila species and in Bactrocera tryoni (Acalyptratae, Diptera). However, interspecific high sequence identity (∼ 88 %) is limited to the first ∼ 30 bp of each tandem repeat, implying that evolutionary constraints operate differently over the monomer length. DINE-TR1 is unevenly distributed across the Drosophila phylogeny. Nevertheless, sequence analysis suggests vertical transmission. We found that CTRs within DINE-TR1 have independently expanded into satellite DNA-like arrays at least twice within Drosophila. By analyzing the genome of Drosophila virilis and Drosophila americana, we show that DINE-TR1 is highly abundant in pericentromeric heterochromatin boundaries, some telomeric regions and in the Y chromosome. It is also present in the centromeric region of one autosome from D. virilis and dispersed throughout several euchromatic sites in both species. We further found that DINE-TR1 is abundant at piRNA clusters, and small DINE-TR1-derived RNA transcripts (∼25 nt) are predominantly expressed in the testes and the ovaries, suggesting active targeting by the piRNA machinery. These features suggest potential piRNA-mediated regulatory roles for DINEs at local and genome-wide scales in Drosophila.

Extraordinary conservation of entire chromosomes in insects over long evolutionary periods

Comparison of the genomes of different Drosophila species has shown that six different chromosomes, the so-called ''Muller elements," constitute the building blocks for all Drosophila species. Here, we confirm previous results suggesting that this conservation of the Muller elements extends far beyond Drosophila, to at least tephritid fruit flies, thought to have diverged from drosophilids 60-70 mYr ago. Less than 10 percent of genes differ in chromosome location between the two insect groups. Within chromosomes, however, the order is highly scrambled, as expected from the comparison between Drosophila species. The data also support the notion that the sex chromosomes of tephritid flies originated from an ancestor of the dot chromosome 4 of Drosophila. Overall, therefore, no new chromosome has been created for perhaps a billion generations over the two evolutionary lines. This stability at the chromosome level, which appears to extend to all Diptera including mosquitoes, is in stark contrast to other groups such as mammals, birds, fish and plants, in which chromosome numbers and organization vary enormously among species that have diverged over much fewer generations.

Resolving cryptic species complexes of major tephritid pests

An FAO/IAEA Co-ordinated Research Project (CRP) on "Resolution of Cryptic Species Complexes of Tephritid Pests to Overcome Constraints to SIT Application and International Trade" was conducted from 2010 to 2015. As captured in the CRP title, the objective was to undertake targeted research into the systematics and diagnostics of taxonomically challenging fruit fly groups of economic importance. The scientific output was the accurate alignment of biological species with taxonomic names; which led to the applied outcome of assisting FAO and IAEA Member States in overcoming technical constraints to the application of the Sterile Insect Technique (SIT) against pest fruit flies and the facilitation of international agricultural trade. Close to 50 researchers from over 20 countries participated in the CRP, using coordinated, multidisciplinary research to address, within an integrative taxonomic framework, cryptic species complexes of major tephritid pests. The following progress was made for the four complexes selected and studied: Anastrepha fraterculus complex - Eight morphotypes and their geographic and ecological distributions in Latin America were defined. The morphotypes can be considered as distinct biological species on the basis of differences in karyotype, sexual incompatibility, post-mating isolation, cuticular hydrocarbon, pheromone, and molecular analyses. Discriminative taxonomic tools using linear and geometric morphometrics of both adult and larval morphology were developed for this complex. Bactrocera dorsalis complex - Based on genetic, cytogenetic, pheromonal, morphometric, and behavioural data, which showed no or only minor variation between the Asian/African pest fruit flies Bactrocera dorsalis, Bactrocera papayae, Bactrocera philippinensis and Bactrocera invadens, the latter three species were synonymized with Bactrocera dorsalis. Of the five target pest taxa studied, only Bactrocera dorsalis and Bactrocera carambolae remain as scientifically valid names. Molecular and pheromone markers are now available to distinguish Bactrocera dorsalis from Bactrocera carambolae. Ceratitis FAR Complex (Ceratitis fasciventris, Ceratitis anonae, Ceratitis rosa) - Morphology, morphometry, genetic, genomic, pheromone, cuticular hydrocarbon, ecology, behaviour, and developmental physiology data provide evidence for the existence of five different entities within this fruit fly complex from the African region. These are currently recognised as Ceratitis anonae, Ceratitis fasciventris (F1 and F2), Ceratitis rosa and a new species related to Ceratitis rosa (R2). The biological limits within Ceratitis fasciventris (i.e. F1 and F2) are not fully resolved. Microsatellites markers and morphological identification tools for the adult males of the five different FAR entities were developed based on male leg structures. Zeugodacus cucurbitae (formerly Bactrocera (Zeugodacus) cucurbitae) - Genetic variability was studied among melon fly populations throughout its geographic range in Africa and the Asia/Pacific region and found to be limited. Cross-mating studies indicated no incompatibility or sexual isolation. Host preference and genetic studies showed no evidence for the existence of host races. It was concluded that the melon fly does not represent a cryptic species complex, neither with regard to geographic distribution nor to host range. Nevertheless, the higher taxonomic classification under which this species had been placed, by the time the CRP was started, was found to be paraphyletic; as a result the subgenus Zeugodacus was elevated to genus level.