Benshui Shu
Jingjing Zhang
Guohua Zhong- Key Laboratory of Crop Integrated Pest Management in South China, Ministry of Agriculture, Key Laboratory of Natural Pesticide and Chemical Biology, Ministry of Education, South China Agricultural University, Guangzhou, China
Azadirachtin, the environmentally friendly botanical pesticide, has been used as an antifeedant and pest growth regulator in integrated pest management for decades. It has shown strong biological activity against Spodoptera litura, but the mechanism of toxicity remains unclear. The present study showed that azadirachtin inhibited the growth of S. litura larvae, which was resulted by structure destroy and size inhibition of the midgut. Digital gene expression (DGE) analysis of midgut suggested that azadirachtin regulated the transcriptional level of multiple unigenes involved in mitogen-activated protein kinase (MAPK) and calcium apoptotic signaling pathways. Simultaneously, the expression patterns of some differentially expressed unigenes were verified by quantitative real time-PCR (qRT-PCR). In addition, the enhanced terminal deoxynucleotidyl transferase biotin-dUTP nick end labeling (TUNEL) staining, the increased expression of caspase family members and apoptosis-binding motif 1 (IBM1) on both gene and protein level and the release of cytochrome c from mitochondria to cytoplasm were induced in midgut after azadirachtin treatment. These results demonstrated that azadirachtin induced structural alteration in S. litura larval midgut by apoptosis activation. These alterations may affect the digestion and absorption of nutrients and eventually lead to the growth inhibition of larvae.
Introduction
The use of synthetic, chemical pesticides is the most effective method for agricultural pest control. However, the conventional pesticides could induce pest resistance (Tilman et al., 2002; Duhan et al., 2017). And the residues of pesticides also confer negative consequences on human health and environment (Cox and Surgan, 2006). Combined, these issues have created an urgent need for sustainable, effective pest control solutions (Helps et al., 2017; Moshi and Matoju, 2017; Patil et al., 2017). The botanical pesticides could be considered to be an alternative tool for integrated pest management due to the advantages of rapid degradation in environment, low toxicity to mammals, and a low risk of resistance development in target pest populations (Schmutterer, 1990; Greenberg et al., 2005; Ahmad et al., 2013).
Currently, the most promising, effective botanical insecticide for integrated pest management is azadirachtin, a tetranortriterpenoid extracted from Azadirachta indica (A. Juss). It has been widely used for decades and affected three aspects of insect biology: feeding behavior, growth, and development (Rembold et al., 1982; Isman et al., 1990; Wang et al., 2015). For example, azadirachtin exhibited a significant antifeedant effect on Drosophila melanogaster, Plutella xylostella, and Galleria mellonella (Sezer and Ozalp, 2011; Huang et al., 2016; Kilani-Morakchi et al., 2017). The antifeedant effect was mediated via the activation of deterrent receptor situated in the medial sensillum styloconicum (Liner et al., 1995). It could also inhibit the transmission of cholinergic nerve signaling and calcium channel in suboesophageal ganglion and reduce the frequency of miniature excitatory postsynaptic currents (mEPSCs). These reductions in nerve cell conductivity resulted in the dysfunction of pest central nervous system and promoting antifeedant behavior (Qiao et al., 2014). In addition to feeding behavior, azadirachtin also inhibited insect growth and development. In D. melanogaster, azadirachtin inhibited the development via blocking the biosynthesis of ecdysteroids and juvenile hormones (Lynn et al., 2012; Lai et al., 2014). It also blocked the release of prothoracicotropic hormone (PTTH) from neuroendocrine cells, which inhibited the development and molting of Trypanosoma cruzi (Garcia et al., 1990; Cortez et al., 2012). Besides, azadirachtin disrupted the proper functioning of endocrine and neuroendocrine systems in Labidura riparia (Sayah et al., 1998).
The tobacco cutworm, S. litura Fabricius (Lepidoptera: Noctuidae), is a polyphagous pest that feeds on more than 150 different host plants and widely distributed throughout tropical and subtropical regions (Gong et al., 2014; Selin-Rani et al., 2016). It has caused devastating destruction to many important field crops and vegetables, such as cotton, soybeans, and cabbage (Kaur et al., 2016; Kiran Gandhi et al., 2016). Azadirachtin has shown significant antifeedant and growth inhibitory action on S. litura at the concentration of 10–100 parts per million (ppm) and 1–10 ppm respectively (Govindachari et al., 1996). In addition, the fecundity was also reduced significantly (Nathan and Kalaivani, 2006). It was also demonstrated that azadirachtin had other significant physiological effects on S. litura, including deformity of larvae, pupae and adult, reduction of the protein synthesis of pupae, hemolymph volume of the last instar larvae, enzyme activities of larvae gut, and cuticular protein level changes of larvae (Sharma et al., 2003; Huang et al., 2004; Nathan et al., 2005; Jeyasankar et al., 2011; Yooboon et al., 2015). Furthermore, the plasma membrane damage and organelle degeneration in plasmatocytes and granular haemocytes were also induced by azadirachtin (Sharma et al., 2003).
The insect midgut functions as the important site for food digestion and nutritional absorption during insect growth and development. It was reported that azadirachtin could exert multiple effects on pest midgut. For example, azadirachtin reduced the levels and activities of midgut digestive enzymes of Glyphodes pyloalis (Khosravi and Sendi, 2013). Additionally, azadirachtin disturbed the serotoninergic system of stomatogastric ganglia and inhibited the peristalsis of midgut in Locusta migratoria (Trumm and Dorn, 2000). However, the mechanism of azadirachtin exposure to the midgut of S. litura has not been clearly defined. Recent study has found that azadirachtin induced apoptosis in SL-1 cell line, which was derived from the ovary of S. litura (Huang et al., 2011). More recently, autophagy-related gene 5 was confirmed to be the molecular switch of autophagy and apoptosis induced by azadirachtin in SL-1 cells (Shao et al., 2016). Therefore, we speculated that the mechanism of azadirachtin could be related to apoptosis in vivo.
In order to elucidate the adverse effects of azadirachtin on S. litura larvae, the present study investigated the morphology and histopathological changes of larvae after azadirachtin expose. Digital gene expression (DGE) analysis of midgut showed that multiple apoptotic signaling pathways involved in the process. In addition, qRT-PCR, TUNEL and western blot had been accomplished and further confirmed that azadirachtin induced apoptosis in the midgut of S. litura larvae. Our data could provide further insight into the mechanism of azadirachtin for growth regulation and it might benefit for the application of these efficient botanic pesticides.
Materials and Methods
Reagents and Antibodies
Azadirachtin (95% purity, #A7430, Sigma-Aldrich, St. Louis, MO, USA) was dissolved into dimethyl sulfoxide (DMSO, #D103272, Aladdin, Shanghai, China) and mixed with artificial diet. The antibody of Sl-IBM1 was prepared by our laboratory. Cleaved Caspase-3 (Asp 175) antibody was obtained from Cell Signaling Technology (#9661, USA) and mouse polyclonal anti-cytochrome c was purchased from Beyotime Biotechnology (#AC909, Shanghai, China). Other chemicals were domestic products.
Insect Culture and Treatment
The S. litura larvae maintained in our laboratory (Guangzhou, China) were fed with an artificial diet upon hatching and maintained in a stable condition of 25 ± 1°C, 60–70% relative humidity and a 16:8 h light: dark cycle (Yi et al., 2017). The adults were maintained with 10% honey water. The third-instar larvae fed with the artificial diet supplemented with 1 μg/g azadirachtin for 7 d were used for the following experiments and defined as larvae with azadirachtin treatment. The larvae fed with the artificial diet containing 200 μL/kg DMSO were used as the control group.
Hematoxylin–Eosin Staining
Azadirachtin-treated and control larvae were dissected and the midgut was washed in cold phosphate buffered saline (PBS) and fixed with 4% paraformaldehyde (#G1101, Servicebio, Wuhan, China) at 4°C for more than 24 h. Then it was embedded in paraffin wax, sectioned at 4 μm slices, mounted on glass slides and stained by hematoxylin and eosin. Slides were visualized on a microscope (Nikon, Japan).
cDNA Library Preparation and Illumina Sequencing
Total RNA was extracted from the midgut samples using TRIzol (#15596026, ThermoFisher Scientific, USA) as described in the manufacturer's instructions. Purity of total RNA was measured by a NanoDrop® spectrophotometer (Thermo Fihser, MA, USA). RNA integrity was determined by an Agilent 2100 (Agilent Technologies, CA, USA). Two libraries were constructed, one for the midgut of azadirachtin-treated larvae and one for the control. A total of 1 μg qualified RNA per sample was used for the library preparation. The sequencing libraries were generated by VAHTS mRNA-seq v2 Library Prep Kit for Illumina® (Catalog NR601, Vazyme, Nanjing, China) following manufacturer's protocol. The clustering of the index-coded samples was performed on a cBot Cluster Generation System (Illumia, USA) according to the manufacturer's instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq X Ten platform with 150 bp paired-end module. The de novo assembly of clean reads was performed by Trinity (settings: –min_contig_length150 –CPU 8 –min_kmer_cov 3 –min_glue 3 –bfly_opts ‘-V 5 –edge-thr = 0.1 –stderr’) (Grabherr et al., 2011). The expression levels of unigenes based on the number of reads uniquely mapped were normalized by Reads Per Kilobase of exon model per Million mapped reads (FPKM) method. The thresholds of differentially expressed unigene were False discovery rate (FDR) ≤ 0.001 and |log2Fold change| ≥1. Differentially expressed unigenes were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis.
Quantitative Real Time PCR (qRT-PCR)
The total RNA of midgut samples was extracted by TRIzol reagent as above and reverse transcribed into cDNA with the PrimeScript™ RT reagent Kit (#RR047A, TaKaRa, Japan) following the manufacturer's recommendations. qRT-PCR was performed using the iTaq™ Universal SYBR® Green Supermix (#1725271, BIO-RAD, USA) with a CFX Connect™ Real-Time System (BIO-RAD, USA) under the following thermal program: one cycle at 95°C for 3 min, 40 cycles at 95°C for 10 s; 60°C for 10 s; 72°C for 15 s, and one step from 65 to 95°C for the dissociation stage. The primers used for qRT-PCR were listed in Supplement Table 1. Gene expression levels were calculated by 2−ΔΔCt method (Livak and Schmittgen, 2001) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as the reference gene (Lu et al., 2013).
Terminal Deoxynucleotidyl Transferase Biotin-Dutp Nick End Labeling (TUNEL) Assay
TUNEL assay is one of the most effective methods for apoptosis detection in tissues (Sarkissian et al., 2014). Paraffin embedded midgut tissue section was prepared as described above. One Step TUNEL Apoptosis Assay Kit (#C1086, Beyotime, Shanghai, China) was used according to the manufacturer's protocol. Briefly, sections were dewaxed in xylene twice for 5 min, fixed with 100% ethanol for 5 min, then dehydrated with 2 min washes in 90% ethanol, 70% ethanol, distilled water. After that, the sections were treated with 20 μg/mL protease K (#1074, Servicebio, Wuhan, China) at 37°C for 30 min and washed three times with PBS. Slides were incubated with the TUNEL reagent in the dark at 37°C for 1 h, washed three times with PBS and counterstained with 4′, 6-diamidino-2-phenylindole (DAPI) prior to mounting in X mounting media. Slides were visualized on a fluorescence microscope (Nikon, Japan).
Western Blot
Whole protein extraction from midgut samples were carried out with the total protein extraction kit (#W034, Nanjing Jiancheng, China) according to the protocol. Cytoplasmic and mitochondrial proteins were extracted using the mitochondrial protein extraction kit (#KGP850, KeyGEN BioTECH, China). Samples were separated on a 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE) gel and transferred to a polyvinylidene fluoride (PVDF) (#IPVH00010, Millipore, USA) with the wet/semi-dry transfer system. The membranes were blocked in tris-buffered saline (TBS: 100 mM Tris–HCl, pH 7.5, 0.9% NaCl) with 5% fat-free milk at room temperature for 2 h and then incubated in primary antibodies (diluted 1:2000 in TBS) at 4°C overnight, washed by TBST (TBS containing 0.1% Tween-20) three times for 5 min and incubated in secondary antibody (1:2000 in TBS) at room temperature for 2 h. Membranes were developed in enhanced chemiluminescence (ECL) reagent (#1705062S, Bio-RAD, USA) and bands were visualized using ECL detection system (Bio-RAD, USA). GAPDH was used as the reference protein to normalize the differences in protein loading.
Caspase-3 Activity Assay
Approximately 3–10 mg of midgut samples dissected from larvae were prepared. The protein concentrations were determined with Bradford method. The caspase-3 activity was detected by the Caspase-3 Activity Assay Kit (#C1115, Beyotime, Shanghai, China) according to the manufacturer's protocol. Enzyme activities were measured by the absorbance of samples at 405 nm with a microplate reader (ThermoFisher Scietific, USA).
Statistical Analysis
Data collected were expressed as the mean ± standard deviation. Statistical significance between control and treatment was determined with SPSS 20.0 using Student's t-test.
Results
Inhibitory Effect of Azadirachtin on the Growth of S. litura Larvae
To confirm the growth inhibition phenotype conferred by azadirachtin, third instar larvae were exposed to azadirachtin added or control diet. After 7 days, the larvae fed with azadirachtin exhibited a growth inhibition phenotype. Specifically, the size of larvae with azadirachtin treatment was obviously smaller (Figure 1A). Additionally, the larvae fed with azadirachtin weighed significantly less, with a 43.4% decrease (1.93 g ± 0.093 vs. 3.41 g ± 0.234, P < 0.05, n = 20) (Figure 1B).
Figure 1. Morphological and dietary changes of S. litura larvae with azadirachtin treatment for 7 days. Control: larvae fed with a DMSO added diet; Azadirachtin: larvae fed with azadirachtin added diet. (A) Morphological change of larvae fed with azadirachtin added diet; (B) dietary change of larvae fed with azadirachtin added diet.
Morphology and Histopathological Examination of Midgut after Azadirachtin Treatment
In order to determine the effect of azadirachtin on structure, the midgut samples of larvae were dissected, measured and stained with hematoxylin-eosin. After azadirachtin exposure, the size of larval midgut was significantly smaller (Figure 2A). In addition to gross structural changes, azadirachtin also conferred histopathological changes. The midgut of control displayed appropriate, tightly arranged cells (Figure 2B). However, it exhibited cells death, abnormal cell structure, and intestinal wall cracking after azadirachtin treatment (Figure 2B). These data demonstrated that azadirachtin disrupted the normal cell structure and induced cell death in the midgut.
Figure 2. Morphological change and histopathological observation of midgut in S. litura larvae fed with azadirachtin added diet. (A) Morphological change of midgut in larvae fed with azadirachtin added diet; (B) Histopathological observation of midgut in larvae fed with azadirachtin added diet. (B1) Control: control group with the normal condition; (B2) Azadirachtin: azadirachtin-treated group with significant histopathological change. Scale bar = 100 μm.
Differentially Expressed Genes Detected by Digital Gene Expression (DGE) Analysis
In order to uncover the molecular mechanisms underlying the histological changes in the midgut of azadirachtin-treated larvae, transcriptional differences were measured by DGE profiling. As shown in Figure 3, a total of 4,091 unigenes were differentially expressed in midgut of azadirachtin treatment. Of these unigenes, 3,108 were up-regulated and 983 were down-regulated.
Figure 3. The statistics of differentially expressed unigenes.
The GO enrichment analysis indicated that the differentially expressed unigenes (DEGs) affected by azadirachtin were classified into three categories and 57 functional groups. Among which 1,181 DEGs were assigned into the category of biological process, 915 DEGs to cellular component and 1,177 DEGs to molecular function (Supplement Figure 1). In biological process category, the DEGs assigned in were divided into 24 functional groups and single-organism process was the biggest group with 881 genes assigned, followed by the group of cellular process (881 DEGs), and metabolic process (740 DEGs). As to cellular component category, the biggest group were cell and cell part, which has the same number of 629 unigenes assigned. The DEGs classified into molecular function category were divided into 16 groups and catalytic activity group has the higher percentage of unigenes and 723 unigenes were concentrated.
The DGE analysis revealed 2,154 differentially expressed unigenes were involved in 247 KEGG pathways. Of these pathways, three pathways with the largest number of annotations were metabolic pathways (358 unigenes, 16.62%), pancreatic secretion (181 unigenes, 8.4%), and protein digestion and absorption (148 unigenes, 6.87%) (Supplement Figure 2).
qRT-PCR Validation
To verify the results of DGE analysis, parts of unigenes with at least a two-fold difference in expression were selected. Because azadirachtin was considered to be a xenobiotic, 17 differentially expressed unigenes involved in xenobiotics biodegradation and metabolism pathway were validated by qRT-PCR (Table 1) and exhibited the same trends as the DGE profiling results (Figure 4A). These results indicated that three major detoxification enzymes (esterase, UDP-glycosyltransferase and cytochrome P450) were activated in midgut after exposure to azadirachtin.
Table 1. The differentially expressed unigenes validated by qRT-PCR.
Figure 4. Validation of differentially expressed unigenes expression levels by qRT-PCR. Data were expressed as arithmetic mean of three replications. (A) 17 differentially expression genes annotated into drug metabolic pathways. (B) 21 differentially expressed genes annotated into MAPK signaling pathway. (C) 20 differentially expressed genes in calcium signaling pathway. (D) 14 differentially expressed genes in other signaling pathway. All data were normalized by the expression level of GAPDH.
In addition, many differentially expressed unigenes were mapped into the pathways involved in apoptosis, such as the mitogen-activated protein kinase (MAPK), calcium, and p53 signaling pathways. A total of 55 unigenes in these apoptotic pathways were also validated by qRT-PCR (Table 1). All the selected unigenes expression profiles were consistent with the results of DGE analysis (Figures 4B–D). These results suggested that azadirachtin activated MAPK, Calcium and other apoptosis signaling pathways within larval midgut.
Detection of Azadirachtin-Induced Apoptosis in the Midgut by TUNEL Assay
TUNEL assay was utilized to directly demonstrate that azadirachtin induced apoptosis and morphological changes in midgut. The TUNEL-positive cells would fluoresce at 488 nm. TUNEL staining was undetectable in midgut of control larvae, while a high level of TUNEL-positive cells appeared in midgut of azadirachtin exposure (Figure 5). These results indicated that azadirachtin induced apoptosis in midgut.
Figure 5. The apoptosis determination by TUNEL assay and DAPI staining in the midgut of S. litura larvae fed with azadirachtin added diet. TUNEL-positive cell which was defined as the cell with the green fluorescence. The cell nuclei stained with DAPI was showed with the blue color.
Azadirachtin Could Up-Regulate the Expression of SL-IBM1
Many unigenes were up-regulated in DGE of azadirachtin-treated midgut. One such differentially expressed unigene IAP-binding motif 1 (Sl-IBM1) (log2 Ratio = 4.11), an IAP antagonist and Lepidoptera order homolog to Reaper, was thought to play an important role in apoptosis process (Wu et al., 2013). Therefore, the expression of IBM1 in midgut after azadirachtin treatment was verified by qRT-PCR and western blot. Azadirachtin exposure caused an up-regulation of Sl-IBM1 mRNA expression in midgut (Figure 6A). In addition, it could also induce an increase of Sl-IBM1 on protein level (Figure 6B). Taken together, azadirachtin exposure up-regulated the pro-apoptotic factor IBM1.
Figure 6. The effects of azadirachtin on mRNA and protein level of IBM1 in larval midgut. (A) The mRNA expression change of IBM1 in midgut of larvae fed with azadirachtin added diet. Data were expressed as arithmetic mean ± SEM of three independent experiments. *Significant difference (P < 0.05), **significant difference (P < 0.01) by ANOVA followed by student's t-test. (B) The Sl-IBM1 protein level in midgut of larvae fed with azadirachtin and DMSO added diet.
Azadirachtin Up-Regulated the Expression of Caspases and Caspase-3-Like Activity
To determine if azadirachtin exposure and the increased Sl-IBM1 expression induced caspase expression within the midgut, mRNA expression of four caspase genes, including caspase-1, caspase-3, caspase-5, and caspase-6 were compared after azadirachtin exposure by qRT-PCR. All four genes were demonstrated to be up-regulated at least 2.75 folds after azadirachtin treatment (Figure 7A). Simultaneously, the cleaved caspase 3 protein was increased significantly (Figure 7B). Additionally, in accordance with gene and protein expression data, caspase-3-like protease activity was obviously increased (Figure 7C). The up-regulation of caspases expression and caspase-3-like activity confirmed that azadirachtin activated apoptosis in midgut by regulating the caspase-dependent apoptotic pathway.
Figure 7. The effects of azadirachtin on mRNA level of caspase family, protein level of cleaved caspase 3 and caspase-3 like activity. (A) The mRNA expression changes of caspase family in midgut of larvae fed with azadirachtin added diet. Data were expressed as arithmetic mean ± SEM of three independent experiments. *Significant difference (P < 0.05), **significant difference (P < 0.01) by ANOVA followed by student's t-test. (B) The cleaved caspase 3 protein level in midgut of larvae fed with azadirachtin and DMSO added diet. (C) The caspase-3 like activity change in midgut of larvae fed with azadirachtin diet.
Azadirachtin Induced the Release of Cytochrome c in Midgut Cells
To reveal whether mitochondria were involved in midgut apoptosis induced by azadirachtin, the protein levels of cytochrome c within cytoplasm and mitochondria of midgut cells were analyzed by western blot. After azadirachtin treatment, cytochrome c in cytoplasmic increased significantly while the mitochondrial reduced (Figure 8). These results indicated that the event of cytochrome c releasing from mitochondria into cytoplasm was occurred in midgut after azadirachtin treatment.
Figure 8. The changes of cytochrome c expression level in cytosolic and mitochondrial of midgut in S. litura larvae fed with azadirachtin diet. Control: larvae fed with DMSO added diet; Azadirachtin: larvae fed with azadirachtin added diet.
Discussion
Because of the botanical efficacy, multiple studies have investigated the toxicity of azadirachtin against various insect species. Prior research analyzed the mechanism of the antifeedant activity and regulation of growth in vitro and apoptosis induction in vivo (Huang et al., 2010; Shu et al., 2015; Xu et al., 2016). Consistent with the prior research, our results demonstrated azadirachtin inhibited the growth of S. litura larvae. The insect midgut plays a crucial role in the physiology of growth, including food intake, digestion, and nutrient absorption (Franzetti et al., 2015). It was reported that azadirachtin immediately and severely diminished the conversion of ingested nutrients in Spodoptera littorali (Martinez and van Emden, 1999). This conversion was showed to be due, in part, to azadirachtin-induced inhibition of the activity of α-amylase, a digestive enzyme produced by midgut epithelial cells (Rharrabe et al., 2008). Our results revealed that azadirachtin also disrupted the structure of larval midgut epithelium. This destruction likely further disturb the digestion of nutrients, which may contribute to explain the growth regulation mechanisms of azadirachtin in S. litura.
Apoptosis is a spontaneous programmed cell death process in response to adverse stimuli (Suganuma et al., 2011; Huang N. et al., 2013). Pesticides could work in part by activating apoptosis. For example, apoptosis occurred as a protective mechanism in midguts, salivary glands and ovaries of honey bee larvae after pesticides treatment (Gregorc and Ellis, 2011). Additionally, apoptosis and autophagy were significantly increased in the brain of worker honey bee (Apis mellifera) after exposure to sublethal doses of imidacloprid (Wu et al., 2015). Multiple botanical pesticides also exhibited the same function of apoptosis induction on a variety of insect species. For example, the cry toxins, produced by Bacillus thuringiensis, induced mitochondrial permeabilization and led to apoptosis in Aedes aegypti larvae (Lemeshko and Orduz, 2013). In addition, previous research of our laboratory has found that camptothecin resulted in midgut epithelial cell apoptosis of S. litura, which was correlated with the increased expression of programmed cell death protein 11 (Gong et al., 2014). Our results confirmed that azadirachtin also induced apoptosis in midgut epithelial cells of S. litura. These results further supported apoptosis induction in response to pesticide exposure.
Both MAPK and calcium signaling pathways regulated multiple biological processes, including cell proliferation, differentiation and apoptosis (Berridge et al., 2000; Kyriakis and Avruch, 2012; Zhu et al., 2017). Neem limonoids contained azadirachtin have previously been reported to target the MAPK signaling pathway (Nagini, 2014; Sui et al., 2014). Additionally, in D. melanogaster S2 cells, azadirachtin induced apoptosis by regulating the Ca2+-Calmodulin signaling pathway and mediating the release of intracellular Ca2+ (Xu et al., 2016; Humeau et al., 2017). In this study, DGE analysis and qRT-PCR indicated that genes involved in MAPK and calcium signaling pathways were differentially expressed in midgut after azadirachtin exposure. Our data suggested that azadirachtin may also act, in part, through MAPK and calcium signaling pathways to induce apoptosis. Further validation is required to link these pathways to apoptosis activation in response to azadirachtin in vivo.
Multiple insect species have been shown to activate apoptosis using homologous pathways in response to negative stimuli. In insects, RHG protein family members and inhibitor of apoptosis (IAP) proteins regulated apoptosis and cell death (Bryant et al., 2009). Additionally, the transcriptional expression of reaper, hid, grim preceded apoptosis (Vasudevan and Ryoo, 2015). In Drosophila, Reaper induced-apoptosis in response to stimuli required both homodimerization and oligimerization with Hid (Sandu et al., 2010). The oligimerization of Reaper and Hid resulted in the recruitment to mitochondrial membrane and activation of apoptosis (Sandu et al., 2010). Besides, in BmN cells and Bombyx mori pupae, the mRNA expression level of Bm-IBM1 (the Reaper homolog) increased in response to infection with baculovirus B. mori Nucleopolyhedrovirus (Wu et al., 2013). In addition, Bm-IBM1 bind to BmIAP1 at the highly conserved IAP-binding motif, which resulted in apoptosis induction (Bryant et al., 2009). Furthermore, Bm-IBM1 induced apoptosis by localizing to mitochondria and activating apoptosis signaling pathway, demonstrating the conservation of Reaper-induced apoptosis across species (Wu et al., 2013). Our data demonstrated that the expression of Sl-IBM1 was also increased in midgut after azadirachtin exposure, which suggested that the transcription level changes of IBM1 in lepidoptera signal apoptosis induction could be one of the indices of apoptosis. Possibly, Sl-IBM1 interacted with IAPs to activate the mitochondrial-associated apoptotic pathways, although no Hid-like protein homolog was found in S. litura.
Azadirachtin induced apoptosis in insect cell lines has been studied more detail over past decades. In Sl-1 cell lines, azadirachtin induced the caspase-dependent mitochondrial pathway via activating the cleavage of caspase-3 (Huang et al., 2011). The release of cytochrome c from mitochondria to cytoplasm is the critical step in mitochondrial apoptotic pathway and the cytochrome c released in cytosol bound to apaf-1 and formed an apoptosome, eventually activated caspases and resulted in apoptosis (Liu et al., 2012; Huang J. et al., 2013; Yang et al., 2017). In this study, the increase of caspases mRNA expression, the cleaved caspase 3 protein expression, the up-regulation of caspase-3 like activity and the release of mitochondrial cytochrome c in vivo further supported that azadirachtin activated apoptosis via the caspase-dependent mitochondrial apoptotic pathway. Additionally, these results indicated that the mitochondrial apoptotic pathway could be prevalent and plays the critical role in apoptosis activation of lepidoptera insects.
In conclusion, azadirachtin inhibited the growth of S. litura larvae by inducing apoptosis and destroying structure of the midgut. DGE analysis suggested that azadirachtin regulated both MAPK and calcium signaling pathways. It also activated the caspase-dependent mitochondrial apoptotic pathway and induced the release of cytochrome c in larval midgut. Our results are the first to demonstrate that azadirachtin caused growth restriction of S. litura larvae via apoptosis activation of midgut in vivo.
Author Contributions
BS: Conceived and designed the experiments; BS, JZ, GC, and RS: Performance of the experiments; BS and JZ: Analysis of the data; BS: Manuscript writing; XY and GZ: Manuscript editing.
Funding
This study was supported by a grant from the National Natural Science Foundation of China (Grants No. 31572335 and 31701812), Public welfare industry (Agriculture) scientific research special fund project, China (Grant No. 201303017) and Science and Technology Planning Project of Guangdong Province (No. 2016A020210090 and No. 2017A010105023).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00137/full#supplementary-material
Supplement Figure 1. Gene Ontology (GO) classification of differentially expressed unigenes. Three different classification represent the three basic categories of Go term (The categories from left to right were biological process, cell composition, molecular functions).
Supplement Figure 2. Top 20 pathways enrichment for differentially expressed unigenes. The y-axis indicates the 20 KEGG pathways and the x-axis indicates the richfactor of each KEGG pathway.
Supplement Table 1. Primers used in the paper.
References
Ahmad, S., Ansari, M. S., and Moraiet, M. A. (2013). Demographic changes in Helicoverpa armigera after exposure to neemazal (1% EC azadirachtin). Crop Prot. 50, 30–36. doi: 10.1016/j.cropro.2013.03.012
Berridge, M. J., Lipp, P., and Bootman, M. D. (2000). The versatility and universality of calcium signalling. Nat. Rev. Mol. Cell Biol. 1, 11–21. doi: 10.1038/35036035
Bryant, B., Zhang, Y., Zhang, C., Santos, C. P., Clem, R. J., and Zhou, L. (2009). A lepidopteran orthologue of IBM1 reveals functional conservation and evolution of IAP antagonists. Insect Mol. Biol. 18, 341–351. doi: 10.1111/j.1365-2583.2009.00878.x
Cortez, M. R., Provençano, A., Silva, C. E., Mello, C. B., Zimmermann, L. T., Schaub, G. A., et al. (2012). Trypanosoma cruzi: effects of azadirachtin and ecdysone on the dynamic development in Rhodnius prolixus larvae. Exp. Parasitol. 131, 363–371. doi: 10.1016/j.exppara.2012.05.005
Cox, C., and Surgan, M. (2006). Unidentified inert ingredients in pesticides: implications for human and environmental health. Environ. Health Persp. 114, 1803–1806. doi: 10.1289/ehp.9374
Duhan, J. S., Kumar, R., Kumar, N., Kaur, P., Nehra, K., and Duhan, S. (2017). Nanotechnology: the new perspective in precision agriculture. Biotechnol. Rep. 15, 11–23. doi: 10.1016/j.btre.2017.03.002
Franzetti, E., Romanelli, D., Caccia, S., Cappellozza, S., Congiu, T., Rajagopalan, M., et al. (2015). The midgut of the silkmoth Bombyx mori is able to recycle molecules derived from degeneration of the larval midgut epithelium. Cell Tissue Res. 361, 509–528. doi: 10.1007/s00441-014-2081-8
Garcia, E. S., Luz, N., Azambuja, P., and Rembold, H. (1990). Azadirachtin depresses release of prothoracicotropic hormone in Rhodnius prolixus larvae: evidence from head transplantations. J. Insect Physiol. 36, 679–682. doi: 10.1016/0022-1910(90)90073-O
Gong, L., Wang, H., Huang, J., Hu, M., Hu, Z., and Zhong, G. (2014). Camptothecin-induced expression of programmed cell death gene 11 in Spodoptera litura. Pest Manag. Sci. 70, 603–609. doi: 10.1002/ps.3588
Govindachari, T. R., Narasimhan, N. S., Suresh, G., Partho, P. D., and Gopalakrishnan, G. (1996). Insect antifeedant and growth-regulating activities of Salannin and other c-seco limonoids from neem oil in relation to Azadirachtin. J. Chem. Ecol. 22, 1453–1461. doi: 10.1007/BF02027724
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Greenberg, S. M., Showler, A. T., and Liu, T. X. (2005). Effects of neem-based insecticides on beet armyworm (Lepidoptera: Noctuidae). Insect Sci. 12, 17–23. doi: 10.1111/j.1672-9609.2005.00003.x
Gregorc, A., and Ellis, J. D. (2011). Cell death localization in situ in laboratory reared honey bee (Apis mellifera L.) larvae treated with pesticides. Pestic. Biochem. Phys. 99, 200–207. doi: 10.1016/j.pestbp.2010.12.005
Helps, J. C., Paveley, N. D., and van den Bosch, F. (2017). Identifying circumstances under which high insecticide dose increases or decreases resistance selection. J. Theor. Biol. 100, 196–202. doi: 10.1016/j.jtbi.2017.06.007
Huang, J., Lv, C., Hu, M., and Zhong, G. (2013). The mitochondria-mediate apoptosis of lepidopteran cells induced by azadirachtin. PLoS ONE 8:e58499. doi: 10.1371/journal.pone.0058499
Huang, J., Shui, K., Li, H., Hu, M., and Zhong, G. (2011). Antiproliferative effect of azadirachtin A on Spodoptera litura Sl-1 cell line through cell cycle arrest and apoptosis induced by up-regulation of p53. Pestic. Biochem. Phys. 99, 16–24. doi: 10.1016/j.pestbp.2010.08.002
Huang, N., Civciristov, S., Hawkins, C. J., and Clem, R. J. (2013). SfDronc, an initiator caspase involved in apoptosis in the fall armyworm Spodoptera frugiperda. Insect Biochem. Mol. Biol. 43, 444–454. doi: 10.1016/j.ibmb.2013.02.005
Huang, T. T., Ali, S., and Long, X. Z. (2016). Joint application of Isaria fumosorsea chitinase and azadirachtin as biopesticides against the diamondback moth, Plutella xylostella L. (Lepidoptera: Plutellidae). Egypt. J. Biol. Pest Control 26, 539–543.
Huang, X. Y., Li, O. W., and Xu, H. H. (2010). Induction of programmed death and cytoskeletal damage on Trichoplusia ni BTI-Tn-5B1-4 cells by azadirachtin. Pestic. Biochem. Phy. 98, 289–295. doi: 10.1016/j.pestbp.2010.06.020
Huang, Z., Shi, P., Dai, J., and Du, J. (2004). Protein metabolism in Spodoptera litura (F.) is influenced by the botanical insecticide azadirachtin. Pestic. Biochem. Phys. 80, 85–93. doi: 10.1016/j.pestbp.2004.07.001
Humeau, J., Bravo-San Pedro, J. M., Vitale, I., Nuñez, L., Villalobos, C., Kroemer, G., et al. (2017). Calcium signaling and cell cycle: progression or death. Cell Calcium 17, S0143–S4160. doi: 10.1016/j.ceca.2017.07.006
Isman, M. B., Koul, O., Luczynski, A., and Kaminski, J. (1990). Insecticidal and antifeedant bioactivities of neem oils and their relationship to azadirachtin content. J. Agric. Food Chem. 38, 1406–1411. doi: 10.1021/jf00096a024
Jeyasankar, A., Raja, N., and Ignacimuthu, S. (2011). Insecticidal compound isolated from Syzygium lineare Wall. (Myrtaceae) against Spodoptera litura (Lepidoptera: Noctuidae). Saudi J. Biol. Sci. 18, 329–332. doi: 10.1016/j.sjbs.2011.01.003
Kaur, M., Kumar, R., Upendrabhai, D. P., Singh, I. P., and Kaur, S. (2016). Impact of sesquiterpenes from Inula racemosa (Asteraceae) on growth, development and nutrition of Spodoptera litura (Lepidoptera: Noctuidae). Pest Manag. Sci. 73, 1031–1038. doi: 10.1002/ps.4429
Khosravi, R., and Sendi, J. J. (2013). Effect of Neem pesticide (Achook) on midgut enzymatic activities and selected biochemical compounds in the hemolymph of lesser mulberry pyralid, Glyphodes pyloalis Walker (Lepidoptera: Pyralidae). J. Plant Prot. Res. 53, 238–247. doi: 10.2478/jppr-2013-0036
Kilani-Morakchi, S., Bezzar-Bendjazia, R., Ferdenache, M., and Aribi, N. (2017). Preimaginal exposure to azadirachtin affects food selection and digestive enzymes in adults of Drosophila melanogaster (Diptera: Drosophilidae). Pestic. Biochem. Physiol. 140, 58–64. doi: 10.1016/j.pestbp.2017.06.004
Kiran Gandhi, B., Patil, R. H., and Srujana, Y. (2016). Field resistance of Spodoptera litura (Fab.) to conventional insecticides in India. Crop Prot. 88, 103–108. doi: 10.1016/j.cropro.2016.06.009
Kyriakis, J. M., and Avruch, J. (2012). Mammalian MAPK signal transduction pathways activated by stress and inflammation: a 10-year update. Physiol. Rev. 92, 689–737. doi: 10.1152/physrev.00028.2011
Lai, D., Jin, X., Wang, H., Yuan, M., and Xu, H. (2014). Gene expression profile change and growth inhibition in Drosophila larvae treated with azadirachtin. J. Biotechnol. 185, 51–56. doi: 10.1016/j.jbiotec.2014.06.014
Lemeshko, V. V., and Orduz, S. (2013). Electrical hypothesis of toxicity of the Cry toxins for mosquito larvae. Biosci. Rep. 33, 125–136. doi: 10.1042/BSR20120101
Liner, L., Vanloon, J. J. A., and Schoonhoven, L. M. (1995). Behavioural and sensory responses to some neem compounds by Pieris brassicae larvae. Physiol. Entomol. 20, 134–140. doi: 10.1111/j.1365-3032.1995.tb00809.x
Liu, K. Y., Yang, H., Peng, J. X., and Hong, H. Z. (2012). Cytochrome c and insect cell apoptosis. Insect Sci. 19, 30–40. doi: 10.1111/j.1744-7917.2011.01431.x
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCt method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Lu, Y., Yuan, M., Gao, X., Kang, T., Zhan, S., Wan, H., et al. (2013). Identification and validation of reference genes for gene expression analysis using quantitative PCR in Spodoptera litura (Lepidoptera: Noctuidae). PLoS ONE 8:e68059. doi: 10.1371/journal.pone.0068059
Lynn, O. M., Kim, J. E., and Lee, K. Y. (2012). Effects of azadirachtin on the development and gene expression of fifth instar larvae of Indianmeal moth, Plodia interpunctella. J. Asia Pac. Entomol. 15, 101–105. doi: 10.1016/j.aspen.2011.08.008
Martinez, S. S., and van Emden, H. F. (1999). Sublethal concentrations of azadirachtin affect food intake, conversion efficiency and feeding behaviour of Spodoptera littoralis (Lepidoptera: Noctuidae). Bull. Entomol. Res. 89, 65–71. doi: 10.1017/S0007485399000085
Moshi, A. P., and Matoju, I. (2017). The status of research on and application of biopesticides in Tanzania. Rev. Crop Prot. 92, 16–28. doi: 10.1016/j.cropro.2016.10.008
Nagini, S. (2014). Neem Limonoids as anticancer agents: modulation of cancer hallmarks and oncogenic signaling. Enzymes 36, 131–147. doi: 10.1016/B978-0-12-802215-3.00007-0
Nathan, S. S., and Kalaivani, K. (2006). Combined effects of azadirachtin and nucleopolyhedrovirus (SpltNPV) on Spodoptera litura Fabricius (Lepidoptera: Noctuidae) larvae. Biol. Control 39, 96–104. doi: 10.1016/j.biocontrol.2006.06.013
Nathan, S. S., Kalaivani, K., and Chung, P. G. (2005). The effects of azadirachtin and nucleopolyhedrovirus on midgut enzymatic profile of Spodoptera litura Fab. (Lepidoptera: Noctuidae). Pestic. Biochem. Phys. 83, 46–57. doi: 10.1016/j.pestbp.2005.03.009
Patil, S. B., Goyal, A., Chitgupekar, S. S., Kumar, S., and El-Bouhssini, M. (2017). Sustainable management of chickpea pod borer. A review. Agron. Sustain. Dev. 37:e20. doi: 10.1007/s13593-017-0428-8
Qiao, J., Zou, X., Lai, D., Yan, Y., Wang, Q., Li, W., et al. (2014). Azadirachtin blocks the calcium channel and modulates the cholinergic miniature synapticcurrent in the central nervous system of Drosophila. Pest Manag. Sci. 70, 1041–1047. doi: 10.1002/ps.3644
Rembold, H., Sharma, G. K., Gzoppelt, C., and Schmutterer, H. (1982). Azadirachtin: apotent insect growth regulator of plant origin. J. Appl. Entomol. 93, 12–17. doi: 10.1111/j.1439-0418.1982.tb03564.x
Rharrabe, K., Amri, H., Bouayad, N., and Sayah, F. (2008). Effects of azadirachtin on post-embryonic development, energy reserves and α-amylase activity of Plodia interpunctella Hübner (Lepidoptera: Pyralidae). J. Stored Prod. Res. 44, 290–294. doi: 10.1016/j.jspr.2008.03.003
Sandu, C., Ryoo, H. D., and Steller, H. (2010). Drosophila IAP antagonists from multimeric complexes to promote cell death. J. Cell Biol. 190, 1039–1052. doi: 10.1083/jcb.201004086
Sarkissian, T., Timmons, A., Arya, R., Abdelwahid, E., and White, K. (2014). Detecting apoptosis in Drosophila tissues and cells. Methods 68, 89–96. doi: 10.1016/j.ymeth.2014.02.033
Sayah, F., Idaomar, M., Soranzo, L., and Karlinsky, A. (1998). Endocrine and neuroendocrine effects of azadirachtin in adult females of the earwig Labidura riparia. Tissue Cell 30, 86–94. doi: 10.1016/S0040-8166(98)80009-8
Schmutterer, H. (1990). Properties and potential of natural pesticides from the neem tree, Azadirachta indica. Ann. Rev. Entomol. 35, 271–297. doi: 10.1146/annurev.en.35.010190.001415
Selin-Rani, S., Senthil-Nathan, S., Thanigaivel, A., Vasantha-Srinivasan, P., Edwin, E. S., Ponsankar, A., et al. (2016). Toxicity and physiological effect of quercetin on generalist herbivore, Spodoptera litura Fab. and a non-target earthworm Eisenia fetida Savigny. Chemosphere 165, 257–267. doi: 10.1016/j.chemosphere.2016.08.136
Sezer, B., and Ozalp, P. (2011). The effects of azadirachtin on the percentage of glycogen contents in larvae of Galleria mellonella. Ekoloji 20, 67–72. doi: 10.5053/ekoloji.2011.8110
Shao, X., Lai, D., Zhang, L., and Xu, H. (2016). Induction of autophagy and apoptosis via PI3K/AKT/TOR pathway by azadirachtin A in Spodoptera litura cells. Sci. Rep. 18:35482. doi: 10.1038/srep35482
Sharma, P. R., Sharma, O. P., and Saxena, B. P. (2003). Effect of Neem gold on haemocytes of the tobacco armyworm, Spodoptera litura (Fabricius) (Lepidoptera; Noctuidae) Curr. Sci. India 84, 690–695.
Shu, B., Wang, W., Hu, Q., Huang, J., Hu, M., and Zhong, G. (2015). A comprehensive study on apoptosis induction by azadirachtin in Spodoptera frugiperda cultured cell line Sf9. Arch. Insect Biochem. Physiol. 89, 153–168. doi: 10.1002/arch.21233
Suganuma, I., Ushiyama, T., Yamada, H., Iwamoto, A., Kobayashi, M., and Ikeda, M. (2011). Cloning and characterization of a dronc homologue in the silkworm, Bombyx mori. Insect Biochem. Mol. Biol. 41, 909–921. doi: 10.1016/j.ibmb.2011.08.005
Sui, X., Kong, N., Ye, L., Han, W., Zhou, J., Zhang, Q., et al. (2014). p38 and JNK MAPK pathways control the balance of apoptosis and autophagy in response to chemotherapeutic agents. Cancer Lett. 344, 174–179. doi: 10.1016/j.canlet.2013.11.019
Tilman, D., Cassman, K. G., Matson, P. A., Naylor, R., and Polasky, S. (2002). Agricultural sustainability and intensive production practices. Nature 418, 671–677. doi: 10.1038/nature01014
Trumm, P., and Dorn, A. (2000). Effects of azadirachtin on the regulation of midgut peristalsis by the stomatogastric nervous system in Locusta migratoria. Phytoparasitica 28, 7–26. doi: 10.1007/BF02994020
Vasudevan, D., and Ryoo, H. D. (2015). Regulation of cell death by IAPs and their antagonists. Curr. Top. Dev. Biol. 114, 185–208. doi: 10.1016/bs.ctdb.2015.07.026
Wang, Z., Cheng, X., Meng, Q., Wang, P., Shu, B., Hu, Q., et al. (2015). Azadirachtin-induced apoptosis involves lysosomal membrane permeabilization and cathepsin L release in Spodoptera frugiperda Sf9 cells. Int. J. Biochem. Cell Biol. 64, 126–135. doi: 10.1016/j.biocel.2015.03.018
Wu, Y., Wu, Y., Hui, T., Wu, H., Wu, Y., and Wang, W. (2013). IBM1 homologue IBM1 in silkworm, Bombyx mori induces apoptosis upon baculovirus infection. FEBS Lett. 587, 600–606. doi: 10.1016/j.febslet.2013.01.072
Wu, Y. Y., Zhou, T., Wang, Q., Dai, P. L., Xu, S. F., Jia, H. R., et al. (2015). Programmed cell death in the Honey Bee (Apis mellifera) (Hymenoptera: Apidae) worker brain induced by Imidacloprid. J. Econ. Entomol. 108, 1486–1494. doi: 10.1093/jee/tov146
Xu, L., Li, S., Ran, X., Liu, C., Lin, R., and Wang, J. (2016). Apoptotic activity and gene responses in Drosophila melanogaster S2 cells, induced by azadirachtin A. Pest Manag. Sci. 72, 1710–1717. doi: 10.1002/ps.4198
Yang, M., Wang, B., Gao, J., Zhang, Y., Xu, W., and Tao, L. (2017). Spinosad induces programmed cell death involves mitochondrial dysfunction and cytochromeC release in Spodoptera frugiperda Sf9 cells. Chemosphere 169, 155–161. doi: 10.1016/j.chemosphere.2016.11.065
Yi, X., Qi, J., Zhou, X., Hu, M. Y., and Zhong, G. H. (2017). Differential expression of chemosensory-protein genes in midguts in response to diet of Spodoptera litura. Sci. Rep. 7:296. doi: 10.1038/s41598-017-00403-5
Yooboon, T., Pluempanupat, W., Koul, O., and Bulangpoti, V. (2015). Effects of azadirachtin on cuticular proteins of Spodoptera litura (Lepidoptera: Noctuidae) vis-à-vis the modes of application. Commun. Agric. Appl. Biol. Sci. 80, 169–177.
Zhu, X., Jiang, X., Li, A., Sun, Y., Liu, Y., Sun, X., et al. (2017). S-allylmercaptocysteine suppresses the growth of human gastric cancer xenografts through induction of apoptosis and regulation of MAPK and PI3K/Akt signaling pathways. Biochem. Biophys. Res. Commun. 491, 821–826. doi: 10.1016/j.bbrc.2017.06.107
Keywords: Spodoptera litura, azadirachtin, midgut, growth regulation, apoptosis
Citation: Shu B, Zhang J, Cui G, Sun R, Yi X and Zhong G (2018) Azadirachtin Affects the Growth of Spodoptera litura Fabricius by Inducing Apoptosis in Larval Midgut. Front. Physiol. 9:137. doi: 10.3389/fphys.2018.00137
Received: 09 October 2017; Accepted: 12 February 2018;
Published: 27 February 2018.
Edited by:
Senthil-Nathan Sengottayan, Manonmaniam Sundaranar University, IndiaReviewed by:
Jalal Jalali Sendi, University of Gilan, IranDong Wei, College of Plant Protection, Southwest University, China
Copyright © 2018 Shu, Zhang, Cui, Sun, Yi and Zhong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guohua Zhong, guohuazhong@scau.edu.cn