The study of the transformer gene from Bactrocera dorsalis and B. correcta with putative core promoter regions
© Laohakieat et al. 2016
Received: 1 October 2015
Accepted: 25 January 2016
Published: 1 February 2016
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).
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.
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.
Eukaryotic sexual reproduction is an important biological process for the continuity of life. Meiosis provides diverse gametes which can be reunited through sexual selection and fertilization. The progeny, endowed with a suitable genome, are naturally selected before the next round of mating. The male and female developmental dimorphism is fundamental for the existence of sexual characteristic .
The sex determination of insects is one of the most well-known system [2, 3]. There are several mechanisms that determine the sexual fate among insect models [3, 4]. The primary signals convey genetic instruction to the sex determination switches in order to designate and sustain sexual identity during development and throughout an organism’s life [3, 4]. In very early embryonic development (approximately before the blastoderm cellularization), an upstream regulator is usually required to be established as a stable genetic switch that serves as a device memory for sex-specific cell fate [5–15].
In Drosophila melanogaster, the primary signal is composed of proteins encoded by a group of genes linked to the X-chromosome, called X-linked signal elements (XSE) [16, 17]. Their encoded proteins interact with the other group of proteins encoded by autosomally localized genes to generate sex-specific dosages (XX/AA or XY/AA) [17, 18]. These dosages regulate the transcription of Sex-lethal (Sxl), upstream regulator, which is also an X-linked gene that presents in both sexes. However, only the XX/AA female dosage can activate the early female-specific promoter of the Sxl gene. This early female-specific SXL protein is a sex-specific splicing regulator of the later promoter Sxl pre-mRNA. The secondary SXL protein can autoregulate the female-specific splicing of its own pre-mRNA [17–20] and the subordinated control gene as transformer (tra) to maintain female development [21, 22]. In contrast, with males, the absence of early SXL protein leads to the production of male-specifically spliced Sxl transcripts which incorporate in-frame stop codons, resulting in the truncated male-specific SXL protein [19, 23].
In Ceratitis capitata (Wiedemann), a tephritid fruit fly, the primary signal is an uncharacterized male-determiner M (M factor) but was only mapped on a Y-chromosome . The M factor was proposed to regulate the tra gene, in C. capitata [6, 7, 25]. It was postulated that the M factor inhibited the maternal TRA protein, which is required for female-specific splicing of tra pre-mRNA and results in a zygotic male-specifically spliced tra mRNA [6, 7, 25]. The tra gene is an upstream regulator of dsx and fru genes in C. capitata because the embryonic RNAi against tra and/or tra-2 experiments led to a stable splicing pattern change in tra, dsx, and fru pre-mRNAs only in XX individuals [7, 25]. In the other tephritid fruit flies, such as the Bactrocera species and Anastrepha species, the tra genes show similar molecular organization and expression of hierarchy patterns to those found in C. capitata [8, 26–29]. These results suggest that the tra orthologues may be the upstream regulators under the M factor control. The alternative sex-specific splicing of the tra gene is under the control of its own autoregulatory loop [6–8, 25–31]. This decision is proposed to be initiated by the signal from the maternal TRA protein, where the M factor is absent in XX females [6–8, 25]. In Drosophila, the functional TRA protein co-operates with the constitutively expressed Transformer-2 (TRA-2) protein, and the complex acts as a splicing regulator [32, 33]. However, in C. capitata, the functional analysis of TRA/TRA-2 complex was indirectly inferred from the RNAi against Cctra-2 effecting tra splicing as well as dsx and fru splicing in XX individuals only . Since the functional analysis of the TRA/TRA-2 complex has been experimentally proven only in D. melanogaster, it was postulated that in C. capitata the TRA/TRA-2 complex may bind to the putative TRA/TRA-2 binding sites located on the male-specific exons and their flanking introns in order to block the strong splice site of male-specific exons. Thereafter, the other SR proteins and spliceosome machinery may be recruited, and the male-specific exons are skipped, resulting in female-specific splicing and the production of a functional TRA protein. On the other hand, in males where the M factor is present, this entire situation may be turned off. Without TRA/TRA-2 regulation, the splicing may be followed by the default which contains in-frame stop codons on male-specific exons and leads to the production of a non-functional truncated TRA protein [6–8, 25–31]. The functional CcTRA proteins may be the sex-specific splicing regulators . They may not only maintain their autoregulatory loop but may regulate the alternative sex-specific splicing of a highly conserved bottom global effector in the sex determination pathway, the so-called doublesex (dsx) gene . The sex determining message carried by this DSX protein can execute sexual differentiation processes in a very broad scheme.
The male- and female-specific DSX (DSXM and DSXF) proteins play transcription factor roles to regulate the downstream genes involved in the sexual differentiation [32, 34]. Most of the mechanism underlining the regulatory roles of the TRA/TRA-2 complex for dsx sex-specifically spliced transcripts was only experimentally proven in D. melanogaster [32–36] while much indirect evidence was gathered from the other tephritid species [37–44]. The presence of functional TRA and TRA-2 proteins acts as a splicing activator that binds to the dsx repeat element (dsx RE) on the female-specific exon. This binding makes a nearby cis-purine rich element (dsx PRE) become a stronger splicing acceptor that facilitates female-specific splicing resulting in the DSXF protein. In contrast, the presence of truncated TRA cancels the strength and potential of this splicing acceptor in males due to no TRA/TRA-2 binding [32–36]. Therefore, the male-specifically spliced dsx transcript contains male-specific exons. This leads to the production of the DSXM protein. The highly conserved dsx genes were also cloned and characterized in many tephritid fruit fly models [37–42] including Bactrocera dorsalis (Hendel) [43, 44] and Bactrocera correcta (Bezzi) . The sex determination pathway would be better understood if tra [7, 8, 26–29] and tra-2 [8, 25, 28–30] counterparts were also characterized. Moreover, it would be advantageous to carry out comparative studies of tra and dsx relationships within the tephritid group to refine the model of sex determination. The role of tra in the sex determination pathway and its related practical applications in pest control programs can be further understood and realized when molecular comparisons are carried out in model organisms belonging to a highly diverse agriculture pest genus such as Bactrocera. The oriental fruit fly (B. dorsalis) and guava fruit fly (B. correcta) belong to the Bactrocera genus, are economically important pests, and are therefore considered targets for control programs .
In this study, we isolated and characterized tra genes from B. dorsalis and B. correcta as well as their putative core promoters. After the full-length cDNA of adult males and females were obtained, RT-PCR was performed with appropriate primers to elucidate the alternative sex-specific splicing in male and female transcripts. The structure and amino acid sequences of both genes were highly conserved among the Bactrocera species, but less so in other tephritid fruit flies. Expression analyses of tra genes (the upstream regulators) were also studied in relation to the dsx genes (the bottom effector) in the wild-type strain (B. dorsalis) and the Salaya1 genetic sexing-mass-reared strain (B. dorsalis). Moreover, the tra RNAi experiment was carried out for the functional analysis of masculinization using the brown and white pupal sexual dimorphism of the Salaya1 strain as a sexual karyotype identifier. The results demonstrated the functions of the tra gene were highly conserved in the sex determination. In addition, the putative core promoter regions of tra genes in B. dorsalis and B. correcta represented the initial discovery of such within the tephritid species. The newly isolated speculative core promoter regions of tra genes from Bactrocera zonata (Saunders), Bactrocera tryoni (Froggatt), Bactrocera carambolae Drew & Hancock (this work); the previously isolated Bactrocera oleae (Gmelin) ; and the dsx putative core promoters from B. dorsalis and B. correcta  were also uncovered. Several molecular features of tra orthologues were compared across Bactrocera species and their related genera to untangle the sex determination mechanism, suggesting further practical genetic tools.
Results and discussion
Molecular organization of tra genes and alternative splicing patterns
In B. correcta, the tra (Bctra) gene was approximately 2 and 2.3 kb long in female- and male-transcripts, respectively [GenBank: KU254116 and KU254117]. There was also conservation of the five common exons and the two male-specific exons (Fig. 1b). However, the novel male-specific exons found in the Bdtra transcripts were not evident. The ORF of the female-specific transcript was 1239 bp and encoded for 413 amino acids. However, the male TRA protein appeared truncated because there were in-frame stop codons located in the male-specific exons.
The male-specific exons of tra genes in tephritid fruit flies located between the common exons 1 and 2 in C. capitata  and Anastrepha species  or 1B and 2A in Bactrocera species [8, 26, 28, 29] contained at least three in-frame stop codons leading to the production of a non-functional TRA polypeptide. Moreover, these male-specific exons and their flanking introns had the putative splicing regulatory elements. These regions of Bdtra and Bctra were analyzed and then compared to the others: Bactrocera species [8, 26, 28, 29], Anastrepha species , C. capitata , and D. melanogaster  (Additional file 1: Figure S1 and Additional file 2: Table S1). Four putative regulatory elements (i.e., TRA/TRA-2, RBP1 binding sites, TRA-2 ISS sequences, and purine-rich elements) were observed [47, 48]. The two male-specific exons found between the two common exons 1B and 2A are also conserved in all Bactrocera species, as in the other tephritids. In this case, many of those putative regulatory elements were found to be clustered [7, 26, 27]. Future functional studies of these TRA/TRA-2 binding sites are required because there is only information regarding dsx and fruitless (fru) sex-specifically spliced genes from D. melanogaster [33, 36, 49] at present.
Sex-specific splicing analysis
5’ upstream regulatory sequences and putative core promoters
The unique finding of this work is the identification of a region that can act as a tra promoter (Bctra). This region may be a core promoter region because it contains many conserved putative core promoter elements. For example, the putative TATA box and putative Inr sequence were 75 and 83.3 % identical, respectively, to their consensus sequences from the Drosophila core promoter database (Fig. 3) [50–52]. In addition, the location of the upstream “T” in the putative TATA box was exactly at the −31 position in relation to the “A” in the putative Inr consensus sequence, concurring with the TATA box location preference with regard to the Inr motif . The other less conserved putative core promoter elements such as the MTE and DPE (60 and 50 % similarity, respectively) were also observed. Whether any of these putative core promoter elements are functional remains to be proven by further promoter analysis experiments.
The 5’ upstream regulatory sequences of tra genes were also amplified from B. tryoni, B. carambolae, and B. zonata using primers Tpro and 1B-1R specific to LOC105232903 and tra CDS, respectively (this work). The upstream sequence of the translation start site (ATG) of Bttra [GenBank: KU254110], Bcatra [GenBank: KU254109], and Bztra [GenBank: KU254108] encompassed 1126, 1271, and 1268 bp, respectively. These regions were also AT rich (67 %), similar to B. dorsalis and B. correcta. The putative core promoter motifs were observed in B. tryoni, B. carambolae, B. zonata, and B. oleae  in similar vicinities. Moreover, similar putative core promoters were also present in the previously characterized dsx genes of B. dorsalis and B. correcta  (Fig. 3). Again, promoter analysis experiments would be required to prove the detailed functions regarding these speculations.
Conservation of tra genes and the newly recognized TRA protein region within tephritid fruit flies
The pair-wise comparison (identity) of tra or dsx ORF (nucleotide) from fruit flies
A newly highly conserved region was discovered among tephritids with 61–83 % similarity (Fig. 5 and Additional file 5: Figure S4). This new region was named “TEP” because it appeared distinctive among tephritids. The TEP region was located approximately 33 amino acids downstream of domain 2 and was encoded by exons 2A and 2B. The TEP region may represent a specific characteristic of the TRA protein in the sex determination pathway of the tephritidae family.
The conservation level of various protein regions (i.e., domains 1 to 4, and the TEP region) within the TRA protein was as restricted as the functional protein domains found in DSX (OD1 and OD2) [44, 56] in Bactrocera. The term “domain” for the DSX protein has a different meaning because the OD domains are well defined protein domains supported by experimental data . This fact implies that protein interactions in the sex determination pathway may be of high fidelity in function. This is especially noted for domain 2 (TRA-CAM) and domain 4 of the TRA protein, which were identical in all tephritids in this study. Although the identification of short regions showing high similarity in all tephritidae TRA proteins is interesting, these putative domains have no similarity to other proteins in metazoans. This fact diminishes the possibility to propose potential functions. In addition, the lack of any structural studies on the TRA proteins and their four putative domains is a major obstacle towards understanding these putative functional domains.
Functional analyses of masculinization in the Salaya1 genetic sexing strain
The microinjected dsRNA containing a 1.1 kb fragment from the 3’ end of Bdtra cDNA can masculinize B. dorsalis (Additional file 6: Figure S5). In this case, the XY male- and XX female-sexual karyotypes are brown- and white-pupae, respectively, because the B. dorsalis Salaya 1 genetic sexing strain was used. The Salaya1 brown-pupae male contains a Y-autosome translocated chromosome which was derived from the wild-type (Phayathai1 strain) during the development of the genetic sexing strain . The Y-autosome translocation chromosome was induced using a low degree of gamma radiation and underwent genetic screening that assumed the chromosomal rearrangement process. The white-pupae pseudomale can be fertile and produce daughter only progenies. However, the abnormal testes (Additional file 7: Figure S6) were also evident, as per Liu et al. 2015 .
For the Bddsx gene, a single band of male-specific product (484 bp) was observed in the wild-type males and the Salaya1 brown-pupae males. Likewise, a female-specific product (677 bp) was in the females. In the RNAi experiment, only the male-specific product was also present in the brown-pupae males. However, both male- and female-specific products were detected in the white-pupae pseudomales (Fig. 6b).
For the functional analysis, RNAi experimentation confirmed that tra and dsx regulated sexual differentiation [7, 25, 26, 28, 59]. The presence of both male and female tra transcripts in the adult stage suggests that the sex determination and development may not be completed, considering the appearance of abnormal testis development and no mating behavior in most of the pseudomales.
The hypothetical Y-linked male determining factor (M factor) is proposed as that which initiates the blocking of tra autoregulatory loops and prohibits female-specific splicing of the zygotic tra pre-mRNA in precellular blastoderm male embryos of C. capitata and Bactrocera jarvisi (Tryon) [6, 8]. During this period, the translation of the inactive TRA protein increases and abolishes the tra autoregulatory loop in males because there is no functional splicing regulation from TRA/TRA-2 complex binding. The abundance of zygotic male-specific, female-specific, and heterogeneous tra mRNA is apparent in both male and female embryos during this primary signalling [6, 28]. Afterward, the male-specific and female-specific tra transcripts are generally and exclusively established in males and females, respectively, presumably because of the stable genetic switch tra.
The coexistence of male- and female-specific tra transcripts during the same period could not be confirmed in the expression analysis of precellularized B. dorsalis embryos because the embryos were not sex sorted . However, the presence of the male-specific transcript similarly appears before cellularization. This suggests that the roles of the non-characterized, but postulated, M factor may be highly conserved in Bactrocera and tephritid fruit flies [6–8, 25–28]. The introgression of the presumably M factor-bearing Y-chromosome from B. jarvisi and B. dorsalis into B. tryoni and B. carambolae, respectively, could interchangeably regulate normal sex determination [8, 60]. However, in this work, the mixed male- and female-specific tra transcripts were observed in all of the brown-pupae adult (XY) males but not in the white-pupae (XX) females of the Salaya1 strain at the adult stage. We hypothesized that the M factor on the Y-chromosome had acquired a leaky mutant stage from the Y-chromosome rearrangement originating from the GSS construction; the M factor cannot completely block the autoregulatory loop in all male karyotype individuals. The leaky M factor function can be rescued by tra RNAi-treatment in the brown-pupae adult males. This results in the disappearance of female-specific tra transcripts. The situation can also be partly imitated in the application of tra RNAi to the white pupae (XX pseudomales). However, a clear difference is that the same coexistence of both male- and female-specifically spliced tra transcripts did not result in the same co-expression of both male- and female-specifically spliced dsx transcripts as in the case of the Salaya1 male (with the M factor) versus the pseudomale (without the M factor). In the pseudomale, the incomplete expression of a functional TRA protein can form a TRA/TRA-2 complex which is effective enough to control the female-specifically spliced dsx. However, in the Salaya1 male that has a leaky M factor, the similarly incomplete expression of a functional TRA protein may also form a TRA/TRA-2 complex, but it seems that the TRA/TRA-2 complex is not effective enough to activate a female-specifically spliced dsx. This fact suggests a possible direct association between the presences of the M factor and TRA-2 gene in the adult stage. In order to provide a justification, we postulate that the M factor may extend its function into adult males by inactivation of not only TRA, but also TRA-2 proteins. A more recent male sex determination model proposed by Salvemini et al. 2009  suggests the possibility of an M factor repressed TRA-2 protein/gene at the embryonic stage. We propose that this mechanism can function during the adult stage. It is proposed that the M factor directly inhibits TRA protein activation in C. capitata [6, 7] and Bactrocera species . If the leaky mutant M factor in the Salaya1 brown-pupae males affects the TRA and TRA-2 protein activities in adult males, the level of TRA/TRA-2 complexes may not reach the threshold to enhance female-specific exonic splicing of dsx pre-mRNA.
Potential application of tra genes for the development of a genetic sexing system in SIT programs
The comparative study of tra in many Bactrocera species suggests a practical application of the tra based-genetic tools. The cross-species utilization of a very highly conserved male-specific exon/intron cassette in the generation of a sex-specific lethality strain is likely among the Bactrocera genera. The TRA/TRA-2 binding sites located in the last male-specific exon of tra genes are identical in B. dorsalis (this work), B. correcta (this work), B. zonata , B. tryoni , B. jarvisi , and B. oleae . Therefore, a single construct for female-specific RIDL (Release of Insects carrying a Dominant Lethal) or fsRIDL may work very well for various Bactrocera pests because it relies on the male-specific exon splicing of the tra gene [62, 63]. The tra RNAi-treatment can be generally used for sexual transformation. However, the derived pseudomales may not be fully masculinized with normal reproductive potential as in this work and elsewhere [7, 25, 26, 28, 59]. Further studies of the sex determination pathway, especially the finding and characterization of the male-determining loci on Y-chromosome and its primary signal pathways would lead to an ideal tool to develop a new masculinization strategy for the genetic sexing strain. The putative core promoters and their upstream regions should contain all necessary promoter/enhancer cis-regulatory elements because these regions are all flanked by the highly conserved LOC 105232903 orthologues. The potential application of the promoter is for it to be used in an early embryonic driver constructions, such as a tetracycline (Tet)-suppressible system [62, 64–66], because the zygotic tra expression is a very early event, before the cellularization stage.
The comparison of Bdtra and Bctra genes to the other Bactrocera and tephritid flies confirms the conserved functional roles of tra in the sex determination pathway. In addition, a new TEP region is proposed because of the conservation among tephritid fruit flies. The generation of new male-specific splicing exons in the Bdtra supports that the initial stage of the alternative splicing process involves multiple exon retention. The evolution of sex-specific splicing is radiating within species from the B. dorsalis model. All of the 5’ regulatory regions of the Bactrocera tra genes contain characteristics of the putative core promoter. The comparative molecular features of Bdtra and Bctra also offer potential applications directly for the B. dorsalis, B. correta, or universally for other Bactrocera genetic sexing systems for SIT.
Fruit fly strains and specimens
The laboratory stocks, B. dorsalis (Phayathai1 strain) and B. correcta (Phayathai2 strain), were used to isolate and characterize Bdtra and Bctra genes, respectively. The B. dorsalis Salaya1 GSS is based on the brown-white pupal color dimorphisms (i.e., males and females that emerged from brown and white pupae, respectively ). They were used for expression and functional analyses of the Bdtra gene. The fruit flies were reared at 25 °C and 13 L: 11D cycles. B. zonata and B. tryoni specimens were from laboratory stocks of the Insect Pest Control Laboratory of the International Atomic Energy Agency, Seiberdorf, Austria. In addition, samples of B. carambolae were from the Jakarta strain.
GenBank Accession Numbers
The GenBank Accession Numbers of tra and dsx genes used in this study were from various tephritid species and D. melanogaster. The tra and dsx orthologues were separately grouped and consolidated in Additional file 8: Table S2 and Additional file 9: Table S3, respectively.
DNA and RNA extractions
Genomic DNA was individually extracted from adult fruit flies according to Baruffi et al. 1995 . Total RNA was isolated from adult male and female fruit flies (2–4 days old) using TRIzol® reagent (Invitrogen, USA) as per the manufacturer’s instructions.
Isolation of Bdtra and Bctra
The primers BoF and BoR designed from Botra [GenBank: AJ715413], were used to firstly isolate the tra orthologues of B. dorsalis and B. correcta, using genomic DNA as templates. All primers used in this work are described in Additional file 10: Table S4. The partial sequences from both species were aligned to derive primers for 3’RACE and nested PCR. Three to five micrograms of total RNA were reverse transcribed to first strand cDNA using the oligo (dT) adapter primer and ImProm-II™ Reverse Transcriptase (Promega, USA), following the manufacturer’s instructions. One-tenth of the initial reverse transcription reaction was used as a template in 3’ RACE using 2B-F and adapter primers. Amplification was performed using a FlexCycler PCR thermal cycler (Analytik Jena, Germany) with the following conditions: 94 °C for 4 min, held at 72 °C for 5 min while Taq polymerase (Vivantis, Malaysia) was added, then one cycle of 60 °C for 2 min and 72 °C for 2 min, followed by five cycles of 94 °C for 1 min, 60 °C for 2 min, 72 °C for 2 min; 28 cycles of 91 °C for 40 s, 58 °C for 2 min, 72 °C for 2 min; and one cycle of final extension at 72 °C for 7 min. The negative controls for each set of RT-PCR experiments were performed excluding reverse transcriptase in order to ensure no genomic DNA contamination.
For the characterization of the 5’ UTR regions and the translational start site of Bdtra, two sets of RT-PCR [(1B-F – 2B-R1) and (5U-BD1 – 2A-2R)] were performed to generate two contiguous fragments to assemble the 5’UTR region and translational start site. Cycling conditions were as follows: 94 °C for 2 min; 29 cycles of 94 °C for 1 min, 60 °C for 30 s, 72 °C for 2 min; and one cycle of final extension at 72 °C for 7 min. PCR was also performed on genomic DNA templates to determine the exon/intron junctions. Three sets of primers (i.e., 5U-BD1 – 2A-2R, 1B-TraF – 2B-R2, and Nested TraF – 3Tra-R) were used to amplify the contiguous fragments with the following cycling conditions: 94 °C for 5 min; 29 cycles of 94 °C for 1 min, 60 °C for 30 s, 72 °C for 1 min; and one cycle of final extension at 72 °C for 7 min. For the same Bctra characterization, the previously described procedures were carried out, except that the primer 5U-BD1 was replaced by 5U-BC1.
Bdtra and Bctra sequences (this work) were inspected to identify the proper restriction enzymes. Subsequently, primers were designed based on the 5’UTR regions. The inverse PCR method was as per Permpoon et al. 2011  except EcoRI restriction digest. Inverse PCR was carried out with 5U-BD1_R – 5U-BD3 primers (for B. dorsalis) and 5U-BC1_R – 5U-BC4 primers (for B. correcta) using the following cycling conditions: 94 °C for 7 min; 29 cycles of 94 °C for 1 min, 58 °C for 30 s, 72 °C for 5 min; and one cycle of final extension at 72 °C for 10 min.
Functional promoter analysis
A 1.4 kb fragment of the 5’ upstream regulatory region including the putative core promoter and translational start codon of Bctra gene was amplified with the BcT-3_NcoI – BcT-4_XbaI primers using genomic DNA as a template. Cycling conditions were as follows: 94 °C for 5 min; 29 cycles of 94 °C for 1 min, 58 °C for 30 s, 72 °C for 2 min; and one cycle of final extension at 72 °C for 7 min. The amplicon was subsequently cloned into the pGEM®-T Easy Vector (Promega, USA). The NcoI-XbaI cut Bctra fragments were ligated in the NcoI-XbaI cut #1215 vector (pSLaf tTA-SV40_af)  to create the #1215-Bctra vector (pSLaf_Bctra–tTA_af). A DsRed restriction fragment was generated from dsRed_XbaIF – dsRed_BamHIR primers using the #1200 vector (pBac[fa PUb-DsRed]) as a template . The Bctra-DsRed construct was created by molecular cloning replacement of the XbaI-BamHI tTA fragment with the XbaI-BamHI DsRed fragment into the #1215-Bctra vector and called #1215-Bctra-DsRed vector (pSLaf_Bctra–DsRed_af).
The transient expression was performed by microinjection of the #1215-Bctra-DsRed constructs (750 ng/µl) into the dechorionated precellular blastoderm embryos from the B. correcta wild-type strain. The microinjected embryos representing developmental times (from 1 to 48 hours after egg laying (h AEL)) were screened every hour by epifluorescence for the expression of the #1215-Bctra-DsRed construct and compared with the uninjected embryos. The screening was performed using the fluorescence stereomicroscope Olympus SZX16 Reflected Fluorescence System with the filters SZX2-FRFP (BP545-580; BA610IF). Photographic documentation was done using an Olympus DP2-BSW version 2.2.
Isolation of additional 5’ upstream putative regulatory region of tra genes from B. zonata, B. tryoni, and B. carambolae
The forward primer Tpro was designed based on putative coding sequences of the LOC105232903 locus [GenBank: XM_011214797.1] from B. dorsalis because the sequence was found to be immediately upstream of the 5’ regulatory region of the Bdtra and Bctra in this work. The reverse primer 1B-1R was designed on the basis of the conserved tra coding sequences approximately 150 bp from the start codon. The genomic DNA of B. zonata, B. tryoni, and B. carambolae were used as templates. The PCR conditions were as per the following: 94 °C for 5 min; 29 cycles of 94 °C for 1 min, 58 °C for 30 s, 72 °C for 2 min; and one cycle of final extension at 72 °C for 7 min. The 5’ and 3’ sequences of the derived specific PCR products were validated for being the corresponding CDS of the LOC105232903 locus and tra orthologues by multiple sequence alignment with the homologous regions in B. dorsalis and B. correcta.
Cloning PCR fragments
The specific amplicon was purified with the QIAquick gel cleanup kit (QIAGEN, Germany). The products were cloned into the pGEM®-T Easy vector (Promega, USA) according to the manufacturer’s instructions. The recombinant plasmids were transformed into DH5α competent cells and screened. Plasmids were extracted as described by Sambrook et al. 1989 . All sequencing reactions were performed on both strands using the sequencing service from Macrogen Inc., Seoul, Korea.
RT-PCR expression analysis
The first stand cDNA of both sexes from each species was separately synthesized. The reverse transcription method and the negative controls were similar to the Bdtra and Bctra isolation procedure. The RT-PCR expression analysis was performed using common and sex-specific primers, represented in Fig. 2a. Cycling conditions were as follows: 94 °C for 2 min; 29 cycles of 94 °C for 1 min, 60 °C for 30 s, 72 °C for 1 min; and one cycle of final extension at 72 °C for 7 min. All of the RT-PCR experiments were repeated three times. The primers used in this study are described in Additional file 10: Table S4.
The 2B-F and adapter primers located on exon 2B and the 3’ end poly (A) tail, respectively, were used to amplify an approximately 1.1 kb fragment of the Bdtra. The amplicon (forward gene fragment) and the inverted sequence (reverse gene fragment) were cloned into pGEM®-T easy vector (Promega, USA) and transformed into E. coli (DH5α). Single-stranded RNA (ssRNA) of Bdtra was synthesized by in vitro bacterial expression using the T7-polymerase system. To produce double-stranded RNA (dsRNA), forward and reverse RNA fragments were mixed in annealing buffer (10 mM Tris pH 8 and 10 mM NaCl). Incubation was done at 75 °C for 5 min (annealing), then cooled down to room temperature (25 °C).
The Salaya1 embryos were harvested and kept for 40 min at 18 °C. Approximately 1000 precellular blastoderm embryos were dechorionated and micro-injected with 1 μg/μl of Bdtra dsRNA at the posterior end. The surviving RNAi-treated embryos were raised to the pupal stage in normal conditions. The brown and white pupae were individually isolated, awaiting sexual identification after adult emergence. The white-pupae males were categorized as pseudomales because the untreated RNAi white pupae were females with the XX sexual karyotype .
Mating behavior and fertility of the sexually mature pseudomales were tested against RNAi-treated and RNAi-untreated brown-pupae males from the Salaya1 strain. The investigation was carried out by pairing individual tested males or pseudomales with three normal sexually mature and virgin white-pupae females. Five RNAi-untreated brown-pupae males, five RNAi-treated brown-pupae males, and six RNAi-treated white-pupae pseudomales were individually tested for mating behavior and fertility when they were sexually mature. Subsequently, their testis development was inspected by microscopic dissection.
RT-PCR analysis of RNAi-treated flies
The previously described RT-PCR conditions were performed to analyze the sex-specific expression patterns of Bdtra and Bddsx from B. dorsalis (Phayathai1 and Salaya1 strains) and their RNAi-treated brown-pupae males versus white-pupae pseudomales. Three sets of primers were used in this experiment: 1B-F – 2B-R1 primers for Bdtra analysis, C3 – Female-specific (f) primers for female-specific Bddsx analysis, and C4 – Male-specific (m) primers for male-specific Bddsx analysis . All of the RT-PCR experiments were repeated three times. The primer sequences are presented in Additional file 10: Table S4.
Sequence alignment and phylogenetic tree reconstruction
Sequences were searched using the BLAST family of programs from the NCBI database for gene identification. Genomic DNA and cDNA sequences were aligned using the ClustalX2 program  and Unipro UGENE Version 1.14.0 . The genetic distance was estimated using Jukes-Cantor measurement. The phylogenetic tree of tra genes were reconstructed using the UPGMA method in the CLC Main Workbench version 7.6.2 (CLC Bio, QIAGEN). The bootstrap values were formulated from 1000 replicates with the 50 % majority rule.
This manuscript describes functional genetics of non-regulated invertebrate (domestic insect pests) thus did not require animal ethics approval.
Availability of data
All supporting data are included as additional files.
The authors gratefully acknowledge all anonymous reviewers for their suggestions and commentaries on the manuscript and Mr. Robert Bachtell Eastland for his English editing services. This research is supported by International Atomic Energy Agency research contract no. 18795 as part of the Agency’s Coordinated Research Project: Comparing Rearing Efficiency and Competitiveness of Males from Salaya Genetic Sexing Strains (Bactrocera spp.) Including the Refinement and Transfer of Existing Technology for the Improvement and Application of Strains for Area-Wide Integrated Pest Management to S. Thanaphum. Part of this study is part of the M.S. thesis of K. Laohakieat under the supervision of S. Thanaphum at the Department of Biotechnology, Faculty of Science, Mahidol University. K. Laohakieat is partially supported by a scholarship from the Faculty of Science, Mahidol University (2013–2015). We would like to thank Dr. Adisak Romsang from the Department of Biotechnology, Faculty of Science, Mahidol University for his help and Prof. Ernst A. Wimmer from the Department of Developmental Biology, Johann-Friedrich-Blumenbach-Institute of Zoology and Anthropology, Georg-August-University Goettingen, Germany for providing us with vectors #1215 and #1200.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Darwin C. Principles of sexual selection. In: The origin of species and the descent of man. 1st ed. 1977. p. 567–611.Google Scholar
- Bachtrog D, Mank JE, Peichel CL, Kirkpatrick M, Otto SP, Ashman T-L, et al. Sex determination: Why so many ways of doing it? PLoS Biol. 2014;12(7):e1001899.PubMed CentralView ArticlePubMedGoogle Scholar
- Sánchez L. Sex-determining mechanisms in insects. Int J Dev Biol. 2008;52(7):837–56.View ArticlePubMedGoogle Scholar
- Bopp D, Saccone G, Beye M. Sex determination in insects: variations on a common theme. Sex Dev. 2014;8(1–3):20–8.View ArticlePubMedGoogle Scholar
- Cline TW. Two closely linked mutations in Drosophila melanogaster that are lethal to opposite sexes and interact with daughterless. Genetics. 1978;90(4):683–97.PubMed CentralPubMedGoogle Scholar
- Gabrieli P, Falaguerra A, Siciliano P, Gomulski L, Scolari F, Zacharopoulou A, et al. Sex and the single embryo: early development in the Mediterranean fruit fly, Ceratitis capitata. BMC Dev Biol. 2010;10(1):1–14.View ArticleGoogle Scholar
- Pane A, Salvemini M, Bovi PD, 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–25.PubMedGoogle Scholar
- Morrow JL, Riegler M, Frommer M, Shearman DCA. Expression patterns of sex-determination genes in single male and female embryos of two Bactrocera fruit fly species during early development. Insect Mol Biol. 2014;23(6):754–67.View ArticlePubMedGoogle Scholar
- Morrow JL, Riegler M, Gilchrist AS, Shearman DCA, Frommer M. Comprehensive transcriptome analysis of early male and female Bactrocera jarvisi embryos. BMC Genet. 2014;15 Suppl 2:S7.PubMed CentralView ArticlePubMedGoogle Scholar
- Hilfiker-Kleiner D, Dubendorfer A, Hilfiker A, Nothiger R. Developmental analysis of two sex-determining genes, M and F, in the housefly, Musca domestica. Genetics. 1993;134(4):1187–94.PubMed CentralPubMedGoogle Scholar
- Hediger M, Henggeler C, Meier N, Perez R, Saccone G, Bopp D. Molecular characterization of the key switch F provides a basis for understanding the rapid divergence of the sex-determining pathway in the housefly. Genetics. 2010;184(1):155–70.PubMed CentralView ArticlePubMedGoogle Scholar
- Shukla JN, Palli SR. Sex determination in beetles: production of all male progeny by parental RNAi knockdown of transformer. Scientific Reports. 2012;2:602.PubMed CentralPubMedGoogle Scholar
- Verhulst EC, Lynch JA, Bopp D, Beukeboom LW, van de Zande L. A new component of the Nasonia sex determining cascade is maternally silenced and regulates transformer expression. PLoS ONE. 2013;8(5):e63618.PubMed CentralView ArticlePubMedGoogle Scholar
- Gempe T, Hasselmann M, Schiøtt M, Hause G, Otte M, Beye M. Sex determination in honeybees: two separate mechanisms induce and maintain the female pathway. PLoS Biol. 2009;7(10):e1000222.PubMed CentralView ArticlePubMedGoogle Scholar
- Kiuchi T, Koga H, Kawamoto M, Shoji K, Sakai H, Arai Y, et al. A single female-specific piRNA is the primary determiner of sex in the silkworm. Nature. 2014;509(7502):633–6.View ArticlePubMedGoogle Scholar
- Erickson JW, Quintero JJ. Indirect effects of ploidy suggest X chromosome dose, not the X:A ratio, signals sex in Drosophila. PLoS Biol. 2007;5(12), e332.PubMed CentralView ArticlePubMedGoogle Scholar
- Salz HK, Erickson JW. Sex determination in Drosophila: the view from the top. Fly. 2010;4(1):60–70.PubMed CentralView ArticlePubMedGoogle Scholar
- Schutt C, Nothiger R. Structure, function and evolution of sex-determining systems in Dipteran insects. Development. 2000;127(4):667–77.PubMedGoogle Scholar
- Keyes LN, Cline TW, Schedl P. The primary sex determination signal of Drosophila acts at the level of transcription. Cell. 1992;68(5):933–43.View ArticlePubMedGoogle Scholar
- Bell LR, Horabin JI, Schedl P, Cline TW. Positive autoregulation of Sex-lethal by alternative splicing maintains the female determined state in Drosophila. Cell. 1991;65(2):229–39.View ArticlePubMedGoogle Scholar
- Inoue K, Hoshijima K, Sakamoto H, Shimura Y. Binding of the Drosophila Sex-lethal gene product to the alternative splice site of transformer primary transcript. Nature. 1990;344(6265):461–3.View ArticlePubMedGoogle Scholar
- Valcárcel J, Singh R, Zamore PD, Green MR. The protein Sex-lethal antagonizes the splicing factor U2AF to regulate alternative splicing of transformer pre-mRNA. Nature. 1993;362(6416):171–5.View ArticlePubMedGoogle Scholar
- Bopp D, Bell LR, Cline TW, Schedl P. Developmental distribution of female-specific Sex-lethal proteins in Drosophila melanogaster. Genes Dev. 1991;5(3):403–15.View ArticlePubMedGoogle Scholar
- Willhoeft U, Franz G. Identification of the sex-determining region of the Ceratitis Capitata Y chromosome by deletion mapping. Genetics. 1996;144(2):737–45.PubMed CentralPubMedGoogle Scholar
- 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–20.View ArticlePubMedGoogle Scholar
- Lagos D, Koukidou M, Savakis C, Komitopoulou K. The transformer gene in Bactrocera oleae: the genetic switch that determines its sex fate. Insect Mol Biol. 2007;16(2):221–30.View ArticlePubMedGoogle Scholar
- Ruiz MF, Milano A, Salvemini M, Eirín-López JM, Perondini ALP, Selivon D, et al. The gene transformer of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects. PLoS ONE. 2007;2(11):e1239.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu G, Wu Q, Li J, Zhang G, Wan F. RNAi-mediated knock-down of transformer and transformer 2 to generate male-only progeny in the oriental fruit fly, Bactrocera dorsalis (Hendel). PLoS ONE. 2015;10(6):e0128892.PubMed CentralView ArticlePubMedGoogle Scholar
- Peng W, Zheng W, Handler A, Zhang H. The role of the transformer gene in sex determination and reproduction in the tephritid fruit fly, Bactrocera dorsalis (Hendel). Genetica. 2015;143(6):717–27.View ArticleGoogle Scholar
- Sarno F, Ruiz MF, Eirín-López JM, Perondini ALP, Selivon D, Sánchez L. The gene transformer-2 of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects. BMC Evol Biol. 2010;10:140.PubMed CentralView ArticlePubMedGoogle Scholar
- Saccone G, Salvemini M, Polito LC. The transformer gene of Ceratitis capitata: a paradigm for a conserved epigenetic master regulator of sex determination in insects. Genetica. 2011;139(1):99–111.View ArticlePubMedGoogle Scholar
- Burtis KC, Baker BS. Drosophila doublesex gene controls somatic sexual differentiation by producing alternatively spliced mRNAs encoding related sex-specific polypeptides. Cell. 1989;56(6):997–1010.View ArticlePubMedGoogle Scholar
- Hoshijima K, Inoue K, Higuchi I, Sakamoto H, Shimura Y. Control of doublesex alternative splicing by transformer and transformer-2 in Drosophila. Science. 1991;252(5007):833–6.View ArticlePubMedGoogle Scholar
- Shukla JN, Nagaraju J. Doublesex: a conserved downstream gene controlled by diverse upstream regulators. J Genet. 2010;89(3):341–56.View ArticlePubMedGoogle Scholar
- Hedley ML, Maniatis T. Sex-specific splicing and polyadenylation of dsx pre-mRNA requires a sequence that binds specifically to tra-2 protein in vitro. Cell. 1991;65(4):579–86.View ArticlePubMedGoogle Scholar
- Tian M, Maniatis T. A splicing enhancer complex controls alternative splicing of doublesex pre-mRNA. Cell. 1993;74(1):105–14.View ArticlePubMedGoogle Scholar
- Saccone G, Peluso I, Testa G, Di Paola F, Pane A, Polito LC. Drosophila Sex-lethal and doublesex homologous genes in Ceratitis capitata: searching for sex-specific genes to develop a medfly transgenic sexing strain, in enhancement of the sterile insect technique through genetic transformation using nuclear techniques. Vienna: IAEA/FAO; 1996.Google Scholar
- Saccone G, Pane A, Testa G, Santoro M, De Martino G, Di Paola F, et al. Sex determination in medfly: a molecular approach. In: Tan KH, editor. Area-wide control of fruit flies and other insect pests. Penung: Penerbit USM; 2000. p. 491–6.Google Scholar
- Shearman DCA, Frommer M. The Bactrocera tryoni homologue of the Drosophila melanogaster sex-determination gene doublesex. Insect Mol Biol. 1998;7(4):355–66.View ArticlePubMedGoogle Scholar
- Lagos D, Ruiz MF, Sánchez L, Komitopoulou K. Isolation and characterization of the Bactrocera oleae genes orthologous to the sex determining Sex-lethal and doublesex genes of Drosophila melanogaster. Gene. 2005;348:111–21.View ArticlePubMedGoogle Scholar
- Ruiz MF, Stefani RN, Mascarenhas RO, Perondini ALP, Selivon D, Sánchez L. The gene doublesex of the fruit fly Anastrepha obliqua (Diptera, Tephritidae). Genetics. 2005;171(2):849–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruiz M, Eirín-López J, Stefani R, Perondini AP, Selivon D, Sánchez L. The gene doublesex of Anastrepha fruit flies (Diptera, Tephritidae) and its evolution in insects. Dev Genes Evol. 2007;217(10):725–31.View ArticlePubMedGoogle Scholar
- Chen S-L, Dai S-M, Lu K-H, Chang C. Female-specific doublesex dsRNA interrupts yolk protein gene expression and reproductive ability in oriental fruit fly, Bactrocera dorsalis (Hendel). Insect Biochem Mol Biol. 2008;38(2):155–65.View ArticlePubMedGoogle Scholar
- Permpoon R, Aketarawong N, Thanaphum S. Isolation and characterization of Doublesex homologues in the Bactrocera species: B. dorsalis (Hendel) and B. correcta (Bezzi) and their putative promoter regulatory regions. Genetica. 2011;139(1):113–27.View ArticlePubMedGoogle Scholar
- White IM, Elson-Harris MM. Fruit flies of economic significance: their identification and bionomics. Oxford: CAB International; 1992. p. 187–92.Google Scholar
- Ast G. How did alternative splicing evolve? Nat Rev Genet. 2004;5(10):773–82.View ArticlePubMedGoogle Scholar
- Qi J, Su S, Mattox W. The doublesex splicing enhancer components Tra2 and Rbp1 also repress splicing through an intronic silencer. Mol Cell Biol. 2007;27(2):699–708.PubMed CentralView ArticlePubMedGoogle Scholar
- Heinrichs V, Baker BS. The Drosophila SR protein RBP1 contributes to the regulation of doublesex alternative splicing by recognizing RBP1 RNA target sequences. EMBO J. 1995;14(16):3987–4000.PubMed CentralPubMedGoogle Scholar
- Heinrichs V, Ryner LC, Baker BS. Regulation of Sex-Specific Selection of fruitless 5′ Splice Sites by transformer and transformer-2. Mol Cell Biol. 1998;18(1):450–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Juven-Gershon T, Kadonaga JT. Regulation of gene expression via the core promoter and the basal transcriptional machinery. Dev Biol. 2010;339(2):225–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Butler JEF, Kadonaga JT. The RNA polymerase II core promoter: a key component in the regulation of gene expression. Genes Dev. 2002;16(20):2583–92.View ArticlePubMedGoogle Scholar
- Kutach AK, Kadonaga JT. The downstream promoter element DPE appears to be as widely used as the TATA box in Drosophila core promoters. Mol Cell Biol. 2000;20(13):4754–64.PubMed CentralView ArticlePubMedGoogle Scholar
- Manley JL, Tacke R. SR proteins and splicing control. Genes Dev. 1996;10(13):1569–79.View ArticlePubMedGoogle Scholar
- Verhulst EC, van de Zande L, Beukeboom LW. Insect sex determination: it all evolves around transformer. Curr Opin Genet Dev. 2010;20(4):376–83.View ArticlePubMedGoogle Scholar
- Geuverink E, Beukeboom LW. Phylogenetic distribution and evolutionary dynamics of the sex determination genes doublesex and transformer in insects. Sex Dev. 2014;8(1–3):38–49.View ArticlePubMedGoogle Scholar
- Permpoon R, Thanaphum S. Isolation and characterization of oligomerization domain I and II coding regions of doublesex genes in agricultural fruit flies (Diptera: Tephritidae). Eur J Entomol. 2010;107(1):121–6.View ArticleGoogle Scholar
- Yang Y, Zhang W, Bayrer JR, Weiss MA. Doublesex and the regulation of sexual dimorphism in Drosophila melanogaster: structure, function, and mutagenesis of a female-specific domain. J Biol Chem. 2008;283(11):7280–92.View ArticlePubMedGoogle Scholar
- Isasawin S, Aketarawong N, Thanaphum S. Characterization and evaluation of microsatellite markers in a strain of the oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae), with a genetic sexing character used in sterile insect population control. Eur J Entomol. 2012;109(3):331–8.View ArticleGoogle Scholar
- Schetelig MF, Milano A, Saccone G, Handler AM. Male only progeny in Anastrepha suspensa by RNAi-induced sex reversion of chromosomal females. Insect Biochem Mol Biol. 2012;42(1):51–7.View ArticlePubMedGoogle Scholar
- Isasawin S, Aketarawong N, Lertsiri S, Thanaphum S. Development of a genetic sexing strain in Bactrocera carambolae (Diptera: Tephritidae) by introgression of sex sorting components from B. dorsalis, Salaya1 strain. BMC Genet. 2014;15 Suppl 2:S2.PubMed CentralView ArticlePubMedGoogle Scholar
- Leftwich PT, Koukidou M, Rempoulakis P, Gong H-F, Zacharopoulou A, Fu G, et al. Genetic elimination of field-cage populations of Mediterranean fruit flies. Proc Roy Soc Lond B Biol Sci. 2014. 281(1792).Google Scholar
- Koukidou M, Alphey L. Practical applications of insects’ sexual development for pest control. Sex Dev. 2014;8(1–3):127–36.View ArticlePubMedGoogle Scholar
- Ant T, Koukidou M, Rempoulakis P, Gong H-F, Economopoulos A, Vontas J, et al. Control of the olive fruit fly using genetics-enhanced sterile insect technique. BMC Biol. 2012;10:51.PubMed CentralView ArticlePubMedGoogle Scholar
- Dafa’alla T, Fu G, Alphey L. Use of a regulatory mechanism of sex determination in pest insect control. J Genet. 2010;89(3):301–5.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Schetelig MF, Handler AM. Strategy for enhanced transgenic strain development for embryonic conditional lethality in Anastrepha suspensa. Proc Natl Acad Sci U S A. 2012;109(24):9348–53.PubMed CentralView ArticlePubMedGoogle Scholar
- 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(4):425–37.View ArticlePubMedGoogle Scholar
- Scolari F, Schetelig MF, Bertin S, Malacrida AR, Gasperi G, Wimmer EA. Fluorescent sperm marking to improve the fight against the pest insect Ceratitis capitata (Wiedemann; Diptera: Tephritidae). New Biotechnol. 2008;25(1):76–84.View ArticleGoogle Scholar
- Sambrook J, Fitsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. 2nd ed. Cold Spring Harbor Press: Cold Spring Harbor; 1989.Google Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.View ArticlePubMedGoogle Scholar
- Okonechnikov K, Golosova O, Fursov M, UTENE team. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. 2012;28(8):1166–7.View ArticlePubMedGoogle Scholar