The regulation of inhibitor of apoptosis proteins (IAPs) during the apoptosis of Cotesia chilonis

Inhibitor of apoptosis proteins (IAPs) are crucial components of apoptosis that perform vital roles in the regulation of caspase activity in organisms. In this study, two IAPs genes were identified from Cotesia chilonis, the dominant parasitic wasp of Chilo suppressalis. CcIAP1 gene is a typical IAP and contains two BIR domains and a RING domain, whereas CcIAP gene is an atypical IAP1 only containing two BIR domains. Phylogenetic analysis indicated that CcIAP1 and CcIAP were grouped with other Hymenopteran IAPs and IAP1 in C. suppressalis. Real-time quantitative PCR revealed that CcIAP1 and CcIAP genes were both highly induced at −6°C and 30°C, and expression was highest at the third instar stage. The expression of CcIAP1 and CcIAP genes were significantly induced during parasitism of C. suppressalis, and the 7-d time point resulted in the highest expression levels for both genes, in which was an advanced stage of larval development of C. chilonis. RNAi experiments showed that CcIAP1 gene was the key IAP in the regulation of apoptosis of C. chilonis and its host. In conclusion, CcIAP1 and CcIAP correlate with the development of C. chilonis and their responses to temperature stress.

1 Introduction

Apoptosis is a basic biological process necessary for the normal development of organisms and the maintenance of the internal environment of tissues (Vaux and Korsmeyer, 1999). This term was first proposed by Kerr, Wyllie and Currie in 1972, and this unique mode of cell death is triggered by internal and external environmental changes and can remove senescent, redundant and damaged cells under the regulation of genes (Kerr et al., 1972). As an important arthropod immune response, apoptosis is an active and orderly process of cell death, which can prevent the proliferation of virus and is regulated by a series of proteins. Among them, the inhibitor of apoptosis protein family (IAPs) is critical in the regulation of apoptosis in arthropods (Everett and McFadden, 1999; Vandergaast et al., 2015). Furthermore, previous studies have shown that apoptosis was a generally regulated cell death program, involving the dynamic balance between the activities of pro- and anti-apoptotic factors (Hamajima et al., 2016).

Inhibitor of apoptosis proteins (IAPs) are the crucial inhibitors in the programmed cell death cascade and are named for their anti-apoptotic activity (Danial and Korsmeyer, 2004). The earliest IAP gene was identified when they screened P³⁵ homologous genes from baculovirus in 1993 (Crook et al., 1993). Up to now, six types of IAP genes have been identified in baculoviruses: IAP1, IAP2, IAP3, IAP4, IAP5, and the newly described IAP6 (Luque et al., 2001; Ikeda et al., 2013; Clem, 2015; Chen et al., 2020). In addition, subsequent studies on other species show that IAPs not only exist in baculoviruses, but also homologs of IAPs have also been found in the other organisms such as fruit flies, nematodes, yeast, mammals and humans (Silke and Vaux, 2001; Zhen and Zhang, 2004). Among them, the relative research mainly focuses on fruit flies, silkworms, mosquitoes and ticks, but rare studies are found in insects (Hay et al., 1995; Huang et al., 2001; Li et al., 2007). Structurally, IAPs family is highly conversed, and there are usually one to three repeating sequences called baculoviral IAP repeat at the N-terminus, which are also called BIR (Baculovirus IAP Repeats) motifs or BIR domains. Additionally, the C-terminus of most IAPs also contains a structure that called RING domain (Williams et al., 2009).

Chilo suppressalis (Walker) is a significant rice pest which is widely distributed in the southern and northern regions of China and is severely damaged. The current prevention situation is still severe (Lu et al., 2013; Luo et al., 2014), which has triggered new biological control methodologies such as the use of natural enemies (He et al., 2013; Gao et al., 2019). Cotesia chilonis (Munakata) is widely distributed in China and has become the dominant parasitic wasp of C. suppressalis as global temperatures have risen (Chen et al., 2002; Pan et al., 2016; Pan et al., 2018). Therefore, the use of C. chilonis to prevent and control C. suppressalis has certain application value (Li et al., 2019; He et al., 2021). However, under the global climate warming, the mechanism of mutual adaptation between C. suppressalis and C. chilonis is not obscure against the background of global warming. And the influence of different factors on IAP genes expression of C. chilonis is well worth exploring. Thus, in this study 2 C. chilonis IAP genes (CcIAP1 and CcIAP) were found and identified and their structural characteristics were described. The expression of these two genes in response to temperature stress along with their expression in different developmental stages were studied. The impact of C. chilonis-mediated endoparasitism on the expression of these two genes in C. suppressalis was also examined. Moreover, RNAi was used to further clarify the function of C. chilonis IAP genes and explore relationships between pro- and anti-apoptotic genes in the host and parasitoid.

2 Materials and methods

2.1 Experimental insects

C. suppressalis and C. chilonis were collected from a suburb of Yangzhou (32.39°N, 119.42°E) and reared under the laboratory conditions at 27 ± 1°C, 60%–70% RH and a 16:8 h (light/dark) photoperiod (Pan et al., 2018). C. suppressalis were reared on an artificial diet (Gao et al., 2019). C. chilonis adults were supplied with a 10% honey/water solution and propagated using fifth instar larvae of C. suppressalis as hosts. Both C. suppressalis and C. chilonis were reared successively for three or more generations. Ethical review and approval was not required for the study on animals in accordance with the local legislation and institutional requirements.

2.2 Sample treatments

For different developmental stages treatment, a single fifth instar of C. suppressalis was placed in a test tube and C. chilonis adults were added in a 2:1 female/male ratio for breeding. Insects were incubated for 6 h at 27°C in darkness to facilitate parasitism of C. suppressalis; once parasitism occurred, C. suppressalis larvae were allowed to feed on artificial diet. Insects were maintained using the 27°C regimes described above until C. chilonis emerged from C. suppressalis. Three parasitized C. suppressalis were dissected from the two treatments on a daily basis, and the development of five randomly-selected C. chilonis individuals was inspected; different developmental stages of C. chilonis (first instar larvae, second instar larvae, third instar larvae, pupae and adults) were collected and stored at −80°C until needed. Treatments either contained 30 (larvae) or five individuals (adults and pupae). All treatments were replicated three times.

For parasitic time treatment, after parasitism for 1 h, 10 h, 24 h, 2 days, 3 days, 4 days, 5 days, 6 days, 7 days, 8 days and 9 days, the fifth instar larvae of C. suppressalis were collected. Four replicate groups for each interval of parasitism were collected and stored at −80°C.

For temperature treatment, 1-day old adults of C. chilonis were subjected to −13, −12, −9, −6, −3, 0, 27, 30, 33, or 36°C for 1 h in a constant-temperature incubator (Pan et al., 2018); samples were then placed in a climate-controlled incubator and allowed to recover at 27°C for 1 h. Each treatment contained 30 1-day old adults, and all treatments were replicated 3 times.

2.3 RNA isolation and synthesis of first strand cDNA

Total RNA was extracted from C. chilonis and C. suppressalis using the SV Total RNA Isolation System (Promega, USA). The integrity of RNA was verified by comparing ribosomal RNA bands in ethidium bromide-stained gels, and RNA purity was examined using spectrophotometric measurements at A₂₆₀ and A₂₈₀ nm (NanoDrop One, Thermo Fisher Scientific, United States). The first strand of cDNA was generated with an oligo (dT) ₁₈ primer (Fermentas, Canada).

2.4 Cloning and genome amplification

Partial gene sequences were obtained from the C. chilonis transcriptome, and specific primers were designed to verify fragments using the first strand of cDNA as template (Supplementary Table S1). Full-length cDNA sequences of genes were obtained with 5′- and 3′- RACE (SMARTer™ RACE, Clontech), and gene-specific primers were designed for verifying full-length cDNA sequences using the 5′-RACE template (Supplementary Table S1). Genomic DNA of C. chilonis adults was extracted using the Axyprep™ Multisource Genomic DNA Kit (Axygen, United States), and primers (Supplementary Table S1) were designed to amplify genomic fragments of CcIAP1 and CcIAP for subsequent cloning or sequencing.

2.5 Sequence analysis of genes

ORFs (Open reading frames) were identified with ORF Finder (https://www.ncbi.nlm.nih.gov/orffinder/), and deduced amino acid sequences were aligned with ClustalX. Sequence analysis tools on the ExPASy Molecular Biology Server including Translate, Compute pI/MW, and Blast (Swiss Institute of Bioinformatics), were used to analyze the predicted sequences. Motif Scan (https://prosite.expasy.org/) and InterPro (http://www.ebi.ac.uk/interpro/) were used to identify motifs characteristic of the IAPs family. Phylogenetic trees were constructed by the neighbor-joining minimum evolution, maximum likelihood and maximum parsimony methods with 1000 bootstrap replicates using MEGA 7.0 (Kumar et al., 2016).

2.6 Real-time qPCR analysis

Total RNA was isolated from the different treatments as described above, and the Bio-Rad iScript™ cDNA Synthesis Kit (Bio-Rad, USA) was used to reverse transcribe 0.5 µg total RNA into first strand cDNA. The primers used for real-time quantitative PCR (Supplementary Table S1) were designed based on the full-length cDNA sequence of genes. Real-time PCR reactions were conducted using SYBR Green I in a 20 μL volume that included 10 μL iTaq™ SYBR^® Green Supermix, 6 μL ddH₂O, 2 μL cDNA template and 1 μL each of the corresponding forward and reverse primers. Reaction conditions for PCR were as follows: 3 min initial denaturation step at 95°C, followed by 40 cycles of 15 s denaturation at 95°C, and 30s annealing at the Tm for each gene (Supplementary Table S1). Melting curve analysis was carried out to evaluate the homogeneity of the amplified PCR products. Each PCR reaction was replicated in triplicate. GAPDH was regarded as the high-temperature reference gene and H3 was regarded as the remaining treatment reference gene (Li et al., 2019).

2.7 Synthesis of dsRNA and RNAi experiments

Two full-length C. chilonis IAP genes were analyzed by the online website (http://sidirect2.rnai.jp); the regions for RNA silencing were determined, and gene-specific primers with an integrated T7 promoter were designed (Supplementary Table S1). A dsRNA specific to gfp encoding green fluorescence protein served as a control (Supplementary Table S1). The amplified products were sequenced and used as template DNA for dsRNA synthesis as recommended for the MEGAscript^®RNAi Kit (Thermo Fisher Scientific, USA). The quality and integrity of dsRNA molecules were evaluated by spectrophotometry and gel electrophoresis. The dsRNAs were stored at −80°C until needed.

Third instar larvae of C. suppressalis parasitized by C. chilonis for 9 days were selected for the RNAi experiments. The dosages of the dsRNA solution and nanomaterial (provided by Prof. Jie Shen, China Agricultural University) were 0.25 μL and 13.71 μL, respectively, for one parasitized third instar of C. suppressalis. The dsRNA solution and nanomaterial were mixed for 10–15 min, and detergent (provided Prof. Jie Shen) was added and incubated for 10–15 min; this solution was then placed on top of parasitized third instar larvae of C. suppressalis for 10–15 min in darkness. At 24 h post-treatment, parasitized C. suppressalis were dissected. C. suppressalis and C. chilonis individuals were collected for RNA extraction, and RNAi efficiency was analyzed by qRT-PCR. Photos were taken with a fluorescence microscope (Nikon ECLIPSE Ts2R) (Deerfield, IL, German) and a KEYENCE VHX-5000 system. The experiment was performed with three replicates per treatment. Non-parasitized C. suppressalis and dsGFP were used as controls.

2.8 Statistical analysis

Gene expression was normalized using the 2^−ΔΔCt method. Differences in mean values were analyzed with the independent-sample t-test and ANOVA. Levene’s test was used to evaluate homogeneity of variances, and Tukey’s test was used to assess significant differences. SPSS v.16.0 (SPSS, Chicago, IL, USA) was used for statistical analysis and presented as means ± SE (standard error).

3 Results

3.1 Characteristics of sequenced genes

The full-length cDNA sequence of CcIAP1 gene was 1404 bp (GenBank accession no. MZ076685) and contained a 210-bp 5′ untranslated region (UTR), a 1155-bp ORF, and 39-bp 3′ UTR. The full-length cDNA of CcIAP gene was 1184 bp (GenBank accession no. MZ076686) and contained a 190-bp 5′ UTR, a 951-bp ORF, and a 43-bp 3’ UTR. Comparison of cDNA and genomic sequences revealed the absence of introns in both CcIAP1 and CcIAP genes. The predicted CcIAP1 protein contained 384 amino acids, with theoretical molecular mass of 43.32 kDa and an isoelectric point (pI) of 6.33. The deduced CcIAP protein contained 316 amino acids with a predicted mass of 35.88 kDa and pI of 5.73. Analysis using Motif Scan and InterPro revealed that CcIAP1 and CcIAP both contained two BIR repeat signatures (residues 59–123 and 202–267 for CcIAP1; residues 59–123 and 202–267 for CcIAP). Furthermore, CcIAP1 also contained a putative zinc finger (residues 337–372). Multiple sequence alignments revealed that CcIAP1 shared 55.53%, 55.61%, 54.63%, 53.73% and 54.22% identity with IAP1 in Dufourea novaeangliae, Melipona quadrifasciata, Habropoda quadrifasciata, Bombus impatiens and Bombus terrestris, respectively. CcIAP shared 46.39%, 44.88%, 45.61%, 43.86% and 44.34% identity with IAP1 in Dufourea novaeangliae, Melipona quadrifasciata, Habropoda quadrifasciata, Bombus impatiens and B. terrestris, respectively (Figure 1).

FIGURE 1

FIGURE 1. Amino acid sequence alignment of CcIAP1 and CcIAP from Cotesia chilonis with orthologous proteins in Dufourea novaeangliae, Melipona quadrifasciata, Habropoda quadrifasciata, Bombus impatiens and Bombus terrestris. Identical amino acids are shaded with the same color. The two BIR sequences are marked by black and red line, respectively. The RING domain is underscored in green, respectively. Accession numbers of species are noted in Supplementary Table S2.

3.2 Phylogenetic analysis of genes

Previous research indicated that insect IAPs clustered with baculovirus IAP3, which suggests acquisition of IAPs via horizontal gene transfer (Chen et al., 2020). To explore evolutionary relationships, BlastP (https://blast.ncbi.nlm.nih.gov/Blast.cgi) was used to search for proteins related to baculovirus IAP3. Phylogenetic analysis of IAPs was performed using neighbor-joining, maximum likelihood, maximum parsimony and minimum evolution methods, and similar results were obtained with the different methods. The dendrogram in Figure 2 is representative of the phylogenetic groupings and was obtained with the neighbor-joining method; it clearly shows that CcIAP1 and CcIAP are closely related to IAP1s in other members of the Hymenoptera and IAP1 in C. suppressalis. Furthermore, phylogenetic analysis revealed that baculovirus IAP3s are closely related to insect IAPs in the orders Lepidoptera, Diptera and Hymenoptera. Analysis of the IAPs using ClustalX and MEGA 7.0 further confirmed that CcIAP1 and CcIAP share high similarity with IAP1s in other hymenopteran insects.

FIGURE 2

FIGURE 2. Phylogenetic tree based on the different species IAPs and baculovirus IAP3 proteins. Trees were generated with MEGA 7.0, and solid circles indicate the location of CcIAP1 and CcIAP Numbers on the branches are bootstrap values (1000 replicates), and only bootstrap values >50 are shown. Accession numbers are provided in Supplementary Table S2.

3.3 Gene expression during different developmental stages

CcIAP1 and CcIAP genes both showed variable expression in the different developmental stages of C. chilonis. CcIAP1 expression was significantly upregulated in the third instar larvae as compared to the second instar larvae (Figure 3A, F_4,11 = 3.365, p = 0.05), and CcIAP expression was significantly higher in the third instar larvae as compared to adults (Figure 3B, F_4,13 = 3.891, p < 0.05). CcIAP1 and CcIAP genes expression levels in the third instar larvae were 3.03- and 2.62-fold higher than expression in the second instar larvae, respectively.

FIGURE 3

FIGURE 3. Relative mRNA expression levels of CcIAP1 (A) and CcIAP (B) at different developmental stages. Statistics are presented as means ± SE. Columns labeled with different letters indicate significance between temperatures using one-way ANOVA followed by Tukey’s multiple comparison analysis (p <0.05).

3.4 Gene expression in response to different parasitic time

When compared for different durations of parasitism, CcIAP1 and CcIAP genes expression patterns were similar but not identical (Figure 4). When compared to the control (1 h of parasitism), CcIAP1 gene was significantly upregulated at 3, 6, 7, 8 and 9 days of parasitism (CcIAP1: F_10,33 = 1869.450, p < 0.001); CcIAP gene was significantly upregulated at 2, 3, 6, 7, 8 and 9 days (CcIAP: F_10,33 = 423.424, p < 0.001). The highest expression levels of CcIAP1 and CcIAP genes were observed at 7 days where levels were 1026.73- and 855.25-fold higher than the control (1 h), respectively.

FIGURE 4

FIGURE 4. Relative mRNA expression levels of CcIAP1 (A) and CcIAP (B) in response to different parasitic time. Statistics represent means ± SE. Columns labeled with different letters indicate significant differences between times using one-way ANOVA followed by Tukey’s multiple comparison analysis (p < 0.05).

3.5 Gene expression in response to different temperatures

The relative mRNA levels of CcIAP1 and CcIAP genes were monitored at temperature gradients ranging from −13°C to 36°C. The expression of both genes were both upregulated by heat and cold stress (CcIAP1: F_9,20 = 132.603, p < 0.001; CcIAP: F_9,20 = 54.613, p < 0.001). CcIAP1 gene was remarkably upregulated at −3°C, −6°C, −9°C, and 30°C (Figure 5A), and CcIAP gene was significantly upregulated at −6°C and 30°C (Figure 5B). The expression levels of both CcIAP1 and CcIAP genes were highest at 30°C and were 24.81- and 39.56-fold greater than the control at 27°C, respectively.

FIGURE 5

FIGURE 5. Relative mRNA expression levels of CcIAP1 (A) and CcIAP (B) at different temperatures. Data represent means ± SE. Columns labeled with different letters indicate significance using one-way ANOVA followed by Tukey’s multiple comparison analysis (p <0.05).

3.6 Knockdown of CcIAP1 and CcIAP genes expression

RNA interference studies were conducted by covering parasitized third instar larvae of C. suppressalis with a solution containing nanocarrier, dsRNA, and detergent, which successfully penetrated the body wall cavity of C. suppressalis (Supplementary Figure S1). The expression of CcIAP1 and CcIAP genes were both significantly reduced when the nanocarrier/dsCcIAP1/detergent or nanocarrier/dsCcIAP/detergent mixture was used for RNAi as compared to the nanomaterial/dsGFP/detergent control (CcIAP1: t = 3.647, p < 0.05; CcIAP: t = 3.357, p < 0.05) (Figure 6). However, the silencing efficiency of dsCcIAP1 was higher than dsCcIAP; treatment with the former reduced the CcIAP1 transcript levels by 78% while the latter reduced CcIAP levels by 55% as compared with the control.

FIGURE 6

FIGURE 6. Relative mRNA expression levels of CcIAP1 (A) and CcIAP (B) in Cotesia chilonis after RNAi for 24 h. Statistics are presented as means ± SE. Data were analyzed using independent samples t-test built in SPSS software. p < 0.05 was considered statistically significant. Asterisks represent significant differences between control and RNAi treatment.

There was a significant upregulation in expression of the apoptosis factor encoded by Cscaspase-1 in parasitized C. suppressalis treated with the nanomaterial/dsCcIAP1/detergent solution as compared to the control (F_3,12 = 5.498, p < 0.05) (Figure 7). Interestingly, expression of the apoptosis effector encoded by Cscaspase-1 was not significantly different from the dsGFP control when samples treated with the nanomaterial/dsCcIAP/detergent mixture.

FIGURE 7

FIGURE 7. Relative mRNA expression levels of Cscaspase-1 in Chilo suppressalis after silencing IAP genes in Cotesia chilonis for 24 h. Statistics represent means ± SE, and columns labeled with different letters indicate significant differences between different treatments using one-way ANOVA followed by Tukey’s multiple comparison analysis (p < 0.05).

4 Discussion

Apoptosis is a genetically and biochemically controlled process and plays crucial roles in the development, issue homeostasis and defense of multicellular organisms. Inhibitor of apoptosis proteins (IAPs) are the class of proteins that can effectively inhibit cell apoptosis (Hamajima et al., 2016). In this study, we cloned and characterized a typical IAP1 gene (CcIAP1) and an atypical IAP1 gene (CcIAP). Multiple sequence alignments and phylogenetic analysis show that CcIAP1 and CcIAP are highly conserved and closely related to orthologues in other hymenopteran insects. CcIAP1 contains two BIR domains and one RING domain while CcIAP only contains two BIR domains, which are characteristic domains of IAPs that can bind one Zn²⁺, and its function is mainly to mediate the interaction between proteins. In addition, The BIR domains are divided into Type 1 and Type 2 according to whether there is a peptide binding groove. The Type 1 lacks the peptide binding groove and can interact with TRAF (tumor necrosis factor receptor-associated factor) and TGF-β (transforming growth factor-β) to regulate their signal pathways (Li et al., 2007). In this study, CcIAP1 and CcIAP both contain Type 2 BIR domains with unique hydrophobic clefts that can be combined with caspase or IAP binding motifs (IBMs) present in IAP antagonists (Verhagen et al., 2007). Most of RING domains that can bind two Zn²⁺ have E3 ubiquitin ligase activities (Feltham et al., 2011; Middleton et al., 2014; Takeda et al., 2014; Clem, 2015). Among them, the combination of BIR domain and caspase plus the participation of the RING domain with ubiquitin ligase (E3) activities can exert the full anti-apoptotic function of IAP (Schile et al., 2008). Furthermore, the two genes both lack introns in C. chilonis. And it has been reported that genes without introns or containing short introns may be expressed at higher levels than genes with multiple or long introns in response to stress (Comeron, 2004; Gao et al., 2019). Previous studies have shown that baculovirus IAP3, which has the most obvious anti-apoptosis function, has a close evolutionary relationship with insect IAPs (Birnbaum et al., 1994; Carpes et al., 2005; Liang et al., 2012). Our phylogenetic results suggest that baculovirus IAP3 genes might be derived from Lepidoptera, Diptera or Hymenoptera through horizontal gene transfer; this speculation is similar to previous results showing that Lepidoptera and Diptera are the primary hosts of baculoviruses (Chen et al., 2020).

The growing development and morphologic changes of organism are in closely connection with cell apoptosis (Elmore, 2007; Kerr et al., 2015). Apoptosis is a normal process in insect development, especially in holometabolous species, and plays a significant role in maintaining homeostasis. IAPs inhibit apoptosis via the BIR domains that mediate protein (Kerr et al., 2015). In different developmental stage experiments, the expression of CcIAP1 and CcIAP were both induced in the third instar larvae stage, which corresponded to predicted times of feeding and growth of C. chilonis. Moreover, C. chilonis that just parasitized have relatively sufficient nutrient support, but with the growth and development, their demand for living space and nutrients increases, and the competition among individuals intensifies. Therefore, it is reasonable to assume that the expression of CcIAP1 and CcIAP genes increased to meet nutritional demands in the third instar larval stage. Parasitism of C. suppressalis by C. chilonis caused a general increase in CcIAP1 and CcIAP genes expression, and mRNA levels at the 3, 6, 7, 8, and 9 days time points were significantly higher than the control, especially at the 7 days time point. Our results showed that the expression levels of CcIAP1 and CcIAP genes were higher in late larval stages as compared to earlier stages, which is consistent with the premise that C. chilonis may inhibit apoptosis in the late larval stages to foster growth and development. One explanation for elevated expression at the 3 days time point is to provide support for the egg-to-larvae transition. As for the downregulation at 4 and 5 days, we speculate that the needs of young C. chilonis larvae have little effect on the development of C. suppressalis and the host-parasite interaction is somewhat balanced at these time points. In organisms, IAPs bind to active caspases that execute the cell death program, and prior research has shown that caspases function in regulation of the mitochondrial apoptotic pathway (Ehrmann et al., 2023). In this study, CcIAP1 and CcIAP genes were both silenced by nanocarrier and dsCcIAP1 or dsCcIAP. These results indicated that nanocarrier could help dsRNA pass through the host silence the objective genes of wasp in the host. Further analyses exhibited Cscaspase-1 expression in C. suppressalis was significantly upregulated when CcIAP1 gene was silenced, which suggests that CcIAP1 could contribute to the inhibition of apoptosis in C. suppressalis. In a related study, the induction of apoptosis could be blocked by the overexpression of Sf-IAP1 in Spodoptera frugiperda (Veeran et al., 2020). Furthermore, our results show that Cscaspase-1 and CcIAP1 play a vital role in the parasitoid-pest interaction, and these complex interactions warrant further study.

Temperature is a critical environmental factor that affects insect growth, development, distribution and abundance (Bale et al., 2002). When subjected to high or low temperature stress, insects may adopt different coping strategies, such as avoidance behaviors or changing physiological functions to tolerate temperature stress (Hoffmann and Parsons, 1991). Some studies shown that external stimulus signals could mediate cell apoptosis, and XIAP play a synergistic role in the mechanism of pro-apoptosis (Igaki et al., 2002; Ludivine et al., 2010). And IAPs are at the center of this system to preserve sudden damages and death (Hanifeh and Ataei, 2022). In this study, we found that the expression of CcIAP1 and CcIAP genes were both significantly upregulated at −6°C and 30°C in response to temperature stress, well indicating that temperature as the abiotic stress plays a part in the apoptosis pathway of cells and the temperature tolerance is also related to the expression and regulation of IAP genes in insects, which had not been shown before. Our experimental results also indicated that temperature tolerance may have evolved in response to long exposure to moderately high or low temperatures rather than to an acute bout of extremely high or low temperatures, which is similar to the results of previous studies (Addo-Bediako et al., 2000; Hance et al., 2007). In addition, heat shock proteins have been reported to regulate apoptosis and cell death (Takayama et al., 2003), suggesting that the mechanism of temperature stress between heat shock proteins and inhibitor of apoptosis proteins is a great direction worthy of study.

In conclusion, two IAP genes (CcIAP1 and CcIAP) in C. chilonis were firstly cloned and characterized, both of them were expressed during different developmental stages, especially in the third instar larvae stage. And the expression of them were both significantly upregulated in response to temperature stress and parasitism. The current study furthers our understanding of the regulation of IAP genes of insect under development and temperature stress, and provides groundwork for future studies aimed at understanding the molecular mechanisms that are driving the evolutionary adaptation of C. chilonis.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

The manuscript presents research on animals that do not require ethical approval for their study.

Author contributions

M-XL: Conceptualization, Investigation, Writing–original draft. F-JH: Investigation, Software, Writing–original draft. FZ: Investigation, Methodology, Writing–review and editing. Y-ZD: Conceptualization, Funding acquisition, Writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Key R&D Program of China (2017YFD0200400), and “Six Talent Peaks” High-level Talent Project (NY-088).

Conflict of interest

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1328167/full#supplementary-material

References

Addo-Bediako A., Chown S. L., Gaston K. J. (2000). Thermal tolerance, climatic variability and latitude. Proc. R. Soc. Lond. B 267, 739–745. doi:10.1098/rspb.2000.1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Bale J. S., Masters G. J., Hodkinson I. D., Awmak C., Bezemer T. M., Brown V., et al. (2002). Herbivory in global climate change research: direct effects of rising temperature on insect herbivores. Glob. Change Biol. 8, 1–16. doi:10.1046/j.1365-2486.2002.00451.x

CrossRef Full Text | Google Scholar

Birnbaum M. J., Clem R. J., Miller L. K. (1994). An apoptosis-inhibiting gene from a nuclear polyhedrosis virus encoding a polypeptide with Cys/His sequence motifs. J. Virol. 68, 2521–2528. doi:10.1128/JVI.68.4.2521-2528.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpes M. P., de Castro M. E., Soares E. F., Villela A. G., Pinedo F. J., Ribeiro B. M. (2005). The inhibitor of apoptosis gene (iap-3) of Anticarsia gemmatalis multicapsid nucleopolyhedrovirus (AgMNPV) encodes a functional IAP. Arch. Virol. 150, 1549–1562. doi:10.1007/s00705-005-0529-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen H. C., Lou Y. G., Cheng J. A. (2002). The research and application of Cotesia chilonis. Chin. J. Biol. 18, 90–93.

Google Scholar

Chen P., Kang T. T., Bao X. Y., Dong Z. Q., Zhu Y., Xiao W. F., et al. (2020). Evolutionary and functional analyses of the interaction between the Bombyx mori inhibitor of apoptosis (IAP) and nucleopolyhedrovirus IAPs. Insect Sci. 27, 463–474. doi:10.1111/1744-7917.12664

PubMed Abstract | CrossRef Full Text | Google Scholar

Clem R. J. (2015). Viral IAPs, then and now. Cell Dev. Biol. 39, 72–79. doi:10.1016/j.semcdb.2015.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Comeron J. M. (2004). Selective and mutational patterns associated with gene expression in humans: influences on synonymous composition and intron presence. Genetics 167, 1293–1304. doi:10.1534/genetics.104.026351

PubMed Abstract | CrossRef Full Text | Google Scholar

Crook N. E., Clem R. J., Miller L. K. (1993). An apoptosis-inhibiting baculovirus gene with a zinc finger-like motif. J. Virol. 67, 2168–2174. doi:10.1128/JVI.67.4.2168-2174.1993

PubMed Abstract | CrossRef Full Text | Google Scholar

Danial N. N., Korsmeyer S. J. (2004). Cell death: critical control points. Cell 116, 205–219. doi:10.1016/s0092-8674(04)00046-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrmann J. F., Grabarczyk D. B., Heinke M., Deszcz L., Kurzbauer R., Hudecz O., et al. (2023). Structural basis for regulation of apoptosis and autophagy by the BIRC6/SMAC complex. Science 379, 1117–1123. doi:10.1126/science.ade8873

PubMed Abstract | CrossRef Full Text | Google Scholar

Elmore S. (2007). Apoptosis: a review of programmed cell death. Toxicol. Pathol. 35, 495–516. doi:10.1080/01926230701320337

PubMed Abstract | CrossRef Full Text | Google Scholar

Everett H., McFadden G. (1999). Apoptosis: an innate immune response to virus infection. Trends Microbiol. 7, 160–165. doi:10.1016/s0966-842x(99)01487-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Feltham R., Bettjeman B., Budhidarmo R., Mace P. D., Shirley S., Condon S. M., et al. (2011). Smac mimetics activate the E3 ligase activity of cIAP1 protein by promoting RING domain dimerization. J. Biol. Chem. 286, 17015–17028. doi:10.1074/jbc.M111.222919

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao P., Lu M. X., Pan D. D., Du Y. Z. (2019). Characterization of an inducible HSP70 gene in Chilo suppressalis and expression in response to environmental and biological stress. Cell Stress Chaperon 25, 65–72. doi:10.1007/s12192-019-01047-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamajima R., Iwamoto A., Tomizaki M., Suganuma I., Kitaguchi K., Kobayashi M., et al. (2016). Functional analysis of inhibitor of apoptosis 1 of the silkworm Bombyx mori. Insect biochem. molec. 79, 97–107. doi:10.1016/j.ibmb.2016.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hance T., Van-Baaren J., Vernon P., Boivin G. (2007). Impact of extreme temperatures on parasitoids in a climate change perspective. Annu. Rev. Entomol. 52, 107–126. doi:10.1146/annurev.ento.52.110405.091333

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanifeh M., Ataei F. (2022). XIAP as a multifaceted molecule in cellular signaling. Apoptosis 27, 441–453. doi:10.1007/s10495-022-01734-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hay B. A., Wassarman D. A., Rubin G. M. (1995). Drosophila homologs of baculovirus inhibitor of apoptosis proteins function to block cell death. Cell 83, 1253–1262. doi:10.1016/0092-8674(95)90150-7

PubMed Abstract | CrossRef Full Text | Google Scholar

He F. J., Zhu F., Lu M. X., Du Y. Z. (2021). Comparison of morphology, development and expression patterns of hsf and hsp11.0 of Cotesia chilonis under normal and high temperature. Peer J. 9, e11353. doi:10.7717/peerj.11353

CrossRef Full Text | Google Scholar

He Y. P., Zhang J. F., Gao C. F., Su J. Y., Chen J. M., Shen J. L. (2013). Regressiom analysis of dynamics of insecticide resistance in field populations of Chilo suppressalis (Lepidoptera: crambidae) during 2002-2011 in China. J. Econ. Entomol. 106, 1832–1837. doi:10.1603/ec12469

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann A. A., Parsons P. A. (1991). Evolutionary genetics and environmental stress. New York: Oxford University Press.

Google Scholar

Huang Q., Deveraux Q. L., Maeda S., Stennicke H. R., Hammock B. D., Ammock B. D., et al. (2001). Cloning and characterization of an inhibitor of apoptosis protein (IAP) from Bombyx mori. Biochim. Biophys. Acta. 1499, 191–198. doi:10.1016/s0167-4889(00)00105-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Igaki T., Yamamoto-goto Y., Tokushige N., Kanda H., Miura M. (2002). Down-regulation of DIAP1 triggers a novel Drosophila cell death pathway mediated by Dark and DRONC. J. Biol. Chem. 277, 23103–23106. doi:10.1074/jbc.C200222200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ikeda M., Yamada H., Hamajima R., Kobayashi M. (2013). Baculovirus genes modulating intracellular innate antiviral immunity of lepidopteran insect cells. Virology 435, 1–13. doi:10.1016/j.virol.2012.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerr J. F., Winterford C. M., Harmon B. V. (2015). Apoptosis. Its significance in cancer and cancer therapy. Cancer 73, 2013–2026. doi:10.1002/1097-0142(19940415)73:8<2013::aid-cncr2820730802>3.0.co;2-j

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerr J. F., Wyllie A. H., Currie A. R. (1972). Apoptosis: a basic biological phenomenon with wide-ranging implications in tissue kinetics. Brit. J. Cancer. 26, 239–257. doi:10.1038/bjc.1972.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar S., Stecher G., Tamura K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi:10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Li Q., Li H., Blitvich B. J., Zhang J. (2007). The Aedes albopictus inhibitor of apoptosis 1 gene protects vertebrate cells from bluetongue virus-induced apoptosis. Insect Mol. Biol. 16, 93–105. doi:10.1111/j.1365-2583.2007.00705.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li Q. Y., Li Z. L., Lu M. X., Cao S. S., Du Y. Z. (2019). Selection of valid reference genes for quantitative real-time PCR in Cotesia chilonis (Hymenoptera: braconidae) exposed to different temperatures. PLoS ONE 14, e0226139. doi:10.1371/journal.pone.0226139

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang C., de Lange J., Chen X., van Oers M. M., Vlak J. M., Westenberg M. (2012). Functional analysis of two inhibitor of apoptosis (iap) orthologs from Helicoverpa armigera nucleopolyhedrovirus. Virus Res. 165, 107–111. doi:10.1016/j.virusres.2012.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu M. X., Du Y. Z., Liu Z. X., Hua J., Liu P. Y., Li J. Y. (2013). Diapause, signal and molecular characteristics of overwintering Chilo suppressalis. Sci. Rep. 3, 3211. doi:10.1038/srep03211

PubMed Abstract | CrossRef Full Text | Google Scholar

Ludivine A., Pradelli M. B., Jean-Ehrland R. (2010). Mitochondrial control of caspase-dependent and -independent cell death. Cell. Mol. Life Sci. 67, 1589–1597. doi:10.1007/s00018-010-0285-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo G. H., Li X. H., Han Z. J., Guo H. F., Yang Q., Wu M., et al. (2014). Molecular characterization of the piggyBac-like element, a candidate marker for phylogenetic research of Chilo suppressalis (Walker) in China. BMC Mol. Biol. 15, 28–29. doi:10.1186/s12867-014-0028-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Luque T., Finch R., Crook N., O'Reilly D. R., Winstanley D. (2001). The complete sequence of the Cydia pomonella granulovirus genome. J. Gen. Virol. 82, 2531–2547. doi:10.1099/0022-1317-82-10-2531

PubMed Abstract | CrossRef Full Text | Google Scholar

Middleton A. J., Budhidarmo R., Day C. L. (2014). Use of E2∼ubiquitin conjugates for the characterization of ubiquitin transfer by RING E3 ligases such as the inhibitor of apoptosis proteins. Method. Enzymol. 545, 243–263. doi:10.1016/B978-0-12-801430-1.00010-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan D. D., Cao S. S., Lu M. X., Hang S. B., Du Y. Z. (2018). Genes encoding heat shock proteins in the endoparasitoid wasp, Cotesia chilonis, and their expression in response to temperatures. J. Integr. Agr. 17, 1012–1022. doi:10.1016/s2095-3119(17)61737-4

CrossRef Full Text | Google Scholar

Pan D. D., Liu Z. X., Lu M. X., Cao S. S., Yan W. F., Du Y. Z. (2016). Species and occurrence dynamics of parasitic wasps of the rice stem borer, Chilo suppressalis (Walker) (Lepidoptera: Pyralidae) in Yangzhou. J. Environ. Entomol. 38, 1106–1113. doi:10.3969/j.issn.1674-0858.2016.06.6

CrossRef Full Text | Google Scholar

Schile A. J., Garcia-Fernandez M., Steller H. (2008). Regulation of apoptosis by XIAP ubiquitin-ligase activity. Genes Dev. 22, 2256–2266. doi:10.1101/gad.1663108

PubMed Abstract | CrossRef Full Text | Google Scholar

Silke J., Vaux D. L. (2001). Two kinds of BIR-containing protein-inhibitors of apoptosis, or required for mitosis. J. Cell Sci. 114, 1821–1827. doi:10.1242/jcs.114.10.1821

PubMed Abstract | CrossRef Full Text | Google Scholar

Takayama S., Reed J. C., Homma S. (2003). Heat-shock proteins as regulators of apoptosis. Oncogene 22, 9041–9047. doi:10.1038/sj.onc.1207114

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeda A. N., Oberoi-Khanuja T. K., Glatz G., Schulenburg K., Scholz R. P., Carpy A., et al. (2014). Ubiquitin-dependent regulation of MEKK2/3-MEK5-ERK5 signaling module by XIAP and cIAP1. EMBO. J. 33, 1784–1801. doi:10.15252/embj.201487808

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandergaast R., Mitchell J. K., Byers N. M., Friesen P. D. (2015). Insect inhibitor-of-apoptosis (IAP) proteins are negatively regulated by signal-induced N-terminal degrons absent within viral IAP proteins. J. Virol. 89, 4481–4493. doi:10.1128/JVI.03659-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaux D. L., Korsmeyer S. J. (1999). Cell death in development. Cell 96, 245–254. doi:10.1016/s0092-8674(00)80564-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Veeran S., Zhang J., Zhong G. (2020). Pro-apoptotic function analysis of the reaper homologue IBM1 in Spodoptera frugiperda. Int. J. Mol. Sci. 21, 27–29. doi:10.3390/ijms21082729

CrossRef Full Text | Google Scholar

Verhagen A. M., Kratina T. K., Hawkins C. J., Silke J., Ekert P. G., Vaux D. L. (2007). Identification of mammalian mitochondrial proteins that interact with IAPs via N-terminal IAP binding motifs. Cell Death Differ. 14, 348–357. doi:10.1038/sj.cdd.4402001

PubMed Abstract | CrossRef Full Text | Google Scholar

Williams T., Chitnis N. S., Bilimoria S. L. (2009). Invertebrate iridovirus modulation of apoptosis. Virol. Sin. 24, 295–304. doi:10.1007/s12250-009-3060-1

CrossRef Full Text | Google Scholar

Zhen H. N., Zhang X. (2004). New inhibitor of apoptosis Livin. J. Fourth Mil. Med. Univ. 25, 1822–1823.

Google Scholar