http://orcid.org/0000-0002-3536-7764
Figures
Abstract
Cnaphalocrocis medinalis (Guenée) is one of the important insect pests in rice field. Bt agents were recommended in the C. medinalis control and Bt rice is bred as a tactic to control this insect. However, the tolerance or resistance of insect to Bt protein is a main threat to the application of Bt protein. In order to investigate the response of C. medinalis transcriptome in defending a Cry1C toxin, high-through RNA-sequencing was carried in the C. medinalis larvae treated with and without Cry1C toxin. A total of 35,586 high-quality unigenes was annotated in the transcriptome of C. medinalis midgut. The comparative analysis identified 6,966 differently expressed unigenes (DEGs) between the two treatments. GO analysis showed that these genes involved in proteolysis and extracellular region. Among these DEGs, carboxylesterase, glutathione S-transferase and P450 were differently expressed in the treated C. medinalis midgut. Furthermore, trypsin, chymotrypsin, and carboxypeptidase were identified in DEGs, and most of them up-regulated. In addition, thirteen ABC transporters were downregulated and three upregulated in Cry1C-treated C. medinalis midgut. Based on the pathway analysis, antigen processing and presentation pathway, and chronic myeloid leukemia pathway were significant in C. medinalis treated with Cry1C toxin. These results indicated that serine protease, detoxification enzymes and ABC transporter, antigen processing and presentation pathway, and chronic myeloid leukemia pathway may involved in the response of C. medinalis to Cry1C toxin. This study provides a transcriptomal foundation for the identification and functional characterization of genes involved in the toxicity of Bt Cry protein against C. medinalis, and provides potential clues to the studies on the tolerance or resistance of an agriculturally important insect pest C. medinalis to Cry1C toxin.
Citation: Yang Y, Xu H, Lu Y, Wang C, Lu Z (2018) Midgut transcriptomal response of the rice leaffolder, Cnaphalocrocis medinalis (Guenée) to Cry1C toxin. PLoS ONE 13(1): e0191686. https://doi.org/10.1371/journal.pone.0191686
Editor: Yulin Gao, Chinese Academy of Agricultural Sciences Institute of Plant Protection, CHINA
Received: September 6, 2017; Accepted: January 9, 2018; Published: January 23, 2018
Copyright: © 2018 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was funded by National Natural Science Foundation of China (31501669), the earmarked fund for China Agriculture Research System (CARS-01-36), National Special Key Project of China (2016ZX08001-001), and Zhejiang Provincial Natural Science Foundation of China (LY15C140003).
Competing interests: The authors have declared that no competing interests exist.
Introduction
The rice leaffolder, Cnaphalocrocis medinalis (Guenée), is one of the important migratory insect pests of rice in China [1–2]. It widely distributes in humid tropical and temperate regions of Oceania, Africa, and Asia [3–4]. It scrapes folds and feeds the rice leaves causing the reduction of photosynthetic activity [1,5]. Heavy occurrence often happens in main rice-growing areas like China, Japan, Korea, Vietnam, Philippines, and so on [3], and frequently results in great rice yield losses [2]. The damage area by C. medinalis were more than 20 million hm2 in China in the most of years ranged from 2003 to 2010 [2]. A report from the Ministry of Agriculture, China showed that the damage area by C. medinalis was up to 15.5 million hm2 in 2015 [4].
So far, the control of C. medinalis mainly relied on the chemicals. However, the misuse and overuse of chemicals resulted in many negative issue such as insect resistance, resurgence, environmental pollution, human health concerned by public [6]. Environment-friendly methods are recommended in the pest control. Bacillus thuringiensis, as one of the important biological agents has been applied in the pest control for more than 60 years due to its crystal protein with insecticidal activity against insect pests [7–9]. Cry genes derived from B. thuringiensis were successfully transferred into plants like cotton, rice and these plants harbored resistance to Lepidopteran insect pests [10–13]. In China, Bt agent was recommended as a friendly bio-agent to use in the control of C. medinalis and Bt rice in China was approved to release for production test in limited area in 2009, but not commercially planted[11]. The application of Bt protein and Bt crop could greatly reduce the use of chemical insecticides. However, insect resistance is a big threat to the application of Bt protein and Bt crop. So far, at least eight insect pests evolved resistance against Bt crop or Bt protein [14–15]. Coincidentally, many insect pests evolved resistance to Bt proteins under selection with Bt protein in the laboratory condition [8,16]. Considering the sustainable use of the Bt protein or Bt rice, the potential of resistance evolution of C. medinalis against Bt protein will become an inescapable issue [17].
Midgut, one of the important organs in the insect body, not only is a place for food ingestion and utilization, but also a place for detoxification to xenobiotics [18]. Interaction of Bt protein and insect mainly occurred in midgut. Several reports revealed that resistance of insect against Bt protein were correlated with the variation of genes in midgut [14,19–21]. Next generation sequencing (NGS) is a strong technology which could generate large volumes of sequence information [22]. NGS provides us with unprecedented high-throughput and low-cost sequencing platforms applied in a variety of manners, including de novo whole-genome sequencing, resequencing of genomes to identify variations, de novo transcriptome and gene expression profiling, and detecting methylation patterns [23]. Now, the Illumina/Solexa sequencing technology is dominated in the NGS market, featuring high data accuracy and a broad range of applications.
In this study, transcriptome of C. medinalis midgut was sequenced and de novo assembled using Illumina sequencing technology. The gene expression patterns in the midgut were compared in C. medinalis larvae treated with and without Bt Cry1C protein. Moreover, the genes of digestion and detoxification system were given greater attention. The results provide clues on the roles of midgut in the response to Bt proteins, information for further study candidate genes involved in interaction of Bt protein and C. medinalis and promote the understanding of potential resistance or tolerance of C. medinalis to Bt proteins.
Materials and methods
Insect rearing
Adults of C. medinalis were collected with sweep net from paddy fields in suburb of Hangzhou, Zhejiang, China in August, 2015. The moths were placed in plastic cup covered with nylon mesh and fed with 10% honey solution. Eggs laid on the mesh were removed and transferred to a box with detached leaf of 45 d-old Taichung Native 1 (TN1). The larvae that hatched from eggs were cultured in the detached leaf then used in the experiment. All insect cultures were kept at 27±1°C with 70–80% RH and a photoperiod of 14: 10 (L: D) h. Insect collection is from our experimental field, and C. medinalis is not endangered or protected species.
Treatment of C. medinalis larvae with Cry1C toxin
Fifth instar larvae of C. medinalis were selected for the feeding experiments. The larvae were reared on rice leaves dipped with Cry1C activated toxin (MP, Cavey, CWRU, US) solution which resulted in a growth inhibition but did not cause visible death of the larvae for 48 hours. Larvae treated with PBS buffer were used as negative control. All treated larvae were incubated for 48 h at 27±1°C with 70–80% RH and a photoperiod of 14: 10 (L: D) h. Approximately 150 larvae were used for each treatment.
RNA isolation, cDNA library construction, transcriptomic sequencing
Total RNA was isolated from the midgut of C. medinalis larvae using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's protocol. The total RNA quantity and purity were analyzed with Bioanalyzer 2100 and RNA 6000 Nano Lab Chip Kit (Agi-lent, CA, USA) with RIN number >7.0.
Approximately 10 μg of total RNA was subjected to isolate Poly (A) mRNA with poly-T oligo attached magnetic beads (Invitrogen). Following purification, the mRNA is fragmented into small pieces using divalent cations under elevated temperature, and the cleaved RNA fragments were reverse-transcribed to create the final cDNA library in accordance with the protocol for the mRNA-Seq sample preparation kit (Illumina, San Diego, USA). The cDNA library was sequenced run with Illumina HiSeqTM2000 sequence platform at LC Sceiences (USA) following the vendor's recommended protocols. Two or three biological replicates for each treatment were prepared.
De novo RNA-seq assembly and annotation
After sequencing, the raw data were filtered against low-quantity reads and adaptor contamination using Trinity. Clean reads were de novo assembled using Trinity with K-mer = 25 [24]. Annotations of all the unigenes were performed by a BLASTx search against the Nr, Swiss-Prot, KOG, KEGG databases and Pfam with a cut-off E-value of ≤1e-5. We obtained the Gene Ontology terms (http://www.geneontology.org) of each C. medinalis unigene with the software Blast2GO (http://www.blast2go.org) using the default parameters [25]. Gene expression level was normalized by RSEM-based lgorithm to get RPKM value [26–27]. Based on the expression levels, thresholds of FDR<0.05 and log2 fold-change (log2FC) ≥1 were set for identifying significant differential expressed unigenes (DEGs) between control and Cry1C treated larvae.
Enrichment analysis
All the DEGs were mapped to terms in KEGG databases to identify significantly enriched metabolic pathway or signal transduction pathways in DEGs. Pathway enrichment analysis provides all terms that significantly enriched in differentially expressed genes in comparison to the control C. medinalis. The P-value was calculated for each pathway by Fisher's exact test with a Q value threshold<0.05.
Quantitative real-time PCR (qRT-PCR) validation
To conform the data, a subset of DEGs was validated by qRT-PCR. qRT-PCR was performed in 20 μL reaction mixtures composed of 2.5 μL of template cDNA,10 μL of 2×SYBR Green PCR Master Mix (Fermentas,Waltham, MA, USA), and 0.25 mM each primer on the CFX96 Real-Time System (Bio-Rad, Hercules, CA, USA). The selected genes were verified with the following cycling conditions: 95°C for 30 s, followed by 40 cycles of 95°C for 30 s, 60°C for 35 s. The melting curve analysis was used to analyze the specificity of the qPCR product. The sequences of the primers used are listed in S1 Table. A ubiquinol-cytochrome c reductase served as an internal control. The relative gene expression values were calculated using the 2−ΔΔCt method [28].
Results
Assembly and annotation of the C. medinalis transcriptome
The transcriptomes of the midgut of C. medinalis larvae treated with and without Cry1C protein were sequenced and compared for each treatment. A total of 37.33G bases of clean reads were obtained from the cDNA libraries, and 35,586 high-quality unigenes and 65,016 transcripts were assembled (Table 1). The average gene and transcript length were 784 and 895 base pair (bp), respectively (Table 1). Mean GC% of gene and transcript were 40.66% and 39.18%, respectively. N50 (the shortest sequence length at 50% of the transcriptome) of gene and transcript were 1,211 and 1,378 bp, respectively. The size distributions of all the unigenes studied were as follows: 18,720 unigenes (17.8%) exhibited lengths of more than 2,000 bp; 17,474 unigenes (16.7%) showed lengths between 500 and 1000 bp; 26,983 unigenes (25.7%) demonstrated lengths between 300 and 500 bp; and 41,789 unigenes (39.8%) were shorter than 300 bp (S1 Fig). The size distributions of all the transcripts studied were as follows: 18,720 unigenes (17.8%) exhibited lengths of more than 2,000 bp; 17,474 unigenes (16.7%) showed lengths between 500 and 1,000 bp; 26,983 unigenes (25.7%) demonstrated lengths between 300 and 500 bp; and 41,789 unigenes (39.8%) were shorter than 300 bp (S2 Fig).
To annotate these unigenes, All unigene sequences were against the protein databases (Nr, SwissProt, KEGG, COG, Pfam, GO) using BLASTx (E-value<e-5). Only 16,741 unigenes (47.04% of all unigenes) displayed Nr annotations. 9,789 unigenes (27.51% of all unigenes) displayed swiss-prot annotations. 12,317 unigenes (34.61% of all unigenes) displayed Pfam annotations. 7,401unigenes (20.80% of all unigenes) displayed KEGG annotations. 11,010 unigenes (30.94% of all unigenes) displayed KOG annotations. 9,003 unigenes (25.30% of all unigenes) displayed KEGG annotations. 11,010 unigenes (30.94% of all unigenes) displayed GO annotations. Only small partial genes were annotated to the databases because of the short nucleotide lengths and the lack of C. medinalis genomic information.
Species-specific distribution revealed that 37.2% of the C. medinalis sequences matched with the silkworm, Bombyx mori. In addition, 27.2% of the C. medinalis sequences resembled those of Danaus plexippus, followed by those of Papilio xuthus (2.8%), Tribolium castaneum (1.3%), Zootermopsis nevadensis (1.2%), Papilio polytes (0.9%) (Fig 1A).
(A) Species distribution of the BLASTX results. (B) GO categories of all unigenes and DEGs. (C) euKaryotic orthologous Groups (KOG) classification.
Classifications for C. medinalis unigenes
GO assignments (level 2) were used to predict the functions of the C. medinalis unigenes. A total of 12,653 unigenes were assigned GO terms of three categories (biological process, cellular component, and molecular function) based on sequence homology. The dominant items of the three categories were nucleus, ATP-binding, cytoplasm, integral to membrane, transcription DNA-dependent. (Fig 1B).
The C. medinalis unigenes were also annotated with COG classifications. Among the 11,010 unigenes with COG annotations, 21% (2,063) belonged to the cluster of general function; followed by 10% (1,079) under the cluster of Signal transduction mechanisms; and 10% (986) under the cluster of Posttranslational modification, protein turnover, chaperones. By contrast, Cytoskeleton corresponded to the smallest clusters (Fig 1C).
The C. medinalis unigenes were also annotated with KEGG classifications. Among the 9,003 unigenes with KEGG annotations, 21% (2,145) belonged to the cluster of Metabolism; followed by 10% (1,014) under the cluster of Organismal Systems; and 10% (991) under the cluster of Human Diseases. Genetic Information Processing, Cellular Processes, Environmental Information Processing.
Overview on complex response of C. medinalis to Cry1C toxin
The number of DEGs in the midgut of C. medinalis treated with Cry1C was 6,966. Of that, 4,621 genes were down-regulated, and 2,345 genes were up-regulated (Fig 2A). The changed ratios of most DEGs accumulated between 2−15 and 213 (Fig 2B). For the top 20 upregulated unigenes in the midgut of C. medinalis treated with Cry1C, three unigenes had explicit annotations, i.e., adenosinetriphosphatase, nucleoprotein TPR, and heat shock 70kDa protein 1/8 (S2 Table). Meanwhile, for the top 20 downregulated unigenes in the midgut of C. medinalis treated with Cry1C, 10 unigenes had explicit annotations, i.e., encoded NADH dehydrogenase I subunit 4, beta-2-microglobulin, glypican 5, CD74 antigen, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, valyl-tRNA synthetase, tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, valyl-tRNA synthetase, protein transport protein SEC61 subunit alpha, saposin (S3 Table).
(A) Number of DEGs. (B) Volcano plot to show the fold change and error rates. The non-DEGs are indicated by black dots, the DEGs up-regulated are indicated by red dots, and the DEGs down-regulated are indicated by blue dots.
From the statistic of GO enrichment, most of DEGs involved in transcription DNA-dependent, regulation of transcription, and proteolysis under biological process; most of DEGs involved in nucleusm, cytoplasm, and integral to membrane under cellular component function item; most of DEGs involved in ATP-binding, metal ion binding, and DNA binding under molecular function (Fig 3A). According to the GO enrichment scatterplot (Fig 3B), large amount of genes involved in proteolysis, extracellular region, extracellular space, and serine-type endopeptidase activity were differently expressed in larval midgut of C. medinalis treat with and without Cry1C toxin.
(A) Histogram. (B) Scatter plot.
To validate the transcriptome data, qRT-PCR analysis was performed for DGEs (S1 Table). The qRT-PCR validation results were presented as fold-changes normalized to the ubiquinol-cytochrome c reductase gene. A comparison of the mean values of qRT-PCR results and the DEGs is summarized in S3 Fig. The consistent expression trend in qRT-PCR results and the original DEGs indicates the reliability of gene expression profiling from transcriptome analysis.
Interest genes and Bt related genes in the transcriptome of C. medinalis midgut
Here, we investigated DEGs potentially involved in insecticidal tolerance/resistance in C. medinalis, including metabolic and Bt-related genes. DEGs sequences encoding enzymes function ingestion and detoxification and receptors of the insecticide and toxins were extracted from the RNA-seq database from the C. medinalis midgut. Here, alkaline phosphatase, cytochrome P450, glutathione S-transferase (GST), carboxylesterase, ABC transporter, and serine protease including transmembrane protease, trypsin and chymotrypsin were (FDR < 0.05 and the absolute value of log2FC ≥ 1) were found significant different in the midgut of C. medinalis larvae treated with Cry1C protein (Table 2). A total of thirteen DEGs encoding carboxylesterases were differentially expressed between treated and untreated C. medinalis, of which 11 of these unigenes were up-regulated and two were down-regulated in the Cry1C-treated C. medinalis (S4 Table). A total of 99 serine protease genes were differently expressed including trypsin, chymotrypsin, carboxypeptidase, aminopeptidase, elastase, cathepsin and so on (Fig 4). Fourty-one unigenes encoding trypsin were differentially expressed between treated and untreated C. medinalis, of which thirty-four unigenes were up-regulated and seven were down-regulated in the Cry1C-treated C. medinalis (Table 3). Eight chymotrypsin were 1.24- to 3.69- fold upregulated. Twenty carboxypeptidase genes were differently expressed with 17 genes upregulated and 3 downregulated in the Cry1C-treated C. medinalis (Table 4). Four aminopeptidase N were upregulated and two downregulated, and two alkaline phosphatase were upregulated in the midgut of C. medinalis larvae treated with Cry1C. One cadherin in Cry1C treated C. medinalis larvae were differently downregulated by 1.36-fold compared with those of the untreated larvae. The DEGs were annotated to the ABC transporter including five subfamilies (ABCB, ABCC, ABCE, ABCF and ABCG), and three upregulated and 13 down regulated (Fig 5).
Significant enriched pathways in the C. medinalis larvae treated with Cry1C toxin
Metabolic pathway enrichment analysis demonstrated that 770 DEGs involved in 232 pathways (S5 Table) might play a role in the response of C. medinalis to Cry1C toxin. DEGs were significantly enriched in two functional categories (antigen processing and presentation, and chronic myeloid leukemia) in C. medinalis treated with Cry1C (Q-value < 0.05) (Table 5). Twenty-five DGEs were found in antigen processing and presentation pathway, and twelve DGEs were found in chronic myeloid leukemia pathway.
Discussion
Despite C. medinalis is one of the important insect pests in rice field of Asia, and also the target insect of Bt agent and Bt rice [3,11]. But the interaction between this insect pest and Bt toxin is unclear. The response of C. medinalis to Bt toxin related with the tolerance and resistance to Bt toxin. In this study, larval midgut transcriptome was sequenced and analyzed, 35,586 high-quality unigenes were obtained, and 6,966 DEGs were found between C. medinalis larvae fed with and without Cry1C toxin. In Plutella xylostella midgut transcriptome, 2,925 and 2,967 unigenes were differently expressed in two resistant strains compared with susceptible strains [29]. In a pink bollworm (Pectinophora gossipiella) larval midgut transcriptome, 39,874 unigenes were obtained from18,623,508 assembled reads [30]. The number of bases in C. medinalis midugt transcriptome in this study was smaller than those in P. xylostella and P. gossipiella larval midgut transcriptome. In the transcriptome of the Lymantria dispar (gypsy moth) larval midgut in response to infection by B. thuringiensis, the DEGs primarily associated with digestive function (eg. a-amylase, lipase and carboxypeptidase), immune response (eg. C-type lectin 4), developmental genes (eg. arylphorin) as well as a variety of binding proteins including cellular retinoic acid binding protein (lipid-binding), insulin-related peptide binding protein (protein binding) and ovary C/EBPg transcription factor (nucleic acid-binding) [31]. Many DGEs involved in the metabolic pathway, ABC transporter pathway, aminopeptidases and cadherins were found between Cry1Ac-resistant and susceptible strains of P. xylostella [29]. Here we focused on the interest genes involved insecticides and Bt resistance or tolerance. These findings may facilitate the understanding of interaction of insect and Bt protein.
Detoxification enzymes
Carboxylesterase, cytochrome P450, and glutathione S-transferase (GST) are important enzymes in insect system and function on the metabolic and detoxification of the xenobiotics. GST in C. medinalis could involved in the detoxification of chlorpyrifos [32]. Chen et al. [33] discovered that two P450 genes CYP6CV1 and CYP9A38 may be involved in detoxification of rice phytochemicals. Veegala and Vemuri [34] reported carboxylesterase and GST mediated in the resistance of C. medinalis to insecticide. In this study, 6 cytochrome P450, 3 GSTs, 13 carboxylesterase were significant different expression in C. medinalis treated with Cry1C, suggesting that detoxification enzymes may be involved the early response of C. medinalis larvae to Cry1C toxin. Similarly, Guo et al. [35] discovered that carboxylesterase were increased in S. exigua as the increased time of exposure to Bt cotton and considered that the allocation of carboxylesterase in body is a metabolic tolerance of herbivorous insects response to a toxic protein.
Protect enzymes
Catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD) are the important protection enzymes in insect [36]. Here, we found 1 CAT, 11 PODs and 4 SODs were found. In midgut of C. medinalis treated with Cry1C, 1 CAT, 9 PODs and 3 SODs were down regulated, and 2 other PODs and one SOD were upregulated. These results showed Cry1C impose environmental stress to C. medinalis larvae, anyhow the upregulated PODs indicated that the body of C. medinalis adjust the physiological activity to reply to Cry1C. Many reports indicated that the increase of POD is the emergence measures of the insect reply to Bt toxin [35,37].
Serine protease
Serine proteases always function on the digestion and utilization of food, and the detoxification of exogenous proteins. Trypsin and chymotrypsin are two kinds of serine proteases and involved in the activation and degradation of Bt toxin [38–40]. The changes of serine proteases may be correlated with the resistance of the insect to Bt toxins. Tetreau et al. [41] reported an increase in larval gut proteolytic activities in a Bt resistance strain of the Dengue fever mosquito. Elleuch et al. [42] found that early degradation of Cry toxins by proteases in mosquito larvae midgut may be one of the reasons why Bt protein toxicity decreased. In our study, 34 trypsin and 8 chymotrypsin were upregulated and 7 other trypsin down-regulated. In addition, twenty carboxypeptidases were differently expressed with 17 genes up-regulated and 3 genes down-regulated. These results indicated that the trypsin, chymotrypsin and carboxypeptidases could play an important role on the tolerance of C. medinalis to Cry1C through degradation.
Bt receptors
In the Bt action, Bt receptors are important on the formation of toxicity [19]. So far, at least cadherin, alkaline phosphatase, and aminopeptidase N (APN) were putative to be the receptor of Bt toxins in insects [15,19]. However, there are no reports on these proteins as Bt receptors in C. medinalis. Unigenes on these proteins could provide invaluable information for the further study to prove these proteins as Bt receptors in C. medinalis. A previous study in the cabbage looper (Trichoplusia ni) indicated that the Cry1Ac resistance was correlated with the concurrent up-regulation of APN6 and down-regulation of APN1 [43]. In our study, four APN were up-regulated (1.68- to 2.54-fold) and two APNs downregulated (2.76- and 5.57- fold) in the C. medinalis treated with Cry1C toxin. The detailed role of these APNs on the tolerance or resistance of C. medinalis to Cry1C remained to be determined. Alkaline phosphatases have also been identified as Bt receptors in many insects [19,44–46]. Chen et al. [47] found a toxin-binding alkaline phosphatase fragment could synergize Bt toxin Cry1Ac against susceptible and resistant Helicoverpa armigera. Here, two alkaline phosphatases were upregulated in C. medinalis larvae treated with Cry1C. According to the role of alkaline phosphatase in Bt action, alkaline phosphatases in tolerant or resistant insects to Bt toxin may do not increase, the results in our studies may indicate these two alkaline phosphatases may be not the receptor of the Cry1C or the results may be unexpected. Flores-Escobar et al. [48] found alkaline phosphatase is more important than APN for Cry1Ab toxicity, while Cry1Ac relied principally on APN1 through RNAi experiment. Cadherin is an important Bt receptor and its mutation was correlated with the insect resistance to Bt toxins [14,19,49]. Here, one unigene annotated as cadherin was downregulated in C. medinalis larvae treated with Cry1C. However, the upregulated cadherin belongs to flamingo-cadherin which is different with the cadherin in other insects related with Bt action. Additionally, the other cadherin in C. medinalis was not differently expressed. Whether this flamingo-cadherin correlated with Cry1C action remains to be determined.
ABC transporters
ATP-binding cassette (ABC) transporters is a super family involved in various functions including the transport of lipids, inorganic ions and especially the detoxification of xenobiotics [50]. ABC transporter plays important role in drug resistance [51]. And in insect ABC transporters were shown functions in glucoside sequestration [52], the transport of eye color pigments [53] and insecticidal resistance [54]. In this study, sixteen ABC transporters attributed to 5 subfamilies were regulated in C. medinalis treated with Cry1C and 3 upregulated and 13 down regulated. Previous studies have linked ABC transporters with Cry1 resistance in seven lepidopterans including Heliothis virescens, H. armigera, P. xylostella, T. ni, Ostrinia nubilalis, H. punctigera and B. mori [29,55–59]. In the transcriptome of P. xylostella Cry1Ac-susceptible and resistant strains, many DGEs of ABC transporters were observed [29]. The role of the ABC transporters found in this study on tolerance or resistance of C. medinalis to Cry1C warrants further investigation.
Enriched pathway
Pathway analysis indicated that two pathways were significant enriched in the C. medinalis treated with Cry1C. The antigen processing and presentation pathway had 25 DEGs in the C. medinalis treated with Cry1C. Antigen processing is an immunological process that prepares antigens for presentation to special cells of the immune system called T lymphocytes [60]. Subsequent presentation of the antigens on class I or class II major histocompatibility complex (MHC) molecules is dependent on the pathway, and MHC I and MHC II antigen presentation typically involves the endogenous and exogenous pathway of antigen processing, respectively [60]. The transporter associated with antigen processing (TAP) represents a key machinery in MHC class I antigen presentation by translocating proteasomal degradation products into the lumen of the endoplasmic reticulum (ER) for loading onto MHC class I molecules, and TAP is a member of the ABC transporters [61]. The ABC transporter pathway was an important enriched one due to previous report that ABC transporter was correlated with Bt resistance in insects [62–63]. Chronic myeloid leukemia (CML) pathway correlated with the apoptosis, and twelve DGEs were found in this pathway. ABC transporter also was important in the CML pathway [64]. In addition, a mechanism of Bt action was reported to be related to apoptosis [65–67]. These two pathways may provide more information on the Bt action and insect response to Bt protein.
In conclusion, multiple aspects were involved in the C. medinalis response to Cry1C protein from the transcriptomal data. Serine protease, detoxification enzymes, ABC transporter, antigen processing and presentation pathway, and chronic myeloid leukemia pathway may be main elements involved in the response of C. medinalis to Cry1C toxin.
Supporting information
S1 Fig. Distribution of unigenes from the midgut transcriptome of Cnaphalocrocis medinalis.
https://doi.org/10.1371/journal.pone.0191686.s001
(TIF)
S2 Fig. Distribution of transcripts from the midgut transcriptome of Cnaphalocrocis medinalis.
https://doi.org/10.1371/journal.pone.0191686.s002
(TIF)
S3 Fig. Validation of transcriptome data by qRT-PCR.
The means of at least three biological replicates are presented as log2FC ± SE.
https://doi.org/10.1371/journal.pone.0191686.s003
(TIF)
S1 Table. List of primers of genes for the validation of transcriptome data by qRT-PCR.
https://doi.org/10.1371/journal.pone.0191686.s004
(DOCX)
S2 Table. Top 20 downregulated unigenes in the midgut of Cnaphalocrocis medinalis treated with Cry1C toxin.
https://doi.org/10.1371/journal.pone.0191686.s005
(DOCX)
S3 Table. Top 20 upregulated unigenes in the midgut of Cnaphalocrocis medinalis treated with Cry1C toxin.
https://doi.org/10.1371/journal.pone.0191686.s006
(DOCX)
S4 Table. Carboxylesterase genes differently expressed in Cnaphalocrocis medinalis larvae treated with Cry1C toxin.
https://doi.org/10.1371/journal.pone.0191686.s007
(DOCX)
S5 Table. KEGG summary of Cnaphalocrocis medinalis response to Cry1C toxin.
https://doi.org/10.1371/journal.pone.0191686.s008
(XLSX)
References
- 1.
Cheng JA. Rice Pests. Beijing: China Agricultural Press; 1996.
- 2. Yang YJ, Xu HX, Zheng XS, Tian JC, Lu YH, Lu ZX. Progresses in management technology of rice leaffolders in China. J Plant Prot. 2015; 42(5): 691–701.
- 3. Barrion AT, Litsinger JA, Medina EB, Aguda RM, Bandong JP, Pantua PC Jr, et al. The rice Cnaphalocrocis and Marasmia (lepidoptera: Pyralidae) leaffolder complex in the Philippines: taxonomy, bionomics, and control. Phil Entomol. 1991; 8(4): 987–1074.
- 4.
Lu ZX. Pictorial handbook on sustainable management of rice leaffolder. Beijing: China Agricultral Press; 2017.
- 5. de Jong PD. Effect of folding and feeding by Cnaphalocrocis medinalis on photosynthesis and transpiration of rice leaves. Entomol Exp Appl. 1992; 63(1): 101–102.
- 6. Heckel DV. Insecticide resistance after Silent Spring. Science 2012; 337(6102): 1612–1614. pmid:23019637
- 7. Ali S, Zafar Y, Ali G.M, Nazir F. Bacillus thuringiensis and its application in agriculture. Afr J Biotechnol. 2010; 9(14): 2022–2031.
- 8. Tabashnik BE, Rensburg JBJV, Carrière Y. Field-evolved insect resistance to Bt crops: definition, theory, and data. J Econ Entomol. 2009; 102(6):2011–2025. pmid:20069826
- 9. Yang YJ, Xu HX, Zheng XS, Lu ZX. Susceptibility and selectivity of Cnaphalocrocis medinalis (Lepidoptera: Pyralidae) to different Cry toxins. J Econ Entomol. 2012; 105(6): 2122–2128. pmid:23356078
- 10.
Cohen MB, Chen M, Bentur JS, Heong KL, Ye G. Bt Rice in Asia: Potential Benefits, Impact, and Sustainability, In: Romeis J, Shelton AM, Kennedy GG, editors. Integration of Insect-Resistant Genetically Modified Crops within IPM Programs. New York: Springer Science + Business Media B.V.; 2008. pp. 223–248.
- 11. Chen M, Shelton A, Ye GY. Insect-resistant genetically modified rice in China: from research to commercialization. Annu Rev Entomol. 2011; 56(1): 81–101.
- 12. Yang Z, Chen H, Tang W, Hua H, Lin Y. Development and characterisation of transgenic rice expressing two Bacillus thuringiensis genes. Pest Manag Sci. 2011; 67(4): 414–422. pmid:21394874
- 13.
James C. Global Status of Commercialized Biotech/GM Crops: 2016. ISAAA Brief No. 52.:, Ithaca: ISAAA; 2016.
- 14. Yang YJ, Chen HY, Wu YD, Yang YH, Wu SW. Mutated cadherin alleles from a field population of Helicoverpa armigera confer resistance to Bacillus thuringiensis toxin Cry1Ac. Appl Environ Microbiol. 2007; 73(21): 6939–6944. pmid:17827322
- 15. Wu YD. Detection and mechanisms of resistance evolved in insects to Cry toxins from Bacillus thuringiensis. Adv Insect Physiol. 2014; 47: 297–342.
- 16. van Frankenhuyzen K. Insecticidal activity of Bacillus thuringiensis crystal proteins. J Insect Physiol. 2009; 101(1): 1–16.
- 17. Wu ZH, Yang YJ, Xu HX, Zheng XS, Tian JC, Lu YH, et al. Changes in growth and development and main enzyme activities in midgut of Cnaphalocrocis medinalis intermittently treated with low amount of Bt rice leaves over generations. Chin J Rice Sci. 2015; 29(4): 417–423.
- 18.
Linser PJ, Dinglasan RR. Chapter One—Insect gut structure, function, development and target of biological toxins, In: Dhadialla TS, Gill SS, editors. Advances in Insect Physiology. New York: Academic Press; 2014; 47: 1–13.
- 19. Pigott CR, Ellar DJ. Role of receptors in Bacillus thuringiensis crystal toxin activity. Microbiol Mol Biol Rev. 2007; 71(2): 255–281. pmid:17554045
- 20. Yamazaki T, Ishikawa T, Pandian GN, Okazaki K, Haginoya K, Tachikawa Y, et al. Midgut juice of Plutella xylostella highly resistant to Bacillus thuringiensis Cry1Ac contains a three times larger amount of glucosinolate sulfatase which binds to Cry1Ac compared to that of susceptible strain. Pestic Biochem Phys. 2011; 101(2): 125–131.
- 21. Brunet JF, Vachon V, Marsolais M, Van Rie J, Schwartz JL, Laprade R. Midgut juice components affect pore formation by the Bacillus thuringiensis insecticidal toxin Cry9Ca. J. Invertebr. Pathol. 2010; 104(3): 203–208. pmid:20399787
- 22. Metzker ML. Sequencing technologies–the next generation. Nat Rev Genet. 2010; 11(1): 31–46 pmid:19997069
- 23. Zhou X, Ren L, Meng Q, Li Y, Yu Y, Yu J. The next-generation sequencing technology and application. Protein Cell 2010; 1(6): 520–536. pmid:21204006
- 24. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29(7): 644–652. pmid:21572440
- 25. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005; 21(18): 3674–3676. pmid:16081474
- 26. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 2008; 5(7): 621–628. pmid:18516045
- 27. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 2011; 12: 323. pmid:21816040
- 28. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001; 25, 402–408. pmid:11846609
- 29. Lei Y, Zhu X, Xie W, Wu Q, Wang S, Guo Z, et al. Midgut transcriptome response to a Cry toxin in the diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae). Gene 2014; 533(1): 180–187. pmid:24120626
- 30. Tassone EE, Zastrowhayes G, Mathis J, Nelson ME, Wu G, Flexner JL, et al. Sequencing, de novo assembly and annotation of a pink bollworm larval midgut transcriptome. GigaScience 2016; 5: 28. pmid:27333791
- 31. Sparks ME, Blackburn MB, Kuhar D, Gundersenrindal DE. Transcriptome of the Lymantria dispar (gypsy moth) larval midgut in response to infection by Bacillus thuringiensis. Plos One 2013; 8(5): e61190. pmid:23658687
- 32. Liu S, Rao XJ, Li MY, Feng MF, He MZ, Li SG. Glutathione S-transferase genes in the rice leaffolder, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae): identification and expression profiles. Arch. Insect Biochem Physiol. 2015;90(1): 1–13. pmid:25917811
- 33. Chen J, Li C, Yang ZF. Identification and expression of two novel cytochrome P450 genes, CYP6CV1 and CYP9A38, in Cnaphalocrocis medinalis (Lepidoptera: Pyralidae). J Insect Sci. 2015; 15(1): 1–7.
- 34.
Veegala RB, Vemuri S. Carboxylesterase and glutathione-S-transferase quantification mediating resistance in populations of rice leaf folder, Cnaphalocrocis medinalis, (Guenee). In: Chakravarthy AK, Editors. New Horizons in Insect Science: Towards Sustainable Pest Management. India: Springer; 2015. pp. 95–100.
- 35. Guo JY, Wu G, Wan FH. Temporal allocation of metabolic tolerance to transgenic Bt cotton in beet armyworm, Spodoptera exigua (Hübner). Sci China Life Sci. 2011; 54(2): 152–158. pmid:21318485
- 36. Shi GQ, Yu QY, Zhang Z. Annotation and evolution of the antioxidant genes in the silkworm, Bombyx mori. Arch Insect Biochem Physiol. 2012;79(2): 87–103. pmid:22392770
- 37. Zhang ZG, Sun JY, Li Y. Effects of Cry1Ab Protein from transgenic Bt rice on the activities of three protective enzymes in larvae of the leaf roller Cnaphalocrocis medinalis. Biotechnol Bull. 2011; 51(1): 148–152.
- 38. Pardo-López L, Soberón M, Bravo A. Bacillus thuringiensis insecticidal three-domain Cry toxins: mode of action, insect resistance and consequences for crop protection. FEMS Microbiol Rev. 2013;37(1): 3–22. pmid:22540421
- 39. Miranda R, Zamudio FZ, Bravo A. Processing of Cry1Ab delta-endotoxin from Bacillus thuringiensis by Manduca sexta and Spodoptera frugiperda midgut proteases: role in protoxin activation and toxin inactivation. Insect Biochem. Mol Biol. 2001;31(12): 1155–1163. pmid:11583928
- 40. Keller M, Sneh B, Strizhov N, Prudovsky E, Regev A, Koncz C, et al. Digestion of δ-entotoxin by gut proteases may explain reduced sensitivity of advanced instar larvae of Spodoptera littoralis to Cry1C. Insect Biochem Mol Biol. 1996;26(4): 365–373. pmid:8814783
- 41. Tetreau G, Stalinski R, David JP, Després L. Increase in larval gut proteolytic activities and Bti resistance in the Dengue fever mosquito. Arch Insect Biochem Physiol. 2013; 82(2): 71–83. pmid:23192850
- 42. Elleuch J, Zribi Zghal R, Lacoix MN, Chandre F, Tounsi S, Jaoua S. Evidence of two mechanisms involved in Bacillus thuringiensis israelensis decreased toxicity against mosquito larvae: Genome dynamic and toxins stability. Microbiol Res. 2015;176: 48–54. pmid:26070692
- 43. Tiewsiri K, Wang P. Differential alteration of two aminopeptidases N associated with resistance to Bacillus thuringiensis toxin Cry1Ac in cabbage looper. PNAS 2011;108(34): 14037–14042. pmid:21844358
- 44. English LH, Readdy TL. δ-Endotoxin inhibits a phosphatase in midgut epithelial membranes of Heliothis virescens. Insect Biochem. 1989;19(2):145–152.
- 45. Fernandez LE, Aimanova KG, Gill SS, Bravo A, Soberón M. A GPI-anchored alkaline phosphatase is a functional midgut receptor of Cry11Aa toxin in Aedes aegypti larvae. Biochem J. 2006; 394(Pt1): 77–84.
- 46. Arenas I, Bravo A, Soberon M, Gomez I. Role of alkaline phosphatase from Manduca sexta in the mechanism of action of Bacillus thuringiensis Cry1Ab toxin. J Biol Chem. 2010; 285(17): 12497–12503. pmid:20177063
- 47. Chen W, Liu C, Xiao Y, Zhang D, Zhang Y, Li X, et al. A toxin-binding alkaline phosphatase fragment synergizes Bt toxin Cry1Ac against susceptible and resistant Helicoverpa armigera. Plos One 2015;10: e0126288. pmid:25885820
- 48. Flores-Escobar B, Rodríguez-Magadan H, Bravo A, Soberón M, Gómez I. Differential role of Manduca sexta aminopeptidase N and alkaline phosphatase in the mode of action of Cry1Aa, Cry1Ab, and Cry1Ac toxins from Bacillus thuringiensis. Appl Environ Microbiol. 2013;79(15): 4543–4550. pmid:23686267
- 49. Xu XJ, Yu LY, Wu YD. Disruption of a cadherin gene associated with resistance to Cry1Ac d-endotoxin of Bacillus thuringiensis in Helicoverpa armigera. Appl Environ Microbiol. 2005; 71(2): 948–954. pmid:15691952
- 50.
Holland IB, Cole SP, Kuchler K, Higgins CF. ABC Proteins: from Bacteria to Man. Academic Press, 2003.
- 51. Shanker M,Willcutts D, Roth JA, Ramesh R. Drug resistance in lung cancer. Lung Cancer: Targets Ther. 2010; 1:23–26.
- 52. Strauss A.S., Peters S., Boland W., Burse A., 2013. ABC transporter functions as a pacemaker for sequestration of plant glucosides in leaf beetles. Elife 2, e01096. pmid:24302568
- 53. Ewart GD, Cannell D, Cox GB, Howells AJ. Mutational analysis of the traffic ATPase (ABC) transporters involved in uptake of eye pigment precursors in Drosophila melanogaster. Implications for structure-function relationships. J Biol Chem. 1994; 269(14): 10370–10377. pmid:8144619
- 54. Aurade RM, Jayalakshmi SK, Sreeramulu K. P-glycoprotein ATPase from the resistant pest, Helicoverpa armigera: purification, characterization and effect of various insecticides on its transport function. Biochim Biophys Acta Biomembr. 2010; 1798(6): 1135–1143.
- 55. Baxter SW, Badenes-Pérez FR, Morrison A, Vogel H, Crickmore N, Kain W, et al. Parallel evolution of Bacillus thuringiensis toxin resistance in Lepidoptera. Genetics 2011;189(2): 675–679. pmid:21840855
- 56. Gahan LJ, Pauchet Y, Vogel H, Heckel DG. An ABC transporter mutation is correlated with insect resistance to Bacillus thuringiensis Cry1Ac toxin. Plos Genet. 2010; 6: e1001248. pmid:21187898
- 57. Atsumi S, Miyamoto K, Yamamoto K, Narukawa J, Kawai S, Sezutsu H, et al. Single amino acid mutation in an ATP-binding cassette transporter gene causes resistance to Bt toxin Cry1Ab in the silkworm, Bombyx mori. PNAS 2012; 109(25): E1591–1598. pmid:22635270
- 58. Xiao Y, Zhang T, Liu CX, Heckel DG, Li XC, Tabashnik BE, et al. Mis-splicing of the ABCC2 gene linked with Bt toxin resistance in Helicoverpa armigera. Sci Rep. 2014; 4: 6184. pmid:25154974
- 59. Coates BS, Siegfried BD. Linkage of an ABCC transporter to a single QTL that controls Ostrinia nubilalis larval resistance to the Bacillus thuringiensis Cry1Fa toxin. Insect Biochem Mol Biol. 2015; 63: 86–96. pmid:26093031
- 60. Blum JS, Wearsch PA,Cresswell P. Pathways of antigen processing. Annu Rev Immunol. 2013; 31: 443–473. pmid:23298205
- 61. Oancea G, O’Mara ML, Bennett WF, Tieleman DP, Abele R, Tampe R. Structural arrangement of the transmission interface in the antigen ABC transport complex TAP. PNAS. 2009; 106(14):5551–5556. pmid:19297616
- 62. Park Y, Gonzálezmartínez RM, Navarrocerrillo G, Chakroun M, Kim Y, Ziarsolo P, et al. ABCC transporters mediate insect resistance to multiple Bt toxins revealed by bulk segregant analysis. BMC Biol. 2014; 12: 1–15.
- 63. Tabashnik BE. ABCs of Insect Resistance to Bt. PLoS Genet. 2015;11: e1005646. pmid:26583776
- 64. Corrêa S, Binato R, Du Rocher B, Morgana TL, Castelo-Branco MTL, Pizzatti L, et al. Wnt/β-catenin pathway regulates ABCB1 transcription in chronic myeloid leukemia. BMC Cancer 2012; 12: 303. pmid:22823957
- 65. Zhang X, Candas M, Griko NB, Rose-Young L, Bulla LA Jr. Cytotoxicity of Bacillus thuringiensis Cry1Ab toxin depends on specific binding of the toxin to the cadherin receptor BT-R1 expressed in insect cells. Cell Death Differ. 2005;12(11): 1407–1416. pmid:15920532
- 66. Zhang X, Candas M, Griko NB, Taussig R, Bulla LA Jr. A mechanism of cell death involving an adenylyl cyclase/PKA signaling pathway is induced by the Cry1Ab toxin of Bacillus thuringiensis. PNAS 2006;103(26): 9897–9902. pmid:16788061
- 67. Zhang X, Griko NB, Corona SK, Bulla L A Jr. Enhanced exocytosis of the receptor BT-R(1) induced by the Cry1Ab toxin of Bacillus thuringiensis directly correlates to the execution of cell death. Comp. Biochem Physiol B Biochem Mol Biol. 2008; 149(4): 581–588. pmid:18230416