Next Article in Journal
Model-Based Tracking of Fruit Flies in Free Flight
Next Article in Special Issue
Identification and Validation of Reference Genes for Expression Analysis Using RT-qPCR in Leptocybe invasa Fisher and La Salle (Hymenoptera: Eulophidae)
Previous Article in Journal
Silkworm and Silk: Traditional and Innovative Applications
Previous Article in Special Issue
Analysis of Long Non-Coding RNA-Mediated Regulatory Networks of Plutella xylostella in Response to Metarhizium anisopliae Infection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis Revealed Genes Related to γ-Irradiation Induced Emergence Failure in Third-Instar Larvae of Bactrocera dorsalis

1
Guizhou Medical University (GMU), Guiyang 550000, China
2
Plant Protection Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou 510640, China
3
Guangdong Provincial Key Laboratory of High Technology for Plant Protection, Guangzhou 510640, China
4
Zanyu Technology Group Co., Ltd., Hangzhou 310000, China
5
Zhejiang Gongzheng Testing Centre Co., Ltd., Hangzhou 310000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Insects 2022, 13(11), 1017; https://doi.org/10.3390/insects13111017
Submission received: 30 August 2022 / Revised: 6 October 2022 / Accepted: 2 November 2022 / Published: 3 November 2022
(This article belongs to the Special Issue Insect Bioinformatics)

Abstract

:

Simple Summary

Bactrocera dorsalis (Diptera: Tephritidae) is a severe insect pest of numerous fruits and crops worldwide. Using the Illumina Hiseq 2000 platform, we generated a comprehensive transcriptome of B. dorsalis in response to radiation. The adult of B. dorsalis could not emerge when third-instar larvae were irradiated with 60Co-γ at 116Gy. Most differentially expressed genes (DEGs) in gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEEG) in enrichment studies were involved in hemolymph coagulation and digestive processes, such as protein digestion and absorption, while pancreatic secretion may be linked to the response to irradiation. Furthermore, the differential expressions of up/down regulated unigenes were confirmed by quantitative reverse transcription-PCR (qRT-PCR). The downregulation of sqd may contribute to the failure of emergence in irradiated samples, while the up-regulation of ENPEP and Ugt are involved in metabolic pathways, and their up-regulation may be considered the results of the irradiation stress reaction. Overall, the current study is a useful resource for the identification of potentially useful biochemical markers that can be used in quarantine.

Abstract

The oriental fruit fly is a polyphagous and highly invasive economically important pest in the world. We proposed the hypothesis that radiation treatment influence RNA expression in the larvae and leads to emergence failure. Therefore, transcriptome analyses of third-instar larvae of B. dorsalis ionizing, irradiated with 60Co-γ at 116Gy, were conducted and compared with the controls; a total of 608 DEGs were identified, including 348 up-regulated genes and 260 down-regulated ones. In addition, 130 SNPs in 125 unigenes were identified. For the DEGs, the most significantly enriched GO item was hemolymph coagulation, and some of the enriched pathways were involved in digestive processes. The subsequent validation experiment confirmed the differential expression of six genes, including sqd, ENPEP, Jhe, mth, Notch, and Ugt. Additionally, the 3401:G->A SNP in the Notch gene was also successfully validated. According to previous research, this was the first comparative transcriptome study to discover the candidate genes involved in insect molt to pupae. These results not only deepen our understanding of the emerging mechanism of B. dorsalis but also provide new insights into the research of biomarkers for quarantine insect treatment with the appropriate dose of radiation.

1. Introduction

The oriental fruit fly, Bactrocera dorsalis (Diptera: Tephritidae), belongs to the B. dorsalis complex and is a severe insect pest of several fruits and vegetables all around the world [1]. It was first detected in Sub-Saharan Africa (SSA) in costal parts of Kenya in 2003 and became widely spread and established successfully across the world [2]. It originates in tropic and sub-tropic regions and is now widespread in China, Southeast Asia, India, and Oceania. Moreover, it has also been found in the Pacific islands, including the Hawaiian Islands, French Polynesia, and Northern Marianas [3]. As an invasive agricultural pest, it causes severe losses to over 250 commercially important plants, including bananas, citrus, star fruit, guava, mango, and eggplants [4]. The female insects oviposit inside the fruit, where the larvae feed until pupation. This usually causes fruit damage and fruit drop [5,6]. Since B. dorsalis is polyphagous and highly invasive, quarantine areas for B. dorsalis have been declared in many areas [7].
Currently, strict phytosanitary treatments are used to prevent the introduction, and/or spread, of exotic tephritidae pests such as B. dorsalis. Chemical fumigants, such as ethylene dibromide (EDB) were first widely used as a pesticide [8]. However, in the late 1980s, EDB was found to be linked to cancer, and in 1984 a ban by the Environmental Protection Agency (EPA), USA, removed EDB from the chemical register for use as an insecticide [9]. Irradiation is an alternative process to chemical fumigation that is easy to apply, quick, and generally safe. Moreover, significant efforts have been made to determine the suitable radiation dose for quarantine treatment. Tephritidae quarantine treatments generally target eggs and larvae because the larvae leave the fruit to pupate in the soil [10]. In addition, the third instar (the last stage of larva) is the most radio-resistant insect stage [11]. Other third-instar tephritid fruit fly larvae have been employed for radiotolerance and substantiative studies, including B. dorsalis [12,13,14,15], Ceratitis capitata Wiedemann [13,15], Zeugodacus cucurbitae Coquillett [13,16], and Z. tau Walker [17]. The International Plant Protection Convention (IPPC) is now evaluating prospective standards for the irradiation treatment of Z. tau, B. dorsalis, and the genus Anastrepha, in order to develop international standards as an appendix to ISPM No. 28 [18] and 210, 225, and 250 Gy, where quarantine doses are currently approved for these insects, respectively. Moreover, for tephritid fruit flies and other insects, the United States Department of Agriculture Animal and Plant Health Inspection Service (USDA-APHIS) approved generic radiation treatment doses of 150 Gy and 400 Gy in 2006, respectively [19,20,21], although the adult and pupae of Lepidoptera may require higher doses. In addition, to encourage global usage irradiation of up to 150 Gy would provide quarantine security against all fruit fly species of the family Tephritidae [22]. To prevent the subsequent infection and re-infestation, the fresh commodities would be required to be sterilized with arthropod-proof packaging [16]. Phytosanitary methods, such as temperature, fumigation, and irradiation treatments are currently employed for exporting commodities. The usage of phytosanitary irradiation (PI) treatment has expanded dramatically, with a global fresh goods volume of around 30,000 tons in 2017 [12,23]. Recently, the emergence of next-generation sequencing (NGS) has provided an opportunity for researchers to efficiently obtain massive quantities of sequence data. Since the transcriptome of B. dorsalis was built by Illumina, the sequencing and transcriptome analyses of B. dorsalis focused on genes involved in reproduction, development, and insecticide resistance have also been performed. The results of previous studies suggested that Illumina transcriptome analysis is a reliable and precise way to study genomic characteristics [24]. Up to now, the comprehensive identification of biochemical markers that can be used in quarantine stations to determine whether the encountering living oriental fruit flies have received an appropriate dose of radiation remains unavailable [25].
In the present study, we present the results from the sequencing and assembly of the transcriptome of two groups (with and without the appropriate dose of radiation with three biological replicates) of the third-instar larvae of B. dorsalis. The readings were first assembled through a de novo genome assembly. DEGs were detected and functionally annotated, and the gene expression patterns for the selected DEGs were validated by qRT-PCR. Irradiated-specific single nucleotide polymorphisms (SNPs) or control-specific SNPs shared by all samples in the group were also identified, and some of them were validated by PCR and Sanger sequencing. This transcriptome is a useful resource for the identification of potentially useful biochemical markers that can be used in quarantine. Furthermore, our data will provide insights into the molecular mechanisms putatively related to the irradiation response of B. dorsalis.

2. Materials and Methods

2.1. Insect Rearing and Sample Preparation

The laboratory colony of B. dorsalis was from a strain originally collected and reared following the method of Zheng et al. [26] in Guangzhou, Guangdong, China. Adults were fed with water and a standard laboratory diet consisting of yeast powder, honey, sugar, and ascorbic acid. All specimens were kept in an incubator at 27 ± 1 °C, 70 ± 5%, at a relative humidity, and a photoperiod of 12:10(L:D) h. Third-instar larvae were collected and divided into two groups: one group was treated and irradiated with 60Co-γ at 116 Gy, according to the Guangzhou Huada Biotechnology Co., Ltd. set manufactured protocol, while the other group was used as controls. Three biological replicates were used for each group: CK (Ck1A, Ck2A, and CK3A) and R (R1A, R2A, and R3A). The Fricke system was employed for both the reference standard and routine dosimetry [27]. This dosimetry system was constructed in accordance with ASTM E1026-13 [27] and ISO/ASTM 51261 [27], and the measured uncertainty value was computed in accordance with ISO/ASTM 51707 [28]. After irradiation, the development from larvae to the pupa of B. dorsalis was delayed for about 3 days, and no adults emerged.

2.2. RNA Isolation and Sequencing

RNA from 18 flies (an irradiated group and control group with three replicates of each group, including 3 flies per replicate) was extracted using a Trizol-reagent (Invitrogen, Carlsbad, CA, USA), following the manufacturer’s procedures. The quality of the total RNA was verified by measuring its absorbance at 260 nm using a NanoVue UV-Vis spectrophotometer (GE Healthcare Bio-Science, Uppsala, Sweden). The integrity of RNA was verified by 1% agarose gel electrophoresis. Poly (A) mRNA was purified using oligo-dT beads. It was fragmented into short fragments and served as a template to synthesize the first-strand cDNA using a random hexamer-primed reverse transcription, followed by the synthesis of the second-strand cDNA. The cDNA libraries were then prepared following Illumina’s protocols [29] and were sequenced on the Illumina Hiseq 2000 platform with 100 bp paired-end reads in BGI (Beijing Genomics Institute, Shenzhen, Guangdong, China).

2.3. De Novo Assembly

Raw data generated by the sequencing were first filtered to remove low-quality reads. Briefly, the following reads were excluded: (1) reads containing adapter sequences; (2) reads with over 5% uncertain nucleotides (Ns); (3) reads with over 10% low-quality bases (Q ≤ 10). The remaining clean reads were used for transcriptome de novo assembly using the Trinity de novo RNAseq assembly pipeline (Inchworm, Chrysalis, and Butterfly) with k-mer sizes set at 25. The output of the Trinity pipeline is a fasta formatted file containing full-length transcripts, erroneous contigs, partial transcript fragments, and non-coding RNA molecules. By using RSEM (RNA-Seq by expectation maximization) software1.2.0 [30], these sequences were then filtered to identify contigs containing full or near full-length transcripts or likely coding regions and isoforms that were represented at a minimum level based on the read abundance. The resulting de novo transcripts (unigenes) were screened against protein databases, such as NR, Swiss-Prot, KEGG, and COG using BLASTX with a cut-off E value of 10−5. For each transcript, the best aligning result was used to determine the sequence direction. If different databases gave conflicting results, a priority order of GenBank non-redundant (NR), Swiss-Prot, KEGG, and COG was followed. If a unigene did not map to any of the entries in the above databases, the software ESTScan [31] was introduced to predict its coding regions and sequence direction.

2.4. Expression Analysis

Clean reads from each sequencing library were independently mapped to the above transcriptome assembly (only contigs containing full and partial ORFs) by using Soapaligner/soap2 [32]. Mismatches of less than three bases were allowed. Gene expression levels were calculated using the RPKM method [33]. If there was more than one transcript for a given gene, the longest transcript remained to calculate its expression level and coverage. The false discovery rate (FDR) method [34] was used to identify differentially expressed genes (DEGs) between the two groups. Genes with FDR < 0.001 and an absolute value of log2Ratio ≥ 1 (fold change ≥ 2) were selected as DEGs in the DESeq package in R (1.18.0) [35]. Then, the DEGs were further annotated by GO function analysis and KEGG pathway analysis [36].

2.5. GO and KEGG Enrichment Analysis

All genes were mapped to the KEGG pathways (http://www.genome.jp/kegg/ (accessed on 1 August 2020)) database and GO database (http://geneontology.org/ (accessed on 1 August 2020)). A hypergeometric distribution test was carried out to identify pathways or GO items overrepresented with DEGs. Multiple testing corrections were performed with the Benjamini-Hochberg procedure [37]. KEEG and GO categories less than p < 0.05 were considered highly enriched.

2.6. Heterozygous SNP Analysis

The general unigene sequences were used as references for SNP identification. The SNP calling was carried out by using SoapSNP. SoapSNP calculates the likelihood of each genotype based on the alignment results and the corresponding sequencing quality scores. It infers the genotype with the highest posterior probability based on Bayes’ theorem. The SNPs from the two groups were classified into three categories: (1) the common SNPs shared by the irradiated and control samples and the SNPs that were heterozygous in each group; (2) the irradiated-specific SNPs, which were derived only from the irradiated samples and the controls that had the same genotype as the reference; and (3) the control-specific SNPs, which were derived from only the controls and irradiated samples which had the same genotype as the reference. Only irradiated-specific SNPs or control-specific SNPs shared by all the samples in its group remained.

2.7. Validation of Gene Expression by qRT-PCR

The differential expressions for the selected genes were validated in two groups (the control and treatment groups), and these genes were picked from the transcriptome data based on their differential expression trends in both groups. In total, we used 18 individuals (from the control and treatment groups) with three biological replicates and three technical replicates. The samples were treated with the same procedure as the samples used for sequencing. RNA was extracted as previously described for sequencing. A total of 2 µg of RNA from each sample was reverse transcribed in a 20 mL reaction system using the PrimerScript TM RT Reagent Kit (Takara Bio Inc., Shiga, Japan). The PCR conditions for both genes were 95 °C for 3 min, followed by 34 cycles of 94 °C for 30 s, 60 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 10 min. The relative fold changes were determined using the 2−ΔΔCt method with the a-tubulin gene (GU269902) of B. dorsalis as an internal control. The GraphPad Prism was used to analyze the data (version 7; GraphPad Software, San Diego, CA, USA, https://www.graphpad.com/scientific-software/prism/ (accessed on 27 August 2018)). For mean comparisons, Student’s t-test (p > 0.05) was employed. All data were displayed as bar charts, along with the means and standard deviations (Mean ± SD).

2.8. SNP Validation

To validate the potential SNP markers, the SNP of the Notch gene was further selected and validated into two groups (the control and treatment groups) by PCR and Sanger sequencing [38]. SNP was validated six times, with three biological replicates and three technical replicates, and each replicate had 3 individuals. The samples were treated with the same procedure as the samples used for sequencing. DNA was extracted for the SNP validation experiment. Primers were designed by the Primer Premier 5 software (PREMIER Biosoft International, Palo Alto, CA, USA). PCR amplification in a 50 μL reaction was performed following the conditions: 95 °C for 2 min, 35 cycles of 95 °C for 15 s, 60 °C for 20 s, 72 °C for 30 s, and 72 °C for 2 min. The PCR product purification was carried out using the E.Z.N.A.® Gel Extraction Kit (Omega Bio-tek, Inc., Norcross, GA, USA). Sanger sequencing was performed in both forward and reverse directions on an ABI 3730 DNA Analyzer. The sequence trace files were analyzed manually.

3. Results

3.1. Sequencing and Transcriptome Assembly

The data generated by Illumina’s Hiseq 2000 sequencing contained a total of 31,987,764,800 bases from 159,938,824 pair-end sequence reads. The average number of reads from the irradiated group and controls were 26,898,794 and 26,414,148, respectively (Table S1). For each sample, the raw data were firstly assembled into contigs separately (Table S2), and then the unigenes were combined into the final assembly (Table 1). In total, we generated 50,268 unigenes with the mean length of 781 bp. As shown in Table 1, the unigene size distribution revealed the following: 59.63% (29,975) were between 100 and 500 bp in length; 19.59% (9845) were between 500 and 1000 bp in length; 12.67% (6370) were between 1000 and 2000 bp; and 8.11% (1367) were more than 2000 bp.

3.2. Functional Annotation

BLASTX alignments between the predicted unigenes and protein databases, including NR, NT, Swiss-Prot, KEGG, COG, and GO, revealed that 28,358 (56.41%) could be annotated with the known biological functions (Table S3). The majority of unigenes (79.09%) had a strong homology with Drosophila (Figure 1). Among these, 13.28% were best matched to the sequences from Drosophila melanogaster, followed by Drosophila virilis (11.78%) and Drosophila willistoni (11.31%). Other unigenes (20.91%) also had hits with other insect species, such as Tribolium castaneum (0.81%) and Bombyx mori (0.17%). Compared to species within Diptera, 1.17% of the unigenes matched sequences from Aedes aegypti, 0.20% from Glossina morsitans morsitans, and 0.14% from Musca domestica.
To classify the functions of the unigenes, the COG database was used. As shown in Figure 2, a total of 21,313 unigenes (42.40%) were annotated and divided into 26 specific categories. The general function category, which contained 16.23% of the mapped unigenes (3459), was the top category, followed by the transcription (2017, 9.46%), carbohydrate transport and metabolism (1683, 7.90%), replication, recombination, and repair (1586, 7.44%), and translation, ribosomal structure, and biogenesis (1515, 7.11%). The nuclear structure was the smallest category, which only contained two unigenes (0.0094%). GO analysis divided the unigenes into three ontologies: molecular function, cellular component, and biological process. As shown in Figure 3, 32.65% (16,412 unigenes) of the unigenes were categorized into 59 groups. The cellular process was the group containing the most unigenes.
Pathway mapping analysis was performed based on the KEGG pathway database. As shown in Table S4, 16,528 unigenes were mapped to 258 pathways. Compared with other, metabolic pathways contained the most unigenes (2421, 14.65%), followed by pathways in cancer (744, 4.5%), focal adhesion (588, 3.56%), and purine metabolism (583, 3.53%).

3.3. Differentially Expressed Genes in Response to Irradiation and Pathway Enrichment Analysis of DEGs

To identify the genes that responded to irradiation, differentially expressed sequences between the two groups were identified (Table S5). There were 608 significantly differentially expressed unigenes detected between the irradiated and control samples, including 348 up-regulated and 260 down-regulated genes. To understand the functions of DEGs, DEGs were mapped to the GO database, and enrichment analysis was performed. A total of 252 DEGs were mapped to the GO database, and 78 GO items were significantly enriched (Table 2). Hemolymph coagulation was the most significantly enriched GO item (p = 4.10 × 10−25), followed by hemostasis (p = 4.83 × 10−22), coagulation (p = 4.83 × 10−22), and the regulation of body fluid levels (p = 1.85 × 10−21). All of the DEGs were also mapped to the KEGG database and compared with the whole transcriptome background. Among all of the genes with KEGG annotation, a total of 380 DEGs were found between the two groups. A total of 32 pathways were notably enriched (Table 3). The most significant pathway was the renin-angiotensin system (p = 5.42 × 10−10), followed by protein digestion and absorption (p = 3.6 × 10−8) and pancreatic secretion (p = 3.6 × 10−8). Most of the enriched pathways (18/32) were metabolic pathways, such as porphyrin and chlorophyll metabolism, retinol metabolism, steroid hormone biosynthesis, and insect hormone biosynthesis.

3.4. Verification of Differentially Expressed Genes

To further evaluate the sequencing results, the expression levels of six genes (eight unigenes) potentially involved in the response to irradiation were analyzed by qRT-PCR. As shown in Table 4 and Figure 4, the validation experiment revealed the same expression tendency as the sequencing data.

3.5. SNP Validation

Among the three pair samples for the SNP validation, one pair of the samples failed the PCR experiment, and only another two pairs of the samples were successfully sequenced. Among them, the 3401:G->A SNP in the Notch gene was successfully validated (Figure 5), implicating its involvement in response to irradiation.

4. Discussion

Irradiation, temperature, and fumigation are currently the most popular phytosanitary treatment techniques used for the disinfestation of postharvest fruit flies in fruits and vegetables [39]. Irradiation has gained popularity as a phytosanitary treatment in the last decade for the eradication of arthropods, ornamentals, and preserved products [40]. Irradiation can cause sub-lethal molecular or biochemical changes, triggering a cascade of physiological changes. In order to eliminate insect pests, gamma radiation is frequently used as a potentially better alternative to toxic fumigants [41].
In this study, with transcriptome sequencing data from two groups (the control and treatment groups with their three biological replicates), we selected molecular markers for the identification of insects that had been treated with gamma irradiation. Previous researchers demonstrated that gamma radiation is environmentally beneficial, with no noticeable environmental adverse effects. Moreover, gamma irradiation prevents the reproduction of several stored grain pest species [42,43].
Our deep sequencing strategy (over 5 G of data for each sample) generated a transcriptome assembly with 50,268 unigenes (a mean length of 781 bp). A total of 20.78% of the unigenes were more than 1000 bp. The previously reported B. dorsalis transcriptome, which was also sequenced on the Illumina platform, generated 49,804 unigenes with a mean size of 456 bp, and less than 10% of the unigenes were over 1000 bp. With 454-pyrosequencing, Zheng et al. [44] generated a B. dorsalis transcriptome with 48,876 unigenes (mean size of 750 bp), and about 18.4% of the unigenes were over 1000 bp. Therefore, compared with the previously reported B. dorsalis transcriptome assembly, we generated more detailed data that facilitated our subsequent analysis. The functional annotation results of the unigenes were generally consistent with the previous transcriptome: the majority of the unigenes had a strong homology with Drosophila; the general function category was the top COG category, and the KEGG metabolic pathways contained the most unigenes.
A total of 608 significantly differentially expressed unigenes between the two groups may be attributed to their response to γ irradiation. GO item enrichment analyses revealed that the most significantly enriched GO item was hemolymph coagulation, which is one of the first responses to injury in insects. KEGG pathway analyses of these genes revealed that some of the pathways were involved in digestive processes, such as protein digestion and absorption and pancreatic secretion. This observation is consistent with the well-established irradiation effects on digestive physiology in other insect pests.
The differential expressions of eight unigenes were further confirmed by qRT-PCR. As shown in Table 4 (Figure 4), the squid gene sqd was downregulated in irradiated samples while ENPEP, Jhe, mth, Notch, and Ugt were highly expressed. The downregulation of sqd may contribute to the failure of emergence in irradiated samples since the initial organization of the Drosophila dorsoventral axis depends on the proteins encoded by sqd. Both ENPEP and Ugt are involved in metabolic pathways, and their up-regulation may be considered the results of the γ irradiation stress reaction. Juvenile hormone esterase (Jhe) plays an important role in the regulation of the juvenile hormone (JH). In lepidopteran hemolymph, decreased JH titers are positively correlated with increased JHE abundance [45]. The high expression of Jhe here may result in the low titers of JH, which are associated with diapause. For example, Monarch adult diapause is instituted by a deficiency in JH production [46]. Therefore, the upregulation of Jhe may contribute to the diapauses of our γ irradiated samples. Individuals with a mutated mth displayed an approximately 35 percent increase in their average lifespan and showed significant resistance to oxidative stress, starvation, and heat stress. We also detected a heterozygous SNP in this gene in the γ irradiated samples, which might be caused by the stress response. Notch (N) signaling is involved in establishing the correct width of various wing veins [47], and the up-regulation of Notch here may also contribute to the failure of emergence in the irradiated samples since the gain of function in Notch showed thinner and incomplete wing veins. In addition to the confirmed up-regulation of Notch, an SNP in this gene was also validated by the Sanger sequencing. It is highly possible that the confirmed SNP was induced by the irradiation treatment. Since the SNP was validated six times in three biological and three technical replicates of the samples, it could be used as DNA markers in quarantine. Other SNPs identified in this study are also useful sources for further studies on selecting DNA markers in quarantine.

5. Conclusions

In conclusion, we have generated a comprehensive transcriptome of the B. dorsalis in response to radiation using the Illumina Hiseq 2000 platform. The genes differentially expressed in the controls and irradiated samples were largely identified and functionally annotated. Our results provide useful biochemical markers that could be used in irradiation quarantine treatment. Furthermore, our data provide insights into the molecular mechanisms putatively related to the irradiation response of the third-instar larvae of B. dorsalis.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/insects13111017/s1, Table S1: Statistics of the transcriptome sequencing data for each sample; Table S2: Statistics of assembled contigs; Table S3: Functional annotation statistics; Table S4: Pathway mapping analyses for all unigenes; Table S5: Differentially expressed sequences between two groups; Table S6: Genes with irradiated-specific SNPs or control-specific SNPs.

Author Contributions

C.S., S.S., W.W. and Q.L. conceived and designed the research; C.G. and Q.L. performed the experiments; Q.L. and S.S. analyzed the data. C.S., S.S., W.W., C.G., Y.G. and Q.L. wrote and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by Natural Science Foundation of Guangdong Province, China (2021A 1515012385); Discipline team building projects of Guangdong Academy of Agricultural Sciences in the 14th Five-Year Period (202105TD);The Project of Collaborative Innovation Center of GDAAS-XTXM202202 (XT202211).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available in the Supplementary Material.

Acknowledgments

This study was supported by the ninth group of the Young Scientist Foundation of Guizhou Province, China (No.40 [2013]); the Science and Technology Foundation of Guizhou Province, China (No. J [2011]2345); and the International Cooperation Foundation of Guizhou Province, China (No.G [2012]7030); the Science and Technology Planning Protection of Guangdong Province (No.2014B070706017). We would like to express our sincere gratitude to Mohsin Nawaz for his scientific support and all types of extended help to this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mutamiswa, R.; Nyamukondiwa, C.; Chikowore, G.; Chidawanyika, F. Overview of oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) in Africa: From invasion, bio-ecology to sustainable management. Crop Prot. 2021, 141, 105492. [Google Scholar] [CrossRef]
  2. Lux, S.A.; Copeland, R.S.; White, I.M.; Manrakhan, A.; Billah, M.K. A new invasive fruit fly species from the Bactrocera dorsalis (Hendel) group detected in East Africa. Int. J. Trop. Insect Sci. 2003, 23, 355–361. [Google Scholar] [CrossRef]
  3. Leblanc, L.; Vueti, E.T.; Drew, R.A.; Allwood, A.J. Host Plant Records for Fruit Flies (Diptera: Tephritidae: Dacini) in the Pacific Islands. Proc. Hawaiian Entomol. Soc. 2012, 44, 11–53. [Google Scholar]
  4. Allwood, A.J.; Chinajariyawong, A.; Kritsaneepaiboon, S.; Drew, R.A.I.; Hamacek, E.L.; Hancock, D.L.; Hengsawad, C.; Jipanin, J.C.; Jirasurat, M.; Krong, C.K.; et al. Host plant records for fruit flies (Diptera: Tephritidae) in Southeast Asia. Raffles Bull. Zool. 1999, 47, 1–92. [Google Scholar]
  5. Fletcher, B.S. The Biology of Dacine Fruit Flies. Annu. Rev. Entomol. 1987, 32, 115–144. [Google Scholar] [CrossRef]
  6. Salmah, M.; Nurul Fatihah, M.; Hailmi, M.; Norhayati, N. Growth and Development of oriental fruit fly, Bactrocera dorsalis Hendel (Diptera: Tephritidae) reared on sweet potatoes (Ipomoea batatas L.) based artificial diet. Serangga 2019, 25, 96–107. [Google Scholar]
  7. Ocitti, P.; Ndlela, S.; Akol, A.; Muyinza, H.; Mohamed, S. Non-chemical post-harvest disinfestation of Bactrocera dorsalis (Hendel)(Diptera: Tephritidae) in Tommy Atkins mango using hot-water immersion treatment. Afr. Entomol. 2021, 29, 238–247. [Google Scholar] [CrossRef]
  8. Arthur, V. Disinfestation of mangoes haden in-fested by Anastrepha fraterculus (WIED.; 1830) (Diptera, Tephritidae), with gamma radiation. Braz. J. Radiat. Sci. 2021, 9. [Google Scholar] [CrossRef]
  9. Anonymous. Ethylene dibromide. Fed. Regist. 1983, 48, 4263. [Google Scholar]
  10. Ladaniya, M. Citrus Fruit: Biology, Technology and Evaluation; Academic Press: Cambridge, MA, USA, 2010. [Google Scholar]
  11. Hallman, G.J. Expanding radiation quarantine treatments beyond fruit flies. Agric. For. Entomol. 2000, 2, 85–95. [Google Scholar] [CrossRef]
  12. Zhao, J.; Ma, J.; Wu, M.; Jiao, X.; Wang, Z.; Liang, F.; Zhan, G. Gamma radiation as a phytosanitary treatment against larvae and pupae of Bactrocera dorsalis (Diptera: Tephritidae) in guava fruits. Food Control 2017, 72, 360–366. [Google Scholar] [CrossRef]
  13. Follett, P.A.; Armstrong, J.W. Revised irradiation doses to control melon fly, Mediterranean fruit fly, and oriental fruit fly (Diptera: Tephritidae) and a generic dose for tephritid fruit flies. J. Econ. Entomol. 2004, 97, 1254–1262. [Google Scholar] [CrossRef] [PubMed]
  14. Srimartpirom, M.; Burikam, I.; Limohpasmanee, W.; Kongratarporn, T.; Thannarin, T.; Bunsiri, A.; Follett, P.A. Low-dose irradiation with modified atmosphere packaging for mango against the oriental fruit fly (Diptera: Tephritidae). J. Econ. Entomol. 2018, 111, 135–140. [Google Scholar] [CrossRef] [PubMed]
  15. Dias, V.S.; Hallman, G.J.; Martínez-Barrera, O.Y.; Hurtado, N.V.; Cardoso, A.A.; Parker, A.G.; Caravantes, L.A.; Rivera, C.; Araújo, A.S.; Maxwell, F.; et al. Modified atmosphere does not reduce the efficacy of phytosanitary irradiation doses recommended for tephritid fruit flies. Insects 2020, 11, 371. [Google Scholar] [CrossRef] [PubMed]
  16. Follett, P.A.; Wall, M.; Bailey, W. Influence of modified atmosphere packaging on radiation tolerance in the phytosanitary pest melon fly (Diptera: Tephritidae). J. Econ. Entomol. 2013, 106, 2020–2026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Guoping, Z.; Lili, R.; Ying, S.; Qiaoling, W.; Daojian, Y.; Yuejin, W.; Tianxiu, L. Gamma irradiation as a phytosanitary treatment of Bactrocera tau (Diptera: Tephritidae) in pumpkin fruits. J. Econ. Entomol. 2015, 108, 88–94. [Google Scholar] [CrossRef]
  18. IIPP. List of Topics for IPPC Standards; 2020. Available online: https://www.ippc.int/en/core-activities/standards-setting/list-topics-ippc-standards/ (accessed on 15 August 2022).
  19. Follett, P.A. Generic radiation quarantine treatments: The next steps. J. Econ. Entomol. 2009, 102, 1399–1406. [Google Scholar] [CrossRef]
  20. Follett, P.A.; Weinert, E.D. Phytosanitary irradiation of fresh tropical commodities in Hawaii: Generic treatments, commercial adoption, and current issues. Radiat. Phys. Chem. 2012, 81, 1064–1067. [Google Scholar] [CrossRef]
  21. UAUSDoAAaPHI. Treatments for Fruits and Vegetables. Fed. Regist. 2006. Available online: https://www.federalregister.gov/documents/2006/01/27/06-746/treatments-for-fruits-and-vegetables (accessed on 15 August 2022).
  22. IPP. International Standards for Phytosanitary Measures (ISPM) No. 28, Phytosanitary Treatments for Regulated Pests; Food and Agricultural Organization Rome: Rome, Italy, 2009. [Google Scholar]
  23. Fang, Y.; Kang, F.; Zhan, G.; Ma, C.; Li, Y.; Wang, L.; Wei, Y.; Gao, X.; Li, Z.; Wang, Y. The effects of a cold disinfestation on Bactrocera dorsalis survival and navel orange quality. Insects 2019, 10, 452. [Google Scholar] [CrossRef] [Green Version]
  24. Jiang, F.; Liang, L.; Wang, J.; Zhu, S. Chromosome-level genome assembly of Bactrocera dorsalis reveals its adaptation and invasion mechanisms. Commun. Biol. 2022, 5, 25. [Google Scholar] [CrossRef]
  25. Chen, X.; Lei, Y.; Li, H.; Xu, L.; Yang, H.; Wang, J.; Jiang, H. CRISPR/Cas9 mutagenesis abolishes odorant-binding protein BdorOBP56f-2 and impairs the perception of methyl eugenol in Bactrocera dorsalis (Hendel). Insect Biochem. Mol. Biol. 2021, 139, 103656. [Google Scholar] [CrossRef] [PubMed]
  26. Zheng, Y.; Song, Z.-W.; Zhang, Y.-P.; Li, D.-S. Ability of Spalangia endius (Hymenoptera: Pteromalidae) to Parasitize Bactrocera dorsalis (Diptera: Tephritidae) after Switching Hosts. Insects 2021, 12, 613. [Google Scholar] [CrossRef] [PubMed]
  27. 51026 IAE; Practice for Using the Frick Dosimetry System. ASTM International: West Conshohocken, PA, USA, 2015.
  28. ISO/ASTM 51707:2002; Standard Guide For Estimating Uncertainties in Dosimetry For Radiation Processing. Available online: https://standards.iteh.ai/catalog/standards/iso/f7a12209-e259-4e42-9245-9e15d7d39daa/iso-astm-51707-2002 (accessed on 15 August 2022).
  29. Shen, G.M.; Dou, W.; Niu, J.Z.; Jiang, H.B.; Yang, W.J.; Jia, F.X.; Hu, F.; Cong, L.; Wang, J.J. Transcriptome analysis of the oriental fruit fly (Bactrocera dorsalis). PLoS ONE 2011, 6, e29127. [Google Scholar] [CrossRef] [PubMed]
  30. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed]
  31. Iseli, C.; Jongeneel, C.V.; Bucher, P. ESTScan: A program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1999, 138–148. [Google Scholar]
  32. Li, R.; Yu, C.; Li, Y.; Lam, T.W.; Yiu, S.M.; Kristiansen, K.; Wang, J. SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics 2009, 25, 1966–1967. [Google Scholar] [CrossRef] [Green Version]
  33. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef]
  34. Audic, S.; Claverie, J.M. The significance of digital gene expression profiles. Genome Res. 1997, 7, 986–995. [Google Scholar] [CrossRef]
  35. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef] [Green Version]
  36. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  37. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  38. Fan, B.; Xie, D.; Li, Y.; Wang, X.; Qi, X.; Li, S.; Meng, Z.; Chen, X.; Peng, J.; Yang, Y.; et al. A single intronic single nucleotide polymorphism in splicing site of steroidogenic enzyme hsd17b1 is associated with phenotypic sex in oyster pompano, Trachinotus anak. Proc. R. Soc. B 2021, 288, 20212245. [Google Scholar] [CrossRef] [PubMed]
  39. Heather, N.W.; Hallman, G.J. Pest Management and Phytosanitary Trade Barriers; CABI: Wallingford, UK, 2008. [Google Scholar]
  40. Joint, F. Irradiation as a phytosanitary treatment of food and agricultural commodities. In Proceedings of A Final Research Coordination Meeting; Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture: Vienna, Austria, 2004. [Google Scholar]
  41. Kiran, R.; Shenoy, K.B.; Venkatesha, M.G. Effect of gamma radiation as a post-harvest disinfestation treatment against life stages of the coffee berry borer, Hypothenemus hampei (Ferrari)(Coleoptera: Curculionidae). Int. J. Radiat. Biol. 2019, 95, 1301–1308. [Google Scholar] [CrossRef] [PubMed]
  42. Seal, D.R.; Tilton, E.W. Effect of gamma radiation on the metamorphic stages of Dermestes maculatus DeGeer (Coleoptera: Dermestidae). Int. J. Radiat. Appl. Instrum. Part A Appl. Radiat. Isot. 1986, 37, 531–535. [Google Scholar] [CrossRef]
  43. Mastrangelo, T.; Walder, J. Use of radiation and isotopes in insects. In Radioisotopes–Applications in Bio-Medical Science; IntechOpen: London, UK, 2011; pp. 67–92. [Google Scholar]
  44. Zheng, W.; Peng, T.; He, W.; Zhang, H. High-throughput sequencing to reveal genes involved in reproduction and development in Bactrocera dorsalis (Diptera: Tephritidae). PLoS ONE 2012, 7, e36463. [Google Scholar] [CrossRef]
  45. Kerkut, G.A.; Gilbert, L.I. Comprehensive Insect Physiology, Biochemistry, and Pharmacology; Pergamon Press: New York, NY, USA, 1985; pp. 431–472. [Google Scholar]
  46. Herman, W.S. Studies on the Adult Reproductive Diapause of the Monarch Butterfly, Danaus plexippus. Biol. Bull. 1981, 160, 89–106. [Google Scholar] [CrossRef]
  47. de Celis, J.; García-Bellido, A. Roles of the Notch gene in Drosophila wing morphogenesis. Mech. Dev. 1994, 46, 109–122. [Google Scholar] [CrossRef]
Figure 1. Species distribution of the BLASTX matches (cutoff value E < 10−5) of the transcriptome unigenes.
Figure 1. Species distribution of the BLASTX matches (cutoff value E < 10−5) of the transcriptome unigenes.
Insects 13 01017 g001
Figure 2. Classification of the clusters of orthologous groups (COG) for the transcriptome.
Figure 2. Classification of the clusters of orthologous groups (COG) for the transcriptome.
Insects 13 01017 g002
Figure 3. Classification of the gene ontology (GO) for the transcriptome.
Figure 3. Classification of the gene ontology (GO) for the transcriptome.
Insects 13 01017 g003
Figure 4. qRT-PCR confirmation of the differentially expressed genes between irradiated samples and controls.
Figure 4. qRT-PCR confirmation of the differentially expressed genes between irradiated samples and controls.
Insects 13 01017 g004
Figure 5. Sanger sequencing results for the SNP identified in the Notch gene.
Figure 5. Sanger sequencing results for the SNP identified in the Notch gene.
Insects 13 01017 g005
Table 1. Statistics of assembled unigenes.
Table 1. Statistics of assembled unigenes.
SampleUnigene CountsN50Mean LengthAll Unigene CountsLength of All Unigene
100–500 bp500–1000 bp1000–1500 bp1500–2000 bp≥2000 bp
CK1A42,012752524811201141966846654,63825,475,845
CK2A37,42164822088986111665646748,09322,476,031
CK3A38,296778028341370179082452552,07027,338,230
R1A32,5635640168377683160845741,49318,942,224
R2A34,811680323121110154284453046,57824,677,137
R3A36,509781628641630224897257151,06729,140,744
All29,9759845400523654078136778150,26839,256,184
Note: Irradiated samples are shown with the prefix of “R” while control samples are shown with the prefix of “CK”.
Table 2. Gene ontology term enrichment analyses results for the differentially expressed genes.
Table 2. Gene ontology term enrichment analyses results for the differentially expressed genes.
CategoryGene Ontology TermAdjusted p-Value
biological processhemolymph coagulation 4.10× 10−25
biological processhemostasis 4.83 × 10−22
biological processcoagulation 4.83 × 10−22
biological processregulation of body fluid levels 1.85 × 10−21
molecular functionstructural constituent of cuticle 1.30 × 10−16
molecular functionstructural constituent of chitin-based larval cuticle 5.55 × 10−16
cellular componentextracellular region 1.35 × 10−15
biological processwound healing 1.70 × 10−15
molecular functionstructural constituent of chitin-based cuticle 3.22 × 10−14
biological processresponse to wounding 1.42 × 10−12
biological processhumoral immune response 5.87 × 10−9
biological processdefense response 8.80 × 10−9
biological processinnate immune response 1.81 × 10−8
molecular functionexopeptidase activity 6.87 × 10−7
molecular functioncatechol oxidase activity 5.89 × 10−6
molecular functionL-DOPA monooxygenase activity 5.89 × 10−6
molecular functiondopamine monooxygenase activity 5.89 × 10−6
biological processcellular lipid catabolic process 1.60 × 10−5
molecular functionpeptidase activity 1.73 × 10−5
cellular componentintracellular ferritin complex 4.74 × 10−5
cellular componentferritin complex 4.74 × 10−5
molecular functionferrous iron binding 6.24 × 10−5
molecular functionpeptidase activity, acting on L-amino acid peptides 7.77 × 10−5
molecular functionoxidoreductase activity, acting on diphenols and related substances as donors, oxygen as acceptor 7.81 × 10−5
biological processregulation of mitochondrial translation 1.10 × 10−4
cellular componentextracellular space 1.60 × 10−4
biological processtriglyceride catabolic process 2.40 × 10−4
biological processneutral lipid catabolic process 2.40 × 10−4
biological processacylglycerol catabolic process 2.40 × 10−4
biological processglycerolipid catabolic process 2.40 × 10−4
biological processtriglyceride metabolic process 2.70 × 10−4
molecular functionaminopeptidase activity 3.50 × 10−4
biological processdopamine metabolic process 5.60 × 10−4
biological processsesquiterpenoid metabolic process 5.70 × 10−4
biological processjuvenile hormone metabolic process 5.70 × 10−4
biological processjuvenile hormone catabolic process 5.70 × 10−4
biological processisoprenoid catabolic process 5.70 × 10−4
biological processsesquiterpenoid catabolic process 5.70 × 10−4
biological processterpenoid catabolic process 5.70 × 10−4
biological processnegative regulation of translational initiation 5.70 × 10−4
molecular functionepoxide hydrolase activity 9.00 × 10−4
biological processlipid catabolic process 9.40 × 10−4
biological processimmune response 1.08 × 10−3
biological processterpenoid metabolic process 1.71 × 10−3
biological processhormone catabolic process 1.71 × 10−3
biological processcatecholamine metabolic process 1.90 × 10−3
biological processcatechol-containing compound metabolic process 1.90 × 10−3
biological processphenol-containing compound metabolic process 1.90 × 10−3
biological processdiol metabolic process 1.90 × 10−3
molecular functionether hydrolase activity 2.07 × 10−3
molecular functionoxidoreductase activity, acting on diphenols and related substances as donors 2.90 × 10−3
molecular functioncarboxypeptidase activity 3.43 × 10−3
biological processthioester biosynthetic process 3.68 × 10−3
biological processacyl-CoA biosynthetic process 3.68 × 10−3
biological processneutral lipid metabolic process 4.01 × 10−3
biological processacylglycerol metabolic process 4.01 × 10−3
biological processimmune system process 4.49 × 10−3
biological processnegative regulation of protein metabolic process 4.73 × 10−3
molecular functionoxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen 7.16 × 10−3
biological processmonocarboxylic acid metabolic process 7.34 × 10−3
molecular functionoxidoreductase activity 9.03 × 10−3
biological processsingle-organism metabolic process 9.08 × 10−3
molecular function4-aminobutyrate transaminase activity 1.11 × 10−2
molecular functionjuvenile hormone epoxide hydrolase activity 1.11 × 10−2
biological processlipid metabolic process 1.13 × 10−2
molecular functionserine-type peptidase activity 1.27 × 10−2
biological processregulation of mitochondrion organization 1.38 × 10−2
biological processacyl-CoA metabolic process 1.49 × 10−2
biological processthioester metabolic process 1.49 × 10−2
molecular functionserine hydrolase activity 1.61 × 10−2
molecular functioneukaryotic initiation factor 4E binding 1.82 × 10−2
biological processdsRNA transport 2.16 × 10−2
biological processregulation of hormone metabolic process 2.61 × 10−2
molecular functiontransferase activity, transferring hexosyl groups 2.82 × 10−2
molecular functionelectron carrier activity 3.32 × 10−2
biological processfatty acid catabolic process 3.38 × 10−2
cellular componentvacuolar proton-transporting V-type ATPase, V1 domain 4.63 × 10−2
molecular functioncatalytic activity 4.79 × 10−2
Table 3. KEGG pathway enrichment analyses results for the differentially expressed genes.
Table 3. KEGG pathway enrichment analyses results for the differentially expressed genes.
Pathway IDPathwayAdjusted p ValueLevel 1Level 2
ko04614Renin-angiotensin system5.42 × 10−10Organismal SystemsEndocrine system
ko04974Protein digestion and absorption3.60 × 10−8Organismal SystemsDigestive system
ko04972Pancreatic secretion3.60 × 10−8Organismal SystemsDigestive system
ko05146Amoebiasis3.60 × 10−8Human DiseasesInfectious diseases: Parasitic
ko00860Porphyrin and chlorophyll metabolism9.70 × 10−8MetabolismMetabolism of cofactors and vitamins
ko00830Retinol metabolism1.60 × 10−7MetabolismMetabolism of cofactors and vitamins
ko05110Vibrio cholerae infection3.06 × 10−6Human DiseasesInfectious diseases: Bacterial
ko00140Steroid hormone biosynthesis5.01 × 10−6MetabolismLipid metabolism
ko00514Other types of O-glycan biosynthesis5.78 × 10−5MetabolismGlycan biosynthesis and metabolism
ko00980Metabolism of xenobiotics by cytochrome P4505.78 × 10−5MetabolismXenobiotics biodegradation and metabolism
ko00982Drug metabolism—cytochrome P4506.99 × 10−5MetabolismXenobiotics biodegradation and metabolism
ko00740Riboflavin metabolism1.09 × 10−4MetabolismMetabolism of cofactors and vitamins
ko04976Bile secretion1.83 × 10−3Organismal SystemsDigestive system
ko00480Glutathione metabolism1.83 × 10−3MetabolismMetabolism of other amino acids
ko00053Ascorbate and aldarate metabolism2.27 × 10−3MetabolismCarbohydrate metabolism
ko05164Influenza A2.36 × 10−3Human DiseasesInfectious diseases: Viral
ko05221Acute myeloid leukemia2.83 × 10−3Human DiseasesCancers: Specific types
ko04640Hematopoietic cell lineage4.31 × 10−3Organismal SystemsImmune system
ko00130Ubiquinone and other terpenoid-quinone biosynthesis8.05 × 10−3MetabolismMetabolism of cofactors and vitamins
ko00592alpha-Linolenic acid metabolism8.05 × 10−3MetabolismLipid metabolism
ko00983Drug metabolism—other enzymes8.05 × 10−3MetabolismXenobiotics biodegradation and metabolism
ko04512ECM—receptor interaction1.21 × 10−2Environmental Information ProcessingSignaling molecules and interaction
ko00040Pentose and glucuronate interconversions1.93 × 10−2MetabolismCarbohydrate metabolism
ko01100Metabolic pathways2.07 × 10−2MetabolismGlobal map
ko00350Tyrosine metabolism2.31 × 10−2MetabolismAmino acid metabolism
ko00030Pentose phosphate pathway2.48 × 10−2MetabolismCarbohydrate metabolism
ko04142Lysosome2.57 × 10−2Cellular ProcessesTransport and catabolism
ko05323Rheumatoid arthritis2.57 × 10−2Human DiseasesImmune diseases
ko00981Insect hormone biosynthesis2.61 × 10−2MetabolismMetabolism of terpenoids and polyketides
ko05130Pathogenic Escherichia coli infection2.61 × 10−2Human DiseasesInfectious diseases: Bacterial
ko04977Vitamin digestion and absorption2.70 × 10−2Organismal SystemsDigestive system
ko00360Phenylalanine metabolism4.50 × 10−2MetabolismAmino acid metabolism
Genes with irradiated-specific SNPs or control-specific SNPs were also considered candidate genes that responded to irradiation. As shown in Table S6, a total of 130 SNPs in 125 unigenes were identified.
Table 4. Validation results of the selected differentially expressed genes.
Table 4. Validation results of the selected differentially expressed genes.
Gene NameUnigeneGene LengthFold Change by Sequencingp ValueFold by qPCR
ENPEPCL1507.contig1_All31324.02 4.90 × 10−31.51
ENPEPunigene5812_All8468.93 1.29 × 10−41.95
JheCL6471.contig1_All171323.44 3.68 × 10−22.76
mthCL6547.contig2_All10446.75 4.54 × 10−410.36
NotchCL4940.contig1_All99775.78 2.27 × 10−42.73
sqdCL816.contig2_All687−4.31 3.14 × 10−3−2.05
Ugtunigene16123_All19114.31 1.07 × 10−31.38
Ugtunigene6574_All10077.14 2.77 × 10−41.55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, C.; Shabbir, S.; Wang, W.; Gao, Y.; Ge, C.; Lin, Q. Transcriptome Analysis Revealed Genes Related to γ-Irradiation Induced Emergence Failure in Third-Instar Larvae of Bactrocera dorsalis. Insects 2022, 13, 1017. https://doi.org/10.3390/insects13111017

AMA Style

Sun C, Shabbir S, Wang W, Gao Y, Ge C, Lin Q. Transcriptome Analysis Revealed Genes Related to γ-Irradiation Induced Emergence Failure in Third-Instar Larvae of Bactrocera dorsalis. Insects. 2022; 13(11):1017. https://doi.org/10.3390/insects13111017

Chicago/Turabian Style

Sun, Chao, Samina Shabbir, Wenxiang Wang, Yan Gao, Cuicui Ge, and Qingsheng Lin. 2022. "Transcriptome Analysis Revealed Genes Related to γ-Irradiation Induced Emergence Failure in Third-Instar Larvae of Bactrocera dorsalis" Insects 13, no. 11: 1017. https://doi.org/10.3390/insects13111017

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