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

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.

The chromosomes of Drosophila suzukii (Diptera: Drosophilidae): detailed photographic polytene chromosomal maps and in situ hybridization data

The spotted wing drosophila, D. suzukii, is a serious agricultural pest attacking a variety of soft fruits and vegetables. Although originating from East Asia it has recently invaded America and Europe raising major concern about its expansion potential and the consequent economic losses. Since cytogenetic information on the species is scarce, we report here the mitotic karyotype and detailed photographic maps of the salivary gland polytene chromosomes of D. suzukii. The mitotic metaphase complement contains three pairs of autosomes, one of which is dot-like, and one pair of heteromorphic (XX/XY) sex chromosomes. The salivary gland polytene complement consists of five long polytene arms, representing the two metacentric autosomes and the acrocentric X chromosome, and one very short polytene element, which corresponds to the dot-like autosome. Banding pattern as well as the most characteristic features and prominent landmarks of each polytene chromosome arm are presented and discussed. Furthermore, twelve gene markers have been mapped on the polytene chromosomes of D. suzukii by in situ hybridization. Their distribution pattern was found quite similar to that of D. melanogaster revealing conservation of synteny although the relative position within each chromosome arm for most of the genes differed significantly between D. suzukii and D. melanogaster. The chromosome information presented here is suitable for comparative cytogenetic studies and phylogenetic exploration, while it could also facilitate the assembly of the genome sequence and support the development of genetic tools for species-specific and environment-friendly biological control applications such as the sterile insect technique.

Comparative Gene Mapping as a Tool to Understand the Evolution of Pest Crop Insect Chromosomes

The extent of the conservation of synteny and gene order in aphids has been previously investigated only by comparing a small subset of linkage groups between the pea aphid Acyrthosiphon pisum and a few other aphid species. Here we compared the localization of eight A. pisum scaffolds (covering more than 5 Mb and 83 genes) in respect to the Drosophila melanogaster Muller elements identifying orthologous loci spanning all the four A. pisum chromosomes. Comparison of the genetic maps revealed a conserved synteny across different loci suggesting that the study of the fruit fly Muller elements could favour the identification of chromosomal markers useful for the study of chromosomal rearrangements in aphids. A. pisum is the first aphid species to have its genome sequenced and the finding that there are several chromosomal regions in synteny between Diptera and Hemiptera indicates that the genomic tools developed in A. pisum will be broadly useful not only for the study of other aphids but also for other insect species.

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

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

Cytogenetic Analysis of the South American Fruit Fly Anastrepha fraterculus (Diptera:Tephritidae) Species Complex: Construction of Detailed Photographic Polytene Chromosome Maps of the Argentinian Af. sp.1 Member

Genetic and cytogenetic studies constitute a significant basis for understanding the biology of insect pests and the design and the construction of genetic tools for biological control strategies. Anastrepha fraterculus is an important pest of the Tephritidae family. It is distributed from southern Texas through eastern Mexico, Central America and South America causing significant crop damage and economic losses. Currently it is considered as a species complex; until now seven members have been described based on multidisciplinary approaches. Here we report the cytogenetic analysis of an Argentinian population characterized as Af. sp.1 member of the Anastrepha fraterculus species complex. The mitotic karyotype and the first detailed photographic maps of the salivary gland polytene chromosomes are presented. The mitotic metaphase complement consists of six (6) pairs of chromosomes, including one pair of heteromorphic sex chromosomes, with the male being the heterogametic sex. The analysis of the salivary gland polytene complement shows a total number of five long chromosomes that correspond to the five autosomes of the mitotic karyotype and a heterochromatic network corresponding to the sex chromosomes. Comparison of the polytene chromosome maps between this species and Anastrepha ludens shows significant similarity. The polytene maps presented here are suitable for cytogenetic studies that could shed light on the species limits within this species complex and support the development of genetic tools for sterile insect technique (SIT) applications.

Extraordinary conservation of entire chromosomes in insects over long evolutionary periods

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

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 Bactrocera oleae genome: localization of nine genes on the polytene chromosomes of the olive fruit fly (Diptera: Tephritidae)

Four homologous and five heterologous gene-specific sequences have been mapped by in situ hybridization on the salivary gland polytene chromosomes of the olive fruit fly, Bactrocera oleae. The nine genes were dispersed on four of the five autosomal chromosomes, thus enriching the available set of chromosome landmarks for this major agricultural pest. Present data further supports the proposed chromosome homologies among B. oleae, Ceratitis capitata, and Drosophila melanogaster and the idea of the conservation of chromosomal element identity throughout dipteran evolution.

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.

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.

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

The olive fruit fly Bactrocera oleae has a unique ability to cope with olive flesh, and is the most destructive pest of olives worldwide. Its control has been largely based on the use of chemical insecticides, however, the selection of insecticide resistance against several insecticides has evolved. The study of detoxification mechanisms, which allow the olive fruit fly to defend against insecticides, and/or phytotoxins possibly present in the mesocarp, has been hampered by the lack of genomic information in this species. In the NCBI database less than 1,000 nucleotide sequences have been deposited, with less than 10 detoxification gene homologues in total. We used 454 pyrosequencing to produce, for the first time, a large transcriptome dataset for B. oleae. A total of 482,790 reads were assembled into 14,204 contigs. More than 60% of those contigs (8,630) were larger than 500 base pairs, and almost half of them matched with genes of the order of the Diptera. Analysis of the Gene Ontology (GO) distribution of unique contigs, suggests that, compared to other insects, the assembly is broadly representative for the B. oleae transcriptome. Furthermore, the transcriptome was found to contain 55 P450, 43 GST-, 15 CCE- and 18 ABC transporter-genes. Several of those detoxification genes, may putatively be involved in the ability of the olive fruit fly to deal with xenobiotics, such as plant phytotoxins and insecticides. In summary, our study has generated new data and genomic resources, which will substantially facilitate molecular studies in B. oleae, including elucidation of detoxification mechanisms of xenobiotic, as well as other important aspects of olive fruit fly biology.

Sex chromosomes and associated rDNA form a heterochromatic network in the polytene nuclei of Bactrocera oleae (Diptera: Tephritidae)

The olive fruit fly, Bactrocera oleae, has a diploid set of 2n = 12 chromosomes including a pair of sex chromosomes, XX in females and XY in males, but polytene nuclei show only five polytene chromosomes, obviously formed by five autosome pairs. Here we examined the fate of the sex chromosomes in the polytene complements of this species using fluorescence in situ hybridization (FISH) with the X and Y chromosome-derived probes, prepared by laser microdissection of the respective chromosomes from mitotic metaphases. Specificity of the probes was verified by FISH in preparations of mitotic chromosomes. In polytene nuclei, both probes hybridized strongly to a granular heterochromatic network, indicating thus underreplication of the sex chromosomes. The X chromosome probe (in both female and male nuclei) highlighted most of the granular mass, whereas the Y chromosome probe (in male nuclei) identified a small compact body of this heterochromatic network. Additional hybridization signals of the X probe were observed in the centromeric region of polytene chromosome II and in the telomeres of six polytene arms. We also examined distribution of the major ribosomal DNA (rDNA) using FISH with an 18S rDNA probe in both mitotic and polytene chromosome complements of B. oleae. In mitotic metaphases, the probe hybridized exclusively to the sex chromosomes. The probe signals localized a discrete rDNA site at the end of the short arm of the X chromosome, whereas they appeared dispersed over the entire dot-like Y chromosome. In polytene nuclei, the rDNA was found associated with the heterochromatic network representing the sex chromosomes. Only in nuclei with preserved nucleolar structure, the probe signals were scattered in the restricted area of the nucleolus. Thus, our study clearly shows that the granular heterochromatic network of polytene nuclei in B. oleae is formed by the underreplicated sex chromosomes and associated rDNA.

Cytogenetic analysis of the Ethiopian fruit fly Dacus ciliatus (Diptera: Tephritidae)

The Ethiopian fruit fly, Dacus ciliatus, is an important pest of cucurbits, which recently invaded the Middle East. The genetics and cytogenetics of D. ciliatus have been scarcely studied. Such information is, however, an essential basis for understanding the biology of insect pests, as well as for the design of modern control strategies. We report here the mitotic karyotype and detailed photographic maps of the salivary gland polytene chromosomes of this species. The mitotic metaphase complement consists of six pairs of chromosomes, including one pair of heteromorphic sex (XX/XY) chromosomes. The heterogametic sex is ascribed to the male. The analysis of the salivary gland polytene complement shows a total number of five long chromosomes (10 polytene arms), which correspond to the five autosomes of the mitotic nuclei, and a heterochromatic mass corresponding to the sex chromosomes. Banding patterns, as well as the most characteristic features and prominent landmarks of each polytene chromosome are presented and discussed. Chromosomal homologies between D. ciliatus and Bactrocera oleae are proposed by comparing chromosome banding patterns and by in situ hybridization of the hsp70 gene.