- Research article
- Open Access
RNA sequencing to characterize transcriptional changes of sexual maturation and mating in the female oriental fruit fly Bactrocera dorsalis
© Zheng et al. 2016
- Received: 2 October 2015
- Accepted: 25 February 2016
- Published: 5 March 2016
Female reproductive potential plays a significant role in the survival and stability of species, and sexual maturation and mating processes are crucial. However, our knowledge of the reproductive genes involved in sexual maturation and mating has been largely limited to model organisms. The oriental fruit fly Bactrocera dorsalis is a highly invasive agricultural pest, known to cause major economic losses; thus, it is of great value to understand the transcriptional changes involved in sexual maturation and mating processes as well as the related genes. Here, we used a high-throughput sequencing method to identify multiple genes potentially involved in sexual maturation and mating in female B. dorsalis.
We sequenced 39,999 unique genes with an average length of 883 bp. In total, 3264 differentially expressed genes (DEGs) were detected between mature virgin and immature Bactrocera dorsalis libraries, whereas only 83 DEGs were identified between flies that had mated or were mature virgins. These DEGs were functionally annotated using the GO and KEGG pathway annotation tools. Results showed that the main GO terms associated with the DEGs from the mature virgin vs. immature groups were primarily assigned to the metabolic and developmental processes, which we focused on, whereas those from the mated vs. mature virgin group largely belonged to the response to stimulus and immune system processes. Additionally, we identified multiple DEGs during sexual maturation that are involved in reproduction, and expression pattern analysis revealed that the majority DEGs detected were highly enriched in those linked to the ovaries or fat bodies. Several mating responsive genes differentially expressed after mating were also identified, and all antimicrobial peptides detected were highly enriched in fat body and significantly up-regulated approximately 2- to 10-fold at 24 h after mating.
This study supplied female reproductive genes involved in sexual maturation and the post-mating response in B. dorsalis, based on RNA-seq. Our data will facilitate molecular research related to reproduction and provide abundant target genes for effective control of this agricultural pest.
- Gene expression
Tephritid fruit flies, especially those belonging to the genus Bactrocera, are globally important due to their destructive impact on agriculture. Their high reproductive ability is an important factor that leads to enhanced survival rates and potential to multiply. Reproductive biology has been well studied in tephritid species, including Bactrocera tryoni (Froggatt), Bactrocera cucurbitae (Coquillett) and Ceratitis capitata (Wiedemann), but most reports have mainly focused on the processes at the physiological and behavioral levels [1–3]. Bactrocera dorsalis (Hendel) (the oriental fruit fly), which is one among the important Bactrocera species, is a highly invasive agricultural pest that is currently distributed across most Asian countries and in a number of Pacific Islands. Adults can lay eggs on various types of host plants, where the hatched larvae subsequently feed and cause crop loss . Some factors that regulate the fecundity and mating behavior of B. dorsalis females have been previously studied, such as methyl eugenol-fed males and multiple mates, etc. . Although the B. dorsalis transcriptome from different developmental stages of the whole body as well as specific tissue in males is available on NCBI, large-scale molecular analysis of reproduction in B. dorsalis females remains limited [6–10]. As an invasive agricultural pest with a wide host range and high fecundity, it is essential to investigate the reproductive biology of B. dorsalis females at the molecular level in order to formulate simple and effective strategies for agricultural pest control.
Recently, transcriptional changes in several genes that may be involved in female sexual maturity and mating were identified in C. capitata based on microarray data ; however, the cDNA library of C. capitata is limited to the head of female adult. Additionally, it is worth noting that the reproductive strategy of B. dorsalis is different from that of C. capitata. For instance, medfly reach sexual maturity rapidly within 2–3 days after emergence, whereas it takes approximately 2 weeks for the oriental fruit fly to complete the transformation after emergence to maturity [12, 13]. Therefore, our study is planned because both species have different mating strategy and that the methodology followed in previous study provided only limited resources.
Female reproductive potential significantly contributed to the multiply and stability of species, which mainly depends on the process of sexual maturation and mating . In most holometabolous insects, adults emerge from the pupal stage as sexually immature individuals. After they ingest a food source, they become sexually mature within several days dependent on the species . During the development of the oocyte, the ovary becomes mature and then when females reach sexual maturity they began to mate and oviposit . Until now, studies on sexual maturity have mainly concentrated on the model species whose genomes have been sequenced except C. capitata, including Drosophila melanogaster, Apis mellifera, Aedes aegypti and Blattella germanica [16–20]. However, these reports on model species mainly focused on the protein vitellogenin and information on other proteins that play roles in reproduction is lacking. While genes involved in female adult reproductive maturation, including EcR, USP and InR, have been identified in Tribolium castaneum using microarrays, they are mainly the early genes in the signal transduction pathway, and little is known about the other genes, especially the late or effect genes .
Successful mating behavior and the following egg-laying is essential for sexual reproduction in insect, especially for its multiply. Female reproductive process has been revealed to be greatly affected by male reproductive gland secretion during copulation. Reports in C. capitata and Bactrocera tryoni showed the effect of mating on remating propensity, refractory period and fecundity of females [22–24]. In D. melanogaster, it was reported that mating regulated ovulation and spawning as well as other processes in mated females, such as, inducing innate immunity, acceleration of reproductive maturation and egg production [25–28]. The sex peptide receptor, which mediates the post-mating switch in females, has been reported in D. melanogaster, Helicoverpa armigera and B. dorsalis [4, 29–31]. However, the other factors involved in this process are still largely unknown. Genome-wide research into the post-mating response in females has been performed in model species such as D. melanogaster, the honeybee queen Apis mellifera, and C. capitata, which have revealed that the post-mating response differs among species as well as the time following mating [11, 25, 28, 32–34]. Regardless, there is no such information available for B. dorsalis, and the transcriptional changes and genes involved in this process is also still unknown.
In the present study, the transcriptome of the whole body of the female B. dorsalis were sequenced at three different physiological stages (immature, mature virgins and mated individuals) using the Illumina HiSeq 2500 system, and assembled based on our previously established transcriptome of B. dorsalis . Furthermore, analysis of the DEGs between mature virgins and immature individuals was performed in order to identify the transcripts involved in the process of sexual maturation. An additional comparison between mated individuals and mature virgins was also carried out to identify the genes that respond to mating. Tissue specific expression profiles of the differentially expressed genes from these two comparisons were analyzed by qRT-PCR. Finally, the expression levels of the genes involved in the post-mating response at different time points were also obtained.
B. dorsalis were cultured in our lab under a 12 h light/dark cycle at 27 °C. Adult flies were reared on artificial diets composed of 25 % yeast extract and 75 % sugar .
Newly emerged adults were collected and sexed within 12 h after emergence, and both sexes were subsequently maintained in separate cages. As male and female oriental fruit flies require 14 days to reach sexual maturity under standard rearing conditions, 1- to 2-day-old (1–2 days) individuals were considered immature while 14-day-old (14 days) individuals were considered sexually mature. Twenty-four hours after emergence, the immature female flies were chilled briefly and RNA was extracted from the whole bodies of 10 individuals using TRIzol regent (Invitrogen, Carlsbad, CA, USA). Similarly, mature virgin samples were collected 14 days after emergence and total RNA was isolated from the whole bodies of 10 mature virgin females.
To obtain mated flies, approximately 100 fourteen-day-old virgin flies of each sex were placed into a 30 cm3 cage shortly before dark as B. dorsalis mating occurs during dusk . When copulation was observed, mating pairs were isolated and placed in small vials about 14 cm in width and 23 cm in height with 11.5 cm diameter. Notably, only pairs that mated for at least 90 min were used in further experiments in order to avoid false mating with little or no sperm transfer. Copulation terminated naturally and the mated pairs were placed in separate cages according to sex. RNA was extracted from the whole bodies of 10 mated females the following day (24 h after mating) with the same procedure as was used for the virgin flies. RNA samples were treated to remove any contaminating DNA with the DNA-free kit (Ambion, Foster City, CA, USA).
Total RNA from each sample was isolated using TRIzol reagent (Invitrogen) following the manufacturer’s instructions. The concentration and integrity of total RNA were determined using the 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Total RNA was then incubated with 10 U DNase I (Takara, Shiga, Japan) at 37 °C for 30 min, and purified using the Dynabeads® Oligo (dT)25 (Thermo Fisher, Waltham, MA, USA) following the manufacturer’s instructions. To ensure the accuracy of the results, we conducted two independent replicates for each sample, specifically immature (1, 2), mature virgin (1, 2) and mated (1, 2) individuals. For each case, the following validation described below was performed.
For cDNA library construction 100 ng purified mRNA was used using NEBNext® UltraTM RNA Library Prep Kit for Illumina sequencing (New England Biolabs, Ipswich, MA, USA). First-strand cDNA was synthesized using ProtoScript II Reverse Transcriptase. After incubation, Second Strand Synthesis Enzyme Mix was added for the synthesis of the second strand of the cDNA. Double-stranded cDNA was then purified using AMPure XP Beads (Beckman Coulter, Brea, CA, USA), followed by an end repair step using End Prep Enzyme Mix. After adaptor ligation and further purification, the cDNA library was obtained. To verify the quality of the library, it was checked using three different methods, including Qubit (Thermo Fisher) quantitative analysis, 2 % agarose gel electrophoresis, and high-sensitivity DNA chip determination.
The cDNA libraries were then used for cluster generation with the TruSeq PE Cluster Kit (Illumina, San Diego, CA, USA) and sequenced on an Illumina HiSeq™ 2500 instrument with paired-end sequencing.
Sequence assembly and annotation
Raw data were first filtered by removing adaptor sequences and low quality reads using FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/), and then assembled into contigs using Trinity (v2013-02-25) with N bases as the default assembly parameters. Finally, redundant sequences were removed while those remaining constituted the transcriptome used in downstream analyses. Furthermore, depth and coverage of sequencing was evaluated. The assembly data have been submitted to NCBI as a TSA under the accession number GEEA00000000.
All contigs were annotated with GetORF of the EMBOSS package . The ORF of each predicted protein was used in a BLASTp search against the Swiss-Prot and the NCBI nr databases, with the e-value threshold set to 10−5. Gene Ontology (GO) annotations were performed with GoPipe . Predicted proteins were first used for a BLASTp search against the Swiss-Prot and TrEMBL databases with an E-value cut-off of 10−5, and then the results were analyzed by GoPipe based with the gene2go software. The KOG (Eukaryotic Orthologous Groups) and KEGG pathway annotations were also analyzed using the Cluster of Orthologous Groups database and the Kyoto Encyclopedia of Genes and Genomes database, respectively .
Differential expression analysis
For the identification of genes that were possibly involved in reproduction, those that were differentially expressed at the various physiological stages of adults were analyzed. Two comparisons of the expression level of paired sets of the three samples were performed between mature virgin vs. immature and mated vs. mature virgin individuals. The number of reads for contigs from each sample was first converted to Reads per Kilobase per Million (RPKM), and then the MARS model (MA-plot-based method with Random Sampling model) in the DEGseq package was used to calculate the expression abundance of each contig . A false discovery rate (FDR) < 0.001 was considered to indicate significant expression abundance. Additionally, differentially expressed genes were further used for GO and KEGG analysis to identify the pathways that the DEGs were predicted to be involved in. Relativity analysis was performed to evaluate the stability and reliability of the replicates. GOEAST was used for GO term enrichment analysis . The relative significant level of enrichment of the differently expressed gene to all genes was calculated with hypergeometric statistical test method, then the GO term whose FDR < 0.001 was considered to be significant. For pathway enrichment analysis, a similar method was used. The KEGG pathway was considered as a unit, and its relative significant level of enrichment of the differently expressed pathway to whole genome was calculated.
Quantitative real-time PCR (qRT-PCR) verification
Selected DEG data was validated by qRT-PCR using the SYBR Premix Ex Taq kit (Takara) according to the manufacturer’s instructions with a real-time thermal cycler (Bio-Rad, Hercules, CA). TRIzol reagent (Invitrogen) was used to extract total RNA from B. dorsalis at three different physiological stages, including the immature (newly emergent within 24 h), mature virgin (14-day-old before mating), and mated (14-day-old post-mating). At least 10 insects were collected for each sample. The first strand cDNA was obtained from 2 μg of total RNA using M-MLV Reverse Transcriptase (Takara) with the primer oligo-anchor R (5′-GACCACGCGTATCGATGTCGACT16 (A/C/G)-3′). The primers used for qRT-PCR detection are listed in Additional file 1: Table S1. The relative gene expression data were analyzed using the 2-ΔΔCt method as described by Zheng et al. . The results were analyzed using a one-way analysis of variance (ANOVA) statistical test. All quantitative PCR experiments were repeated in three biological replicates.
Tissue specific expression pattern analysis of DEGs involved in reproduction
Total RNA was extracted from different tissues of the mature female adults, including the head, thorax, ovary, fat body and midgut. The tissue samples were dissected from 30 individuals. After the cDNA was obtained, qRT-PCR was performed according to the methods described above. A total of eleven genes differentially expressed during maturation and six AMPs whose transcript abundance changed after mating were used in this study.
Change in transcriptional abundance of mating-responsive genes
To investigate the expression level of mating-responsive genes before and after mating, newly emergent individuals (1 day after emergence) of each sex were separated and the mated pairs were obtained as described above. RNA was extracted from the head or fat body of mated females 1, 12, 24 h after mating. At least 30 individuals were used for each sample and qRT-PCR experiments were repeated in three biological replicates.
De novo sequence assembly
The cDNA samples from different stages were sequenced, and these raw data (Additional file 2: Table S2) were assembled into 39,999 unique genes (Additional file 3: Table S3). The mean contig size was 883 bp with lengths ranging from 201 to 27,791 bp (Additional file 3: Table S3 and Additional file 4: Figure S1). The contig size distribution revealed that more than half of the contigs (23,645; 59.11 %) were between 200 and 500 bp in length, whereas 34.88 % (13,951) were between 500 and 3000 bp in length (Additional file 4: Figure S1).
Gene ontology and clusters of orthologous group classification
Gene ontology (GO) analysis was performed for functional categorization of unigenes (Additional file 5: Figure S2). Within the biological processes, the main GO terms were grouped into cellular (3362; 19 %) and metabolic (2696; 15 %; Additional file 5: Figure S2A). As for the molecular function category, binding (3663; 50 %) constituted the largest group, followed by catalytic activity (2096; 29 %; Additional file 5: Figure S2B). In the previously sequenced B. dorsalis transcriptome obtained from different developmental stages, the metabolic process group (35 %) was much larger than the cellular process group (16 %) . The primary reason for this difference might be that the metabolic activities of insects during larval development is greater than in the adult stage. Many complicated physiological processes occurs during molting and metamorphosis, such as a molting cascade similar to the larval molt, histolysis of larval tissues, remodeling and formation of adult tissues, etc. [43, 44].
Assignments of clusters of orthologous groups (COG) were performed to further evaluate the completeness of the transcriptome and the effectiveness of our annotation process. The annotated sequences were grouped into 25 major functional classes (Additional file 6: Figure S3). Among them, the majority of the clusters were “Signal transduction mechanisms” (795; 13.01 %), “General function prediction only” (742; 12.14 %), whereas “Cell motility” (17; 0.28 %) and “Nuclear structure” (28; 0.46 %) represented the smallest groups (Additional file 6: Figure S3). Our transcriptome will aid further research focused on signal transduction in tephritid files.
Transcriptional changes during female maturation
Transcriptional changes of mating-responsive genes in females
To identify transcripts that may have been involved in female mating, differentially expressed genes between the mated and mature virgin females were identified. Out of 83 transcripts that displayed significant transcriptional changes between virgin mature and mated females 24 h post-mating, 65 (78 %) were more abundant in mated females, while only 18 (22 %) were more abundant in virgin mature females (FDR < 0.001, Fig. 1 & Additional file 7: File S1). This is similar with the result found in C. capitata and D. melanogaster [11, 27, 28]. Only 32 and 28 transcripts were altered in abundance 24 h after mating in C. capitata and D. melanogaster, respectively. However, in a study of Drosophila using DNA microarray analysis, 539 transcripts are differentially expressed in the lower reproductive tracts . One reason for the variable changes in transcription could be that like previous studies by Lawniczak and Begun  in Drosophila, our study focused on broader changes of gene expression in the whole body, while a study by Mack et al.  in Drosophila concentrated on the lower reproductive tracts, where significant changes take place following mating. Furthermore, the time point for sample collection after mating in both our study and that of the medfly was 24 h , that in the Lawniczak and Begun  analysis was 1–3 h, whereas the Mack et al.  examined several time points (0, 3, 6 and 24 h post-mating). It was demonstrated in previous studies that post-mating transcriptional changes are highly variable at different time points, thus possibly explaining transcriptional differences among these studies .
Statistically significant differences (p < 0.05) among enriched biological process GO terms during mating included those involved in the response to stimulus, as well as immune system, developmental, cellular, biological regulation and metabolic processes. Of these transcripts, those that are related to the response to stimulus and immune system process are the largest group (Fig. 2b). A table of these genes is provided in the supplementary data (Additional file 10: Table S5). Similar results were found in A. mellifera where the genes involved in stress response were significantly regulated following mating . Additionally, an immune response was stimulated by mating in Drosophila, and it might it arise from sexually antagonistic interaction between the sexes [46, 47]. Furthermore, enrichment analysis of KEGG pathways for the mating responsive genes was also performed (Additional file 9: File S2). Our results showed that the metabolism pathways were the main groups, including carbohydrate metabolism, transport and catabolism, glycan biosynthesis and metabolism, lipid metabolism, metabolism of cofactors and vitamins. The different results between KEGG enrichment and GO enrichment might be caused by the different database used in the analysis.
Genes involved in oogenesis and ovary maturation
hts was reported to modulate actin polymerization in neurons, sharing a conserved function with mammalian adducing in actin-capping activity . One isoform of hts, Ovhts-RC, was shown to be distributed along a cortical anterior-posterior gradient during Drosophila late-oogenesis to modulate spatially restricted actin filament growth at the oocyte cortex while alterations in Ovhts-RC led to actin overgrowth in the oocyte . A previous report indicated that oskar participated in regulating pole plasm assembly, and a long form of oskar played an important role in yolk endocytosis and F-actin projection at the posterior pole . Recently, it was found that oskar functions upstream of D-EndoB in the polarized endocytic activity of the yolk uptake process in Drosophila oocytes . It was also found that mago nashi interacts with Tsunagi/Y14 and other proteins, forming a multiprotein complex to play an essential role in localization of oskar mRNA within the oocyte in Drosophila . Recent analyses demonstrated that mago nashi forms a complex with Tsunagi/Y14 and Ranshi to influence oocyte differentiation within the posterior pole of the presumptive oocyte in Drosophila . Vitellogenin, also known as yolk protein, has been well studied in insects, especially in Drosophila, Aedes aegypti, Anopheles gambiae and C. capitata [54, 59–62]. In these species, it is expressed in the fat body, secreted into the hemolymph and subsequently sequestered by oocytes to facilitate the transport of carbohydrates, lipids and other nutrients to the ovaries . In the transcriptome we obtained, three vitellogenin genes were identified, including vitellogenin 1, vitellogenin 2 and vitellogenin 2 precursor, all of which are up-regulated in sexually mature females, indicating that they are essential for Bactrocera sexual maturation. The gene yolkless, commonly known as the vitellogenin receptor, belongs to the apolipoproteins receptor (LDLR) superfamily . Drosophila yolkless is expressed very early during the development of the oocyte at both transcriptional and translational level, long before vitellogenesis begins. Previous studies found that yolkless mutants failed to incorporate yolk proteins during oogenesis . Our findings are validated in the qRT-PCR results, revealing that these genes were enriched in the ovary or fat body.
Furthermore, we found disheveled and axin (isoform A) in our transcriptome, both of which are core components of Wnt signaling pathways (Additional file 8: Table S4). They are also required for cell movement and FAK regulation during ovarian morphogenesis, which may explain why both genes are abundantly expressed in mature B. dorsalis females  (Fig. 3).
Notably, two genes involved in sex determination were also identified, including sex-determining region y protein (sry) and transformer-2 (tra-2). It has been reported that sry plays a role in mammalian sex determination, reduced or delayed sry expression impairs testis development. This indicates a potential role of sry dysregulation in human intersex disorders . Additionally, tra-2 was also reported to be an essential switch for the sex determination cascade in dipteran insect species [50, 51]. Knockdown of tra-2 leads to the production of male-only progeny in B. dorsalis . In the present study, we identified these two genes and found they were up-regulated in sexually mature females compared to their newly eclosed counterparts, which is the first evidence that sex determination related genes may also function during sex maturation.
Genes involved in the post-mating response of females
This study significantly advances the characterization of the sexual maturation and mating process at the molecular level in the oriental fruit fly, a major invasive agricultural tephritid pest. Numerous transcripts displayed significant changes in abundance during maturation of the female oriental fruit fly and generally mirrored the physiological and behavioral alterations in female flies. Multiple reproductive genes were differentially expressed between mature virgin and immature individuals, notably including two genes involved in sex determination, possibly indicating they have important roles during sexual maturation. However, modest post-mating changes were also observed in females. The difference in transcriptional changes pertaining to sexual maturation and mating, compared to C. capitata, is perhaps a result of different reproductive strategies employed by these species.
Our data reveals the existence of additional potential molecular targets for B. dorsalis control, and may aid in future studies of reproductive mechanisms in fruit flies. Furthermore, our findings provide a comprehensive resource for exploring the complex molecular mechanisms driving the evolution of reproductive processes, and will likely facilitate the development of effective and eco-friendly pest control strategies.
This Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GEEA00000000. Other data sets supporting the results of this article are included in the additional files.
This research was supported by grants from National Natural Science Foundation of China (No: 31201516 and 31471774), project No. 2662015PY069 supported by the Fundamental Research Funds for the Central Universities, and the earmarked fund for the China Agriculture Research System (No. CARS-27).
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.
- Kuba H, Itô Y. Remating inhibition in the melon fly, Bactrocera (= Dacus) cucurbitae (Diptera: Tephritidae): copulation with spermless males inhibits female remating. J Ethol. 1993;11(1):23–8.View ArticleGoogle Scholar
- Harmer AM, Radhakrishnan P, Taylor PW. Remating inhibition in female Queensland fruit flies: effects and correlates of sperm storage. J Insect Physiol. 2006;52(2):179–86.View ArticlePubMedGoogle Scholar
- Miyatake T, Chapman T, Partridge L. Mating-induced inhibition of remating in female Mediterranean fruit flies Ceratitis capitata. J Insect Physiol. 1999;45(11):1021–8.View ArticlePubMedGoogle Scholar
- Zhang H, Li H. Photographic guide to key control techniques for citrus disease and insect pests. Beijing: Chinese Agricultural Press; 2012.Google Scholar
- Shelly T. Fecundity of female oriental fruit flies (Diptera: Tephritidae): effects of methyl eugenol-fed and multiple mates. Ann Entomol Soc Am. 2000;93(3):559–64.View ArticleGoogle Scholar
- Zheng W, Peng T, He W, Zhang H. High-throughput sequencing to reveal genes involved in reproduction and development in Bactrocera dorsalis (Diptera: Tephritidae). PLoS One. 2012;7(5):e36463.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen SL, Dai SM, Lu KH, 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
- Wei D, Tian CB, Liu SH, Wang T, Smagghe G, Jia FX, Dou W, Wang JJ. Transcriptome analysis to identify genes for peptides and proteins involved in immunity and reproduction from male accessory glands and ejaculatory duct of Bactrocera dorsalis. Peptides. 2015. http://dx.doi.org/10.1016/j.peptides.2015.08.007.
- Wei D, Li HM, Yang WJ, Wei DD, Dou W, Huang Y, Wang JJ. Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis. Insect Mol Biol. 2015;24(1):41–57.Google Scholar
- Suganya R, Chen SL, Lu KH. Target of rapamycin in the oriental fruit fly Bactrocera dorsalis (Hendel): its cloning and effect on yolk protein expression. Arch Insect Biochem Physiol. 2010;75(1):45–56.View ArticlePubMedGoogle Scholar
- Gomulski LM, Dimopoulos G, Xi Z, Scolari F, Gabrieli P, Siciliano P, Clarke AR, Malacrida AR, Gasperi G. Transcriptome profiling of sexual maturation and mating in the Mediterranean fruit fly, Ceratitis capitata. PLoS One. 2012;7(1):e30857.Google Scholar
- Hong X, Ding J. Agricultural entomology. 2nd ed. Beijing: China Agriculture Press; 2007.Google Scholar
- Papadopoulos N, Katsoyannos B, Kouloussis N, Economopoulos A, Carrey J. Effect of adult age, food, and time of day on sexual calling incidence of wild and mass-reared Ceratitis capitata males. Entomol Exp Appl. 1998;89(2):175–82.View ArticleGoogle Scholar
- Lei CL, Rong XL. General entomology. 3rd ed. Beijing: Chinese Agricultural Press; 2007.Google Scholar
- Riddiford LM. How does juvenile hormone control insect metamorphosis and reproduction? Gen Comp Endocrinol. 2012;179(3):477–84.View ArticlePubMedGoogle Scholar
- Hyun S. Body size regulation and insulin-like growth factor signaling. Cell Mol Life Sci. 2013;70(13):2351–65.View ArticlePubMedGoogle Scholar
- Corona M, Velarde RA, Remolina S, Moran-Lauter A, Wang Y, Hughes KA, Robinson GE. Vitellogenin, juvenile hormone, insulin signaling, and queen honey bee longevity. Proc Natl Acad Sci USA. 2007;104(17):7128–33.Google Scholar
- Sun G, Zhu J, Li C, Tu Z, Raikhel AS. Two isoforms of the early E74 gene, an Ets transcription factor homologue, are implicated in the ecdysteroid hierarchy governing vitellogenesis of the mosquito, Aedes aegypti. Mol Cell Endocrinol. 2002;190(1–2):147–57.View ArticlePubMedGoogle Scholar
- Cruz J, Mane-Padros D, Zou Z, Raikhel AS. Distinct roles of isoforms of the heme-liganded nuclear receptor E75, an insect ortholog of the vertebrate Rev-erb, in mosquito reproduction. Mol Cell Endocrinol. 2012;349(2):262–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Suren-Castillo S, Abrisqueta M, Maestro JL. FoxO inhibits juvenile hormone biosynthesis and vitellogenin production in the German cockroach. Insect Biochem Mol Biol. 2012;42(7):491–8.View ArticlePubMedGoogle Scholar
- Parthasarathy R, Palli SR. Molecular analysis of nutritional and hormonal regulation of female reproduction in the red flour beetle, Tribolium castaneum. Insect Biochem Mol Biol. 2011;41(5):294–305.View ArticlePubMedPubMed CentralGoogle Scholar
- Marchini D, Del Bene G, Cappelli L, Dallai R. Ultrastructure of the male reproductive accessory glands in the medfly Ceratitis capitata (Diptera:Tephritidae) and preliminary characterization of their secretions. Arthropod Struct Dev. 2003;31(4):313–27.View ArticlePubMedGoogle Scholar
- Radhakrishnan P, Nair S, Raftos D, Taylor PW. Transfer and fate of male ejaculate in female Queensland fruit flies. Physiol Entomol. 2008;33:302–9.View ArticleGoogle Scholar
- Radhakrishnan P, Taylor PW. Ability of male Queensland fruit flies to inhibit receptivity in multiple mates, and the associated recovery of accessory glands. J Insect Physiol. 2008;54(2):421–8.View ArticlePubMedGoogle Scholar
- Mack PD, Kapelnikov A, Heifetz Y, Bender M. Mating-responsive genes in reproductive tissues of female Drosophila melanogaster. Proc Natl Acad Sci U S A. 2006;103(27):10358–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Gillott C. Male accessory gland secretions: modulators of female reproductive physiology and behavior. Annu Rev Entomol. 2003;48:163–84.View ArticlePubMedGoogle Scholar
- McGraw LA, Gibson G, Clark AG, Wolfner MF. Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster. Curr Biol. 2004;14(16):1509–14.View ArticlePubMedGoogle Scholar
- Lawniczak MK, Begun DJ. A genome-wide analysis of courting and mating responses in Drosophila melanogaster females. Genome. 2004;47(5):900–10.View ArticlePubMedGoogle Scholar
- Zheng WP, Liu YR, Zheng WW, Xiao YL, Zhang HY. Influence of the silencing sex-peptide receptor on Bactrocera dorsalis adults and offspring by feeding with ds-spr. J Asia-Pac Entomol. 2015;18(3):477–81.View ArticleGoogle Scholar
- Hanin O, Azrielli A, Zakin V, Applebaum S, Rafaeli A. Identification and differential expression of a sex-peptide receptor in Helicoverpa armigera. Insect Biochem Mol Biol. 2011;41(8):537–44.View ArticlePubMedGoogle Scholar
- Hanin O, Azrielli A, Applebaum SW, Rafaeli A. Functional impact of silencing the Helicoverpa armigera sex-peptide receptor on female reproductive behaviour. Insect Mol Biol. 2012;21(2):161–7.View ArticlePubMedGoogle Scholar
- Manfredini F, Brown MJ, Vergoz V, Oldroyd BP. RNA-sequencing elucidates the regulation of behavioural transitions associated with the mating process in honey bee queens. BMC Genomics. 2015;16:563.View ArticlePubMedPubMed CentralGoogle Scholar
- Kocher SD, Richard FJ, Tarpy DR, Grozinger CM. Genomic analysis of post-mating changes in the honey bee queen (Apis mellifera). BMC Genomics. 2008;9:232.View ArticlePubMedPubMed CentralGoogle Scholar
- Kocher SD, Tarpy DR, Grozinger CM. The effects of mating and instrumental insemination on queen honey bee flight behaviour and gene expression. Insect Mol Biol. 2010;19(2):153–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Li XX, Zhang MY, Zhang HY. RNA interference of four genes in adult Bactrocera dorsalis by feeding their dsRNAs. PLoS One. 2011;6(3):e17788.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng WW, Peng W, Zhu CP, Zhang Q, Saccone G, Zhang HY. Identification and expression profile analysis of odorant binding proteins in the oriental fruit fly Bactrocera dorsalis. Int J Mol Sci. 2013;14(7):14936–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000;16(6):276–7.View ArticlePubMedGoogle Scholar
- Gene Ontology C, Blake JA, Dolan M, Drabkin H, Hill DP, Li N, Sitnikov D, Bridges S, Burgess S, Buza T et al. Gene Ontology annotations and resources. Nucleic Acids Res. 2013;41(D1):D530–5.Google Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010;38 Suppl 1:D355–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.View ArticlePubMedGoogle Scholar
- Zheng Q, Wang XJ. GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008;36(s2):W358–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng WW, Zhu CP, Peng T, Zhang HY. Odorant receptor co-receptor Orco is upregulated by methyl eugenol in male Bactrocera dorsalis (Diptera: Tephritidae). J Insect Physiol. 2012;58(8):1122–7.View ArticlePubMedGoogle Scholar
- Mesce KA, Fahrbach SE. Integration of endocrine signals that regulate insect ecdysis. Front Neuroendocrinol. 2002;23(2):179–99.View ArticlePubMedGoogle Scholar
- Lee CY, Wendel DP, Reid P, Lam G, Thummel CS, Baehrecke EH. E93 directs steroid-triggered programmed cell death in Drosophila. Mol Cell. 2000;6(2):433–43.View ArticlePubMedGoogle Scholar
- Chapman T, Liddle LF, Kalb JM, Wolfner MF, Partridge L. Cost of mating in Drosophila melanogaster females is mediated by male accessory gland products. Nature. 1995;373(6511):241–4.View ArticlePubMedGoogle Scholar
- Peng J, Zipperlen P, Kubli E. Drosophila sex-peptide stimulates female innate immune system after mating via the Toll and Imd pathways. Curr Biol. 2005;15(18):1690–4.View ArticlePubMedGoogle Scholar
- Innocenti P, Morrow EH. Immunogenic males: a genome-wide analysis of reproduction and the cost of mating in Drosophila melanogaster females. J Evol Biol. 2009;22(5):964–73.View ArticlePubMedGoogle Scholar
- Cohen ED, Mariol MC, Wallace RM, Weyers J, Kamberov YG, Pradel J, Wilder EL. DWnt4 regulates cell movement and focal adhesion kinase during Drosophila ovarian morphogenesis. Dev Cell. 2002;2(4):437–48.Google Scholar
- Larney C, Bailey TL, Koopman P. Switching on sex: transcriptional regulation of the testis-determining gene Sry. Development. 2014;141(11):2195–205.View ArticlePubMedPubMed CentralGoogle Scholar
- 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–25.PubMedGoogle 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
- Pokrywka NJ, Zhang H, Raley-Susman K. Distinct roles for hu li tai shao and swallow in cytoskeletal organization during Drosophila oogenesis. Dev Dynam. 2014;243(7):906–16.View ArticleGoogle Scholar
- Lewandowski JP, Sheehan KB, Bennett Jr PE, Boswell RE. Mago Nashi, Tsunagi/Y14, and Ranshi form a complex that influences oocyte differentiation in Drosophila melanogaster. Dev Biol. 2010;339(2):307–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Rina M, Savakis C. A cluster of vitellogenin genes in the Mediterranean fruit fly Ceratitis capitata: sequence and structural conservation in dipteran yolk proteins and their genes. Genetics. 1991;127(4):769–80.PubMedPubMed CentralGoogle Scholar
- Pielage J, Bulat V, Zuchero JB, Fetter RD, Davis GW. Hts/Adducin controls synaptic elaboration and elimination. Neuron. 2011;69(6):1114–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Vanzo N, Oprins A, Xanthakis D, Ephrussi A, Rabouille C. Stimulation of endocytosis and actin dynamics by Oskar polarizes the Drosophila oocyte. Dev Cell. 2007;12(4):543–55.View ArticlePubMedGoogle Scholar
- Tsai YC, Chiang W, Liou W, Lee WH, Chang YW, Wang PY, Li YC, Tanaka T, Nakamura A, Pai LM. Endophilin B is required for the Drosophila oocyte to endocytose yolk downstream of Oskar. Development. 2014;141(3):563–73.Google Scholar
- Hachet O, Ephrussi A. Drosophila Y14 shuttles to the posterior of the oocyte and is required for oskar mRNA transport. Curr Biol. 2001;11(21):1666–74.View ArticlePubMedGoogle Scholar
- Goldman TD, Arbeitman MN. Genomic and functional studies of Drosophila sex hierarchy regulated gene expression in adult head and nervous system tissues. PLoS Genet. 2007;3(11):e216.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujii S, Amrein H. Genes expressed in the Drosophila head reveal a role for fat cells in sex-specific physiology. EMBO J. 2002;21(20):5353–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Roy SG, Hansen IA, Raikhel AS. Effect of insulin and 20-hydroxyecdysone in the fat body of the yellow fever mosquito, Aedes aegypti. Insect Biochem Mol Biol. 2007;37(12):1317–26.View ArticlePubMedPubMed CentralGoogle Scholar
- Rogers DW, Whitten MM, Thailayil J, Soichot J, Levashina EA, Catteruccia F. Molecular and cellular components of the mating machinery in Anopheles gambiae females. Proc Natl Acad Sci U S A. 2008;105(49):19390–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Tufail M, Takeda M. Molecular characteristics of insect vitellogenins. J Insect Physiol. 2008;54(12):1447–58.View ArticlePubMedGoogle Scholar
- Schonbaum CP, Lee S, Mahowald AP. The Drosophila yolkless gene encodes a vitellogenin receptor belonging to the low density lipoprotein receptor superfamily. Proc Natl Acad Sci U S A. 1995;92(5):1485–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Schonbaum CP, Perrino JJ, Mahowald AP. Regulation of the vitellogenin receptor during Drosophila melanogaster oogenesis. Mol Biol Cell. 2000;11(2):511–21.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar