Sex and the single embryo: early deveiopment in the Mediterranean fruit fly, Ceratitis capitata

Gene drive and genetic sex conversion in the global agricultural pest Ceratitis capitata

Homing-based gene drives are recently proposed interventions promising the area-wide, species-specific genetic control of harmful insect populations. Here we characterise a first set of gene drives in a tephritid agricultural pest species, the Mediterranean fruit fly Ceratitis capitata (medfly). Our results show that the medfly is highly amenable to homing-based gene drive strategies. By targeting the medfly transformer gene, we also demonstrate how CRISPR-Cas9 gene drive can be coupled to sex conversion, whereby genetic females are transformed into fertile and harmless XX males. Given this unique malleability of sex determination, we modelled gene drive interventions that couple sex conversion and female sterility and found that such approaches could be effective and tolerant of resistant allele selection in the target population. Our results open the door for developing gene drive strains for the population suppression of the medfly and related tephritid pests by co-targeting female reproduction and shifting the reproductive sex ratio towards males. They demonstrate the untapped potential for gene drives to tackle agricultural pests in an environmentally friendly and economical way.

In toto light sheet fluorescence microscopy live imaging datasets of Ceratitis capitata embryonic development

The Mediterranean fruit fly (medfly), Ceratitis capitata, is an important model organism in biology and agricultural research with high economic relevance. However, information about its embryonic development is still sparse. We share nine long-term live imaging datasets acquired with light sheet fluorescence microscopy (484.5 h total recording time, 373 995 images, 256 Gb) with the scientific community. Six datasets show the embryonic development in toto for about 60 hours at 30 minutes intervals along four directions in three spatial dimensions, covering approximately 97% of the entire embryonic development period. Three datasets focus on germ cell formation and head involution. All imaged embryos hatched morphologically intact. Based on these data, we suggest a two-level staging system that functions as a morphogenetic framework for upcoming studies on medfly. Our data supports research on wild-type or aberrant morphogenesis, quantitative analyses, comparative approaches to insect development as well as studies related to pest control. Further, they can be used to test advanced image processing approaches or to train machine learning algorithms and/or neuronal networks.

Early embryonic transcriptomes of Zeugodacus tau provide insight into sex determination and differentiation genes

Zeugodacus tau (Walker) is an invasive pest. The sterile insect technique is an environment-friendly method for pest control. Understanding the mechanism of sex determination will contribute to improving efficiency of this technique. In this study, we identified the transformer (tra) gene in Z. tau. One female-specific and two male-specific isoforms of tra were found in Z. tau, and the male-specific splicing pattern of tra was found to occur 5 h after egg laying. We performed transcriptome sequencing at 1 h (E1), 5 h (E5), and 9 h (E9) after egg laying and obtained high-quality transcriptome libraries of early embryo development. We identified 13 297 and 11 713 differentially expressed genes (DEGs) from E5 versus E1 and E9 versus E1 comparisons, respectively. To explore the potential functions of the DEGs during embryonic development, Gene Ontology, Clusters of Orthologous Groups of proteins, and Kyoto Encyclopedia of Genes and Genomes analyses were performed. Twenty-six genes potentially involved in sex determination or differentiation, including Maleness-on-the-Y (MoY), were identified in Z. tau. To verify the transcriptome results, 15 genes were selected for quantitative real-time PCR validation. The results were consistent with the transcriptome sequencing results. Moreover, U2 small nuclear riboprotein auxiliary factor (U2AF-50), female lethal d (fl(2)d), and virilizer (vir) were highly expressed at E5, indicating that they may be related to the sex-specific splicing of tra. Further functional analysis is needed to confirm this speculation. Our data provide an insight into the mechanism underlying sex determination and differentiation in tephritid species.

Conditional knockdown of transformer in sheep blow fly suggests a role in repression of dosage compensation and potential for population suppression

The transformer (tra) gene is essential for female development in many insect species, including the Australian sheep blow fly, Lucilia cuprina. Sex-specific tra RNA splicing is controlled by Sex lethal (Sxl) in Drosophila melanogaster but is auto-regulated in L. cuprina. Sxl also represses X chromosome dosage compensation in female D. melanogaster. We have developed conditional Lctra RNAi knockdown strains using the tet-off system. Four strains did not produce females on diet without tetracycline and could potentially be used for genetic control of L. cuprina. In one strain, which showed both maternal and zygotic tTA expression, most XX transformed males died at the pupal stage. RNAseq and qRT-PCR analyses of mid-stage pupae showed increased expression of X-linked genes in XX individuals. These results suggest that Lctra promotes somatic sexual differentiation and inhibits X chromosome dosage compensation in female L. cuprina. However, XX flies homozygous for a loss-of-function Lctra knockin mutation were fully transformed and showed high pupal eclosion. Two of five X-linked genes examined showed a significant increase in mRNA levels in XX males. The stronger phenotype in the RNAi knockdown strain could indicate that maternal Lctra expression may be essential for initiation of dosage compensation suppression in female embryos.

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.

Transcribed sex-specific markers on the Y chromosome of the oriental fruit fly, Bactrocera dorsalis

BACKGROUND

The Oriental fruit fly, Bactrocera dorsalis, is a highly polyphagous invasive species with a high reproductive potential. In many tropical and subtropical parts of the world it ranks as one of the major pests of fruits and vegetables. Due to its economic importance, genetic, cytogenetic, genomic and biotechnological approaches have been applied to understand its biology and to implement the Sterile Insect Technique, currently a part of area-wide control programmes against this fly. Its chromosome complement includes five pairs of autosomes and the sex chromosomes. The X and Y sex chromosomes are heteromorphic and the highly heterochromatic and degenerate Y harbours the male factor BdMoY. The characterization of the Y chromosome in this fly apart from elucidating its role as primary sex determination system, it is also of crucial importance to understand its role in male biology. The repetitive nature of the Y chromosome makes it challenging to sequence and characterise.

RESULTS

Using Representational Difference Analysis, fluorescent in situ hybridisation on mitotic chromosomes and in silico genome resources, we show that the B. dorsalis Y chromosome harbours transcribed sequences of gyf, (typo-gyf) a homologue of the Drosophila melanogaster Gigyf gene, and of a non-LTR retrotransposon R1. Similar sequences are also transcribed on the X chromosome. Paralogues of the Gigyf gene are also present on the Y and X chromosomes of the related species B. tryoni. Another identified Y-specific repetitive sequence linked to BdMoY appears to be specific to B. dorsalis.

CONCLUSIONS

Our random scan of the Y chromosome provides a broad picture of its general composition and represents a starting point for further applicative and evolutionary studies. The identified repetitive sequences can provide a useful Y-marking system for molecular karyotyping of single embryos. Having a robust diagnostic marker associated with BdMoY will facilitate studies on how BdMoY regulates the male sex determination cascade during the embryonic sex-determination window. The Y chromosome, despite its high degeneracy and heterochromatic nature, harbours transcribed sequences of typo-gyf that may maintain their important function in post-transcriptional mRNA regulation. That transcribed paralogous copies of Gigyf are present also on the X and that this genomic distribution is maintained also in B. tryoni raises questions on the evolution of sex chromosomes in Bactrocera and other tephritids.

Transcriptome analysis of Anastrepha fraterculus sp. 1 males, females, and embryos: insights into development, courtship, and reproduction

BACKGROUND

Anastrepha fraterculus sp. 1 is considered a quarantine pest in several American countries. Since chemical control applied in an integrated pest management program is the only strategy utilized against this pest, the development of pesticide-free methods, such as the Sterile Insect Technique, is being considered. The search for genes involved in sex-determination and differentiation, and in metabolic pathways associated with communication and mating behaviour, contributes with key information to the development of genetic control strategies. The aims of this work were to perform a comprehensive analysis of A. fraterculus sp. 1 transcriptome and to obtain an initial evaluation of genes associated with main metabolic pathways by the expression analysis of specific transcripts identified in embryos and adults.

RESULTS

Sexually mature adults of both sexes and 72 h embryos were considered for transcriptome analysis. The de novo transcriptome assembly was fairly complete (62.9% complete BUSCO orthologs detected) with a total of 86,925 transcripts assembled and 28,756 GO annotated sequences. Paired-comparisons between libraries showed 319 transcripts differently expressed between embryos and females, 1242 between embryos and males, and 464 between sexes. Using this information and genes searches based on published studies from other tephritid species, we evaluated a set of transcripts involved in development, courtship and metabolic pathways. The qPCR analysis evidenced that the early genes serendipity alpha and transformer-2 displayed similar expression levels in the analyzed stages, while heat shock protein 27 is over-expressed in embryos and females in comparison to males. The expression of genes associated with courtship (takeout-like, odorant-binding protein 50a1) differed between males and females, independently of their reproductive status (virgin vs mated individuals). Genes associated with metabolic pathways (maltase 2-like, androgen-induced gene 1) showed differential expression between embryos and adults. Furthermore, 14,262 microsatellite motifs were identified, with 11,208 transcripts containing at least one simple sequence repeat, including 48% of di/trinucleotide motifs.

CONCLUSION

Our results significantly expand the available gene space of A. fraterculus sp. 1, contributing with a fairly complete transcript database of embryos and adults. The expression analysis of the selected candidate genes, along with a set of microsatellite markers, provides a valuable resource for further genetic characterization of A. fraterculus sp. 1 and supports the development of specific genetic control strategies.

The transformer-2 and fruitless characterisation with developmental expression profiles of sex-determining genes in Bactrocera dorsalis and B. correcta

Sex determination in tephritid fruit flies involves a signaling cascade of alternatively spliced genes. The Transformer (TRA) and Transformer-2 (TRA-2) complex establishes an autoregulatory loop switching sex-specific splicing of tra pre-mRNA in females. The TRA/TRA-2 complex also regulates the sex-specific splicing of downstream effector genes, doublesex (dsx) and fruitless (fru). In Ceratitis capitata, a Maleness-on the-Y (MoY) gene modulates sex-specifically spliced Cctra pre-mRNA and results in the breakdown of the Cctra autoregulatory loop in males. In this study, the tra-2 and fru genes were characterised in two key pests, Bactrocera dorsalis and B. correcta. The tra-2 genes showed high degrees of conservation among tephritids. The complex gene organisation for each of Bdfru and Bcfru were identified. There are sex-specific and non sex-specific transcripts generated by alternative promoters as found in Drosophila melanogaster and other insects. RNAi knockdown of Bdtra transcripts showed that BdTRA controls the sex-specific splicing of Bddsx and Bdfru pre-mRNAs. Developmental expression analysis shows that multiple splice variants of Bdtra and Bctra RNAs are present before and during cellular blastoderm formation and that the mature sex-specific variants become fixed later in embryogenesis. Furthermore, the Bddsx^M splice variants are found in early embryos at the beginning of gastulation, but Bdfru^M does not appear until the larval stage. We proposed that the zygotic tra loop is initiated in both female and male embryos before becoming automatised or abolished by MoY, respectively.

Characterization of the First W-Specific Protein-Coding Gene for Sex Identification in Helicoverpa armigera

Helicoverpa armigera is a globally-important crop pest with a WZ (female)/ZZ (male) sex chromosome system. The absence of discernible sexual dimorphism in its egg and larval stages makes it impossible to address any sex-related theoretical and applied questions before pupation unless a W-specific sequence marker is available for sex diagnosis. To this end, we used one pair of morphologically pre-sexed pupae to PCR-screen 17 non-transposon transcripts selected from 4855 W-linked candidate reads identified by mapping a publicly available egg transcriptome of both sexes to the male genome of this species and detected the read SRR1015458.67499 only in the female pupa. Subsequent PCR screenings of this read and the previously reported female-specific RAPD (random amplified polymorphic DNA) marker AF18 with ten more pairs of pre-sexed pupae and different annealing positions and/or temperatures as well as its co-occurrence with the female-specific transcript splicing isoforms of doublesex gene of H. armigera (Hadsx) and amplification and sequencing of their 5′ unknown flanking sequences in three additional pairs of pre-sexed pupae verified that SRR1015458.67499 is a single copy protein-coding gene unique to W chromosome (named GUW1) while AF18 is a multicopy MITE transposon located on various chromosomes. Test application of GUW1 as a marker to sex 30 neonates of H. armigera yielded a female/male ratio of 1.14: 1.00. Both GUW1 and Hadsx splicing isoforms assays revealed that the H. armigera embryo cell line QB-Ha-E-1 is a male cell line. Taken together, GUW1 is not only a reliable DNA marker for sexing all stages of H. armigera and its cell lines, but also represents the first W-specific protein-coding gene in lepidopterans.

Transcriptome Analysis of the Oriental Fruit Fly Bactrocera dorsalis Early Embryos

The oriental fruit fly, Bactrocera dorsalis (Hendel), is one of the most devastating and highly invasive agricultural pests world-wide, resulting in severe economic loss. Thus, it is of great interest to understand the transcriptional changes that occur during the activation of its zygotic genome at the early stages of embryonic development, especially the expression of genes involved in sex determination and the cellularization processes. In this study, we applied Illumina sequencing to identify B. dorsalis sex determination genes and early zygotic genes by analyzing transcripts from three early embryonic stages at 0–1, 2–4, and 5–8 h post-oviposition, which include the initiation of sex determination and cellularization. These tests generated 13,489 unigenes with an average length of 2185 bp. In total, 1683, 3201 and 3134 unigenes had significant changes in expression levels at times after oviposition including at 2–4 h versus 0–1 h, 5–8 h versus 0–1 h, and 5–8 h versus 2–4 h, respectively. Clusters of gene orthology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were performed throughout embryonic development to better understand the functions of differentially expressed unigenes. We observed that the RNA binding and spliceosome pathways were highly enriched and overrepresented during the early stage of embryogenesis. Additionally, transcripts for 21 sex-determination and three cellularization genes were identified, and expression pattern analysis revealed that the majority of these genes were highly expressed during embryogenesis. This study is the first assembly performed for B. dorsalis based on Illumina next-generation sequencing technology during embryogenesis. Our data should contribute significantly to the fundamental understanding of sex determination and early embryogenesis in tephritid fruit flies, and provide gene promoter and effector gene candidates for transgenic pest-management strategies for these economically important species.

miRNA-1-3p is an early embryonic male sex-determining factor in the Oriental fruit fly Bactrocera dorsalis

Regulation of male sexual differentiation by a Y chromosome-linked male determining factor (M-factor) is one of a diverse array of sex determination mechanisms found in insects. By deep sequencing of small RNAs from Bactrocera dorsalis early embryos, we identified an autosomal-derived microRNA, miR-1-3p, that has predicted target sites in the transformer gene (Bdtra) required for female sex determination. We further demonstrate by both in vitro and in vivo tests that miR-1-3p suppresses Bdtra expression. Injection of a miR-1-3p mimic in early embryos results in 87-92% phenotypic males, whereas knockdown of miR-1-3p by an inhibitor results in 67-77% phenotypic females. Finally, CRISPR/Cas9-mediated knockout of miR-1-3p results in the expression of female-specific splice variants of Bdtra and doublesex (Bddsx), and induced sex reversal of XY individuals into phenotypic females. These results indicate that miR-1-3p is required for male sex determination in early embryogenesis in B. dorsalis as an intermediate male determiner.

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.

The autoregulatory loop: A common mechanism of regulation of key sex determining genes in insects

Sex determination in most insects is structured as a gene cascade, wherein a primary signal is passed through a series of sex-determining genes, culminating in a downstream double-switch known as doublesex that decides the sexual fate of the embryo. From the literature available on sex determination cascades, it becomes apparent that sex determination mechanisms have evolved rapidly. The primary signal that provides the cue to determine the sex of the embryo varies remarkably, not only among taxa, but also within taxa. Furthermore, the upstream key gene in the cascade also varies between species and even among closely related species. The order Insecta alone provides examples of astoundingly complex diversity of upstream key genes in sex determination mechanisms. Besides, unlike key upstream genes, the downstream double-switch gene is alternatively spliced to form functional sex-specific isoforms. This sex-specific splicing is conserved across insect taxa. The genes involved in the sex determination cascade such as Sex-lethal (Sxl) in Drosophila melanogaster, transformer (tra) in many other dipterans, coleopterans and hymenopterans, Feminizer (fem) in Apis mellifera, and IGF-II mRNA-binding protein (Bmimp) in Bombyx mori are reported to be regulated by an autoregulatory positive feedback loop. In this review, by taking examples from various insects, we propose the hypothesis that autoregulatory loop mechanisms of sex determination might be a general strategy. We also discuss the possible reasons for the evolution of autoregulatory loops in sex determination cascades and their impact on binary developmental choices.

Isolation and Molecular Characterization of the Transformer Gene From Bactrocera cucurbitae (Diptera: Tephritidae)

transformer (tra) is a switch gene of sex determination in many insects, particularly in Dipterans. However, the sex determination pathway in Bactrocera cucurbitae (Coquillett), a very destructive pest on earth, remains largely uncharacterized. In this study, we have isolated and characterized one female-specific and two male-specific transcripts of the tra gene (Bcutra) of B. cucurbitae. The genomic structure of Bcutra has been determined and the presence of multiple conserved Transformer (TRA)/TRA-2 binding sites in Bcutra has been found. BcuTRA is highly conservative with its homologues in other tephritid fruit flies. Gene expression analysis of Bcutra at different developmental stages demonstrates that the female transcript of Bcutra appears earlier than the male counterparts, indicating that the maternal TRA is inherited in eggs and might play a role in the regulation of TRA expression. The conservation of protein sequence and sex-specific splicing of Bcutra and its expression patterns during development suggest that Bcutra is probably the master gene of sex determination of B. cucurbitae. Isolation of Bcutra will facilitate the development of a genetic sexing strain for its biological control.

Delivery of Nucleic Acids through Embryo Microinjection in the Worldwide Agricultural Pest Insect, Ceratitis capitata

The Mediterranean fruit fly (medfly) Ceratitis capitata (Wiedemann) (Diptera: Tephritidae) is a pest species with extremely high agricultural relevance. This is due to its reproductive behavior: females damage the external surface of fruits and vegetables when they lay eggs and the hatched larvae feed on their pulp. Wild C. capitata populations are traditionally controlled through insecticide spraying and/or eco-friendly approaches, the most successful being the Sterile Insect Technique (SIT). The SIT relies on mass-rearing, radiation-based sterilization and field release of males that retain their capacity to mate but are not able to generate fertile progeny. The advent and the subsequent rapid development of biotechnological tools, together with the availability of the medfly genome sequence, has greatly boosted our understanding of the biology of this species. This favored the proliferation of new strategies for genome manipulation, which can be applied to population control. In this context, embryo microinjection plays a dual role in expanding the toolbox for medfly control. The ability to interfere with the function of genes that regulate key biological processes, indeed, expands our understanding of the molecular machinery underlying medfly invasiveness. Furthermore, the ability to achieve germ-line transformation facilitates the production of multiple transgenic strains that can be tested for future field applications in novel SIT settings. Indeed, genetic manipulation can be used to confer desirable traits that can, for example, be used to monitor sterile male performance in the field, or that can result in early life-stage lethality. Here we describe a method to microinject nucleic acids into medfly embryos to achieve these two main goals.

The whole genome sequence of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), reveals insights into the biology and adaptive evolution of a highly invasive pest species

Background

The Mediterranean fruit fly (medfly), Ceratitis capitata, is a major destructive insect pest due to its broad host range, which includes hundreds of fruits and vegetables. It exhibits a unique ability to invade and adapt to ecological niches throughout tropical and subtropical regions of the world, though medfly infestations have been prevented and controlled by the sterile insect technique (SIT) as part of integrated pest management programs (IPMs). The genetic analysis and manipulation of medfly has been subject to intensive study in an effort to improve SIT efficacy and other aspects of IPM control.

Results

The 479 Mb medfly genome is sequenced from adult flies from lines inbred for 20 generations. A high-quality assembly is achieved having a contig N50 of 45.7 kb and scaffold N50 of 4.06 Mb. In-depth curation of more than 1800 messenger RNAs shows specific gene expansions that can be related to invasiveness and host adaptation, including gene families for chemoreception, toxin and insecticide metabolism, cuticle proteins, opsins, and aquaporins. We identify genes relevant to IPM control, including those required to improve SIT.

Conclusions

The medfly genome sequence provides critical insights into the biology of one of the most serious and widespread agricultural pests. This knowledge should significantly advance the means of controlling the size and invasive potential of medfly populations. Its close relationship to Drosophila, and other insect species important to agriculture and human health, will further comparative functional and structural studies of insect genomes that should broaden our understanding of gene family evolution.

Electronic supplementary material

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

The role of the transformer gene in sex determination and reproduction in the tephritid fruit fly, Bactrocera dorsalis (Hendel)

Transformer (tra) is a switch gene in the somatic sex-determination hierarchy that regulates sexual dimorphism based on RNA splicing in many insects. In tephritids, a Y-linked male determining gene (M) controls sex in the sex-determination pathway. Here, homologues of Drosophila tra and transformer-2 (tra-2) genes were isolated and characterized in Bactrocera dorsalis (Hendel), one of the most destructive agricultural insect pests in many Asian countries. Two male-specific and one female-specific isoforms of B. dorsalis transformer (Bdtra) were identified. The presence of multiple TRA/TRA-2 binding sites in Bdtra suggests that the TRA/TRA-2 proteins are splicing regulators promoting and maintaining, epigenetically, female sex determination by a tra positive feedback loop in XX individuals during development. The expression patterns of female-specific Bdtra transcripts during early embryogenesis shows that a peak appears at 15 h after egg laying. Using dsRNA to knock-down Bdtra expression in the embryo and adult stages, we showed that sexual formation is determined early in the embryo stage and that parental RNAi does not lead to the production of all male progeny as in Tribolium castaneum. RNAi results from adult abdominal dsRNA injections show that Bdtra has a positive influence on female yolk protein gene (Bdyp1) expression and fecundity.

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.

RNAi-Mediated Knock-Down of transformer and transformer 2 to Generate Male-Only Progeny in the Oriental Fruit Fly, Bactrocera dorsalis (Hendel)

The transformer (tra) gene appears to act as the genetic switch that promotes female development by interaction with the transformer2 (tra-2) gene in several dipteran species including the Medfly, housefly and Drosophila melanogaster. In this study, we describe the isolation, expression and function of tra and tra-2 in the economically important agricultural pest, the oriental fruit fly, Bactrocera dorsalis (Hendel). Bdtra and Bdtra-2 are similar to their homologs from other tephritid species. Bdtra demonstrated sex-specific transcripts: one transcript in females and two transcripts in males. In contrast, Bdtra-2 only had one transcript that was common to males and females, which was transcribed continuously in different adult tissues and developmental stages. Bdtra-2 and the female form of Bdtra were maternally inherited in eggs, whereas the male form of Bdtra was not detectable until embryos of 1 and 2 h after egg laying. Function analyses of Bdtra and Bdtra-2 indicated that both were indispensable for female development, as nearly 100% males were obtained with embryonic RNAi against either Bdtra or Bdtra-2. The fertility of these RNAi-generated males was subsequently tested. More than 80% of RNAi-generated males could mate and the mated females could lay eggs, but only 40-48.6% males gave rise to progeny. In XX-reversed males and intersex individuals, no clear female gonadal morphology was observed after dissection. These results shed light on the development of a genetic sexing system with male-only release for this agricultural pest.

De novo assembly and transcriptome analysis of the Mediterranean fruit fly Ceratitis capitata early embryos

The agricultural pest Ceratitis capitata, also known as the Mediterranean fruit fly or Medfly, belongs to the Tephritidae family, which includes a large number of other damaging pest species. The Medfly has been the first non-drosophilid fly species which has been genetically transformed paving the way for designing genetic-based pest control strategies. Furthermore, it is an experimentally tractable model, in which transient and transgene-mediated RNAi have been successfully used. We applied Illumina sequencing to total RNA preparations of 8–10 hours old embryos of C. capitata, This developmental window corresponds to the blastoderm cellularization stage. In summary, we assembled 42,614 transcripts which cluster in 26,319 unique transcripts of which 11,045 correspond to protein coding genes; we identified several hundreds of long ncRNAs; we found an enrichment of transcripts encoding RNA binding proteins among the highly expressed transcripts, such as CcTRA-2, known to be necessary to establish and, most likely, to maintain female sex of C. capitata. Our study is the first de novo assembly performed for Ceratitis capitata based on Illumina NGS technology during embryogenesis and it adds novel data to the previously published C. capitata EST databases. We expect that it will be useful for a variety of applications such as gene cloning and phylogenetic analyses, as well as to advance genetic research and biotechnological applications in the Medfly and other related Tephritidae.

Comprehensive transcriptome analysis of early male and female Bactrocera jarvisi embryos

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.

Subtractive and differential hybridization molecular analyses of Ceratitis capitata XX/XY versus XX embryos to search for male-specific early transcribed genes

The agricultural pest Ceratitis capitata, also known as the Mediterranean fruit fly or Medfly, is a fruit crop pest of very high economic relevance in different continents. The strategy to separate Ceratitis males from females (sexing) in mass rearing facilities is a useful step before the sterilization and release of male-only flies in Sterile Insect Technique control programs (SIT). The identification of genes having early embryonic male-specific expression, including Y-linked genes, such as the Maleness factor, could help to design novel and improved methods of sexing in combination with transgenesis, aiming to confer conditional female-specific lethality or female-to-male sexual reversal.

We used a combination of Suppression Subtractive Hybrydization (SSH), Mirror Orientation Selection (MOS) and differential screening hybridization (DSH) techniques to approach the problem of isolating corresponding mRNAs expressed in XX/XY embryos versus XX-only embryos during a narrow developmental window (8-10 hours after egg laying, AEL ). Here we describe a novel strategy we have conceived to obtain relatively large amounts of XX-only embryos staged at 8-10 h AEL and so to extract few micrograms of polyA+ required to apply the complex technical procedure. The combination of these 3 techniques led to the identification of a Y-linked putative gene, CcGm2, sharing high sequence identity to a paralogous gene, CcGm1, localized either on an autosome or on the X chromosome.

We propose that CcGm2 is a first interesting putative Y-linked gene which could play a role in sex determination. The function exterted by this gene should be investigated by novel genetic tools, such as CRISPR-CAS9, which will permit to target only the Y-linked paralogue, avoiding to interfere with the autosomal or X-linked paralogue function.

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

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

Expression patterns of sex-determination genes in single male and female embryos of two Bactrocera fruit fly species during early development

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.

Conservation and sex-specific splicing of the transformer gene in the calliphorids Cochliomyia hominivorax, Cochliomyia macellaria and Lucilia sericata

Transformer (TRA) promotes female development in several dipteran species including the Australian sheep blowfly Lucilia cuprina, the Mediterranean fruit fly, housefly and Drosophila melanogaster. tra transcripts are sex-specifically spliced such that only the female form encodes full length functional protein. The presence of six predicted TRA/TRA2 binding sites in the sex-specific female intron of the L. cuprina gene suggested that tra splicing is auto-regulated as in medfly and housefly. With the aim of identifying conserved motifs that may play a role in tra sex-specific splicing, here we have isolated and characterized the tra gene from three additional blowfly species, L. sericata, Cochliomyia hominivorax and C. macellaria. The blowfly adult male and female transcripts differ in the choice of splice donor site in the first intron, with males using a site downstream of the site used in females. The tra genes all contain a single TRA/TRA2 site in the male exon and a cluster of four to five sites in the male intron. However, overall the sex-specific intron sequences are poorly conserved in closely related blowflies. The most conserved regions are around the exon/intron junctions, the 3′ end of the intron and near the cluster of TRA/TRA2 sites. We propose a model for sex specific regulation of tra splicing that incorporates the conserved features identified in this study. In L. sericata embryos, the male tra transcript was first detected at around the time of cellular blastoderm formation. RNAi experiments showed that tra is required for female development in L. sericata and C. macellaria. The isolation of the tra gene from the New World screwworm fly C. hominivorax, a major livestock pest, will facilitate the development of a “male-only” strain for genetic control programs.

Transgenic sexing system for Ceratitis capitata (Diptera: Tephritidae) based on female-specific embryonic lethality

Fruit fly pest species have been successfully controlled and managed via the Sterile Insect Technique (SIT), a control strategy that uses infertile matings of sterile males to wild females to reduce pest populations. Biological efficiency in the field is higher if only sterile males are released in SIT programs and production costs are also reduced. Sexing strains developed in the Mediterranean fruit fly Ceratitis capitata (medfly) through classical genetics are immensely beneficial to medfly SIT programs but exhibit reduced fertility and fitness. Moreover, transfer of such classical genetic systems to other tephritid species is difficult. Transgenic approaches can overcome this limitation of classical genetic sexing strains (GSSs), but had resulted so far in transgenic sexing strains (TSSs) with dominant lethality at late larval and pupal stages. Here we present a transgene-based female-specific lethality system for early embryonic sexing in medfly. The system utilizes the sex-specifically spliced transformer intron to restrict ectopic mRNA translation of the pro-apoptotic gene hid(Ala5) to females only. The expression of this lethal effector gene is driven by a tetracycline-repressible transactivator gene tTA that is under the control of promoters/enhancers of early-acting cellularization genes. Despite observed position effects on the sex-specific splicing, we could effectively establish this early-acting transgenic sexing system in the medfly C. capitata. After satisfactory performance in large scale tests, TSSs based on this system will offer cost-effective sexing once introduced into SIT programs. Moreover, this approach is straight forward to be developed also for other insect pest and vector species.

Transcriptional profiles of mating-responsive genes from testes and male accessory glands of the Mediterranean fruit fly, Ceratitis capitata

BACKGROUND

Insect seminal fluid is a complex mixture of proteins, carbohydrates and lipids, produced in the male reproductive tract. This seminal fluid is transferred together with the spermatozoa during mating and induces post-mating changes in the female. Molecular characterization of seminal fluid proteins in the Mediterranean fruit fly, Ceratitis capitata, is limited, although studies suggest that some of these proteins are biologically active.

METHODOLOGY/PRINCIPAL FINDINGS

We report on the functional annotation of 5914 high quality expressed sequence tags (ESTs) from the testes and male accessory glands, to identify transcripts encoding putative secreted peptides that might elicit post-mating responses in females. The ESTs were assembled into 3344 contigs, of which over 33% produced no hits against the nr database, and thus may represent novel or rapidly evolving sequences. Extraction of the coding sequences resulted in a total of 3371 putative peptides. The annotated dataset is available as a hyperlinked spreadsheet. Four hundred peptides were identified with putative secretory activity, including odorant binding proteins, protease inhibitor domain-containing peptides, antigen 5 proteins, mucins, and immunity-related sequences. Quantitative RT-PCR-based analyses of a subset of putative secretory protein-encoding transcripts from accessory glands indicated changes in their abundance after one or more copulations when compared to virgin males of the same age. These changes in abundance, particularly evident after the third mating, may be related to the requirement to replenish proteins to be transferred to the female.

CONCLUSIONS/SIGNIFICANCE

We have developed the first large-scale dataset for novel studies on functions and processes associated with the reproductive biology of Ceratitis capitata. The identified genes may help study genome evolution, in light of the high adaptive potential of the medfly. In addition, studies of male recovery dynamics in terms of accessory gland gene expression profiles and correlated remating inhibition mechanisms may permit the improvement of pest management approaches.

Molecular evolution of Drosophila Sex-lethal and related sex determining genes

Background

Sex determining mechanisms are evolutionarily labile and related species often use different primary signals and gene regulatory networks. This is well illustrated by the sex determining cascade of Drosophila fruitflies, which have recruited Sex-lethal as the master switch and cellular memory of sexual identity, a role performed in other insects by the gene transformer. Here we investigate the evolutionary change in the coding sequences of sex determining genes associated with the recruitment of Sex-lethal. We analyze sequences of Sex-lethal itself, its Drosophila paralogue sister-or-Sex-lethal and downstream targets transformer and doublesex.

Results

We find that the recruitment of sister-or-Sex-lethal was associated with a number of adaptive amino acid substitutions, followed by a tightening of purifying selection within the Drosophila clade. Sequences of the paralogue sister-or-Sex-lethal, in contrast, show a signature of rampant positive selection and relaxation of purifying selection. The recruitment of Sex-lethal as top regulator and memory gene is associated with a significant release from purifying selection in transformer throughout the Drosophila clade. In addition, doublesex shows a signature of positive selection and relaxation of purifying selection in the Drosophila clade. A similar pattern is seen in sequences from the sister Tephritidae clade.

Conclusions

The pattern of molecular evolution we observe for Sex-lethal and its paralogue sister-or-Sex-lethal is not characteristic of a duplication followed by neo-functionalization. Rather, evidence suggests a sub-functionalization scenario achieved through the evolution of sophisticated splicing. As expected, we find that transformer evolves under relaxed purifying selection after the recruitment of Sex-lethal in Drosophila. Finally, the observation of doublesex adaptation in both Drosophila and Tephritidae suggests that these changes are due to ongoing adaptation of downstream sex-specific regulation, rather than being associated the recruitment of Sex-lethal and the resulting change in the topology of the sex determining cascade.

Sex determination in insects: a binary decision based on alternative splicing

The gene regulatory networks that control sex determination vary between species. Despite these differences, comparative studies in insects have found that alternative splicing is reiteratively used in evolution to control expression of the key sex-determining genes. Sex determination is best understood in Drosophila where activation of the RNA binding protein-encoding gene Sex-lethal is the central female-determining event. Sex-lethal serves as a genetic switch because once activated it controls its own expression by a positive feedback splicing mechanism. Sex fate choice in is also maintained by self-sustaining positive feedback splicing mechanisms in other dipteran and hymenopteran insects, although different RNA binding protein-encoding genes function as the binary switch. Studies exploring the mechanisms of sex-specific splicing have revealed the extent to which sex determination is integrated with other developmental regulatory networks.

Interchromosomal duplications on the Bactrocera oleae Y chromosome imply a distinct evolutionary origin of the sex chromosomes compared to Drosophila

BACKGROUND

Diptera have an extraordinary variety of sex determination mechanisms, and Drosophila melanogaster is the paradigm for this group. However, the Drosophila sex determination pathway is only partially conserved and the family Tephritidae affords an interesting example. The tephritid Y chromosome is postulated to be necessary to determine male development. Characterization of Y sequences, apart from elucidating the nature of the male determining factor, is also important to understand the evolutionary history of sex chromosomes within the Tephritidae. We studied the Y sequences from the olive fly, Bactrocera oleae. Its Y chromosome is minute and highly heterochromatic, and displays high heteromorphism with the X chromosome.

METHODOLOGY/PRINCIPAL FINDINGS

A combined Representational Difference Analysis (RDA) and fluorescence in-situ hybridization (FISH) approach was used to investigate the Y chromosome to derive information on its sequence content. The Y chromosome is strewn with repetitive DNA sequences, the majority of which are also interdispersed in the pericentromeric regions of the autosomes. The Y chromosome appears to have accumulated small and large repetitive interchromosomal duplications. The large interchromosomal duplications harbour an importin-4-like gene fragment. Apart from these importin-4-like sequences, the other Y repetitive sequences are not shared with the X chromosome, suggesting molecular differentiation of these two chromosomes. Moreover, as the identified Y sequences were not detected on the Y chromosomes of closely related tephritids, we can infer divergence in the repetitive nature of their sequence contents.

CONCLUSIONS/SIGNIFICANCE

The identification of Y-linked sequences may tell us much about the repetitive nature, the origin and the evolution of Y chromosomes. We hypothesize how these repetitive sequences accumulated and were maintained on the Y chromosome during its evolutionary history. Our data reinforce the idea that the sex chromosomes of the Tephritidae may have distinct evolutionary origins with respect to those of the Drosophilidae and other Dipteran families.

Genomic organization and splicing evolution of the doublesex gene, a Drosophila regulator of sexual differentiation, in the dengue and yellow fever mosquito Aedes aegypti

Background

In the model system Drosophila melanogaster, doublesex (dsx) is the double-switch gene at the bottom of the somatic sex determination cascade that determines the differentiation of sexually dimorphic traits. Homologues of dsx are functionally conserved in various dipteran species, including the malaria vector Anopheles gambiae. They show a striking conservation of sex-specific regulation, based on alternative splicing, and of the encoded sex-specific proteins, which are transcriptional regulators of downstream terminal genes that influence sexual differentiation of cells, tissues and organs.

Results

In this work, we report on the molecular characterization of the dsx homologue in the dengue and yellow fever vector Aedes aegypti (Aeadsx). Aeadsx produces sex-specific transcripts by alternative splicing, which encode isoforms with a high degree of identity to Anopheles gambiae and Drosophila melanogaster homologues. Interestingly, Aeadsx produces an additional novel female-specific splicing variant. Genomic comparative analyses between the Aedes and Anopheles dsx genes revealed a partial conservation of the exon organization and extensive divergence in the intron lengths. An expression analysis showed that Aeadsx transcripts were present from early stages of development and that sex-specific regulation starts at least from late larval stages. The analysis of the female-specific untranslated region (UTR) led to the identification of putative regulatory cis-elements potentially involved in the sex-specific splicing regulation. The Aedes dsx sex-specific splicing regulation seems to be more complex with the respect of other dipteran species, suggesting slightly novel evolutionary trajectories for its regulation and hence for the recruitment of upstream splicing regulators.

Conclusions

This study led to uncover the molecular evolution of Aedes aegypti dsx splicing regulation with the respect of the more closely related Culicidae Anopheles gambiae orthologue. In Aedes aegypti, the dsx gene is sex-specifically regulated and encodes two female-specific and one male-specific isoforms, all sharing a doublesex/mab-3 (DM) domain-containing N-terminus and different C-termini. The sex-specific regulation is based on a combination of exon skipping, 5' alternative splice site choice and, most likely, alternative polyadenylation. Interestingly, when the Aeadsx gene is compared to the Anopheles dsx ortholog, there are differences in the in silico predicted default and regulated sex-specific splicing events, which suggests that the upstream regulators either are different or act in a slightly different manner. Furthermore, this study is a premise for the future development of transgenic sexing strains in mosquitoes useful for sterile insect technique (SIT) programs.

The transformer gene of Ceratitis capitata: a paradigm for a conserved epigenetic master regulator of sex determination in insects

The transformer gene in Ceratitis capitata (Cctra(ep)) is the founding member of a family of related SR genes that appear to act as the master epigenetic switch in sex determination in insects. A functional protein seems to be produced only in individuals with a female XX karyotype where it is required to maintain the productive mode of expression through a positive feedback loop and to direct female development by instructing the downstream target genes accordingly. When zygotic activation of this loop is prevented, male development follows. Recently, tra(ep) orthologues were isolated in more distantly related dipteran species including Musca domestica, Glossina morsitans and Lucilia cuprina and in the Hymenopterans Apis mellifera and Nasonia vitripennis. All of these tra(ep) orthologues seem to act as binary switches that govern all aspects of sexual development. Transient silencing leads to complete masculinization of individuals with a female karyotype. Reciprocally, in some systems it has been shown that transient expression of the functional TRA product is sufficient to transactivate the endogenous gene and implement female development in individuals with a male karyotype. Hence, a mechanism based on tra(ep) epigenetic autoregulation seems to represent a common and presumably ancestral single principle of sex determination in Insecta. The results of these studies will not only be important for understanding divergent evolution of basic developmental processes but also for designing new strategies to improve genetic sexing in different insect species of economical or medical importance.