In embryos the maternal-to-zygotic transition (MTZ) integrates post-transcriptional regulation of maternal transcripts with transcriptional activation of the zygotic genome. Although the molecular mechanisms underlying this event are being clarified in Drosophila melanogaster, little is know about the embryogenic processes in other insect species. The recent publication of expressed sequence tags (ESTs) from embryos of the global pest species Ceratitis capitata (medfly) has enabled the investigation of embryogenesis in this species and has allowed a comparison of the embryogenic processes in these two related dipteran species, C. capitata and D. melanogaster, that shared a common ancestor 80-100 mya.
Using a novel PCR-based sexing method, which takes advantage of a putative LTR retrotransposon MITE insertion on the medfly Y chromosome, the transcriptomes of individual early male and female embryos were analysed using RT-PCR. This study is focused on two crucial aspects of the onset of embryonic development: sex determination and cellular blastoderm formation. Together with the three known medfly genes (Cctransformer, Cctransformer2 and Ccdoublesex), the expression patterns of other medfly genes that are similar to the D. melanogaster sex-determination genes (sisterlessA, groucho, deadpan, Sex-lethal, female lethal d, sans fille and intersex) and four cellular blastoderm formation genes (Rho1, spaghetti squash, slow-as-molasses and serendipity-α) were analyzed, allowing us to sketch a preliminary outline of the embryonic process in the medfly. Furthermore, a putative homologue of the Zelda gene has been considered, which in D. melanogaster encodes a DNA-binding factor responsible for the maternal-to-zygotic transition.
Our novel sexing method facilitates the study of i) when the MTZ transition occurs in males and females of C. capitata, ii) when and how the maternal information of "female-development" is reprogrammed in the embryos and iii) similarities and differences in the regulation of gene expression in C. capitata and D. melanogaster. We suggest a new model for the onset of the sex determination cascade in the medfly: the maternally inherited Cctra transcripts in the female embryos are insufficient to produce enough active protein to inhibit the male mode of Cctra splicing. The slow rate of development and the inefficiency of the splicing mechanism in the pre-cellular blastoderm facilitates the male-determining factor (M) activity, which probably acts by inhibiting CcTRA protein activity.
The onset of embryonic development is the result of a complex interaction between maternal and zygotic genetic information: during the maternal-to-zygotic transition (MTZ) embryos integrate post-transcriptional regulation of maternal transcripts with transcriptional activation of the zygotic genome .
In the model insect Drosophila melanogaster maternal transcripts and proteins are deposited in the egg during oogenesis and two waves of zygotic gene transcription occur in the embryos: the first between 1 and 2 h of development and the second major burst between 2 and 3 h of development . For at least one third of zygotically active genes, the expression profile is coupled to the degradation of the corresponding maternal mRNA. It has been hypothesized that the degradation of maternal mRNA could be due to the activity of Pumilio, of Smaug [3, 4] and of miRNAs  or to the presence of a general mediator of mRNA degradation in the mRNA itself, i.e. the AU-rich element . The activation of the earliest zygotic genes is linked to the presence of TAGteam sites in their promoter regions; the zinc-finger protein Zelda (zelda, zld) binds to these sites and triggers the expression of the downstream genes. TAGteam sites have been found in the promoters of genes involved in early developmental processes, including cellular blastoderm formation, sex-determination and the dorsoventral-expression pattern .
Although the molecular mechanisms involved in the MTZ transition during the first stages of embryogenesis are being clarified in D. melanogaster , little is known about embryogenic processes including sex determination in other insect species. This lack of knowledge does not allow us to determine whether there are features common to all embryogenic processes.
The recent publication of ESTs from medfly embryos has facilitated the investigation of embryogenesis in this species . The medfly is a highly invasive agricultural pest species and belongs to the family Tephritidae that is estimated to have diverged from the family Drosophilidae between 80 and 100 mya [8, 9]. Even though the medfly genome is more than three times larger than that of D. melanogaster, a recent analysis has shown that functionally non-coding sequences are conserved between these two species and that homologous genes, such as the even-skipped gene (eve), have the same expression profile in both species . Thus, it is reasonable to hypothesize that medfly embryogenesis may follow a similar sequence of events as in D. melanogaster. If that is the case, C. capitata could be used as a reference species to study the molecular mechanisms that trigger embryogenesis, as a big genome facilitates the identification of cis-regulator elements through comparative genomic methods .
It is known that in the medfly the cellular blastoderm formation occurs between 9 and 11 h after oviposition  and one of the major events during early embryogenesis, the molecular mechanism that controls sex-determination, is well studied in C. capitata . In order to detect sex-specific differences in the timing of gene expression and/or mRNA splicing in early developmental stages in the medfly, a primary requirement is the identification of the sex of the early life stages. Male and female embryos are not morphologically distinguishable ; however the sexes may be recognized in the early life stages using PCR-based approaches. These sexing assays are based on the repetitive nature of the medfly Y chromosome. Several repetitive DNA sequences have been isolated in medfly that are Y-specific or enriched on the Y chromosome . Copies of these repetitive sequences, which appear to occur in tandemly linked arrays, are distributed over 90% of the length of the long arm of the Y chromosome . The basic structure of the repeats is similar: each unit (1.3-1.7 kb) contains an AT-rich (83-87%) internal region (200 bp or more), called the A-T element. Although three previous approaches [16–18], based on the repetitive nature of the Y chromosome, may be suitable for determining the sex of embryos, an improved method has been developed and is reported in this paper. Using this molecular sexing approach it is possible to study the sex-specific expression of genes in early embryos; we have addressed four principal points: i) when the MTZ transition occurs in C. capitata; ii) when the sexual fate is established at the molecular level; iii) how the maternal information of "female-development" is reprogrammed in the embryos; and iv) how similar or different is the regulation of gene expression in early embryogenesis between C. capitata and D. melanogaster.
Individual samples of C. capitata from five laboratory strains and one natural population were considered in this study. The main laboratory strain used in this work was the ISPRA strain, established in 1968 at the European Community Joint Research Centre (Ispra, Italy), with wild flies from Sicily and Greece. The strain has been maintained in the quarantine facility at the Dept. of Animal Biology, University of Pavia (Italy) since 1979. Other laboratory strains used, maintained in the same insectary, were: Egypt II, obtained from the FAO/IAEA Agriculture and Biotechnology Laboratory (Seibersdorf, Vienna, Austria); Guatemala, obtained in 1990 from MOSCAMED (Guatemala); Israel and Hawaii established from flies collected in the corresponding countries, maintained in our laboratory since 1990 and 1992, respectively. In addition, a sample of wild flies collected in May 2007 from coffee drupes in Ruiru near Nairobi, Kenya, part of the ancestral home range of the species, was used in the study.
Egg and embryo collection
To obtain fertilized and unfertilized eggs, two cages were set up using 200 mated or virgin ISPRA females, respectively. When the flies were five-days old, eggs were collected over a 5 min period. The eggs were maintained at 25°C and 65% humidity for set time intervals. After dechorionation (using 1.5-2% hypochlorite solution), the eggs were repeatedly washed in distilled water and individually transferred to 1.5 ml microcentrifuge tubes.
Nucleic acid preparations
DNA and RNA were isolated from individual embryos using TRIzol® Reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol. The extracted RNA was resuspended in RNAse-free water. Genomic DNA was individually extracted from 3rd instar larvae, 7-day old pupae and 3-day old adult flies using the phenol-chloroform extraction method . Following treatment with RNAse A, the DNA was extracted with phenol/chloroform, precipitated with ethanol and resuspended in TE (10 mM Tris-HCl, pH 8, 1 mM EDTA). The DNA concentration was quantified using a Nanodrop ND-1000 spectrophotometer (Nanodrop Technologies Inc., Wilmington, DE, USA).
PCR amplifications and agarose gel electrophoresis
Oligonucleotide primers CcYf (5'-gctcgaagacatgcattgaa-3') and CcYr (5'-gacggtaagtgccattcgtt-3') were designed on a known Y-specific sequence of C. capitata [GenBank:AF115330] using Primer3 . PCR amplifications were performed in 10 μl reaction volumes using ~50 ng DNA, 1.5 mM MgCl2, Reaction Buffer (10 mM Tris, 50 mM KCl; pH 8.3), 0.2 mM dNTPs mixture, 10 pmol of each primer and 1 unit Taq DNA polymerase (Invitrogen). Amplification was achieved on an Eppendorf Mastercycler Gradient using the following programme: an initial denaturing step at 94°C for 2 min; 30 cycles of denaturation at 94°C for 30 sec, annealing at 55.8°C for 1 min, extension at 72°C for 45 sec; and a final extension at 72°C for 5 min. The PCR products were analysed on 14 × 10 cm 1.5% agarose gel slabs in 1 × TAE buffer together with a 1 kb and a 100 bp ladder standard (Invitrogen). The bands were visualized by ethidium bromide staining and exposure to UV light. Bands selected for further analysis were eluted from the gel using the PureLink Quick Gel Extraction Kit (Invitrogen).
EST bioinformatic analysis and primer design
ESTs  were subjected to bioinformatic analysis via BLASTX against the D. melanogaster genome and sequences similar to genes involved in sex determination and cellular blastoderm formation were selected (Table 1). Data on the D. melanogaster gene expression profiles were obtained from the flybase server . The primers used to amplify gene-specific sequences by reverse transcription PCR (RT-PCR) were designed using Primer3 and may span one or more exon sequences (Figure 1).
Medfly homologues of genes involved in sex-determination and cellularization and their function and embryonic expression in D. melanogaster
Gene, (symbol), GenBank Acc.
Drosophilagene, FlyBase ID
Alignment length (aa)
D. melanogastergene functions
D. melanogastergene expression during early embryogenesis
zinc-finger transcription factor involved in MTZ process and neuronal development
transcription factor forming the primary signal of the sex-determination cascade (XSEs) and involved in endoderm development
expressed in female zygote at the end of nuclear cycle 8 (1:12 h)
transcription factor forming the primary signal of the sex-determination cascade (autosomal) and involved in dendrite morphogenesis
expressed in the zygote starting from nuclear cycle 12 (2 h)
corepressor of transcription forming the primary signal of the sex-determination cascade (autosomal) and involved in dendrite morphogenesis
maternal origin, gro transcripts are detectable until stage 11 after complete gastrulation (5:20 h)
Sex lethal (Sxl)
RNA-binding protein, member of the sex-determination molecular cascade
expressed in the zygote starting from nuclear cycle 12-13 (2 h)
RNA-binding protein, member of the sex-determination molecular cascade
maternal origin; an increased amount of tra transcripts is detectable 2-3 h after oviposition
RNA-binding protein, member of the sex-determination molecular cascade
maternal origin. No change in tra2 transcript amount is detectable after 2-3 h and in chromosome 2R deleted embryos
zinc-finger transcription factor involved in sex differentiation process
dsx transcripts are first detectable during stage 11 (9:20 h) in the fat bodies and gonad primordia
sans fille (snf)
U1 snRNA binding protein involved in sex-determination and oogenesis
the SNF protein localizes to the nucleus in all tissues during development
female lethal d (fl(2)d
RNA splicing factor involved in sex-determination and female-germline sex-determination
transcription factor binding protein involved in sex-differentiation process
maternal origin; an increased amount of ix transcripts is detectable 2-3 h after oviposition
protein kinase coupled with GTPase activity involved in cellular-blastoderm formation
spaghetti squash (sqh)
myosin heavy chain binding protein coupled with ATPase activity involved in cellularization and ovarian follicle cell development
slow as molasses (slam)
protein of unknown function involved in cellularization
slam expression rapidly increases during cycle 11 (1:45 h) and peaks during cycle 14 (2:30 h)
Similar to serendipity-α gene product, involved in cellularization
amale specific isoform; bfemale-specific isoform
cDNA preparation and RT-PCR
cDNA was generated using the cloned AMV first-strand synthesis kit (Invitrogen) according to the manufacturer's instructions. The RT-PCR mixture and the cycling conditions were as reported above, but the annealing temperature was 59°C. Established primers and cycling conditions were used to amplify the Cctra sequences . For each set of RT-PCR amplifications control reactions were performed excluding reverse transcriptase. RT-PCR amplifications were considered only when these negative controls produced no amplification products. Electrophoresis and visualization of amplification products were performed as described above.
Cloning and Sequencing of PCR and RT-PCR Products
Gel eluted PCR and RT-PCR products were cloned into the PCR®2.1-TOPO® vector using the TOPO TA cloning kit (Invitrogen). Positive colonies were selected and the size of the insert quantified by Eco RI digestion and electrophoresis. Clones were sequenced using an ABI-310 automatic sequencer and the ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction Kit v. 3.1 (Applied Biosystems, Foster City, CA, USA).
Nucleotide sequences were aligned using T-coffee . Sequence comparisons were performed using the BLAST family of programs from the National Centre for Biotechnology Information . Identity values were calculated using PAUP 4.0b10 .
Southern hybridization analyses
Genomic DNA from individual flies (4 μg) was digested with SspI endonuclease. The digested DNAs were electrophoresed on a 20 × 14 cm 1% agarose gel in 1 × TBE and transferred to a positively charged nylon membrane according to Southern . The membrane was hybridized at 55°C with 200 ng of probe DNA labelled with the Gene Images Alkphos Direct labelling system (GE Healthcare, Little Chalfont, UK) using the random primer method. The hybridization and detection protocols were those described by the manufacturer. Signal detection was performed using CDP-star followed by exposure to autoradiographic film (X-OMAT AR, Kodak).
Chromosome preparation and in situ hybridization
Mitotic chromosome spreads were obtained using medfly third instar larvae . Briefly, brain tissue was incubated in 1% sodium citrate for 10 min at room temperature and transferred to methanol-acetic acid 3:1 solution for 4 min. The material was disrupted in 100 μl 60% acetic acid and dropped onto clean slides and dried. Pre-hybridization was performed according to . In situ hybridization was performed using the following protocol: the probe DNA was labelled using the Biotin High Prime kit (Roche, Basel, Switzerland) and detection of hybridization signals was performed using the Vectastain ABC elite kit (Vector Laboratories, Burlingame, CA, USA) and DAB (3,3'-Diaminobenzidine tetrahydrochloride, Sigma-Aldrich, St. Louis, MO, USA) .
Results and Discussion
A putative MITE insertion on the Y chromosome permits a fast sex-assignment
The PCR-based sexing method reported here was designed to easily assign sexual identity to the earlier medfly life stages in order to study sex-specific expression and/or mRNA splicing of early, zygotically transcribed genes. This method follows previous sexing assays [16–18] that used primers to amplify a Y-specific repetitive sequence (Y114). In the previous assays, amplification was limited to the male: the presence of a sex-specific pattern only in males may give a false sex assignment should the amplification fail.
The CcYf and CcYr primers were designed outside of the AT-rich element of the Y-specific sequence to amplify a 250 bp sequence (Figure 2A). Genomic DNA from individual male and female adults from the ISPRA laboratory strain was subjected to PCR amplification and, unexpectedly, amplification products were obtained in both sexes. The amplification pattern was different in males and females: only one distinct band of 242 bp was amplified in females, whereas two distinct bands of 250 bp and 727 bp were amplified from the male samples (Figure 2B). These male- and female-specific amplification patterns remain constant in all medfly laboratory strains analyzed (Egypt II, Guatemala, Israel and Hawaii) and from wild flies collected in Kenya, the home range of C. capitata , indicating that there is no pattern variability, irrespective of the geographic origin of the samples. Moreover they are constant in all developmental stages, including early embryos (Figure 2C).
The 242, 250 and 727 bp amplicons were sequenced [Genbank: GU122238-Genbank: GU122240] and BLASTN analyses showed that they share high identity with the Y-specific repetitive DNA sequence (ranging from 93% to 95%). Alignments among the three sequences indicated a high degree of identity (89.6 to 94.5%) (Figure 3). The male-specific 727 bp sequence, used as probe in Southern blot analysis on male and female genomic DNA, resulted in multiple intense signals with male-derived DNA and at least one weaker signal with female-derived DNA samples (Figure 4A). In situ hybridization on mitotic chromosomes again using the 727 bp sequence as a probe clearly indicated that this is a repetitive sequence located on the long arm of the Y chromosome (Figure 4B).
The female 242 bp sequence lacks eight base pairs in the 3' region, that are present in both male sequences, while the male 727 bp sequence contains an additional 477-bp sequence, which is absent from the other two fragments. This indel has no homology to the Y-specific repetitive DNA sequence, but it shares 86-89% identity with two medfly ESTs , the functional/biological roles of which are unknown. The presence of terminal inverted repeats in these sequences, their small size, the absence of coding sequence, the potential to form a stable secondary structure, the presence of duplicated target regions in the ESTs and the presence of multiple copies in the medfly genome, as confirmed by Southern hybridization (Additional file 1, Figure S1; Additional file 2, Figure S2), led us to hypothesize that the 467 bp indel sequence is a miniature inverted repeat transposable element (MITE) . It is the insertion of the putative MITE into one or more repetitive units on the Y-chromosome that results in the amplification of the male-specific 727 bp sequence (Figure 2A).
The dynamics of the sex determination molecular cascade in early embryos
The MTZ transition process integrates the post-transcriptional regulation mechanisms of the maternally inherited transcripts and the transcriptional activation and control of the zygotic genome. Particular attention has to be focused on the events that regulate the sex determination process, because the zygote has to reset the information of "female development" inherited from its mother to its own sex-development information. This process involves the activities of transcriptional and of mRNA-splicing regulators.
Three of the medfly genes involved in the molecular sex-determination cascade have previously been identified, Cctra [22, 31], Cctra2 [7, 32], and Ccdsx . Cctra is the key-gene of the sex-determination cascade: it generates mRNAs encoding full-length active proteins only in females and displays an autocatalytic function, which guarantees the female-specific development of cell memory. Cctra, in cooperation with Cctra2, determines the sex-specific splicing of Ccdsx, the transcription factor that is the regulator of the sex-differentiation process.
Cctra and Ccdsx mRNAs are maternally inherited as female-specific splicing variants; the maternal information for female-specific development is reset in embryos through the reprogramming of Cctra mRNA splicing and the degradation of the maternal Ccdsx mRNAs (Figure 5A).
Cctra expression in the zygote starts between 4 and 5 h after oviposition. In 5 to 8 h embryos of both sexes, the Cctra mRNA population is heterogeneous, composed of transcripts of different molecular weights (Figure 5A). These patterns are not attributable to contaminating genomic DNA as controls showed that there are no amplification products in the absence of retro-transcriptase. Notably, in 6 h embryos the Cctra amplification patterns appear to be almost identical in both sexes. Sequencing of these transcripts showed that both male and female embryos contain six transcripts including the mature male and female Cctra mRNA and partially or unprocessed Cctra transcripts [GenBank: GU122232-Genbank: GU122237]. Moreover, the splicing mechanism does not appear to act in an orderly manner as transcripts with different combinations of intron/exon sequences are detected (Figure 6).
The presence of different splicing transcripts, including male-specific transcripts, in female embryos indicate that the maternal Cctra transcripts do not produce enough active CcTRA protein to prevent the male mode of splicing in female embryos. Concentration-dependant splicing activity is a general property of proteins belonging to the SR-superfamily of splicing factors  which includes the products of Cctra and Cctra2. Consequently, we suggest that in female embryos, the Cctra autoregulation loop produces more and more active CcTRA female protein, reaching the threshold concentration that results in the female-specific splicing of all the Cctra transcripts at 8-9 h (Figure 5A). By contrast, in male embryos the postulated Y-linked male-determining factor, M [15, 22, 32], directly or indirectly blocks the Cctra autoregulation loop, inducing the translation of more and more inactive CcTRA protein, resulting in the male-specific splicing of all the Cctra transcripts.
Apart from these sex-specific differences in splicing, as we cannot detect clearcut differences in Cctra transcription levels between male and female embryos, it is unlikely that the M factor interferes with Cctra transcription. Even if we cannot completely exclude this hypothesis, we presume that the M factor acts by affecting Cctra mRNA splicing or CcTRA protein activity (e.g. reducing its phosphorylation).
Sequencing of the Cctra transcripts from embryos also permitted the characterization of the genomic organization of the gene: the exon/intron structure is different from that reported previously . In particular, we found that intron 1 is 50 bp longer than expected and we found no evidence of the presence of a third 52 bp long intron within the genomic portion that we analyzed. This discrepancy could be due to the general genetic variability that occurs between populations of different geographical origins and between different laboratory strains.
In contrast to Cctra, the Cctra2 transcripts appeared in RT-PCR experiments as a single band in both sexes at all the embryonic stages analysed, even though the primers spanned five introns (Figure 1). The presence of only the mature transcript, as ascertained by comparison with the published genomic sequence [Genbank:EU437408], could be due to Cctra2 not being expressed in the zygote during early embryogenesis and that only maternal Cctra2 mRNA supplies the required CcTRA2 protein. This hypothesis is supported by the observation in D. melanogaster that the amount of tra2 transcripts is lower in 2-3 h than in 0-1 h embryos, suggesting that there is no zygotic expression or that the zygotic expression cannot compensate for the concomitant degradation of the maternal mRNAs. Furthermore, in embryos in which the tra2-carrying chromosome 2 is deleted, the quantity of tra2 mRNAs is similar to that in normal embryos, indicating that the maternal tra2 mRNA is stable and that there is no zygotic contribution .
As far as it concerns the third medfly sex determination gene, Ccdsx expression begins in 10 h male and female embryos and its transcripts are sex-specifically spliced, as by this stage only female embryos express the functional female-specific CcTRA protein (Figure 5A).
In D. melanogaster sex determination occurs before cellularization as it is coupled with dosage compensation . Our data show that sex determination in the medfly also occurs before cellularization. This is surprising as the medfly X chromosome is largely heterochromatic and contains very few genes, suggesting that there is no need for dosage compensation . In both species, the establishment of this cascade during the first stages of embryogenesis requires the activity of proteins that act in a concentration-dependent manner. These molecular events could be facilitated in a syncitial stage rather than in a cellularized organism.
Our results show that the medfly sex determination cascade is complete at 10 h, before the end of cellularization (9-11 h ), whereas in D. melanogaster it ends after gastrulation (stage 13) . A possible explanation for this phenomenon may be the shorter cascade in C. capitata and also the late onset of cellularization (9-11 h) compared to 2.10-2.45 h in D. melanogaster.
Possible medfly sex determination-related genes
In D. melanogaster the ratio of X-linked and autosomal gene products represents the primary signal of the sex determination cascade. Among these genes, sisterless A (sisA) is one of the X-linked genes and deadpan (dpn) and groucho (gro) are autosomal genes. Early during development, the ratio of their product concentrations, together with those of other genes, determines the expression status of the key gene Sex-lethal (Sxl) starting from an early promoter (SxlPe) [37, 38]. In females the Sxl early mRNA encodes an active SXL protein. At later developmental stages a maintenance promoter (SxlPm) triggers Sxl transcription in both males and females, but this pre-mRNA requires the presence of an active SXL protein before it can mature into full-length SXL-encoding mRNA. Only in female embryos does the SXL translated from SxlPe mRNA act as the memory device for female sexual development via its auto-regulatory function. The function of SXL requires the activity of cofactors, such as sans-fille (snf) and female lethal (d) (fl(2)d). When Sxl is activated, it sets in motion a cascade of regulatory genes, transformer (tra), transformer-2 (tra2) and doublesex (dsx) . In female embryos the DSXF protein cooperates with the product of the intersex (ix) gene to drive the expression of the genes responsible for sexual differentiation .
Even though the C. capitata sex determination cascade differs from that of D. melanogaster, several medfly EST transcripts with homology to the sisA, dpn and gro genes were identified (Table 1). Their expression in the embryos comply with the classification of their D. melanogaster homologues: CcsisA and Ccdpn are expressed by the zygote earlier than 4 h after oviposition while Ccgro is a maternal gene, as its transcripts are inherited from the mother (Figure 5B). There were no apparent differences in expression between male and female embryos, apart from an apparent lower amount of Ccdpn transcripts in the 8-9 h female embryos, compared to their male counterparts. This difference has to be confirmed with further analyses and its biological significance has to be investigated. As the medfly primary sex determination signal is not based on the ratio of X and autosomes gene products, the presence of transcripts of these genes during embryonic development opens the interesting prospect of identifying the roles that they play during early medfly development.
As expected, the CcSxl gene has no sex-specific splicing variants and it is inherited maternally as a mature transcript. This difference with respect to D. melanogaster was also reported for MdSxl in Musca domestica , a species that appears to have a sex determination system similar to that of the medfly . This is indicative of a fundamental difference in transcriptional regulation of this gene and shows that this gene has a different function in drosophilid and non-drosophilid species.
The medfly embryos also inherit transcripts of the Ccsnf, Ccfl(2)d and Ccix genes. These genes could have the same functions in the medfly as their Drosophila homologues, as co-regulators of Cctra splicing (Ccfl(2)d and Ccsnf) or as a co-factor of Ccdsx (Ccix).
Like Cctra, the Ccgro, CcSxl and Ccfl(2)d maternal genes present at least two PCR products in 5 to 8 h old embryos in both sexes (Figure 5B). Sequencing demonstrated that the larger fragment corresponds to the pre-mRNA, while the smaller fragment is the spliced/mature mRNA. This suggests that the newly transcribed mRNAs are not immediately processed into mature forms.
To conclude, we can hypothesize that two waves of zygotic expression occur in the medfly: the earliest genes, which include CcsisA and Ccdpn, are expressed before 4 h after oviposition. The second major burst of transcription occurs at 5 h after oviposition and the Ccgro, CcSxl, Cctra and Ccfl(2)d genes belong to this second wave.
Towards the cellular blastoderm formation
A second major event during the first stage of embryogenesis is the cellular blastoderm formation, which in C. capitata occurs between 9 h and 11 h after oviposition .
In D. melanogaster a set of genes is involved in this process that includes maternal, maternal-zygotic and zygotic genes. In the medfly EST database we identified two sequences that are similar to the D. melanogaster maternal genes, Rho1 and spaghetti-squash (sqh). These genes are also maternal in the medfly and the immature form of CcRho1 is present in 5-8 h old embryos, indicating that this gene is expressed during the second major wave of transcription in the embryo (Figure 5C).
The slow as molasses (slam) and serendipity-α (sry-α) genes are expressed only in the zygote in D. melanogaster. The Ccslam gene complies with this classification, whereas the Ccsry-α gene is maternally inherited in the embryos as a mature transcript. Again, in 5 to 8 h old embryos, the immature form of the Ccsry-α mRNA was detected. The recent embryonic in situ detection of Ccsry-α did not reveal the presence of maternal transcripts of this gene . This discrepancy may be due to the higher sensibility of a PCR- compared to a hybridization-based approach. What is more, the Ccsry-α gene showed greater amino-acid similarity to the D. melanogaster CG8247 gene than to sry-α. CG8247 is a sry-α-like gene that is also involved in cellular blastoderm formation but is, unlike sry-α, inherited maternally in D. melanogaster
The putative CcZelda gene
In D. melanogaster the presence of TAGteam sites in promoter-sequences is a feature of genes that are expressed early during embryogenesis ; these sites are bound by the zinc-finger transcription factor Zelda . zelda (zld) transcripts are present in the germ line cells of the ovary, in unfertilized eggs and throughout early embryonic development, while later on in development zld becomes restricted to the nervous system. The Zld protein is a general activator of early genes and it is responsible for the MTZ transition. It regulates the activation of genes involved in the three early processes that occur in the embryos, sex determination (including sisA, Sxl and dpn), cellular blastoderm formation (including slam and sry-α), and the dorsoventral genes.
In the medfly EST database we identified a sequence that encodes for a protein with a high degree of identity with Zelda (Table 1). Like its putative homologue zld, Cczld transcripts are maternally derived and they are present throughout the syncitial period (Figure 5D). The promoter regions of two of the medfly genes putatively regulated by Zelda, Ccslam and Ccsry-α, are known, but it was not possible to identify conserved TAGteam sites between the medfly and Drosophila homologues of these genes. In the D. melanogaster slam gene the TAG sites are present only in the first intron, as we do not have the genomic sequence of Ccslam we cannot confirm that these sites are conserved in the medfly gene. In the sry-α gene the TAG sites are located in the promoter, whereas in Ccsry-α the promoter lacks these sites. This lack of TAGteam sequences in Ccsry-α could be a further demonstration that this gene has a different transcription-control molecular mechanism with respect to sry-α in Drosophila. The sequencing of the entire medfly genome, currently in progress, will provide a powerful tool to understand the molecular mechanisms of the MTZ transition process in this species.
The novel PCR-based sexing method reported in this paper takes advantage of a putative MITE insertion on the medfly Y chromosome and it allows the analysis of the MTZ transition at very early stages of embryogenesis in both sexes. During this period, as in D. melanogaster, two main events take place: sex determination and the cellular blastoderm formation. The zygotic transcriptional activation of genes involved in these processes follows two waves (Figure 7). The first wave starts before 4 h after oviposition and includes the zygotic genes CcsisA, Ccdpn and Ccslam. The second major burst of expression activation begins 5 h after oviposition and includes the maternal-zygotic genes Ccgro, CcSxl, Cctra, Ccfl(2)d, CcRho1 and Ccsry-α. The C. capitata sex determination cascade is shorter than that of D. melanogaster, probably due to the fact that Cctra, instead of CcSxl, is the master switch gene. Further differences in C. capitata, compared to Drosophila early development, concern the late onset of cellularization and the accomplishment of sex determination before the cellular blastoderm formation. This last issue proves that the transcription of Ccdsx begins during the cellularization process, while that of D. melanogaster occurs after the cellular blastoderm formation.
The maternal information of female development is reset in the medfly embryos through the reprogramming of Cctra mRNA splicing and the degradation of the maternal Ccdsx mRNA. The heterogeneity of the Cctra mRNA population during the "splicing-resetting" phase (5-8 h after oviposition) is indicative of a threshold-dependent activity of the CcTRA protein. Thus we suggest that the maternally-inherited Cctra transcripts in female embryos are not sufficient to produce enough active protein to inhibit the male mode of Cctra splicing. These data suggest that the M-factor  acts by inhibiting the CcTRA protein activity, so that it does not reach the threshold concentration required for female-specific splicing of all the Cctra transcripts. The slow rate of development and the apparent inefficacy of the splicing mechanism in the pre-cellular blastoderm, exemplified by the presence of immature mRNA precursors of many genes (Ccgro, CcSxl, Ccfl(2)d, CcRho1 and Ccsry-α), could facilitate the M-factor activity. This model of sex determination in the medfly differs from that previously proposed .
Operationally, the slow rate of early development in the medfly may facilitate the study of the expression and splicing mechanisms that occur in slow motion compared to D. melanogaster, thus making C. capitata a suitable candidate to deepen our knowledge of the molecular mechanisms involved in the MTZ transition in Diptera.
As a practical consequence, promoter and enhancer sequences that are active in early stages of development will become available as tools for population control strategies against pest insects.
We thank the two anonymous referees for their constructive comments and suggestions. We thank Saverio Brogna for critically reading the manuscript. This work was supported by the grant PRIN 2006, #2006050870_004 (Programmi di Ricerca Scientifica di Rilevante Interesse Nazionale) from the Italian Ministry of University and Research.
Department of Animal Biology, University of Pavia
Entomology Unit, FAO/IAEA Agriculture and Biotechnology Laboratory, Joint FAO/IAEA Programme, International Atomic Energy Agency
Division of Genetics, Cell and Developmental Biology, Department of Biology, University of Patras
Gouw JW, Pinkse MW, Vos HR, Moshkin Y, Verrijzer CP, Heck AJ, Krijgsveld J: In vivo stable isotope labeling of fruit flies reveals post-transcriptional regulation in the maternal-to-zygotic transition. Mol Cell Proteomics. 2009, 8 (7): 1566-1578. 10.1074/mcp.M900114-MCP200.PubMed CentralView ArticlePubMed
Liang HL, Nien CY, Liu HY, Metzstein MM, Kirov N, Rushlow C: The zinc-finger protein Zelda is a key activator of the early zygotic genome in Drosophila. Nature. 2008, 456 (7220): 400-403. 10.1038/nature07388.PubMed CentralView ArticlePubMed
Tadros W, Goldman AL, Babak T, Menzies F, Vardy L, Orr-Weaver T, Hughes TR, Westwood JT, Smibert CA, Lipshitz HD: SMAUG is a major regulator of maternal mRNA destabilization in Drosophila and its translation is activated by the PAN GU kinase. Dev Cell. 2007, 12 (1): 143-155. 10.1016/j.devcel.2006.10.005.View ArticlePubMed
Benoit B, He CH, Zhang F, Votruba SM, Tadros W, Westwood JT, Smibert CA, Lipshitz HD, Theurkauf WE: An essential role for the RNA-binding protein Smaug during the Drosophila maternal-to-zygotic transition. Development. 2009, 136 (6): 923-932. 10.1242/dev.031815.PubMed CentralView ArticlePubMed
Bushati N, Stark A, Brennecke J, Cohen SM: Temporal reciprocity of miRNAs and their targets during the maternal-to-zygotic transition in Drosophila. Curr Biol. 2008, 18 (7): 501-506. 10.1016/j.cub.2008.02.081.View ArticlePubMed
De Renzis S, Elemento O, Tavazoie S, Wieschaus EF: Unmasking activation of the zygotic genome using chromosomal deletions in the Drosophila embryo. PLoS Biol. 2007, 5 (5): e117-10.1371/journal.pbio.0050117.PubMed CentralView ArticlePubMed
Gomulski LM, Dimopoulos G, Xi Z, Soares MB, Bonaldo MF, Malacrida AR, Gasperi G: Gene discovery in an invasive tephritid model pest species, the Mediterranean fruit fly, Ceratitis capitata. BMC Genomics. 2008, 9: 243-10.1186/1471-2164-9-243.PubMed CentralView ArticlePubMed
Beverley SM, Wilson AC: Molecular evolution in Drosophila and the higher Diptera II. A time scale for fly evolution. J Mol Evol. 1984, 21 (1): 1-13. 10.1007/BF02100622.View ArticlePubMed
Kwiatowski J, Skarecky D, Bailey K, Ayala FJ: Phylogeny of Drosophila and related genera inferred from the nucleotide sequence of the Cu, Zn Sod gene. J Mol Evol. 1994, 38 (5): 443-454. 10.1007/BF00178844.View ArticlePubMed
Peterson BK, Hare EE, Iyer VN, Storage S, Conner L, Papaj DR, Kurashima R, Jang E, Eisen MB: Big genomes facilitate the comparative identification of regulatory elements. PLoS ONE. 2009, 4 (3): e4688-10.1371/journal.pone.0004688.PubMed CentralView ArticlePubMed
Schetelig MF, Horn C, Handler AM, Wimmer EA: Development of an embryonic lethality system for SIT in Ceratitis capitata. Area-wide Control of Insect Pests: from Research to Field Implementation. Edited by: Vreysen MJB, Robinson AS, Hendrichs J. 2007, Dordrecht, Netherlands Springer, 85-93. full_text.View Article
Saccone G, Pane A, Polito LC: Sex determination in flies, fruitflies and butterflies. Genetica. 2002, 116 (1): 15-23. 10.1023/A:1020903523907.View ArticlePubMed
White I, Elson-Harris MM: Fruit flies of economic significance: their identification and bionomics. 1992, Wallingford, Oxon, UK: CAB International in association with ACIAR
Anleitner JE, Haymer DS: Y enriched and Y specific DNA sequences from the genome of the Mediterranean fruit fly, Ceratitis capitata. Chromosoma. 1992, 101 (5-6): 271-278. 10.1007/BF00346005.View ArticlePubMed
Willhoeft U, Franz G: Identification of the sex-determining region of the Ceratitis capitata Y chromosome by deletion mapping. Genetics. 1996, 144 (2): 737-745.PubMed CentralPubMed
San Andrés V, Urbaneja A, Sabater-Muñoz B, Castañera P: A novel molecular approach to assess mating success of sterile Ceratitis capitata (Diptera: Tephritidae) males in sterile insect technique programs. J Econ Entomol. 2007, 100 (4): 1444-1449. 10.1603/0022-0493(2007)100[1444:ANMATA]2.0.CO;2.View ArticlePubMed
Zhou Q, Untalan PM, Haymer DS: Repetitive A-T rich DNA sequences from the Y chromosome of the Mediterranean fruit fly, Ceratitis capitata. Genome. 2000, 43 (3): 434-438. 10.1139/gen-43-3-434.View ArticlePubMed
Douglas LJ, Untalan PM, Haymer DS: Molecular sexing in the Mediterranean fruit fly, Ceratitis capitata. Insect Biochem Mol Biol. 2004, 34 (2): 159-165. 10.1016/j.ibmb.2003.08.006.View ArticlePubMed
Baruffi L, Damiani G, Guglielmino CR, Bandi C, Malacrida AR, Gasperi G: Polymorphism within and between populations of Ceratitis capitata: comparison between RAPD and multilocus enzyme electrophoresis data. Heredity. 1995, 74 (Pt 4): 425-437. 10.1038/hdy.1995.60.View ArticlePubMed
Rozen S, Skaletsky H: Primer3 on the www for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMed
Tweedie S, Ashburner M, Falls K, Leyland P, McQuilton P, Marygold S, Millburn G, Osumi-Sutherland D, Schroeder A, Seal R, et al: FlyBase: enhancing Drosophila Gene Ontology annotations. Nucleic Acids Res. 2009, D555-559. 10.1093/nar/gkn788. 37 Database
Pane A, Salvemini M, Delli Bovi P, Polito C, Saccone G: The transformer gene in Ceratitis capitata provides a genetic basis for selecting and remembering the sexual fate. Development. 2002, 129 (15): 3715-3725.PubMed
Poirot O, O'Toole E, Notredame C: Tcoffee@igs: A web server for computing, evaluating and combining multiple sequence alignments. Nucleic Acids Res. 2003, 31 (13): 3503-3506. 10.1093/nar/gkg522.PubMed CentralView ArticlePubMed
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (and Other Methods). 1998, Sinauer Associates, Sunderland, MA
Southern EM: Detection of specific sequences among DNA fragments separated by gel electrophoresis. J Mol Biol. 1975, 98 (3): 503-517. 10.1016/S0022-2836(75)80083-0.View ArticlePubMed
Zacharopoulou A: Polytene chromosome maps in the Medfly Ceratitis capitata. Genome. 1990, 33 (2): 184-197.
Zacharopoulou A, Frisardi M, Savakis C, Robinson AS, Tolias P, Konsolaki M, Komitopoulou K, Kafatos FC: The genome of the Mediterranean fruitfly Ceratitis capitata: localization of molecular markers by in situ hybridization to salivary gland polytene chromosomes. Chromosoma. 1992, 101 (7): 448-455. 10.1007/BF00582839.View ArticlePubMed
Malacrida AR, Gomulski LM, Bonizzoni M, Bertin S, Gasperi G, Guglielmino CR: Globalization and fruitfly invasion and expansion: the medfly paradigm. Genetica. 2007, 131 (1): 1-9. 10.1007/s10709-006-9117-2.View ArticlePubMed
Wessler SR, Bureau TE, White SE: LTR-retrotransposons and MITEs: important players in the evolution of plant genomes. Curr Opin Genet Dev. 1995, 5 (6): 814-821. 10.1016/0959-437X(95)80016-X.View ArticlePubMed
Ruiz MF, Milano A, Salvemini M, Eirín-López JM, Perondini AL, Selivon D, Polito C, Saccone G, Sánchez L: The gene transformer of anastrepha fruit flies (Diptera, tephritidae) and its evolution in insects. PLoS ONE. 2007, 2 (11): e1239-10.1371/journal.pone.0001239.PubMed CentralView ArticlePubMed
Salvemini M, Robertson M, Aronson B, Atkinson P, Polito LC, Saccone G: Ceratitis capitata transformer-2 gene is required to establish and maintain the autoregulation of Cctra, the master gene for female sex determination. Int J Dev Biol. 2009, 53 (1): 109-120. 10.1387/ijdb.082681ms.View ArticlePubMed
Saccone G, Salvemini M, Pane A, Polito LC: Masculinization of XX Drosophila transgenic flies expressing the Ceratitis capitata DoublesexM isoform. Int J Dev Biol. 2008, 52 (8): 1051-1057. 10.1387/ijdb.082657gs.View ArticlePubMed
Fu XD: The superfamily of arginine/serine-rich splicing factors. RNA. 1995, 1 (7): 663-680.PubMed CentralPubMed
Graham P, Penn JK, Schedl P: Masters change, slaves remain. Bioessays. 2003, 25 (1): 1-4. 10.1002/bies.10207.View ArticlePubMed
Hempel LU, Oliver B: Sex-specific DoublesexM expression in subsets of Drosophila somatic gonad cells. BMC Dev Biol. 2007, 7: 113-10.1186/1471-213X-7-113.PubMed CentralView ArticlePubMed
Harrison DA: Sex determination: controlling the master. Curr Biol. 2007, 17 (9): R328-330. 10.1016/j.cub.2007.03.012.View ArticlePubMed
Lu H, Kozhina E, Mahadevaraju S, Yang D, Avila FW, Erickson JW: Maternal Groucho and bHLH repressors amplify the dose-sensitive X chromosome signal in Drosophila sex determination. Dev Biol. 2008, 323 (2): 248-260. 10.1016/j.ydbio.2008.08.012.PubMed CentralView ArticlePubMed
Schütt C, Nöthiger R: Structure, function and evolution of sex-determining systems in Dipteran insects. Development. 2000, 127 (4): 667-677.PubMed
Chase BA, Baker BS: A genetic analysis of intersex, a gene regulating sexual differentiation in Drosophila melanogaster females. Genetics. 1995, 139 (4): 1649-1661.PubMed CentralPubMed
Meise M, Hilfiker-Kleiner D, Dübendorfer A, Brunner C, Nöthiger R, Bopp D: Sex-lethal, the master sex-determining gene in Drosophila, is not sex-specifically regulated in Musca domestica. Development. 1998, 125 (8): 1487-1494.PubMed
Burghardt G, Hediger M, Siegenthaler C, Moser M, Dübendorfer A, Bopp D: The transformer2 gene in Musca domestica is required for selecting and maintaining the female pathway of development. Dev Genes Evol. 2005, 215 (4): 165-176. 10.1007/s00427-004-0464-7.View ArticlePubMed
Schetelig MF, Caceres C, Zacharopoulou A, Franz G, Wimmer EA: Conditional embryonic lethality to improve the sterile insect technique in Ceratitis capitata (Diptera: Tephritidae). BMC Biol. 2009, 7: 4-10.1186/1741-7007-7-4.PubMed CentralView ArticlePubMed
ten Bosch JR, Benavides JA, Cline TW: The TAGteam DNA motif controls the timing of Drosophila pre-blastoderm transcription. Development. 2006, 133 (10): 1967-1977. 10.1242/dev.02373.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.