BACKGROUND
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 , bactrocera cucurbitae  and ceratitis capitata , but most reports have mainly focused on the processes at the physiological and behavioral levels .
bactrocera dorsalis  , 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  <cit> .
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.
<cit> .
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 .
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  <cit> ; 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– <dig> days after emergence, whereas it takes approximately 2 weeks for the oriental fruit fly to complete the transformation after emergence to maturity  <cit> .
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  <cit> .
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  <cit> .
during the development of the oocyte, the ovary becomes mature and then when females reach sexual maturity they began to mate and oviposit  <cit> .
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 .
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  <cit> .
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 .
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 .
the sex peptide receptor, which mediates the post-mating switch in females, has been reported in d. melanogaster, helicoverpa armigera and b. dorsalis .
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 .
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  using the illumina hiseq  <dig> system, and assembled based on our previously established transcriptome of b. dorsalis  <cit> .
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.
methods
insect culture
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  <cit> .
sample preparation
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  individuals were considered immature while 14-day-old  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  <dig> individuals using trizol regent .
similarly, mature virgin samples were collected 14 days after emergence and total rna was isolated from the whole bodies of  <dig> mature virgin females.
to obtain mated flies, approximately  <dig> fourteen-day-old virgin flies of each sex were placed into a 30 cm <dig> cage shortly before dark as b. dorsalis mating occurs during dusk  <cit> .
when copulation was observed, mating pairs were isolated and placed in small vials about 14 cm in width and 23 cm in height with  <dig>  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  <dig> mated females the following day  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 .
sequencing
total rna from each sample was isolated using trizol reagent  following the manufacturer’s instructions.
the concentration and integrity of total rna were determined using the  <dig> bioanalyzer .
total rna was then incubated with  <dig> u dnase i  at 37 °c for 30 min, and purified using the dynabeads® oligo  <dig>  following the manufacturer’s instructions.
to ensure the accuracy of the results, we conducted two independent replicates for each sample, specifically immature , mature virgin  and mated  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 .
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 , 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  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  and sequenced on an illumina hiseq™  <dig> 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 , and then assembled into contigs using trinity  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 geea <dig>
all contigs were annotated with getorf of the emboss package  <cit> .
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− <dig>  gene ontology  annotations were performed with gopipe  <cit> .
predicted proteins were first used for a blastp search against the swiss-prot and trembl databases with an e-value cut-off of 10− <dig>  and then the results were analyzed by gopipe based with the gene2go software.
the kog  and kegg pathway annotations were also analyzed using the cluster of orthologous groups database and the kyoto encyclopedia of genes and genomes database, respectively  <cit> .
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 , and then the mars model  in the degseq package was used to calculate the expression abundance of each contig  <cit> .
a false discovery rate  <  <dig>  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  <cit> .
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 <  <dig>  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  verification
selected deg data was validated by qrt-pcr using the sybr premix ex taq kit  according to the manufacturer’s instructions with a real-time thermal cycler .
trizol reagent  was used to extract total rna from b. dorsalis at three different physiological stages, including the immature , mature virgin , and mated .
at least  <dig> insects were collected for each sample.
the first strand cdna was obtained from 2 μg of total rna using m-mlv reverse transcriptase  with the primer oligo-anchor r -3′).
the primers used for qrt-pcr detection are listed in additional file 1: table s <dig>  the relative gene expression data were analyzed using the 2-ΔΔct method as described by zheng et al.
<cit> .
the results were analyzed using a one-way analysis of variance  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  <dig> 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  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  <dig>   <dig>  24 h after mating.
at least  <dig> individuals were used for each sample and qrt-pcr experiments were repeated in three biological replicates.
RESULTS
de novo sequence assembly
the cdna samples from different stages were sequenced, and these raw data  were assembled into  <dig>  unique genes .
the mean contig size was 883 bp with lengths ranging from  <dig> to  <dig>  bp .
the contig size distribution revealed that more than half of the contigs  were between  <dig> and 500 bp in length, whereas  <dig>  %  were between  <dig> and 3000 bp in length .
gene ontology and clusters of orthologous group classification
gene ontology  analysis was performed for functional categorization of unigenes .
within the biological processes, the main go terms were grouped into cellular  and metabolic .
as for the molecular function category, binding  constituted the largest group, followed by catalytic activity .
in the previously sequenced b. dorsalis transcriptome obtained from different developmental stages, the metabolic process group  was much larger than the cellular process group   <cit> .
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.
<cit> .
assignments of clusters of orthologous groups  were performed to further evaluate the completeness of the transcriptome and the effectiveness of our annotation process.
the annotated sequences were grouped into  <dig> major functional classes .
among them, the majority of the clusters were “signal transduction mechanisms” , “general function prediction only” , whereas “cell motility”  and “nuclear structure”  represented the smallest groups .
our transcriptome will aid further research focused on signal transduction in tephritid files.
transcriptional changes during female maturation
to identify transcripts that may play a role in female maturation, differentially expressed gene sequences between the mature and immature libraries were identified with a total of  <dig> degs .
among them,  <dig> transcripts were more abundant in the mature females, while  <dig> were more abundant in immature females.
our results are different from those revealed by the adult transcriptome of c. capitata obtained with a microarray, in which only  <dig> transcripts displayed significant changes during female maturation.
furthermore, in the transcriptome of c. capitata, the transcripts with greater abundance were found in the immature females  instead of the mature individuals   <cit> .
these differences may be attributed to the use of different sequencing methods and the different sexual maturation time utilized by these two species.
for example, b. dorsalis is sexually mature after 14 days, while c. capitata can mate after only 4 days  <cit> .
therefore, more genes are possibly implicated in the sexual maturation in b. dorsalis.fig.
<dig> comparisons of differentially expressed genes  between the mature virgin and immature samples, as well as between the mated and mature virgin samples.
the number of up-regulated genes is represented in the red column, while down-regulated genes are shown in green
analysis of enriched biological process go terms with significant differences  during female maturation was performed.
we focused on the physiological processes that appear to play essential roles in maturation, including reproductive, signaling, developmental, metabolic, immune processes as well as response to stimulus .
a table of the identified genes is provided in the supplementary data .
results revealed that differentially expressed genes in the metabolic and developmental processes are a much larger group than in other processes, which is similar to the results in c. capitata that indicated the most enriched terms in mature females were categorized in the metabolic process group  <cit> .
these data indicate that a variety of metabolic activities and developmental affairs usually occur during oogenesis and ovary maturation.
additionally, enrichment analysis of the kegg pathways for the differentially expressed genes during maturation was also performed.
results showed that the number of differentially expressed genes in the metabolism pathways are the largest groups, including carbohydrate metabolism, energy metabolism, nucleotide metabolism, lipid metabolism, metabolism of cofactors and vitamins, amino acid metabolism, etc.
.
this is similar with the result of go enrichment analysis.fig.
<dig> go classification of the degs to demonstrate changes in transcript abundance with maturation and mating.
goeast was used for go term enrichment analysis, and 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 <  <dig>  was considered to be significant enriched.
the number of enriched transcripts for go annotated sequences was revealed in the various samples of different biological categories.
a mature virgin females compared to immature females, and b mated females compared to mature virgin females
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  <dig> transcripts that displayed significant transcriptional changes between virgin mature and mated females 24 h post-mating,  <dig>  were more abundant in mated females, while only  <dig>  were more abundant in virgin mature females .
this is similar with the result found in c. capitata and d. melanogaster  <cit> .
only  <dig> and  <dig> 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,  <dig> transcripts are differentially expressed in the lower reproductive tracts  <cit> .
one reason for the variable changes in transcription could be that like previous studies by lawniczak and begun  <cit>  in drosophila, our study focused on broader changes of gene expression in the whole body, while a study by mack et al.
<cit>  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  <cit> , that in the lawniczak and begun  <cit>  analysis was 1–3 h, whereas the mack et al.
<cit>  examined several time points .
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  <cit> .
statistically significant differences  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 .
a table of these genes is provided in the supplementary data .
similar results were found in a. mellifera where the genes involved in stress response were significantly regulated following mating  <cit> .
additionally, an immune response was stimulated by mating in drosophila, and it might it arise from sexually antagonistic interaction between the sexes  <cit> .
furthermore, enrichment analysis of kegg pathways for the mating responsive genes was also performed .
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
from deg analysis of the samples from immature  vs. mature virgin  individuals, we identified seven genes that reportedly participate in oogenesis and ovary maturation in insects .
interestingly, two genes that were previously found to play important roles in wnt signaling pathways and an additional two involved in sex determination were also identified in this deg library .
the subsequent qrt-pcr validated this result .
all of these genes were significantly highly expressed at the sexually mature and post-mating stage, which is consistent with our transcriptome data, indicating that they are involved in physiological processes of mature females as is the case in other insects  .
additionally, we analyzed expression of these genes in several tissues to elucidate their roles and found that most are highly enriched in the ovary, including hu-li tai shao , oskar, mago nashi, yolkless , disheveled, axin, sry and transformer- <dig> .
the expression level of oskar in the ovary is approximately  <dig> fold higher than that in the fat body and  <dig> fold in the midgut.
notably, three vitellogenin genes can be detected in the head, thorax, ovary, fat body, while they are more abundantly expressed in fat body .
these results are explained below by analyzing their roles determined for other insects.fig.
<dig> qrt-pcr confirmation of the differentially expressed genes between mature virgin females and the immature individuals.
total rnas from b. dorsalis at three different physiological stages were isolated, including the immature , virgin  and mated .
relative transcript levels are calculated by real-time pcr using the 16 s rrna gene as the standard control.
different letters represent significant differences  among samples.
three independent biological replicates were performed for each sample
fig.
<dig> differentially expressed genes between mature virgin females and immature individuals in various tissues.
qrt-pcr was used to analyze mrna levels of differentially expressed genes, including those found in the head, thorax, ovary, fat body and midgut.
relative transcript levels are calculated by real-time pcr using the 16 s rrna gene as the standard control.
various letters indicate significant differences in the expression level .
three biological replicates were performed
hts was reported to modulate actin polymerization in neurons, sharing a conserved function with mammalian adducing in actin-capping activity  <cit> .
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  <cit> .
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  <cit> .
recently, it was found that oskar functions upstream of d-endob in the polarized endocytic activity of the yolk uptake process in drosophila oocytes  <cit> .
it was also found that mago nashi interacts with tsunagi/y <dig> and other proteins, forming a multiprotein complex to play an essential role in localization of oskar mrna within the oocyte in drosophila  <cit> .
recent analyses demonstrated that mago nashi forms a complex with tsunagi/y <dig> and ranshi to influence oocyte differentiation within the posterior pole of the presumptive oocyte in drosophila  <cit> .
vitellogenin, also known as yolk protein, has been well studied in insects, especially in drosophila, aedes aegypti, anopheles gambiae and c. capitata .
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  <cit> .
in the transcriptome we obtained, three vitellogenin genes were identified, including vitellogenin  <dig>  vitellogenin  <dig> and vitellogenin  <dig> 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  superfamily  <cit> .
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  <cit> .
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  in our transcriptome, both of which are core components of wnt signaling pathways .
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  <cit>  .
notably, two genes involved in sex determination were also identified, including sex-determining region y protein  and transformer- <dig> .
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  <cit> .
additionally, tra- <dig> was also reported to be an essential switch for the sex determination cascade in dipteran insect species  <cit> .
knockdown of tra- <dig> leads to the production of male-only progeny in b. dorsalis  <cit> .
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
upon mating, significant expression changes occurred in genes involved in the immune or stress response .
the immune response stimulated by mating might arise from sexually antagonistic interaction between the sexes  <cit> .
several up-regulated genes from the deg analysis showed significant increases in transcript abundance in females, including two prrs  and six antimicrobial peptides   .
to investigate the expression pattern of these genes, the expression levels of amps in different tissues were analyzed by qrt-pcr.
the results showed that all were highly enriched in the fat body, which is usually an important organ involved in the immune response in insects .
the expression pattern of these genes indicates that they may participate in the immune response after mating in b. dorsalis.fig.
<dig> expression pattern of the immune-related genes between mated females and the mature virgin individuals in different tissues, including the head, thorax, ovary, fat body and midgut.
relative transcript levels of six antimicrobial peptide are calculated by real-time pcr using the 16 s rrna gene as the standard control.
different letters indicate significant differences  among samples and three independent biological replicates were performed for each sample
to further validate, the changes in expression at different time points after mating were analyzed in both the fat body and head.
all genes except attacin-c-like and phormicin-like were significantly up-regulated 24 h after mating in the fat body approximately 2– <dig> fold .
this is partly supported by the previous studies in d. melanogaster, which showed that the immune response is activated in mated females by mating and associated sex peptides transfer  <cit> .
upon mating, potent signal molecules from the males enter the females and the behavior and physiology of the mated females shifted from reproductive receptivity to fecundity.
the immune response was induced in the females and it might arise from sexually antagonistic interaction between the sexes  <cit> .
however, this is different from the results of the previous report on c. capitata, which revealed no significant changes of transcript abundance were detected in immune-related genes after mating, except for a large reduction in defensin in the female abdomen  <cit> .
this data suggests that post mating changes differ among species  <cit> .
additionally, the expression level of most amps did not significantly change during the time course in the female head 24 h post-mating, except for sapecin and cecropin- <dig>  both of which significantly increased in transcript abundance, by  <dig>  and  <dig>  fold, respectively .
it is clear that the fat body is the main tissue where the immune response takes place and not the head.fig.
<dig> the effect of mating on the expression level of immune-related genes in the fat body  and head  of b. dorsalis females.
rna was extracted from the head or fat body of mature vigin females as well as mated females at  <dig>   <dig>  24 h after mating.
relative transcript levels of six antimicrobial peptide are calculated by real-time pcr using the 16 s rrna gene as the standard control.
various letters represent significant differences of the expression level of genes .
three biological replicates were performed
CONCLUSIONS
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.
availability of supporting data
this transcriptome shotgun assembly project has been deposited at ddbj/embl/genbank under the accession geea <dig>  other data sets supporting the results of this article are included in the additional files.
additional files
additional file 1: table s <dig>  primers used for qrt-pcr.
additional file 2: table s <dig>  summary of the solexa sequencing statistics of b. dorsalis adult transcriptome.
additional file 3: table s <dig>  summary of the b. dorsalis adult transcriptome.
additional file 4: figure s <dig>  assembled contig length distribution of the transcriptome.
the x-axis shows contig size and the y-axis indicates the number of contigs for each given size.
additional file 5: figure s <dig>  gene ontology  analyses of transcriptome sequences predicted to participate in:  biological processes and  molecular functions.
data are presented as level  <dig> go categorization for molecular function and level  <dig> go categorization for biological process.
classified gene objects are showed as percentages  of the total number according to the go assignments.
additional file 6: figure s <dig>  histogram presentation of clusters of orthologous groups  classification.
a total of  <dig> predicted unigenes were classified in  <dig> cog categories.
additional file 7: file s <dig>  genes differentially expressed between mature virgin females and immature individuals, and between mated females and mature virgins.
reads, rpkm and the fold values for each comparison group  is shown.
additional file 8: table s <dig>  oriental fruit fly assembled sequences of degs  with best-hit matches to dipteran genes involved in the sexual maturation.
additional file 9: file s <dig>  enrichment analysis of kegg pathways for the differentially expressed genes.
additional file 10: table s <dig>  oriental fruit fly assembled sequences of degs  with best-hit matches to dipteran genes involved in the immune system process or response to stimulus.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
the study was conceived and designed by wwz.
dyl and fyw prepared the oriental fruit fly samples.
fyw conducted the rna extraction.
wwz performed the qrt-pcr.
wwz, jlw and hyz analyzed the transcriptomic data.
wwz prepared the initial draft of the manuscript and hyz modified the final version.
all authors have read and approved the final version of the manuscript.

 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, bactrocera cucurbitae and ceratitis capitata, but most reports have mainly focused on the processes at the physiological and behavioral levels. bactrocera dorsalis, 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. 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– days after emergence, whereas it takes approximately 2 weeks for the oriental fruit fly to complete the transformation after emergence to maturity. 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. 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. 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. the sex peptide receptor, which mediates the post-mating switch in females, has been reported in d. melanogaster, helicoverpa armigera and b. dorsalis. 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. 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 using the illumina hiseq 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. methods insect culture 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. sample preparation 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 individuals were considered immature while 14-day-old 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 individuals using trizol regent. similarly, mature virgin samples were collected 14 days after emergence and total rna was isolated from the whole bodies of mature virgin females. to obtain mated flies, approximately fourteen-day-old virgin flies of each sex were placed into a 30 cm 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  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 mated females the following day 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. sequencing total rna from each sample was isolated using trizol reagent following the manufacturer’s instructions. the concentration and integrity of total rna were determined using the bioanalyzer. total rna was then incubated with u dnase i at 37 °c for 30 min, and purified using the dynabeads® oligo following the manufacturer’s instructions. to ensure the accuracy of the results, we conducted two independent replicates for each sample, specifically immature, mature virgin and mated 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. 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, 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 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 and sequenced on an illumina hiseq™ 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, and then assembled into contigs using trinity 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 geea 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− gene ontology 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− and then the results were analyzed by gopipe based with the gene2go software. the kog 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, and then the mars model in the degseq package was used to calculate the expression abundance of each contig. a false discovery rate  <  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 <  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 verification selected deg data was validated by qrt-pcr using the sybr premix ex taq kit according to the manufacturer’s instructions with a real-time thermal cycler. trizol reagent was used to extract total rna from b. dorsalis at three different physiological stages, including the immature, mature virgin, and mated. at least insects were collected for each sample. the first strand cdna was obtained from 2 μg of total rna using m-mlv reverse transcriptase with the primer oligo-anchor r -3′). the primers used for qrt-pcr detection are listed in additional file 1: table s 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 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 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 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 24 h after mating. at least 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 were assembled into unique genes. the mean contig size was 883 bp with lengths ranging from to  bp. the contig size distribution revealed that more than half of the contigs were between and 500 bp in length, whereas  % were between and 3000 bp in length. gene ontology and clusters of orthologous group classification gene ontology analysis was performed for functional categorization of unigenes. within the biological processes, the main go terms were grouped into cellular and metabolic. as for the molecular function category, binding constituted the largest group, followed by catalytic activity. in the previously sequenced b. dorsalis transcriptome obtained from different developmental stages, the metabolic process group was much larger than the cellular process group. 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.. assignments of clusters of orthologous groups were performed to further evaluate the completeness of the transcriptome and the effectiveness of our annotation process. the annotated sequences were grouped into major functional classes. among them, the majority of the clusters were “signal transduction mechanisms”, “general function prediction only”, whereas “cell motility” and “nuclear structure” represented the smallest groups. our transcriptome will aid further research focused on signal transduction in tephritid files. transcriptional changes during female maturation to identify transcripts that may play a role in female maturation, differentially expressed gene sequences between the mature and immature libraries were identified with a total of degs. among them, transcripts were more abundant in the mature females, while were more abundant in immature females. our results are different from those revealed by the adult transcriptome of c. capitata obtained with a microarray, in which only transcripts displayed significant changes during female maturation. furthermore, in the transcriptome of c. capitata, the transcripts with greater abundance were found in the immature females instead of the mature individuals. these differences may be attributed to the use of different sequencing methods and the different sexual maturation time utilized by these two species. for example, b. dorsalis is sexually mature after 14 days, while c. capitata can mate after only 4 days. therefore, more genes are possibly implicated in the sexual maturation in b. dorsalis.fig. comparisons of differentially expressed genes between the mature virgin and immature samples, as well as between the mated and mature virgin samples. the number of up-regulated genes is represented in the red column, while down-regulated genes are shown in green analysis of enriched biological process go terms with significant differences during female maturation was performed. we focused on the physiological processes that appear to play essential roles in maturation, including reproductive, signaling, developmental, metabolic, immune processes as well as response to stimulus. a table of the identified genes is provided in the supplementary data. results revealed that differentially expressed genes in the metabolic and developmental processes are a much larger group than in other processes, which is similar to the results in c. capitata that indicated the most enriched terms in mature females were categorized in the metabolic process group. these data indicate that a variety of metabolic activities and developmental affairs usually occur during oogenesis and ovary maturation. additionally, enrichment analysis of the kegg pathways for the differentially expressed genes during maturation was also performed. results showed that the number of differentially expressed genes in the metabolism pathways are the largest groups, including carbohydrate metabolism, energy metabolism, nucleotide metabolism, lipid metabolism, metabolism of cofactors and vitamins, amino acid metabolism, etc.. this is similar with the result of go enrichment analysis.fig. go classification of the degs to demonstrate changes in transcript abundance with maturation and mating. goeast was used for go term enrichment analysis, and 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 <  was considered to be significant enriched. the number of enriched transcripts for go annotated sequences was revealed in the various samples of different biological categories. a mature virgin females compared to immature females, and b mated females compared to mature virgin females 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 transcripts that displayed significant transcriptional changes between virgin mature and mated females 24 h post-mating, were more abundant in mated females, while only were more abundant in virgin mature females. this is similar with the result found in c. capitata and d. melanogaster. only and 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, 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. 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 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. a table of these genes is provided in the supplementary data. 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. furthermore, enrichment analysis of kegg pathways for the mating responsive genes was also performed. 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 from deg analysis of the samples from immature vs. mature virgin individuals, we identified seven genes that reportedly participate in oogenesis and ovary maturation in insects. interestingly, two genes that were previously found to play important roles in wnt signaling pathways and an additional two involved in sex determination were also identified in this deg library. the subsequent qrt-pcr validated this result. all of these genes were significantly highly expressed at the sexually mature and post-mating stage, which is consistent with our transcriptome data, indicating that they are involved in physiological processes of mature females as is the case in other insects. additionally, we analyzed expression of these genes in several tissues to elucidate their roles and found that most are highly enriched in the ovary, including hu-li tai shao, oskar, mago nashi, yolkless, disheveled, axin, sry and transformer-. the expression level of oskar in the ovary is approximately fold higher than that in the fat body and fold in the midgut. notably, three vitellogenin genes can be detected in the head, thorax, ovary, fat body, while they are more abundantly expressed in fat body. these results are explained below by analyzing their roles determined for other insects.fig. qrt-pcr confirmation of the differentially expressed genes between mature virgin females and the immature individuals. total rnas from b. dorsalis at three different physiological stages were isolated, including the immature, virgin and mated. relative transcript levels are calculated by real-time pcr using the 16 s rrna gene as the standard control. different letters represent significant differences among samples. three independent biological replicates were performed for each sample fig. differentially expressed genes between mature virgin females and immature individuals in various tissues. qrt-pcr was used to analyze mrna levels of differentially expressed genes, including those found in the head, thorax, ovary, fat body and midgut. relative transcript levels are calculated by real-time pcr using the 16 s rrna gene as the standard control. various letters indicate significant differences in the expression level. three biological replicates were performed 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/y 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/y 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. 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 vitellogenin and vitellogenin 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 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 in our transcriptome, both of which are core components of wnt signaling pathways. 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. notably, two genes involved in sex determination were also identified, including sex-determining region y protein and transformer-. 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- was also reported to be an essential switch for the sex determination cascade in dipteran insect species. knockdown of tra- 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 upon mating, significant expression changes occurred in genes involved in the immune or stress response. the immune response stimulated by mating might arise from sexually antagonistic interaction between the sexes. several up-regulated genes from the deg analysis showed significant increases in transcript abundance in females, including two prrs and six antimicrobial peptides. to investigate the expression pattern of these genes, the expression levels of amps in different tissues were analyzed by qrt-pcr. the results showed that all were highly enriched in the fat body, which is usually an important organ involved in the immune response in insects. the expression pattern of these genes indicates that they may participate in the immune response after mating in b. dorsalis.fig. expression pattern of the immune-related genes between mated females and the mature virgin individuals in different tissues, including the head, thorax, ovary, fat body and midgut. relative transcript levels of six antimicrobial peptide are calculated by real-time pcr using the 16 s rrna gene as the standard control. different letters indicate significant differences among samples and three independent biological replicates were performed for each sample to further validate, the changes in expression at different time points after mating were analyzed in both the fat body and head. all genes except attacin-c-like and phormicin-like were significantly up-regulated 24 h after mating in the fat body approximately 2– fold. this is partly supported by the previous studies in d. melanogaster, which showed that the immune response is activated in mated females by mating and associated sex peptides transfer. upon mating, potent signal molecules from the males enter the females and the behavior and physiology of the mated females shifted from reproductive receptivity to fecundity. the immune response was induced in the females and it might arise from sexually antagonistic interaction between the sexes. however, this is different from the results of the previous report on c. capitata, which revealed no significant changes of transcript abundance were detected in immune-related genes after mating, except for a large reduction in defensin in the female abdomen. this data suggests that post mating changes differ among species. additionally, the expression level of most amps did not significantly change during the time course in the female head 24 h post-mating, except for sapecin and cecropin- both of which significantly increased in transcript abundance, by and fold, respectively. it is clear that the fat body is the main tissue where the immune response takes place and not the head.fig. the effect of mating on the expression level of immune-related genes in the fat body and head of b. dorsalis females. rna was extracted from the head or fat body of mature vigin females as well as mated females at 24 h after mating. relative transcript levels of six antimicrobial peptide are calculated by real-time pcr using the 16 s rrna gene as the standard control. various letters represent significant differences of the expression level of genes. three biological replicates were performed  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. availability of supporting data this transcriptome shotgun assembly project has been deposited at ddbj/embl/genbank under the accession geea other data sets supporting the results of this article are included in the additional files. additional files additional file 1: table s primers used for qrt-pcr. additional file 2: table s summary of the solexa sequencing statistics of b. dorsalis adult transcriptome. additional file 3: table s summary of the b. dorsalis adult transcriptome. additional file 4: figure s assembled contig length distribution of the transcriptome. the x-axis shows contig size and the y-axis indicates the number of contigs for each given size. additional file 5: figure s gene ontology analyses of transcriptome sequences predicted to participate in: biological processes and molecular functions. data are presented as level go categorization for molecular function and level go categorization for biological process. classified gene objects are showed as percentages of the total number according to the go assignments. additional file 6: figure s histogram presentation of clusters of orthologous groups classification. a total of predicted unigenes were classified in cog categories. additional file 7: file s genes differentially expressed between mature virgin females and immature individuals, and between mated females and mature virgins. reads, rpkm and the fold values for each comparison group is shown. additional file 8: table s oriental fruit fly assembled sequences of degs with best-hit matches to dipteran genes involved in the sexual maturation. additional file 9: file s enrichment analysis of kegg pathways for the differentially expressed genes. additional file 10: table s oriental fruit fly assembled sequences of degs with best-hit matches to dipteran genes involved in the immune system process or response to stimulus. competing interests the authors declare that they have no competing interests. authors’ contributions the study was conceived and designed by wwz. dyl and fyw prepared the oriental fruit fly samples. fyw conducted the rna extraction. wwz performed the qrt-pcr. wwz, jlw and hyz analyzed the transcriptomic data. wwz prepared the initial draft of the manuscript and hyz modified the final version. all authors have read and approved the final version of the manuscript. 
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. in total, differentially expressed genes were detected between mature virgin and immature bactrocera dorsalis libraries, whereas only degs were identified between flies that had mated or were mature virgins. 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. 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. 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. this study supplied female reproductive genes involved in sexual maturation and the post-mating response in b. dorsalis, based on rna-seq. we sequenced unique genes with an average length of 883 bp. our data will facilitate molecular research related to reproduction and provide abundant target genes for effective control of this agricultural pest.
 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 unique genes with an average length of 883 bp. in total, differentially expressed genes were detected between mature virgin and immature bactrocera dorsalis libraries, whereas only 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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords rnaseqtranscriptometephritidmaturitypost-matingreproductiongene expressionhttp://dx.doi.org/ /501100001809national natural science foundation of chinaid0esnai8701zheng weiwei http://dx.doi.org/ /501100001809national natural science foundation of chinaid0e1tai8702zheng weiwei the fundamental research funds for the central universitiesid0ec1ai8703zheng weiwei china agriculture research systemid0eiabi8704zhang hongyu issue-copyright-statement© the author 2016 