Molecular mechanism by which Apis cerana cerana MKK6 (AccMKK6)-mediated MAPK cascades regulate the oxidative stress response

Mitogen-activated protein kinase kinases (MKKs) are important components of the MAPK signaling pathways, which play a key role in responding to stress and inflammatory stimuli. Here, a new MKK gene, AccMKK6, was identified and functionally analyzed in Apis cerana cerana. Real-time quantitative PCR (qPCR) and Western blot analysis demonstrated that the AccMKK6 expression level was up-regulated by several environmental stresses. Moreover, the knockdown of AccMKK6 by RNA interference technology altered the expression levels of some antioxidant genes. In addition, the knockdown of AccMKK6 resulted in increased malonyldialdehyde (MDA) concentration and decreased antioxidant-related enzymes activity in honeybees. To explore the MAPK signaling pathways involved in AccMKK6, we identified the transcription factor kayak in A. cerana cerana. We analyzed the interactions of AccMKK6, Accp38b, and Acckayak using the yeast two-hybrid system. AccMKK6 and Acckayak showed similar expression profiles after several stress treatments. In addition, the expression level of Acckayak was significantly increased when AccMKK6 was silenced. Therefore, we speculate that AccMKK6 may be involved in the MAPK cascades, which play a crucial role in counteracting oxidative stress caused by external stimuli.


Introduction
Insects face a variety of adverse environmental pressures, including global warming, heavy metal pollution and the excessive use of pesticides. These environmental stresses can enhance the generation of reactive oxygen species (ROS), which cause oxidative damage [1][2][3]. ROS are formed as byproducts of aerobic metabolism, including superoxide anions, hydrogen peroxide, and hydroxyl radicals [4]. Generally, ROS are the second messenger in intracellular signaling pathways and maintain a dynamic equilibrium under normal conditions [5]. However, this equilibrium can be disrupted by adverse environmental stresses, which result in lipid peroxidation and in protein and DNA damage and can adversely affect cell viability by causing membrane damage and enzyme inactivation [6]. To maintain ROS homeostasis and prevent oxidative damage, organisms have developed complex antioxidant mechanisms to avoid oxidative damage [7]. Mitogen-activated protein kinase (MAPK) signaling pathways are critical in responding to oxidative stress [8,9].
MAPK signaling pathways transfer extracellular information to the interior of the cell, regulating cell activities in response to external stimuli [10,11]. These signaling pathways are activated by a variety of stimuli including cytokines, hormones, and diverse cellular stressors such as oxidative stress and endoplasmic reticulum (ER) stress [12]. MAPK signaling pathways have been implicated in cancer, neurodegenerative diseases, and the pathogenesis of several chronic diseases, as the accumulation of ROS and the activation of MAPK signaling pathways are considered to lead to disease in humans via the induction of

Insects and treatments
The honeybees (A. cerana cerana) used in this study were obtained from the experimental apiary of the College of Animal Science and Technology, Shandong Agricultural University (Taian, Shandong, China). Larvae [first (L1) and fourth (L4) instar larvae], pupae [prepupae (PP), white-eyed (Pw), pink-eyed (Pp), brown-eyed (Pb), and dark-eyed (Pd) pupae] and adult workers were collected. The experimental honeybees were kept in bee boxes at a constant temperature of 32 • C and 80% humidity under a 24 h dark regimen [27]. The adult workers were divided into nine groups of 48 individuals each. Each group was exposed to different harmful conditions. Group 1 was injected with 1 μl of 20 mM H 2 O 2 solution and analyzed after 0, 0.5, 1, 1.5, 2, or 2.5 h. The honeybees in groups 2 and 3 were injected with 1 μl of HgCl 2 and CdCl 2 (3 mg/ml), respectively, and analyzed after 0, 0.5, 1, 1.5, 2, or 2.5 h. Groups 4-6 of honeybees were injected with 1 μl of pesticide (paraquat, emamectin benzoate, acetamiprid) at concentrations of 20 μl/ml. Treated at 0 h honeybees were used as controls. The treated materials were immediately snap-frozen in liquid nitrogen and stored at −70 • C.

RNA extraction and cDNA synthesis
Total RNA was extracted from the honeybees using RNAiso Plus (TaKaRa, Japan) and then stored at −70 • C. The RNA was reverse transcribed into first-strand complementary DNA (cDNA) using 5× All-In-One RT MasterMix with the AccuRT Genomic DNA Removal Kit (Applied Biological Materials Inc., Richmond, BC, Canada), which facilitates the elimination of contaminating genomic DNA from RNA.

Isolation of the AccMKK6 and Acckayak open reading frame sequences
To clone the open reading frame (ORF) sequences of AccMKK6 (GenBank accession no. XP 016921098.1) and Acckayak (GenBank accession no. XP 016911846.1), specific primers (as shown in Supplementary Table S1) were designed based on the AccMKK6 and Acckayak gene sequences and synthesized by Biosune Biotechnological Company (Shanghai, China). The PCR mixtures were as described by Zhao et al. [28], and the amplification conditions are shown in Supplementary Table S2. The PCR products were purified using a Gel Extraction Kit (Solarbio, Beijing, China), ligated with pEASY-T1 vectors (TransGen, Biotech, Beijing, China), and then transformed into Escherichia coli cells (DH5α) for sequencing. The primers used in this study are listed in Supplementary Table S1.

Bioinformatics and phylogenetic analysis of AccMKK6 and Acckayak
The amino acid sequences of AccMKK6 and Acckayak were obtained from GenBank (https://www.ncbi.nlm.nih. gov/genbank/), and the isoelectric point and molecular mass were predicted by ExPASy (https://www.expasy.org/). Homologous AccMKK6 and Acckayak amino acid sequences from other species were obtained using the BLAST search from the NCBI (https://blast.ncbi.nlm.nih.gov/Blast.cgi), and multiple amino acid sequence alignments were performed using DNAMAN version 6.0.3. A neighbor-joining (NJ) phylogenetic tree analysis was performed on the homologous sequences of MKKs from other species using MEGA version 4.1.

Quantitative real-time PCR
To quantify the gene expression, real-time quantitative PCR was performed via the Bestar R One-step RT qPCR Kit (SyBR Green) (DBI Bioscience, China) and a CFX96TM Real-time PCR Detection System (Bio-Rad, Hercules, CA, U.S.A.). The specific qPCR primers for AccMKK6 and Acckayak were designed based on the cDNA sequence. The housekeeping gene β-actin (GenBank accession no. HM 640276) was selected as an internal control [29]. The qPCR was performed in a total volume of 20 μl with 8 μl of double-distilled water, 1 μl of cDNA template, 0.5 μl of each primer, and 10 μl of Bestar R qPCR Mastermix. The qPCR program was as follows: 95 • C for 30 s; 40 cycles of 95 • C for 5 s, 55 • C for 15 s, and 72 • C for 15 s; and a final melt cycle from 55 to 96 • C. All experimental samples were evaluated in triplicate, and the data were analyzed using CFX Manager Software version 1.1. The significant differences between the samples were identified by SPSS software version 17.0.

The protein expression of AccMKK6 and antibody preparation
Primers (M2F, M2R) containing BamHI and SacI restriction enzyme sites were used to insert the ORF of AccMKK6 into the expression vector pET-30a(+) (Novagen, Darmstadt, Germany). The recombinant plasmid pET-30a(+)-AccMKK6 was transformed into E. coli Transetta (DE3) (TransGen Biotech, Beijing, China). The bacterial solution was cultured in Luria Bertani (LB) with 50 μg/ml kanamycin at 37 • C for 1-2 h until the cell density reached 0.4-0.6 at 600 nm. The expression of the recombinant AccMKK6 protein was induced with 2 mM isopropyl-1-thio-β-d-galactopyranoside (IPTG) at 28 • C for 6-8 h, and the protein was then separated by 12% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). To prepare an antibody against AccMKK6, the recombinant AccMKK6 protein band was excised from the gel, milled with 1 ml of normal saline and injected into white mice. The mice were injected once a week for a total of four injections. Blood was collected on the fourth day after the last injection and stored at 37 • C for 1 h and then 4 • C for 6 h. After centrifugation at 3000 g for 15 min, the supernatant antibodies were stored at −70 • C.

Protein extraction and Western blot analysis
The total protein was extracted from the honeybee samples using a Tissue Protein Extraction Kit (CWbiotech, Beijing, China). Western blotting was performed according to the procedure described by Li et al. [30] with some modifications. The protein samples were separated by 12% SDS-PAGE and electrotransferred onto a polyvinylidene difluoride (PVDF) membrane. The membrane was incubated overnight with the primary antibody (anti-AccMKK6) at 1:1000 (v/v) dilution. Tubulin was used as the reference antibody at a 1:5000 (v/v). The results were observed using the Super-Signal West Femto Trial Kit (Thermo Scientific, U.S.A.). The Western blot results were analyzed using ImageJ 1.51j8 (National Institutes of Health, U.S.A.).

RNAi of AccMKK6
We used an RNAi experiment to knock down AccMKK6 in adult workers. Primers were designed to select regions of lower homology, and 23 bp of the T7 promoter sequence (TAATACGACTCACTATAGGGCGA) was added to the start of the AccMKK6 forward primer (M3F) and the end of the AccMKK6 reverse primer (M3R). The double-stranded RNAs (dsRNA) were prepared using the T7 RiboMAX TM Express RNAi System (Promega, U.S.A.) according to the manufacturer's instructions. We selected the green fluorescent protein gene (GFP, GenBank accession no. U87974), which has no homolog in the honeybees, as a control [31]. The dsRNA (6 μg/individual) was injected into adult honeybees between the first and second abdominal segments using a microsyringe [32]. In addition, no injected honeybees and honeybees injected with 1 μl of water were used as control groups. The silencing efficiency was determined by qPCR and Western blot.

Transcription level analysis of antioxidant genes and enzymatic activities after AccMKK6 knockdown
qPCR was performed to analyze the expression profiles of AccSOD1 (GenBank ID: JN700517), AccSOD2 (GenBank ID: JN637476), AccGSTS4 (GenBank ID: JN008721), AccCYP4G11 (GenBank ID: kc243984), AccGSTO2 (Gen-Bank ID: JX434029), and AccGSTD (GenBank ID: JF798573) after AccMKK6 was knocked down. The total protein was extracted from the whole adult bees 36 h after they had been injected with dsRNA and quantified using the BCA Protein Assay Kit (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The MDA concentration was measured using an MDA Assay Kit (Nanjing JianCheng Bioengineering Institute, Nanjing, China). The activity levels of the enzymes SOD and POD were measured using the SOD and POD Assay Kits (Nanjing JianCheng Bioengineering Institute, Nanjing, China), respectively, according to the manufacturer's instructions.

Matchmaker Gold yeast two-hybrid system
The Matchmaker Gold yeast two-hybrid system was used to detect the interactions between AccMKK6, Accp38b, and Acckayak, as described in a previous study [33]. The typical MAPK member Accp38b (GenBank accession no. GU321334) was isolated as previously described [34] and ligated into the pGBKT7 vector. The ORFs of AccMKK6 and Acckayak were cloned and fused to the pGADT7 vector. The plasmids Accp38b-pGBKT7 and AccMKK6-pGADT7, Accp38b-pGBKT7, and Acckayak-pGADT7 were co-transformed into the Y2H Gold yeast strain, respectively, and plated on DDO (SD/-Leu/-Trp) medium. Then, the positive clones were grown on selective SD (QDO/X, SD/-Ade/-His/-Leu/-Trp with X-α-gal and Aureobasidin A) medium.

Glutathione S-transferase pull-down experiments
The ORF of Accp38b was ligated into the pGEX-4T-1 vector with the glutathione S-transferase (GST)-tag. The ORF of AccAP2m and AccMKK6 was ligated into the expression vector pET-30a(+) with a His-tag. The recombinant plasmids were then transformed into E. coli Transetta (DE3). The pull-down assay was performed using Glutathione sepharose TM 4B (GE Healthcare Life Sciences, Sweden) and eluted with elution buffer. The proteins samples were separated by 12% SDS-PAGE, and the separated proteins were blotted onto PVDF membrane. Western blot analyses were performed using His or GST antibodies (Abmart, Shanghai, China).

Bioinformatics and phylogenetic analysis of AccMKK6
MKK6 is a novel redox-regulated Rac-binding protein that mediates MAPK cascades and plays important roles in ROS production and disease resistance in many animals [35]. To determine the functions of MKK6 in A. cerana cerana, the ORF of AccMKK6 was obtained by RT-PCR using the specific primers M1F and M1R. The AccMKK6 ORF consists of 1002 bp and encodes a 333 amino acid protein with a predicted average mass of 37.63 kDa and a theoretical isoelectric point (pI) of 7.15. Multiple sequence analysis of several MKK6s from different species revealed that the amino acid sequence of AccMKK6 has strong similarity to those of AmMKK6 (93.7%), CfMKK6 (91%), PxMKK6 (73.5%), BmMKK6 (73.1%), and AaMKK6 (70.7%). As shown in Figure 1A, AccMKK6 contains typical features of MKKs, including a conserved SxxxT motif and serine/threonine protein kinase catalytic (S-TKc) domain. Based on the amino acid sequences of MKKs from A. cerana cerana and other species, an NJ phylogenetic tree was constructed to identify the position of AccMKK6 in the evolutionary history of MKK proteins ( Figure 1B). The results showed that AccMKK6 lies within the distinct co-orthologous pattern of insect MKK6s and has a close evolutionary relationship with AmMKK3.

Developmental stage-specific and tissue-specific expression patterns of AccMKK6
To analyze the characteristics of AccMKK6, we used qPCR to detect the expression patterns of AccMKK6 in various tissues and in different developmental stages of honeybees. The total RNA was collected from the head, thorax, abdomen, epidermis and midgut. As shown in Figure 2A, AccMKK6 was expressed in all the tissues that we isolated. Among the tissues examined, the most well-defined expression was observed in the thorax, followed by the epidermis and head. We also analyzed AccMKK6 expression during different developmental stages and found that AccMKK6 was highly expressed in the pupal stages. As shown in Figure 2B, AccMKK6 was highest in the Pw phase. During the larvae stage, the expression levels in L4 were higher than in L1. The expression of AccMKK6 was lower in the adult stages than in the earlier stages.

Expression profiles of AccMKK6 under a variety of environmental stresses
Honeybees are inevitably affected by various external stimuli, including climate change, ecological deterioration, and pathogenic microbial invasion. Therefore, we simulated several abiotic stimuli (H 2 O 2 , heavy metals, and pesticides) to investigate the expression profile of AccMKK6, as examined by qPCR. Surprisingly, the relative expression of AccMKK6 was up-regulated in all of the tested treatments. As shown in Figure 2C, the expression of Ac-cMKK6 increased with H 2 O 2 treatment. HgCl 2 also up-regulated AccMKK6 expression by 1.6-fold at 1 h, and CdCl 2 up-regulated AccMKK6 expression by 3.8-fold at 0.5 h, after which expression returned to the basal level ( Figure  2D,E). Honeybee poisoning by pesticides is a serious problem in beekeeping. As shown in Figure 2F-H, AccMKK6 responded robustly to treatment with a variety of pesticides. The expression of AccMKK6 increased by 29-fold after 1.5 h of exposure to paraquat. Under exposure to emamectin benzoate, the relative expression level of AccMKK6 increased by 8.8-fold after 3 h. Moreover, the expression of AccMKK6 increased more drastically after exposure to acetamiprid. All of these conditions are believed to induce the formation of ROS. These data confirmed that AccMKK6 may play an important role in the resistance of honeybees to oxidative stress.

Protein expression levels of AccMKK6 under environmental stresses
To further explore the expression pattern of AccMKK6 in response to various types of oxidative damage, Western blot analysis was used to assess the changes in AccMKK6 after exposure to adverse stress. AccMKK6 was detected using anti-AccMKK6. Density analysis showed that increased AccMKK6 expression was induced after HgCl 2 treatment at 0.5 h ( Figure 3A). As shown in Figure 3B, CdCl 2 treatment also increased AccMKK6 protein expression levels. Treatment with paraquat and emamectin benzoate up-regulated AccMKK6, and the peak levels were observed at 1 and 0.5 h, respectively ( Figure 3C,D). The protein expression of AccMKK6 was reduced after acetamiprid treatment ( Figure 3E).

Effects of AccMKK6 knockdown on oxidative stress
RNAi mediated by dsRNA has been used to study gene function in recent years. Here, adult workers were injected with dsAccMKK6-RNA, dsGFP-RNA, or water. The expression of AccMKK6 was detected by qPCR and Western blot to determine the optimal silencing time. The transcription level of AccMKK6 in the dsAccMKK6-RNA injection group was significantly inhibited compared with those of the control groups, especially at 12, 24, and 36 h ( Figure 4A). In the control groups, the level of AccMKK6 mRNA fluctuated and then stabilized at 36 h after injection. To clarify the effect of AccMKK6 silencing at the protein level, we performed a Western blot. As shown in Figure 4B, the AccMKK6 expression in the dsAccMKK6-RNA injection group was lower than in the control groups. The above results indicate that dsRNA-mediated gene silencing is successful, and honeybees at 36 h after the dsRNA injection were selected for further study.
To further explore the role of AccMKK6 in combating oxidative damage, we analyzed the transcription levels of some antioxidant genes at 36 h after AccMKK6 silencing. The qPCR results in Figure 4C-H showed that Acc-SOD1, AccSOD2, AccGSTS4, AccCYP4G11, AccGSTO2, and AccGSTD were down-regulated when AccMKK6 was knocked down. To confirm the effect of RNAi with AccMKK6 on the antioxidant capacity of A. cerana cerana, we measured the enzymatic activities of SOD and POD and the levels of MDA 36 h after RNAi treatment. The results proved that the MDA concentration in the honeybees in which AccMKK6 was knocked down were higher than those in the control honeybees ( Figure 4I). Moreover, the activities of SOD and POD were lower in the silenced samples than in the control groups ( Figure 4J,K).

AccMKK6 interacted with Accp38b in vitro
Previous studies have reported that MKK6 mediates the activation of p38 MAPK in the response to oxidative stress [36]. To investigate the downstream component of AccMKK6, based on our previous research and the honeybee genome, we detected the interactions between AccMKK6 and Accp38b using yeast two-hybrid assays. The yeast two-hybrid assay results showed that the positive control clone and the clone co-transformed with AccMKK6 and Accp38b grew well on SD medium DDO and QDO plates ( Figure 5A). To verify the yeast two-hybrid results, we performed GST pull-down assays. As shown in Figure 5C, the AccMKK6 protein was captured and specifically bound to p38b protein, according to an assay using the antibody against the His epitope. These results indicated that AccMKK6 interacts with Accp38b.

Cloning and characterization of Acckayak
The transcription factor kayak is among the most important substrates for MAPK signal pathways [37]. Therefore, kayak was cloned to investigate the involvement of AccMKK6 in the MAPK cascades in A. cerana cerana. The ORF of Acckayak is 747 bp in length and encodes a putative polypeptide of 248 amino acid residues with a predicted molecular mass of 27.73 kDa and a theoretical pI of 6.8. Multiple sequence analysis of several kayaks from different species revealed that the deduced amino acid sequence of Acckayak is highly homologous to those of Amkayak, Btkayak, Oakayak, Tckayak, and Wakayak. Acckayak shared an average similarity of 70.59% with kayak sequences from other eukaryotic organisms, which suggested that kayak is highly conserved across species. As shown in Figure  6A, Acckayak contained a conserved basic leucine zipper protein (bZIP) domain. An NJ phylogenetic tree was constructed to understand the evolutionary relationships among kayaks ( Figure 6B). Phylogenetic analysis revealed that the homologous protein Amkayak shares the highest similarity with Acckayak.

Accp38b interacted with Acckayak, and both were induced by AccMKK6 silencing
To investigate the interaction between Accp38 and Acckayak, a yeast two-hybrid assay was performed. As shown in Figure 5B, Accp38b and Acckayak were co-transformed into the Y2H Gold yeast strain, which grew well on the QDO plates. The interaction between Accp38b and Acckayak was further confirmed using a GST pull-down assay. As shown in Figure 5D, the Acckayak protein was captured and specifically bound to p38b protein, according to the assay using the antibody against the His epitope. These results indicated that Accp38b interacts with Acckayak. In addition, the expression levels of Acckayak and Accp38b were examined after AccMKK6 silencing. Both Acckayak and Accp38b had higher expression levels in AccMKK6-silenced honeybees than in control honeybees ( Figure 5E,F).

Expression profiles of Acckayak under a variety of environmental stresses
To explore the relationship between AccMKK6 and Acckayak, we analyzed the expression patterns of Acckayak after several stress treatments. As shown in Figure 7, the expression of Acckayak was induced after the H 2 O 2 , HgCl 2 ,  CdCl 2 , and paraquat treatments. Notably, the expression patterns of Acckayak were approximately consistent with those of AccMKK6 after the same treatments.

Discussion
The MAPK signaling pathway is one of the important signal transduction systems in organisms, and many MKKs have been identified. MKKs are a family of serine/threonine protein kinases that are involved in many cellular responses including cell growth, nervous system development and immunity [17,38]. Previous studies demonstrated that MKKs may play a role in combatting oxidative stress [12]. Thus, we guessed that AccMKK6 might participate in the oxidative stress response and have conducted a series of experiments to test our prediction.
MKK6 is a member of the dual-specificity protein kinase family and is activated upon the phosphorylation of serine and threonine residues by upstream MKKKs [39]. The serine/threonine protein kinase catalytic domain (S-TKc) was found in the polypeptide sequence of AccMKK6 subject to protein phosphorylation (residues 48-309), and Ser 202 and Thr 206 sites were identified. Sequence analysis revealed that AccMKK6 contained sequences that are highly conserved among other typical MKK6s in insects. As shown in the phylogenetic tree, different MKKs from insects were clearly clustered into four clades, and AccMKK6 belongs to the MKK3/6 clade, sharing the highest degree of homology with the MKK3 protein from Apis mellifera.
Understanding the stage and tissue distributions of AccMKK6 mRNA may help to better understand its functional mechanisms. Therefore, we analyzed AccMKK6 mRNA expression at multiple developmental stages and examined its tissue distribution. The transcription level of AccMKK6 was strikingly higher in the pupal stages than in the other stages. Because ROS accumulation can cause oxidative damage, particularly in fast-growing organisms [40], we hypothesized that AccMKK6 is essential for honeybee development and may prevent ROS-mediated oxidative damage during the pupae growth stages. Tissue-specific expression analysis showed that AccMKK6 is expressed in various tissues but is most highly expressed in the thorax and epidermis. Similarly, the highest expression of Accp38b has been detected in the thorax [34]. The epidermis acts as the exoskeleton, conferring physical stability, and participates in the responses to various environmental stresses [41]. These data indicated that AccMKK6 might be involved in protecting honeybees against harm from oxidative stress.
Honeybees suffer constantly from a variety of adverse environmental stresses, and their number has plummeted [42]. H 2 O 2 is a typical oxidant that induces ROS elevation and cell death [43]. Heavy metal ions play essential roles in many physiological processes, and excess heavy metals can activate the MAPK cascades to interfere with cellular oxygen metabolism [44]. Moreover, Cd 2+ and Hg 2+ exposure enhances lipid peroxidation and ROS generation and alters the level of antioxidant enzymes [45,46]. Pesticides are significant threats to a honeybee's life and can leading to impaired biochemical and physiological functions [47]. For instance, emamectin benzoate is a novel macrocyclic lactone insecticide that can induce intracellular ROS accumulation and cell apoptosis [48]. Paraquat is one of the most widely used herbicides in agriculture and is a strong redox agent that yields paraquat monocation radicals, which then react with molecular oxygen to produce superoxide anions [49]. In A. mellifera, the neonicotinoid insecticide acetamiprid acts on acetylcholine nicotinic receptors and affects behavior [50]. After the exposure of Cydia pomonella to acetamiprid, the activity of CYP450 activity increased, suggesting activation of the antioxidant system to maintain redox homeostasis [51]. Therefore, we simulated several adverse life-threatening environmental conditions that could cause ROS damage in A. cerana cerana and explored the resulting changes in AccMKK6 expression at the transcription levels to look for evidence of their antioxidant functionality under environmental stress. We found that the expression level of AccMKK6 was up-regulated after H 2 O 2 , CdCl 2 and HgCl 2 treatment. Our results also revealed that the pesticide treatments affected the expression profiles of AccMKK6. As we all know, the gene expression pattern is usually an indicator of its function. Studies have shown that organisms repair oxidative damage by altering the expression of related genes, including by up-regulating antioxidant genes [52]. These results collectively indicated that AccMKK6 might protect honeybees from ROS-induced damage under oxidative stress.
We also detected the expression of AccMKK6 at the protein level after HgCl 2 , CdCl 2 , emamectin benzoate, and acetamiprid treatments. The results were not completely consistent with the mRNA expression data. The protein-to-mRNA ratio is different for different genes; additionally, the ratio might change after different treatments [53]. The stability of the transcript itself and the efficiency of translation and the stability of the translated product will affect the expression of the final product, as will problems with post-translational processing.
The RNA-mediated inhibition of endogenous target gene expression has become a popular strategy for determining gene function. An RNAi experiment was performed to further investigate the function of AccMKK6. AccMKK6 was successfully knocked down by RNAi. Organisms have many systems to control oxidative damage caused by ROS. Cellular antioxidant systems consist of multifunctional proteins and enzymes [54]. The expression of antioxidant genes changed after AccMKK6 knockdown, indicating that AccMKK6 may response to oxidative stress. By investigating the transcription levels of most of the antioxidant genes, such as AccSOD1, AccSOD2, AccGSTS4, AccCYP4G11, AccGSTO2, and AccGSTD, we found that they were reduced when AccMKK6 was knocked down. Previous studies reported that AccYB-1 may play an important role in response to oxidative damage [55]. Nrf can restrain the expression of antioxidant enzymes to maintain the activity of cellular defenses and/or to rapidly restore induced enzymes to normal levels [56]. We analyzed the transcription levels of AccYB-1 and AccNrf in AccMKK6-silenced honeybees and found that AccYB-1 and AccNrf were induced when AccMKK6 was silenced (Supplementary Figure S1). These results indicate that the AccMKK6 signal cascade plays a role in oxidative stress.
Lipid oxidation results in the generation of many harmful secondary products, and MDA can act as a reliable indicator of lipid peroxidation, indirectly reflecting cellular damage [57]. When AccMKK6 was silenced in honeybees, lipid peroxides began to accumulate faster, resulting in higher MDA concentrations. SOD and POD, as antioxidant enzymes, are the first line of defense that acts directly on ROS attacks [40]. SODs convert superoxide to oxygen into H 2 O 2 , and then PODs convert the H 2 O 2 into H 2 O [58]. The activity levels of these enzymes were changed by AccMKK6 knockdown, which may result in elevated ROS levels. Based on these findings, we conclude that the AccMKK6 may possess potent antioxidant properties.
Previous studies have shown that in mammals, MKK6 regulates a variety of fundamental cellular processes involved in the response to environmental stimuli by participating in the MAPK cascades [59,60]. In the p38 MAPK cascade, MKK3, MKK4, and MKK6 serve as upstream MAPK kinases responsible for p38 activation, and while downstream targets of p38 are either other kinases or transcription factors such as ATF-2 and MEF2 [14,61]. The qPCR results showed that the transcriptional factors ATF2 and MEF2 were responsive to environmental and chemical stresses (Supplementary Figure S2). These results suggest that MAPK cascade pathways play an important role in the response to various stresses in A. cerana cerana. To verify whether AccMKK6 participates in the MAPK cascades in A. cerana cerana, we cloned the transcription factor kayak from A. cerana cerana and named it Acckayak. Yeast two-hybrid analysis indicated that AccMKK6 interacts with Accp38b, and Accp38b interacts with Acckayak. At the same time, we also found that the expression levels of Accp38b and Acckayak were up-regulated after AccMKK6 was silenced, which supports the participation of AccMKK6 in the MAPK cascades in A. cerana cerana. A recent study reported CAR binding and inhibition of the GADD45B-MKK6 scaffold to repress the phosphorylation of p38 MAPK [62].