Search

The damage caused by Bradysia odoriphaga is the main factor threatening the production of vegetables in the Liliaceae family. However, few genetic studies of B. odoriphaga have been conducted because of a lack of genomic resources. Many long-read sequencing technologies have been developed in the last decade; therefore, in this study, the transcriptome including all development stages of B. odoriphaga was sequenced for the first time by Pacific single-molecule long-read sequencing. Here, 39,129 isoforms were generated, and 35,645 were found to have annotation results when checked against sequences available in different databases. Overall, 18,473 isoforms were distributed in 25 various Clusters of Orthologous Groups, and 11,880 isoforms were categorized into 60 functional groups that belonged to the three main Gene Ontology classifications. Moreover, 30,610 isoforms were assigned into 44 functional categories belonging to six main Kyoto Encyclopedia of Genes and Genomes functional categories. Coding DNA sequence (CDS) prediction showed that 36,419 out of 39,129 isoforms were predicted to have CDS, and 4319 simple sequence repeats were detected in total. Finally, 266 insecticide resistance and metabolism-related isoforms were identified as candidate genes for further investigation of insecticide resistance and metabolism in B. odoriphaga.

Aphids cause serious losses to the production of wheat. The grain aphid, Sitobion avenae, which is the dominant species of aphid in all wheat regions of China, is resistant to a variety of insecticides, including imidacloprid and chlorpyrifos. However, the resistance and mechanism of insecticide tolerance of S. avenae are still unclear. Therefore, this study employed transcriptome analysis to compare the expression patterns of stress response genes under imidacloprid and chlorpyrifos treatment for 15 min, 3 h, and 36 h of exposure. S. avenae adult transcriptome was assembled and characterized first, after which samples treated with insecticides for different lengths of time were compared with control samples, which revealed 602267 differentially expressed unigenes (DEUs). Among these DEUs, 31&ndash, 790 unigenes were classified into 66&ndash, 786 categories of gene ontology (GO) functional groups, and 24&ndash, 760 DEUs could be mapped into 54&ndash, 268 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Finally, 11 insecticide-tolerance-related unigenes were chosen to confirm the relative expression by quantitative real-time polymerase chain reaction (qRT-PCR) in each treatment. Most of the results between qRT-PCR and RNA sequencing (RNA-Seq) are well-established. The results presented herein will facilitate molecular research investigating insecticide resistance in S. avenae, as well as in other wheat aphids.