De novo sequencing-based transcriptome and digital gene expression analysis reveals insecticide resistance-relevant genes in Propylaea japonica (Thunberg) (Coleoptea: Coccinellidae)

Omics approaches to unravel insecticide resistance mechanism in Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae)

Bemisia tabaci (Gennadius) whitefly (BtWf) is an invasive pest that has already spread worldwide and caused major crop losses. Numerous strategies have been implemented to control their infestation, including the use of insecticides. However, prolonged insecticide exposures have evolved BtWf to resist these chemicals. Such resistance mechanism is known to be regulated at the molecular level and systems biology omics approaches could shed some light on understanding this regulation wholistically. In this review, we discuss the use of various omics techniques (genomics, transcriptomics, proteomics, and metabolomics) to unravel the mechanism of insecticide resistance in BtWf. We summarize key genes, enzymes, and metabolic regulation that are associated with the resistance mechanism and review their impact on BtWf resistance. Evidently, key enzymes involved in the detoxification system such as cytochrome P450 (CYP), glutathione S-transferases (GST), carboxylesterases (COE), UDP-glucuronosyltransferases (UGT), and ATP binding cassette transporters (ABC) family played key roles in the resistance. These genes/proteins can then serve as the foundation for other targeted techniques, such as gene silencing techniques using RNA interference and CRISPR. In the future, such techniques will be useful to knock down detoxifying genes and crucial neutralizing enzymes involved in the resistance mechanism, which could lead to solutions for coping against BtWf infestation.

Behavioral and Transcriptomic Analyses in the Indoxacarb Response of a Non-Target Damselfly Species

Ischnura senegalensis, which widely spreads in paddy fields, has the potential to be used as a natural predator of insect pests. However, the application of insecticides in the field could pose a threat to the survival of I. senegalensis. Among these pesticides, indoxacarb, an oxadiazine insecticide, is renowned for its broad-spectrum efficacy against numerous insect pests. In this study, we examined the toxicity of indoxacarb towards the larvae of I. senegalensis. Behavioral experiments and transcriptome analyses were conducted under indoxacarb treatments. Results revealed that indoxacarb induced abnormal body gestures and significant locomotory impairments, which could ultimately reduce the survival rate of the larvae in their natural habitat. Moreover, transcriptome analyses indicated that genes related to muscle function were significantly affected. Interestingly, at lower concentrations of indoxacarb (0.004 mg/L), the larvae seem to detoxify the indoxacarb with the aid of the cytochrome P450 gene. However, under higher concentrations (0.4 mg/L), the sensory abilities of the larvae were significantly diminished, and they were unable to degrade the toxicity of indoxacarb. Our study underscores the importance of carefully evaluating the impact of insecticides on non-target predatory insects before their widespread application.

Sublethal toxicity, transgenerational effects, and transcriptome expression of the neonicotinoid pesticide cycloxaprid on demographic fitness of Coccinella septempunctata

Evaluating side effects of new neonicotinoids in terms of sublethal doses and transcriptome expression is a crucial but challenging part of integrated pest management (IPM) approaches. To this end, a study of lethal and sublethal effects on Coccinella septempunctata larvae was conducted, and an age-stage, two-sex life table procedure was performed to investigate life-table parameters. Cycloxaprid (CYC) was shown to have adverse effects on survival, development, total longevity, reproductive capacity, and predation ability in C. septempunctata. In addition, demographic growth parameters of the F1 generation such as net reproductive rate, and the intrinsic and finite rates of increase were significantly decreased under sublethal dosage LR₃₀ (1.91 g ai/hm²). These results demonstrated that the population growth of C. septempunctata was impacted by a sublethal dosage of CYC. For transcriptome expression, 544 up- and 338 down-regulated significantly differentially expressed genes (DEGs), were observed between LR₃₀ treatment and control groups. Moreover, pathways related to metabolism of retinol, carcinogenesis, biosynthesis of steroid hormone, P450 metabolism, and metabolism of xenobiotics were identified in KEGG pathway analysis. Ten DEGs were chosen and confirmed with quantitative real-time PCR analysis. Based on these findings, CYC should be considered as a component of IPM strategies in the field.

Transcriptome analysis reveals differential effects of beta-cypermethrin and fipronil insecticides on detoxification mechanisms in Solenopsis invicta

Insecticide resistance poses many challenges in insect pest control, particularly in the control of destructive pests such as red imported fire ants (Solenopsis invicta). In recent years, beta-cypermethrin and fipronil have been extensively used to manage invasive ants, but their effects on resistance development in S. invicta are still unknown. To investigate resistance development, S. invicta was collected from populations in five different cities in Guangdong, China. The results showed 105.71- and 2.98-fold higher resistance against fipronil and beta-cypermethrin, respectively, in the Guangzhou population. The enzymatic activities of acetylcholinesterase, carboxylases, and glutathione S-transferases significantly increased with increasing beta-cypermethrin and fipronil concentrations. Transcriptomic analysis revealed 117 differentially expressed genes (DEGs) in the BC-ck vs. BC-30 treatments (39 upregulated and 78 downregulated), 109 DEGs in F-ck vs. F-30 (33 upregulated and 76 downregulated), and 499 DEGs in BC-30 vs. F-30 (312 upregulated and 187 downregulated). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that DEGs associated with insecticide resistance were significantly enriched in metabolic pathways, the AMPK signaling pathway, the insulin signaling pathway, carbon metabolism, peroxisomes, fatty acid metabolism, drug metabolism enzymes and the metabolism of xenobiotics by cytochrome P450. Furthermore, we found that DEGs important for insecticide detoxification pathways were differentially regulated under both insecticide treatments in S. invicta. Comprehensive transcriptomic data confirmed that detoxification enzymes play a significant role in insecticide detoxification and resistance development in S. invicta in Guangdong Province. Numerous identified insecticide-related genes, GO terms, and KEGG pathways indicated the resistance of S. invicta workers to both insecticides. Importantly, this transcriptome profile variability serves as a starting point for future research on insecticide risk evaluation and the molecular mechanism of insecticide detoxification in invasive red imported fire ants.

Emamectin benzoate induced enzymatic and transcriptional alternation in detoxification mechanism of predatory beetle Paederus fuscipes (Coleoptera: Staphylinidae) at the sublethal concentration

In this study, the detoxification enzyme activity and the transcriptional profile changes in the second instar through RNA-sequencing technology due to emamectin benzoate (EB) were assessed. The cytochrome P450 monooxygenases (P450) enzyme activity was not altered by EB due to the change in concentration and exposure time in all treatments. The glutathione S-transferase (GST) enzyme was not considerably varying in all treatments, while exposure time significantly changed the enzyme activity. Results showed that the esterase (Ests) activity was elevated with the increasing concentrations and exposure time. Two libraries were generated, containing 107,767,542 and 108,142,289 clean reads for the samples treated with LC₃₀ of EB and control. These reads were grouped into 218,070 transcripts and 38,097 unigenes. A total of 2257 differentially expressed genes (DEGs) were identified from these unigenes, of which 599 up-regulated and 1658 were down-regulated. The majority of these DEGs related to pesticides resistance were identified in numerous Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, e.g., steroid hormone biosynthesis, glutathione metabolism, drug metabolism-other enzymes, chemical carcinogenesis, pathways of cancer, metabolism of xenobiotics by cytochrome P450, drug metabolism of cytochrome P450, linoleic acid metabolism, retinol metabolism, and insect hormone biosynthesis. These pathways also shared the same genes as cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), Esterase (Ests), UDP-glucosyltransferases (UGTs), and ATP-binding cassettes (ABCs). A heatmap analysis also showed that regulation of genes in a pathway causes a series of gene expression regulation in subsequent pathways. Our quantitative reverse transcription-PCR (qRT-PCR) results were consistent with the DEG's data of transcriptome analysis. The comprehensive transcriptome sequence resource attained through this study evidence that the EB induces significant modification in enzyme activity and transcriptome profile of Paederus fuscipes, which may enable more significant molecular underpinnings behind the insecticide-resistance mechanism for further investigations.

Next-generation biological control: the need for integrating genetics and genomics

Biological control is widely successful at controlling pests, but effective biocontrol agents are now more difficult to import from countries of origin due to more restrictive international trade laws (the Nagoya Protocol). Coupled with increasing demand, the efficacy of existing and new biocontrol agents needs to be improved with genetic and genomic approaches. Although they have been underutilised in the past, application of genetic and genomic techniques is becoming more feasible from both technological and economic perspectives. We review current methods and provide a framework for using them. First, it is necessary to identify which biocontrol trait to select and in what direction. Next, the genes or markers linked to these traits need be determined, including how to implement this information into a selective breeding program. Choosing a trait can be assisted by modelling to account for the proper agro-ecological context, and by knowing which traits have sufficiently high heritability values. We provide guidelines for designing genomic strategies in biocontrol programs, which depend on the organism, budget, and desired objective. Genomic approaches start with genome sequencing and assembly. We provide a guide for deciding the most successful sequencing strategy for biocontrol agents. Gene discovery involves quantitative trait loci analyses, transcriptomic and proteomic studies, and gene editing. Improving biocontrol practices includes marker-assisted selection, genomic selection and microbiome manipulation of biocontrol agents, and monitoring for genetic variation during rearing and post-release. We conclude by identifying the most promising applications of genetic and genomic methods to improve biological control efficacy.

Bt, Not a Threat to Propylea japonica

Given the ever-increasing commercial planting of transgenic plants across the world, an evaluation of their impacts on non-target organisms is as an important part of the risk assessment process. Propylea japonica is a dominant non-target predator and pollen feeder insect that is prevalent in Bt cotton fields, and it is thus in direct contact with Bt proteins. However, the effect of Bt proteins on P. japonica has not received much attention. In this study, the effects of Cry1Ac and/or Cry2Ab proteins on P. japonica were investigated from three aspects. First, no significant differences in the diversity of the microbiota nor change in species composition and community structure were observed among Cry protein treatments. Firmicutes are the most abundant bacterial phylum present in P. japonica, followed by Proteobacteria and Actinobacteria. The most abundant genus was Staphylococcus. Second, the expression levels of the detoxification and digestion-related genes did not change significantly in any Cry protein treatment. Third, none of the Cry proteins affected the population fitness of P. japonica. These results indicated that P. japonica was not sensitive to Bt proteins, suggesting that growing Bt cotton expressing Cry1Ac and/or Cry2Ab will pose negligible risks to P. japonica.

Bacillus thuringiensis protein Vip3Aa does not harm the predator Propylea japonica: A toxicological, histopathological, biochemical and molecular analysis

The ladybeetle Propylea japonica is a widely distributed natural enemy in many agricultural systems. P. japonica is often used as a test organism for safety assessments of transgenic Bacillus thuringiensis crops. Plant varieties expressing the Vip3Aa insecticidal protein are not currently commercially available in China. In this study, protease inhibitor E-64 was used as a positive control to examine the responses of P. japonica larvae to a high concentration of Vip3Aa proteins. Larvae that were fed E-64 had increased mortality and prolonged developmental period, but these parameters were unaffected when larvae were fed Vip3Aa. The epithelial cells of midguts were intact and closely connected with the basal membrane when larvae were fed Vip3Aa, but the epithelial cells degenerated in the E-64 treatment. The activities of antioxidative enzymes and expression levels of detoxification-related genes in P. japonica larvae were not altered after exposure to Vip3Aa; however, these biochemical and molecular parameters were significantly changed in the E-64 treatment. The results demonstrate that Vip3Aa protein is not harmful to the predator P. japonica.

Evaluation of a cuticle protein gene as a potential RNAi target in aphids

BACKGROUND

RNA interference (RNAi) has potential as a pest insect control technique. One possible RNAi target is the cuticle protein, which is important in insect molting and development. As an example, here we evaluate the possibility of designing double-stranded RNA (RNA) that is effective for silencing the cuticle protein 19 gene (CP19) in aphids but is harmless to non-target predator insects.

RESULTS

The sequences of CP19s were similar (86.6-94.4%) among the tested aphid species (Aphis citricidus, Acyrthosiphon pisum, and Myzus persicae) but different in the predator Propylaea japonica. Ingestion of species-specific dsRNAs of CP19 by the three aphids produced 39.3-64.2% gene silencing and 45.8-55.8% mortality. Ingestion of non-species-specific dsRNA (dsAcCP19) by Ac. pisum and M. persicae gave gene silencing levels ranging from 40.4% to 50.3% and 43.3-50.8% mortality. The dsApCP19 did not affect PjCP19 expression or developmental duration in P. japonica.

CONCLUSION

The results demonstrate that CP19 is a promising RNAi target for aphid control via one dsRNA design. The targeting of genes that are conserved in insect pests but not present in beneficial insects is a useful RNAi-based pest control strategy. © 2019 Society of Chemical Industry.

Selection of appropriate reference genes for RT-qPCR analysis in Propylea japonica (Coleoptera: Coccinellidae)

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique commonly used in molecular biology to analyze RNA expression. The selection of suitable reference genes for data normalization is a precondition for credible measurements of gene expression levels using RT-qPCR. Propylea japonica is one of the most common pests of many crop systems throughout East Asia, and has often been used in the testing of non-target impacts during environmental risk assessments of genetically engineered plants. The present study assessed the suitability of nine frequently used reference genes for comparisons of P. japonica gene expression. Expression stability was compared across developmental stages, sex, a range of tissues, and following exposure to different temperatures. Data were analyzed using RefFinder, which integrated the results obtained using NormFinder, geNorm, BestKeeper, and the ΔCt method. This led to the identification of unique sets of reference genes for each experimental condition: ribosomal protein S18 (RPS18) and elongation factor 1 α (EF1A) for developmental stage comparisons, RPS18 and EF1A for sex comparisons, EF1A and ribosomal protein L4 for tissue comparisons, and RPS18 and EF1A for analyses of temperature-mediated effects. These reference genes will help to enhance the accuracy of RT-qPCR analyses of P. japonica gene expression. This work represents an initial move towards building a standardized system for RT-qPCR analysis of P. japonica, providing a basis for the ecological risk assessment of RNAi-based insect control products.

Vitality and Stability of Insecticide Resistance in Adult Propylaea japonica (Coleoptera: Coccinellidae)

Propylaea japonica (Thunberg) was a dominant species among the predacious ladybirds in the fields and active from March to November during a year in Fuzhou, China. Stability of insecticide resistance and vitality in adult P. japonica were investigated. The field ladybird P. japonica in Fuzhou, China, showed 9- to 16-fold resistance ratios to chlorpyrifos, 13- to 2,083-fold to methamidophos, 32- to 230-fold to fenvalerate, and 4- to 49-fold to avermectins, respectively, based on the field monitoring during 2004, 2009, and 2012, as compared with insecticide-susceptible F39 progenies. The resistance levels in the field P. japonica were high during May and November but low during summer. The population growth tendency index in field P. japonica was 0.8-fold as high as that in insecticide-susceptible P. japonica. The field P. japonica also showed high tolerance to the insecticide as compared with pest Lipaphis erysimi Kaltenbach and two parasitoids Diaeretiella rapae and Pachyneuron aphidis. Stable insecticide resistance levels and high vitality were found first in adult P. japonica with 1-, 30-, or 60-d-old adults, or among the adults developed form the eggs produced by newly emerged adults or by 60-d-old adults. Increased activity of glutathione S-transferases, carboxylesterases, and cytochrome P450 monooxygenases might be involved in the resistance of P. japonica. The results indicated that, in certain areas, inclusion of P. japonica for pest control in the integrated pest management would be highly recommended.

Prey Foraging Under Sublethal Lambda-Cyhalothrin Exposure on Pyrethroid-Susceptible and -Resistant Lady Beetles (Eriopis connexa (Coleoptera: Coccinelidae))

Sublethal insecticide exposure may affect foraging of insects, including natural enemies, although the subject is usually neglected. The lady beetle Eriopis connexa (Germar, 1824) (Coleoptera: Coccinelidae) is an important predator of aphids with existing pyrethroid-resistant populations that are undergoing scrutiny for potential use in pest management systems characterized by frequent insecticide use. However, the potential effect of sublethal pyrethroid exposure on this predator's foraging activity has not yet been assessed and may compromise its use in biological control. Therefore, our objective was to assess the effect of sublethal lambda-cyhalothrin exposure on three components of the prey foraging activity (i.e., walking, and prey searching and handling), in both pyrethroid-susceptible and -resistant adults of E. connexa. Both lady beetle populations exhibited similar walking patterns without insecticide exposure in noncontaminated arenas, but in partially contaminated arenas walking differed between strains, such that the resistant insects exhibited greater walking activity. Behavioral avoidance expressed as repellence to lambda-cyhalothrin was not observed for either the susceptible or resistant populations of E. connexa, but the insecticide caused avoidance by means of inducing irritability in 40% of the individuals, irrespective of the strain. Insects remained in the insecticide-contaminated portion of the arena for extended periods resulting in greater exposure. Although lambda-cyhalothrin exposure did not affect prey searching by susceptible lady beetles, prey searching was extended for exposed resistant predators. In contrast, prey handling was not affected by population or by lambda-cyhalothrin exposure. Thus, sublethal exposure to the insecticide in conjunction with the insect resistance profile can affect prey foraging with pyrethroid-exposed resistant predators exhibiting longer prey searching time associated with higher walking activity reducing its predatory performance.

Leukemia inhibitory factor regulates marmoset induced pluripotent stem cell proliferation via a PI3K/‌Akt‑dependent Tbx‑3 activation pathway

Leukemia inhibitory factor (LIF) is the most pleiotropic cytokine of the interleukin-6 family, and is widely used to establish and maintain pluripotent stem cells, particularly mouse pluripotent stem cells. However, no reports have fully elucidated the application of LIF in marmoset induced pluripotent stem cell (iPSC) culture, particularly the underlying mechanisms. To demonstrate the feasibility of the application of LIF to marmoset iPSCs, the present study assessed these cells in the presence of LIF. Cell proliferation was measured using MTT assay, cell apoptosis was determined by flow cytometric analysis of fluorescein isothiocyanate Annexin V staining and the differentially expressed genes were analysed using Digital Gene Expression (DGE) analysis. The altered expression of pluripotency-associated genes was confirmed by reverse transcription-quantitative polymerase chain reaction and western blot analysis. Furthermore, following treatment with LY294002, cell proliferation was measured by MTT assay and protein levels were confirmed by western blot analysis. The results showed that LIF significantly promoted the number of proliferating cells, but had no effect on apoptosis. Digital Gene Expression analysis was used to examine the differentially expressed genes of marmoset iPSCs in the presence of LIF. The results showed that the pluripotency-associated transcription factor-encoding gene T-box 3 (Tbx-3) was activated by LIF. Notably, LIF increased the levels of phosphorylated (p-)AKT and Tbx-3 in the marmoset iPSCs. Furthermore, pretreatment with LY294002, an inhibitor of phosphoinositide 3-kinase (PI3K), significantly impaired the LIF-induced upregulation of p-AKT and Tbx-3 in the marmoset iPSCs, suggesting that the PI3K/Akt signaling pathway is involved in this regulation. Taken together, the results suggested that LIF is effective in maintaining marmoset iPSCs in cultures, which is associated with the activation of Tbx-3 through regulation of the PI3K/Akt signaling pathway.

Identification and characterization of ace1-type acetylcholinesterase in insecticide-resistant and -susceptible Propylaea japonica (Thunberg)

Characterization and gene cloning of acetylecholinesterase (AChE) in the insecticide-resistant (R) and -susceptible (S) insects have been reported in the past. However, the studies focused mostly on herbivorous pests, rather than predacious species, such as ladybird beetles. Using R and S Propylaea japonica (thunberg), a full-length cDNA sequence (2928 bp) of the ace1-type AChE gene was determined for the first time. The ace1 encoding a protein of 645 amino acids contained typical conserved motifs, such as FGESAG domains, catalytic triad, acyl pocket, oxyanino hole, choline binding site, peripheral anionic site, omega loop and conserved aromatic residues. R P. japonica displayed 50-times greater resistance to chlorpyrifos or mathamidophos with a significantly lower AChE sensitivity to paraoxon, malaoxon, chlorpyrifos or methamidophos than its S counterpart. Five amino acids in the ace1 of R P. japonica differed from those found in S P. japonica. One of them, F358S, located in the acyl-binding pocket, might play a crucial role in the resistance of the insect to organophosphates (OPs). Whereas, K493E and I538V, which were close to some of the conserved aromatic amino acids (i.e., H509, Y511, and W499) in the gorge, and G571R and T576A near C593 that formed the disulfide bonds with C471, might also involve in the change of insecticide resistance in P. japonica. AChE insensitivity and amino acid replacements, particularly F358S, might be the determining factors in the alteration of OPs-resistance in P. japonica.

Adaptation to toxic hosts as a factor in the evolution of insecticide resistance

Insecticide resistance is a serious economic problem that jeopardizes sustainability of chemical control of herbivorous insects and related arthropods. It can be viewed as a specific case of adaptation to toxic chemicals, which has been driven in large part, but not exclusively, by the necessity for insect pests to tolerate defensive compounds produced by their host plants. Synthetic insecticides may simply change expression of specific sets of detoxification genes that have evolved due to ancestral associations with host plants. Feeding on host plants with more abundant or novel secondary metabolites has even been shown to prime insect herbivores to tolerate pesticides. Clear understanding of basic evolutionary processes is important for achieving lasting success in managing herbivorous arthropods.

Transcriptional responses of detoxification genes to four plant allelochemicals in Aphis gossypii

Aphis gossypii Glover (Hemiptera: Aphididae) can damage a variety of agricultural crops, so it is very important for cotton aphids to evolve adaptive mechanisms to various allelochemicals from host plants. Our results aim to provide a fundamental and rich resource for exploring aphid functional genes in A. gossypii. A transcriptome data set and five expression profile data sets of A. gossypii samples were analyzed by Illumina sequencing platform. In total, 53,763,866 reads were assembled into 1,963,516 contigs and 28,555 unigenes. Compared with the control, 619 genes were significantly up- or downregulated in the treatment group by 2-tridecanone. There were 516, 509, and 717 of differential expression genes in tannic acid, quercetin, and gossypol treatment groups, respectively. Furthermore, there were 4 of 54 putative cytochrome P450 genes and 1 of 7 putative carboxylesterases downregulated in all treatment groups by four plant allelochemicals. When aphids fed on 2-tridecanone, tannic acid, and quercetin, only one P450 gene was upregulated. These results show that plant allelochemical stress can induce differential gene expression in A. gossypii. The differential response information of gene expression based on a large-scale sequence would be useful to reveal molecular mechanisms of adaptation for A. gossypii to plant allelochemicals.

Modification of Flight and Locomotion Performances, Respiratory Metabolism, and Transcriptome Expression in the Lady Beetle Harmonia axyridis through Sublethal Pesticide Exposure

Biological control is usually used in combination with chemical control for practical agricultural applications. Thus, the influence of insecticides on the natural predators used for biological control should be investigated for integrated pest management. The ladybird Harmonia axyridis is an effective predator on aphids and coccids. Beta-cypermethrin is a broad-spectrum insecticide used worldwide for controlling insect pests. H. axyridis is becoming increasingly threatened by this insecticide. Here, we investigated the effect of a sublethal dose of beta-cypermethrin on flight, locomotion, respiration, and detoxification system of H. axyridis. After exposure to beta-cypermethrin, succinic female adults flew more times, longer distances, and during longer time periods. Exposure to a sublethal dose of beta-cypermethrin also promoted an increase in walking rate, walking distance, walking duration, and also an increase in respiratory quotient and respiratory rate. To investigate the effects of beta-cypermethrin on H. axyridis detoxification system, we analyzed the transcriptome of H. axyridis adults, focusing on genes related to detoxification systems. De novo assembly generated 65,509 unigenes with a mean length of 799 bp. From these genes, 26,020 unigenes (40.91% of all unigenes) exhibited clear homology to known genes in the NCBI non-redundant database. In addition, 10,402 unigenes were annotated in the Cluster of Orthologous Groups database, 12,088 unigenes were assigned to the Gene Ontology database and 12,269 unigenes were in the Kyoto Encyclopedia of Genes and Genome (KEGG) database. Exposure to beta-cypermethrin had significant effects on the transcriptome profile of H. axyridis adult. Based on uniquely mapped reads, 3,296 unigenes were differentially expressed, 868 unigenes were up-regulated and 2,248 unigenes were down-regulated. We identified differentially-expressed unigenes related to general detoxification systems in H. axyridis. This assembled, annotated transcriptome provides a valuable genomic resource for further understanding the molecular basis of detoxification mechanisms in H. axyridis.

De novo transcriptome and expression profile analyses of the Asian corn borer (Ostrinia furnacalis) reveals relevant flubendiamide response genes

Background

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), has become the most damaging insect pest of corn in Asia. However, the lack of genome or transcriptome information heavily hinders our further understanding of ACB in every aspect at a molecular level and on a genome-wide scale. Here, we used the Ion Torrent Personal Genome Machine (PGM) Sequencer to explore the ACB transcriptome and to identify relevant genes in response to flubendiamide, showing high selective activity against ACB.

Results

We obtained 35,430 unigenes, with an average length of 716 bp, representing a dramatic expansion of existing cDNA sequences available for ACB. These sequences were annotated with Non-redundant Protein (Nr), Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to better understand their functions. A total of 31 cytochrome P450 monooxygenases (P450s), 27 carboxyl/cholinesterases (CCEs) and 19 glutathione S-transferases (GSTs) were manually curated to construct phylogenetic trees, and 25 unigenes encoding target proteins (acetylcholinesterase, nicotinic acetylcholine receptor, gamma-aminobutyric acid receptor, glutamate-gated chloride channel, voltage-gated sodium channel and ryanodine receptor) were identified. In addition, we compared and validated the differentially expressed unigenes upon flubendiamide treatment, revealing that the genes for detoxification enzymes (P450s and esterase), calcium signaling pathways and muscle control pathways (twitchin and tropomyosin), immunoglobulin (hemolin), chemosensory protein and heat shock protein 70 were significantly overexpressed in response to flubendiamide, while the genes for cuticular protein, protease and oxidoreductase showed much lower expression levels.

Conclusion

The obtained transcriptome information provides large genomic resources available for further studies of ACB. The differentially expressed gene data will elucidate the molecular mechanisms of ACB in response to the novel diamide insecticide, flubendiamide. In particular, these findings will facilitate the identification of the genes involved in insecticide resistance and the development of new compounds to control the ACB.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3431-6) contains supplementary material, which is available to authorized users.

Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, Toxoptera citricida

Winged and wingless morphs in insects represent a trade-off between dispersal ability and reproduction. We studied key genes associated with apterous and alate morphs in Toxoptera citricida (Kirkaldy) using RNAseq, digital gene expression (DGE) profiling, and RNA interference. The de novo assembly of the transcriptome was obtained through Illumina short-read sequencing technology. A total of 44,199 unigenes were generated and 27,640 were annotated. The transcriptomic differences between alate and apterous adults indicated that 279 unigenes were highly expressed in alate adults, whereas 5,470 were expressed at low levels. Expression patterns of the top 10 highly expressed genes in alate adults agreed with wing bud development trends. Silencing of the lipid synthesis and degradation gene (3-ketoacyl-CoA thiolase, mitochondrial-like) and glycogen genes (Phosphoenolpyruvate carboxykinase [GTP]-like and Glycogen phosphorylase-like isoform 2) resulted in underdeveloped wings. This suggests that both lipid and glycogen metabolism provide energy for aphid wing development. The large number of sequences and expression data produced from the transcriptome and DGE sequencing, respectively, increases our understanding of wing development mechanisms.

Variation in life history traits and transcriptome associated with adaptation to diet shifts in the ladybird Cryptolaemus montrouzieri

BACKGROUND

Despite the broad diet range of many predatory ladybirds, the mechanisms involved in their adaptation to diet shifts are not completely understood. Here, we explored how a primarily coccidophagous ladybird Cryptolaemus montrouzieri adapts to feeding on aphids.

RESULTS

Based on the lower survival rate, longer developmental time, and lower adult body weight and reproduction rate of the predator, the aphid Megoura japonica proved being less suitable to support C. montrouzieri as compared with the citrus mealybug Planococcus citri. The results indicated up-regulation of genes related to ribosome and translation in fourth instars, which may be related to their suboptimal development. Also, several genes related to biochemical transport and metabolism, and detoxification were up-regulated as a result of adaptation to the changes in nutritional and non-nutritional (toxic) components of the prey.

CONCLUSION

Our results indicated that C. montrouzieri succeeded in feeding on aphids by regulation of genes related to development, digestion and detoxification. Thus, we argue that these candidate genes are valuable for further studies of the functional evolution of ladybirds led by diet shifts.

Bt proteins Cry1Ah and Cry2Ab do not affect cotton aphid Aphis gossypii and ladybeetle Propylea japonica

Plant varieties expressing the Bt (Bacillus thuringiensis) insecticidal proteins Cry1Ah and Cry2Ab have potential commercialization prospects in China. However, their potential effects on non-target arthropods (NTAs) remain uncharacterized. The cotton aphid Aphis gossypii is a worldwide pest that damages various important crops. The ladybeetle Propylea japonica is a common and abundant natural enemy in many cropping systems in East Asia. In the present study, the effects of Cry1Ah and Cry2Ab proteins on A. gossypii and P. japonica were assessed from three aspects. First, neither of the Cry proteins affected the growth or developmental characteristics of the two test insects. Second, the expression levels of the detoxification-related genes of the two test insects did not change significantly in either Cry protein treatment. Third, neither of the Cry proteins had a favourable effect on the expression of genes associated with the amino acid metabolism of A. gossypii and the nutrition utilization of P. japonica. In conclusion, the Cry1Ah and Cry2Ab proteins do not appear to affect the cotton aphid A. gossypii or the ladybeetle P. japonica.

De novo transcriptome sequencing and analysis of Coccinella septempunctata L. in non-diapause, diapause and diapause-terminated states to identify diapause-associated genes

BACKGROUND

The most common ladybird beetle, Coccinella septempunctata L., is an excellent predator of crop pests such as aphids and white flies, and it shows a wide range of adaptability, a large appetite and a high reproductive ability. Diapause research plays an important role in the artificial propagation and shelf-life extension of insect products. Although this lady beetle's regulatory, physiological and biochemical characteristics in the diapause period are well understood, the molecular mechanism of diapause remains unknown. Therefore, we collected female adults in three different states, i.e., non-diapause, diapause and diapause termination, for transcriptome sequencing.

RESULTS

After transcriptome sequencing using the Illumina HiSeq 2500 platform with pretreatment, a total of 417.6 million clean reads from nine samples were filtered using the program FASTX (version 0.0). Additionally, 106,262 contigs were assembled into 82,820 unigenes with an average length of 921 bp and an N50 of 1,241 bp. All of the unigenes were annotated through BLASTX alignment against the Nr or UniProt database, and 37,872 unigenes were matched. We performed further analysis of these unigenes using the Clusters of Orthologous Groups of proteins (COG), Gene Ontology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Through pairwise comparisons of the non-diapause (ND), diapause (D), and diapause-terminated (DT) groups, 3,501 and 1,427 differentially expressed genes (DEGs) were identified between D and ND and between DT and D, respectively. Moreover, 443 of the DEGs were specifically expressed during the diapause period (i.e., DEGs that were expressed at the highest or lowest levels during diapause compared with the other stages). GO function and KEGG pathway enrichment were performed on all DEGs and showed that RNA-directed DNA polymerase activity and fatty acid metabolism were significantly affected. Furthermore, eight specific expressed genes were selected for validation using qRT-PCR. Among these eight genes, seven genes were up-regulated, and one gene was down-regulated; the change trends of the eight genes were the same between the qRT-PCR and RNA-seq analysis results.

CONCLUSIONS

In this study, a new method for collecting and identifying diapause insects was described. We generated a vast quantity of transcriptome data from C. septempunctata L., providing a resource for gene function research. The diapause-associated genes that we identified establish a foundation for future studies on the molecular mechanisms of diapause.

Genome-Wide Expression Profiling of Anoxia/Reoxygenation in Rat Cardiomyocytes Uncovers the Role of MitoKATP in Energy Homeostasis

Mitochondrial ATP-sensitive potassium channel (mitoK(ATP)) is a common end effector of many protective stimuli in myocardial ischemia-reperfusion injury (MIRI). However, the specific molecular mechanism underlying its myocardial protective effect is not well elucidated. We characterized an anoxia/reoxygenation (A/R) model using freshly isolated adult rat cardiomyocytes. MitoK(ATP) status was interfered with its specific opener diazoxide (DZ) or blocker 5-hydroxydecanote (5-HD). Digital gene expression (DGE) and bioinformatic analysis were deployed. Three energy metabolism related genes (MT-ND6, Idh2, and Acadl) were upregulated when mitoK(ATP) opened. In addition, as many as 20 differentially expressed genes (DEGs) were significantly enriched in five energy homeostasis correlated pathways (PPAR, TCA cycle, fatty acid metabolism, and peroxisome). These findings indicated that mitoK(ATP) opening in MIRI resulted in energy mobilization, which was confirmed by measuring ATP content in cardiomyocytes. These causal outcomes could be a molecular mechanism of myocardial protection of mitoKATP and suggested that the mitoK(ATP) opening plays a physiologic role in triggering cardiomyocytes' energy homeostasis during MIRI. Strategies of modulating energy expenditure during myocardial ischemia-reperfusion may be promising approaches to reduce MIRI.