De novo characterization of the Dialeurodes citri transcriptome: mining genes involved in stress resistance and simple sequence repeats (SSRs) discovery

The metabolism and detoxification effects of lead exposure on Aleurolyphus ovatus (Acari: Acaridae) via transcriptome analysis

Aleurolyphus ovatus Troupeau is one of the most predominant species of the Acaridae family worldwide. Recent reports have demonstrated that the accumulation of lead in stored grains and dietary items exceeds the required standards. However, the molecular mechanism of heavy metal stress on mites has not been reported. To understand the mechanism underlying the heavy metal response of A. ovatus, comparative transcriptome analysis was performed in this study using an Illumina high throughput mRNA sequencing (RNA-seq) platform. A. ovatus was fed on artificial diets containing two different concentrations of lead, namely, a low concentration of 12.5 mg/kg (LAO) and a high concentration of 100 mg/kg (HAO), while the mites in the control (NAO) group were not exposed to lead. A total of 44,362 unigenes, with an average length of 1547 bp, were identified. Of these, 996 unigenes were successfully annotated in seven functional databases. The number of differentially expressed genes (DEGs) in A. ovatus under different lead concentrations was compared. In NAO versus LAO group, including 310 up-regulated and 1580 down-regulated DEGs. In NAO versus HAO group, including 3928 up-regulated and 1761 down-regulated DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that detoxification enzyme genes were significantly expressed in pathways, such as cytochrome P450 foreign body metabolism, glutathione metabolism and drug metabolism-cytochrome pathway. The results of gene annotation and quantitative real-time PCR showed that high concentration of lead significantly stimulated the expression of metabolic detoxification enzyme genes such as glutathione S transferase (GST) and superoxide dismutase (SOD), while low concentration inhibited their expression. This study will provide a basis for the molecular mechanism of A. ovatus in response to heavy metal lead stimulation in stored grain.

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

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

RNA-seq Analysis Reveals Mitochondrial and Cuticular Protein Genes Are Associated with Phosphine Resistance in the Rusty Grain Beetle (Coleoptera:Laemophloeidae)

The rusty grain beetle, Cryptolestes ferrugineus (Stephens), is a serious pest of stored grain, which has developed high levels of resistance to phosphine. In this study, five geographically distant populations of C. ferrugineus had been collected in China, specifically in granaries where phosphine fumigant is used for pest control, and they showed a high resistance ratio up to 1,907 (LC50 = 21.0 mg/liter). Then, a reference transcriptome was constructed to use as a basis for investigating the molecular mechanisms of phosphine resistance in this species, which consisted of 47,006 unigenes with a mean length of 1,090. Subsequently, the RNA-Seq analysis of individuals from the most susceptible and resistant populations led to the identification of 54 genes that are differentially expressed. GO and KEGG analysis demonstrated that genes associated with mitochondrial and respiration functions were significantly enriched. Also, the 'structural constituent of cuticle' term was annotated in the GO enrichment analysis and further qRT-PCR confirmed that the expression levels of nine cuticular protein genes were significantly increased in the resistant population. In conclusion, we present here a transcriptome-wide overview of gene expression changes between resistant and susceptible populations of C. ferrugineus, and this in turn documents that mitochondria and cuticular protein genes may play together a crucial role in phosphine resistance. Further gene function analysis should enable the provision of advice to expedite resistance management decisions.

A De Novo Transcriptomics Approach Reveals Genes Involved in Thrips Tabaci Resistance to Spinosad

The onion thrip, Thrips tabaci (Thysanoptera: Thripidae) is a major polyphagous pest that attacks a wide range of economically important crops, especially Allium species. The thrip's damage can result in yield loss of up to 60% in onions (Allium cepa). In the past few decades, thrip resistance to insecticides with various modes of actions have been documented. These include resistance to spinosad, a major active compound used against thrips, which was reported from Israel. Little is known about the molecular mechanisms underlying spinosad resistance in T. tabaci. We attempted to characterize the mechanisms involved in resistance to spinosad using quantitative transcriptomics. Susceptible (LC50 = 0.6 ppm) and resistant (LC50 = 23,258 ppm) thrip populations were collected from Israel. An additional resistant population (LC50 = 117 ppm) was selected in the laboratory from the susceptible population. De novo transcriptome analysis on the resistant and susceptible population was conducted to identify differently expressed genes (DGEs) that might be involved in the resistance against spinosad. In this analysis, 25,552 unigenes were sequenced, assembled, and functionally annotated, and more than 1500 DGEs were identified. The expression levels of candidate genes, which included cytochrome P450 and vittelogenin, were validated using quantitative RT-PCR. The cytochrome P450 expression gradually increased with the increase of the resistance. Higher expression levels of vitellogenin in the resistant populations were correlated with higher fecundity, suggesting a positive effect of the resistance on resistant populations. This research provides a novel genetic resource for onion thrips and a comprehensive molecular examination of resistant populations to spinosad. Those resources are important for future studies concerning thrips and resistance in insect pests regarding agriculture.

Detoxification enzymes associated with butene-fipronil resistance in Epacromius coerulipes

BACKGROUND

Epacromius coerulipes is a widely distributed locust pest species. Chemical control is the main method used to kill locusts; however, this can result in the selection of locusts with resistance to chemical pesticides. Therefore, the study of resistance is of great significance for the sustainable management of locusts.

RESULTS

In this study, to investigate the relationship between detoxification enzymes and butene-fipronil resistance in E. coerulipes, resistant strains of the locust were compared with sensitive strains. The synergism of synergistic agents was significantly enhanced, and the activities of multifunctional oxidase, carboxylesterase, and glutathione sulfur transferase were significantly increased. Transcriptome sequencing revealed 226 detoxification enzyme genes and 23 upregulated genes. Neighbor-joining was used to construct a phylogenetic tree of related gene families, which included 59 P450 genes, 52 carboxylesterases (CarE) genes, and 25 glutathione S-transferase (GST) genes. Reverse transcription polymerase chain reaction (RT-PCR) analysis results of overexpressed genes in the resistant population combined with a phylogenetic tree showed that four P450 genes belonged to the CYP6, CYP4, CYP18 and CYP302 families, two CarE genes belonged to Clade A families, and one GST gene belonged to the Sigma family. These family members were annotated as detoxification enzyme genes of metabolic insecticide in the transcriptome databases.

CONCLUSIONS

This study showed that P450, CarE and GST together resulted in moderate resistance to butene-fipronil in locusts. The analysis revealed several overexpressed detoxification enzyme genes that will be the focus of future studies on the mechanism of resistance to butene-fipronil. © 2019 Society of Chemical Industry.

Single-Molecule Long-Read Sequencing of Zanthoxylum bungeanum Maxim. Transcriptome: Identification of Aroma-Related Genes

Zanthoxylum bungeanum Maxim. is an economically important tree species that is resistant to drought and infertility, and has potential medicinal and edible value. However, comprehensive genomic data are not yet available for this species, limiting its potential utility for medicinal use, breeding programs, and cultivation. Transcriptome sequencing provides an effective approach to remedying this shortcoming. Herein, single-molecule long-read sequencing and next-generation sequencing approaches were used in parallel to obtain transcript isoform structure and gene functional information in Z. bungeanum. In total, 282,101 reads of inserts (ROIs) were identified, including 134,074 full-length non-chimeric reads, among which 65,711 open reading frames (ORFs), 50,135 simple sequence repeats (SSRs), and 1492 long non-coding RNAs (lncRNAs) were detected. Functional annotation revealed metabolic pathways related to aroma components and color characteristics in Z. bungeanum. Unexpectedly, 30 transcripts were annotated as genes involved in regulating the pathogenesis of breast and colorectal cancers. This work provides a comprehensive transcriptome resource for Z. bungeanum, and lays a foundation for the further investigation and utilization of Zanthoxylum resources.

RNA interference-mediated silencing of genes involved in the immune responses of the soybean pod borer Leguminivora glycinivorella (Lepidoptera: Olethreutidae)

RNA interference (RNAi) technology may be useful for developing new crop protection strategies against the soybean pod borer (SPB; Leguminivora glycinivorella), which is a critical soybean pest in northeastern Asia. Immune-related genes have been recently identified as potential RNAi targets for controlling insects. However, little is known about these genes or mechanisms underlying their expression in the SPB. In this study, we completed a transcriptome-wide analysis of SPB immune-related genes. We identified 41 genes associated with SPB microbial recognition proteins, immune-related effectors or signalling molecules in immune response pathways (e.g., Toll and immune deficiency pathways). Eleven of these genes were selected for a double-stranded RNA artificial feeding assay. The down-regulated expression levels of LgToll-5-1a and LgPGRP-LB2a resulted in relatively high larval mortality rates and abnormal development. Our data represent a comprehensive genetic resource for immune-related SPB genes, and may contribute to the elucidation of the mechanism regulating innate immunity in Lepidoptera species. Furthermore, two immune-related SPB genes were identified as potential RNAi targets, which may be used in the development of RNAi-mediated SPB control methods.

Functional Analysis of RNA Interference-Related Soybean Pod Borer (Lepidoptera) Genes Based on Transcriptome Sequences

RNA interference (RNAi) is useful for controlling pests of agriculturally important crops. The soybean pod borer (SPB) is the most important soybean pest in Northeastern Asia. In an earlier study, we confirmed that the SPB could be controlled via transgenic plant-mediated RNAi. Here, the SPB transcriptome was sequenced to identify RNAi-related genes, and also to establish an RNAi-of-RNAi assay system for evaluating genes involved in the SPB systemic RNAi response. The core RNAi genes, as well as genes potentially involved in double-stranded RNA (dsRNA) uptake were identified based on SPB transcriptome sequences. A phylogenetic analysis and the characterization of these core components as well as dsRNA uptake related genes revealed that they contain conserved domains essential for the RNAi pathway. The results of the RNAi-of-RNAi assay involving Laccase 2 (a critical cuticle pigmentation gene) as a marker showed that genes encoding the sid-like (Sil1), scavenger receptor class C (Src), and scavenger receptor class B (Srb3 and Srb4) proteins of the endocytic pathway were required for SPB cellular uptake of dsRNA. The SPB response was inferred to contain three functional small RNA pathways (i.e., miRNA, siRNA, and piRNA pathways). Additionally, the SPB systemic RNA response may rely on systemic RNA interference deficient transmembrane channel-mediated and receptor-mediated endocytic pathways. The results presented herein may be useful for developing RNAi-mediated methods to control SPB infestations in soybean.

Molecular cloning, mRNA expression and alternative splicing of a ryanodine receptor gene from the citrus whitefly, Dialeurodes citri (Ashmead)

Insect ryanodine receptors are the main targets of diamide insecticides that have highly selective insecticidal activity but are less toxic to mammals. Therefore, these insecticides are ideal for pest control. Ryanodine receptors (RyRs) play a critical role in Ca²⁺ signaling in muscle and non-muscle cells. In this study, we cloned the complete cDNA (DcRyR) of the RyR from the citrus whitefly, Dialeurodes citri, a serious pest of citrus orchards in China. The open reading frame of RyR is 15,378bp long and encodes a protein with 5126 amino acids with a computed molecular weight of 579.523kDa. DcRyR shows a high amino acid sequence identity to RyRs from other insects (76%-95%) and low identity to those from nematodes and mammals (44%-52%). DcRyR shares many features of insect and vertebrate RyRs, including a MIR domain, two RIH domains, three SPRY domains, four copies of RyR repeat domain, RIH-associated domain at the N-terminus, two consensus calcium-binding EF-hands and six transmembrane domains at the C-terminus. The expression of DcRyR mRNA was the highest in the nymphs and lowest in eggs; DcRyR mRNA was 1.85-fold higher in the nymphs than in the eggs. Among the tissues, DcRyR mRNA expression was 4.18- and 4.02-fold higher in the adult head and thorax than in the abdomen. DcRyR had three alternative splice sites and the splice variants showed body part-specific expression and were developmentally regulated. These results may help investigate target-based resistance to diamide insecticides in D. citri.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Transcriptome Profiling Analysis of Wolf Spider Pardosa pseudoannulata (Araneae: Lycosidae) after Cadmium Exposure

Pardosa pseudoannulata is one of the most common wandering spiders in agricultural fields and a potentially good bioindicator for heavy metal contamination. However, little is known about the mechanisms by which spiders respond to heavy metals at the molecular level. In the present study, high-throughput transcriptome sequencing was employed to characterize the de novo transcriptome of the spiders and to identify differentially expressed genes (DEGs) after cadmium exposure. We obtained 60,489 assembled unigenes, 18,773 of which were annotated in the public databases. A total of 2939 and 2491 DEGs were detected between the libraries of two Cd-treated groups and the control. Functional enrichment analysis revealed that metabolism processes and digestive system function were predominately enriched in response to Cd stress. At the cellular and molecular levels, significantly enriched pathways in lysosomes and phagosomes as well as replication, recombination and repair demonstrated that oxidative damage resulted from Cd exposure. Based on the selected DEGs, certain critical genes involved in defence and detoxification were analysed. These results may elucidate the molecular mechanisms underlying spiders' responses to heavy metal stress.

Identification of Immunity-Related Genes in Dialeurodes citri against Entomopathogenic Fungus Lecanicillium attenuatum by RNA-Seq Analysis

Dialeurodes citri is a major pest in citrus producing areas, and large-scale outbreaks have occurred increasingly often in recent years. Lecanicillium attenuatum is an important entomopathogenic fungus that can parasitize and kill D. citri. We separated the fungus from corpses of D. citri larvae. However, the sound immune defense system of pests makes infection by an entomopathogenic fungus difficult. Here we used RNA sequencing technology (RNA-Seq) to build a transcriptome database for D. citri and performed digital gene expression profiling to screen genes that act in the immune defense of D. citri larvae infected with a pathogenic fungus. De novo assembly generated 84,733 unigenes with mean length of 772 nt. All unigenes were searched against GO, Nr, Swiss-Prot, COG, and KEGG databases and a total of 28,190 (33.3%) unigenes were annotated. We identified 129 immunity-related unigenes in transcriptome database that were related to pattern recognition receptors, information transduction factors and response factors. From the digital gene expression profile, we identified 441 unigenes that were differentially expressed in D. citri infected with L. attenuatum. Through calculated Log₂Ratio values, we identified genes for which fold changes in expression were obvious, including cuticle protein, vitellogenin, cathepsin, prophenoloxidase, clip-domain serine protease, lysozyme, and others. Subsequent quantitative real-time polymerase chain reaction analysis verified the results. The identified genes may serve as target genes for microbial control of D. citri.

Antennal transcriptome and differential expression of olfactory genes in the yellow peach moth, Conogethes punctiferalis (Lepidoptera: Crambidae)

The yellow peach moth (YPM), Conogethes punctiferalis (Guenée), is a multivoltine insect pest of crops and fruits. Antennal-expressed receptors are important for insects to detect olfactory cues for host finding, mate attraction and oviposition site selection. However, few olfactory related genes were reported in YPM until now. In the present study, we sequenced and characterized the antennal transcriptomes of male and female YPM. In total, 15 putative odorant binding proteins (OBPs), 46 putative odorant receptors (ORs) and 7 putative ionotropic receptors (IRs) were annotated and identified as olfactory-related genes of C. punctiferalis. Further analysis of RT-qPCR revealed that all these olfactory genes are primarily or uniquely expressed in male and female antennae. Among which, 3 OBPs (OBP4, OBP8 and PBP2) and 4 ORs (OR22, OR26, OR44 and OR46) were specially expressed in male antennae, whereas 4 ORs (OR5, OR16, OR25 and OR42) were primarily expressed in female antennae. The predicted protein sequences were compared with homologs in other lepidopteran species and model insects, which showed high sequence homologies between C. punctiferalis and O. furnacalis. Our work allows for further functional studies of pheromone and general odorant detection genes, which might be meaningful targets for pest management.

Transcriptome and Difference Analysis of Fenpropathrin Resistant Predatory Mite, Neoseiulus barkeri (Hughes)

Several fenpropathrin-resistant predatory mites have been reported. However, the molecular mechanism of the resistance remains unknown. In the present study, the Neoseiulus barkeri (N. barkeri) transcriptome was generated using the Illumina sequencing platform, 34,211 unigenes were obtained, and 15,987 were manually annotated. After manual annotation, attentions were attracted to resistance-related genes, such as voltage-gated sodium channel (VGSC), cytochrome P450s (P450s), and glutathione S-transferases (GSTs). A polymorphism analysis detected two point mutations (E1233G and S1282G) in the linker region between VGSC domain II and III. In addition, 43 putative P450 genes and 10 putative GST genes were identified from the transcriptome. Among them, two P450 genes, NbCYP4EV2 and NbCYP4EZ1, and four GST genes, NbGSTd01, NbGSTd02, NbGSTd03 and NbGSTm03, were remarkably overexpressed 3.64-46.69-fold in the fenpropathrin resistant strain compared to that in the susceptible strain. These results suggest that fenpropathrin resistance in N. barkeri is a complex biological process involving many genetic changes and provide new insight into the N. barkeri resistance mechanism.

Differential Gene Expression Analysis of the Epacromius coerulipes (Orthoptera: Acrididae) Transcriptome

Epacromius coerulipes (Ivanov) is one of the most widely distributed locusts. To date, the main methods to kill locusts still rely on chemical controls, which can result in the selection of locusts with resistance to chemical pesticides. Butene-fipronil is a new pesticide that was discovered by the structural modification of fipronil. This pesticide has been used to control various agricultural pests and has become an important pesticide product to control pests that exhibit resistance to other pesticides, including locusts. To extend its useful half-life, studies of the initiation and progression of resistance to this pesticide are needed. Herein, two E. coerulipes strains, a pesticide-sensitive (PS) and a pesticide-resistant (PR) strain, were chosen to undergo de novo assembly by paired-end transcriptome Illumina sequencing. Overall, 63,033 unigenes were detected; the average gene length was 772 bp and the N50 was 1,589 bp. Among these unigenes, ∼ 25,132 (39.87% of the total) could be identified as known proteins in bioinformatic databases from national centers. A comparison of the PR and PS strains revealed that 2,568 genes were differentially expressed, including 1,646 and 922 genes that were up- and down-regulated, respectively. According to the Gene Ontology (GO) database, among biological processes the metabolic process group was the largest group (6,900 genes, 22.47%) and contained a high frequency of differentially expressed genes (544 genes, 27.54%). According to the Clusters of Orthologous Groups (COG) categories, 28 genes, representing 2.98% of all genes, belonged to the group of genes involved in the biosynthesis, transportation, and catabolism of secondary metabolites. The differentially expressed genes that we identified are involved in 50 metabolic pathways. Among these pathways, the metabolism pathway was the most represented. After enrichment analysis of differential gene expression pathways, six pathways--ribosome; starch, and sucrose metabolism; ascorbate and aldarate metabolism; drug metabolism-cytochrome P450; metabolism of xenobiotics by cytochrome P450; and glutathione metabolism--showed a high degree of enrichment. Among these pathways, drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, and glutathione metabolism have been associated with pesticide metabolism. Furthermore, 316 unigenes in the E. coerulipes transcriptome encode detoxifying enzymes and 76 unigenes encode target proteins of pesticides. Among these genes, 23 genes that encode detoxifying enzymes in the resistance group were found to be up-regulated. The transcriptome sequencing results of E. coerulipes established a genomics database of E. coerulipes for the first time. This study also establishes a molecular basis for gene function analysis of E. coerulipes Moreover, it provides a theoretical resource for mechanistic studies on pesticide resistance through the screening and investigation of resistance genes.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.

Transcriptome Analysis of the Carmine Spider Mite, Tetranychus cinnabarinus (Boisduval, 1867) (Acari: Tetranychidae), and Its Response to β-Sitosterol

Tetranychus cinnabarinus (Acari: Tetranychidae) is a worldwide polyphagous agricultural pest that has the title of resistance champion among arthropods. We reported previously the identification of the acaricidal compound β-sitosterol from Mentha piperita and Inula japonica. However, the acaricidal mechanism of β-sitosterol is unclear. Due to the limited genetic research carried out, we de novo assembled the transcriptome of T. cinnabarinus using Illumina sequencing and conducted a differential expression analysis of control and β-sitosterol-treated mites. In total, we obtained >5.4 G high-quality bases for each sample with unprecedented sequencing depth and assembled them into 22,941 unigenes. We identified 617 xenobiotic metabolism-related genes involved in detoxification, binding, and transporting of xenobiotics. A highly expanded xenobiotic metabolic system was found in mites. T. cinnabarinus detoxification genes-including carboxyl/cholinesterase and ABC transporter class C-were upregulated after β-sitosterol treatment. Defense-related proteins, such as Toll-like receptor, legumain, and serine proteases, were also activated. Furthermore, other important genes-such as the chloride channel protein, cytochrome b, carboxypeptidase, peritrophic membrane chitin binding protein, and calphostin-may also play important roles in mites' response to β-sitosterol. Our results demonstrate that high-throughput-omics tool facilitates identification of xenobiotic metabolism-related genes and illustration of the acaricidal mechanisms of β-sitosterol.

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

The testis is a highly specialized tissue that plays a vital role in ensuring fertility by producing spermatozoa, which are transferred to the female during mating. Spermatogenesis is a complex process, resulting in the production of mature sperm, and involves significant structural and biochemical changes in the seminiferous epithelium of the adult testis. The identification of genes involved in spermatogenesis of Bactrocera dorsalis (Hendel) is critical for a better understanding of its reproductive development. In this study, we constructed a cDNA library of testes from male B. dorsalis adults at different ages, and performed de novo transcriptome sequencing to produce a comprehensive transcript data set, using Illumina sequencing technology. The analysis yielded 52 016 732 clean reads, including a total of 4.65 Gb of nucleotides. These reads were assembled into 47 677 contigs (average 443 bp) and then clustered into 30 516 unigenes (average 756 bp). Based on BLAST hits with known proteins in different databases, 20 921 unigenes were annotated with a cut-off E-value of 10(-5). The transcriptome sequences were further annotated using the Clusters of Orthologous Groups, Gene Orthology and the Kyoto Encyclopedia of Genes and Genomes databases. Functional genes involved in spermatogenesis were analysed, including cell cycle proteins, metalloproteins, actin, and ubiquitin and antihyperthermia proteins. Several testis-specific genes were also identified. The transcripts database will help us to understand the molecular mechanisms underlying spermatogenesis in B. dorsalis. Furthermore, 2913 simple sequence repeats and 151 431 single nucleotide polymorphisms were identified, which will be useful for investigating the genetic diversity of B. dorsalis in the future.