Analysis of the Olive Fruit Fly Bactrocera oleae Transcriptome and Phylogenetic Classification of the Major Detoxification Gene Families

Insights on reproduction-related genes in the striped fruit fly, Zeugodacus scutellata (Hendel) (Diptera: Tephritidae)

The striped fruit fly, Zeugodacus scutellata is a significant pest in East and Southeast Asia by damaging Cucurbitaceae blossoms and fruits. To control this pest, a novel strategy to suppress the gene(s) associated with sexually dimorphic phenotypes has been devised and implemented in a laboratory scale. However, comprehensive transcriptomic analysis related to this sex differentiation of Z. scutellata was necessary to determine effective target genes for the genetic control. We performed de novo assembly of the transcript obtained by paired-end sequencing using an Illumina HiSeq platform and let to 217,967 unigenes (i.e., unique genes) with a minimum length of 200 bp. The female produced 31, 604, 442 reads with 97.93% of Q20, 94.76% of Q30, and the male produced 130, 592, 828 reads with 97.93% of Q20 and 94.76 of Q30%. The differentially expressed genes were used to predict genetic factors associated with sex differentiation, which included Rho1, extra-macrochaetae (emc), hopscotch (hop), doublesex (dsx), sex-lethal (sxl), transformer-2 (tra-2), testis-specific serine/threonine-protein kinase (tssk1), tektin1 (tkt1) and 2 (tkt2), odorant binding proteins (OBPs), fruitless (fru), vitellogenin receptor, and hormone receptors in Z. scutellata. In addition, this transcriptome analysis provides the additional gene associated with sex determination and mating behaviors, which would be applied to develop a novel sterile insect technique against Z. scutellata.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

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.

Nanopore long-read RNA-seq and absolute quantification delineate transcription dynamics in early embryo development of an insect pest

The olive fruit fly, Bactrocera oleae, is the most important pest for the olive fruit but lacks adequate transcriptomic characterization that could aid in molecular control approaches. We apply nanopore long-read RNA-seq with internal RNA standards allowing absolute transcript quantification to analyze transcription dynamics during early embryo development for the first time in this organism. Sequencing on the MinION platform generated over 31 million reads. Over 50% of the expressed genes had at least one read covering its entire length validating our full-length approach. We generated a de novo transcriptome assembly and identified 1768 new genes and a total of 79,810 isoforms; a fourfold increase in transcriptome diversity compared to the current NCBI predicted transcriptome. Absolute transcript quantification per embryo allowed an insight into the dramatic re-organization of maternal transcripts. We further identified Zelda as a possible regulator of early zygotic genome activation in B. oleae and provide further insights into the maternal-to-zygotic transition. These data show the utility of long-read RNA in improving characterization of non-model organisms that lack a fully annotated genome, provide potential targets for sterile insect technic approaches, and provide the first insight into the transcriptome landscape of the developing olive fruit fly embryo.

Comparative Microbiomics of Tephritid Frugivorous Pests (Diptera: Tephritidae) From the Field: A Tale of High Variability Across and Within Species

The family Tephritidae includes some of the most notorious insect pests of agricultural and horticultural crops in tropical and sub-tropical regions. Despite the interest in the study of their gut microbiome, our present knowledge is largely based on the analysis of laboratory strains. In this study, we present a first comparative analysis of the gut microbiome profiles of field populations of ten African and Mediterranean tephritid pests. For each species, third instar larvae were sampled from different locations and host fruits and compared using 16S rRNA amplicon sequencing and a multi-factorial sampling design. We observed considerable variation in gut microbiome diversity and composition both between and within fruit fly species. A “core” microbiome, shared across all targeted species, could only be identified at most at family level (Enterobacteriaceae). At genus level only a few bacterial genera (Klebsiella, Enterobacter, and Bacillus) were present in most, but not all, samples, with high variability in their relative abundance. Higher relative abundances were found for seven bacterial genera in five of the fruit fly species considered. These were Erwinia in Bactrocera oleae, Lactococcus in B. zonata, Providencia in Ceratitis flexuosa, Klebsiella, and Rahnella in C. podocarpi and Acetobacter and Serratia in C. rosa. With the possible exception of C. capitata and B. dorsalis (the two most polyphagous species considered) we could not detect obvious relationships between fruit fly dietary breadth and microbiome diversity or abundance patterns. Similarly, our results did not suggest straightforward differences between the microbiome profiles of species belonging to Ceratitis and the closely related Bactrocera/Zeugodacus. These results provide a first comparative analysis of the gut microbiomes of field populations of multiple economically relevant tephritids and provide base line information for future studies that will further investigate the possible functional role of the observed associations.

Knocking down the sex peptide receptor by dsRNA feeding results in reduced oviposition rate in olive fruit flies

Insect pests can cause crop damage in yield or quality, resulting in profit losses for farmers. The primary approach to control them is still the use of chemical pesticides resulting in significant hazards to the environment and human health. Biological control and the sterile insect technique are alternative strategies to improve agriculture protection. However, both strategies have significant limitations. A newly introduced approach that could be both effective and species-specific is the RNA interference mechanism. One key point for the success of this strategy is the delivery method of double-strand RNA (dsRNA) to the insects. A method of dsRNA delivery to insects with potential use in the field is the oral delivery, feeding the insects engineered microorganisms that produce dsRNA. Here, we present the first protocol for dsRNA feeding using modified bacteria, in the olive fruit fly, the most important insect pest of cultivated olives. We chose to target the sex peptide receptor gene. The sex peptide receptor interacts with the sex peptide, a peptide that is responsible for the postmating behavior in the model organism, Drosophila melanogaster. Feeding the female olive fruit fly with bacteria that produced dsRNA for the sex peptide receptor gene resulted in the development of female insects with significantly lower oviposition rates. Administration of dsRNA producing bacteria in insect diet against target genes that lead to genetic sexing or female-specific lethality could be added in the armory of control methods.

De novo assembly of the olive fruit fly (Bactrocera oleae) genome with linked-reads and long-read technologies minimizes gaps and provides exceptional Y chromosome assembly

BACKGROUND

The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly's biology and proposing alternative control methods to pesticide use.

RESULTS

Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR.

CONCLUSIONS

The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome's organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Characterisation of GST genes from the Hyphantria cunea and their response to the oxidative stress caused by the infection of Hyphantria cunea nucleopolyhedrovirus (HcNPV)

The fall webworm, Hyphantria cunea (Drury) (Lepidoptera: Noctuidae), is a major pest found in forests. In this study, the effects of Hyphantria cunea nucleopolyhedrovirus (HcNPV) infection on the transcription levels and activities of glutathione S-transferases (GSTs) in H. cunea were determined. In the present study, 18 GST family genes were identified from the H. cunea transcriptome dataset by using bioinformatic analyses. These GST genes were classified into cytosolic (15 genes) and microsomal (three genes) classes. The 15 cytosolic GST genes belonged to four different subclasses (epsilon, sigma and delta). The all GST genes, especially GSTe4, showed high expression levels in egg and 1st~4th instar larval stage while their low expression levels in 5th~7th instar larvae using real-time quantitative PCR analysis. However, the expression levels of the 18 GST genes were varied after exposure to sublethal doses of HcNPV. The expression levels of most GSTs were downregulated and upregulated at low and high concentrations of HcNPV, respectively. The corresponding total GST activities also showed similar patterns. In H. cunea, changes in the expression levels and enzymatic activities of GSTs after exposure to HcNPV indicated that they may have important functions in the defense against HcNPV, and the stress, which may be reflected by the high GST enzymatic activities.

Symbiosis in Sustainable Agriculture: Can Olive Fruit Fly Bacterial Microbiome Be Useful in Pest Management?

The applied importance of symbiosis has been gaining recognition. The relevance of symbiosis has been increasing in agriculture, in developing sustainable practices, including pest management. Insect symbiotic microorganisms' taxonomical and functional diversity is high, and so is the potential of manipulation of these microbial partners in suppressing pest populations. These strategies, which rely on functional organisms inhabiting the insect, are intrinsically less susceptible to external environmental variations and hence likely to overcome some of the challenges posed by climate change. Rates of climate change in the Mediterranean Basin are expected to exceed global trends for most variables, and this warming will also affect olive production and impact the interactions of olives and their main pest, the obligate olive fruit fly (Bactrocera oleae). This work summarizes the current knowledge on olive fly symbiotic bacteria towards the potential development of symbiosis-based strategies for olive fruit fly control. Particular emphasis is given to Candidatus Erwinia dacicola, an obligate, vertically transmitted endosymbiont that allows the insect to cope with the olive-plant produced defensive compound oleuropein, as a most promising target for a symbiosis disruption approach.

Impact of Bactrocera oleae on the fungal microbiota of ripe olive drupes

The olive fruit fly (OFF), Bactrocera oleae is the most devastating pest affecting olive fruit worldwide. Previous investigations have addressed the fungal microbiome associated with olive drupes or B. oleae, but the impact of the insect on fungal communities of olive fruit remains undescribed. In the present work, the fungal microbiome of olive drupes, infested and non-infested by the OFF, was investigated in four different localities and cultivars. Olive fruit fly infestations caused a general reduction of the fungal diversity, a higher quantity of the total DNA and an increase in taxa that remained unidentified or had unknown roles. The infestations led to imbalanced fungal communities with the growth of taxa that are usually outcompeted. While it was difficult to establish a cause-effect link between fly infestation and specific fungi, it is clear that the fly alters the natural microbial balance, especially the low abundant taxa. On the other hand, the most abundant ones, were not significantly influenced by the insect. In fact, despite the slight variation between the sampling locations, Aureobasidium, Cladosporium, and Alternaria, were the dominant genera, suggesting the existence of a typical olive fungal microbiome.

Residual degradation and toxicity of insecticides against Bactrocera oleae

Field and lab trials took place in Crete (July to September 2016), concerning the residual degradation and toxicity of seven active ingredients applied as bait sprays against the olive fruit fly. Highest residues were recorded in olive leaves for dimethoate and phosmet (~ 60 mg/kg) immediately after application (day 1+), while a threefold and fivefold reduction was observed 1 week later, respectively. Residues of pyrethroids were determined at lower levels (< 10 mg/kg) but remained almost stable for a longer period of time. Finally, thiacloprid and spinosad residues were determined at 5.81 and 0.19 mg/kg respectively (day 1+), and rapidly decreased below the LOQ. Highest toxicity against the olive fruit fly was observed just right after the application of dimethoate (100%), a-cypermethrin (80%), and L-cyhalothrin (72.92%). Although the toxicity of dimethoate was significantly reduced 1 week after the application (80%) and then minimized, toxicity of pyrethroids remained almost stable (> 60%) for the first 2 weeks and then decreased to 30-40%, which remained stable up to the end of the study (8 weeks). Concerning phosmet, its toxicity ranged from 35 to 56% for 3 weeks with no significant reduction, while spinosad presented a lower toxicity profile (50% only for 1 week). The benefits of these results in the knowledge of insecticide residues and their toxicity against olive fruit fly can be used for improving olive fruit fly control.

Transcriptional responses of detoxification genes to four plant allelochemicals in Aphis gossypii

Aphis gossypii Glover (Hemiptera: Aphididae) can damage a variety of agricultural crops, so it is very important for cotton aphids to evolve adaptive mechanisms to various allelochemicals from host plants. Our results aim to provide a fundamental and rich resource for exploring aphid functional genes in A. gossypii. A transcriptome data set and five expression profile data sets of A. gossypii samples were analyzed by Illumina sequencing platform. In total, 53,763,866 reads were assembled into 1,963,516 contigs and 28,555 unigenes. Compared with the control, 619 genes were significantly up- or downregulated in the treatment group by 2-tridecanone. There were 516, 509, and 717 of differential expression genes in tannic acid, quercetin, and gossypol treatment groups, respectively. Furthermore, there were 4 of 54 putative cytochrome P450 genes and 1 of 7 putative carboxylesterases downregulated in all treatment groups by four plant allelochemicals. When aphids fed on 2-tridecanone, tannic acid, and quercetin, only one P450 gene was upregulated. These results show that plant allelochemical stress can induce differential gene expression in A. gossypii. The differential response information of gene expression based on a large-scale sequence would be useful to reveal molecular mechanisms of adaptation for A. gossypii to plant allelochemicals.

Transcriptomic responses of the olive fruit fly Bactrocera oleae and its symbiont Candidatus Erwinia dacicola to olive feeding

The olive fruit fly, Bactrocera oleae, is the most destructive pest of olive orchards worldwide. The monophagous larva has the unique capability of feeding on olive mesocarp, coping with high levels of phenolic compounds and utilizing non-hydrolyzed proteins present, particularly in the unripe, green olives. On the molecular level, the interaction between B. oleae and olives has not been investigated as yet. Nevertheless, it has been associated with the gut obligate symbiotic bacterium Candidatus Erwinia dacicola. Here, we used a B.oleae microarray to analyze the gene expression of larvae during their development in artificial diet, unripe (green) and ripe (black) olives. The expression profiles of Ca. E. dacicola were analyzed in parallel, using the Illumina platform. Several genes were found overexpressed in the olive fly larvae when feeding in green olives. Among these, a number of genes encoding detoxification and digestive enzymes, indicating a potential association with the ability of B. oleae to cope with green olives. In addition, a number of biological processes seem to be activated in Ca. E. dacicola during the development of larvae in olives, with the most notable being the activation of amino-acid metabolism.

Transcriptome Analysis and Identification of Major Detoxification Gene Families and Insecticide Targets in Grapholita Molesta (Busck) (Lepidoptera: Tortricidae)

The oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae), is an important pest of most stone and pome fruits and causes serious damage to the fruit industry worldwide. This insect pest has been primarily controlled through the application of insecticides; as a result, G. molesta has developed resistance to many different types of insecticides. To identify detoxification genes, we have, de novo, sequenced the transcriptome of G. molesta (Lepidoptera: Tortricidae) and yielded 58,970 unigenes of which 26,985 unigenes matched to known proteins. In total, 2,040 simple sequence repeats have been identified. The comprehensive transcriptome data set has permitted us to identify members of important gene families related to detoxification in G. molesta, including 77 unigenes of putative cytochrome P450s, 28 of glutathione S-transferases, 46 of Carboxylesterases, and 31 of insecticide targets. Orthologs of some of these unigenes have shown to play a pivotal role in insecticide resistance in other insect species and those unigenes likely have similar functions in G. molesta.

Transcriptome Analysis of an Insecticide Resistant Housefly Strain: Insights about SNPs and Regulatory Elements in Cytochrome P450 Genes

Background

Insecticide resistance in the housefly, Musca domestica, has been investigated for more than 60 years. It will enter a new era after the recent publication of the housefly genome and the development of multiple next generation sequencing technologies. The genetic background of the xenobiotic response can now be investigated in greater detail. Here, we investigate the 454-pyrosequencing transcriptome of the spinosad-resistant 791spin strain in relation to the housefly genome with focus on P450 genes.

Results

The de novo assembly of clean reads gave 35,834 contigs consisting of 21,780 sequences of the spinosad resistant strain. The 3,648 sequences were annotated with an enzyme code EC number and were mapped to 124 KEGG pathways with metabolic processes as most highly represented pathway. One hundred and twenty contigs were annotated as P450s covering 44 different P450 genes of housefly. Eight differentially expressed P450s genes were identified and investigated for SNPs, CpG islands and common regulatory motifs in promoter and coding regions. Functional annotation clustering of metabolic related genes and motif analysis of P450s revealed their association with epigenetic, transcription and gene expression related functions. The sequence variation analysis resulted in 12 SNPs and eight of them found in cyp6d1. There is variation in location, size and frequency of CpG islands and specific motifs were also identified in these P450s. Moreover, identified motifs were associated to GO terms and transcription factors using bioinformatic tools.

Conclusion

Transcriptome data of a spinosad resistant strain provide together with genome data fundamental support for future research to understand evolution of resistance in houseflies. Here, we report for the first time the SNPs, CpG islands and common regulatory motifs in differentially expressed P450s. Taken together our findings will serve as a stepping stone to advance understanding of the mechanism and role of P450s in xenobiotic detoxification.

Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype

Background

Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids.

Results

We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated.

A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains.

DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases.

Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively.

Conclusion

Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2161-5) contains supplementary material, which is available to authorized users.

Sequence Analysis of Insecticide Action and Detoxification-Related Genes in the Insect Pest Natural Enemy Pardosa pseudoannulata

The pond wolf spider Pardosa pseudoannulata, an important natural predatory enemy of rice planthoppers, is found widely distributed in paddy fields. However, data on the genes involved in insecticide action, detoxification, and response are very limited for P. pseudoannulata, which inhibits the development and appropriate use of selective insecticides to control insect pests on rice. We used transcriptome construction from adult spider cephalothoraxes to analyze and manually identify genes enconding metabolic enzymes and target receptors related to insecticide action and detoxification, including 90 cytochrome P450s, 14 glutathione S-transferases (GSTs), 17 acetylcholinesterases (AChEs), 17 nicotinic acetylcholine receptors (nAChRs), and 17 gamma-aminobutyric acid (GABA) receptors, as well as 12 glutamate-gated chloride channel (GluCl) unigenes. Sequence alignment and phylogenetic analysis revealed the different subclassifications of P450s and GSTs, some important sequence diversities in nAChRs and GABA receptors, polymorphism in AChEs, and high similarities in GluCls. For P450s in P. pseudoannulata, the number of unigenes belonging to the CYP2 clade was much higher than that in CYP3 and CYP4 clades. The results differed from insects in which most P450 genes were in CYP3 and CYP4 clades. For GSTs, most unigenes belonged to the delta and sigma classes, and no epsilon GST class gene was found, which differed from the findings for insects and acarina. Our results will be useful for studies on insecticide action, selectivity, and detoxification in the spider and other related animals, and the sequence differences in target genes between the spider and insects will provide important information for the design of selective insecticides.

The molecular biology of the olive fly comes of age

Background

Olive cultivation blends with the history of the Mediterranean countries since ancient times. Even today, activities around the olive tree constitute major engagements of several people in the countryside of both sides of the Mediterranean basin. The olive fly is, beyond doubt, the most destructive pest of cultivated olives. The female fly leaves its eggs in the olive fruit. Upon emergence, the larvae feed on the olive sap, thus destroying the fruit. If untreated, practically all olives get infected. The use of chemical insecticides constitutes the principal olive fly control approach. The Sterile Insect Technique (SIT), an environmentally friendly alternative control method, had been tried in pilot field applications in the 1970's, albeit with no practical success. This was mainly attributed to the low, non-antagonistic quality of the mixed-sex released insects. Many years of experience from successful SIT applications in related species, primarily the Mediterranean fruit fly, Ceratitis capitata, demonstrated that efficient SIT protocols require the availability of fundamental genetic and molecular information.

Results

Among the primary systems whose understanding can contribute towards novel SIT approaches (or its recently developed alternative RIDL: Release of Insects carrying a Dominant Lethal) is the reproductive, since the ability to manipulate the reproductive system would directly affect the insect's fertility. In addition, the analysis of early embryonic promoters and apoptotic genes would provide tools that confer dominant early-embryonic lethality during mass-rearing. Here we report the identification of several genes involved in these systems through whole transcriptome analysis of female accessory glands (FAGs) and spermathecae, as well as male testes. Indeed, analysis of differentially expressed genes in these tissues revealed higher metabolic activity in testes than in FAGs/spermathecae. Furthermore, at least five olfactory-related genes were shown to be differentially expressed in the female and male reproductive systems analyzed. Finally, the expression profile of the embryonic serendipity-α locus and the pre-apoptotic head involution defective gene were analyzed during embryonic developmental stages.

Conclusions

Several years of molecular studies on the olive fly can now be combined with new information from whole transcriptome analyses and lead to a deep understanding of the biology of this notorious insect pest. This is a prerequisite for the development of novel embryonic lethality female sexing strains for successful SIT efforts which, combined with improved mass-reared conditions, give new hope for efficient SIT applications for the olive fly.

How functional genomics will impact fruit fly pest control: the example of the Mediterranean fruit fly, Ceratitis capitata

The highly invasive agricultural insect pest Ceratitis capitata (Diptera: Tephritidae) is the most thoroughly studied tephritid fruit fly at the genetic and molecular levels. It has become a model for the analysis of fruit fly invasions and for the development of area-wide integrated pest management (AW-IPM) programmes based on the environmentally-friendly Sterile Insect Technique (SIT). Extensive transcriptome resources and the recently released genome sequence are making it possible to unravel several aspects of the medfly reproductive biology and behaviour, opening new opportunities for comparative genomics and barcoding for species identification. New genes, promotors and regulatory sequences are becoming available for the development/improvement of highly competitive sexing strains, for the monitoring of sterile males released in the field and for determining the mating status of wild females. The tools developed in this species have been transferred to other tephritids that are also the subject of SIT programmes.

Australian endemic pest tephritids: genetic, molecular and microbial tools for improved Sterile Insect Technique

Among Australian endemic tephritid fruit flies, the sibling species Bactrocera tryoni and Bactrocera neohumeralis have been serious horticultural pests since the introduction of horticulture in the nineteenth century. More recently, Bactrocera jarvisi has also been declared a pest in northern Australia. After several decades of genetic research there is now a range of classical and molecular genetic tools that can be used to develop improved Sterile Insect Technique (SIT) strains for control of these pests. Four-way crossing strategies have the potential to overcome the problem of inbreeding in mass-reared strains of B. tryoni. The ability to produce hybrids between B. tryoni and the other two species in the laboratory has proved useful for the development of genetically marked strains. The identification of Y-chromosome markers in B. jarvisi means that male and female embryos can be distinguished in any strain that carries a B. jarvisi Y chromosome. This has enabled the study of homologues of the sex-determination genes during development of B jarvisi and B. tryoni, which is necessary for the generation of genetic-sexing strains. Germ-line transformation has been established and a draft genome sequence for B. tryoni released. Transcriptomes from various species, tissues and developmental stages, to aid in identification of manipulation targets for improving SIT, have been assembled and are in the pipeline. Broad analyses of the microbiome have revealed a metagenome that is highly variable within and across species and defined by the environment. More specific analyses detected Wolbachia at low prevalence in the tropics but absent in temperate regions, suggesting a possible role for this endosymbiont in future control strategies.

Transcriptome-based identification of ABC transporters in the western tarnished plant bug Lygus hesperus

ATP-binding cassette (ABC) transporters are a large superfamily of proteins that mediate diverse physiological functions by coupling ATP hydrolysis with substrate transport across lipid membranes. In insects, these proteins play roles in metabolism, development, eye pigmentation, and xenobiotic clearance. While ABC transporters have been extensively studied in vertebrates, less is known concerning this superfamily in insects, particularly hemipteran pests. We used RNA-Seq transcriptome sequencing to identify 65 putative ABC transporter sequences (including 36 full-length sequences) from the eight ABC subfamilies in the western tarnished plant bug (Lygus hesperus), a polyphagous agricultural pest. Phylogenetic analyses revealed clear orthologous relationships with ABC transporters linked to insecticide/xenobiotic clearance and indicated lineage specific expansion of the L. hesperus ABCG and ABCH subfamilies. The transcriptional profile of 13 LhABCs representative of the ABCA, ABCB, ABCC, ABCG, and ABCH subfamilies was examined across L. hesperus development and within sex-specific adult tissues. All of the transcripts were amplified from both reproductively immature and mature adults and all but LhABCA8 were expressed to some degree in eggs. Expression of LhABCA8 was spatially localized to the testis and temporally timed with male reproductive development, suggesting a potential role in sexual maturation and/or spermatozoa protection. Elevated expression of LhABCC5 in Malpighian tubules suggests a possible role in xenobiotic clearance. Our results provide the first transcriptome-wide analysis of ABC transporters in an agriculturally important hemipteran pest and, because ABC transporters are known to be important mediators of insecticidal resistance, will provide the basis for future biochemical and toxicological studies on the role of this protein family in insecticide resistance in Lygus species.

Characterizing the developmental transcriptome of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae) through comparative genomic analysis with Drosophila melanogaster utilizing modENCODE datasets

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis, is an important pest of fruit and vegetable crops throughout Asia, and is considered a high risk pest for establishment in the mainland United States. It is a member of the family Tephritidae, which are the most agriculturally important family of flies, and can be considered an out-group to well-studied members of the family Drosophilidae. Despite their importance as pests and their relatedness to Drosophila, little information is present on B. dorsalis transcripts and proteins. The objective of this paper is to comprehensively characterize the transcripts present throughout the life history of B. dorsalis and functionally annotate and analyse these transcripts relative to the presence, expression, and function of orthologous sequences present in Drosophila melanogaster.

RESULTS

We present a detailed transcriptome assembly of B. dorsalis from egg through adult stages containing 20,666 transcripts across 10,799 unigene components. Utilizing data available through Flybase and the modENCODE project, we compared expression patterns of these transcripts to putative orthologs in D. melanogaster in terms of timing, abundance, and function. In addition, temporal expression patterns in B. dorsalis were characterized between stages, to establish the constitutive or stage-specific expression patterns of particular transcripts. A fully annotated transcriptome assembly is made available through NCBI, in addition to corresponding expression data.

CONCLUSIONS

Through characterizing the transcriptome of B. dorsalis through its life history and comparing the transcriptome of B. dorsalis to the model organism D. melanogaster, a database has been developed that can be used as the foundation to functional genomic research in Bactrocera flies and help identify orthologous genes between B. dorsalis and D. melanogaster. This data provides the foundation for future functional genomic research that will focus on improving our understanding of the physiology and biology of this species at the molecular level. This knowledge can also be applied towards developing improved methods for control, survey, and eradication of this important pest.

Olive fly transcriptomics analysis implicates energy metabolism genes in spinosad resistance

Background

The olive fly, Bactrocera oleae, is the most devastating pest of cultivated olives. Its control has been traditionally based on insecticides, mainly organophosphates and pyrethroids. In recent years, the naturalyte spinosad is used against the olive fly. As with other insecticides, spinosad is subject to selection pressures that have led to resistance development. Mutations in the α6 subunit of the nicotinic acetylcholine receptor (nAChR) have been implicated in spinosad resistance in several species (e.g., Drosophila melanogaster) but excluded in others (e.g., Musca domestica). Yet, additional mechanisms involving enhanced metabolism of detoxification enzymes (such as P450 monooxygenases or mixed function oxidases) have also been reported. In order to clarify the spinosad resistance mechanisms in the olive fly, we searched for mutations in the α6-subunit of the nAChR and for up-regulated genes in the entire transcriptome of spinosad resistant olive flies.

Results

The olive fly α6-subunit of the nAChR was cloned from the laboratory sensitive strain and a spinosad selected resistant line. The differences reflected silent nucleotide substitutions or conserved amino acid changes. Additionally, whole transcriptome analysis was performed in the two strains in order to reveal any underlying resistance mechanisms. Comparison of over 13,000 genes showed that in spinosad resistant flies nine genes were significantly over-expressed, whereas ~40 were under-expressed. Further functional analyses of the nine over-expressed and eleven under-expressed loci were performed. Four of these loci (Yolk protein 2, ATP Synthase FO subunit 6, Low affinity cationic amino acid transporter 2 and Serine protease 6) showed consistently higher expression both in the spinosad resistant strain and in wild flies from a resistant California population. On the other side, two storage protein genes (HexL1 and Lsp1) and two heat-shock protein genes (Hsp70 and Hsp23) were unfailingly under-expressed in resistant flies.

Conclusion

The observed nucleotide differences in the nAChR-α6 subunit between the sensitive and spinosad resistant olive fly strains did not advocate for the involvement of receptor mutations in spinosad resistance. Instead, the transcriptome comparison between the two strains indicated that several immune system loci as well as elevated energy requirements of the resistant flies might be necessary to lever the detoxification process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-714) contains supplementary material, which is available to authorized users.

De novo cloning and annotation of genes associated with immunity, detoxification and energy metabolism from the fat body of the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis, is a destructive pest in tropical and subtropical areas. In this study, we performed transcriptome-wide analysis of the fat body of B. dorsalis and obtained more than 59 million sequencing reads, which were assembled into 27,787 unigenes with an average length of 591 bp. Among them, 17,442 (62.8%) unigenes matched known proteins in the NCBI database. The assembled sequences were further annotated with gene ontology, cluster of orthologous group terms, and Kyoto encyclopedia of genes and genomes. In depth analysis was performed to identify genes putatively involved in immunity, detoxification, and energy metabolism. Many new genes were identified including serpins, peptidoglycan recognition proteins and defensins, which were potentially linked to immune defense. Many detoxification genes were identified, including cytochrome P450s, glutathione S-transferases and ATP-binding cassette (ABC) transporters. Many new transcripts possibly involved in energy metabolism, including fatty acid desaturases, lipases, alpha amylases, and trehalose-6-phosphate synthases, were identified. Moreover, we randomly selected some genes to examine their expression patterns in different tissues by quantitative real-time PCR, which indicated that some genes exhibited fat body-specific expression in B. dorsalis. The identification of a numerous transcripts in the fat body of B. dorsalis laid the foundation for future studies on the functions of these genes.

De novo assembly, gene annotation, and marker discovery in stored-product pest Liposcelis entomophila (Enderlein) using transcriptome sequences

Background

As a major stored-product pest insect, Liposcelis entomophila has developed high levels of resistance to various insecticides in grain storage systems. However, the molecular mechanisms underlying resistance and environmental stress have not been characterized. To date, there is a lack of genomic information for this species. Therefore, studies aimed at profiling the L. entomophila transcriptome would provide a better understanding of the biological functions at the molecular levels.

Methodology/Principal Findings

We applied Illumina sequencing technology to sequence the transcriptome of L. entomophila. A total of 54,406,328 clean reads were obtained and that de novo assembled into 54,220 unigenes, with an average length of 571 bp. Through a similarity search, 33,404 (61.61%) unigenes were matched to known proteins in the NCBI non-redundant (Nr) protein database. These unigenes were further functionally annotated with gene ontology (GO), cluster of orthologous groups of proteins (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of genes potentially involved in insecticide resistance were manually curated, including 68 putative cytochrome P450 genes, 37 putative glutathione S-transferase (GST) genes, 19 putative carboxyl/cholinesterase (CCE) genes, and other 126 transcripts to contain target site sequences or encoding detoxification genes representing eight types of resistance enzymes. Furthermore, to gain insight into the molecular basis of the L. entomophila toward thermal stresses, 25 heat shock protein (Hsp) genes were identified. In addition, 1,100 SSRs and 57,757 SNPs were detected and 231 pairs of SSR primes were designed for investigating the genetic diversity in future.

Conclusions/Significance

We developed a comprehensive transcriptomic database for L. entomophila. These sequences and putative molecular markers would further promote our understanding of the molecular mechanisms underlying insecticide resistance or environmental stress, and will facilitate studies on population genetics for psocids, as well as providing useful information for functional genomic research in the future.