Shan Zhao1,2,3
Ruijuan Wang1,2,3Yan Liu1,2,3Long Su1,2,3Xiaoyan Dai1,2,3Dongyun Qin1,2Hao Chen1,2,3Zhenjuan Yin2,4Li Zheng1,2,3
Yifan Zhai1,2,3*- 1Institute of Plant Protection, Shandong Academy of Agricultural Sciences, Jinan, China
- 2Key Laboratory of Natural Enemies Insects, Ministry of Agriculture and Rural Affairs, Jinan, China
- 3China Ministry of Agriculture and Rural Affairs (MARA) - Centre for Agriculture and Bioscience International (CABI) Joint Laboratory for Bio-Safety Shandong Sub-Center, Jinan, China
- 4College of Agriculture, Guizhou University, Guiyang, China
Forkhead box O (FoxO), a key transcription factor in many species, participates in numerous physiological and pathological processes of organisms through a variety of signaling pathways. In the present study, we established DsFoxO knockout (DsFoxO-KO) strain using CRISPR/Cas9, and the influence on development and fecundity of mutant strain were evaluated. To clarify the corresponding mechanism, a transcriptome analysis was conducted subsequently. The results showed that the survival rates of the DsFoxO-KO strain in larval, pupal, and adult stages were all significantly lower than those of control. The duration of the pupal stage was similar between the two strains; however, durations of egg, larva, adult preoviposition period (APOP), and total APOP (TPOP) in the DsFoxO-KO strain were all significantly longer compared to those of the control strain. The fecundity of the DsFoxO-KO strain was 20.31 eggs/female, which was significantly lower than that of the control strain (430.47 eggs/female). With the transcriptome analysis, 612 differentially expressed genes (DEGs) were identified. Following COG and GO analyses, we found that most of the DEGs were associated with the metabolic process. According to the KEGG database, the mTOR signaling, MAPK signaling, Wnt signaling, and Toll and Imd signaling pathways; insect hormone biosynthesis; autophagy; and apoptosis were altered in the DsFoxO-KO strain. These results demonstrated that knockout of DsFoxO in D. suzukii significantly influenced its development and fecundity, while transcriptome analysis provided insights to explore the corresponding molecular mechanism. These findings highlighted the critical role of FoxO in D. suzukii and might contribute to the development of novel management strategies for these flies in the future.
Introduction
Drosophila suzukii, commonly known as spotted wing Drosophila, was first reported in Japan in 1931 as a cause of damage to fruit crops (Matsumura, 1931). Currently, the fly has spread across the world to more than 52 countries via commercial exchanges and tourism (Andreazza et al., 2017; Dos Santos et al., 2017; Little et al., 2017). D. suzukii is now considered a major agricultural pest species of thin-skinned fruits (Lee et al., 2011; Burrack et al., 2013; Mazzi et al., 2017). Unlike most Drosophila species that prefer to lay eggs on fermenting fruits, D. suzukii is a more serious threat as its special oviposition site selection (Atallah et al., 2014). D. suzukii females are able to pierce the skin of ripening healthy fruits with their enlarged ovipositor and deposit eggs inside the flesh (Poyet et al., 2015; Muto et al., 2018), causing evident damage before harvest. Additionally, other pathogens, such as bacteria, yeasts, and fungi, could infect the fruits via damaged skin, leading to further deterioration (Asplen et al., 2015). The primary and effective approach to control D. suzukii is the application of chemical insecticides; however, the emergence of insecticide resistance and the environmental unsustainability have hampered the control efficiency of chemical insecticides (Shaw et al., 2018; Shaw et al., 2019; Mastore et al., 2021). Thus, researchers are actively seeking new methods and searching for potential target genes that might provide insights to develop alternative control strategies (Alphey and Bonsall, 2018; Ni et al., 2021).
The forkhead box O (FoxO) protein, a highly conserved transcription factor, has been shown to participate in a variety of physiological and pathological processes, including longevity, growth, stress resistance, and metabolism (Greer et al., 2007; Puig and Mattila, 2011). Unlike mammals that have four kinds of FoxO genes (FoxO1, FoxO3a, FoxO4, and FoxO6), Drosophila has only one (dFoxO) (Kaestner et al., 2000; de la Torre-Ubieta et al., 2010). The DNA-binding domain of dFoxO, located in the N-terminal, has 45% identity with FoxO4, whereas 84% identity in the three α-helix regions. In the C-terminal, dFoxO has a higher Ser and Gln content than mammals (Puig et al., 2003). In Drosophila melanogaster, activated FoxO had been shown to extend the lifespan of the organisms (Alic et al., 2014) and regulate RNA interference, providing protection against RNA virus infection (Spellberg and Marr, 2015). In mammals, FoxO proteins are widely expressed in the brain, where they could mediate both neuroprotection and neurodegeneration (Yuan et al., 2009; Siegrist et al., 2010). Numerous studies have also been conducted on mice, and the results showed that constant expression of FoxO proteins would regulate the activity of follicular, suggesting that the proteins are also involved in fecundity (Pelosi et al., 2013; Barilovits et al., 2014). Similar functions of FoxO were also found in Caenorhabditis elegans (Michaelson et al., 2010; Hibshman et al., 2016). Given that FoxO plays such important roles in a variety of organisms, it emerges as a potential candidate for targeting genes in pest control strategies.
A detailed molecular mechanism of FoxO has been gradually uncovered through the dedicated efforts of researchers, and several signaling pathways have been proved to be involved in its function. Among these pathways, insulin/insulin-like growth factor (insulin/IGF) signaling seems to be studied the most for FoxO functions, and through that, FoxO has been found to regulate a number of target genes involved in metabolism and cell progression (Tatar et al., 2001; Yamamoto and Tatar, 2011; Martins et al., 2016). Nam et al. (2022) revealed that translational controlled tumor protein (Tctp) regulates cell growth by reducing cytoplasmic FoxO levels in Drosophila. In Tribolium castaneum, Blattella germanica, and D. melanogaster, the effects of FoxO on fecundity have been linked to the downregulation of the expression of vitellogenin (Sheng et al., 2011; Süren-Castillo et al., 2012; Yang et al., 2013). After a deep research on Helicoverpa armigera, Zhang et al. (2022) suggested that pupal diapause is induced by the decreases of transforming growth factor-beta (TGFβ), which was stimulated by the FoxO-activated ubiquitin-proteasome system (UPS). All the aforementioned results demonstrate that FoxO is involved in different signaling pathways in different species, and further studies are still required to uncover the specific pathways in which FoxO participates in D. suzukii.
In the present study, to explore the influence of FoxO on development and fecundity of D. suzukii, we knocked out the DsFoxO gene through a 17-bp deletion in this species using the method of CRISPR/Cas9. Subsequent evaluation revealed that the development and fecundity of flies were largely influenced in this knockout strain. To explore the corresponding molecular mechanism, transcriptome analysis was conducted to screen the related genes and signaling pathways. The findings of this study provide further evidence of FoxO’s influence on D. suzukii and offer potential insights for the development of sustainable strategies in the management of this pest species.
Materials and methods
Insect rearing
The D. suzukii colony was originally collected from cherry orchards in Tai’an (Shandong Province, China) in 2012 and identified as the TA2012 strain. The flies were maintained in our laboratory in a climate-controlled incubator at 25°C ± 1°C under the condition of the 16L:8D daylight cycle and 70%–80% relative humidity (Zhai et al., 2016).
Preparation of sgRNA
The single guide (sg)RNA target site of DsFoxO (XM_036818322.1) was identified with the principle of 5′-GGN18NGG-3’ (the protospacer-adjacent motif (PAM) sequence is marked as bold). Gene structures are sketched in Figure 1A. The sgRNA used in this study was designed as 5′-GGCGACCTACCCCTGGACGTGGG-3’. The sgRNA template was synthesized with two oligonucleotides, with one containing the T7 promoter and target sequence (5′-TAATACGACTCACTATAGGN18 GTTTTAGAGCTAGAAATAGC-3′) and the other was the universal oligonucleotide encoding the remaining sgRNA sequences (5′-AAAAGCACCGACTCGGTGCCACTTTTTCAAGTTGATAACGGACTAGCCTTATTTTAACTTGCTATTTCTAGCTCTAAAAC-3′). PCR, as well as the PCR product purification, was performed as described previously (Wang et al., 2016). Using purified DNA as templates, sgRNA was synthesized with the GeneArt™ Precision gRNA Synthesis Kit (Thermo Fisher Scientific, Vilnius, Lithuania), according to the manufacturer’s instruction.
FIGURE 1. CRISPR/Cas9-mediated knockout of DsFoxO in D. suzukii. (A) Genomic structure of DsFoxO including protein-coding regions (black polygon) and untranslated regions (white polygon), and the position for sgRNA (red triangle). (B) Alignment analysis of PCR products from the wild-type (WT) strain and mutant (Mut) strain. Target sequences of sgRNA are in green with bold, while the PAM sequences in red with bold. Lowercase letters represent inserted bases. The amino acids of FoxO are in blue, while mutant amino acids are in brown. * indicates the terminator. The numbers represent the location of the base in the FoxO gene.
Embryo microinjection
Eggs, produced within half an hour, were collected and lined up on a microscope slide. Using the FemtoJet and InjectMan NI2 Microinjection system (Eppendorf, Hamburg, Germany), each egg was injected with 1 nl mixture of sgRNA (500 ng/μl) and Cas9 protein (250 ng/μl, Thermo Fisher, Shanghai, China). A total of 400 eggs were implemented with injection; after that, injected eggs were maintained in a climate-controlled incubator, as described previously.
Identification of the DsFoxO-KO strain
Adults that hatched from injected eggs or resulted from the following generations were collected, and their genomic DNA was extracted using the AxyPrep™ DNA Extraction Kit (Axygen Biosciences, Union City, CA, United States) for genotyping. Using the extracted DNA as a template, PCR was performed with the primer pair (F: 5′-ATGATGGACGGCTTCGCGCAGGACT-3′ and R: 5′-CTCTGTGGATGAGAACCTACCTC-3′). PCR products were sequenced by Tsingke Biological Company (Qingdao, China) to detect the indel mutation on the target site.
Detection of life history traits
A total of 270 eggs of the DsFoxO-KO strain were collected and reared individually in Petri dishes (3.5 cm in diameter). The dishes were divided into three groups with each group containing 90 dishes. The duration of the egg incubation period, and larval and pupal stages were recorded daily, and the stage-specific survival rates (i.e., the proportion of individuals that survived from the previous stage) were calculated. Immediately after adult emergence, 15 male–female pairs were placed individually for mating in bottles. The bottles were divided into three groups with each group containing five bottles. The oviposition was observed over a period of 30 days, and the adult preoviposition period (APOP), total APOP (TPOP), and fecundity were recorded and calculated. The TA2012 strain was implemented with the same treatments and used as control.
RNA sequencing
Total RNA of six adult females (15 days post-emergence) was extracted using TRIzol (Vazyme, Nanjing, China), following manufacturer’s instruction, and digested with DNase I (Vazyme, Nanjing, China). The mRNA with polyA tails was enriched with oligo (dT) magnetic beads and dealt with fragmentation buffer. cDNA was synthesized with random primers using mRNA as a template. After purification, the cDNA was repaired and A bases were added to 3′ ends. The cDNA sequences were screened with AMPure XP beads, and the cDNA library was obtained through PCR. After the library was confirmed as qualified, sequencing was performed with the Illumina NovaSeq6000 platform (BioMarker Technologies, Qingdao, China). Three biological replicates were used for RNA sequencing, and the TA2012 strain of D. suzukii was used as control.
qRT-PCR
Six genes were selected randomly to verify the results of RNA sequencing. Total RNA of adults was extracted and treated with DNase I. cDNA was synthesized using the PrimeScript™ IV 1st strand cDNA Synthesis Mix (TaKaRa, Dalian, China) in a 20 μl reaction, containing 5×Mix 4 μl, Random 6-mers 2 μl, RNA 5 μg, and an appropriate amount of RNase-free dH2O, with a procedure of 30°C 10 min, 42°C 20 min, and 95°C 5 min. qRT-PCR was carried out using the 2×SYBR Green qPCR Mix (SparkJade, Qingdao, China) in a 20 μl reaction containing 2×SYBR qPCR Mix 10 μl, cDNA 2 μl, forward primer (10 μM) 0.4 μl, reverse primer (10 μM) 0.4 μl, ROX reference dye II 0.4 μl, and RNase-free H2O 6.8 μl, with a procedure of 95°C 20 s, 40 cycles of 95°C 3 s, and 60°C 30 s, and a dissociation stage of 95°C 15 s and 60°C 1 min. qRT-PCR was performed using the 7500 Fast Real-Time PCR System (ABI). The relative expression levels of genes were calculated using the 2−ΔΔCt method, and DsActin was used as a comparator with the primer pair (F: 5′-CTACGAGGGTTATGCCCTGC-3′ and R: 5′-CGGTGGTGGTGAACGAGTAA-3’).
Statistical analysis
Data were analyzed by IBM SPSS Statistics 24, followed by separation of means ± SE using Student’s t-test. Differences were significant at p < 0.05.
Results
Establishment of the DsFoxO-KO strain of D. suzukii
To knockout DsFoxO, 400 eggs of the TA2012 strain were injected with sgRNA and Cas9, and 90 G0 embryos survived to adulthood. Of these, 72 were divided into eight groups for mating with each group containing six females and three males. After G1 eggs were collected, we genotyped the G0 adults for the expected CRISPR-mediated indel mutation with PCR and confirmed by direct sequencing of the PCR products. The results showed that 38 of these 72 adults were edited. Then, 20 G1 adults were collected and a single pair was carried out with the TA2012 strain. After G2 eggs were collected, the G1 adults were genotyped and the results showed that four adults were mutated. From the four mutants, we selected a heterozygous mutant with a 17-bp deletion (19 bp deleted and 2 bp inserted, Figure 1B) for further study. From the single-pair G1 family that had this one parent heterozygous, we obtained 65 adults (G2). G2 adults were genotyped with 1/4 wings, and the results showed that, among the 65 adults, seven females and six males were heterozygous. The 13 adults were pooled for mating to produce G3. G3 adults were also genotyped with 1/4 wings, and two females and one male were homozygous. These three adults were pooled to generate the DsFoxO-KO strain.
Life history traits of the DsFoxO-KO strain
The effects of FoxO on life history traits in D. suzukii were detected in this study. The results showed that the survival rates of the DsFoxO-KO strain in the stages of larva, pupa, and adult were 54.07%, 63.69%, and 63.63%, respectively, which were all significantly lower than those in the TA2012 strain displayed as 91.33% (larva), 76.94% (pupa), and 95.27% (adult) (p < 0.05) (Figure 2).
FIGURE 2. Survival rate of different life history stages of DsFoxO-KO and TA2012 strains in different stages. A total of 270 eggs were used and divided into three groups in this experiment. The data were shown as means ± SE. *** represents p < 0.005, and **** represents p < 0.001.
The developmental periods of each stage were also detected in our study. The duration of the pupa stage was similar in the two strains (p > 0.05), whereas the egg and larval stages of the DsFoxO-KO strain (1.53 d and 6.07 d) were significantly longer than those in the TA2012 strain (1.17 d and 4.78 d) (p < 0.05). APOP and TPOP in the DsFoxO-KO strain were 12.08 d and 24.44 d, respectively, which were significantly longer than those in the TA2012 strain (3.47 d and 14.60 d) (p < 0.05). In 30 days after pairing, the fecundity in the DsFoxO-KO strain was 20.31 eggs/female, which was significantly lower than that in the TA2012 strain (430.47 eggs/female) (p < 0.05) (Figure 3).
FIGURE 3. Comparison of life stage and fecundity parameters of DsFoxO-KO and TA2012 strains. The data were shown as means ± SE. *** represents p < 0.005, and **** represents p < 0.001.
Differentially expressed genes (DEGs) resulting from DsFoxO knockout
Six samples (three replicates in each group) were used for RNA sequencing analysis, and low-quality and adapter sequences were removed from raw data using SOAPnuke software (v1.4.0) and Trimmomatic (v0.36). About 46.22 Gb clean data were obtained, with each sample containing clean data over 5.75 Gb. The clean reads in the DsFoxO-KO groups were 28882942, 26731905, and 29757666, respectively, and the percentages of GC contents were 54.30%, 54.39%, and 54.22%, respectively, while in the TA2012 strain, the clean reads were 20892031, 19266328, and 29112744 with the percentages of GC contents of 52.78%, 53.33%, and 54.62%, respectively. The percentages of bases with quality value ≥30 were all over 91.87% (Table 1).
TABLE 1. Statistics of RNA sequencing data.
The expression levels of transcripts were quantified using the RSEM package (Bo and Dewey, 2011) and analyzed based on read counts with DESeq software (Love et al., 2014). The fold change (FC) ≥ 2 and false discovery rate (FDR) < 0.01 were used as standards, and FDR was classified by the p-value. In total, compared to the TA2012 group, 612 DEGs were identified, of which 404 were upregulated and 208 were downregulated (Figure 4A).
FIGURE 4. Transcriptome analysis of DsFoxO-KO and TA2012 strains. (A) Volcano plot of DEGs. Each dot represents a gene. The X-axis shows log2 (FC), and the Y-axis shows −log10 (FDR). All dots were filtered by FC ≥ 2 and FDR < 0.01. The red dots indicate upregulated DEGs. The blue dots indicate downregulated DEGs. The gray dots indicate normal genes. (B) Gene function classification based on the COG annotation. (C) Gene function classification according to GO assignments.
DEG analysis
For functional annotation, all DEGs were analyzed with the methods of cluster of orthologous groups of proteins (COG) (Tatusov et al., 2000), the Gene Ontology (GO) consortium (Ye et al., 2006), and Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Kanehisa et al., 2008).
The results of COG analysis showed that 248 sequences were classified to 26 categories with seven categories containing no sequences. Among these categories, the cluster for “Posttranslational modification, protein turnover, and chaperones” constituted the largest group (48, 19.35%), followed by “Carbohydrate transport and metabolism” with a number of 32 (12.90%) and “Secondary metabolites biosynthesis, transport, and catabolism” (24, 9.68%). On the contrary, “Cell cycle control, cell division, and chromosome partitioning” (1, 4.03%) and “Coenzyme transport and metabolism” (1, 4.03%) represented the smallest groups (Figure 4B).
The DEG enrichment was classified into 32 different groups, according to GO assignments, including three categories of biological process (20 groups), cellular component (three groups), and molecular function (nine groups). The biological process category accounted for the largest number of DEGs, among which the cellular process (211 DEGs) and metabolic process (196 DEGs) were the two most abundant terms. The cellular anatomical entity (264 DEGs) was the most abundant term in the cellular component category. In the molecular function category, catalytic activity (213 DEGs) and binding (180 DEGs) were the two most abundant terms (Figure 4C).
The KOBAS database and clusterProfiler software were used to test the statistical enrichment of DEGs in KEGG pathways (Mao et al., 2005). In addition, 350 annotated genes were analyzed according to the KEGG database, and six biological pathways were identified, including cellular process, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems. Of these, the pathways represented by most genes were metabolism (155, 44.29%), followed by environmental information processing (60, 17.14%) (Figure 5).
FIGURE 5. Biological pathways according to the KEGG database. The X-axis represents the percentage of genes. The Y-axis represents KEGG terms. Numbers on the right margin of each bar represent numbers of genes in the corresponding term.
In addition, the mTOR signaling pathway (six DEGs), MAPK signaling pathway (three DEGs), Wnt signaling pathway (five DEGs), Toll and Imd signaling pathways (14 DEGs), insect hormone biosynthesis (six DEGs), autophagy (five DEGs), and apoptosis (two DEGs) were all identified based on the KEGG database (Figure 5). The key genes in these pathways were analyzed, and we found that most genes involved in insect hormone biosynthesis were upregulated in the DsFoxO-KO strain, whereas genes in other pathways were mainly downregulated (Figure 6).
FIGURE 6. Expression heatmap of key genes in several pathways identified by KEGG. Red and blue indicate upregulated and downregulated genes relative to mean expression levels of the TA2012 strain, respectively, and the ratios were plotted on a color scale. Gene expression levels are represented by fragments per kilobase of transcript per million mapped reads (FPKM) and analyzed with log2 and Z-score.
Validation with qRT-PCR
To validate the results of RNA sequencing, six genes were selected for qRT-PCR. The information about the genes and the primers is listed in Table 2. The results showed that the expression levels of LOC108012461, LOC108008525, LOC108013851, and LOC108013219 were upregulated, whereas those of LOC108008512 and LOC108020151 were downregulated compared with each gene in the TA2012 strain (Figure 7). The expression levels of all six genes were consistent with the results of RNA sequencing.
TABLE 2. Sequences of the primers used for the qRT-PCR analysis of six randomly selected genes.
FIGURE 7. Expression level detection of six genes with qRT-PCR. The results were calculated using DsActin as control and were expressed as mean ± SE. * represents p < 0.05, and **** represents p < 0.001.
Discussion
As a key transcription factor, the functions of FoxO have been well studied in various species (Bai et al., 2013; Riedel et al., 2013). In Drosophila, FoxO has been demonstrated to play a role in metabolic adaptation, neuromuscular junction homeostasis, cytoskeletal dynamics, and microtubule functioning (Nechipurenko and Broihier, 2012; McLaughlin et al., 2016; Molaei et al., 2019; Birnbaum et al., 2021). In the present study, we aimed to explore the potential effects of FoxO on development and fecundity of D. suzukii. To obtain reliable experimental data and materials for further study, we knocked out the DsFoxO gene of D. suzukii using the method of CRISPR/Cas9. After selection for three generations, the mutant strain was genotyped as homozygous with 17 bp deletion, which changed the amino acid sequences from the indel site and terminated the mRNA at the location of 400 bp approximately (Figure 1B), indicating that the protein was knocked out completely. The establishment of the DsFoxO-KO strain laid the foundation for subsequent research.
FoxO has been previously suggested to affect the reproduction of female Drosophila (Yang et al., 2013). In our study, we assessed the reproduction of the DsFoxO-KO strain. In the DsFoxO-KO strain, fecundity was significantly reduced from 430.47 eggs/female in the TA2012 strain to 20.80 eggs/female, indicating that the reproduction of D. suzukii was influenced by FoxO. Similarly, FoxO was also reported to affect the number of eggs laid by female Aedes aegypti and Nilaparvata lugens (Hansen et al., 2007; Dong et al., 2021). In addition, the survival rates of DsFoxO-KO larvae, pupae, and adults were all significantly lower than those in the TA2012 strain. The duration of different life stages was also tested, and we found that the duration of egg, larval, APOP, and TPOP stages in the DsFoxO-KO strain were all significantly longer than those in TA2012 strain. These results suggested that FoxO would perform a huge impact on the development of D. suzukii, which were consistent with previous research studies (Alic et al., 2014; Dobson et al., 2019; Shi et al., 2019). Therefore, knocking out DsFoxO leads to a significant retardation in the development and fecundity of D. suzukii. These findings highlight the potential functions of FoxO in D. suzukii.
Cell homeostasis, which is determined by the metabolic steady state, is quite important for the development and fecundity of organisms. In Drosophila, as well as in mammals, FoxO can influence growth by regulating the cell cycle and energy metabolism (Jünger et al., 2003; Greer et al., 2007). Moreover, in differentiating cells, non-normal expression of FoxO would lead to cellular atrophy and promote a catabolic state (Accili and Arden, 2004). In the present study, after the transcriptome analysis was implemented, the data were analyzed with the methods of GO and KEGG. The results of GO analysis showed that DEGs were enriched in the cellular process, metabolic process, and catalytic activity, whereas most genes were presented in the metabolism pathway with KEGG analysis, suggesting that normal cellular and metabolic processes might be affected when FoxO was knocked out. These results were consistent with previous studies and suggested that knockout of FoxO in Drosophila might disrupt the cellular homeostasis, and thus, the development and fecundity of the flies were influenced. The development and fecundity of Drosophila would also be affected by the utilization of nutrition directly or indirectly as the energy for the organism was generated by various forms of nutrients (Ratnaparkhi and Sudhakaran, 2022). In COG analysis, numerous DEGs were classified into the category of “Carbohydrate transport and metabolism,” and these genes might restrict the availability of nutrients, which would also result in the alteration of development and fecundity of the flies.
Several pathways in which FoxO might be involved were also identified based on DEG analysis, including the mTOR signaling pathway, MAPK signaling pathway, Wnt signaling pathway, Toll and Imd signaling pathways, insect hormone biosynthesis, autophagy, and apoptosis (Figure 5). These aforementioned pathways all have been suggested to play important roles in insect growth and development (Zeng et al., 2017; Zhang and Zhang, 2019; Hughson et al., 2021). For example, alpha-ketoglutarate (AKG) would extend Drosophila lifespan by inhibiting the mTOR pathway (Su et al., 2019), whereas the development of Drosophila larvae was altered through the MAPK signaling pathway at the condition of ancestral dietary change (Towarnicki et al., 2022). Particularly, various DEGs were associated with the Toll and Imd signaling pathway, which is an important pathway in the innate immune response. DEGs in this pathway mainly encoded peptidoglycan recognition proteins (PGRPs) and nuclear factor kappa-B (NF-κB) proteins, which are the two key signaling molecules of the Toll and Imd pathways (Kawai and Akira, 2007; Royet and Dziarski, 2007), and both have also been proposed to play important roles in development processes (Maillet et al., 2008; Morgan and Liu, 2011). With a deep exploration of previous reports, we found that these pathways also regulated several physiological processes, including cell cycle, cell proliferation, apoptosis, and autophagy (Xing et al., 2015; Zeng et al., 2017; Maiese, 2021; Yan et al., 2022), all of which could influence the development and fecundity of the organism.
In terms of insect hormone biosynthesis classification, we observed that the expression levels of the ecdysone 20-monooxygenase isoform (E20-MO), juvenile hormone acid O-methyltransferase (JHAMT), and cytochrome P450 were all upregulated in the DsFoxO-KO strain. From the previous reports, we learned that the ecdysone biosynthesis mediated the Drosophila body size through a FoxO-Ultraspiracle interaction in case of nutritional control (Koyama et al., 2014). Similarly, Mirth et al. (2014) indicated that the effects of juvenile hormone (JH) on the growth rate in Drosophila are mainly dependent on FoxO. Although few studies have been reported on the relationship between FoxO and cytochrome P450, this protein had been reported to be essential for the synthesis and degradation of ecdysone and JH in insects (Bergé et al., 1998). Based on the results of our study, we hypothesized that, in D. suzukii, these hormones might impact the development and fecundity via their interactions with FoxO, and in the DsFoxO-KO strain, the growth of the flies were retarded to some extent, which warranted further investigation in future studies.
Based on the aforementioned information, we hypothesized that the functions of some hormones probably were limited without FoxO and that FoxO might influence the development and fecundity directly or via one or multiple pathways with the alteration of normal physiological processes. Perhaps, some of the pathways might modulate the development and fecundity by the reaction with FoxO. This hypothesis would be explored in future mechanistic studies. On the other hand, FoxO has also been implicated in diapause regulation in various species, and insulin signaling seems to be a key developmental pathway involved in this regulation (Tatar and Yin, 2001; Sim and Denlinger, 2013; Sim et al., 2015; Zhang et al., 2022). In our study, we also discovered that the expression level of insulin-like peptide 2 was upregulated in the DsFoxO-KO strain. Thus, further work is required on the relationship between diapause and FoxO in Drosophila.
In conclusion, we successfully established a DsFoxO-KO strain of D. suzukii using CRISPR/Cas9 in the present study. Life history trait analysis indicated that the development and fecundity were significantly influenced without FoxO. Through RNA sequencing, several DEGs were identified as their expression levels change or classification of pathways. Further studies with these DEGs might elucidate the molecular mechanism underlying the influence of FoxO on the development and fecundity of D. suzukii. Such insights will be useful for the development of strategies to control this economically important pest species.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA1013261.
Author contributions
SZ: data curation, formal analysis, investigation, and writing–original draft. RW: data curation and writing–review and editing. YL: formal analysis and writing–review and editing. LS: data curation, investigation, and writing–review and editing. XD: formal analysis and writing–review and editing. DQ: methodology and writing–review and editing. HC: conceptualization, formal analysis, and writing–review and editing. ZY: methodology and writing–review and editing. LZ: conceptualization, funding acquisition, and writing–review and editing. YZ: conceptualization, funding acquisition, supervision, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was funded by the National Natural Science Foundation of China (31972273 and 32202313) and the Shandong Provincial Natural Science Foundation (ZR2021YQ21 and ZR2021QC218).
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.
References
Accili D., Arden K. C. (2004). Foxos at the crossroads of cellular metabolism, differentiation, and transformation. Cell 117, 421–426. doi:10.1016/s0092-8674(04)00452-0
Alic N., Giannakou M. E., Papatheodorou I., Hoddinott M. P., Andrews T. D., Bolukbasi E., et al. (2014). Interplay of dFOXO and two ETS-family transcription factors determines lifespan in Drosophila melanogaster. PLoS Genet. 10, e1004619. doi:10.1371/journal.pgen.1004619
Alphey N., Bonsall M. B. (2018). Genetics-based methods for agricultural insect pest management. Agric. For. Entomology 20, 131–140. doi:10.1111/afe.12241
Andreazza F., Bernardi D., Dos Santos R. S. S., Garcia F. R. M., Oliveira E. E., Botton M., et al. (2017). Drosophila suzukii in southern neotropical region: current status and future perspectives. Neotrop. Entomol. 46, 591–605. doi:10.1007/s13744-017-0554-7
Asplen M. K., Anfora G., Biondi A., Choi D. S., Chu D., Daane K. M., et al. (2015). Invasion biology of spotted wing Drosophila (Drosophila suzukii): a global perspective and future priorities. J. Pest. Sci. 88, 469–494. doi:10.1007/s10340-015-0681-z
Atallah J., Teixeira L., Salazar R., Zaragoza G., Kopp A. (2014). The making of a pest: the evolution of a fruit-penetrating ovipositor in Drosophila suzukii and related species. Proc. Biol. Sci. 281, 20132840. doi:10.1098/rspb.2013.2840
Bai H., Kang P., Hernandez A. M., Tatar M. (2013). Activin signaling targeted by insulin/dFOXO regulates aging and muscle proteostasis in Drosophila. PLoS Genet. 9, e1003941. doi:10.1371/journal.pgen.1003941
Barilovits S. J., Newsom K. J., Bickford J. S., Beachy D. E., Rhoton-Vlasak A., Nick H. S. (2014). Characterization of a mechanism to inhibit ovarian follicle activation. Fertil. Steril. 101, 1450–1457. doi:10.1016/j.fertnstert.2014.01.025
Bergé J. B., Feyereisen R., Amichot M. (1998). Cytochrome P450 monooxygenases and insecticide resistance in insects. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 353, 1701–1705. doi:10.1098/rstb.1998.0321
Birnbaum A., Sodders M., Bouska M., Chang K., Kang P., McNeill E., et al. (2021). FOXO regulates neuromuscular junction homeostasis during Drosophila aging. Front. Aging Neurosci. 12, 567861. doi:10.3389/fnagi.2020.567861
Bo L., Dewey C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf 12, 323. doi:10.1186/1471-2105-12-323
Burrack H. J., Fernandez G. E., Spivey T., Kraus D. A. (2013). Variation in selection and utilization of host crops in the field and laboratory by Drosophila suzukii matsumara (Diptera: drosophilidae), an invasive frugivore. Pest Manag. Sci. 69, 1173–1180. doi:10.1002/ps.3489
de la Torre-Ubieta L., Gaudillière B., Yang Y., Ikeuchi Y., Yamada T., DiBacco S., et al. (2010). A FOXO-Pak1 transcriptional pathway controls neuronal polarity. Genes Dev. 24, 799–813. doi:10.1101/gad.1880510
Dobson A. J., Boulton-McDonald R., Houchou L., Svermova T., Ren Z., Subrini J., et al. (2019). Longevity is determined by ETS transcription factors in multiple tissues and diverse species. PLoS Genet. 15, e1008212. doi:10.1371/journal.pgen.1008212
Dong Y., Chen W., Kang K., Pang R., Dong Y., Liu K., et al. (2021). FoxO directly regulates the expression of TOR/S6K and vitellogenin to modulate the fecundity of the brown planthopper. Sci. China Life Sci. 64, 133–143. doi:10.1007/s11427-019-1734-6
Dos Santos L. A., Mendes M. F., Krüger A. P., Blauth M. L., Gottschalk M. S., Garcia F. R. (2017). Global potential distribution of Drosophila suzukii (Diptera, Drosophilidae). PLoS One 12, e0174318. doi:10.1371/journal.pone.0174318
Greer E. L., Dowlatshahi D., Banko M. R., Villen J., Hoang K., Blanchard D., et al. (2007). An AMPK-FOXO pathway mediates longevity induced by a novel method of dietary restriction in C. elegans. Curr. Biol. 17, 1646–1656. doi:10.1016/j.cub.2007.08.047
Hansen I. A., Sieglaff D. H., Munro J. B., Shiao S. H., Cruz J., Lee I. W., et al. (2007). Forkhead transcription factors regulate mosquito reproduction. Insect biochem. Mol. Biol. 37, 985–997. doi:10.1016/j.ibmb.2007.05.008
Hibshman J. D., Hung A., Baugh L. R. (2016). Maternal diet and insulin-like signaling control intergenerational plasticity of progeny size and starvation resistance. PLoS Genet. 12, e1006396. doi:10.1371/journal.pgen.1006396
Hughson B. N., Shimell M., O'Connor M. B. (2021). AKH Signaling in D. melanogaster alters larval development in a nutrient-dependent manner that influences adult metabolism. Front. Physiol. 12, 619219. doi:10.3389/fphys.2021.619219
Jünger M. A., Rintelen F., Stocker H., Wasserman J. D., Végh M., Radimerski T., et al. (2003). The Drosophila forkhead transcription factor FOXO mediates the reduction in cell number associated with reduced insulin signaling. J. Biol. 2, 20. doi:10.1186/1475-4924-2-20
Kaestner K. H., Knochel W., Martinez D. E. (2000). Unified nomenclature for the winged helix/forkhead transcription factors. Genes Dev. 14, 142–146. doi:10.1101/gad.14.2.142
Kanehisa M., Araki M., Goto S., Hattori M., Hirakawa M., Itoh M., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, D480–D484. doi:10.1093/nar/gkm882
Kawai T., Akira S. (2007). Signaling to NF-kappaB by toll-like receptors. Trends Mol. Med. 13, 460–469. doi:10.1016/j.molmed.2007.09.002
Koyama T., Rodrigues M. A., Athanasiadis A., Shingleton A. W., Mirth C. K. (2014). Nutritional control of body size through FoxO-Ultraspiracle mediated ecdysone biosynthesis. elife 3, e03091. doi:10.7554/eLife.03091
Lee J. C., Bruck D. J., Curry H., Edwards D., Haviland D. R., Van Steenwyk R. A., et al. (2011). The susceptibility of small fruits and cherries to the spotted-wing Drosophila, Drosophila suzukii. Pest Manag. Sci. 67, 1358–1367. doi:10.1002/ps.2225
Little C. M., Chapman T. W., Moreau D. L., Hillier N. K. (2017). Susceptibility of selected boreal fruits and berries to the invasive pest Drosophila suzukii (Diptera: drosophilidae). Pest Manag. Sci. 73, 160–166. doi:10.1002/ps.4366
Love M. I., Huber W., Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Maiese K. (2021). Targeting the core of neurodegeneration: FoxO, mTOR, and SIRT1. Neural Regen. Res. 16, 448–455. doi:10.4103/1673-5374.291382
Maillet F., Bischoff V., Vignal C., Hoffmann J., Royet J. (2008). The Drosophila peptidoglycan recognition protein PGRP-LF blocks PGRP-LC and IMD/JNK pathway activation. Cell Host Microbe 3, 293–303. doi:10.1016/j.chom.2008.04.002
Mao X., Cai T., Olyarchuk J. G., Wei L. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi:10.1093/bioinformatics/bti430
Martins R., Lithgow G. J., Link W. (2016). Long live FOXO: unraveling the role of FOXO proteins in aging and longevity. Aging Cell 15, 196–207. doi:10.1111/acel.12427
Mastore M., Caramella S., Quadroni S., Brivio M. F. (2021). Drosophila suzukii susceptibility to the oral administration of Bacillus thuringiensis, Xenorhabdus nematophila and its secondary metabolites. Insects 12, 635. doi:10.3390/insects12070635
Mazzi D., Bravin E., Meraner M., Finger R., Kuske S. (2017). Economic impact of the introduction and establishment of Drosophila suzukii on sweet cherry production in Switzerland. Insects 8, 18. doi:10.3390/insects8010018
McLaughlin C. N., Nechipurenko I. V., Liu N., Broihier H. T. (2016). A toll receptor-foxo pathway represses Pavarotti/MKLP1 to promote microtubule dynamics in motoneurons. J. Cell Biol. 214, 459–474. doi:10.1083/jcb.201601014
Michaelson D., Korta D. Z., Capua Y., Hubbard E. J. A. (2010). Insulin signaling promotes germline proliferation in C. elegans. Development 137, 671–680. doi:10.1242/dev.042523
Mirth C. K., Tang H. Y., Makohon-Moore S. C., Salhadar S., Gokhale R. H., Warner R. D., et al. (2014). Juvenile hormone regulates body size and perturbs insulin signaling in Drosophila. Proc. Natl. Acad. Sci. USA. 111, 7018–7023. doi:10.1073/pnas.1313058111
Molaei M., Vandehoef C., Karpac J. (2019). NF-κB shapes metabolic adaptation by attenuating Foxo-mediated lipolysis in Drosophila. Dev. Cell 49, 802–810. doi:10.1016/j.devcel.2019.04.009
Morgan M. J., Liu Z. G. (2011). Crosstalk of reactive oxygen species and NF-κB signaling. Cell Res. 21, 103–115. doi:10.1038/cr.2010.178
Muto L., Kamimura Y., Tanaka K. M., Takahashi A. (2018). An innovative ovipositor for niche exploitation impacts genital coevolution between sexes in a fruit-damaging Drosophila. Proc. Biol. Sci. 285, 20181635. doi:10.1098/rspb.2018.1635
Nam S., Le T. P., Chung S., Choi K. W. (2022). Tctp regulates the level and localization of Foxo for cell growth in Drosophila. Cell Death Discov. 8, 146. doi:10.1038/s41420-022-00937-2
Nechipurenko I. V., Broihier H. T. (2012). Foxo limits microtubule stability and is itself negatively regulated by microtubule disruption. J. Cell Biol. 196, 345–362. doi:10.1083/jcb.201105154
Ni X. Y., Lu W. J., Qiao X., Huang J. (2021). Genome editing efficiency of four Drosophila suzukii endogenous U6 promoters. Insect Mol. Biol. 30, 420–426. doi:10.1111/imb.12707
Pelosi E., Omari S., Michel M., Ding J., Amano T., Forabosco A., et al. (2013). Constitutively active Foxo3 in oocytes preserves ovarian reserve in mice. Nat. Commun. 4, 1843. doi:10.1038/ncomms2861
Poyet M., Le Roux V., Gibert P., Meirland A., Prévost G., Eslin P., et al. (2015). The wide potential trophic niche of the asiatic fruit fly Drosophila suzukii: the key of its invasion success in temperate Europe? PLoS One 10, e0142785. doi:10.1371/journal.pone.0142785
Puig O., Marr M. T., Ruhf M. L., Tjian R. (2003). Control of cell number by Drosophila FOXO: downstream and feedback regulation of the insulin receptor pathway. Genes Dev. 17, 2006–2020. doi:10.1101/gad.1098703
Puig O., Mattila J. (2011). Understanding forkhead box class O function: lessons from Drosophila melanogaster. Antioxid. Redox Sign 14, 635–647. doi:10.1089/ars.2010.3407
Ratnaparkhi A., Sudhakaran J. (2022). Neural pathways in nutrient sensing and insulin signaling. Front. Physiol. 13, 1002183. doi:10.3389/fphys.2022.1002183
Riedel C. G., Dowen R. H., Lourenco G. F., Kirienko N. V., Heimbucher T., West J. A., et al. (2013). DAF-16 employs the chromatin remodeller SWI/SNF to promote stress resistance and longevity. Nat. Cell Biol. 15, 491–501. doi:10.1038/ncb2720
Royet J., Dziarski R. (2007). Peptidoglycan recognition proteins: pleiotropic sensors and effectors of antimicrobial defences. Nat. Rev. Microbiol. 5, 264–277. doi:10.1038/nrmicro1620
Shaw B., Brain P., Wijnen H., Fountain M. T. (2018). Reducing Drosophila suzukii emergence through inter-species competition. Pest Manag. Sci. 74, 1466–1471. doi:10.1002/ps.4836
Shaw B., Fountain M., Wijnen H. (2019). Control of daily locomotor activity patterns in Drosophila suzukii by the circadian clock, light, temperature and social interactions. J. Biol. Rhythms 34, 463–481. doi:10.1177/0748730419869085
Sheng Z., Xu J., Bai H., Zhu F., Palli S. R. (2011). Juvenile hormone regulates vitellogenin gene expression through insulin-like peptide signaling pathway in the red flour beetle, Tribolium castaneum. J. Biol. Chem. 286, 41924–41936. doi:10.1074/jbc.M111.269845
Shi C., Booth L. N., Murphy C. T. (2019). Insulin-like peptides and the mTOR-TFEB pathway protect Caenorhabditis elegans hermaphrodites from mating-induced death. eLife 8, e46413. doi:10.7554/eLife.46413
Siegrist S. E., Haque N. S., Chen C. H., Hay B. A., Hariharan I. K. (2010). Inactivation of both Foxo and reaper promotes long-term adult neurogenesis in Drosophila. Curr. Biol. 20, 643–648. doi:10.1016/j.cub.2010.01.060
Sim C., Denlinger D. L. (2013). Insulin signaling and the regulation of insect diapause. Front. Physiol. 4, 189. doi:10.3389/fphys.2013.00189
Sim C., Kang D. S., Kim S., Bai X., Denlinger D. L. (2015). Identification of FOXO targets that generate diverse features of the diapause phenotype in the mosquito Culex pipiens. Proc. Natl. Acad. Sci. USA. 112, 3811–3816. doi:10.1073/pnas.1502751112
Spellberg M. J., Marr M. T. (2015). FOXO regulates RNA interference in Drosophila and protects from RNA virus infection. Proc. Natl. Acad. Sci. USA. 112, 14587–14592. doi:10.1073/pnas.1517124112
Su Y., Wang T., Wu N., Li D., Fan X., Xu Z., et al. (2019). Alpha-ketoglutarate extends Drosophila lifespan by inhibiting mTOR and activating AMPK. Aging (Albany NY) 11, 4183–4197. doi:10.18632/aging.102045
Süren-Castillo S., Abrisqueta M., Maestro J. L. (2012). FoxO inhibits juvenile hormone biosynthesis and vitellogenin production in the German cockroach. Insect biochem. Mol. Biol. 42, 491–498. doi:10.1016/j.ibmb.2012.03.006
Tatar M., Kopelman A., Epstein D., Tu M. P., Yin C. M., Garofalo R. S. (2001). A mutant Drosophila insulin receptor homolog that extends life-span and impairs neuroendocrine function. Science 292, 107–110. doi:10.1126/science.1057987
Tatar M., Yin C. (2001). Slow aging during insect reproductive diapause: why butterflies, grasshoppers and flies are like worms. Exp. Gerontol. 36, 723–738. doi:10.1016/s0531-5565(00)00238-2
Tatusov R. L., Galperin M. Y., Natale D. A. (2000). The COG database: a tool for genome scale analysis of protein functions and evolution. Nucleic Acids Res. 28, 33–36. doi:10.1093/nar/28.1.33
Towarnicki S. G., Youngson N. A., Corley S. M., St John J. C., Melvin R. G., Turner N., et al. (2022). Ancestral dietary change alters the development of Drosophila larvae through MAPK signaling. Fly. (Austin) 16, 299–311. doi:10.1080/19336934.2022.2088032
Wang J., Zhang H. N., Wang H. D., Zhao S., Zuo Y. Y., Yang Y. H., et al. (2016). Functional validation of cadherin as a receptor of Bt toxin Cry1Ac in Helicoverpa armigera utilizing the CRISPR/Cas9 system. Insect biochem. Mol. Biol. 76, 11–17. doi:10.1016/j.ibmb.2016.06.008
Xing Y., Su T. T., Ruohola-Baker H. (2015). Tie-mediated signal from apoptotic cells protects stem cells in Drosophila melanogaster. Nat. Commun. 6, 7058. doi:10.1038/ncomms8058
Yamamoto R., Tatar M. (2011). Insulin receptor substrate chico acts with the transcription factor FOXO to extend Drosophila lifespan. Aging Cell 10, 729–732. doi:10.1111/j.1474-9726.2011.00716.x
Yan F., Zhao Q., Li Y., Zheng Z., Kong X., Shu C., et al. (2022). The role of oxidative stress in ovarian aging: a review. J. Ovarian Res. 15, 100. doi:10.1186/s13048-022-01032-x
Yang S. A., Wang W. D., Chen C. T., Tseng C. Y., Chen Y. N., Hsu H. J. (2013). FOXO/Fringe is necessary for maintenance of the germline stem cell niche in response to insulin insufficiency. Dev. Biol. 382, 124–135. doi:10.1016/j.ydbio.2013.07.018
Ye J., Fang L., Zheng H. K., Zhang Y., Chen J., Zhang Z. J., et al. (2006). WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 34, W293–W297. doi:10.1093/nar/gkl031
Yuan Z., Lehtinen M. K., Merlo P., Villén J., Gygi S., Bonni A. (2009). Regulation of neuronal cell death by MST1-FOXO1 signaling. J. Biol. Chem. 284, 11285–11292. doi:10.1074/jbc.M900461200
Zeng B., Huang Y., Xu J., Shiotsuki T., Bai H., Palli S. R., et al. (2017). The FOXO transcription factor controls insect growth and development by regulating juvenile hormone degradation in the silkworm, Bombyx mori. J. Biol. Chem. 292, 11659–11669. doi:10.1074/jbc.M117.777797
Zhai Y. F., Lin Q. C., Zhang J. P., Zhang F., Zheng L., Yu Y. (2016). Adult reproductive diapause in Drosophila suzukii females. J. Pest Sci. 89, 679–688. doi:10.1007/s10340-016-0760-9
Zhang M., Zhang X. (2019). The role of PI3K/AKT/FOXO signaling in psoriasis. Arch. Dermatol. Res. 311, 83–91. doi:10.1007/s00403-018-1879-8
Keywords: D. suzukii, DsFoxO, development, fecundity, transcriptome analysis
Citation: Zhao S, Wang R, Liu Y, Su L, Dai X, Qin D, Chen H, Yin Z, Zheng L and Zhai Y (2023) DsFoxO knockout affects development and fecundity of Drosophila suzukii. Front. Physiol. 14:1290732. doi: 10.3389/fphys.2023.1290732
Received: 08 September 2023; Accepted: 26 October 2023;
Published: 10 November 2023.
Edited by:
Natraj Krishnan, Mississippi State University, United StatesReviewed by:
Kui Kang, Zunyi Normal University College, ChinaJustin Flaven-Pouchon, University of Tübingen, Germany
Copyright © 2023 Zhao, Wang, Liu, Su, Dai, Qin, Chen, Yin, Zheng and Zhai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yifan Zhai, enlpZmFuQHRvbS5jb20=