Transcriptome differences between Cry1Ab resistant and susceptible strains of Asian corn borer

MiR-2b-3p Downregulated PxTrypsin-9 Expression in the Larval Midgut to Decrease Cry1Ac Susceptibility of the Diamondback Moth, Plutella xylostella (L.)

Crystal (Cry) toxins, produced by Bacillus thuringiensis, are widely used as effective biological pesticides in agricultural production. However, insects always quickly evolve adaptations against Cry toxins within a few generations. In this study, we focused on the Cry1Ac protoxin activated by protease. Our results identified PxTrypsin-9 as a trypsin gene that plays a key role in Cry1Ac virulence in Plutella xylostella larvae. In addition, P. xylostella miR-2b-3p, a member of the micoRNA-2 (miR-2) family, was significantly upregulated by Cry1Ac protoxin and targeted to PxTrypsin-9 downregulated its expression. The mRNA level of PxTrypsin-9, regulated by miR-2b-3p, revealed an increased tolerance of P. xylostella larvae to Cry1Ac at the post-transcriptional level. Considering that miR-2b and trypsin genes are widely distributed in various pest species, our study provides the basis for further investigation of the roles of miRNAs in the regulation of the resistance to Cry1Ac and other insecticides.

CeRNA network identified hsa-miR-17-5p, hsa-miR-106a-5p and hsa-miR-2355-5p as potential diagnostic biomarkers for tuberculosis

This study aims to analyze the regulatory non-coding RNAs in the pathological process of tuberculosis (TB), and identify novel diagnostic biomarkers. A longitudinal study was conducted in 5 newly diagnosed pulmonary tuberculosis patients, peripheral blood samples were collected before and after anti-TB treatment for 6 months, separately. After whole transcriptome sequencing, the differentially expressed RNAs (DE RNAs) were filtrated with |log2 (fold change) | > log2(1.5) and P value < .05 as screening criteria. Then functional annotation was actualized by gene ontology enrichment analysis, and enrichment pathway analysis was conducted by Kyoto Encyclopedia of Genes and Genomes database. And finally, the competitive endogenous RNA (ceRNA) regulatory network was established according to the interaction of ceRNA pairs and miRNA-mRNA pairs. Five young women were recruited and completed this study. Based on the differential expression analysis, a total of 1469 mRNAs, 996 long non-coding RNAs, 468 circular RNAs, and 86 miRNAs were filtrated as DE RNAs. Functional annotation demonstrated that those DE-mRNAs were strongly involved in the cellular process (n = 624), metabolic process (n = 513), single-organism process (n = 505), cell (n = 651), cell part (n = 650), organelle (n = 569), and binding (n = 629). Enrichment pathway analysis revealed that the differentially expressed genes were mainly enriched in HTLV-l infection, T cell receptor signaling pathway, glycosaminoglycan biosynthesis-heparan sulfate/heparin, and Hippo signaling pathway. CeRNA networks revealed that hsa-miR-17-5p, hsa-miR-106a-5p and hsa-miR-2355-5p might be regarded as potential diagnostic biomarkers for TB. Immunomodulation-related genes are differentially expressed in TB patients, and hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-2355-5p might serve as potential diagnostic biomarkers.

Critical Analysis of Multi-Omic Data from a Strain of Plutella xylostella Resistant to Bacillus thuringiensis Cry1Ac Toxin

Rapid evolution of resistance in crop pests to Bacillus thuringiensis (Bt) products threatens their widespread use, especially as pests appear to develop resistance through a range of different physiological adaptations. With such a diverse range of mechanisms reported, researchers have resorted to multi-omic approaches to understand the molecular basis of resistance. Such approaches generate a lot of data making it difficult to establish where causal links between physiological changes and resistance exist. In this study, a combination of RNA-Seq and iTRAQ was used with a strain of diamondback moth, Plutella xylostella (L.), whose resistance mechanism is well understood. While some of the causal molecular changes in the resistant strain were detected, other previously verified changes were not detected. We suggest that while multi-omic studies have use in validating a proposed resistance mechanism, they are of limited value in identifying such a mechanism in the first place.

Toxicity of Bacillus thuringiensis Strains Derived from the Novel Crystal Protein Cry31Aa with High Nematicidal Activity against Rice Parasitic Nematode Aphelenchoides besseyi

The plant parasitic nematode, Aphelenchoides besseyi, is a serious pest causing severe damage to various crop plants and vegetables. The Bacillus thuringiensis (Bt) strains, GBAC46 and NMTD81, and the biological strain, FZB42, showed higher nematicidal activity against A. besseyi, by up to 88.80, 82.65, and 75.87%, respectively, in a 96-well plate experiment. We screened the whole genomes of the selected strains by protein-nucleic acid alignment. It was found that the Bt strain GBAC46 showed three novel crystal proteins, namely, Cry31Aa, Cry73Aa, and Cry40ORF, which likely provide for the safe control of nematodes. The Cry31Aa protein was composed of 802 amino acids with a molecular weight of 90.257 kDa and contained a conserved delta-endotoxin insecticidal domain. The Cry31Aa exhibited significant nematicidal activity against A. besseyi with a lethal concentration (LC₅₀) value of 131.80 μg/mL. Furthermore, the results of in vitro experiments (i.e., rhodamine and propidium iodide (PI) experiments) revealed that the Cry31Aa protein was taken up by A. besseyi, which caused damage to the nematode's intestinal cell membrane, indicating that the Cry31Aa produced a pore-formation toxin. In pot experiments, the selected strains GBAC46, NMTD81, and FZB42 significantly reduced the lesions on leaves by up to 33.56%, 45.66, and 30.34% and also enhanced physiological growth parameters such as root length (65.10, 50.65, and 55.60%), shoot length (68.10, 55.60, and 59.45%), and plant fresh weight (60.71, 56.45, and 55.65%), respectively. The number of nematodes obtained from the plants treated with the selected strains (i.e., GBAC46, NMTD81, and FZB42) and A. besseyi was significantly reduced, with 0.56, 0.83., 1.11, and 5.04 seedling mL^-1 nematodes were achieved, respectively. Moreover, the qRT-PCR analysis showed that the defense-related genes were upregulated, and the activity of hydrogen peroxide (H₂O₂) increased while malondialdehyde (MDA) decreased in rice leaves compared to the control. Therefore, it was concluded that the Bt strains GBAC46 and NMTD81 can promote rice growth, induce high expression of rice defense-related genes, and activate systemic resistance in rice. More importantly, the application of the novel Cry31Aa protein has high potential for the efficient and safe prevention and green control of plant parasitic nematodes.

CRISPR/Cas9-Mediated Mutagenesis of Abdominal-A and Ultrabithorax in the Asian Corn Borer, Ostrinia furnacalis

Simple Summary

Homeotic genes encode transcription factors that coordinated the anatomical structure formation during the early embryonic development of organisms. In this study, we functionally characterized two homeotic genes, Abdominal-A (Abd-A) and Ultrabithorax (Ubx), in the Asian corn borer, Ostrinia furnacalis (a maize pest that has devastated the Asia-Pacific region) by using a CRISPR/Cas9 genome editing system. Our results show that the mutagenesis of OfAbd-A and OfUbx led to severe morphological defects in O. furnacalis, which included fused segments and segmental twist during the larval stage, and hollowed and incision-like segments during the pupal stage in OfAbd-A mutants, as well as defects in the wing-pad development in pupal and adult OfUbx mutants. Overall, knocking out Abd-A and Ubx in O. furnacalis resulted in the embryonic lethality to, and pleiotropic impact on, other homeotic genes. This study not only confirms the conserved body planning functions in OfAbd-A and OfUbx, but it also strengthens the control implications of these homeotic genes for lepidopteran pests.

Abstract

(1) Background: Abdominal-A (Abd-A) and Ultrabithorax (Ubx) are homeotic genes that determine the identity and morphology of the thorax and abdomen in insects. The Asian corn borer, Ostrinia furnacalis (Guenée) (Lepidoptera: Pyralidae), is a devastating maize pest throughout Asia, the Western Pacific, and Australia. Building on previous knowledge, we hypothesized that the knockout of Abd-A and Ubx would disrupt the abdominal body planning in O. furnacalis. (2) Methods: CRISPR/Cas9-targeted mutagenesis was employed to decipher the functions of these homeotic genes. (3) Results: Knockout insects demonstrated classical homeotic transformations. Specifically, the mutagenesis of OfAbd-A resulted in: (1) Fused segments and segmental twist during the larval stage; (2) Embryonic lethality; and (3) The pleiotropic upregulation of other homeotic genes, including Lab, Pd, Dfd, Antp, and Abd-B. The mutagenesis of OfUbx led to: (1) Severe defects in the wing pads, which limited the ability of the adults to fly and mate; (2) Female sterility; and (3) The pleiotropic upregulation of other homeotic genes, including Dfd, Abd-B, and Wnt1. (4) Conclusions: These combined results not only support our hypothesis, but they also strengthen the potential of using homeotic genes as molecular targets for the genetic control of this global insect pest.

The stability and sequence cleavage preference of dsRNA are key factors differentiating RNAi efficiency between migratory locust and Asian corn borer

We compared the stability of double-stranded RNA (dsRNA) in each of two body fluids (hemolymph, midgut fluid) and in each of two tissues (integument, midgut), and the uptake of dsRNA in each of two cultured tissues (integument, midgut) between the migratory locust (Locusta migratoria) and the Asian corn borer (Ostrinia furnacalis). We further compared the abundance of putative small interfering RNAs (siRNAs) generated from each of two dsRNAs (dsβ-actin, dsEf1α) and the preference of dsRNA cleavages between the two insect species. Our studies showed a rapid degradation of dsRNA in the midgut fluids of both insect species and in O. furnacalis hemolymph. However, dsRNA remained reasonably stable in L. migratoria hemolymph. When nuclease degradation of dsRNA in cultured tissues was inhibited, dsRNA uptake was not significantly different between the two species. We further showed that the silencing efficiency against target genes was consistent with the abundance of putative siRNAs processed from the dsRNA. In addition, O. furnacalis showed a strong preference in cleaving dsRNA when the nucleotide G was in the position of "1" at 5'-end whereas L. migratoria showed broad spectrum in cleavage sites to generate siRNA. Taken together, our study revealed that silencing efficiency of a target gene by RNAi was directly related to the dsRNA degradation by nucleases and the abundance of siRNAs generated from the dsRNA.

Increased RNAi Efficiency by dsEGFP-Induced Up-Regulation of Two Core RNAi Pathway Genes (OfDicer2 and OfAgo2) in the Asian Corn Borer (Ostrinia furnacalis)

Simple Summary

RNA interference (RNAi) has shown great potentials as a novel technology for insect pest management. However, numerous studies have shown that the efficiency of RNAi varies substantially among different insect species. For example, as a major insect pest of corn, the Asian corn borer (Ostrinia furnacalis) showed very low RNAi efficiency. Therefore, it is necessary to develop new strategies for enhancing RNAi efficiency in insects with low RNAi efficiency. In this study, six core RNAi pathway genes were identified and characterized from O. furnacalis transcriptome database. After dsEGFP was injected into O. furnacalis, the expression of the core RNAi pathway genes (OfDicer2 and OfAgo2) was significantly up-regulated in response to the exposure of dsEGFP. As a result, the RNAi efficiency against the target genes in certain tissues of O. furnacalis was significantly improved. These results suggest that RNAi efficiency can be improved by inducing the expression of key RNAi pathway genes in O. furnacalis.

Abstract

RNA interference (RNAi) is a sequence-specific gene silencing mechanism that holds great promise for effective management of agricultural pests. Previous studies have shown that the efficacy of RNAi varies among different insect species, which limits its wide spread application in the field of crop protection. In this study, we identified and characterized six core RNAi pathway genes including OfDicer1, OfDicer2, OfR2D2, OfAgo1, OfAgo2, and OfAgo3 from the transcriptomic database of the Asian corn borer (Ostrinia furnacalis). Domain analysis showed that the six deduced proteins contained the necessary functional domains. Insect developmental stage- and tissue-specific expression analysis showed that five genes were expressed in all the stages and tissues examined except OfAgo3, which showed low expression in larvae, and high expression in pupae and adults and in the midgut. RT-qPCR was performed to examine the response of these six genes to exogenous double-stranded RNA (dsRNA). Interestingly, the transcript levels of OfDicer2 and OfAgo2 were significantly enhanced after the injection of dsEGFP at different time points and tissues investigated. Consequently, the RNAi efficiency in targeting the insect endogenous genes can be greatly enhanced in the hemolymph or midgut. Taken together, our investigations suggest that RNAi efficiency can be enhanced by pre-injection of dsRNA to induce the RNAi core machinery genes, which could be a useful strategy to improving RNAi efficiency for studying gene functions under laboratory conditions.

A Novel Reference for Bt-Resistance Mechanism in Plutella xylostella Based on Analysis of the Midgut Transcriptomes

Simple Summary

Plutella xylostella is a very serious pest to cruciferous vegetables. At present, the control methods used are mainly traditional insecticides and the cultivation of Bt crops. However, with the long-term and large-scale use of insecticides, the diamondback moth has developed strong resistance to many kinds of insecticides and Bt crops. The Cry1S1000 strain of P. xylostella used here is a strain with more than 8000 times resistance to Bt Cry1Ac protoxin. In this paper, we used transcriptome sequencing to determine the midgut transcriptome of the G88-susceptible strain, Cry1S1000-resistant strain and its corresponding toxin-induced strains to find more genes related to Bt resistance. Our results will provide a reference for optimizing the control strategy of diamondback moth resistance and improving the control efficiency of biopesticides and Bt crops.

Abstract

The diamondback moth, Plutella xylostella, is a lepidopteran insect that mainly harms cruciferous vegetables, with strong resistance to a variety of agrochemicals, including Bacillus thuringiensis (Bt) toxins. This study intended to screen genes associated with Bt resistance in P. xylostella by comparing the midgut transcriptome of Cry1Ac-susceptible and -resistant strains together with two toxin-treated strains 24 h before sampling. A total of 12 samples were analyzed by BGISEQ-500, and each sample obtained an average of 6.35 Gb data. Additionally, 3284 differentially expressed genes (DEGs) were identified in susceptible and resistant strains. Among them, five DEGs for cadherin, 14 for aminopeptidase, zero for alkaline phosphatase, 14 for ATP binding cassette transport, and five heat shock proteins were potentially involved in resistance to Cry1Ac in P. xylostella. Furthermore, DEGs associated with “binding”, “catalytic activity”, “cellular process”, “metabolic process”, and “cellular anatomical entity” were more likely to be responsible for resistance to Bt toxin. Thus, together with other omics data, our results will offer prospective genes for the development of Bt resistance, thereby providing a brand new reference for revealing the resistance mechanism to Bt of P. xylostella.

The draft genome of the Asian corn borer yields insights into ecological adaptation of a devastating maize pest

The Asian corn borer (ACB) is the most devastating pest on maize in the western Pacific region of Asia. Despite broad interests in insecticide resistance, seasonal adaptation, and larval color mimicry regarding the ACB system, lacking of reference genomic information and a powerful gene editing approach have hindered the in-depth studies of these aspects. Here we present a 455.7 Mb draft genome of ACB with 98.4% completeness. Comparative genomics analysis showed an evident expansion in gene families of gustatory receptors (105), which is related to polyphagous characteristics. Based on the comparative transcriptome analysis of resistant and susceptible ACB against Bt Cry1Ab toxin, we identified 26 genes related to Cry1Ab resistance. Additionally, transcriptomics of insects exposed to conditions of low temperature and diapause (LT) vs. room temperature and diapause (RT) provided insights into the genetic mechanisms of cold adaptation. We also successfully developed an efficient CRISPR/Cas9-based genome editing system and applied it to explore the role of color pattern genes in the ecological adaptation of ACB. Taken together, our study provides a fully annotated high-quality reference genome and efficient gene editing system to realize the potential of ACB as a study system to address important biological questions such as insecticide resistance, seasonal adaptation, and coloration. These valuable genomic resources will also benefit the development of novel strategies for maize pest management.

Toxicological and molecular profiling of insecticide resistance in a Brazilian strain of fall armyworm resistant to Bt Cry1 proteins

BACKGROUND

Spodoptera frugiperda, fall armyworm (FAW) is the major pest of maize in Brazil and has readily acquired field resistance to a broad range of synthetic insecticides and to Bacillus thuringiensis (Bt) insecticidal proteins expressed in important crops. This study aims to understand patterns of cross-resistance in FAW by investigating the toxicological profile of a Bt-resistant Brazilian strain (Sf_Des) in comparison to a Bt-susceptible strain (Sf_Bra).

RESULTS

Laboratory bioassays with 15 active substances of nine mode of action classes revealed that Sf_Des has a medium level of resistance to deltamethrin and chlorpyrifos. Very high cross-resistance was observed among Cry1 toxins, but high susceptibility against Vip3A. Strain Sf_Des exhibited - depending on the substrate - up to 19-fold increased cytochrome P450 activity in comparison to Sf_Bra. RNA-Seq data support a major role of P450 enzymes in the detoxification of insecticides because we detected 85 P450 transcripts upregulated in Sf_Des. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) analysis confirmed that CYP9A-like and CYP6B39 are significantly upregulated (>200-fold) in Sf_Des in comparison to Sf_Bra strain. No target-site mutation linked to pyrethroid resistance was detected, but mutations in the AChE linked to organophosphate resistance were observed in Sf_Des. A Gene Ontology (GO) analysis of differentially expressed genes (DEG) categorized most of them into the biological process category, involved in oxidation-reduction and metabolic processes.

CONCLUSION

Our results indicate that multiple/cross-resistance mechanisms may have developed in the Sf_Des strain to conventional insecticides and Bt insecticidal proteins. The systematic toxicological analysis presented will help to guide recommendations for an efficient resistance management. © 2020 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Transcriptome reveal the response to Cry1Ac toxin in susceptible Bombyx mori

Bombyx mori as a representative in Lepidoptera is an important economic insect in agriculture production. Bacillus thuringiensis (Bt) is a bacterial pathogen in silkworm production. Understanding how silkworm respond to Bt-toxin can provide guidance to cultivate resistant silkworm strains. Cry1Ac is one type of Bt-toxin. In current research, Dazao, a susceptible B. mori strain to Bt-toxin, was treated by Cry1Ac toxin and compared its transcriptome with untreated samples. This analysis detected 1234 differentially expressed genes (DEGs). Gene Ontology, KEGG, and UniProt keyword enrichment analysis showed that DEGs include ATP-binding cassette (ABC) transporter, stress response, cuticle, and protein synthesis, and folding process. Five ABC genes were upregulated after Cry1Ac treatment including ABCA2, ABCA3, and ABCC4. They are also known as the transporters of Bt-toxin in lepidopteran insect. Expression of cuticle proteins was significantly increased at 6 h after Cry1Ac treatment. Sex-specific storage-proteins and heat shock protein were also upregulated in Cry1Ac treated samples. Our data provide an expression profile about the response of Cry1Ac toxin in susceptible B. mori strain.

Which Is Stronger? A Continuing Battle Between Cry Toxins and Insects

In this article, we review the latest works on the insecticidal mechanisms of Bacillus thuringiensis Cry toxins and the resistance mechanisms of insects against Cry toxins. Currently, there are two models of insecticidal mechanisms for Cry toxins, namely, the sequential binding model and the signaling pathway model. In the sequential binding model, Cry toxins are activated to bind to their cognate receptors in the mid-intestinal epithelial cell membrane, such as the glycophosphatidylinositol (GPI)-anchored aminopeptidases-N (APNs), alkaline phosphatases (ALPs), cadherins, and ABC transporters, to form pores that elicit cell lysis, while in the signaling pathway model, the activated Cry toxins first bind to the cadherin receptor, triggering an extensive cell signaling cascade to induce cell apoptosis. However, these two models cannot seem to fully describe the complexity of the insecticidal process of Cry toxins, and new models are required. Regarding the resistance mechanism against Cry toxins, the main method insects employed is to reduce the effective binding of Cry toxins to their cognate cell membrane receptors by gene mutations, or to reduce the expression levels of the corresponding receptors by trans-regulation. Moreover, the epigenetic mechanisms, host intestinal microbiota, and detoxification enzymes also play significant roles in the insects' resistance against Cry toxins. Today, high-throughput sequencing technologies like transcriptomics, proteomics, and metagenomics are powerful weapons for studying the insecticidal mechanisms of Cry toxins and the resistance mechanisms of insects. We believe that this review shall shed some light on the interactions between Cry toxins and insects, which can further facilitate the development and utilization of Cry toxins.

Transcriptional profiling analysis of susceptible and resistant strains of Anticarsia gemmatalis and their response to Bacillus thuringiensis

Anticarsia gemmatalis is one of the main defoliators of soybean in Brazil. Bacillus thuringiensis (Bt) transgenic crops are used for their management. In this paper we used RNA-seq to explore the response of A. gemmatalis to Bt HD73, as well as to detect transcriptional differences after Bt infection between resistant and susceptible strains. A total of 3853 and 6224 differentially expressed genes (DGEs) were identified in susceptible and resistant larvae after Bt exposure, respectively. We identified 2143 DEGs between susceptible and resistant larvae and 1991 between susceptible and resistant larvae Bt exposed. Immunity-related genes, Bt toxins receptors, proteases, genes involved in metabolic processes, transporters, cuticle proteins and mobile elements have been identified. qRT-PCR data demonstrated upregulation of five genes in susceptible strain after Bt exposure. These results provide insights to understand the molecular and cellular mechanisms of response to Bt that could be used in strategies to control agricultural pests.

Assessing the Single and Combined Toxicity of Chlorantraniliprole and Bacillus thuringiensis (GO33A) against Four Selected Strains of Plutella xylostella (Lepidoptera: Plutellidae), and a Gene Expression Analysis

Concerns about resistance development to conventional insecticides in diamondback moth (DBM) Plutella xylostella (L.), the most destructive pest of Brassica vegetables, have stimulated interest in alternative pest management strategies. The toxicity of Bacillus thuringiensis subsp. aizawai (Bt GO33A) combined with chlorantraniliprole (Chl) has not been documented. Here, we examined single and combined toxicity of chlorantraniliprole and Bt to assess the levels of resistance in four DBM strains. Additionally, enzyme activities were tested in field-original highly resistant (FOH-DBM), Bt-resistant (Bt-DBM), chlorantraniliprole-resistant (CL-DBM), and Bt + chlorantraniliprole-resistant (BtC-DBM) strains. The Bt product had the highest toxicity to all four DBM strains followed by the mixture of insecticides (Bt + Chl) and chlorantraniliprole. Synergism between Bt and chlorantraniliprole was observed; the combination of Bt + (Bt + Chl) (1:1, LC50:LC50) was the most toxic, showing a synergistic effect against all four DBM strains with a poison ratio of 1.35, 1.29, 1.27, and 1.25. Glutathione S-transferase (GST) and carboxyl-esterase (CarE) activities showed positive correlations with chlorantraniliprole resistance, but no correlation was observed with resistance to Bt and Bt + Chl insecticides. Expression of genes coding for PxGST, CarE, AChE, and MFO using qRT-PCR showed that the PxGST and MFO were significantly overexpressed in Bt-DBM. However, AChE and CarE showed no difference in the four DBM strains. Mixtures of Bt with chlorantraniliprole exhibited synergistic effects and may aid the design of new combinations of pesticides to delay resistance in DBM strains substantially.

Response Mechanisms of Invertebrates to Bacillus thuringiensis and Its Pesticidal Proteins

Extensive use of chemical insecticides adversely affects both environment and human health. One of the most popular biological pest control alternatives is bioinsecticides based on Bacillus thuringiensis This entomopathogenic bacterium produces different protein types which are toxic to several insect, mite, and nematode species. Currently, insecticidal proteins belonging to the Cry and Vip3 groups are widely used to control insect pests both in formulated sprays and in transgenic crops. However, the benefits of B. thuringiensis-based products are threatened by insect resistance evolution. Numerous studies have highlighted that mutations in genes coding for surrogate receptors are responsible for conferring resistance to B. thuringiensis Nevertheless, other mechanisms may also contribute to the reduction of the effectiveness of B. thuringiensis-based products for managing insect pests and even to the acquisition of resistance. Here, we review the relevant literature reporting how invertebrates (mainly insects and Caenorhabditis elegans) respond to exposure to B. thuringiensis as either whole bacteria, spores, and/or its pesticidal proteins.

Differences in midgut transcriptomes between resistant and susceptible strains of Chilo suppressalis to Cry1C toxin

Background

Chilo suppressalis is a widespread rice pest that poses a major threat to food security in China. This pest can develop resistance to Cry toxins from Bacillus thuringiensis (Bt), threatening the sustainable use of insect-resistant transgenic Bt rice. However, the molecular basis for the resistance mechanisms of C. suppressalis to Cry1C toxin remains unknown. This study aimed to identify genes associated with the mechanism of Cry1C resistance in C. suppressalis by comparing the midgut transcriptomic responses of resistant and susceptible C. suppressalis strains to Cry1C toxin and to provide information for insect resistance management.

Results

A C. suppressalis midgut transcriptome of 139,206 unigenes was de novo assembled from 373 million Illumina HiSeq and Roche 454 clean reads. Comparative analysis identified 5328 significantly differentially expressed unigenes (DEGs) between C. suppressalis Cry1C-resistant and -susceptible strains. DEGs encoding Bt Cry toxin receptors, aminopeptidase-P like protein, the ABC subfamily and alkaline phosphatase were downregulated, suggesting an association with C. suppressalis Cry1C resistance. Additionally, Cry1C resistance in C. suppressalis may be related to changes in the transcription levels of enzymes involved in hydrolysis, digestive, catalytic and detoxification processes.

Conclusion

Our study identified genes potentially involved in Cry1C resistance in C. suppressalis by comparative transcriptome analysis. The assembled and annotated transcriptome data provide valuable genomic resources for further study of the molecular mechanisms of C. suppressalis resistance to Cry toxins.

Multiple Known Mechanisms and a Possible Role of an Enhanced Immune System in Bt-Resistance in a Field Population of the Bollworm, Helicoverpa zea: Differences in Gene Expression with RNAseq

Several different agricultural insect pests have developed field resistance to Bt (Bacillus thuringiensis) proteins (ex. Cry1Ac, Cry1F, etc.) expressed in crops, including corn and cotton. In the bollworm, Helicoverpa zea, resistance levels are increasing; recent reports in 2019 show up to 1000-fold levels of resistance to Cry1Ac, a major insecticidal protein in Bt-crops. A common method to analyze global differences in gene expression is RNA-seq. This technique was used to measure differences in global gene expression between a Bt-susceptible and Bt-resistant strain of the bollworm, where the differences in susceptibility to Cry1Ac insecticidal proteins were 100-fold. We found expected gene expression differences based on our current understanding of the Bt mode of action, including increased expression of proteases (trypsins and serine proteases) and reduced expression of Bt-interacting receptors (aminopeptidases and cadherins) in resistant bollworms. We also found additional expression differences for transcripts that were not previously investigated, i.e., transcripts from three immune pathways-Jak/STAT, Toll, and IMD. Immune pathway receptors (ex. PGRPs) and the IMD pathway demonstrated the highest differences in expression. Our analysis suggested that multiple mechanisms are involved in the development of Bt-resistance, including potentially unrecognized pathways.

Midgut de novo transcriptome analysis and gene expression profiling of Spodoptera exigua larvae exposed with sublethal concentrations of Cry1Ca protein

Spodoptera exigua (Hübner) is a polyphagous pest on agricultural crops, whose control is based mainly on the application of chemical insecticides. Bacillus thuringiensis (Bt) is one of the most important biological agents that have been successfully applied as a biological control, and Cry1Ca protein is considered to be active against S. exigua. Therefore, to understand the response of S. exigua to Cry1Ca protein, high-throughput sequencing was used to analyse the S. exigua larval midgut after treatment with sublethal concentrations of Cry1Ca protein. Transcriptome data showed that a total of 98,571 unigenes with an N50 value of 1135 bp and a mean length of 653 bp were obtained. Furthermore, 2962 differentially expressed genes (DEGs) were identified after Cry1Ca challenge, including 1508 up-regulated and 1454 down-regulated unigenes. Among these DEGs, detoxification (CYP, CarE, and GST) and Bt resistance (ALP, APN, and ABC transporter)-related genes were differentially expressed in the midgut of S. exigua after Cry1Ca treatment. However, most DEGs of protective enzymes were down-regulated, while most DEGs related with serine protease and REPAT were up-regulated. Furthermore, almost all DEGs related to the immune signaling pathway, antimicrobial protein, and lysozyme were up-regulated by Cry1Ca treatment. These results indicated that the detoxification enzyme, protective enzymes, Bt resistance-related genes, serine protease, REPAT, and the immune response might have been involved in the response of S. exigua to Cry1Ca protein. In summary, analysis of the transcriptomal expression of genes involved in Cry1Ca protein against S. exigua provided potential clues for elucidating the host response processes and defensive mechanisms underlying Cry1Ca toxicity.

Identification of Cry1Ah-binding proteins through pull down and gene expression analysis in Cry1Ah-resistant and susceptible strains of Ostrinia furnacalis

Bacillus thuringiensis produces insecticidal Cry toxins used in the control of multiple insect pests. Evolution of insect resistance to Bt toxins endangers the use of Cry toxins for pest control. Analysis of the Cry1Ah-binding proteins from brush border membrane vesicles (BBMV) of Ostrinia furnacalis, Asian corn borer (ACB) from the Cry1Ah-resistant (ACB-AhR) and susceptible (ACB-BtS) strains was performed by an improved pull down assay that includes coupling Cry1Ah to NHS-activated Sepharose combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS). Our data show that Cry1Ah bound to alkaline phosphatase (ALP), cadherin-like (CAD), actin, aminopeptidase-N (APN), prophenoloxidase (proPO), serine proteinase inhibitor (SPI), immulectin, and V-ATPase and to other proteins that were not previously characterized as Cry-binding proteins in ACB-BtS strain. Analysis of Cry1Ah-pulled down proteins of the BBMV from ACB-AhR revealed that Cry1Ah toxin did not bind to ALP in ACB-AhR strain, suggesting that this protein may correlate with the resistant phenotype of this strain. Additionally, we analyzed the expression of representative genes coding for Cry1Ah-binding proteins such as ALP, APN, CAD, proPO, SPI, and immulectin by qRT-PCR. ACB-AhR showed increased expression levels of proPO (7.5 fold), ALP (6.2 fold) and APN (1.4 fold) in comparison to ACB-BtS strain. In contrast, the cad gene showed slight decreased expression in ACB-AhR strain (0.7 fold) compared with ACB-BtS strain. Our data suggest that differences in the susceptibility to Cry1Ah toxin in the ACB-AhR strain may be associated with reduced ALP binding sites and with an increased immune response. This study also brings evidence of a possible binding interaction of Cry1Ah toxin to immune related proteins like proPO.

HT-SuperSAGE of the gut tissue of a Vip3Aa-resistant Heliothis virescens (Lepidoptera: Noctuidae) strain provides insights into the basis of resistance

Multitoxin Bt-crops expressing insecticidal toxins with different modes of action, for example, Cry and Vip, are expected to improve resistance management in target pests. While Cry1A resistance has been relatively well characterized in some insect species, this is not the case for Vip3A, for which no mechanism of resistance has yet been identified. Here we applied HT-SuperSAGE to analyze the transcriptome of the gut tissue of tobacco budworm Heliothis virescens (F.) laboratory-selected for Vip3Aa resistance. From a total of 1 324 252 sequence reads, 5 895 126-bp tags were obtained representing 17 751 nonsingleton unique transcripts (UniTags) from genetically similar Vip3Aa-resistant (Vip-Sel) and susceptible control (Vip-Unsel) strains. Differential expression was significant (≥2.5 fold or ≤0.4; P < 0.05) for 1989 sequences (11.2% of total UniTags), where 420 represented overexpressed (OE) and 1569 underexpressed (UE) genes in Vip-Sel. BLASTN searches mapped 419 UniTags to H. virescens sequence contigs, of which, 416 (106 OE and 310 UE) were unambiguously annotated to proteins in NCBI nonredundant protein databases. Gene Ontology distributed 345 of annotated UniTags in 14 functional categories with metabolism (including serine-type hydrolases) and translation/ribosome biogenesis being the most prevalent. A UniTag homologous to a particular member of the REsponse to PAThogen (REPAT) family was found among most overexpressed, while UniTags related to the putative Vip3Aa-binding ribosomal protein S2 (RpS2) were underexpressed. qRT-PCR of a subset of UniTags validated the HT-SuperSAGE data. This study is the first providing lepidopteran gut transcriptome associated with Vip3Aa resistance and a foundation for future attempts to elucidate the resistance mechanism.

Transcriptome and Proteome Alternation With Resistance to Bacillus thuringiensis Cry1Ah Toxin in Ostrinia furnacalis

Background: Asian corn borer (ACB), Ostrinia furnacalis can develop resistance to transgenic Bacillus thuringiensis (Bt) maize expressing Cry1Ah-toxin. However, the mechanisms that regulate the resistance of ACB to Cry1Ah-toxin are unknown.

Objective: In order to understand the molecular basis of the Cry1Ah-toxin resistance in ACB, “omics” analyses were performed to examine the difference between Cry1Ah-resistant (ACB-AhR) and susceptible (ACB-BtS) strains of ACB at both transcriptional and translational levels.

Results: A total of 7,007 differentially expressed genes (DEGs) and 182 differentially expressed proteins (DEPs) were identified between ACB-AhR and ACB-BtS and 90 genes had simultaneous transcription and translation profiles. Down-regulated genes associated with Cry1Ah resistance included aminopeptidase N, ABCC3, DIMBOA-induced cytochrome P450, alkaline phosphatase, glutathione S-transferase, cadherin-like protein, and V-ATPase. Whereas, anti-stress genes, such as heat shock protein 70 and carboxylesterase were up-regulated in ACB-AhR, displaying that a higher proportion of genes/proteins related to resistance was down-regulated compared to up-regulated. The Kyoto encyclopedia of genes and genomes (KEGG) analysis mapped 578 and 29 DEGs and DEPs, to 27 and 10 pathways, respectively (P < 0.05). Furthermore, real-time quantitative (qRT-PCR) results based on relative expression levels of randomly selected genes confirmed the “omics” response.

Conclusion: Despite the previous studies, this is the first combination of a study using RNA-Seq and iTRAQ approaches on Cry1Ah-toxin binding, which led to the identification of longer length of unigenes in ACB. The DEGs and DEPs results are valuable for further clarifying Cry1Ah-mediated resistance.

Phylogeny and evolution of the cholesterol transporter NPC1 in insects

Sterols are essential nutrients for eukaryotes. Insects are obligate sterol auxotrophs and must acquire this key nutrient from their diets. The digestive tract is the organ for absorbing nutrients as well as sterols from food. In mice, the Niemann-Pick type C1 Like 1 (NPC1L1) gene is highly expressed in the intestine and is critical for cholesterol absorption. In contrast, the molecular mechanisms for the absorption of dietary sterols in insects have not been well studied. We annotated NPC1 genes in 39 insects from 10 orders using available genomic and transcriptomic information and inferred phylogenetic relationships. Insect NPC1 genes were grouped into two sister-clades, NPC1a and NPC1b, suggesting a likely duplication in the ancestor of insects. The former exhibited weaker gut-biased expression or a complete lack of tissue-biased expression, depending on the species, while the latter was highly enriched in the gut of three lepidopteran species. This result is similar to previous findings in Drosophila melanogaster. In insects, NPC1a accumulated non-synonymous substitutions at a lower rate than NPC1b. This pattern was consistent across orders, indicating that NPC1a evolved under stronger molecular constraint than NPC1b.

Does Bt maize expressing Cry1Ac protein have adverse effects on the parasitoid Macrocentrus cingulum (Hymenoptera: Braconidae)?

The potential effects of insect-resistant, genetically engineered (GE) crops on non-target organisms, especially on predators and parasitoids, must be evaluated before their commercial cultivation. The effects of GE maize that produces Cry1Ac toxin on the parasitoid Macrocentrus cingulum were assessed by direct bioassay and indirect bioassay. In the indirect bioassay, parasitism rate, cocoon weight and the number of M. cingulum progeny produced per host were significantly reduced when M. cingulum-parasitized Cry1Ac-susceptible Ostrinia furnacalis were fed a diet containing purified Cry1Ac; however, life-table parameters of M. cingulum were not adversely affected when the same assay was performed with Cry1Ac-resistant O. furnacalis. These results indicated that the detrimental effects detected with a Cry1Ac-susceptible host were mediated by poor host quality. In a direct bioassay, no difference in life-table parameters were detected when M. cingulum adults were directly fed a 20% honey solution with or without Cry1Ac; however, survival and longevity were significantly reduced when M. cingulum adults were fed a honey solution containing potassium arsenate, which was used as a positive control. The stability and bioactivity of Cry1Ac toxin in the food sources and Cry1Ac toxin uptake by the host insect and parasitoid were confirmed by enzyme-linked immunosorbent assay and sensitive-insect bioassays. Our results demonstrate that M. cingulum is not sensitive to Cry1Ac toxin at concentrations exceeding those encountered in Bacillus thuringiensis maize fields. This study also demonstrates the power of using resistant hosts when assessing the risk of genetically modified plants on non-target organisms and will be useful for assessing other non-target impacts.

Feeding regulation by neuropeptide Y on Asian corn borer Ostrinia furnacalis

The Asian Corn Borer Ostrinia furnacalis is a major agricultural pest. In this study, a full-length neuropeptide Y (npy) gene in O. furnacalis was sequenced and cloned from cDNA library, which contains an ORF of 273 bp by encoding 90 amino acid residues. The mature OfurNPY is composed of 29 amino acids with amidation in C-terminal. The spatiotemporal expression analysis showed that npy highest expression level was in the midgut of the fifth instar larvae (the gluttony period). When the expression of npy was knocked down by feeding or injecting dsNPY, larval food consumption, body size, and body weight were significantly inhibited compared to controls. These results indicate that NPY is an important regulator in the control of feeding of O. furnacalis.

Annual post-market environmental monitoring (PMEM) report on the cultivation of genetically modified maize MON 810 in 2015 from Monsanto Europe S.A

Following a request from the European Commission, the Panel on Genetically Modified Organisms of the European Food Safety Authority (GMO Panel) assessed the annual post-market environmental monitoring (PMEM) report for the 2015 growing season of the Cry1Ab-expressing maize event MON 810 provided by Monsanto Europe S.A. The GMO Panel concludes that the insect resistance monitoring data submitted to EFSA do not indicate a decrease in susceptibility of field Iberian populations of corn borers to the Cry1Ab protein during the 2015 season. However, since the methodology for insect resistance monitoring remained unchanged compared to previous PMEM reports, the GMO Panel reiterates its previous recommendations on resistance monitoring to provide sufficient detection sensitivity. Although the farmer alert system to report complaints about product performance could complement the information obtained from the laboratory bioassays, the GMO Panel is currently not in a position to appraise its usefulness, and therefore encourages the consent holder to provide more information on this complementary resistance monitoring tool. The data on general surveillance (GS) do not indicate any unanticipated adverse effects on human and animal health or the environment arising from the cultivation of maize MON 810. The GMO Panel reiterates its previous recommendations on the analysis of farmer questionnaires, and advises the consent holder to provide more detailed information on the conducting and reporting of the literature search in future annual PMEM reports. Moreover, the GMO Panel encourages relevant parties to continue developing a methodological framework to use existing networks in the broader context of environmental monitoring. The GMO Panel concludes that the case-specific monitoring (CSM) and GS activities of maize MON 810 as carried out by the consent holder do not provide evidence that would invalidate previous GMO Panel evaluations on the safety of maize MON 810.

Patterns of Gene Expression in Western Corn Rootworm (Diabrotica virgifera virgifera) Neonates, Challenged with Cry34Ab1, Cry35Ab1 and Cry34/35Ab1, Based on Next-Generation Sequencing

With Next Generation Sequencing technologies, high-throughput RNA sequencing (RNAseq) was conducted to examine gene expression in neonates of Diabrotica virgifera virgifera (LeConte) (Western Corn Rootworm, WCR) challenged with individual proteins of the binary Bacillus thuringiensis insecticidal proteins, Cry34Ab1 and Cry35Ab1, and the combination of Cry34/Cry35Ab1, which together are active against rootworm larvae. Integrated results of three different statistical comparisons identified 114 and 1300 differentially expressed transcripts (DETs) in the Cry34Ab1 and Cry34/35Ab1 treatment, respectively, as compared to the control. No DETs were identified in the Cry35Ab1 treatment. Putative Bt binding receptors previously identified in other insect species were not identified in DETs in this study. The majority of DETs (75% with Cry34Ab1 and 68.3% with Cry34/35Ab1 treatments) had no significant hits in the NCBI nr database. In addition, 92 DETs were shared between Cry34Ab1 and Cry34/35Ab1 treatments. Further analysis revealed that the most abundant DETs in both Cry34Ab1 and Cry34/35Ab1 treatments were associated with binding and catalytic activity. Results from this study confirmed the nature of these binary toxins against WCR larvae and provide a fundamental profile of expression pattern of genes in response to challenge of the Cry34/35Ab1 toxin, which may provide insight into potential resistance mechanisms.

A potential and novel type transgenic corn plant for control of the Corn Borer

The corn borer is a world-wide agricultural pest. In this study, a full-length neuropeptide F (npf) gene in Ostrinia furnacalis was sequenced and cloned from a cDNA library, in which the npf gene produces two splicing mRNA variants - npf1 and npf2 (with a 120 bp segment inserted into the npf1 sequence to generate npf2). A spatio-temporal expression analysis showed that the highest expression level of npf was in the midgut of 5^th instar larvae (the gluttony period), and their npf expression and food consumption were significantly promoted after food deprivation for 6 h. When npf was knocked down by double-stranded RNA for NPF, larval food intake, weight and body size were effectively inhibited through changes of a biosynthesis and metabolism pathway; i.e. gene silencing of NPF causes decreases of total lipid and glycogen and increases of trehalose production. Moreover, we produced transgenic corn plants with stably expressed dsNPF. Results showed that O. furnacalis larvae fed on these transgenic leaves had lower food consumption and smaller body size compared to controls. These results indicate that NPF is important in the feeding control of O. furnacalis and valuable for production of potential transgenic corn.

Current knowledge of detoxification mechanisms of xenobiotic in honey bees

The western honey bee Apis mellifera is the most important managed pollinator species in the world. Multiple factors have been implicated as potential causes or factors contributing to colony collapse disorder, including honey bee pathogens and nutritional deficiencies as well as exposure to pesticides. Honey bees' genome is characterized by a paucity of genes associated with detoxification, which makes them vulnerable to specific pesticides, especially to combinations of pesticides in real field environments. Many studies have investigated the mechanisms involved in detoxification of xenobiotics/pesticides in honey bees, from primal enzyme assays or toxicity bioassays to characterization of transcript gene expression and protein expression in response to xenobiotics/insecticides by using a global transcriptomic or proteomic approach, and even to functional characterizations. The global transcriptomic and proteomic approach allowed us to learn that detoxification mechanisms in honey bees involve multiple genes and pathways along with changes in energy metabolism and cellular stress response. P450 genes, is highly implicated in the direct detoxification of xenobiotics/insecticides in honey bees and their expression can be regulated by honey/pollen constitutes, resulting in the tolerance of honey bees to other xenobiotics or insecticides. P450s is also a key detoxification enzyme that mediate synergism interaction between acaricides/insecticides and fungicides through inhibition P450 activity by fungicides or competition for detoxification enzymes between acaricides. With the wide use of insecticides in agriculture, understanding the detoxification mechanism of insecticides in honey bees and how honeybees fight with the xenobiotis or insecticides to survive in the changing environment will finally benefit honeybees' management.

Transcriptome Comparison Analysis of Ostrinia furnacalis in Four Developmental Stages

The Asian corn borer, Ostrinia furnacalis, is one of the most destructive pests of maize and causes huge losses in maize yield each year. In order to characterize the different developmental stages, a high-throughput sequencing platform was employed to perform de novo transcriptome assembly and gene expression analysis for the egg, larva, pupa and adult stages. Approximately 185 million reads were obtained, trimmed, and assembled into 42,638 unigenes with an average length of 801.94 bp and an N50 length of 1,152 bp. These unigene sequences were annotated and classified by performing Gene Ontology (GO), Cluster of Orthologous Groups (KOG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional classifications. Comparison of the gene expression profiles of the two transitional stages revealed dramatic differences. Some differentially expressed genes are associated with digestion, cuticularization olfactory recognition and wing formation as well as growth and development. In total, 12 putative insect development-related genes were identified. Real-time quantitative PCR (RT-qPCR) results and sequencing based on relative expression levels of randomly selected genes confirmed these expression patterns. These data represent the most comprehensive transcriptomic resource currently available for O. furnacalis and will facilitate the study of developmental pathways, cuticularization, wing formation and olfactory recognition.

Potential resistance management for the sustainable use of insect-resistant genetically modified corn and rice in China

Many lines of insect-resistant genetically modified (IRGM) corn and rice containing Bacillus thuringiensis (Bt) insecticidal genes have been developed and undergone different environmental biosafety assessments stages in China, showing robust application prospects. The potential of targeted pests to develop resistance to Bt crops is widespread, which threatens the sustainable utility of IRGM corn and rice. In this study, the potential risks of target pest complexes developing resistance to IRGM corn and rice are evaluated. Theoretical and empirical studies implementing precautionary insect resistance management (IRM) strategies to delay resistance evolution are summarized and challenges to IRM are discussed. Additionally, solutions facing these challenges are proposed. Finally, directions for future studies in developing IRGM corn and rice and IRM plans are discussed.

Identification of differentially expressed microRNAs between Bacillus thuringiensis Cry1Ab-resistant and -susceptible strains of Ostrinia furnacalis

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), can develop strong resistance to Cry1Ab, the most widely commercialized Cry toxin for Bt maize worldwide. It is essential to understand the mechanism of resistance for management of this species, but information on the post-transcriptional regulation of Bt resistance in this target insect is limited. In the present study, RNA was extracted from the ACB in various larval stages (1–5 instar) from Cry1Ab-sensitive (ACB-BtS) and -resistant (ACB-AbR) strains, each of which included two biological replicates. Using Illumina sequencing, a total of 23,809,890 high-quality reads were collected from the four ACB libraries. The numbers of known microRNAs (miRNAs) were 302 and 395 for ACB-BtS and 268 and 287 for ACB-AbR. Using Mireap software, we identified 32 and 16 potential novel miRNAs for ACB-BtS and 18 and 22 for ACB-AbR. Among them, 21 known and 1 novel miRNAs had significantly different expression between ACB-BtS and ACB-AbR. Several miRNAs were observed to target potential Bt receptor genes, such as aminopeptidase N and cadherin-like protein. The glycosylphosphatidylinositol-anchor biosynthetic process and ABC transporters pathway were identified through Gene Ontology and KEGG pathway analysis of target genes of the differentially expressed miRNAs.

Transcriptional analysis of susceptible and resistant European corn borer strains and their response to Cry1F protoxin

Background

Despite a number of recent reports of insect resistance to transgenic crops expressing insecticidal toxins from Bacillus thuringiensis (Bt), little is known about the mechanism of resistance to these toxins. The purpose of this study is to identify genes associated with the mechanism of Cry1F toxin resistance in European corn borer (Ostrinia nubilalis Hübner). For this, we compared the global transcriptomic response of laboratory selected resistant and susceptible O. nubilalis strain to Cry1F toxin. We further identified constitutive transcriptional differences between the two strains.

Results

An O. nubilalis midgut transcriptome of 36,125 transcripts was assembled de novo from 106 million Illumina HiSeq and Roche 454 reads and used as a reference for estimation of differential gene expression analysis. Evaluation of gene expression profiles of midgut tissues from the Cry1F susceptible and resistant strains after toxin exposure identified a suite of genes that responded to the toxin in the susceptible strain (n = 1,654), but almost 20-fold fewer in the resistant strain (n = 84). A total of 5,455 midgut transcripts showed significant constitutive expression differences between Cry1F susceptible and resistant strains. Transcripts coding for previously identified Cry toxin receptors, cadherin and alkaline phosphatase and proteases were also differentially expressed in the midgut of the susceptible and resistant strains.

Conclusions

Our current study provides a valuable resource for further molecular characterization of Bt resistance and insect response to Cry1F toxin in O. nubilalis and other pest species.

Electronic supplementary material

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

Companion cropping with potato onion enhances the disease resistance of tomato against Verticillium dahliae

Intercropping could alleviate soil-borne diseases, however, few studies focused on the immunity of the host plant induced by the interspecific interactions. To test whether or not intercropping could enhance the disease resistance of host plant, we investigated the effect of companion cropping with potato onion on tomato Verticillium wilt caused by Verticillium dahliae (V. dahliae). To investigate the mechanisms, the root exudates were collected from tomato and potato onion which were grown together or separately, and were used to examine the antifungal activities against V. dahliae in vitro, respectively. Furthermore, RNA-seq was used to examine the expression pattern of genes related to disease resistance in tomato companied with potato onion compared to that in tomato grown alone, under the condition of infection with V. dahliae. The results showed that companion cropping with potato onion could alleviate the incidence and severity of tomato Verticillium wilt. The further studies revealed that the root exudates from tomato companied with potato onion significantly inhibited the mycelia growth and spore germination of V. dahliae. However, there were no significant effects on these two measurements for the root exudates from potato onion grown alone or from potato onion grown with tomato. RNA-seq data analysis showed the disease defense genes associated with pathogenesis-related proteins, biosynthesis of lignin, hormone metabolism and signal transduction were expressed much higher in the tomato companied with potato onion than those in the tomato grown alone, which indicated that these defense genes play important roles in tomato against V. dahliae infection, and meant that the disease resistance of tomato against V. dahliae was enhanced in the companion copping with potato onion. We proposed that companion cropping with potato onion could enhance the disease resistance of tomato against V. dahliae by regulating the expression of genes related to disease resistance response. This may be a potential mechanism for the management of soil-borne plant diseases in the intercropping system.