Comprehensive transcriptome analysis of early male and female Bactrocera jarvisi embryos

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

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

Sexual Dimorphism in Growth Rate and Gene Expression Throughout Immature Development in Wild Type Chrysomya rufifacies (Diptera: Calliphoridae) Macquart

Reliability of forensic entomology analyses to produce relevant information to a given case requires an understanding of the underlying arthropod population(s) of interest and the factors contributing to variability. Common traits for analyses are affected by a variety of genetic and environmental factors. One trait of interest in forensic investigations has been species-specific temperature-dependent growth rates. Recent work indicates sexual dimorphism may be important in the analysis of such traits and related genetic markers of age. However, studying sexual dimorphic patterns of gene expression throughout immature development in wild-type insects can be difficult due to a lack of genetic tools, and the limits of most sex-determination mechanisms. Chrysomya rufifacies, however, is a particularly tractable system to address these issues as it has a monogenic sex determination system, meaning females have only a single-sex of offspring throughout their life. Using modified breeding procedures (to ensure single-female egg clutches) and transcriptomics, we investigated sexual dimorphism in development rate and gene expression. Females develop slower than males (9 h difference from egg to eclosion respectively) even at 30°C, with an average egg-to-eclosion time of 225 h for males and 234 h for females. Given that many key genes rely on sex-specific splicing for the development and maintenance of sexually dimorphic traits, we used a transcriptomic approach to identify different expression of gene splice variants. We find that 98.4% of assembled nodes exhibited sex-specific, stage-specific, to sex-by-stage specific patterns of expression. However, the greatest signal in the expression data is differentiation by developmental stage, indicating that sexual dimorphism in gene expression during development may not be investigatively important and that markers of age may be relatively independent of sex. Subtle differences in these gene expression patterns can be detected as early as 4 h post-oviposition, and 12 of these nodes demonstrate homology with key Drosophila sex determination genes, providing clues regarding the distinct sex determination mechanism of C. rufifacies. Finally, we validated the transcriptome analyses through qPCR and have identified five genes that are developmentally informative within and between sexes.

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

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

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

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

Evidence of Adaptive Evolution and Relaxed Constraints in Sex-Biased Genes of South American and West Indies Fruit Flies (Diptera: Tephritidae)

Several studies have demonstrated that genes differentially expressed between sexes (sex-biased genes) tend to evolve faster than unbiased genes, particularly in males. The reason for this accelerated evolution is not clear, but several explanations have involved adaptive and nonadaptive mechanisms. Furthermore, the differences of sex-biased expression patterns of closely related species are also little explored out of Drosophila. To address the evolutionary processes involved with sex-biased expression in species with incipient differentiation, we analyzed male and female transcriptomes of Anastrepha fraterculus and Anastrepha obliqua, a pair of species that have diverged recently, likely in the presence of gene flow. Using these data, we inferred differentiation indexes and evolutionary rates and tested for signals of selection in thousands of genes expressed in head and reproductive transcriptomes from both species. Our results indicate that sex-biased and reproductive-biased genes evolve faster than unbiased genes in both species, which is due to both adaptive pressure and relaxed constraints. Furthermore, among male-biased genes evolving under positive selection, we identified some related to sexual functions such as courtship behavior and fertility. These findings suggest that sex-biased genes may have played important roles in the establishment of reproductive isolation between these species, due to a combination of selection and drift, and unveil a plethora of genetic markers useful for more studies in these species and their differentiation.

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

Background

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

Results

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

Conclusion

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

The study of the transformer gene from Bactrocera dorsalis and B. correcta with putative core promoter regions

Background

The transformer (tra) is a sex determining switch in different orders of insects, including Diptera, as in the family Tephritidae. The lifelong autoregulatory loop of tra female-specific splicing can be reset by the intervention of male-specific primary signals (M factor). In early development, the functional female and truncated male TRA proteins relay the sexual fates to the alternative splicing of a bisexual switch gene, doublesex (dsx) cascading the sexual differentiation processes. Bactrocera dorsalis (Hendel) and Bactrocera correcta (Bezzi) are among the Bactrocera model worldwide key pests. Area-wide integrated pest management using the male-only Sterile Insect Technique (SIT) relying on genetic sexing systems is effective in control programs. We undertook the molecular characterization and comparative studies of the tra orthologues in the Bactrocera species, including the Salaya1 genetic sexing strain (GSS).

Results

RT-PCR revealed that B. dorsalis tra (Bdtra) and B. correcta tra (Bctra) transcripts contained conservation of both constitutive exons and male-specific exons as in other Bactrocera. However, new Bdtra male-specific exons were retained, diversifying the pattern of the male-specifically spliced transcripts. The coding sequences of tra were highly conserved in Bactrocera (86–95 %) but less so among related genera (61–65 %) within the same Tephritidae family. A conservation of deduced amino acid sequences (18 residues), called the TEP region, was identified to be distinctive among tephritids. The 5’ regulatory sequence containing many structural characteristics of the putative core promoter was discovered in B. correcta. The expression patterns of Bdtra and Bctra were sex-specifically spliced and the signals relayed to the dsx genes in the adult wild-types. However, the coexistence of male- and female-specifically spliced transcripts (980 and 626 bp, respectively) of the B. dorsalis wild-type strain was found in the Salaya1 GSS adult males. The Bdtra RNA interference masculinized the XX karyotype females into pseudomales, but their testes were mostly not well developed.

Conclusions

Bdtra and Bctra have sex-specific splicing, similar to Bactroceras, Ceratitis capitata (Wiedemann), and Anastrephas. A newly identified TEP region is proposed in tephritids. A putative core promoter has been discovered in Bctra.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0342-0) contains supplementary material, which is available to authorized users.

The Microbiome of Field-Caught and Laboratory-Adapted Australian Tephritid Fruit Fly Species with Different Host Plant Use and Specialisation

Tephritid fruit fly species display a diversity of host plant specialisation on a scale from monophagy to polyphagy. Furthermore, while some species prefer ripening fruit, a few are restricted to damaged or rotting fruit. Such a diversity of host plant use may be reflected in the microbial symbiont diversity of tephritids and their grade of dependency on their microbiomes. Here, we investigated the microbiome of six tephritid species from three genera, including species that are polyphagous pests (Bactrocera tryoni, Bactrocera neohumeralis, Bactrocera jarvisi, Ceratitis capitata) and a monophagous specialist (Bactrocera cacuminata). These were compared with the microbiome of a non-pestiferous but polyphagous tephritid species that is restricted to damaged or rotting fruit (Dirioxa pornia). The bacterial community associated with whole fruit flies was analysed by 16S ribosomal DNA (rDNA) amplicon pyrosequencing to detect potential drivers of taxonomic composition. Overall, the dominant bacterial families were Enterobacteriaceae and Acetobacteraceae (both Proteobacteria), and Streptococcaceae and Enterococcaceae (both Firmicutes). Comparisons across species and genera found different microbial composition in the three tephritid genera, but limited consistent differentiation between Bactrocera species. Within Bactrocera species, differentiation of microbial composition seemed to be influenced by the environment, possibly including their diets; beyond this, tephritid species identity or ecology also had an effect. The microbiome of D. pornia was most distinct from the other five species, which may be due to its ecologically different niche of rotting or damaged fruit, as opposed to ripening fruit favoured by the other species. Our study is the first amplicon pyrosequencing study to compare the microbiomes of tephritid species and thus delivers important information about the turnover of microbial diversity within and between fruit fly species and their potential application in pest management strategies.

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

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