- Research article
- Open Access
DnaJ homolog subfamily A member1 (DnaJ1) is a newly discovered anti-apoptotic protein regulated by azadirachtin in Sf9 cells
BMC Genomics volume 19, Article number: 413 (2018)
Azadirachtin, one of the most promising botanical insecticides, has been widely used for pest control. Azadirachtin induces apoptosis in insect cell lines, including Sf9, SL-1 and BTI-Tn-5B1–4. Mitochondrial and lysosomal pathways are likely involved in the azadirachtin-induced apoptosis, however, detailed molecular mechanisms remain largely undefined.
Azadirachtin-induced apoptosis in Sf9 cells was verified by morphological observation, Hoechst 33258 staining, and a Caspase-3-based analysis. Comparative two-dimensional gel electrophoresis (2-DE) coupled with a linear ion trap quadrupole (LTQ)-MS/MS analysis identified 12 prominent, differentially expressed proteins following azadirachtin treatment. These differentially expressed genes are involved in regulating cytoskeleton development, signal transduction, gene transcription, and cellular metabolism. Knockdown gene expression of a gene encoding a DnaJ homolog enhanced apoptosis induced by azadirachtin in Sf9 cells.
Azadirachtin treatment induces apoptosis in Sf9 cells and affects expression of multiple genes with functions in cytoskeleton development, signal transduction, gene regulation, and cellular metabolisms. Azadirachtin induces apoptosis at least partially by down-regulation of Sf-DnaJ in Sf9 cells.
Azadirachtin, a prototypical botanical tetranortriterpenoid isolated from neem trees (Azadirachta indica, A. Juss), is one of the most potent botanical chemicals and has been used extensively in pest control [1,2,3]. Azadirachtin is effective against more than 550 species of insect pests, including insects from Lepidoptera, Hemiptera, Diptera and Orthoptera. The mode of action of azadirachtin against insects include antifeedant effect and disruption of insect growth and development. On the other hand, azadirachtin has little toxicity to mammals and decays fast in the environment [4,5,6], which makes it a preferred choice for pest management in the field. Azadirachtin is also toxic to cultured insect cells. Inhibition of cell proliferation has been observed in Sf9 cells derived from the ovaries of Spodoptera frugiperda, SL-1 cells derived from Spodoptera litura, BTI-Tn-5B1–4 cells derived from Trichoplusia ni, and C6/36 cells derived from Aedes albopictus [7,8,9,10,11]. Treatments of these cells with 10 to 100 nM azadirachtin result in completely inhibition of cell proliferation [7, 8]. Studies with some of the insect cell lines suggest that apoptosis is the cause of cell death based on observed morphological, physiological, biochemical, and toxicological changes [9,10,11,12].
The high efficacy of azadirachtin against cultured cells and insects has attracted a great deal of attention to reveal the molecular pathways for its mode of action. However, so far most information on molecular mechanisms associated with azadirachtin toxicity has been obtained from cancer cell lines. Apoptotic signaling pathways are activated in cancer cells following azadirachtin treatments, including the caspase-dependent pathway, AIF-mediated pathway, p38 and JNK1/2 pathway, ROS-dependent MAPK pathway and death receptor pathway [13,14,15]. In insect cells, the p53 gene is induced in azadirachtin-treated SL-1 cells, resulting in cell cycle arrest and the induction of apoptosis . Using insect Sf9 cells, our group has previously demonstrated that both mitochondrial and lysosomal pathways are involved in apoptosis after azadirachtin treatments [17, 18]. Specifically, we found that cathepsin L released from lysosome to cytosol was induced in azadirachtin-treated Sf9 cells, resulting in the activation of caspase-3 . Despite significant progress has been made, our knowledge on molecular components and pathways leading to apoptosis in azadirachtin-treated cells remains fragmented.
Comparative proteomic analyses are powerful and effective tools for large-scale identification of proteins involved in a specific biological process. Two-dimensional gel electrophoresis (2-DE) combined with mass spectrometry (MS) has commonly used for proteomics and has been extensively applied to analyze the differentially expressed proteins in identical biological samples that are treated differently [19, 20]. For example, 10 proteins of S. litura (Fabricius) affected by azadirachtin significantly have been identified using 2-DE, and six of them are functionally assigned based on matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF-MS) . Two induced hemolymph proteins with functions in lipid metabolism have also been identified using 2-DE coupled with MS/MS from azadirachtin-treated Ostrinia furnacalis (Lepidoptera: Crambidae) . Twenty-one differentially expressed proteins have been identified using the 2-DE/MS/MS method in azadirachtin-treated Drosophila melanogaster larvae, with results indicating that heat shock protein 23 is the potential target of azadirachtin action .
The objective of this study is to conduct a systematic study on proteins expressed differentially in Sf9 cells after treatments with azadirachtin via comparative proteomic analyses. Twelve most differentially expressed proteins were identified by linear ion trap quadrupole (LTQ)-MS/MS. Among these 12 proteins, Sf-DnaJ1 (DnaJ homolog to subfamily A member 1) was down-regulated significantly in azadirachtin-treated Sf9 cells. Knockdown of Sf-DnaJ1 via RNA interference (RNAi) resulted in increased apoptosis in Sf9 cells after azadirachtin treatments. Our results suggest that Sf-DnaJ1 is a regulator of azadirachtin-induced apoptosis.
Hyclone SFX-insect cell culture medium was purchased from Thermo Scientific (USA), and fetal bovine serum (FBS) was purchased from Gibco (USA). IPG drystrip and IPG buffer were purchased from Amersham Biosciences (Uppsala, Sweden). Azadirachtin (95%) was obtained from Sigma (St.Louis, MO, USA). Other chemicals were domestic products with analytical grade. Rabbit polyclonal antibodies against HSP40, TCTP and GAPDH, respectively, were obtained from BOSTER (Wuhan, China).
Sf9 cells were obtained from the State Key Laboratory for Biocontrol/Institute of Entomology, Sun Yat-Sen University (Guangzhou, China), and were maintained at 27 °C in 25 cm2 culture flasks (Corning, USA) containing 3 mL Hyclone SFX-insect cell culture medium supplemented with 5% fetal bovine serum. The doubling time under optimum conditions was 18–24 h and cells were subcultured every 2 days.
Cell viability assay
Sf9 cells were seeded onto a 96-well plate (5 × 103/well) and incubated for 24 h, then exposed to a series of concentrations of azadirachtin for 24 h and 0.1% DMSO was included as a control. Fifty μL of methylthiazoletetrazolium (MTT) solution was added to each well and cells were incubated in dark for another 4 h. After removing the supernatant, 150 μL DMSO was added and mixed thoroughly with the pipette. Cell viability was measured based on absorbance at 490 nm using a microplate reader (Thermo Scientific, Waltham, MA, USA).
Morphological observation by inverted phase contrast microscopy
Cells seeded in 6-well plates were treated with 0.75 μg/mL azadirachtin and 0.1% DMSO was used as control. Morphological characteristics of cells at 0, 24, 48, and 72 h after treatments were recorded by an inverted phase contrast microscope (Lecia, Japan), respectively.
Hoechst 33258 analysis
Hoechst 33258 is a blue fluorescent dye that could penetrate cell membranes and stain the cell nuclei with blue color. Sf9 cells treated with azadirachtin for 24, 36 and 48 h were stained with 0.5 mL Hoechst 33258 solution for 5 min and then washed with phosphate-buffered saline (PBS) twice for 3 min each time. The stained cells were observed under the fluorescent microscope (Nikon, Japan).
Analysis of caspase-3-like enzymatic activity
Sf9 cells treated with azadirachtin were collected and the caspase-3-like proteolytic activity was measured using a Caspase-3 Colorimetric Assay Kit (KeyGEN BioTECH, Nanjing, China). The cells were then washed with PBS twice and collected by centrifugation at 2000 rpm for 5 min. Total cellular proteins were extracted using a cold lysis buffer on ice for 20–60 min. Protein concentrations were determined following the Bradford approach . The solution containing 150 mg proteins with the Caspase-3 substrate (integrating specific luminescence substrate) was incubated in dark at 37 °C for 4 h. Caspase-3-like enzymatic activity was measured based on the absorbance of samples measured at 405 nm using a microplate reader (Thermo Scientific, USA). The Cacpase-3 inhibitor Z-VAD-FMK was also used with the final concentration of 20 μM.
Preparation of protein samples
The adherent Sf9 cells treated with 0.75 μg/mL azadirachtin for 24 h and control group were washed with PBS twice and then mixed with 1 mL lysis buffer containing 40 mM Tris-base, 7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 2% (v/v) carrier ampholytes pH 3–10 and 65 mM DTT. The homogenates were shaken for 15 min in an ice-water bath and centrifuged at 14000 rpm for 15 min at 4 °C. Protein concentrations of supernatants were determined by the Bradford method.
2-DE, gel staining and image analysis
Before loading for 2-DE, samples were dissolved in 350 μL rehydration buffer containing 7 M urea, 2 M thiourea, 4% (w/v) CHAPS, 2% (v/v) IPG buffer, 20 mM DTT, and a trace of bromophenol blue, then centrifuged at 14000 rpm for 5 min. Total protein extracts from control and treated samples were separated through 2-DE. Two protein samples (140 μg) were loaded onto analytical and preparative gels, respectively. Isoelectric focusing (IEF) was carried out on an IPGphor system (Amersham Biosciences) with pH 4–7 IPG strips (18 cm, linear) according to the manufacturer’s instructions. A total of 60 kVh was applied. Then the IPG strips were equilibrated in 3 mL equilibration buffer twice for 15 min. The first equilibration was performed in a buffer containing 50 mM Tris-HCl (pH 8.8), 6 M urea, 30% (v/v) glycerol, 2% (w/v) SDS, 1% (w/v) DTT, and a trace amount of bromophenol blue. The second equilibration was performed in a buffer modified by 2.5% w/v IAA instead of DTT. The strips were placed on the top of 12.5% SDS-polyacrylamide gels and sealed with 0.5% agarose. Electrophoresis was carried out on a Hoefer SE 600 apparatus (Amersham Biosciences) at 20 °C with the current of 15 mA/gel for 40 min, and then 45 mA/gel for 6 h. The protein spots in gels were visualized by staining with silver nitrate . At least three replicates were performed for each sample. Images of each gel were acquired using Lab-Scan version 3.0 software (GEHealthcare) on an Image-Scanner. Images were analyzed by ImageMaster 2-DE platinum version 5.0 software. The intensity of the protein spots was calculated with PDQuest 8.0 software.
Protein digestion, LTQ-MS/MS and database searching
In order to locate protein spots with different intensity in gels, Coomassie Brilliant Blue G-250 was used to stain the gels for mass spectrometric analysis. The gels were fixed with a buffer containing 40% ethanol and 10% glacial acetic acid for 1 h, washed with double distilled water three times, stained with Coomassie Brilliant Blue G-250 staining solution overnight, and decolorized in a destaining solution for at least 4 h. The significantly altered protein spots were located by comparing gels with control and treated samples side by side. The identified protein spots were cut out from the gel, further destained with 30 mM potassium ferricyanide/100 mM sodium thiosulfate (1:1, v/v) for 20 min, and washed in Milli-Q water until the gels were completely destained. The spots were kept in 0.2 M NH4HCO3 for 20 min and then lyophilized. Each spot was digested in 12.5 ng/mL trypsin with 0.1 M NH4HCO3 overnight. The peptides were extracted with 50% Acetonitrile, and 0.1% TFA three times.
Separation and identification of the digested proteins were conducted on a Finnigan LTQ mass spectrometer (ThermoQuest, San Jose, CA, USA) coupled with a Surveyor HPLC system (ThermoQuest, San Jose, CA, USA). A Microcore RP column (C18 0.15 mm × 120 mm; ThermoHypersil, San Jose, CA, USA) was used to separate the protein digests. Solvent A was 0.1% (v/v) formic acid, and solvent B was 0.1% (v/v) formic acid in 100% (v/v) ACN. The gradient was held at 2% solvent B for 15 min, and increased linearly to 98% solvent B for 90 min. The peptides were eluted from the C18 microcapillary column at a flow rate of 150 μL/min and then electrosprayed directly into an LCQ-Deca mass spectrometer with the application of spray voltage of 3.2 kV and capillary temperature at 200 °C. The full scan was ranged from M/Z 400 to 2000. Protein identification based on MS/MS data was performed with SEQUEST software (University of Washington, licensed to Thermo Finnigan) based on the database of Swiss Port. The species for sequence search is Lepidoptera. Protein identification results were filtered with a stringent filter condition of Xcorr (1 + ≥ 1.9, 2 + ≥ 2.2, 3 + ≥ 3.75) and DelCn (≥ 0.1).
Identification and sequencing of sf-DnaJ1 cDNA
To obtain a full length Sf-DnaJ1 cDNA, total RNA was isolated from Sf9 cells with an E.Z.N.A.™ Total RNA Kit II (OMEGA, USA) according to the manufacturer’s instructions. First strand cDNAs were synthesized using a PrimeScript® 1st Strand cDNA Synthesis Kit (TaKaRa) according to the provided protocol. cDNA of Sf-DnaJ1 was amplified by PCR with the degenerate primers Sf-DnaJ1-F and Sf-DnaJ1-R (Additional file 1: Table S1) and 50 μL reaction mixture contained 0.5 μL template, 1 μL of each 10 mM primer, 0.5 μL Taq DNA polymerase (TIANGEN, Beijing, China), 4 μL of 2.5 mM dNTP mixture (TIANGEN, Beijing, China), and 5 μL 10× Taq Buffer. The PCR program was performed with 32 cycles of 30 s at 94 °C, 30 s at 55 °C and 30 s at 72 °C. The 5’-RACE (SMART RACE, Clontech) and 3’-RACE (TaKaRa) methods were used to fulfill the full-length cDNA of Sf-DnaJ1. To ensure the 5′ and 3′ fragments were cloned from the same gene, specific primers were designed and PCR was used to amplify the coding region of the transcript encoding Sf-DnaJ1.
Quantitative real time PCR
In order to confirm the expression profiles of six identified proteins, quantitative real-time PCR (qRT-PCR) was performed. Total RNA was extracted from Sf9 cells treated with azadirachtin for 24 h and control cells using a Total RNA Kit II (OMEGA, USA). The cDNA for qRT-PCR was synthesized using a PrimeScript™ RT reagent Kit (TaKaRa, Japan), which has a gDNA Eraser to eliminate DNA contamination. qRT-PCR was performed on CFX Connect™ Real-Time System (Bio-Rad, USA) using SsoAdvanced™ SYBR® Green Supermix (Bio-Rad, USA). The PCR was carried out as follows: 95 °C for 3 min for denaturation, 40 cycles of 95 °C for 10 s, 60 °C for 10 s, 72 °C for 30 s, and a dissociation step at the end. GAPDH was used as a reference for normalization. Relative expression levels were calculated by the 2−ΔΔCT method. Primers used in the experiments are listed in Additional file 1: Table S1.
Western blot assays
Cells were collected and washed with PBS. Total cellular proteins were extracted using the CytoBuster™ Protein Extraction Reagent (Novagen, USA) according to the manufacturer’s protocol. Protein concentrations were determined by the Bradford method. Equal amounts of proteins from different samples were separated on a 12% SDS-PAGE gel. Proteins in the gel were then transferred to a polyvinylidene difluoride membrane (PVDF, Millipore, USA). The membrane was washed with TBS for 3 times, incubated with TBS supplemented with 5% fat-free milk at 4 °C overnight, and incubated with HSP 40 antibody, TCTP antibody or GAPDH antibody at room temperature for 2 h. Subsequently, the membrane was washed and incubated with the peroxidase-conjugated secondary antibody at room temperature for more than 2 h. The protein bands were detected by the enhanced chemiluminescence western blot kit (CW0049, CWBIO, Beijing, China) and detected by exposure to X-ray film a dark room.
Double-stranded RNA synthesis and transfection
dsRNA against the Sf-DnaJ1 transcript was synthesized using a T7 RiboMAX™ Express RNAi System (Promega, USA) according to the provided protocol. Template DNA was amplified by PCR with the RNAi primers listed in Additional file 1: Table S1. The dsRNA against egfp, which used as the negative control, was similarly synthesized by the template pEGFP-C and primers in Additional file 1: Table S1. The size and integrity of dsRNAs were checked by agarose gel electrophoresis.
Transfection of Sf9 cells was performed based on the Lipofectin transfection method. Monolayer cultures of Sf9 cells were prepared in 35-mm cell culture dishes (Corning, USA). Transfection was carried out by incubation with 2 mL Hyclone SFX-insect cell culture medium (without FBS) containing 5 μg dsRNAs overnight at 27 °C, followed by incubation in 10 μL lipofectamine 2000 (Invitrogen) for 6 h. The medium was then replaced with a medium containing FBS. After 24 h treatment, the RNAi efficiency was examined based on qRT-PCR and western blot results.
Annexin V-FITC/propidium iodide double-staining and flow cytometry
Anchorage-dependent Sf9 cells treated with 0.75 μg/mL azadirachtin for 24 h were collected by centrifugation at 2000 rpm for 5 min at 4 °C. The cells were then resuspended and washed twice with PBS. The cells were then fixed in 500 μL binding buffer. Prior to cytometry analysis, 5 μL Annexin V-FITC and 5 μL PI were added to the fixed cells which were then incubated in dark for 15 min at room temperature. The cells were analyzed through a flow cytometry with an Ar laser with excitation and emission wavelengths 488 nm and 530 nm, respectively. At least 2.0 × 104 cells were counted in each assay. The Sf9 cells with Sf-DnaJ1 knocked down and GFP control cells were treated with 0.75 μg/mL azadirachtin for 24 h and used for analyses.
Each treatment had three replicates and data were expressed as the mean values ± SEM. One-way ANOVA followed with Duncan’s new multiple range test (DMRT) and student’s t test were conducted during statistical analyses (P < 0.05).
Azadirachtin inhibited cell viability and proliferation
As shown in Fig. 1a, the inhibition rates of Sf9 cell proliferation were 14.6 ± 2.21%, 23.1 ± 3.32%, 28.7 ± 2.44%, 36.1 ± 1.56%, 43.8 ± 2.25% and 45.5 ± 1.35%, respectively, after cells treated with azadirachtin at the concentrations of 2, 5, 10, 20, 40, 50 μg/mL for 24 h. These results revealed that azadirachtin had a strong inhibition effect on cells proliferation and decreased cell viability in a dose-dependent manner.
Morphological changes associated with apoptosis
Our previous study showed that a low concentration of azadirachtin induced apoptosis of Sf9 cells . As shown in Fig. 1b, the morphological changes of Sf9 cells treated with 0.75 μg/mL azadirachtin could be observed clearly under the inverted phase contrast microscopy. Typical morphological characteristics of apoptosis were displayed in Sf9 cells treated with azadirachtin for 24 h, including cell shrinkage, increased gaps, membrane blebbing and apoptotic bodies. After treatment for 48 h, gaps of cells increased and apoptotic bodies occurred widely. The number of viable cells and apoptotic bodies reduced greatly after treatment for 72 h, and few cells kept the normal morphological appearance. These results suggest that typical morphological characteristics of apoptosis are induced in cells treated with 0.75 μg/mL azadirachtin in Sf9 cells.
Azadirachtin induced nuclear condensation
As shown in Fig. 1c, Hoechst 33258 staining revealed nuclear changes induced by azadirachtin in Sf9 cells. The nuclei of control cells exhibited uniform sizes and were stained homogeneous. In comparison, nuclei of azadirachtin-treated cells exhibited deeper blue staining with nuclear and chromatin condensation after treatment for 24, 36 and 48 h, respectively. The number of live cells decreased gradually and nuclear condensation became more obvious with longer treatment time.
Azadirachtin increased caspase-3-like proteolytic activity
In order to examine if caspases were activated in Sf9 cells treated with azadirachtin, caspase-3-like proteolytic activity was determined. Compared with that in control cells, the caspase-3-like proteolytic activity increased 2.66, 7.08, 9.26 and 9.91 fold in cells treated with azadirachtin for 12, 24, 36 and 48 h, respectively (Fig. 1d). Apoptosis induced by azadirachtin was inhibited completely by Z-VAD-FMK and the caspase-3-like proteolytic activity in cells treated with both azadirachtin and Z-VAD-FMK was similar to normal cells (Additional file 2: Figure S1). These results indicate that azadirachtin induce apoptosis via activating caspase-3 activity in a time-dependent manner.
Identification of differential expressed proteins
Comparative proteomic analyses were performed to identify azadirachtin-responsive proteins. Approximately 800 protein spots were detected on a 2-D gel using the Image Master 5.0 software. The criterion for proteins with significant changes in abundance was at least a 1.5-fold increase or decrease between control cells and cells treated with azadirachtin for 24 h. Thirteen protein spots indicated by arrows satisfied this criterion and were considered differentially expressed proteins (Fig. 2). Relative intensity of these proteins was showed in Additional file 3: Figure S2. Abundance of proteins corresponding to spots D1 to D11 decreased in treated cells whilst abundance of proteins corresponding to spots U1 and U2 increased in treated cells (Fig. 3).
The identified azadirachtin-responsive proteins are listed in Table 1 and the matched peptide sequences are listed in Additional file 4: Table S2. Proteins up-regulated by azadirachtin were ribosomal protein L9 (Sf-RpL9, U1) and abnormal wing disc-like protein (Sf-AWD, U2). Proteins down-regulated were DnaJ homolog subfamily A member1 (Sf-DnaJ1, D1), ribosomal protein SA (Sf-RpSA, D2), actin (D3-D5), beta-tubulin (D6), proteasomezeta subunit (Sf-PS, D7), P27BBP/eIF6-like (Sf-P27BBP/eIF6, D8), translationally-controlled tumor protein homolog (Sf-TCTP, D9), Proteasome subunit alpha type 6-A (Sf-PMSA6, D11), and an unknown protein (D10).
Characterization of sf-DnaJ1transcript
A full length cDNA encoding Sf-DnaJ1 was cloned and sequenced as described in Materials and Methods. The nucleotide sequence was deposited in GenBank with the accession number KF562156. The Sf-DnaJ1-encoding cDNA is 1757 bp long with 148 bp at 5′- and 397 bp at the 3′-untranslated regions. The predicted protein has 404 amino acid residues with calculated molecular mass 45.46 kDa. The predicted Sf-DnaJ1 share 89, 77 and 67% sequence identity with DnaJ family proteins from Bombyx mori, Tribolium castaneum and Aedes aegypti, respectively. The phylogenetic relationship of Sf-DnaJ1 with paralogs from other species is shown in Additional file 5: Figure S3.
Validation of differentially expressed proteins by qRT-PCR and western blot
To examine the expression of the six genes affected by azadirachtin at transcriptional level, qRT-PCR was conducted. Consistent with the proteomic analysis, the levels of the transcripts encoding Sf-DnaJ1, Sf-PS, Sf-P27BBP/eIF6, Sf-TCTP and Sf-PMSA6 decreased, whereas the level of the transcript encoding Sf-AWD increased in azadirachtin-treated Sf9 cells. Specifically, the level of transcripts encoding Sf-DnaJ1, Sf-PS, Sf-P27BBP/eIF6, Sf-TCTP and Sf-PMSA6 decreased by 39.2, 43, 47, 57.1and 23.7%, respectively (Fig. 4a), while the level of the transcript encoding Sf-AWD increased 107.6%. Western blot analysis revealed that DnaJ1 and TCTP decreased in Sf9 cells treated with azadirachtin for 24 h.
Knockdown of sf-DnaJ enhanced apoptosis induced by azadirachtin
To further confirm the role of Sf-DnaJ1 in azadirachtin-induced apoptosis, dsRNA against Sf-DnaJ1 transcripts was used for gene silencing. The transcript and protein levels of Sf-DnaJ1 were reduced significantly (98.2%) in dsRNA-treated cells compared to controls (Fig. 5a and b). In addition, the cells treated with ds-DnaJ1 could not induced apoptosis and the apoptotic rate had no significant difference with normal cells (Additional file 6: Figure S4). However, silencing the ds-DnaJ1-encoding gene resulted in the enhancing of azadirachtin’s effect on apoptosis (Fig. 5c, d).
Azadirachtin has been proven to be an effective insecticidal ingredient for pest management and has the ability to induce apoptosis in cultured cell lines . In this study, Sf9 cells treated with 0.75 μg/mL azadirachtin exhibited morphological changes typical of apoptosis, as have been previously observed in other insect cells in vitro [11, 17]. Apoptosis is a form of physiological cell death with important roles in various biological processes . Although it has been well documented that apoptosis is associated with azadirachtin treatments, the information on the molecules and pathways affected by azadirachtin remained fragmented. In order to identify key regulators and uncover the molecular mechanism of apoptosis induced by azadirachtin, proteomic analyses including 2-DE coupled with MS-based protein identification was carried out. Twelve proteins were identified with significant changes after azadirachtin treatment (Fig. 2).
One of the down-regulated proteins was identified as DnaJ homolog subfamily A member 1 (alternative name: heat shock 40 kDa protein 4, hsp40–4). Heat shock proteins (HSPs) have been recognized as common molecular chaperones induced under various stress conditions, including heat shock, heavy metals, exposure to radiation, ethanol, oxidative stress and so on . Previous studies have provided substantial evidence that HSPs may regulate apoptosis. Hsp27 and Hsp70 have effects as anti-apoptotic proteins [29, 30]. Hsp60 and Hsp10 can enhance the proteolytic maturation of caspase-3 . Inhibition of Hsp70 and Hsp40 by HSP inhibitor KNK437 results in enhanced PS-341-induced cell death . Hsp40–Hsp70 pair plays a counter role through interacting with Bax, resulting in the inhibition of its translocation to mitochondria in NO-induced apoptosis . HSP 90 has been identified as a potential target of azadirachtin effect through reversing docking . Hsp40 has been characterized as a co-chaperone involved in regulation of Hsp70 chaperone activity, but it remains unclear whether it can regulate apoptosis independently [35, 36]. In this study, we provided two pieces of evidence that, for the first time, demonstrated that a Hsp40 member (Sf-DnaJ1) regulates apoptosis. First, we found that the protein abundance of Sf-DnaJ1 decreased significantly in Sf9 cells treated with azadirachtin, indicating that Sf-DnaJ1 was one of the regulators of azadirachtin action. Second, silencing Sf-DnaJ1-encoding gene enhanced the capability of azadirachtin on inducing apoptosis. Our results indicated that Hsp40 can be an independent apoptosis regulator. Further studies are required to reveal the molecular pathway between azadirachtin and Sf-DnaJ1, and downstream steps.
Previous studies have shown that azadirachtin could induce the depolymerization of actin, leading to cell arrest and subsequently apoptosis in a caspase-independent manner [37, 38]. Our results provide further support to this observation. Compared with untreated control cells, the abundance of actins (three protein spots corresponding to D3, D4 and D5) and β-tubulin (D6) decreased dramatically in Sf9 cells treated with azadirachtin, suggesting that azadirachtin exerted its effects by inhibiting actin and tubulin to remodel apoptotic cells.
Proteasomes are ubiquitous and abundant multi-catalytic enzyme complexes in charge of the degradation of most intracellular proteins . A previous study has showed that a ubiquitin-proteasome complex is involved in the apoptotic process . Proteasomes can stabilize pro-apoptotic proteins including p53 and Bax in various cell types, but they also could suppress apoptosis induced by some stimuli [41,42,43,44], suggesting that the roles of proteasomes in apoptosis are complex and can go to either direction. In our results, the abundance of proteasome subunits zeta (D7) and alpha type 6-A (D11) decreased in Sf9 cells after azadirachtin treatment, indicating that azadirachtin could be a new proteasome inhibitor.
Translationally-controlled tumor protein homolog (TCTP) is generally regarded as a highly conserved protein with multiple functions. It is widely expressed in various tissues or cell types and plays an important role in maintaining survival of a variety of cell types . Several lines of evidence indicate that TCTP has a role as an anti-apoptotic protein in cultured cells and it protects cells from apoptotic cell death by inserting into the mitochondrial membrane and inhibiting the dimerization of Bax . TCTP probably interacts selectively with actin and microtubule cytoskeleton to regulate cell shape during interphase and mitosis . TCTP inhibits apoptosis by controlling the stability of tumor suppressor p53, which in turns represses the transcription of the TCTP-encoding gene [48, 49]. In this study, the abundance of TCTP (D9) decreased in Sf9 cells treated with azadirachtin, indicating that azadirachtin could regulate the proliferation and apoptosis of Sf9 cells via inhibiting the gene expression of TCTP.
The abundance of three ribosome proteins was affected by azadirachtin. Azadirachtin reduced the abundance of ribosomal protein SA (D2) and P27BBP/eIF6-like (D8), but elevated the abundance of the ribosomal protein L9 (U1). Ribosomal proteins are generally involved in protein synthesis. The impact of azadirachtin on the abundance of ribosomal proteins could suggest that this chemical plays a role in protein biosynthesis.
Several proteins including Bcl-2, NF-κB, P53 and PI3K have been previously suggested to be targets of azadirachtin [37,38,39]. However, these proteins were not identified in our proteomic analyses. One possible explanation is that the abundance of these regulatory proteins was less affected by azadirachtin, and therefore, were not picked up in our analysis since we did focused on the proteins with most differences. Alternatively, the pH range for our gel strip was pH 4–7, which might not be optimal for separation of these proteins. More studies with different pH ranges should be used to identify other differentially affected proteins after azadirachtin treatments.
In summary, the present study confirmed that azadirachtin induced apoptosis in Sf9 cells. Azadirachtin-responsive proteins in Sf9 cells were analyzed and 12 differentially expressed protein spots were identified through 2-DE and LTQ-MS/MS analyses. Transcription of the genes encoding Sf-DnaJ1, Sf-PS, Sf-P27BBP/eIF6, Sf-TCTP, Sf-PMSA6 and Sf-AWD was also affected by azadirachtin. In addition, changes in protein levels of Sf-DnaJ1 and Sf-TCTP were confirmed by western blot. Knockdown of the gene encoding Sf-DnaJ1 by RNAi resulted in increased apoptosis induced by azadirachtin, suggesting that Sf-DnaJ1 played an anti-apoptotic role in azadirachtin-induced apoptosis of Sf9 cells. Overall, our results indicated that azadirachtin triggered apoptosis by down-regulating Sf-DnaJ1-encoding gene. Further experiments are needed to obtain more evidence for a more detailed picture on the molecular processes of azadirachtin toxicity on cultured cells and field insects.
Linear ion trap quadrupole-mass spectromtry
Two-dimensional gel electrophoresis
Apoptosis inducing factor
DnaJ homolog subfamily A member1
c-Jun N-terminal kinase
Reactive oxygen species
Mitogen-activated protein kinase
Matrix-assisted laser desorption/ionization time-of-flight
Fetal bovine serum
3-(4, 5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide
Translationally-controlled tumor protein homolog
Ribosomal protein L9
Abnormal wing disc-like protein Sf-RpSA: ribosomal protein SA
Proteasome subunit alpha type 6-A
Mordue AJ, Blackwell A. Azadirachtin: an update. J Insect Physiol. 1993;39:903–24.
Isman MB. Botanical insecticides, deterrents, and repellents in modern agriculture and an increasingly regulated world. Annu Rev Entomol. 2006;51:45–66.
Morgan ED. Azadirachtin, a scientific gold mine. Bioorgan Med Chem. 2009;17:4096–105.
Thomas LD, Murray BI. Insect growth regulating effects of neem extracts and Azadirachtin on aphids. Entomol Exp Appl. 1994;72:77–84.
Shi P, Huang Z, Chen G, Zhou L, Tan X. Effects of Azadirachtin on six inorganic cation distributions in Ostrinia furnacalis (G.). Biol Trace Elem Res. 2006;113(1):105–12.
Khosravi R, Sendi JJ. 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 Protection Res. 2013;53(3):238–47.
Rembold H, Annadurai RS. Azadirachtin inhibits proliferation SF-9 cells in monolayer culture. Z Naturforsch C. 1993;48:495–9.
Salehzadeh A, Jabbar A, Jennens L, Ley SV, Annadurai RS, Adams R, Strang RH. The effects of phytochemical pesticides on the growth of cultured invertebrate and vertebrate cells. Pest Manag Sci. 2002;58:268–76.
Zhong G, Shui K, Huang J, Jia JW, Hu M. Induction of apoptosis by botanical components in Spodoptera litura cultured cell line. Acta Entomol Sin. 2008;51(4):449–53.
Zhong G, Shui K, Lv C, Jia JW, Ren T, Hu M. Induction of apoptosis by azadirachtin, a botanical insecticidal component, in Spodoptera litura cultured cell line SL-1. Acta Entomol Sin. 2008;51(6):618–27.
Huang XY, Li OW, Xu HH. Induction of programmed death and cytoskeletal damage on Trichoplusia ni BTI-Tn-5B1-4 cells by azadirachtin. Pestic Biochem Phys. 2010;98:289–95.
Shu B, Wang W, Hu Q, Huang J, Hu M, Zhong G. A comprehensive study on apoptosis induction by azadirachtin in Spodoptera frugiperda cultured cell line Sf9. Arch Insect Biochem. 2015;89(3):153–68.
Srivastava P. Neem oil limonoids induces p53-independent apoptosis and autophagy. Carcinogenesis. 2012;33(11):2199–207.
Babykutty S. Nimbolide retards tumor cell migration, invasion, and angiogenesis by downregulating MMP-2/9 expression via inhibiting ERK1/2 and reducing DNA-binding activity of NF-kB in colon cancer cells. Mol Carcinogen. 2012;51(6):475–90.
Gupta SC, Reuter S, Phromnoi K, Park B, Hema PS, Nair M, Aggarwal BB. Nimbolide sensitizes human colon cancer cells to TRAIL through reactive oxygen species- and ERK-dependent up-regulation of death receptors, p53, and Bax. J Biol Chem. 2011;286(2):1134–46.
Huang JF, Shui KJ, Li HY, Hu M, Zhong G. 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. 2011;99:16–24.
Huang J, Lv C, Hu M, Zhong G. The mitochondria-mediate apoptosis of lepidopteran cells induced by Azadirachtin. PLoS One. 2013;8(3):e58499. doi: 10.1371.
Wang Z, Cheng X, Meng Q, Wang P, Shu B, Hu Q, Hu M, Zhong G. Azadirachtin-induced apoptosis involves lysosomal membrane permeabilization and cathepsin L release in Spodoptera frugiperda Sf9 cells. Int J Biochem Cell Bio. 2015;64:126–35.
Gorg A, Weiss W, Dunn MJ. Current two-dimensional electrophoresis technology for proteomics. Proteomics. 2009;4:3665–85.
Wang YQ, Zhang LY, Lai D, Xu HH. The nematicidal and proteomic effects of Huanong AVM (analog of avermectin) on the pine-wilt nematode, Bursaphelenchus xylophilus. Pestic Biochem Phys. 2011;98:224–30.
Huang ZW, Shi P, Dai JQ, Du JW. Protein metabolism in Spodoptera litura (F.) is influenced by the botanical insecticide azadirachtin. Pestic Biochem Phys. 2004;80:85–93.
Huang ZW, Shi P, Chen GC, Du JW. Effects of Azadirachtin on hemolymph protein expression in Ostrinia furnacalis (Lepidoptera: Crambidae). Ann Entomol Soc Am. 2007;100(2):245–50.
Wang H, Lai D, Yuan M, Xu HH. Growth inhibition and differences in protein profiles in azadirachtin-treated Drosophila melanogaster larvae. Electrophoresis. 2014;35(8):1122–9.
Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principles of protein-dye binding. Anal Biochem. 1976;72:248–54.
Jungblut PR, Seifert R. Analysis by high-resolution two-dimensional electrophoresis of differentiation-dependent alterations in cytosolic protein pattern of HL-60 leukemic cells. J Biochem Bioph Meth. 1990;21:47–58.
Ahmad S, Ansari MS, Moraiet MA. Demographic changes in Helicoverpa armigera after exposure to neemazal (1% EC azadirachtin). Crop Prot. 2013;50:30–6.
Kerr JFR, Wyllie AH, Curie AR. Apoptosis: a basic biological phenomenon with wide-ranging implications in tissue kinetics. Brit J Cancer. 1972;26:239–57.
Qian SB, McDonough H, Boellmann F, Cyr DM, Patterson C. CHIP-mediated stress recovery by sequential ubiquitination of substrates and Hsp70. Nature. 2006;440:551–5.
Pandey P, Farber R, Nakazawa A, Kumar S, Bharti A, Nalin C, Weichselbaum R, Kufe D, Kharbanda S. Hsp27 functions as a negative regulator of cytochrome c-dependent activation of procaspase-3. Oncogene. 2000;19:1975–81.
Bienemann AS, Lee YB, Howarth J, Uney JB. Hsp70 suppresses apoptosis in sympathetic neurones by preventing the activation of c-Jun. J Neurochem. 2008;104:271–8.
Samali A, Cai J, Zhivotovsky B, Jones DP, Orrenius S. Presence of a pre-apoptotic complex of pro-caspase-3, Hsp60 and Hsp10 in the mitochondrial fraction of jurkat cells. EMBO J. 1999;18:2040–8.
Liu Y, Zheng T, Zhao S, Liu H, Han D, Zhen Y, Xu D, Wang Y, Yang H, Zheng G, Wang C, Wu J, Ye Y. Inhibition of heat shock protein response enhances PS-341-mediated glioma cell death. Ann Surg Oncol. 2012;19:S421–9.
Gotoh T, Terada K, Oyadomari S, Mori M. Hsp70-DnaJ chaperone pair prevents nitric oxide- and CHOP-induced apoptosis by inhibiting translocation of Bax to mitochondria. Cell Death Differ. 2004;11:390–402.
Chitta K, Paulus A, Caulfield TR, Akhtar S, Blake MK, Ailawadhi S, Knight J, Heckman MG, Pinkerton A, Chanan-Khan A. Nimbolide targets BCL2 and induces apoptosis in preclinical models of Waldenströms macroglobulinemia. Blood Cancer J. 2014;4:e260.
Borges JC, Fischer H, Craievich AF, Ramos CH. Low resolution structural study of two human HSP40 chaperones in solution. DJA1 from subfamily a and DJB4 from subfamily B have different quaternary structures. J Biol Chem. 2005;280:13671–81.
Qiu XB, Shao YM, Miao S, Wang L. The diversity of the DnaJ/Hsp40 family, the crucial partners for Hsp70 chaperones. Cell Mol Life Sci. 2006;63:2560–70.
Anuradha A, Annadurai RS, Shashidhara LS. Actin cytoskeleton as a putative target of the neem limonoid Azadirachtin A. Insect Biochem Molec. 2007;37:627–34.
Pravin Kumar R, Manoj MN, Kush A, Annadurai RS. In silico approach of azadirachtin binding with actins. Insect Biochem Molec. 2007;37:635–40.
Rock KL, Gramm C, Rothstein L, Clark K, Stein R, Dick L, Hwang D, Goldberg AL. Inhibitors of the proteasome block the degradation of most cell proteins and the generation of peptides presented on MHC class I molecules. Cell. 1994;78:761–71.
Strous GJ, Govers R. The ubiquitin-proteasome system and endocytosis. J Cell Sci. 1999;112:1417–23.
Wójcik C. Regulation of apoptosis by the ubiquitin and proteasome pathway. J Cell Mol Med. 2002;6(1):25–48.
Nawrocki ST, Carew JS, Pino MS, Highshaw RA Jr, Dunner K, Huang P, Abbruzzese JL, McConkey DJ. Bortezomib sensitizes pancreatic cancer cells to endoplasmic reticulum stress-mediated apoptosis. Cancer Res. 2005;65:11658–66.
Perry DK, Burns JM, Pollinger HS, Amiot BP, Gloor JM, Gores GJ, Stegall MD. Proteasome inhibition causes apoptosis of normal human plasma cells preventing alloantibody production. Am J Transplant. 2009;9:201–9.
Orlowski RZ, Kuhn DJ. Proteasome inhibitors in Cancer therapy: lessons from the first decade. Clin Cancer Res. 2008;14:1649–57.
Thaw P, Baxter NJ, Hounslow AM, Price C, Waltho JP, Craven CJ. Structure of TCTP reveals unexpected relationship with guanine nucleotide-free chaperones. Nat Struct Biol. 2001;8:701–4.
Susini L, Besse L, Duflaut D, Lespagnol A, Beekman C, Fiucci G, Atkinson AR, Busso D, Poussin P, Marine JC, Martinou JC, Cavarelli J, Moras D, Amson R, Telerman A. TCTP protects from apoptotic cell death by antagonizing bax function. Cell Death Differ. 2008;15:1211–20.
Bazile F, Pascal A, Arnal I, Le Clainche C, Chesnel F, Kubiak JZ. Complex relationship between TCTP, microtubules and actin microfilaments regulates cell shape in normal and cancer cells. Carcinogenesis. 2010;30:555–65.
Rho SB, Lee JH, Park MS, Byun HJ, Kang S, Seo SS, Kim JY, Park SY. Anti-apoptotic protein TCTP controls the stability of the tumor suppressor p53. FEBS Lett. 2010;585:29–35.
Amson R, Pece S, Lespangnol A, Vyas R, Mazzarol G, Tosoni D, Colaluca I, Viale G, Rodrigues-Ferreira S, Wynendaele J, Chaloin O, Hoebeke J, Marine JC, Di Fiore PP, Telerman A. Reciprocal repression between P53 and TCTP. Nat Med. 2012;18(1):91–9.
The authors sincerely thank Prof. Ming-Shun Chen of Kansas State University for revising this manuscript.
The work and publication costs were supported by grants from National Nature Science Foundation of China (No. 31071713).
Availability of data and materials
All data generated or analyzed during this study are available from this published article and its Additional files.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Details of the primer pairs used for genes cloning, RT-qPCR, RNAi. (DOCX 14 kb)
Figure S1. Z-VAD-FMK inhibited the apoptosis induced by azadirachtin in Sf9 cells. A: Morphological changes induced by different treatments in Sf9 cells. 1, 2, 3 and 4 were represented as normal cells, cells treated with Z-VAD-FMK, cells treated with azadirachtin and cells treated with azadirachtin and Z-VAD-FMK, respectively. B: Caspase-3 like activity induced by different treatments in Sf9 cells. (TIF 3988 kb)
Figure S2. Quantitative analysis of the azadirachtin-responsive proteins in Sf9 cells. The data are expressed as arithmetic mean ± SEM of protein intensity on gels from three independent experiments. Statistical analysis was carried out using the SPSS software and different letters above bars indicate significant differences between different treatments at the same time (P < 0.05) by ANOVA followed by DMRT. (TIF 792 kb)
Table S2. The matched peptide sequences of the identified proteins. (DOCX 14 kb)
Figure S3. The coding region sequence, deduced amino acid sequences and phylogenetic analysis of Sf-DnaJ1. A: The coding region sequence and the deduced amino acid sequence of Sf-DnaJ1. B: Phylogenetic analysis of selected DnaJ1. The sequences participate in the phylogenetic tree are: Linepithema humile (XP_012222446.1); Harpegnathos saltator (XP_011150153.1); Pogonomyrmex barbatus (XP_011644113.1); Apis mellifera (XP_006566003.1); Melipona quadrifasciata (KOX75023.1); Megachile rotundata (XP_003699212.1); Athalia rosae (XP_012260458.1); Orussus abietinus (XP_012274281.1); Nasonia vitripennis (XP_008205330.1); Microplitis demolitor (XP_008552261.1); Tribolium castaneum (XP_971446.1); Zootermopsis nevadensis (KDR22500.1); Plutella xylostella (XP_011557028.1); Bombyx mori (NP_001040292.1); Amyelois transitella (XP_013190352.1); Papilio xuthus (XP_013165050.1); Culex quinquefasciatus (XP_001844792.1); Aedes aegypti (ABF18277.1). (PNG 695 kb)
Figure S4. Silence of DnaJ1 didn’t induce apoptosis in Sf9 cells. Fig A-B. A means control, B shows morphological characteristics of Sf9 cells with ds-DnaJ1 treatment for 24 h. Fig C-D. C means control, D shows the apoptosis of Sf9 cells with ds-DnaJ1 treatment for 24 h, ten thousand cells were counted for each sample. Fig E: Apoptotic rate of Sf9 cells with different treatments. The data represent the mean values± S.E.M of three independent experiments. The apoptotic rate of cells with dsDnaJ1 treatment had no significant difference with normal cells. (TIF 1126 kb)
About this article
Cite this article
Shu, B., Jia, J., Zhang, J. et al. DnaJ homolog subfamily A member1 (DnaJ1) is a newly discovered anti-apoptotic protein regulated by azadirachtin in Sf9 cells. BMC Genomics 19, 413 (2018). https://doi.org/10.1186/s12864-018-4801-z