Search

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

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

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.

Three Australian tephritid fruit flies (Bactrocera tryoni - Q-fly, Bactrocera neohumeralis - NEO, and Bactrocera jarvisi - JAR) are promising models for genetic studies of pest status and invasiveness. The long history of ecological and physi- ological studies of the three species has been augmented by the development of a range of genetic and genomic tools, including the capacity for forced multigeneration crosses between the three species followed by selection experiments, a draft genome for Q-fly, and tissue- and stage-specific transcriptomes. The Q-fly and NEO species pair is of particular interest. The distribution of NEO is contained entirely within the wider distribution of Q-fly and the two species are ecologically extremely similar, with no known differences in pheromones, temperature tolerance, or host-fruit utilisation. However there are three clear differences be- tween them: humeral callus colour, complete pre-mating isolation based on mating time-of-day, and invasiveness. NEO is much less invasive, whereas in historical times Q-fly has invaded southeastern Australia and areas of Western Australia and the North- ern Territory. In southeastern fruit-growing regions, microsatellites suggest that some of these outbreaks might derive from ge- netically differentiated populations overwintering in or near the invaded area. Q-fly and NEO show very limited genome diffe- rentiation, so comparative genomic analyses and QTL mapping should be able to identify the regions of the genome controlling mating time and invasiveness, to assess the genetic bases for the invasive strains of Q-fly, and to facilitate a variety of improve- ments to current sterile insect control strategies for that species (Current Zoology 61 (3): 477-487, 2015). Keywords Bactrocera, Genomics, STERILE-insect-technique, Interspecific-hybridisation

In tephritids, the sex-determination pathway follows the sex-specific splicing of transformer (tra) mRNA, and the cooperation of tra and transformer-2 (tra-2) to effect the sex-specific splicing of doublesex (dsx), the genetic double-switch responsible for male or female somatic development. The Dominant Male Determiner (M) is the primary signal that controls this pathway. M, as yet uncharacterized, is Y-chromosome linked, expressed in the zygote and directly or indirectly diminishes active TRA protein in male embryos. Here we first demonstrated the high conservation of tra, tra-2 and dsx in two Australian tephritids, Bactrocera tryoni and Bactrocera jarvisi. We then used quantitative reverse transcription PCR on single, sexed embryos to examine expression of the key sex-determination genes during early embryogenesis. Individual embryos were sexed using molecular markers located on the B. jarvisi Y-chromosome that was also introgressed into a B. tryoni line. In B. jarvisi, sex-specific expression of tra transcripts occurred between 3 to 6 h after egg laying, and the dsx isoform was established by 7 h. These milestones were delayed in B. tryoni lines. The results provide a time frame for transcriptomic analyses to identify M and its direct targets, plus information on genes that may be targeted for the development of male-only lines for pest management.

BackgroundMaternally inheritedWolbachiabacteria infect many insect species. They can also be transferred horizontally into uninfected host lineages. AWolbachiaspillover from an infected source population must occur prior to the establishment of heritable infections, but this spillover may be transient. In a previous study of tephritid fruit fly species of tropical Australia we detected a high incidence of identicalWolbachiastrains in several species as well asWolbachiapseudogenes. Here, we have investigated this further by analysing field specimens of 24 species collected along a 3,000 km climate gradient of eastern Australia.ResultsWolbachiasequences were detected in individuals of nine of the 24 (37 %) species. Seven (29 %) species displayed four distinctWolbachiastrains based on characterisation of full multi locus sequencing (MLST) profiles; the strains occurred as single and double infections in a small number of individuals (2–17 %). For the two remaining species all individuals had incomplete MLST profiles andWolbachiapseudogenes that may be indicative of lateral gene transfer into host genomes. The detection ofWolbachiawas restricted to northern Australia, including in five species that only occur in the tropics. Within the more widely distributedBactrocera tryoniandBactrocera neohumeralis,Wolbachiaalso only occurred in the north, and was not linked to any particular mitochondrial haplotypes.ConclusionsThe presence ofWolbachiapseudogenes at high prevalence in two species in absence of complete MLST profiles may represent footprints of historic infections that have been lost. The detection of identical low prevalence strains in a small number of individuals of seven species may question their role as reproductive manipulator and their vertical inheritance. Instead, the findings may be indicative of transient infections that result from spillover events from a yet unknown source. These spillover events appear to be restricted to northern Australia, without proliferation in host lineages further south. Our study highlights that tropical fruit fly communities containWolbachiapseudogenes and may be exposed to frequent horizontalWolbachiatransfer. It also emphasises that global estimates ofWolbachiafrequencies may need to consider lateral gene transfer andWolbachiaspillover that may be regionally restricted, transient and not inherited.

The DNA sequence of the coding region of the wild-type and mutant alleles of the white eye gene from the melon fly, Bactrocera cucurbitae Coquillett, was obtained. The mutant white eye allele had a single base pair mutation at the 5′ end of intron 4b in the RNA splice recognition site. Due to the disrupted splice recognition site, intron 4b was not removed from the mutant RNA transcript. The resulting mRNA transcript was 68 bp longer than that of the wild type, containing a frameshift and premature stop codon. Transport of pigment precursors would be blocked, resulting in a lack of pigmentation deposition consistent with the known mutant phenotype.

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.

Background Developing embryos are provided with maternal RNA transcripts and proteins, but transcription from the zygotic nuclei must be activated to control continuing embryonic development. Transcripts are generated at different stages of early development, and those involved in sex determination and cellularisation are some of the earliest to be activated. The male sex in tephritid fruit flies is determined by the presence of a Y chromosome, and it is believed that a transcript from the Y-chromosome sets in motion a cascade that determines male development, as part of the greater maternal to zygotic transition (MTZ). Here we investigate the poly(A+) transcriptome in early male and female embryos of the horticultural pest Bactrocera jarvisi (Diptera: Tephritidae). Results Bactrocera jarvisi embryos were collected over two pre-blastoderm time periods, 2-3h and 3-5h after egg laying. Embryos were individually sexed using a Y-chromosome marker, allowing the sex-specific poly(A+) transcriptome of single-sex embryo pools to be deep-sequenced and assembled de novo. Transcripts for sixteen sex-determination and two cellularisation gene homologues of Drosophila melanogaster (Diptera: Drosophilidae) were identified in early embryos of B. jarvisi, including transcripts highly upregulated prior to cellularisation. No strong candidates for transcripts derived solely from the Y chromosome were recovered from the poly(A+) fraction. Conclusions Bactrocera jarvisi provides an excellent model for embryonic studies due to available Y-chromosome markers and the compact time frame for zygotic transcription and the sex-determined state. Our data contribute fundamental information to sex-determination research, and provide candidates for the sourcing of gene promoters for transgenic pest-management strategies of tephritid fruit flies.

Numts are an integral component of many eukaryote genomes offering a snapshot of the evolutionary process that led from the incorporation of an α-proteobacterium into a larger eukaryotic cell some 1.8 billion years ago. Although numt sequence can be harnessed as molecular marker, these sequences often remain unidentified and are mistaken for genuine mtDNA leading to erroneous interpretation of mtDNA data sets. It is therefore indispensable that during the process of amplifying and sequencing mitochondrial genes, preventive measures are taken to ensure the exclusion of numts to guarantee the recovery of genuine mtDNA. This applies to mtDNA analyses in general but especially to studies where mtDNAs are sequenced de novo as the launch pad for subsequent mtDNA-based research. By using a combination of dilution series and nested rolling circle amplification (RCA), we present a novel strategy to selectively amplify mtDNA and exclude the amplification of numt sequence. We have successfully applied this strategy to de novo sequence the mtDNA of the Black Field Cricket Teleogryllus commodus, a species known to contain numts. Aligning our assembled sequence to the reference genome of Teleogryllus emma (GenBank EU557269.1) led to the identification of a numt sequence in the reference sequence. This unexpected result further highlights the need of a reliable and accessible strategy to eliminate this source of error.

Background The tephritid fruit flies include a number of economically important pests of horticulture, with a large accumulated body of research on their biology and control. Amongst the Tephritidae, the genus Bactrocera, containing over 400 species, presents various species groups of potential utility for genetic studies of speciation, behaviour or pest control. In Australia, there exists a triad of closely-related, sympatric Bactrocera species which do not mate in the wild but which, despite distinct morphologies and behaviours, can be force-mated in the laboratory to produce fertile hybrid offspring. To exploit the opportunities offered by genomics, such as the efficient identification of genetic loci central to pest behaviour and to the earliest stages of speciation, investigators require genomic resources for future investigations. Results We produced a draft de novo genome assembly of Australia’s major tephritid pest species, Bactrocera tryoni. The male genome (650 -700 Mbp) includes approximately 150Mb of interspersed repetitive DNA sequences and 60Mb of satellite DNA. Assessment using conserved core eukaryotic sequences indicated 98% completeness. Over 16,000 MAKER-derived gene models showed a large degree of overlap with other Dipteran reference genomes. The sequence of the ribosomal RNA transcribed unit was also determined. Unscaffolded assemblies of B. neohumeralis and B. jarvisi were then produced; comparison with B. tryoni showed that the species are more closely related than any Drosophila species pair. The similarity of the genomes was exploited to identify 4924 potentially diagnostic indels between the species, all of which occur in non-coding regions. Conclusions This first draft B. tryoni genome resembles other dipteran genomes in terms of size and putative coding sequences. For all three species included in this study, we have identified a comprehensive set of non-redundant repetitive sequences, including the ribosomal RNA unit, and have quantified the major satellite DNA families. These genetic resources will facilitate the further investigations of genetic mechanisms responsible for the behavioural and morphological differences between these three species and other tephritids. We have also shown how whole genome sequence data can be used to generate simple diagnostic tests between very closely-related species where only one of the species is scaffolded. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-1153) contains supplementary material, which is available to authorized users.