Next Article in Journal
Mayflies (Ephemeroptera) and Their Contributions to Ecosystem Services
Previous Article in Journal
Identification and Virulence Characterization of Two Akanthomyces attenuatus Isolates Against Megalurothrips usitatus (Thysanoptera: Thripidae)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 10
Issue 6
10.3390/insects10060169
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae)

by
Jia Wang
*,
Huan Fan
,
Pan Wang
and
Ying-Hong Liu
College of Plant Protection, Southwest University, Chongqing 400716, China
*
Author to whom correspondence should be addressed.
Insects 2019, 10(6), 169; https://doi.org/10.3390/insects10060169
Submission received: 29 April 2019 / Revised: 3 June 2019 / Accepted: 11 June 2019 / Published: 13 June 2019
Browse Figures
Versions Notes

Abstract

:
The Chinese citrus fly, Bactrocera minax, is a devastating pest of citrus, which enters the obligatory diapause in overwintering pupae to resist harsh environmental conditions. However, little is known about the molecular mechanisms underlying pupal diapause. The previous transcriptomic analysis revealed that a large number of genes were regulated throughout the pupal stage. Of these genes, 12 and six ones that are remarkably up- and downregulated, respectively, specifically in intense diapause were manually screened out in present study. To validate the expression of these genes throughout the pupal stage, the quantitative real-time PCR (qRT-PCR) was conducted, and the genes displaying different expression patterns with those of previous study were excluded. Then, the expressions of remaining genes were compared between diapause-destined and non-diapause-destined pupae to reveal their association with diapause using qRT-PCR and semiquantitative PCR. Finally, five genes, TTLL3B, Cyp6a9, MSTA, Fru, and UC2, and two genes, KSPI and LYZ1, were demonstrated to be positively and negatively associated with diapause, respectively. These findings provide a solid foundation for the further investigation of molecular mechanisms underlying B. minax pupal diapause.

1. Introduction

The Chinese citrus fly, Bactrocera minax (Enderlein) (Diptera: Tephritidae), an oligophagous pest of citrus fruits, has been considered as a devastating pest of citrus plants in the temperate areas of Asia [1,2]. Usually, the female adults lay eggs into the citrus fruit where the larvae hatch and feed, causing the fruits to ripen prematurely and drop down. Subsequently, the mature larvae drill out of the fallen fruits and burrow into soil for pupation [2]. Recently, the outbreaks of B. minax in several provinces of China have caused serious economic losses to the citrus industry [2,3]. Therefore, this pest has aroused great concerns in citrus-growing regions in China and the research on this pest has been widely conducted [4,5,6,7,8,9,10,11].
Diapause is a genetically programmed developmental arrest, accompanied by suppressed metabolism and enhanced stress tolerance [12,13], in response to environmental stimuli or internal cues to resist adverse environmental conditions in unfavorable season [14]. Many univoltine insects enter obligatory diapause at a specific developmental stage in each generation but require no token stimuli for diapause induction and preparation [12]. Likewise, the univoltine B. minax also enters the pupal diapause to overwinter. The overwintering pupal stage lasts for 160–170 d, synchronizing the adult emergence with citrus fruit bearing [15]. Therefore, understanding the molecular mechanisms of pupal diapause of B. minax will contribute to elucidating the inherent mechanisms underlying pupal development and adaptation to harsh environment. In the past decade, several relevant studies on pupal diapause of B. minax had been carried out in China. For example, the previous studies have demonstrated that the chilling in diapause is critical for pupal survival and adult emergence [10], and the 20E application was able to terminate pupal diapause and significantly accelerate adult emergence [7,16]. Moreover, the respiratory rate throughout the pupal stages was measured, and the fitted respiratory rate trajectory indicated that B. minax gradually entered intense diapause approximately 40 d after pupation, and the intense diapause lasted for two months (40–100 d after pupation) [17]. Although substantial progress has been made, the understanding of mechanisms underlying pupal diapause of B. minax is still limited.
Mining the diapause-associated genes could be an effective way in an effort to reveal the mechanisms underlying diapause. A number of studies have been conducted in this endeavor and many diapause-associated genes have been discovered in insects [18,19,20,21]. Recently, the next-generation sequencing technology has widely been used to study the mechanisms underlying a lot of biological processes in organisms by revealing the gene expression profiles in high-throughput and large-scale manners [22,23,24]. Given the importance of insect diapause, the underlying mechanisms have been studied using next-generation sequencing technology in many insects [25,26,27]. In our previous study, the transcriptome of B. minax was obtained using next-generation sequencing technology, based on which the gene expression profiles throughout pupal stage were compared, and a large number of genes were found to be regulated specifically in diapause [17]. In order to reveal the physiological variation in diapause, genes regulated most remarkably during diapause were manually screened out in this study, according to the reads per kb of exon model per million mapped reads (RPKM) values obtained previously [17]. Then, the quantitative real-time PCR (qRT-PCR) was conducted to verify the expression patterns of these genes throughout the pupal stage. Those genes which displayed inconsistent expression patterns between qRT-PCR and transcriptomic analysis were excluded. To further verify the association of remaining genes with diapause, the qRT-PCR and semiquantitative PCR were performed to compare the gene expression levels between diapause-destined and non-diapause-destined pupae. The findings of this study laid basis for further investigations on the mechanisms underlying B. minax diapause.

2. Materials and Methods

2.1. Insect Rearing and Sample Collection

Fallen oranges infested with maggots were taken back to the lab from an orchard in Wulong County, Chongqing Municipality, China. Third-instar larvae were collected from oranges and placed over sand in plastic dishes to pupate. All pupae were reared at 18 ± 2 °C, 70 ± 10% relative humidity, and photoperiod of 14 L:10 D. Twenty pupae were collected and stored in liquid nitrogen every 10 days until adult emergence for subsequent gene expression analysis.

2.2. Acquisition of Non-Diapause-Destined Pupae and Sample Collection

Non-diapause-destined pupae were acquired by injecting 20E into newly-formed pupae. The 20E (Sigma, St. Louis, MO, USA) was dissolved in 10% ethanol to concentration of 1 μg/μL. Each newly-formed pupa with the same size (≈9 mm) was injected with 1 μL 20E solution. Then all pupae were reared in plastic dishes under the same condition described above. Twenty pupae were collected and stored in liquid nitrogen every 10 days for subsequent gene expression analysis.

2.3. RNA Isolation and First-Strand cDNA Synthesis

Total RNA was isolated from whole body of pupa using TRIzol Reagent (Life Technologies, Carlsbad, CA, USA) according to manufacturer’s protocols. The integrity of RNA was detected by 1% agarose gel electrophoresis. Reverse transcription of 500 ng total RNA into the first-strand cDNA was performed using a PrimeScriptTM RT Master Mix (Perfect Real Time) Kit (Takara, Shiga, Japan). Three biological replicates were set for each time point.

2.4. Selection of Putative Diapause-Associated Genes

The previous transcriptomic study discovered significantly differentially expressed genes (DEGs) throughout pupal stage according to the RPKM values. The cluster analysis indicated that 1434 and 848 DEGs were up- and downregulated, respectively, specifically in diapause (40–100 d after pupation) [17]. In this study, the top hit DEGs with stable expression among biological replicates were manually selected as putative highly diapause-associated genes. These genes were annotated by alignment to NCBI Nr database using BLASTx (http://www.ncbi.nlm.nih.gov/) (Table 1).

2.5. Identification of Highly Diapause-Associated Genes

The expression profiles of 18 selected genes throughout pupal stage were analyzed by qRT-PCR. The genes displaying different expression patterns between qRT-PCR analysis and previous transcriptomic study were excluded. Then, the expression levels of remaining genes were compared between diapause-destined and non-diapause-destined pupae, using qRT-PCR and semiquantitative PCR, to determine the highly diapause-associated genes. The specific primers for 18 putative genes were designed (Table S1), and UBQ was used as the reference gene [28]. The qRT-PCR was performed with a CFX96TM Real-Time PCR Detection System thermal cycler (BIO-RAD, Hercules, CA, USA) in a reaction volume of 10 μL, including 4 μL SYBR® Premix Ex Taq II (Takara), 1 μL cDNA (50 ng/μL), and 0.5 μL forward and reverse primers (10 μM). Relative expression level was calculated according to the 2−ΔΔCt method. Three technique replicates were set for each sample. For semiquantitative PCR, the reaction volume was 25 µL, including 2 µL 10× PCR reaction buffer, 2.5 μL dNTP Mix (10 mM), 2.5 μL MgCl2 (25 mM), 0.25 μL Taq polymerase (2.5 U), 1 μL forward primer (10 μM) and reverse primer (10 μM), 1 μL cDNA template, and 14.8 μL ddH2O. The PCR reaction was carried out under conditions: 94 °C for 5 min, 28–32 cycles of 94 °C for 30 s followed by 55–60 °C 30 s, 72 °C 30 s, 72 °C 5 min (Table S1). To confirm the identity of highly diapause-associated genes, phylogenetic analysis was performed by Molecular Evolutionary Genetics Analysis version 5 (MEGA 5, https://www.megasoftware.net/) with the neighbor-joining method on the basis of the amino acid sequences available in GenBank (http://www.ncbi.nlm.nih.gov/). The reliability of the branching was tested using bootstrap method (1000 replications).

2.6. Statistical Analysis

The relative expression level of each gene throughout pupal stages was compared using one-way analysis of variance (ANOVA) with Tukey’s test after verification of normality and homogeneity of variances. The relative expression level of each gene was compared between diapause-destined and non-diapause-destined pupae using independent t-test. All data analyses were performed with SPSS 22.0 software (SPSS Inc., Chicago, IL, USA).

3. Results

3.1. Expression Profiles of Selected Genes throughout Pupal Stage

A total of 12 upregulated and six downregulated genes during diapause were manually selected as putative diapause-associated genes according to the RPKM values (Figure 1). All eighteen genes were significantly regulated throughout the pupal stage. However, only 13 genes, Arrdc3 (F13,41 = 12.931, p < 0.05), Cyp6a9 (F13,41 = 37.071, p < 0.05), Fru (F13,41 = 19.433, p < 0.05), GEF (F13,41 = 14.4, p < 0.05), KSPI (F13,41 = 28.247, p < 0.05), LYZ1 (F13,41 = 524.051, p < 0.05), MSTA (F13,41 = 205.801, p < 0.05), MYBLI (F13,41 = 18.417, p < 0.05), setmar (F13,41 = 11.244, p < 0.05), TTLL3B (F13,41 = 25.601, p < 0.05), twk-7 (F13,41 = 66.57, p < 0.05), UC1 (F13,41 = 19.150, p < 0.05), and UC2 (F13,41 = 23.198, p < 0.05), showed similar expression patterns with those of previous transcriptomic study. The remaining five genes, ECE-2L (F13,41 = 251.969, p < 0.05), grk (F13,41 = 24.785, p < 0.05), MEGF10 (F13,41 = 7.35, p < 0.05), Pcd (F13,41 = 2.672, p < 0.05), and wbl (F13,41 = 17.428, p < 0.05), were excluded. The expression levels of Arrdc3, Cyp6a9, Fru, GEF, MSTA, MYBLI, setmar, TTLL3B, UC1, and UC2 were upregulated along with intensification of diapause and downregulated after diapause termination, although the differences were not significant at certain time points. The expression levels of KSPI, LYZ1, and twk-7 were remarkably suppressed 10 days after pupation and upregulated after diapause termination (Figure 2).

3.2. Comparison of Gene Expression Levels between Diapause-Destined and Non-Diapause-Destined Pupae

The newly-formed pupae were injected with 20E to acquire non-diapause-destined pupae. To validate the association of 13 genes mentioned above with diapause, their expression levels were compared between diapause-destined and non-diapause-destined pupae using qRT-PCR. Nine genes, Arrdc3, Cyp6a9, Fru, KSPI, LYZ1, MSTA, setmar, TTLL3B, and UC2, were found to be putatively diapause-associated. The expression of Arrdc3, Cyp6a9, Fru, MSTA, setmar, TTLL3B, and UC2 were suppressed, while that of KSPI and LYZ1 were activated in non-diapause-destined pupae (Figure 3). In addition, the MYBLI, twk-7, UC1, and GEF were not regulated as expected (Figure 3).
Subsequently, the semiquantitative PCR was conducted to further validate the association of these nine genes with diapause. The expression levels of 7 genes, Cyp6a9, Fru, KSPI, LYZ1, MSTA, TTLL3B, and UC2 were consistent with the qRT-PCR data, displaying significant differences between diapause-destined and non-diapause-destined pupae. Nevertheless, Arrdc3 and setmar only exhibited slight differences between diapause-destined and non-diapause-destined pupae (Figure 4). Taken together, Cyp6a9, Fru, MSTA, TTLL3B, and UC2 are highly positively associated with diapause, while KSPI and LYZ1 are highly negatively associated with diapause. The identities of these diapause-associated genes, except UC2, were confirmed by phylogenetic analyses (Figures S1–S6).

4. Discussion

According to our previous transcriptomic study, over 4000 genes were regulated throughout the pupal stage of B. minax, and several physiological pathways were deployed in diapause programming [17]. Despite the high sensitivity and accuracy in revealing the gene expression profile, the transcriptomic analysis still has errors. Therefore, the subsequent qRT-PCR validation is usually adopted [17,25,26]. In this study, 18 DEGs of top hits and displaying stable expression among biological replicates were manually selected as putative diapause-associated genes. The expression patterns of these 18 genes throughout the pupal stage were evaluated by qRT-PCR analysis, and 13 of them were consistent with that of previous transcriptomic study. Then, the newly-formed pupae were injected with 20E and solvent to acquire non-diapause-destined and diapause-destined pupae, respectively. Along with pupal development, the expression levels of positively diapause-associated genes ought to be increased with intensification of diapause in diapause-destined pupae, but not in non-diapause-destined pupae. In contrast, the negatively diapause-associated genes were expected to be activated in non-diapause-destined pupae but suppressed in diapause-destined pupae. Therefore, the qRT-PCR and semiquantitative PCR were performed to assess the differences in expression levels of remaining 13 genes between diapause-destined and non-diapause-destined pupae. The findings of these two methods were not exactly consistent with each other as Arrdc3 and setmar displayed a significant difference in qRT-PCR analysis, but not in semiquantitative PCR analysis, probably due to the relatively higher sensitivity of qRT-PCR. Only seven genes which showed a significant difference in both analyses were deemed to be highly diapause-associated genes.
TTLL3B, a member of tubulin-tyrosine ligase (TTLL) family, catalyzes the polyglycylation of proteins, especially tubulin [29]. Polyglycylation is the long-chain glycines modification at a single or multiple glutamate sites on the C-terminal tail of protein, and the long-chain glycine is linked to the γ-carboxyl group of the glutamic acid residue via isopeptide bond. This modification consists of two steps, initiation and extend [30,31,32,33]. To date, glycylation has been found solely in axonemes of cilia and flagella [30,31,32,33,34,35], playing roles in stability of cell cilia and male sterility. However, little is known about its function in insects. To our knowledge, this is the first time that glycylation has been found associated with diapause, while its functions necessitate further investigation.
P450s are involved in the metabolism of exogenous substances such as insecticides and plant secondary metabolites, and also endogenous compounds like juvenile hormones, 20-hydroxyecdysone, ecdysteroids, and pheromones [36]. The CYP6a9 identified in this study belongs to P450 CYP6 family which is specific to insect [37]. Likewise, some members in P450 CYP6 family were also found to be associated with diapause of other insects, such as Drosophila montana [26]. It has been established that P450 genes are related with insecticide resistance, stress resistance, and immunity in organisms [38,39,40]. However, the roles CYP6a9 plays in B. minax diapause remains unknown.
MSTA is a member of SET (suppressor of variegation, enhancer of zeste, trithorax) and MYND (myeloid-Nervy-DEAF1) domain-containing proteins family (SMYD) which are special protein lysine methyltransferases involved in methylation of histones and nonhistone targets [41,42,43]. The SET domain is an evolutionarily conserved catalytic unit for lysine methylation [44,45], and the MYND domain is a conserved zinc-binding domain primarily involved in protein–protein interaction [46]. SMYD family consists of six groups, SMYD1–5 in mammal [47], and arthropod-specific SMYDA [48]. SMYD proteins are related to epigenetic control, development, and cell proliferation through posttranslational modifications in histones and nonhistones. In insects, the studies on SMYD proteins are still limited and mainly focused on fruit flies. Knockdown of Drosophila Smyd4 leads to abnormal eclosion, suggesting the vital role for Smyd4 in development [49]. That may also be the case in diapausing B. minax.
Fru gene encodes a sex-specifically spliced transcription factor containing the zinc finger and BTB (Broad-complex, Tramtrack, and Bric-à-brac) domain [50,51,52]. It plays a crucial role in the courtship behavior and neuronal sex differentiation of male fruit flies [51,53]. For example, the female fruit flies expressing male-specific Fru in the nervous system showed the courtship behavior to other females, while the male flies expressing female-specific Fru developed into sterile males, lost the intrinsic properties of males, and exhibited female characteristics [54,55]. Knockdown of Fru gene in Schistocerca gregaria interfered mating and also affected fertility and fecundity of males [56]. Therefore, the Fru gene has been considered as a switch gene for courtship and thus contributing to the survival and reproduction of insects. However, the role Fru plays in insect diapause has not yet been studied. It is likely that certain unknown downstream genes other than courtship-related ones were activated by Fru in diapausing B. minax.
KSPI identified in this study is a serine protease inhibitor which modulates the bioactivity of target proteins through proteolytic cleavage. KSPI belongs to Kunitz-type protease inhibitor family, also known as bovine pancreatic trypsin inhibitor (BPTI)-like proteases [57]. As a classic representation of this family, BPTI has a broad inhibitory activity and mainly inhibits trypsin, chymotrypsin, and elastase-like serine protease [58]. The structures of all known Kunitz-type protease inhibitors of invertebrates are all similar to BPTI [59]. Serine protease inhibitors play important roles in the innate immunity of insects. It has been demonstrated that Kunitz-type protease inhibitors can be deployed to defense against microorganism invasion in insects [60,61]. Therefore, the downregulation of KSPI may suppress the immune response in diapausing B. minax, which is consistent with the findings that the immunity was attenuated in diapausing Samia cynthia pryeri pupae [62].
Lysozyme is a vital bacteriolytic enzyme in insect innate immune system, catalyzing the hydrolysis of 1,4-beta-linkages between N-acetylmuramic acid and N-acetylglucosamine (NAG) residues in the peptidoglycan of the bacterial cell wall [63]. The lysozymes fall into six different classes: c-type, g-type, i-type, bacterial, T4 phage-type, and plant lysozymes based on the source and mechanism [64]. The LYZ1 found in this study belongs to c-type. Insect lysozymes are produced in fat body and released into hemolymph. The activity of lysozyme in hemolymph increases dramatically when insects are infected with pathogens [65,66]. Like KSPI, the LYZ1 was downregulated during B. minax diapause as well, supporting the speculation that the pupal immunity was suppressed during diapause.

5. Conclusions

The univoltine B. minax enters obligatory diapause in pupal stage, during which a large number of genes are regulated. In this study, 5 genes, TTLL3B, Cyp6a9, MSTA, Fru, and UC2, and 2 genes, KSPI and LYZ1, were demonstrated to be positively and negatively associated with diapause, respectively. Nevertheless, there could be other genes that are involved in the B. minax diapause programming were not analyzed in the present study. Therefore, more effort can be exerted to explore the relevant genes. To elucidate the mechanisms underlying B. minax diapause, the functions of these diapause-associated genes ought to be further investigated in the future.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/10/6/169/s1, Table S1: Primers used for qRT-PCR and Semiquantitative PCR, Figure S1: Phylogenetic analysis of the Cyp6a9, Figure S2: Phylogenetic analysis of the Fru, Figure S3: Phylogenetic analysis of the KSPI, Figure S4: Phylogenetic analysis of the LYZ1, Figure S5: Phylogenetic analysis of the MSTA, Figure S6: Phylogenetic analysis of the TTLL3.

Author Contributions

Conceptualization, J.W. and Y.-H.L.; methodology, J.W.; formal analysis, J.W. and H.F.; writing—original draft preparation, J.W. and H.F.; writing—review and editing, J.W.; supervision, J.W.; project administration, J.W.; funding acquisition, J.W.

Funding

This study was financially supported by the National Natural Science Foundation of China (31401742), the Chongqing Research Program of Basic Research and Frontier Technology (cstc2016jcyjA0203 and cstc2018jcyjAX0516), and the Fundamental Research Funds for the Central Universities (XDJK2018C092).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Dorji, C.; Clarke, A.; Drew, R.; Fletcher, B.; Loday, P.; Mahat, K.; Raghu, S.; Romig, M. Seasonal phenology of Bactrocera minax (Diptera: Tephritidae) in western Bhutan. Bull. Entomol. Res. 2006, 96, 531–538. [Google Scholar] [PubMed]
  2. Xia, Y.L.; Ma, X.L.; Hou, B.H.; Ouyang, G.C. A review of Bactrocera minax (Diptera: Tephritidae) in China for the purpose of safeguarding. Adv. Entomol. 2018, 6, 35–61. [Google Scholar] [CrossRef]
  3. Zhang, Y.A. Citrus fruit flies of Sichuan Province (China). EPPO Bull. 2010, 19, 649–654. [Google Scholar]
  4. Chen, E.H.; Dou, W.; Hu, F.; Tang, S.; Zhao, Z.M.; Wang, J.J. Purification and biochemical characterization of glutathione s-transferases in Bactrocera minax (Diptera: Tephritidae). Fla. Entomol. 2012, 95, 593–601. [Google Scholar] [CrossRef]
  5. Dong, Y.; Wan, L.; Rui, P.; Desneux, N.; Niu, C. Feeding and mating behaviour of Chinese citrus fly Bactrocera minax (Diptera, Tephritidae) in the field. J. Pest Sci. 2014, 87, 647–657. [Google Scholar] [CrossRef]
  6. Liu, H.; Jiang, G.; Zhang, Y.; Chen, F.; Li, X.; Yue, J.; Ran, C.; Zhao, Z. Effect of six insecticides on three populations of Bactrocera (Tetradacus) minax (Diptera: Tephritidae). Curr. Pharm. Biotechnol. 2015, 16, 77–83. [Google Scholar] [CrossRef]
  7. Wang, J.; Zhou, H.Y.; Zhao, Z.M.; Liu, Y.H. Effects of juvenile hormone analogue and ecdysteroid on adult eclosion of the fruit fly Bactrocera minax (Diptera: Tephritidae). J. Econ. Entomol. 2014, 107, 1519–1525. [Google Scholar] [CrossRef]
  8. Wang, X.L.; Zhang, R.J. Review on biology, ecology and control of Bactrocera (Tetradacus) minax Enderlein. J. Environ. Entomol. 2009, 31, 73–79. [Google Scholar]
  9. Zhang, B.; Nardi, F.; Hull-Sanders, H.; Wan, X.; Liu, Y. The complete nucleotide sequence of the mitochondrial genome of Bactrocera minax (Diptera: Tephritidae). PLoS ONE 2014, 9, e100558. [Google Scholar] [CrossRef]
  10. Dong, Y.C.; Wang, Z.J.; Clarke, A.R.; Rui, P.; Desneux, N.; Niu, C.Y. Pupal diapause development and termination is driven by low temperature chilling in Bactrocera minax. J. Pest Sci. 2013, 86, 429–436. [Google Scholar] [CrossRef]
  11. Lü, Z.C.; Wang, L.H.; Zhang, G.F.; Wan, F.H.; Guo, J.Y.; Yu, H.; Wang, J.B. Three heat shock protein genes from Bactrocera (tetradacus) minax Enderlein: Gene cloning, characterization, and association with diapause. Neotrop. Entomol. 2014, 43, 362–372. [Google Scholar] [CrossRef] [PubMed]
  12. Denlinger, D.L. Regulation of diapause. Annu. Rev. Entomol. 2002, 47, 93–122. [Google Scholar] [CrossRef] [PubMed]
  13. Koštál, V. Eco-physiological phases of insect diapause. J. Insect Physiol. 2006, 52, 113–127. [Google Scholar] [CrossRef] [PubMed]
  14. Ragland, G.J.; Egan, S.P.; Feder, J.L.; Berlocher, S.H.; Hahn, D.A. Developmental trajectories of gene expression reveal candidates for diapause termination: A key life-history transition in the apple maggot fly Rhagoletis pomonella. J. Exp. Biol. 2011, 214, 3948–3960. [Google Scholar] [CrossRef]
  15. Wang, X.J.; Luo, L.Y. Research progress in the Chinese citrus fruit fly. Entomol. Knowl. 1995, 32, 310–315. [Google Scholar]
  16. Chen, Z.; Dong, Y.; Wang, Y.; Andongma, A.A.; Rashid, M.A.; Krutmuang, P.; Niu, C. Pupal diapause termination in Bactrocera minax: An insight on 20-hydroxyecdysone induced phenotypic and genotypic expressions. Sci. Rep. 2016, 6, 27440. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, J.; Fan, H.; Xiong, K.C.; Liu, Y.H. Transcriptomic and metabolomic profiles of Chinese citrus fly, Bactrocera minax (Diptera: Tephritidae), along with pupal development provide insight into diapause program. PLoS ONE 2017, 12, e0181033. [Google Scholar] [CrossRef]
  18. Fremdt, H.; Amendt, J.; Zehner, R. Diapause-specific gene expression in Calliphora vicina (Diptera: Calliphoridae)—A useful diagnostic tool for forensic entomology. Int. J. Leg. Med. 2014, 128, 1001–1011. [Google Scholar] [CrossRef]
  19. Sasibhushan, S.; Ponnuvel, K.M.; Vijayaprakash, N.B. Diapause specific gene expression in the eggs of multivoltine silkworm Bombyx mori, identified by suppressive subtractive hybridization. Comp. Biochem. Phys. B 2012, 161, 371–379. [Google Scholar] [CrossRef]
  20. Hayward, S.A.L.; Pavlides, S.C.; Tammariello, S.P.; Rinehart, J.P.; Denlinger, D.L. Temporal expression patterns of diapause-associated genes in flesh fly pupae from the onset of diapause through post-diapause quiescence. J. Insect Physiol. 2005, 51, 631–640. [Google Scholar] [CrossRef]
  21. Chen, B.; Kayukawa, T.; Monteiro, A.; Ishikawa, Y. The expression of the HSP90 gene in response to winter and summer diapauses and thermal-stress in the onion maggot, Delia antiqua. Insect Mol. Biol. 2005, 14, 697–702. [Google Scholar] [CrossRef] [PubMed]
  22. Mardis, E.R. Next-generation DNA sequencing methods. Annu. Rev. Genom. Hum. Genet. 2008, 9, 387–402. [Google Scholar] [CrossRef] [PubMed]
  23. Margulies, M.; Egholm, M.; Altman, W.E.; Attiya, S.; Bader, J.S.; Bemben, L.A.; Berka, J.; Braverman, M.S.; Chen, Y.J.; Chen, Z. Genome sequencing in microfabricated high-density picolitre reactors. Nature 2005, 437, 376–380. [Google Scholar] [CrossRef] [PubMed]
  24. Morozova, O.; Marra, M.A. Applications of next-generation sequencing technologies in functional genomics. Genomics 2008, 92, 255–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Deng, Y.; Li, F.; Rieske, L.K.; Sun, L.L.; Sun, S.H. Transcriptome sequencing for identification of diapause-associated genes in fall webworm, Hyphantria cunea Drury. Gene 2018, 668, 229–236. [Google Scholar] [CrossRef] [PubMed]
  26. Maaria, K.; Parker, D.J.; Mikko, M.; Salminen, T.S.; Anneli, H. Transcriptional differences between diapausing and non-diapausing D. Montana females reared under the same photoperiod and temperature. PLoS ONE 2016, 11, e0161852. [Google Scholar]
  27. Qi, X.; Zhang, L.; Han, Y.; Ren, X.; Huang, J.; Chen, H. De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes. BMC Genom. 2015, 16, 1086. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, J.; Zhao, J.; Liu, Y.H. Evaluation of endogenous reference genes in Bactrocera minax (Diptera:Tephritidae). Acta Entomol. Sin. 2014, 57, 1375–1380. [Google Scholar]
  29. Redeker, V.; Levilliers, N.; Schmitter, J.M.; Le, C.J.; Rossier, J.; Adoutte, A.; Bré, M.H. Polyglycylation of tubulin: A posttranslational modification in axonemal microtubules. Science 1994, 266, 1688–1691. [Google Scholar] [CrossRef] [PubMed]
  30. Garnham, C.P.; Yu, I.; Li, Y.; Rollmecak, A. Crystal structure of tubulin tyrosine ligase-like 3 reveals essential architectural elements unique to tubulin monoglycylases. Proc. Natl. Acad. Sci. USA 2017, 114, 6545–6550. [Google Scholar] [CrossRef] [Green Version]
  31. Ikegami, K.; Setou, M. TTLL10 can perform tubulin glycylation when co-expressed with TTLL8. FEBS Lett. 2009, 583, 1957–1963. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Rogowski, K.; Juge, F.; Dijk, J.V.; Wloga, D.; Strub, J.M.; Levilliers, N.; Thomas, D.; Bré, M.; Dorsselaer, A.V.; Gaertig, J. Evolutionary divergence of enzymatic mechanisms for posttranslational polyglycylation. Cell 2009, 137, 1076–1087. [Google Scholar] [CrossRef] [PubMed]
  33. Wloga, D.; Webster, D.M.; Rogowski, K.; Bré, M.H.; Levilliers, N.; Jerkadziadosz, M.; Janke, C.; Dougan, S.T.; Gaertig, J. TTLL3 Is a tubulin glycine ligase that regulates the assembly of cilia. Dev. Cell 2009, 16, 867–876. [Google Scholar] [CrossRef] [PubMed]
  34. Grau, M.B.; Curto, G.G.; Rocha, C.; Magiera, M.M.; Sousa, P.M.; Giordano, T.; Spassky, N.; Janke, C. Tubulin glycylases and glutamylases have distinct functions in stabilization and motility of ependymal cilia. J. Cell Biol. 2013, 202, 441–451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Pathak, N.; Austin, C.A.; Drummond, I.A. Tubulin tyrosine ligase-like genes ttll3 and ttll6 maintain zebrafish cilia structure and motility. J. Biol. Chem. 2011, 286, 11685–11695. [Google Scholar] [CrossRef] [PubMed]
  36. Feyereisen, R. Insect P450 enzymes. Annu. Rev. Entomol. 1999, 44, 507–533. [Google Scholar] [CrossRef]
  37. Scott, J.G. Cytochromes P450 and insecticide resistance. Insect Biochem. Mol. Biol. 1999, 29, 757–777. [Google Scholar] [CrossRef]
  38. Le, G.G.; Boundy, S.; Daborn, P.J.; Yen, J.L.; Sofer, L.; Lind, R.; Sabourault, C.; Madiravazzi, L.; Ffrenchconstant, R.H. Microarray analysis of cytochrome P450 mediated insecticide resistance in Drosophila. Insect Biochem. Mol. Biol. 2003, 33, 701–708. [Google Scholar]
  39. Pletcher, S.D.; Macdonald, S.J.; Marguerie, R.; Certa, U.; Stearns, S.C.; Goldstein, D.B.; Partridge, L. Genome-wide transcript profiles in aging and calorically restricted Drosophila melanogaster. Curr. Biol. 2002, 12, 712–723. [Google Scholar] [CrossRef]
  40. Riddell, C.E.; Garces, J.D.L.; Adams, S.; Barribeau, S.M.; Twell, D.; Mallon, E.B. Differential gene expression and alternative splicing in insect immune specificity. BMC Genom. 2014, 15, 1031. [Google Scholar] [CrossRef]
  41. Gottlieb, P.D.; Pierce, S.A.; Sims, R.J.; Yamagishi, H.; Weihe, E.K.; Harriss, J.V.; Maika, S.D.; Kuziel, W.A.; King, H.L.; Olson, E.N. Bop encodes a muscle-restricted protein containing MYND and SET domains and is essential for cardiac differentiation and morphogenesis. Nat. Genet. 2002, 31, 25–32. [Google Scholar] [CrossRef] [PubMed]
  42. Huang, J.; Perez-Burgos, L.; Placek, B.J.; Sengupta, R.; Richter, M.; Dorsey, J.A.; Kubicek, S.; Opravil, S.; Jenuwein, T.; Berger, S.L. Repression of p53 activity by Smyd2-mediated methylation. Nature 2006, 444, 629–632. [Google Scholar] [CrossRef] [PubMed]
  43. Mazur, P.K.; Reynoird, N.; Khatri, P. SMYD3 links lysine methylation of MAP3K2 to Ras-driven cancer. Nature 2014, 510, 283–287. [Google Scholar] [CrossRef] [PubMed]
  44. Marmorstein, R. Structure of SET domain proteins: A new twist on histone methylation. Trends Biochem. Sci. 2003, 28, 59–62. [Google Scholar] [CrossRef]
  45. Yeates, T.O. Structures of SET domain proteins: Protein lysine methyltransferases make their mark. Cell 2002, 111, 5–7. [Google Scholar] [CrossRef]
  46. Berger, S.L. The complex language of chromatin regulation during transcription. Nature 2007, 447, 407–412. [Google Scholar] [CrossRef] [PubMed]
  47. Du, S.J.; Tan, X.; Zhang, J. SMYD proteins: Key regulators in skeletal and cardiac muscle development and function. Anat. Rec. 2014, 297, 1650–1662. [Google Scholar] [CrossRef]
  48. Calpena, E.; Palau, F.; Espinós, C.; Galindo, M.I. Evolutionary history of the Smyd gene family in metazoans: A framework to identify the orthologs of human Smyd genes in Drosophila and other animal species. PLoS ONE 2015, 10, e0134106. [Google Scholar] [CrossRef]
  49. Thompson, E.C.; Travers, A.A. A Drosophila Smyd4 homologue is a muscle-specific transcriptional modulator involved in development. PLoS ONE 2008, 3, e3008. [Google Scholar] [CrossRef]
  50. Bjorum, S.M.; Simonette, R.A.; Jr, A.R.; Wang, J.E.; Lewis, B.M.; Trejo, M.H.; Hanson, K.A.; Beckingham, K.M. The Drosophila BTB domain protein Jim Lovell has roles in multiple larval and adult behaviors. PLoS ONE 2013, 8, e61270. [Google Scholar] [CrossRef]
  51. Clynen, E.; Ciudad, L.; Bellés, X.; Piulachs, M.D. Conservation of fruitless’ role as master regulator of male courtship behaviour from cockroaches to flies. Dev. Genes Evol. 2011, 221, 43–48. [Google Scholar] [CrossRef] [PubMed]
  52. Sadler, A.J.; Rossello, F.J.; Yu, L.; Deane, J.A.; Yuan, X.; Wang, D.; Irving, A.T.; Kaparakis-Liaskos, M.; Gantier, M.P.; Ying, H. BTB-ZF transcriptional regulator PLZF modifies chromatin to restrain inflammatory signaling programs. Proc. Natl. Acad. Sci. USA 2015, 112, 1535–1540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Dalton, J.E.; Fear, J.M.; Knott, S.; Baker, B.S.; Mcintyre, L.M.; Arbeitman, M.N. Male-specific Fruitless isoforms have different regulatory roles conferred by distinct zinc finger DNA binding domains. BMC Genom. 2013, 14, 659. [Google Scholar] [CrossRef] [PubMed]
  54. Demir, E.; Dickson, B.J. Fruitless splicing specifies male courtship behavior in Drosophila. Cell 2005, 121, 785–794. [Google Scholar] [CrossRef] [PubMed]
  55. Latham, K.L.; Liu, Y.S.; Taylor, B.J. A small cohort of FRUM and Engrailed-expressing neurons mediate successful copulation in Drosophila melanogaster. BMC Neurosci. 2013, 14, 57. [Google Scholar] [CrossRef] [PubMed]
  56. Boerjan, B.; Tobback, J.; De, L.A.; Schoofs, L.; Huybrechts, R. Fruitless RNAi knockdown in males interferes with copulation success in Schistocerca gregaria. Insect Biochem. Mol. Biol. 2011, 41, 340–347. [Google Scholar] [CrossRef] [PubMed]
  57. Lima, C.A.; Torquato, R.J.; Sasaki, S.D.; Justo, G.Z.; Tanaka, A.S. Biochemical characterization of a Kunitz type inhibitor similar to dendrotoxins produced by Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) hemocytes. Vet. Parasitol. 2010, 167, 279–287. [Google Scholar] [CrossRef] [PubMed]
  58. Ascenzi, P.; Bocedi, A.; Bolognesi, M.; Spallarossa, A.; Coletta, M.; De, C.R.; Menegatti, E. The bovine basic pancreatic trypsin inhibitor (Kunitz inhibitor): A milestone protein. Curr. Protein Pept. Sci. 2003, 4, 231–251. [Google Scholar] [CrossRef]
  59. Bayrhuber, M.; Vijayan, V.; Ferber, M.; Graf, R.; Korukottu, J.; Imperial, J.; Garrett, J.E.; Olivera, B.M.; Terlau, H.; Zweckstetter, M. Conkunitzin-S1 is the first member of a new Kunitz-type neurotoxin family. J. Biol. Chem. 2005, 280, 23766–23770. [Google Scholar] [CrossRef]
  60. Kurioka, A.; Yamazaki, M.; Hirano, H. Primary structure and possible functions of a trypsin inhibitor of Bombyx mori. FEBS J. 1999, 259, 120–126. [Google Scholar]
  61. Nirmala, X.; Kodrík, D.; Zurovec, M.; Sehnal, F. Insect silk contains both a Kunitz-type and a unique Kazal-type proteinase inhibitor. FEBS J. 2001, 268, 2064–2073. [Google Scholar] [CrossRef] [Green Version]
  62. Nakamura, A.; Miyado, K.; Takezawa, Y.; Ohnami, N.; Sato, M.; Ono, C.; Harada, Y.; Yoshida, K.; Kawano, N.; Kanai, S.; et al. Innate immune system still works at diapause, a physiological state of dormancy in insects. Biochem. Biophys. Res. Commun. 2011, 410, 351–357. [Google Scholar] [CrossRef] [PubMed]
  63. Jollès, P.; Jollès, J. What’s new in lysozyme research? Mol. Cell. Biochem. 1984, 63, 165–189. [Google Scholar] [CrossRef] [PubMed]
  64. Qasba, P.K.; Kumar, S.; Brew, D.K. Molecular divergence of lysozymes and α-lactalbumin. Crit. Rev. Biochem. Mol. Biol. 1997, 32, 255–306. [Google Scholar] [CrossRef] [PubMed]
  65. Elmogy, M.; Bassal, T.T.M.; Yousef, H.A.; Dorrah, M.A.; Mohamed, A.A.; Duvic, B. Isolation, characterization, kinetics, and enzymatic and nonenzymatic microbicidal activities of a novel c-type lysozyme from plasma of Schistocerca gregaria (Orthoptera: Acrididae). J. Insect Sci. 2015, 15, 57. [Google Scholar] [CrossRef] [PubMed]
  66. Tateishi, K.; Kasahara, Y.; Watanabe, K.; Hosokawa, N.; Doi, H.; Nakajima, K.; Adachi, H.; Nomoto, A. A new cell line from the fat body of Spodoptera litura (Lepidoptera, Noctuidae) and detection of lysozyme activity release upon immune stimulation. In Vitro Cell. Dev. Biol. Anim. 2015, 51, 15–18. [Google Scholar] [CrossRef] [PubMed]
Insects 10 00169 g001 550
Figure 1. Heat map of expression levels of 18 putative diapause-associated genes throughout the pupal stage of B. minax according to the reads per kb of exon model per million mapped reads (RPKM) values. Red and green color gradients indicate the high and low expression levels, respectively. PreD, pre-diapause, 1 d after pupation; ED, early-diapause, 57 d after pupation; MD, middle-diapause, 75 d after pupation; LD, late-diapause, 95 d after pupation. PD, post-diapause, 135 d after pupation.
Figure 1. Heat map of expression levels of 18 putative diapause-associated genes throughout the pupal stage of B. minax according to the reads per kb of exon model per million mapped reads (RPKM) values. Red and green color gradients indicate the high and low expression levels, respectively. PreD, pre-diapause, 1 d after pupation; ED, early-diapause, 57 d after pupation; MD, middle-diapause, 75 d after pupation; LD, late-diapause, 95 d after pupation. PD, post-diapause, 135 d after pupation.
Insects 10 00169 g001
Insects 10 00169 g002 550
Figure 2. Relative expression levels of 18 putative diapause-associated genes throughout pupal stage of B. minax. For each gene, the relative expression level was calculated against a replicate of newly-formed pupae (1 d after pupation) according to the 2−ΔΔCt method. Different letters indicate significant differences (p < 0.05).
Figure 2. Relative expression levels of 18 putative diapause-associated genes throughout pupal stage of B. minax. For each gene, the relative expression level was calculated against a replicate of newly-formed pupae (1 d after pupation) according to the 2−ΔΔCt method. Different letters indicate significant differences (p < 0.05).
Insects 10 00169 g002
Insects 10 00169 g003 550
Figure 3. Relative gene expression levels between diapause-destined and non-diapause-destined B. minax pupae measured by qRT-PCR. Asterisks indicate significant differences (p < 0.05).
Figure 3. Relative gene expression levels between diapause-destined and non-diapause-destined B. minax pupae measured by qRT-PCR. Asterisks indicate significant differences (p < 0.05).
Insects 10 00169 g003
Insects 10 00169 g004 550
Figure 4. Relative gene expression levels between diapause-destined and non-diapause-destined B. minax pupae measured by semiquantitative PCR.
Figure 4. Relative gene expression levels between diapause-destined and non-diapause-destined B. minax pupae measured by semiquantitative PCR.
Insects 10 00169 g004
Table
Table 1. Information of selected putative diapause-associated genes.
Table 1. Information of selected putative diapause-associated genes.
GeneIDAccession No.Annotation
Arrdc3BmUnigene27028.coGBEY01019016.1Arrestin domain-containing protein 3
Cyp6a9BmUnigene26209.coGBEY01017833.1Cytochrome P450 6a9
FruBmUnigene22870.coGBEY01013664.1Sex determination fruitless protein type C
GEFBmUnigene27885.c5GBEY01020590.1Ras guanine nucleotide exchange factor P
grkBmUnigene18469.coGBEY01008416.1Protein gurken
MEGF10BmUnigene26874.coGBEY01018777.1Multiple epidermal growth factor-like domains protein 10
MSTABmUnigene27964.coGBEY01020759.1Protein msta
MYBLIBmUnigene27469.coGBEY01019748.1Myb-like protein I
PcdBmUnigene23179.coGBEY01014042.1Pterin-4-alpha-carbinolamine dehydratase
setmarBmUnigene15667.coGBEY01005630.1Histone-lysine N-methyltransferase SETMAR
TTLL3BBmUnigene14579.coGBEY01004593.1Tubulin glycylase 3B
UC1BmUnigene24070.coGBEY01015150.1Uncharacterized protein
UC2BmUnigene30159.coGBEY01022960.1Uncharacterized protein
ECE-2BmUnigene27012.coGBEY01018994.1Endothelin-converting enzyme 2
KSPIBmUnigene25112.coGBEY01016431.1Kunitz-type serine protease inhibitor
LYZ1BmUnigene21628.coGBEY01012185.1Lysozyme 1
twk-7BmUnigene25271.coGBEY01016631.1TWiK family of potassium channels protein 7
wblBmUnigene21615.coGBEY01012168.1Protein windbeutel

Share and Cite

MDPI and ACS Style

Wang, J.; Fan, H.; Wang, P.; Liu, Y.-H. Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae). Insects 2019, 10, 169. https://doi.org/10.3390/insects10060169

AMA Style

Wang J, Fan H, Wang P, Liu Y-H. Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae). Insects. 2019; 10(6):169. https://doi.org/10.3390/insects10060169

Chicago/Turabian Style

Wang, Jia, Huan Fan, Pan Wang, and Ying-Hong Liu. 2019. "Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae)" Insects 10, no. 6: 169. https://doi.org/10.3390/insects10060169

APA Style

Wang, J., Fan, H., Wang, P., & Liu, Y. -H. (2019). Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae). Insects, 10(6), 169. https://doi.org/10.3390/insects10060169

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop