Long non-coding RNAs regulate Aedes aegypti vector competence for Zika virus and reproduction

Long non-coding RNAs (lncRNAs) play critical regulatory roles in various cellular and metabolic processes in mosquitoes and all other organisms studied thus far. In particular, their involvement in essential processes such as reproduction makes them potential targets for the development of novel pest control approaches. However, their function in mosquito biology remains largely unexplored. To elucidate the role of lncRNAs in mosquitoes' reproduction and vector competence for arboviruses, we have implemented a computational and experimental pipeline to mine, screen, and characterize lncRNAs related to these two biological processes. Through analysis of publicly available Zika virus (ZIKV) infection-regulated Aedes aegypti transcriptomes, at least six lncRNAs were identified as being significantly upregulated in response to infection in various mosquito tissues. The roles of these ZIKV-regulated lncRNAs (designated Zinc1, Zinc2, Zinc3, Zinc9, Zinc10 and Zinc22), were further investigated by dsRNA-mediated silencing studies. Our results show that silencing of Zinc1, Zinc2, and Zinc22 renders mosquitoes significantly less permissive to ZIKV infection, while silencing of Zinc22 also reduces fecundity, indicating a potential role for Zinc22 in trade-offs between vector competence and reproduction. We also found that silencing of Zinc9 significantly increases fecundity but has no effect on ZIKV infection, suggesting that Zinc9 may be a negative regulator of oviposition. Our work demonstrates that some lncRNAs play host factor roles by facilitating viral infection in mosquitoes. We also show that lncRNAs can influence both mosquito reproduction and permissiveness to virus infection, two biological systems with important roles in mosquito vectorial capacity.

Development of Toehold Switches as a Novel Ribodiagnostic Method for West Nile Virus

West Nile virus (WNV) is an emerging neurotropic RNA virus and a member of the genus Flavivirus. Naturally, the virus is maintained in an enzootic cycle involving mosquitoes as vectors and birds that are the principal amplifying virus hosts. In humans, the incubation period for WNV disease ranges from 3 to 14 days, with an estimated 80% of infected persons being asymptomatic, around 19% developing a mild febrile infection and less than 1% developing neuroinvasive disease. Laboratory diagnosis of WNV infection is generally accomplished by cross-reacting serological methods or highly sensitive yet expensive molecular approaches. Therefore, current diagnostic tools hinder widespread surveillance of WNV in birds and mosquitoes that serve as viral reservoirs for infecting secondary hosts, such as humans and equines. We have developed a synthetic biology-based method for sensitive and low-cost detection of WNV. This method relies on toehold riboswitches designed to detect WNV genomic RNA as transcriptional input and process it to GFP fluorescence as translational output. Our methodology offers a non-invasive tool with reduced operating cost and high diagnostic value that can be used for field surveillance of WNV in humans as well as in bird and mosquito populations.

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.

Decoding the Reproductive System of the Olive Fruit Fly, Bactrocera oleae

In most diploid organisms, mating is a prerequisite for reproduction and, thus, critical to the maintenance of their population and the perpetuation of the species. Besides the importance of understanding the fundamentals of reproduction, targeting the reproductive success of a pest insect is also a promising method for its control, as a possible manipulation of the reproductive system could affect its destructive activity. Here, we used an integrated approach for the elucidation of the reproductive system and mating procedures of the olive fruit fly, Bactrocera oleae. Initially, we performed a RNAseq analysis in reproductive tissues of virgin and mated insects. A comparison of the transcriptomes resulted in the identification of genes that are differentially expressed after mating. Functional annotation of the genes showed an alteration in the metabolic, catalytic, and cellular processes after mating. Moreover, a functional analysis through RNAi silencing of two differentially expressed genes, yellow-g and troponin C, resulted in a significantly reduced oviposition rate. This study provided a foundation for future investigations into the olive fruit fly's reproductive biology to the development of new exploitable tools for its control.

Anosmic flies: what Orco silencing does to olive fruit flies

Background

The olive fruit fly (Bactrocera oleae) is the most destructive pest of the olive cultivation worldwide causing significant production losses and olive fruit impoverishment, as its larvae feed exclusively on the olive fruit. Reproductive and sexual behavior, as well as host-plant recognition of the fly, are highly dependent on its chemosensory system. Therefore, exploring the role of genes that play a critical role in olfaction, could reveal potential molecular targets that determine species-specific features on chemical communication and could be used to impair sexual behavior.

Results

In this study we identified the gene that encodes the conserved olfactory co-receptor Orco (Odorant receptor co-receptor), which interacts with all divergent insect odorant receptors, and investigated how disruption of its expression affects chemoreception. We initially searched the expression profile of Bo-Orco in both sexes during sexual maturation, as well as pre- and post-mating communication by relative quantitative real time PCR (qRT-PCR) analysis suggesting that Bo-Orco was abundantly expressed in sexually mature adults.

We further investigated the functional role of Bo-Orco in mating and oviposition behavior via transient gene silencing that was performed through in vivo dsRNA hemolymph injections in sexually mature flies 7 days after eclosion. Orco-knockdown phenotypes in both sexes showed reduced copulation rates in mating competitiveness tests, possibly through impaired olfactory-mediated detection of sex pheromone. In addition, oviposition was significantly inhibited in dsRNA-Orco injected females in a post-mating behavior test.

Conclusions

Our results demonstrate that Orco plays a crucial role in the reproductive behavior of the olive fruit fly, since pre- and post-mating processes were affected. This is the first report in the olive fruit fly that links the chemosensory pathway with the mating behavior and the reproductive potential at a molecular basis, rendering this gene a potential target for the improvement of the olive fruit fly population control techniques.

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.

Targeted somatic mutagenesis through CRISPR/Cas9 ribonucleoprotein complexes in the olive fruit fly, Bactrocera oleae

The olive fruit fly, Bactrocera oleae (Diptera: Tephritidae), is the most destructive insect pest of olive cultivation, causing significant economic and production losses. Here, we present the establishment of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 methodology for gene disruption in this species. We performed targeted mutagenesis of the autosomal gene white (Bo-we), by injecting into early embryos in vitro preassembled and solubilized Cas9 ribonucleoprotein complexes loaded with two gene-specific single-guide RNAs. Gene disruption of Bo-we led to somatic mosaicism of the adult eye color. Large eye patches or even an entire eye lost the iridescent reddish color, indicating the successful biallelic mutagenesis in somatic cells. Cas9 induced either indels in each of the two simultaneously targeted Bo-we sites or a large deletion of the intervening region. This study demonstrates the first efficient implementation of the CRISPR/Cas9 technology in the olive fly, providing new opportunities towards the development of novel genetic tools for its control.

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.

Developing CRISPR-based sex-ratio distorters for the genetic control of fruit fly pests: A how to manual

Agricultural pest control using genetic-based methods provides a species-specific and environmentally harmless way for population suppression of fruit flies. One way to improve the efficiency of such methods is through self-limiting, female-eliminating approaches that can alter an insect populations' sex ratio toward males. In this microreview, we summarize recent advances in synthetic sex ratio distorters based on X-chromosome shredding that can induce male-biased progeny. We outline the basic principles to guide the efficient design of an X-shredding system in an XY heterogametic fruit fly species of interest using CRISPR/Cas gene editing, newly developed computational tools, and insect genetic engineering. We also discuss technical aspects and challenges associated with the efficient transferability of this technology in fruit fly pest populations, toward the potential use of this new class of genetic control approaches for pest management purposes.

Maleness-on-the-Y (MoY) orchestrates male sex determination in major agricultural fruit fly pests

In insects, rapidly evolving primary sex-determining signals are transduced by a conserved regulatory module controlling sexual differentiation. In the agricultural pest Ceratitis capitata (Mediterranean fruit fly, or Medfly), we identified a Y-linked gene, Maleness-on-the-Y (MoY), encoding a small protein that is necessary and sufficient for male development. Silencing or disruption of MoY in XY embryos causes feminization, whereas overexpression of MoY in XX embryos induces masculinization. Crosses between transformed XY females and XX males give rise to males and females, indicating that a Y chromosome can be transmitted by XY females. MoY is Y-linked and functionally conserved in other species of the Tephritidae family, highlighting its potential to serve as a tool for developing more effective control strategies against these major agricultural insect pests.

Housekeeping in Tephritid insects: the best gene choice for expression analyses in the medfly and the olive fly

Real-time quantitative-PCR has been a priceless tool for gene expression analyses. The reaction, however, needs proper normalization with the use of housekeeping genes (HKGs), whose expression remains stable throughout the experimental conditions. Often, the combination of several genes is required for accurate normalization. Most importantly, there are no universal HKGs which can be used since their expression varies among different organisms, tissues or experimental conditions. In the present study, nine common HKGs (RPL19, tbp, ubx, GAPDH, α-TUB, β-TUB, 14-3-3zeta, RPE and actin3) are evaluated in thirteen different body parts, developmental stages and reproductive and olfactory tissues of two insects of agricultural importance, the medfly and the olive fly. Three software programs based on different algorithms were used (geNorm, NormFinder and BestKeeper) and gave different ranking of HKG stabilities. This confirms once again that the stability of common HKGs should not be taken for granted and demonstrates the caution that is needed in the choice of the appropriate HKGs. Finally, by estimating the average of a standard score of the stability values resulted by the three programs we were able to provide a useful consensus key for the choice of the best HKG combination in various tissues of the two insects.

Achilles, a New Family of Transcriptionally Active Retrotransposons from the Olive Fruit Fly, with Y Chromosome Preferential Distribution

Sex chromosomes have many unusual features relative to autosomes. The in depth exploration of their structure will improve our understanding of their origin and divergence (degeneration) as well as the evolution of genetic sex determination pathways which, most often are attributed to them. In Tephritids, the structure of Y chromosome, where the male-determining factor M is localized, is largely unexplored and limited data concerning its sequence content and evolution are available. In order to get insight into the structure and organization of the Y chromosome of the major olive insect pest, the olive fly Bactrocera oleae, we characterized sequences from a Pulse Field Gel Electrophoresis (PFGE)-isolated Y chromosome. Here, we report the discovery of the first olive fly LTR retrotransposon with increased presence on the Y chromosome. The element belongs to the BEL-Pao superfamily, however, its sequence comparison with the other members of the superfamily suggests that it constitutes a new family that we termed Achilles. Its ~7.5 kb sequence consists of the 5'LTR, the 5'non-coding sequence and the open reading frame (ORF), which encodes the polyprotein Gag-Pol. In situ hybridization to the B. oleae polytene chromosomes showed that Achilles is distributed in discrete bands dispersed on all five autosomes, in all centromeric regions and in the granular heterochromatic network corresponding to the mitotic sex chromosomes. The between sexes comparison revealed a variation in Achilles copy number, with male flies possessing 5-10 copies more than female (CI range: 18-38 and 12-33 copies respectively per genome). The examination of its transcriptional activity demonstrated the presence of at least one intact active copy in the genome, showing a differential level of expression between sexes as well as during embryonic development. The higher expression was detected in male germline tissues (testes). Moreover, the presence of Achilles-like elements in different species of the Tephritidae family suggests an ancient origin of this element.

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.

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.

Detection and geographical distribution of the organophosphate resistance-associated Δ3Q ace mutation in the olive fruit fly, Bactrocera oleae (Rossi)

BACKGROUND

The olive fruit fly, Bactrocera oleae (Rossi) (Diptera: Tephritidae), is the most important pest of olives. Its control is based mostly on organophosphate (OP) insecticides, a practice that has led to resistance development. OP resistance in B. oleae has been associated with three mutations in the acetylcholinesterase (AChE), the product of ace gene. The current study presents new diagnostic tests for the detection of the ace mutations and aims at monitoring the frequency of the Δ3Q mutation, which appears associated with resistance at higher OP doses in natural olive fly populations.

RESULTS

An allele-specific polymerase chain reaction (PCR), a PCR-RFLP (restriction fragment length polymorphism) and a Taq-Man test were developed for the Δ3Q mutation detection and a new duplex quantitative PCR assay was designed for the G488S and I214V mutations. Moreover, the frequency of Δ3Q mutation was examined in ten populations of eight countries around the Mediterranean basin. The highest frequencies (10%) were found in Greece and Italy, whereas a gradual decrease of Δ3Q frequency towards the western Mediterranean was noted.

CONCLUSION

Robust tests for insecticide resistance mutations at their incipient levels are essential tools to monitor the increase and geographical spread of such mutations. Three different tests were developed for AChE-Δ3Q that indicated its association with OP applications across the Mediterranean.

Molecular characterization and chromosomal distribution of a species-specific transcribed centromeric satellite repeat from the olive fruit fly, Bactrocera oleae

Satellite repetitive sequences that accumulate in the heterochromatin consist a large fraction of a genome and due to their properties are suggested to be implicated in centromere function. Current knowledge of heterochromatic regions of Bactrocera oleae genome, the major pest of the olive tree, is practically nonexistent. In our effort to explore the repetitive DNA portion of B. oleae genome, a novel satellite sequence designated BoR300 was isolated and cloned. The present study describes the genomic organization, abundance and chromosomal distribution of BoR300 which is organized in tandem, forming arrays of 298 bp-long monomers. Sequence analysis showed an AT content of 60.4%, a CENP-B like-motif and a high curvature value based on predictive models. Comparative analysis among randomly selected monomers demonstrated a high degree of sequence homogeneity (88%-97%) of BoR300 repeats, which are present at approximately 3,000 copies per haploid genome accounting for about 0.28% of the total genomic DNA, based on two independent qPCR approaches. In addition, expression of the repeat was also confirmed through RT-PCR, by which BoR300 transcripts were detected in both sexes. Fluorescence in situ hybridization (FISH) of BoR300 on mitotic metaphases and polytene chromosomes revealed signals to the centromeres of two out of the six chromosomes which indicated a chromosome-specific centromeric localization. Moreover, BoR300 is not conserved in the closely related Bactrocera species tested and it is also absent in other dipterans, but it's rather restricted to the B. oleae genome. This feature of species-specificity attributed to BoR300 satellite makes it a good candidate as an identification probe of the insect among its relatives at early development stages.

Genomic structure, organization and localization of the acetylcholinesterase locus of the olive fruit fly, Bactrocera oleae

Acetylcholinesterase (AChE), encoded by the ace gene, is a key enzyme of cholinergic neurotransmission. Insensitive acetylcholinesterase (AChE) has been shown to be responsible for resistance to OPs and CBs in a number of arthropod species, including the most important pest of olives trees, the olive fruit fly Bactrocera oleae. In this paper, the organization of the B. oleae ace locus, as well as the structural and functional features of the enzyme, are determined. The organization of the gene was deduced by comparison to the ace cDNA sequence of B. oleae and the organization of the locus in Drosophila melanogaster. A similar structure between insect ace gene has been found, with conserved exon-intron positions and junction sequences. The B. oleae ace locus extends for at least 75 kb, consists of ten exons with nine introns and is mapped to division 34 of the chromosome arm IIL. Moreover, according to bioinformatic analysis, the Bo AChE exhibits all the common features of the insect AChE. Such structural and functional similarity among closely related AChE enzymes may implicate similarities in insecticide resistance mechanisms.

Altered GPI modification of insect AChE improves tolerance to organophosphate insecticides

The olive fruit fly Bactrocera oleae is the most destructive and intractable pest of olives. The management of B. oleae has been based on the use of organophosphate (OP) insecticides, a practice that induced resistance. OP-resistance in the olive fly was previously shown to be associated with two mutations in the acetylcholinesterase (AChE) enzyme that, apparently, hinder the entrance of the OP into the active site. The search for additional mutations in the ace gene that encodes AChE revealed a short deletion of three glutamines (Δ3Q) from a stretch of five glutamines, in the C-terminal peptide that is normally cleaved and substituted by a GPI anchor. We verified that AChEs from B. oleae and other Dipterans are actually GPI-anchored, although this is not predicted by the "big-PI" algorithm. The Δ3Q mutation shortens the unusually long hydrophilic spacer that follows the predicted GPI attachment site and may thus improve the efficiency of GPI anchor addition. We expressed the wild type B. oleae AChE, the natural mutant Δ3Q and a constructed mutant lacking all 5 consecutive glutamines (Δ5Q) in COS cells and compared their kinetic properties. All constructs presented identical K(m) and k(cat) values, in agreement with the fact that the mutations did not affect the catalytic domain of the enzyme. In contrast, the mutants produced higher AChE activity, suggesting that a higher proportion of the precursor protein becomes GPI-anchored. An increase in the number of GPI-anchored molecules in the synaptic cleft may reduce the sensitivity to insecticides.

Isolation, annotation and applications of expressed sequence tags from the olive fly, Bactrocera oleae

The olive fruit fly, Bactrocera oleae, is the major pest of the olive tree. Despite its importance, very little genetic and molecular knowledge is available. The present study is a first attempt to identify and characterize B. oleae expressed sequence tags (ESTs). One hundred and ninety-five randomly selected cDNA clones were isolated and the obtained sequences were annotated through BLASTX similarity searches. A set of 159 unique putative transcripts were functionally assigned using Gene Ontology terms in broad categories of biological process, molecular function and cellular component based on D. melanogaster matches. Moreover, the cytogenetic location of 35 ESTs was determined by in situ hybridization to B. oleae polytene chromosomes. The resulting low-resolution EST map more than doubles the available entry points to the insect's genome and can assist syntenic comparisons with other distant species. The deduced codon usage of the isolated ESTs suggested a conserved pattern of B. oleae with its closest relatives. Additionally, the comparative analysis of B. oleae ESTs with the homologous D. melanogaster genes led to the development of 17 nuclear EPIC-PCR markers for the amplification of intron sequences of 11 Tephritidae species. Sequencing analysis of several cross-amplified intron sequences revealed a high degree of conservation among Bactrocera species and a varying transferability of the generated markers across the examined genera, suggesting that this method can provide a useful tool for the clarification of phylogenetic relationships among different species, particularly in cases of species complexes.

Spinosad resistance development in wild olive fruit fly Bactrocera oleae (Diptera: Tephritidae) populations in California

BACKGROUND

Among target pests of the insecticide spinosad is the olive fruit fly, Bactrocera oleae (Rossi) (Diptera: Tephritidae). In Cyprus, spinosad has been sporadically used since its registration in 2002, whereas in Greece its use has been very limited since its registration in 2004, particularly in biological olive cultivars in Crete. By contrast, in California it has been the only insecticide used against the olive fruit fly since its registration in 2004. This study aimed at examining the resistance status of the olive fruit fly to spinosad.

RESULTS

Populations from California, Greece and Cyprus, plus a laboratory population, were tested. Bioassays were performed by oral or topical application of different concentrations of the insecticide. Cypriot populations demonstrated no resistance as compared with that of the laboratory population. Among the Greek populations, only one from Crete demonstrated a fourfold increase in resistance, whereas five populations from California demonstrated a 9-13-fold increase.

CONCLUSION

The observed resistance increase was associated with spinosad applications in the respective areas. These values are relatively low and do not yet pose a serious control problem in the field. However, the observed variation documents that spinosad tolerance has increased in areas where the insecticide has been more extensively used.

Molecular detection of the entomopathogenic bacterium Pseudomonas entomophila using PCR

AIMS

To develop a specific, fast and simple molecular method useful to detect the entomopathogenic bacterium Pseudomonas entomophila.

METHODS AND RESULTS

The use of bioinformatics tools allowed the identification of unique genes present in P. entomophila genome. Using such genes, we designed primers aiming to detect specifically P. entomophila by PCR. Furthermore, a pair of primers specifically designed to amplify the 16S rRNA gene in Pseudomonas species was used. Primer specificity was checked using environmental pseudomonad and nonpseudomonad species. A 618 -bp fragment was amplified only in Pseudomonas using the 16S rDNA primers. Primers (PSEEN1497) designed to detect P. entomophila amplified a 570 -bp fragment only in P. entomophila. A duplex PCR was developed combining 16S rDNA and PSEEN1497 primers that allowed the detection of P. entomophila present in experimentally infected Drosophila melanogaster.

CONCLUSIONS

We developed a molecular method useful to detect P. entomophila present in bacterial cultures or directly from infected insects.

SIGNIFICANCE OF THE STUDY

To the best of our knowledge, this is the first molecular method aiming to detect P. entomophila in environmental samples. The use of our method will facilitate studies related to ecology and insect host range of this entomopathogenic bacterium.

Isolation and characterization of microsatellite markers from the Mediterranean fruit fly, Ceratitis capitata: cross-species amplification in other Tephritidae species reveals a varying degree of transferability

The Mediterranean fruit fly, Ceratitis capitata, is a pest of major economic importance and has become a model for the development of SIT control programs for insect pests. Significant information has been accumulated on classical and population genetics of this species during the past 2 decades. However, the availability of molecular markers is limited. Here, we present the isolation and characterization of 159 microsatellite clones and the development of 108 polymorphic microsatellite markers for this insect pest. Mapping by in situ hybridization to polytene chromosomes of 21 microsatellite clones enriched the cytogenetic map that was previously constructed by our group. The enriched map provides a large number of STSs for future genome mapping projects. Cross-species amplification of these microsatellite loci in 12 Tephritidae species and sequence analysis of several amplification products indicated a varying degree of transferability and their possible usefulness as molecular and genetic markers in these species where genetic and molecular tools are limited.

Analysis of olive fly invasion in California based on microsatellite markers

The olive fruit fly, Bactrocera oleae, is the main pest of the olive fruit and its expansion is exclusively restricted to the cultivation zone of the olive tree. Even though olive production has a century-old history in California, the olive fly was first detected in the Los Angeles area in 1998. Within 5 years of the first observation, the insect was reported from all olive cultivation areas of the state. Field-collected flies from five locations in California and another from Israel were analyzed on the basis of microsatellite polymorphisms in 10 microsatellite loci. These results were integrated with those of a previous study of olive fly populations around the European part of the Mediterranean basin. The analysis pointed to the eastern part of the Mediterranean as the putative source of the observed invasion. Moreover, samples from California were quite different from Mediterranean samples implying the participation of phenomena such as genetic drift during the invasion and expansion of the olive fly in California.

Isolation and characterization of microsatellite markers from the olive fly, Bactrocera oleae, and their cross-species amplification in the Tephritidae family

Background

The Tephritidae family of insects includes the most important agricultural pests of fruits and vegetables, belonging mainly to four genera (Bactrocera, Ceratitis, Anastrepha and Rhagoletis). The olive fruit fly, Bactrocera oleae, is the major pest of the olive fruit. Currently, its control is based on chemical insecticides. Environmentally friendlier methods have been attempted in the past (Sterile Insect Technique), albeit with limited success. This was mainly attributed to the lack of knowledge on the insect's behaviour, ecology and genetic structure of natural populations. The development of molecular markers could facilitate the access in the genome and contribute to the solution of the aforementioned problems. We chose to focus on microsatellite markers due to their abundance in the genome, high degree of polymorphism and easiness of isolation.

Results

Fifty-eight microsatellite-containing clones were isolated from the olive fly, Bactrocera oleae, bearing a total of sixty-two discrete microsatellite motifs. Forty-two primer pairs were designed on the unique sequences flanking the microsatellite motif and thirty-one of them amplified a PCR product of the expected size. The level of polymorphism was evaluated against wild and laboratory flies and the majority of the markers (93.5%) proved highly polymorphic. Thirteen of them presented a unique position on the olive fly polytene chromosomes by in situ hybridization, which can serve as anchors to correlate future genetic and cytological maps of the species, as well as entry points to the genome. Cross-species amplification of these markers to eleven Tephritidae species and sequencing of thirty-one of the amplified products revealed a varying degree of conservation that declines outside the Bactrocera genus.

Conclusion

Microsatellite markers are very powerful tools for genetic and population analyses, particularly in species deprived of any other means of genetic analysis. The presented set of microsatellite markers possesses all features that would render them useful in such analyses. This could also prove helpful for species where SIT is a desired outcome, since the development of effective SIT can be aided by detailed knowledge at the genetic and molecular level. Furthermore, their presented efficacy in several other species of the Tephritidae family not only makes them useful for their analysis but also provides tools for phylogenetic comparisons among them.

A small deletion in the olive fly acetylcholinesterase gene associated with high levels of organophosphate resistance

Organophosphate resistance in the olive fly was previously shown to associate with two point mutations in the ace gene. The frequency of these mutations was monitored in Bactrocera oleae individuals of increasing resistance. In spite of the difference in resistance among the individuals, there was no correlation between mutation frequencies and resistance level, indicating that other factors may contribute to this variation. The search for additional mutations in the ace gene of highly resistant insects revealed a small deletion at the carboxyl terminal of the protein (termed Delta3Q). Significant correlation was shown between the mutation frequency and resistance level in natural populations. In addition, remaining activity of acetylcholinesterase enzyme (AChE) after dimethoate inhibition was higher in genotypes carrying the mutation. These results strongly suggest a role of Delta3Q in high levels of organophosphate (OP) resistance. Interestingly, the carboxyl terminal of AChE is normally cleaved and substituted by a glycosylphosphatidylinositol (GPI) anchor. We hypothesize that Delta3Q may improve GPI anchoring, thus increasing the amount of AChE that reaches the synaptic cleft. In this way, despite the presence of insecticide, enough enzyme would remain in the cleft for its normal role of acetylcholine hydrolysis, allowing the insect to survive. This provides a previously un-described mechanism of resistance.

Organophosphate resistance in olive fruit fly, Bactrocera oleae, populations in Greece and Cyprus

The olive fruit fly Bactrocera oleae (Gmelin) (Diptera: Tephritidae) is the most important pest of olives in countries around the Mediterranean basin. Its control has been based mostly on bait sprays with organophosphate insecticides (usually dimethoate or fenthion) for about 40 years. In the present study, the resistance status of olive fruit fly populations to dimethoate was examined in Greece and Cyprus over 2 years. Thirty-one populations from various regions of Greece, nine from Cyprus and one laboratory susceptible strain, which served as a control, were assayed by topical application of dimethoate. Considerable variation in the resistance levels to dimethoate was recorded in the populations of B. oleae, with resistance ratios ranging from 6.3 to 64.4 (ED(50) values 12.5-128.7 ng dimethoate per insect). The highest resistance ratios were found in populations from Crete, and the lowest in those from Cyprus. This variation could be attributed to different selection pressures from insecticidal applications among populations from the various regions. Migration of resistant genotypes, either autonomous or via commerce, may also be involved.

Microsatellite analysis of olive fly populations in the Mediterranean indicates a westward expansion of the species

Bactrocera oleae is the major insect pest of the olive fruit. Twelve microsatellite loci isolated from the genome of this insect were used in a Mediterranean-wide population analysis. These loci were highly polymorphic with a mean number of alleles per locus of 10.42 and a mean effective number of alleles of 2.76. The analysis was performed on a sample of 671 flies collected from nineteen locations around the European part of the Mediterranean basin. Despite the high level of gene flow across the Mediterranean, results support the notion of a differentiation of three subpopulations: one of the Iberian Peninsula, one of Greece and Italy and one of Cyprus. In addition, the gradual decrease of heterozygosity from the Eastern to the Western part of the Mediterranean indicates a westward expansion of the species.

Identification and partial characterization of a new Ceratitis capitata-specific 44-bp pericentromeric repeat

Tandem satellite DNA repeats are often associated with centromeres. In spite of their importance in the organization of the centromere, they do not seem to be broadly conserved among species and their role is still unclear. Here we report the identification of a new 44-bp tandem pericentromeric repeat from the medfly, Ceratitis capitata. The repeat is specific to this insect and is not found in any of the other closely related species tested. It localizes in four out of its five autosomes and in the X chromosome. It is organized in long arrays, interspersed by transposable elements and other less well-defined sequence motifs.

Genetic markers in ribosomal DNA for the identification of members of the genus Anisakis (Nematoda: ascaridoidea) defined by polymerase-chain-reaction-based restriction fragment length polymorphism

Polymerase-chain-reaction-based restriction fragment length polymorphism analysis was performed to establish genetic markers in rDNA, for the identification of the three sibling species of the Anisakis simplex complex and morphologically differentiated Anisakis species, i.e. Anisakis physeteris, Anisakis schupakovi, Anisakis typica and Anisakis ziphidarum. Different restriction patterns were found between A. simplex sensu stricto and Anisakis pegreffii with two of the restriction endonucleases used (HinfI and TaqI), between A. simplex sensu stricto and A. simplex C with one endonuclease (HhaI), and between A. simplex C and Aniskis pegreffii with three endonucleases (HhaI, HinfI and TaqI), while no variation in patterns was detected among individuals within each species. The species A. physeteris, A. schupakovi, A. typica and A. ziphidarum were found to be different from each other and different from the three sibling species of the A. simplex complex by distinct fragments using 10-12 of the endonucleases tested. The polymorphisms obtained by restriction fragment length polymorphisms have provided a new set of genetic markers for the accurate identification of sibling species and morphospecies.

Genomic structure and ecdysone regulation of the prophenoloxidase 1 gene in the malaria vector Anopheles gambiae

Prophenoloxidase, a melanin-synthesizing enzyme, is considered to be an important arthropod immune protein. In mosquitoes, prophenoloxidase has been shown to be involved in refractory mechanisms against malaria parasites. In our study we used Anopheles gambiae, the most important human malaria vector, to characterize the first arthropod prophenoloxidase gene at the genomic level. The complete nucleotide sequence, including the immediate 5' flanking sequence (-855 bp) of the prophenoloxidase 1 gene, was determined. The gene spans 10 kb and is composed of five exons and four introns coding for a 2.5-kb mRNA. In the 5' flanking sequence, we found several putative regulatory motifs, two of which were identified as ecdysteroid regulatory elements. Electrophoretic mobility gel-shift assays and supershift assays demonstrated that the Aedes aegypti ecdysone receptor/Ultraspiracle nuclear receptor complex, and, seemingly, the endogenous Anopheles gambiae nuclear receptor complex, was able to bind one of the ecdysteroid response elements. Furthermore, 20-hydroxyecdysone stimulation was shown to up-regulate the transcription of the prophenoloxidase 1 gene in an A. gambiae cell line.

Diagnosis of a case of gastric anisakidosis by PCR-based restriction fragment length polymorphism analysis

A set of genetic markers, based on PCR-RFLPs of three diagnostic restriction enzymes (Hhal, Hinfl and Taql), which proved to be suitable for the identification of the species of the genus Anisakis, was used for the first molecular identification of a larva obtained by endoscopy in a case of gastric anisakidosis, in a 51 year old woman from Southern Italy. The analysis of the restriction profiles obtained allowed the larva to be identified as Anisakis pegreffii, one of the three sibling species of the A. simplex complex. PCR-RFLP proved to be a cost-effective and reliable tool for the exact identification of Anisakis larvae recovered from infected humans.

Are chromosomal inversions induced by transposable elements? A paradigm from the malaria mosquito Anopheles gambiae

Chromosomal rearrangements abound in nature and can be studied in detail in organisms with polytene chromosomes. In Drosophila and in Anopheline mosquitoes most speciation processes seem to be associated with the establishment of chromosomal rearrangements, particularly of paracentric inversions. It is not known what triggers inversions in natural populations. In the laboratory inversions are commonly generated by X-rays, mutagens or after the activity of certain transposable elements (TEs). The Anopheles gambiae complex is comprised of six sibling species, each one characterized by the presence of fixed paracentric inversions on their chromosomes. Two of these, An. gambiae s.s. and An. arabiensis, are the most important vectors of human malaria and are structured into sub-populations, each carrying a characteristic set of polymorphic chromosomal inversions. We have cloned the breakpoints of the naturally occurring polymorphic inversion In(2R)d' of An. arabiensis. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a transposable element that we called Odysseus. Characteristics of Odysseus' terminal region and its cytological distribution in different strains as well as within the same strain indicate that Odysseus is an actively transposing element. The presence of Odysseus at the junction of the naturally occurring inversion In(2R)d' suggests that the inversion may be the result of the TEs activity. Cytological evidence from Drosophila melanogaster has also implicated the hobo transposable element in the generation of certain Hawaiian endemic inversions. This picture supports the hypothesis of the important role of TEs in generating natural inversions.

Characterization of the soluble guanylyl cyclase beta-subunit gene in the mosquito Anopheles gambiae

Genomic DNA corresponding to the soluble guanylyl cyclase beta-subunit (GCSbeta) gene was cloned and sequenced from Anopheles gambiae. The sequence was 8103 bp long and presumably included the entire coding region. The deduced amino acid sequence was 71% and 62% similar to previously known Drosophila and vertebrate GCSbeta, while the C-terminus of A. gambiae GCSbeta was shorter. Because of the conserved characteristics in each functional domain, the high G+C% in the third codon positions compared to the introns, the lack of internal stop codons, and the fact that we identified the gene from a cDNA, we conclude that this A. gambiae gene is functional. This is the first detailed description of a guanylyl cyclase gene structure (e.g. intron-exon boundaries). Interestingly, within the fifth intron we found high similarity to the flanking regions of the Pegasus-27 transposable element and other noncoding regions of the A. gambiae genome.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd' of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd' suggests that the inversion may be the result of the transposable element's activity. Characteristics of Odysseus' terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

Molecular cloning and expression of a hexamerin cDNA from the malaria mosquito, Anopheles gambiae

During the last larval instar, dipteran insects synthesize two hexamerins rich in aromatic residues, typified by the larval serum proteins 1 and 2 (LSP-1 and LSP-2) of Drosophila melanogaster. We report here the characterization of a complete cDNA sequence encoding a LSP-1-like protein from a lower dipteran insect, the malaria mosquito Anopheles gambiae. The cDNA encodes the subunit of a homohexamer, A. gambiae hexamerin-1.1 (AgHex-1.1), which is a major pupal protein but only a minor constituent of late larval hemolymph. AgHex-1.1 is moderately rich in methionine (3.9%) and particularly rich in aromatic residues (21% Phe+Tyr). Cytogenetic analysis reveals AgHex-1.1 to be encoded by a single-copy gene localized to division 22F within the proximal 2La inversion breakpoint of chromosome 2 of A. gambiae. The AgHex-1.1 transcript is first detected in fourth-instar larvae (L4) and disappears abruptly in early pupae. In situ hybridization shows accumulation of the transcript uniquely in the larval fat body. AgHex-1.1 mRNA is re-expressed in male and female adults at about 10% of the L4 level, with no effect of bloodfeeding in females. The potential roles of AgHex-1.1 in Anopheles development and reproductive maturation are discussed.

Inversion monophyly in African anopheline malaria vectors

The African Anopheles gambiae complex of six sibling species has many polymorphic and fixed paracentric inversions detectable in polytene chromosomes. These have been used to infer phylogenetic relationships as classically done with Drosophila. Two species, A. gambiae and A. merus, were thought to be sister taxa based on a shared X inversion designated Xag. Recent DNA data have conflicted with this phylogenetic inference as they have supported a sister taxa relationship of A. gambiae and A. arabiensis. A possible explanation is that the Xag is not monophyletic. Here we present data from a gene (soluble guanylate cyclase) within the Xag that strongly supports the monophyly of the Xag. We conjecture that introgression may be occurring between the widely sympatric species A. gambiae and A. arabiensis and that the previous DNA phylogenies have been detecting the introgression. Evidently, introgression is not uniform across the genome, and species-specific regions, like the X-chromosome inversions, do not introgress probably due to selective elimination in hybrids and backcrosses.

Physical map of the malaria vector Anopheles gambiae

Random cDNA clones, cosmid clones and RAPD polymorphic fragments have been localized by in situ hybridization to the ovarian nurse cell polytene chromosomes of the malaria vector Anopheles gambiae. We thus established 85 molecular markers for 110 sites within the whole A. gambiae polytene chromosome complement. The cDNA clones analyzed were isolated at random, and their exact localizations were determined by in situ hybridization. For 15 of the cDNA clones, a partial nucleotide sequence has been obtained; for nine of them sequence searches in the GenBank database revealed high degrees of similarity with published sequences. The cosmid clones analyzed were obtained as the result of screening with a few of the aforementioned cDNA clones of particular interest, or taken from a small set of randomly isolated cosmid clones. The RAPD clones are polymorphic fragments, potentially diagnostic for the various chromosomal forms of A. gambiae that are currently being analyzed.

Distribution of genetic diversity in relation to chromosomal inversions in the malaria mosquito Anopheles gambiae

The epidemiology of malaria in Africa is complicated by the fact that its principal vector, the mosquito Anopheles gambiae, constitutes a complex of six sibling species. Each species is characterized by a unique array of paracentric inversions, as deduced by karyotypic analysis. In addition, most of the species carry a number of polymorphic inversions. In order to develop an understanding of the evolutionary histories of different parts of the genome, we compared the genetic variation of areas inside and outside inversions in two distinct inversion karyotypes of A. gambiae. Thirty-five cDNA clones were mapped on the five arms of the A. gambiae chromosomes with divisional probes. Sixteen of these clones, localized both inside and outside inversions of chromosome 2, were used as probes in order to determine the nucleotide diversity of different parts of the genome in the two inversion karyotypes. We observed that the sequence diversity inside the inversion is more than three-fold lower than in areas outside the inversion and that the degree of divergence increases gradually at loci at increasing distance from the inversion. To interpret the data we present a selectionist and a stochastic model, both of which point to a relatively recent origin of the studied inversion and may suggest differences between the evolutionary history of inversions in Anopheles and Drosophila species.

An anchored restriction-mapping approach applied to the genetic analysis of the Anopheles gambiae malaria vector complex 1

We introduce here a simple approach for rapidly determining restriction maps for a number of regions of a genome; this involves "anchoring" a map with a rare restriction site (in this case the seldom-cutting EagI) followed by partial digestion of a frequent-cutting enzyme (e.g., Sau 3A). We applied this technology to five species of the Anopheles gambiae complex. In a single Southern blot we obtained about a 15-kb restriction map each for the mtDNA, rRNA gene, and a scnDNA region for each of five species. Phylogenetic analyses of these regions yield trees at odds with the more traditional chromosome inversion-based trees. The value of the approach for systematic purposes is the ease with which several large, independent regions of the genome can be quickly assayed for molecular variation.