Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis

Comparative analysis of testicular fusion in Spodoptera litura (cutworm) and Bombyx mori (silkworm): Histological and transcriptomic insights

Spodoptera litura commonly known as the cutworm, is among the most destructive lepidopteran pests affecting over 120 plants species. The powerful destructive nature of this lepidopteran is attributable to its high reproductive capacity. The testicular fusion that occurs during metamorphosis from larvae to pupa in S.litura positively influences the reproductive success of the offspring. In contrast, Bombyx mori, the silkworm, retains separate testes throughout its life and does not undergo this fusion process. Microscopic examination reveals that during testicular fusion in S.litura, the peritoneal sheath becomes thinner and more translucent, whereas in B.mori, the analogous region thickens. The outer basement membrane in S.litura exhibits fractures, discontinuity, and uneven thickness accompanied by a significant presence of cellular secretions, large cell size, increased vesicles, liquid droplets, and a proliferation of rough endoplasmic reticulum and mitochondria. In contrast, the testicular peritoneal sheath of B.mori at comparable developmental stage exhibits minimal change. Comparative transcriptomic analysis of the testicular peritoneal sheath reveals a substantial difference in gene expression between the two species. The disparity in differential expressed genes (DEGs) is linked to an enrichment of numerous transcription factors, intracellular signaling pathways involving Ca2+ and GTPase, as well as intracellular protein transport and signaling pathways. Meanwhile, structural proteins including actin, chitin-binding proteins, membrane protein fractions, cell adhesion, extracellular matrix proteins are predominantly identified. Moreover, the study highlights the enrichment of endopeptidases, serine proteases, proteolytic enzymes and matrix metalloproteins, which may play a role in the degradation of the outer membrane. Five transcription factors-Slforkhead, Slproline, Slcyclic, Slsilk, and SlD-ETS were identified, and their expression pattern were confirmed by qRT-PCR. they are candidates for participating in the regulation of testicular fusion. Our findings underscore significant morphological and trancriptomic variation in the testicular peritoneal sheath of S.litura compared to the silkworm, with substantial changes at the transcriptomic level coinciding with testicular fusion. The research provides valuable clues for understanding the complex mechanisms underlying this unique phenomenon in insects.

BdTTLL3B-mediated polyglycylation is involved in the spermatogenesis in Bactrocera dorsalis

Polyglycylation is a post-translational modification that generates glycine side chains in the C-terminal domains of both α- and β-tubulins. To date, the patterns and significance of polyglycylation across insect species remain largely unknown. The TTLL3B was thought to be a polyglycylase and be essential for polyglycylation in dipteran insects. In this study, the TTLL3B of Bactrocera dorsalis (BdTTLL3B) was identified and characterized. The BdTTLL3B expressed remarkably higher in adult males, especially in testes. The spatio-temporal patterns of polyglycylation were consistent with that of BdTTLL3B. Along with spermatogenesis, the intensity of polyglycylation was enhanced steadily and concentrated in elongated flagella. The expression of recombinant BdTTLL3B in Hela cells, which are genetically deficient in polyglycylation, catalyzed intracellular polyglycylation, validating the identity of BdTTLL3B as a polyglycylase. Knockout of BdTTLL3B significantly suppressed polyglycylation in testes and impaired male fertility, probably due to abnormal morphology of mitochondrial derivatives and over-accumulation of paracrystalline. Taken together, these findings indicated that the BdTTLL3B-mediated polyglycylation is involved in the spermatogenesis and play an important role in fertility of adult B. dorsalis. Therefore, the BdTTLL3B can be considered as a candidate target gene for the management of B. dorsalis, such as developing gene silencing/knockout-based sterile insect technology (SIT).

CRISPR/Cas9 mediated validation of spermatogenesis-related gene, tssk2 as a component of genetic pest management of fall armyworm, Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae)

Invasive insect pests, currently, pose a serious economic threat to several staple crops all over the world, one such being the fall armyworm, Spodoptera frugiperda. It was first observed in Africa since 2016, outside of its natural habitat in the Americas. Subsequently, it invaded several countries in South and South East Asia and also very recently in Australia. In all the newly invaded regions, maize is the principal crop attacked causing a serious economic concern to the poor farmers, particularly in the developing countries. Owing to the innate genetic ability, it defies many of the management options that include insecticides, Bt transgenics, and so forth. This is due to its high mobility, polyphagy and ability for quick development of resistance to several classes of insecticides. At this critical juncture, CRISPR/Cas9 mediated genome editing has shown a lot of promise in developing a novel area-wide pest management strategy called precision-guided sterile insect technique (pgSIT). pgSIT was initially demonstrated in Drosophila melanogaster which holds a greater promise for the environmentally friendly management of several globally significant agricultural pests such as S. frugiperda. Therefore, before developing both sgRNA and Cas9 transgenic lines, we have validated the target gene such as tssk2 through a non-transgenic approach by microinjecting ribo nucleo protein complex (Cas9 protein and tssk2 sgRNA) into G0 eggs of S. frugiperda. In the current investigation, we have obtained five edited males with distinct mutations which were further used for crossing studies to ascertain the effect of tssk2 editing affecting egg hatchability.

Identification of testis development-related genes by combining Iso-Seq and RNA-Seq in Zeugodacus tau

Introduction: Zeugodacus tau (Walker) is an invasive pest. An effective method to control this pest is the sterile insect technique (SIT). To better apply this technique, it is necessary to understand testis development progression. Methods: Differentially expressed genes (DEGs) during testis development were analyzed by PacBio Iso-Seq and RNA-seq. Results: RNA-Seq library of Z. tau testes on day 1, 6, and 11 post eclosion were constructed. We identified 755 and 865 differentially expressed genes in the comparisons of T6 (testes on day 6) vs. T1 and T11 vs. T1, respectively. The KEGG pathway analysis showed that the DEGs were significantly enriched in retinol metabolism, vitamin B6 metabolism, and ascorbate and aldarate metabolism pathways. Knockdown of retinol dehydrogenase 12-like (rdh12-like), pyridoxal kinase (pdxk) and regucalcin (rgn), the representative gene in each of the above 3 pathways, reduced the hatching rate of Z. tau offspring. In addition, we identified 107 Drosophila spermatogenesis-related orthologous genes in Z. tau, of which innexin 2 (inx2) exhibited significantly up-regulated expression throughout testis development, and the knockdown of this gene reduced offspring hatching rate. Discussion: Our data indicated that rdh12-like, pdxk, rgn, and inx2 genes were related to testis development, and they were conserved in tephritid species. These results suggested that this gene might have the same function in tephritid. The findings provide an insight into testis development and spermatogenesis in tephritid species.

Candidate target genes of the male-specific expressed Doublesex in the termite Reticulitermes speratus

Eusocial insects such as termites, ants, bees, and wasps exhibit a reproductive division of labor. The developmental regulation of reproductive organ (ovaries and testes) is crucial for distinguishing between reproductive and sterile castes. The development of reproductive organ in insects is regulated by sex-determination pathways. The sex determination gene Doublesex (Dsx), encoding transcription factors, plays an important role in this pathway. Therefore, clarifying the function of Dsx in the developmental regulation of sexual traits is important to understand the social evolution of eusocial insects. However, no studies have reported the function of Dsx in hemimetabolous eusocial group termites. In this study, we searched for binding sites and candidate target genes of Dsx in species with available genome information as the first step in clarifying the function of Dsx in termites. First, we focused on the Reticulitermes speratus genome and identified 101 candidate target genes of Dsx. Using a similar method, we obtained 112, 39, and 76 candidate Dsx target genes in Reticulitermes lucifugus, Coptotermes formosanus, and Macrotermes natalensis, respectively. Second, we compared the candidate Dsx target genes between species and identified 37 common genes between R. speratus and R. lucifugus. These included several genes probably involved in spermatogenesis and longevity. However, only a few common target genes were identified between R. speratus and the other two species. Finally, Dsx dsRNA injection resulted in the differential expression of several target genes, including piwi-like protein and B-box type zinc finger protein ncl-1 in R. speratus. These results provide valuable resource data for future functional analyses of Dsx in termites.

RNAi-mediated knockdown of α-Spectrin depresses reproductive performance in female Bactrocera dorsalis

The female reproductive potential plays a crucial role in reproduction, population dynamics and population maintenance. However, the function of endogenous genes in undifferentiated germ cells has been largely unknown in Bactrocera dorsalis. In this study, the conservative analysis showed that α-Spectrin shared a similarity in B. dorsalis and other dipteral flies. Further, the differential expression of α-Spectrin was examined in B. dorsalis by RT-qPCR, and the expression pattern of α-Spectrin protein was identified in female adult ovaries by using immunostaining. During the development of ovary, the change on the number of undifferentiated germ cells was also characterized and analyzed. To understand the function of α-Spectrin in B. dorsalis ovary, the RNAi-based knockdown was conducted, and the RNAi efficiency was examined by RT-qPCR, western blot and bioassay. The results revealed that the α-Spectrin dsRNA could strikingly decrease the expression level of α-Spectrin in ovaries and diminish oviposition and ovary size as a consequence of downregulation of α-Spectrin. Overall, our study facilitates reproductive research on the function of conservative genes in B. dorsalis ovary, which may provide a new insight into seeking novel target genes for pest management control.

Diallyl Trisulfide Causes Male Infertility with Oligoasthenoteratospermia in Sitotroga cerealella through the Ubiquitin-Proteasome Pathway

Essential oils extracted from plant sources along with their biologically active components may have negative effects on insects. Diallyl trisulfide (DAT) is an active component of garlic essential oil, and it exhibits multi-targeted activity against many organisms. Previously we reported that DAT induces male infertility and leads to apyrene and eupyrene sperm dysfunction in Sitotroga cerealella. In this study, we conducted an analysis of testis-specific RNA-Seq data and identified 449 downregulated genes and 60 upregulated genes in the DAT group compared to the control group. The downregulated genes were significantly enriched in the ubiquitin-proteasome pathway. Furthermore, DAT caused a significant reduction in mRNA expression of proteasome regulatory subunit particles required for ATP-dependent degradation of ubiquitinated proteins as well as decreased the expression profile of proteasome core particles, including β1, β2, and β5. Sperm physiological analysis showed that DAT decreased the chymotrypsin-like activity of the 20S proteasome and formed aggresomes in spermatozoa. Overall, our findings suggest that DAT impairs the testis proteasome, ultimately causing male infertility characterized by oligoasthenoteratospermia due to disruption in sperm proteasome assembly in S. cerealella.

Comparative multi-tissue analyses identify testis-specific serine/threonine protein kinase (TSSK) genes involved in male fertility in the melon fly Zeugodacus cucurbitae

BACKGROUND

Zeugodacus cucurbitae is an agricultural pest species with robust reproductive capabilities capable of causing extensive damage. The advent of novel male fertility-related pest control strategies has been an area of active entomological research focused on the sterile insect technique (SIT) strategy.

RESULTS

RNA-sequencing analyses were conducted using 16 tissue samples from adult male Z. cucurbitae, leading to the identification of 5338 genes that were differentially expressed between the testes and three other analyzed tissue types. Of these genes, 808 exhibited high levels of testis expression. A quantitative polymerase chain reaction (qPCR) approach was used to validate the expression of ten of these genes selected at random, including ZcTSSK1 and ZcTSSK3, which are testis-specific serine/threonine protein kinase (TSSK) genes. Evaluation via a loss-of-function-based knockdown assay showed that the down-regulation of either of these two genes in males was associated with significantly decreased egg hatching rates. In situ hybridization analyses revealed the expression of both of these transcripts in the transformation zone, and significant decreases in Z. cucurbitae sperm numbers were observed following double-stranded RNA treatment. Together, these results suggested that inhibiting ZcTSSK1 and ZcTSSK3 expression was sufficient to alter male fertility in Z. cucurbitae.

CONCLUSION

These transcriptional sequencing results provide a foundation for further efforts to clarify the regulators of Z. cucurbitae male fertility. These preliminary analyses of the functions of ZcTSSK family genes as regulators of spermatogenesis underscore their importance in the processes integral to male fecundity and their potential as targets for pest control efforts centered on the genetic manipulation of males. © 2023 Society of Chemical Industry.

Identification, expression profiles and involvement in insecticides tolerance and detoxification of carboxylesterase genes in Bactrocera dorsalis

Carboxylesterases (CarEs) are a multifunctional superfamily of enzymes and play an important role in detoxification of various insecticides in insects. The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive agricultural pests and has developed different degrees of resistance to organophosphates in field. However, the involvement of BdCarEs in tolerance or resistance to other alternative insecticides are still unclear. In the present study, 33 BdCarEs genes were identified based on the genome database of B. dorsalis. Phylogenetic analysis demonstrated that they were classified into nine clades, with abundance of α-esterases. Meanwhile, the sequence characterization and the chromosome distribution were also analyzed. The spatiotemporal expression analysis of BdCarEs genes suggested that the diversity of potential function in different physiological processes. With the exception of BdCarE21, all BdCarEs genes responded to at least one insecticide exposure, and BdCarE20 was found to be up-regulated after exposure to all five tested insecticides individually. Eight BdCarEs genes were overexpressed in MR strain when compared to that in SS strain. Subsequently, knockdown the expression of representative BdCarEs genes significantly increased the susceptibility of the oriental fruit fly to corresponding insecticides, which indicated that the tested BdCarEs genes contributed to one or multiple insecticide detoxification. These findings provide valuable insights into the potential role in respond to tolerance or resistance to insecticides with different mode of action, and will facilitate development of efficiency management strategy for B. dorsalis.

Chromosome-level genome assembly reveals potential epigenetic mechanisms of the thermal tolerance in the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis (Hendel), has very strong ecological adaptability and phenotypic plasticity. Here, the genome of B. dorsalis was assembled into 549.45 Mb sequences with a contig N50 length of 12.81 Mb. Among, 95.67 % assembled genome sequences were anchored on six chromosomes with an N50 length of 94.63 Mb. According to the basic characteristics of the sex chromosomes of Tephritidae, the X chromosome of B. dorsalis was identified. Significant gene expansions were detected in several important gene families related to adaptability. In particular, we annotated 50 histone modification enzymes (HMEs) in this genome. A comparative transcriptome analysis indicated that 12 HME genes were differentially expressed in two thermo-tolerant strains (heat and cold). Interestingly, four and seven of the 12 HME genes responded to heat shock or cold hardening, respectively. These evidences suggested that the histone modification as an epigenetic modification may be involved in the thermal tolerance of B. dorsalis, but with different regulation mechanisms in thermal acclimation and hardening. The high quality genome of B. dorsalis provides an invaluable resource for further functional genomic study. Moreover, comparative genomic analysis will shed insights on revealing the mechanisms of adaptive evolution in this fly.

Knockdown of a Nicotinic Acetylcholine Receptor Subunit Gene Bdorβ1 Decreases Susceptibility to Oxa-Bridged trans- instead of cis-Nitromethylene Neonicotinoid Insecticides in Bactrocera dorsalis

Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that mediate the fast action of acetylcholine in synaptic cholinergic transmissions. Insect nAChRs are the target of several classes of insecticides. Here, the full-length cDNA encoding a nAChR beta1 subunit (Bdorβ1) was identified and characterized from a destructive pest, Bactrocera dorsalis. The amino acid sequence of Bdorβ1 shows high identities to other insect nAChRs β1 subunits. Double injection of dsBdorβ1 reduced the expression of Bdorβ1 and in turn significantly decreased susceptibility to oxa-bridged trans- instead of cis-nitromethylene neonicotinoids. Our results support the involvement of Bdorβ1 in the susceptibility of B. dorsalis to oxa-bridged trans- instead of cis-nitromethylene neonicotinoids and imply that these two classes of neonicotinoids might be acting at different nAChR subtypes.

Comparative analysis of testis transcriptome between a genetic male sterile line (GMS) and its wild-type 898WB in silkworm, Bombyx mori

The silkworm, Bombyx mori, is an important model organism of lepidopteran insects, and its testis is a main male reproductive organ and spermatogenesis place. Studying the testis helps to understand the mechanisms of genetic sterility of lepidopteran insects and to achieve sterile insect technique (SIT) for pest control. Herein, we performed a comparative transcriptome analysis of testes between three biological replicates of the GMS mutant and wild strain 898WB, respectively. In total, 1872 up-regulated genes and 1823 down-regulated genes were identified in the testis of the GMS mutant. Several genes contribute significantly to spermatogenesis and testis development, such as "serine/threonine protein kinase", "organic cation transporter protein", "tyrosine protein kinase", "lncRNAs" and "immune-associated genes". The KEGG pathway analysis shows that the DEGs were annotated to 123 pathways, and 10 pathways were significantly enriched, such as "metabolic pathway", "biosynthesis of amino acids", and "phagosome-lysosome pathway", which are associated with testis development and spermatogenesis. The results of the qPCR expression were consistent with the RNA-seq data, which shows that the RNA-seq results were accurate. The DEGs of the testes between GMS mutant and 898WB were screened by RNA-Seq technology, which provides a reliable reference to understand the molecule mechanism of male sterility of the GMS mutant.

Testes morphology and the identification of transcripts of the hormonal pathways of the velvetbean caterpillar Anticarsia gemmatalis Hübner, 1818 (Lepidoptera: Erebidae)

Anticarsia gemmatalis is one of the main defoliating pests of soybeans in Brazil. In the current study, we characterized the histomorphology of the testes and the spermatogenesis process in A. gemmatalis. We also identified transcripts involved in the biosynthesis, metabolism, and signaling of juvenile and ecdysteroid hormones, in order to provide information about potential mechanisms of regulation of hormonal pathways in this species. Our analyses revealed that the A. gemmatalis larvae have a pair of kidney-shaped testicles. These are divided into four testicular follicles, where there are germ cell cysts at different stages of development. In the pupal stage, the testicles are fused, so adults have a single spherical testis, with a variable number of follicles. The A. gemmatalis has centripetal spermatogenesis and exhibits spermatic dimorphism. We identified 31 transcripts that encode proteins involved in juvenile hormone and ecdysteroid pathways, such as mevalonate kinase, CYP14A1, ecdysone receptor, among others. Our results on the morphology of the testes and spermatogenesis process, as well as identification of the genes involved in hormonal pathways in A. gemmatalis, provide important data for understanding the biology of this agricultural pest, which can be used as a basis for further research in other economically important lepidopterans.

Molecular Characterization and Transcriptional Expression Analysis of ABC Transporter H Subfamily Genes in the Oriental Fruit Fly

The oriental fruit fly, Bactrocera dorsalis Hendel (Diptera: Tephretidae), is a serious pest of fruits and vegetables and has developed high levels of insecticide resistance. ATP-binding cassette transporter genes (ABC transporters) are involved in mediating the energy-driven transport of many substances across membranes and are closely associated with development and insecticide detoxification. In this study, three ABC transporters in the H subfamily were identified, and the possible roles of these genes in B. dorsalis are discussed. Bioinformatics analysis revealed that those genes are conserved, typical of half-transporters. The expression profiles of BdABCH genes (BdABCHs) in the developmental stages, tissues, and following insecticide exposure, extreme temperature, warm- and cold-acclimated strain, starvation, and desiccation stress were determined by quantitative real-time PCR. Expression of BdABCHs can be detected in various tissues and in different developmental stages. They were most highly expressed in the hindgut and in newly emerged adults. The mRNA levels of BdABCHs in males (including most tissues and body segments) were higher than in females. The expression of BdABCH1 was significantly upregulated 3.8-fold in the cold-acclimated strain, and was significantly upregulated by 1.9-, 3.8- and 4.1-fold in the 0°C, starvation, and desiccation treatments, respectively. Treatment with malathion and avermectin at LD20 and LD30 concentrations produced no obvious changes in the levels of BdABCHs. BdABCHs may be involved in the transport of related hormones during eclosion, as well as water and inorganic salts. BdABCH1 also demonstrated that it is related to the ability to cope with adverse environments.

The adipokinetic hormone signaling system regulates the sensitivity of Bactrocera dorsalis to malathion

The neuropeptide adipokinetic hormone (AKH) binds to the AKH receptor (AKHR) to regulate carbohydrate and lipid metabolism. It also participates in the insect anti-stress response. We used RT-qPCR to detect the expression levels of 39 neuropeptides in malathion-susceptible (MS) and malathion-resistant (MR) strains of Bactrocera dorsalis. AKH and AKHR were highly expressed in the MR strain. Using a malathion bioassay and RNA interference (RNAi), we demonstrated that AKHR is involved in the susceptibility of B. dorsalis to malathion. We found significantly reduced expression of two detoxification enzyme genes (glutathione-S-transferase, GST and α-esterase, CarE) after AKHR RNAi. Based on our previous data, GSTd10 and CarE6 participate the direct metabolism of malathion in this fly, which is also verified by a malathion metabolism assay by HPLC using the crude enzymes in the current study. These results suggest that AKHR plays an important role in affecting malathion susceptibility via detoxification enzyme genes.

Genome-wide identification of long non-coding RNAs (lncRNAs) associated with malathion resistance in Bactrocera dorsalis

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in the regulation of biological processes and have been identified in many species including insects. However, the association between lncRNAs and pesticide resistance in insect species such as Bactrocera dorsalis is unknown.

RESULTS

RNA-seq was performed on malathion resistant (MR1) and susceptible (MS) strains of B. dorsalis and a total of 6171 lncRNAs transcripts were identified. These included 3728 lincRNAs, 653 antisense lncRNAs, 1402 intronic lncRNAs, and 388 sense lncRNAs. A total of 40 and 52 upregulated lncRNAs were found in females and males of the MR1 strain compared to 54 and 49 in the same sexes of the MS strain, respectively. Twenty-seven of these lncRNAs showed the same trend of expression in both females and males in the MR1 strain, in which 15 lncRNAs were upregulated and 12 were downregulated. RT-qPCR results indicated that the differentially expressed lncRNAs were associated with malathion resistance. The lnc15010.10 and lnc3774.2 were highly expressed in the cuticle of the MR1 strain, indicating that these two lncRNAs may be related to malathion resistance. RNAi of lnc3774.2 and a bioassay showed that malathion resistance was possibly influenced by changes in the B. dorsalis cuticle.

CONCLUSION

LncRNAs of B. dorsalis potentially related to the malathion resistance were identified. Two lncRNAs appear to influence malathion resistance via modulating the structure, or components, of the cuticle. © 2021 Society of Chemical Industry.

Ovary-Specific Transcriptome and Essential Role of Nanos in Ovary Development in the Oriental Fruit Fly (Diptera: Tephritidae)

We used transcriptome analysis to research ovary development in Bactrocera dorsalis (Hendel). The ovary transcriptome of B. dorsalis yielded 66,463,710 clean reads that were assembled into 23,822 unigenes. After aligning to the Nr database in NCBI, 15,473 (64.95%) of the unigenes were matched to identified proteins. As determined by BLAST search, 11,043 (46.36%), 6,102 (25.61%), and 12,603 (52.90%) unigenes were each allocated to clusters via gene ontology, orthologous groups, and SwissProt, respectively. The Kyoto encyclopedia database of genes and genomes (KEGG) was further used to annotate these sequences, and 11,068 unigenes were mapped to 255 known pathways. Afterward, the genes that were possibly involved in oogenesis and ovary development were obtained from the transcriptome data and analyzed. Interestingly, seven ovary-specific genes were identified, including a Nanos gene that is involved in maintaining the primordial germ cells in many insects. Therefore, we further focused on the function of the BdNanos gene, and the gene was injected into B. dorsalis. As expected, the knocking down of Nanos gene expression led to significant inhibition of ovary development, suggesting an important role of this gene in the reproductive process of B. dorsalis. In summary, the present study provides an important reference for identifying the molecular mechanisms of oogenesis and ovary development in B. dorsalis. The BdNanos gene is crucial for ovary development in B. dorsalis and is therefore a potential new pest control target.

Decoding the Reproductive System of the Olive Fruit Fly, Bactrocera oleae

In most diploid organisms, mating is a prerequisite for reproduction and, thus, critical to the maintenance of their population and the perpetuation of the species. Besides the importance of understanding the fundamentals of reproduction, targeting the reproductive success of a pest insect is also a promising method for its control, as a possible manipulation of the reproductive system could affect its destructive activity. Here, we used an integrated approach for the elucidation of the reproductive system and mating procedures of the olive fruit fly, Bactrocera oleae. Initially, we performed a RNAseq analysis in reproductive tissues of virgin and mated insects. A comparison of the transcriptomes resulted in the identification of genes that are differentially expressed after mating. Functional annotation of the genes showed an alteration in the metabolic, catalytic, and cellular processes after mating. Moreover, a functional analysis through RNAi silencing of two differentially expressed genes, yellow-g and troponin C, resulted in a significantly reduced oviposition rate. This study provided a foundation for future investigations into the olive fruit fly's reproductive biology to the development of new exploitable tools for its control.

Genome-wide gene expression profiling of the melon fly, Zeugodacus cucurbitae, during thirteen life stages

The melon fly, Zeugodacus cucurbitae (Coquillett), is an important destructive pest worldwide. Functional studies of the genes associated with development and reproduction during different life stages are limited in Z. cucurbitae. There have yet to be comprehensive transcriptomic resources for genetic and functional genomic studies to identify the molecular mechanisms related to its development and reproduction. In this study, we comprehensively sequenced the transcriptomes of four different developmental stages: egg, larva, pupa, and adults. Using the Illumina RNA-Seq technology, we constructed 52 libraries from 13 stages with four biological replicates in each and generated 435.61 Gb clean reads. We comprehensively characterized the transcriptomes with high-coverage mapping to the reference genome. A total of 13,760 genes were mapped to the reference genome, and another 4481 genes were characterized as new genes. Finally, 14,931 genes (81.85%) were functionally annotated against six annotation databases. This study provides the first comprehensive transcriptome data of all developmental stages of Z. cucurbitae, and will serve as a valuable resource for future genetic and functional studies.

Comparative proteomic analysis of sex-biased proteins in ovary and testis at different stages of Spodoptera litura

Sex-biased protein is thought to be able to drive the phenotypic differences in males and females in insects. In this study, 1385 and 1727 proteins were identified as differentially accumulated proteins (DAPs) by comparing the protein abundances at pupae stage with those at adult stage in ovary and testis of S.litura, respectively. And among which, 548 DAPs were showed to be expressed in both ovary and testis, and 837 and 1179 proteins were considered as ovary-specific and testis-specific DAPs, respectively. To further identify DAPs related to gonad development and sex dimorphism, a total of 320 DAPs were selected and defined as "proteins of specific interest" based on several selecting criteria. Sex dimorphism is a complex and dynamic developmental progress, and these identified DAPs were suggested to be involved in multiple functions such as organonitrogen compound catabolic process, glycosylation, proteasome, N-Glycan biosynthesis and other reproduction-related processes. Overall, our results highlighted these sexual-biased, gonad development related and sexual dimorphism related DAPs, and their abundance variations along with development were also examined, which could provide important information for their functional analysis in reproduction and potential biomarkers for developing useful strategies against S. litura and other orthologous pests. BIOLOGICAL SIGNIFICANCE: Sex dimorphism entails the differentiation of two sexual functions, resulting in sexually phenotypic differences and leading to the development of female and male morphologies and behaviors. However, sex dimorphism related proteins remain to be identified in many non-model insects. In this study, iTRAQ-based proteomic analysis was applied to examine the variations of protein abundances at pupae stage and adult stage in ovary and testis of S.litura, respectively. Reproduction and sex dimorphism related proteins were further identified as "proteins of specific interest". These identified candidate proteins provided valuable information for their further functional analysis in reproduction and could serve as potential biomarkers for developing useful strategies against S. litura and other orthologous pests.

De Novo Transcriptome Analysis Reveals Abundant Gonad-specific Genes in the Ovary and Testis of Henosepilachna vigintioctopunctata

Henosepilachna vigintioctopunctata (Coleoptera: Coccinellidae) is a major pest affecting Solanaceae plants in Asian countries. In this study, we sequenced the ovary and testis transcriptomes of H. vigintioctopunctata to identify gonad-related genes. Comparison of the unigene sequences in ovary and testis libraries identified 1,421 and 5,315 ovary- and testis-specific genes, respectively. Among the ovary-specific genes, we selected the RC2-like and PSHS-like genes to investigate the effects of gene silencing on the mortality, percentage infertility, pre-oviposition period, fecundity, daily number of eggs laid, and hatching rate in female adults. Although the percentage mortality and infertility of females did not differ significantly among dsRNA treatments, fecundity was significantly reduced in the dsRC2-like and dsPSHS-like treatment groups. Moreover, the pre-oviposition period was markedly prolonged in response to dsPSHS-like treatment. This is the first reported RNA sequencing of H. vigintioctopunctata. The transcriptome sequences and gene expression profiles of the ovary and testis libraries will provide useful information for the identification of gonad-related genes in H. vigintioctopunctata and facilitate further research on the reproductive biology of this species. Moreover, the gonad-specific genes identified may represent candidate target genes for inhibiting the population growth of H. vigintioctopunctata.

Reduced Glutamine Synthetase Activity Alters the Fecundity of Female Bactrocera dorsalis (Hendel)

Glutamine synthetase (GS) is a key enzyme in glutamine synthesis and is associated with multiple physiological processes in insects, such as embryonic development, heat shock response, and fecundity regulation. However, little is known about the influence of GS on female fecundity in the oriental fruit fly, Bactrocera dorsalis. Based on the cloning of BdGSs, mitochondrial BdGSm and cytoplasmic BdGSc, we determined their expressions in the tissues of adult B. dorsalis. BdGSm was highly expressed in the fat body, while BdGSc was highly expressed in the head and midgut. Gene silencing by RNA interference against two BdGSs isoforms suppressed target gene expression at the transcriptional level, leading to a reduced ovarian size and lower egg production. The specific inhibitor L-methionine S-sulfoximine suppressed enzyme activity, but only the gene expression of BdGSm was suppressed. A similar phenotype of delayed ovarian development occurred in the inhibitor bioassay. Significantly lower expression of vitellogenin and vitellogenin receptor was observed when GS enzyme activity was suppressed. These data illustrate the effects of two GS genes on adult fecundity by regulating vitellogenin synthesis in different ways.

RNAi-Mediated Knockdown of Tssk1 and Tektin1 Genes Impair Male Fertility in Bactrocera dorsalis

The genetic-based sterile insect technique (SIT) is an effective and environmentally safe strategy to diminish populations of agricultural and horticultural insect pests. Functional characterization of genes related to male fertility can enhance the genetic-based SIT. Tssk1 has been involved to control male fertility in both mammals and insects. Moreover, Tektin1 has also been revealed to influence male fertility in both human and mammals. These findings suggested that Tssk1 and Tektin1 identified from Bactrocera dorsalis could be required for male fertility in B. dorsalis. In this study, expression profiles of these two genes were studied at different developmental stages and in various tissues of adult males. Remarkably, it was found that Tssk1 and Tektin1 were highly expressed in the testis of mature adult males of B. dorsalis. Furthermore, Tssk1 and Tektin1 genes were downregulated by using the RNA interference (RNAi) method. Fertility assays including egg laying, hatching, and spermatozoa count were also performed to investigate male fertility of B. dorsalis. Results showed that knockdown of Tssk1 and Tektin1 caused male sterility up to 58.99% and 64.49%, respectively. As expected, the total numbers of spermatozoa were also significantly reduced by 65.83% and 73.9%, respectively. These results suggested that male sterility was happened wing to the low number of spermatozoa. In conclusion, we demonstrate that Tssk1 and Tektin1 are the novel agents that could be used to enhance the genetic-based SIT, or their double-stranded RNA (dsRNA) can be used as biopesticides to control the population of B. dorsalis.

Gene expression profiling of ovary identified eggshell proteins regulated by 20-hydroxyecdysone in Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive pests worldwide. The frequent use of chemical insecticides has led B. dorsalis to develop resistance to many insecticides in recent decades. New high-throughput-sequenced transcriptomes, as well as genomes, have revealed a large number of reference genes for functional target identification. Here, we performed digital gene expression profiling of ovary and testis of B. dorsalis adults. Various genes were identified to be highly expressed in B. dorsalis ovary. The genes encoding components of eggshell, vitelline membrane proteins (Vmps) and chorion-related proteins, were identified to be tissue-specifically expressed in ovary. Five cytochrome P450 genes were also identified to be highly expressed in ovary. Three of them were ecdysone synthesis pathway genes indicating the ovary as a potential synthesis site of female. The up-regulated expression of Vmps by exogenous 20-hydroxyecdysone implied the hormonal regulation of eggshell formation during ovarian development. Many other genes with potential functions in ovarian development were also identified, including vitellogenin receptor, insulin receptor, NASP protein, and odorant binding protein. These findings should promote our understanding of the regulation of vitellogenesis and eggshell formation and enable exploration of potentially novel pest control targets.

Genome-Wide Analysis of MicroRNAs in Relation to Pupariation in Oriental Fruit Fly

Insect metamorphosis is a complex process involving drastic morphological and physiological changes. microRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play key roles in regulating various biological processes, including metamorphosis, by post-transcriptional repression of mRNAs. The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive insect pests in many Asian countries and the Pacific Islands. The regulatory role of miRNAs in B. dorsalis metamorphosis is unclear. To better understand the molecular regulatory mechanisms of miRNAs in pupariation, Illumina sequencing of the wandering stage (WS), the late WS and the white puparium stage of B. dorsalis were performed. Two hundred forty-nine miRNAs, including 184 known miRNAs and 65 novel miRNAs, were obtained. Among these miRNAs, 19 miRNAs were differentially expressed in pupariation, and eight miRNAs showed relative high expression levels (>50 TPM), of which five differentially expressed miRNAs (DEMs) had target differentially expressed genes (DEGs) predicted by the expected miRNA-mRNA negative regulation pattern using the Illumina HiSeq data. Four sets of DEMs and their predicted target DEGs were confirmed by qPCR. Of the four miRNAs, two miRNAs were down-regulated: miR-981, which may target pdpc, and Bdo-novel-mir-55, which potentially regulates spsX1, psB/C, and chit3. The other two miRNAs were up-regulated: let-7a-3p, which possibly controls lap, and Bdo-novel-mir-24, which may regulate ipc and sp1/2. This study provides a useful resource to elucidate the regulatory role of miRNAs and understand the molecular mechanisms of metamorphosis.

RNA Interference-Mediated Knockdown of Male Fertility Genes in the Queensland Fruit Fly Bactrocera tryoni (Diptera: Tephritidae)

The Queensland fruit fly, Bactrocera tryoni, is Australia’s most important horticultural pest. The Sterile Insect Technique (SIT) has been used to control this species for decades, using radiation to sterilize males before field-release. This method of sterilization can potentially reduce the insects’ abilities to compete for mates. In this study, RNA interference (RNAi) techniques were examined for their potential to sterilize male B. tryoni without adversely affecting mating competitiveness. B. tryoni adults were injected or fed double-stranded RNAs (dsRNAs) targeting spermatogenesis genes (tssk1, topi and trxt); quantitative reverse-transcriptase PCR analyses confirmed that transcript levels were reduced 60–80% for all three genes following injections. Feeding produced a significant gene knockdown for tssk1 and trxt after three days, but interestingly, two genes (trxt and topi) produced an excess of transcripts after 10 days of feeding. Despite these fluctuations in transcript levels, all three dsRNAs impacted the fecundity of treated males, with tssk1- and topi-dsRNA-treated males producing 75% fewer viable offspring than the negative controls. Mating competition assays demonstrated that dsRNA-treated males can actively compete with untreated males. These findings suggest that RNAi technology could serve as an alternative to radiation as a means of sterilizing these insects in an SIT program.

Genome-wide annotation of cuticular proteins in the oriental fruit fly (Bactrocera dorsalis), changes during pupariation and expression analysis of CPAP3 protein genes in response to environmental stresses

Cuticular proteins (CPs) are essential components of the insect cuticle as they create a structural and protective shield and may have a role in insect development. In this paper, we studied the CPs in the oriental fruit fly (Bactrocera dorsalis), one of the most economically important pests in the Tephritidae family around the world. The availability of a complete genome sequence (NCBI Assembly: ASM78921v2) allowed the identification of 164 CP genes in B. dorsalis. Comparative analysis of the CPs in B. dorsalis with those in the model insect Drosophila melanogaster and the closely related Ceratitis capitata, and CPs from mosquitoes, Lepidoptera, Hymenoptera and Coleoptera identified Diptera-specific genes and cuticle development patterns. Analysis of their evolutionary relationship revealed that some CP families had evolved according to the phylogeny of the different insect species, while others shared a closer relationship based on domain architecture. Subsequently, transcriptome analysis showed that while most of the CPs (60-100% of the family members) are expressed in the epidermis, some were also present in internal organs such as the fat body and the reproductive organs. Furthermore, the study of the expression profiles throughout development revealed a profound change in the expression of CPs during the formation of the puparium (pupariation). Further analysis of the expression profiles of the CPAP3 genes under various environmental stresses revealed them to be involved in the response to pesticides and arid and extreme temperatures conditions. In conclusion, the data provide a particular overview of CPs and their evolutionary and transcriptional dynamics, and in turn they lay a molecular foundation to explore their roles in the unique developmental process of insect metamorphosis and stress responses.

Proteomic analysis of watery saliva secreted by white-backed planthopper, Sogatella furcifera

The white-backed planthopper, Sogatella furcifera, is a phloem sap feeder that secretes watery and gelling saliva during feeding. In this study, we identified the major proteins in watery saliva of S. furcifera by shotgun LC-MS/MS analysis combined with transcriptomic analysis. A total of 161 proteins were identified, which were divided into 8 function categories, including enzymes, transporter, calcium ion binding protein, salivary sheath protein, cytoskeleton protein, DNA-, RNA-, and protein-binding or regulating proteins, other non-enzyme proteins and unknown proteins. Gene expression pattern of 11 secretory proteins were analyzed by real time quantitative-PCR. We detected the mucin-like protein, which had a unique expression level in salivary gland, most likely to be a candidate effector involved in regulation of plant defense. This study identified the watery saliva component of S. furcifera and it provided a list of proteins which may play a role in interaction between S. furcifera and rice.

Involvement of Met and Kr-h1 in JH-Mediated Reproduction of Female Bactrocera dorsalis (Hendel)

Juvenile hormone (JH) prevents metamorphosis during insect larval stages and promotes adult reproductive processes. Krüppel-homolog 1 (Kr-h1), a zinc finger transcription factor assumed to be induced by JH via the JH receptor methoprene-tolerant (Met), mediates the antimetamorphic effect of JH in insects, but its function in JH-mediated reproductive processes has not been fully explored. In this study, Met and Kr-h1 involved in the JH signaling pathway were first cloned and identified from the oriental fruit fly, Bactrocera dorsalis, an important pest infesting fruit and vegetables worldwide. Subsequent spatiotemporal expression analysis revealed that Met and Kr-h1 were both highly expressed in 7-day-old adults and fat body of female adults, respectively. Treatment with a JH analog (methoprene) significantly induced the expression of JH signaling and vitellogenin (Vg) genes and accelerated ovary development. RNA interference (RNAi) further revealed that either Met or Kr-h1 depletion at the adult stage of B. dorsalis impeded ovary development, with significantly lower egg production noted as well. In addition, rescue through methoprene application after RNAi stimulated the expression of JH signaling and Vg genes. Although there were still differences in ovary phenotype between rescued insects and the pre-RNAi control, ovary redevelopment with a larger surface area was observed, consistent with the spatiotemporal expression and phenotypes recorded in the original methoprene experiment. Our data reveal the involvement of Met and Kr-h1 in insect vitellogenesis and egg production, thus indicating the crucial role of the JH signaling pathway in insect reproduction.

The neuropeptides and protein hormones of the agricultural pest fruit fly Bactrocera dorsalis: What do we learn from the genome sequencing and tissue-specific transcriptomes?

Neuropeptides and protein hormones are very important signaling molecules, and are involved in the regulation and coordination of various physiological processes in invertebrates and vertebrates. Using a bioinformatics approach, we screened the recently sequenced genome and six tissue-specific transcriptome databases (central nervous system, fat body, ovary, testes, male accessory glands, antennae) of the oriental fruit fly (Bactrocera dorsalis) that is economically one of the most important pest insects of tropical and subtropical fruit. Thirty-nine candidate genes were found to encode neuropeptides or protein hormones. These include most of the known insect neuropeptides and protein hormones, with the exception of adipokinetic hormone-corazonin-related peptide, allatropin, diuretic hormone 34, diuretic hormone 45, IMFamide, inotocin, and sex peptide. Our results showed the neuropeptides and protein hormones of Diptera insects appear to have a reduced repertoire compared to some other insects. Moreover, there are also differences between B. dorsalis and the super-model of Drosophila melanogaster. Interesting features of the oriental fruit fly are the absence of genes coding for sex peptide and the presence of neuroparsin and two genes coding neuropeptide F. The majority of the identified neuropeptides and protein hormones is present in the central nervous system, with only a limited number of these in the other tissues. Moreover, we predicted their physiological functions via comparing with data of FlyBase and FlyAtlas. Taken together, owing to the large number of identified peptides, this study can be used as a reference about structure, tissue distribution and physiological functions for comparative studies in other model and important pest insects.

Transcriptome profiling of the Plutella xylostella (Lepidoptera: Plutellidae) ovary reveals genes involved in oogenesis

BACKGROUND

As a specialized organ, the insect ovary performs valuable functions by ensuring fecundity and population survival. Oogenesis is the complex physiological process resulting in the production of mature eggs, which are involved in epigenetic programming, germ cell behavior, cell cycle regulation, etc. Identification of the genes involved in ovary development and oogenesis is critical to better understand the reproductive biology and screening for the potential molecular targets in Plutella xylostella, a worldwide destructive pest of economically major crops.

RESULTS

Based on transcriptome sequencing, a total of 7.88Gb clean nucleotides was obtained, with 19,934 genes and 1861 new transcripts being identified. Expression profiling indicated that 61.7% of the genes were expressed (FPKM≥1) in the P. xylostella ovary. GO annotation showed that the pathways of multicellular organism reproduction and multicellular organism reproduction process, as well as gamete generation and chorion were significantly enriched. Processes that were most likely relevant to reproduction included the spliceosome, ubiquitin mediated proteolysis, endocytosis, PI3K-Akt signaling pathway, insulin signaling pathway, cAMP signaling pathway, and focal adhesion were identified in the top 20 'highly represented' KEGG pathways. Functional genes involved in oogenesis were further analyzed and validated by qRT-PCR to show their potential predominant roles in P. xylostella reproduction.

CONCLUSIONS

Our newly developed P. xylostella ovary transcriptome provides an overview of the gene expression profiling in this specialized tissue and the functional gene network closely related to the ovary development and oogenesis. This is the first genome-wide transcriptome dataset of P. xylostella ovary that includes a subset of functionally activated genes. This global approach will be the basis for further studies on molecular mechanisms of P. xylostella reproduction aimed at screening potential molecular targets for integrated pest management.

Comparative Proteomic Profiling Reveals Molecular Characteristics Associated with Oogenesis and Oocyte Maturation during Ovarian Development of Bactrocera dorsalis (Hendel)

Time-dependent expression of proteins in ovary is important to understand oogenesis in insects. Here, we profiled the proteomes of developing ovaries from Bactrocera dorsalis (Hendel) to obtain information about ovarian development with particular emphasis on differentially expressed proteins (DEPs) involved in oogenesis. A total of 4838 proteins were identified with an average peptide number of 8.15 and sequence coverage of 20.79%. Quantitative proteomic analysis showed that a total of 612 and 196 proteins were differentially expressed in developing and mature ovaries, respectively. Furthermore, 153, 196 and 59 potential target proteins were highly expressed in early, vitellogenic and mature ovaries and most tested DEPs had the similar trends consistent with the respective transcriptional profiles. These proteins were abundantly expressed in pre-vitellogenic and vitellogenic stages, including tropomyosin, vitellogenin, eukaryotic translation initiation factor, heat shock protein, importin protein, vitelline membrane protein, and chorion protein. Several hormone and signal pathway related proteins were also identified during ovarian development including piRNA, notch, insulin, juvenile, and ecdysone hormone signal pathways. This is the first report of a global ovary proteome of a tephritid fruit fly, and may contribute to understanding the complicate processes of ovarian development and exploring the potentially novel pest control targets.

Discovery and functional identification of fecundity-related genes in the brown planthopper by large-scale RNA interference

Recently, transcriptome and proteome data have increasingly been used to identify potential novel genes related to insect phenotypes. However, there are few studies reporting the large-scale functional identification of such genes in insects. To identify novel genes related to fecundity in the brown planthopper (BPH), Nilaparvata lugens, 115 genes were selected from the transcriptomic and proteomic data previously obtained from high- and low-fecundity populations in our laboratory. The results of RNA interference (RNAi) feeding experiments showed that 91.21% of the genes were involved in the regulation of vitellogenin (Vg) expression and may influence BPH fecundity. After RNAi injection experiments, 12 annotated genes were confirmed as fecundity-related genes and three novel genes were identified in the BPH. Finally, C-terminal binding protein (CtBP) was shown to play an important role in BPH fecundity. Knockdown of CtBP not only led to lower survival, underdeveloped ovaries and fewer eggs laid but also resulted in a reduction in Vg protein expression. The novel gene resources gained from this study will be useful for constructing a Vg regulation network and may provide potential target genes for RNAi-based pest control.

Characterization of a β-Adrenergic-Like Octopamine Receptor in the Oriental Fruit Fly, Bactrocera dorsalis (Hendel)

The biogenic amine octopamine plays a critical role in the regulation of many physiological processes in insects. Octopamine transmits its action through a set of specific G-protein coupled receptors (GPCRs), namely octopamine receptors. Here, we report on a β-adrenergic-like octopamine receptor gene (BdOctβR1) from the oriental fruit fly, Bactrocera dorsalis (Hendel), a destructive agricultural pest that occurs in North America and the Asia-Pacific region. As indicated by RT-qPCR, BdOctβR1 was highly expressed in the central nervous system (CNS) and Malpighian tubules (MT) in the adult flies, suggesting it may undertake important roles in neural signaling in the CNS as well as physiological functions in the MT of this fly. Furthermore, its ligand specificities were tested in a heterologous expression system where BdOctβR1 was expressed in HEK-293 cells. Based on cyclic AMP response assays, we found that BdOctβR1 could be activated by octopamine in a concentration-dependent manner, confirming that this receptor was functional, while tyramine and dopamine had much less potency than octopamine. Naphazoline possessed the highest agonistic activity among the tested agonists. In antagonistic assays, mianserin had the strongest activity and was followed by phentolamine and chlorpromazine. Furthermore, when the flies were kept under starvation, there was a corresponding increase in the transcript level of BdOctβR1, while high or low temperature stress could not induce significant expression changes. The above results suggest that BdOctβR1 may be involved in the regulation of feeding processes in Bactrocera dorsalis and may provide new potential insecticide leads targeting octopamine receptors.

Testis-specific RNA-Seq of Calligrapha (Chrysomelidae) as a transcriptomic resource for male-biased gene inquiry in Coleoptera

We report the architecture of testis transcriptomes of four closely related species of Calligrapha (Chrysomelidae) beetles, which diverged during the last 3 million years. Five cDNA libraries were sequenced using Illumina HiSeq technology, retrieving 102 884-176 514 assembled contigs, of which ~33-45% of these longer than 499 nt were functionally annotated. Annotation and sequence similarity comparisons of these libraries revealed high homogeneity in gene composition and the presence of several functional candidates related to reproduction or reproductive processes (0.72-1.08% of annotated sequences). Stringent sequence similarity analyses of these transcriptomes against empirically demonstrated male-biased genes in Drosophila melanogaster and Tribolium castaneum allowed the identification of 77 homologues in Calligrapha, possible candidates of male-biased expression. Some of these genes - including CG9313, Tektin-A or tomboy40 - were confirmed as orthologs of these male-biased genes using phylogenetic inference and available model insect data, increasing our confidence that they represent functional homologues too. Our transcriptomes are a valuable transcriptomic resource for the analysis of male-biased genes in Calligrapha, which has the added interest of including several female-only species. But it simultaneously represents a landmark for similar studies in Coleoptera, broadening the taxonomic diversity currently represented by the model species T. castaneum, and incipient genomic data in other herbivorous lineages, including weevils, longhorn beetles and leaf beetles.

Transcriptome analysis to identify genes for peptides and proteins involved in immunity and reproduction from male accessory glands and ejaculatory duct of Bactrocera dorsalis

In the male reproductive system of insects, the male accessory glands and ejaculatory duct (MAG/ED) are important organs and their primary function is to enhance the fertility of spermatozoa. Proteins secreted by the MAG/ED are also known to induce post-mating changes and immunity responses in the female insect. To understand the gene expression profile in the MAG/ED of the oriental fruit fly Bactrocera dorsalis (Hendel), that is an important pest in fruits, we performed an Illumina-based deep sequencing of mRNA. This yielded 54,577,630 clean reads corresponding to 4.91Gb total nucleotides that were assembled and clustered to 30,669 unigenes (average 645bp). Among them, 20,419 unigenes were functionally annotated to known proteins/peptides in Gene Orthology, Clusters of Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes pathway databases. Typically, many genes were involved in immunity and these included microbial recognition proteins and antimicrobial peptides. Subsequently, the inducible expression of these immunity-related genes was confirmed by qRT-PCR analysis when insects were challenged with immunity-inducible factors, suggesting their function in guaranteeing fertilization success. Besides, we identified some important reproductive genes such as juvenile hormone- and ecdysteroid-related genes in this de novo assembly. In conclusion, this transcriptomic sequencing of B. dorsalis MAG/ED provides insights to facilitate further functional research of reproduction, immunity and molecular evolution of reproductive proteins in this important agricultural pest.

RNA sequencing to characterize transcriptional changes of sexual maturation and mating in the female oriental fruit fly Bactrocera dorsalis

BACKGROUND

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.

RESULTS

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.

CONCLUSION

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.

Deciphering the main venom components of the ectoparasitic ant-like bethylid wasp, Scleroderma guani

Similar to venom found in most venomous animals, parasitoid venoms contain a complex cocktail of proteins with potential agrichemical and pharmaceutical use. Even though parasitoids are one of the largest group of venomous animals, little is known about their venom composition. Recent few studies revealed high variated venom composition existing not only in different species but also between closely related strains, impling that increasing information on the venom proteins from more greater diversity of species of different taxa is key to comprehensively uncover the complete picture of parasitoid venom. Here, we explored the major protein components of the venom of ectoparasitic ant-like bethylid wasp, Scleroderma guani by an integrative transcriptomic-proteomic approach. Illumina deep sequencing of venom apparatus cDNA produced 49,873 transcripts. By mapping the peptide spectral data derived from venom reservoir against these transcripts, mass spectrometry analysis revealed ten main venom proteins, including serine proteinase, metalloprotease, dipeptidyl peptidase IV, esterase, antithrombin-III, acid phosphatase, neural/ectodermal development factor IMP-L2 like protein, venom allergen 3, and unknown protein. Interestingly, one serine proteinase was firstly identified with rarely high molecular weight about 200 kDa in parasitoid venom. The occurrence of abundant acid phosphatase, antithrombin-III and venom allergen 3 demonstrated that S. guani venom composition is similar to that of social wasp venoms. All identified venom genes showed abundantly biased expression in venom apparatus, indicating their virulent functions involved in parasitization. This study shed light on the more better understanding of parasitoid venom evolution across species and will facilitate the further elucidation of function and toxicity of these venom proteins.

Identification, characterization and target gene analysis of testicular microRNAs in the oriental fruit fly Bactrocera dorsalis

MicroRNAs (miRNAs) are small noncoding RNAs that regulate various diverse biological processes including insect spermatogenesis. The oriental fruit fly, Bactrocera dorsalis, is one of the most destructive horticultural pests in East Asia and the Pacific region. Although developmental miRNA profiles of B. dorsalis have enriched our knowledge, specific testicular miRNAs in this dipteran species are unexplored. In this study, we identified miRNAs from B. dorsalis testes by deep sequencing, which provided an overview of miRNA expression during spermatogenesis. Small RNA libraries were constructed from the testes of fully mature (FM), immature (IM) and middle-aged (MA) adult flies of B. dorsalis. Small RNA sequencing and data analysis revealed 172 known and 78 novel miRNAs amongst these libraries. Pairwise comparisons of libraries led to the identification of 24, 15 and 14 differentially expressed miRNAs in FM vs. IM, FM vs. MA and IM vs. MA insects, respectively. Using a bioinformatic approach, we predicted 124 target genes against the 13 most differentially expressed miRNAs. Furthermore, the expression patterns of six randomly selected miRNAs (from the 13 most differentially expressed miRNAs) and their putative target genes (from the 124 predicted target genes) were analysed in the testis of B. dorsalis by quantitative real-time PCR, which showed that out of six, four tested miRNAs-mRNAs had an inverse expression pattern and are probably co-regulated. This study is the first comparative profile of the miRNA transcriptome in three developmental stages of the testis, and provides a useful resource for further studies on the role of miRNAs in spermatogenesis in B. dorsalis.