Phenotypes, antioxidant responses, and gene expression changes accompanying a sugar-only diet in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)
BMC Evolutionary Biology volume 17, Article number: 194 (2017)
Diet composition (yeast:carbohydrate ratio) is an important determinant of growth, development, and reproduction. Recent studies have shown that decreased yeast intake elicits numerous transcriptomic changes and enhances somatic maintenance and lifespan, which in turn reduces reproduction in various insects. However, our understanding of the responses leading to a decrease in yeast ratio to 0% is limited.
In the present study, we investigated the effects of a sugar-only diet (SD) on the gene expression patterns of the oriental fruit fly, Bactrocera dorsalis (Hendel), one of the most economically important pests in the family Tephritidae. RNA sequencing analyses showed that flies reared on an SD induced significant changes in the expression levels of genes associated with specific metabolic as well as cell growth and death pathways. Moreover, the observed upregulated genes in energy production and downregulated genes associated with reproduction suggested that SD affects somatic maintenance and reproduction in B. dorsalis. As expected, we observed that SD altered B. dorsalis phenotypes by significantly increasing stress (starvation and desiccation) resistance, decreasing reproduction, but did not extend lifespan compared to those that received a normal diet (ND) regime. In addition, administration of an SD resulted in a reduction in antioxidant enzyme activities and an increase in MDA concentrations, thereby suggesting that antioxidants cannot keep up with the increase in oxidative damage induced by SD regime.
The application of an SD diet induces changes in phenotypes, antioxidant responses, and gene expressions in B. dorsalis. Previous studies have associated extended lifespan with reduced fecundity. The current study did not observe a prolongation of lifespan in B. dorsalis, which instead incurred oxidative damage. The findings of the present study improve our understanding of the molecular, biochemical, and phenotypic response of B. dorsalis to an SD diet.
Nutrients are critical environmental signals that influence growth and development in insects . Macronutrient balance is an important determinant of fitness, and previous studies have shown that altering the concentrations of yeast and sugar has a profound impact on the lifespan of Drosophila melanogaster, thereby suggesting that dietary protein:carbohydrate (P:C) balance is the main driver of lifespan and ageing processes [2,3,4]. In addition, yeast is a complex mixture of nutrients, and lifespan differences between diets with small amounts of yeast or no yeast are explained by deficiencies in both protein and micronutrients in Bactrocera tryoni . In many tephritid flies, a diet combining sugar with yeast hydrolysate, which is a rich source of carbohydrates and proteins, has been shown to enhance reproductive development, sexual performance, fecundity, and lifespan in comparison with a diet of sugar alone . Moderate restriction of protein in the adult diet extends the lifespan of various species ranging from yeast to primates, and the explanations for lifespan extension by dietary protein restriction assumed an adaptive allocation of limited resources between reproductive effort and somatic repair . A higher ratio of protein to carbohydrate in adult diets has been associated with shorter lifespan in insects, thereby suggesting that protein overconsumption is toxic . Dietary proteins are essential for the growth and development of organisms, and a previous study has shown that D. melanogaster fed on a sugar-only diet (SD) can lead to the inhibition of cell growth and cell cycle progression .
Life-history traits are often negatively associated with each other and thus play an important role in insect development . Research investigations have described a trade-off between fecundity and longevity and qualitative and quantitative changes in the diets of multicellular organisms, including Anastrepha ludens, Trichopria drosophilae, and D. melanogaster [11,12,13]. Different diet sources vary in nutritional content, thereby influencing diverse biochemical pathways. Studies have been conducted to underlie the genetic and molecular mechanisms involved in different dietary status, which will assist in bettering understand the relationship among phenotype, gene expression, and the environment [14, 15].
Reactive oxygen species (ROS) are by-products of aerobic metabolism that occurs in the mitochondria that induce oxidative damage to cell, which in turn causes aging and eventually death . Therefore, the activities of different antioxidant enzymes may have a significant effect on the lifespan of organisms. For example, the overexpression of superoxide dismutase (SOD), a major cytosolic enzyme responsible for scavenging highly toxic superoxide radicals in certain tissues, increases the lifespan of Drosophila . On the contrary, knocking out the SOD gene dramatically shortens lifespan of D. melanogaster to a few days by inducing higher levels of oxidative damage . Protein restriction has been demonstrated to reduce the amount of oxidative damage of DNA and other major cellular components. For example, the activities of SOD and catalase (CAT) in liver homogenates of 12- and 24-month-old mice increase with protein restriction . The mechanisms underlying lifespan extension by protein restriction may be the upregulation of SOD activity, which is associated with the insulin signaling pathway . Insulin belongs to a superfamily of peptide hormones, which include relaxins, insulin-like peptides (ILPs) and insulin-like growth factors (IGFs). The insulin signaling pathway plays a central role in nutrition regulated growth . In response to dietary proteins, ILPs function in an endocrine manner and trigger growth by binding to insulin receptors (InRs) and subsequently activating a conserved PI3 kinase (PI3K) and Akt kinase signaling pathways in all tissues of D. melanogaster .
The oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), is one of the most important pests of fruits and vegetables in South East Asia and the Pacific region . It is highly polyphagous, infesting a wide variety of fruit crops such as citrus, mandarin, peach, and mango , and induces significant economic losses through direct fruit damage, fruit drop, and export limitations associated with quarantine restrictions. Moreover, diet is one of the most important environmental variables that affects its growth, reproduction, and distribution; therefore, the ability to adapt to different dietary conditions is critical for their survival . Several studies have investigated the phenotypic response of flies to different nutritional diets; however, the effects of SD on the phenotype and gene expression in B. dorsalis remain unclear. In the present study, gene expression, ovary development, lifespan, somatic maintenance (starvation and desiccation resistance), antioxidant enzyme activities (SOD, glutathione S-transferase (GST), peroxidase (POD), and CAT), the total antioxidant capacity (T-AOC), and malondialdehyde (MDA) content were examined in B. dorsalis after exposure to an SD or a normal diet (ND) regimen. The aim of the present study was to facilitate a better understanding the phenotypes, antioxidant responses, and gene expression patterns in response to the SD regimen in B. dorsalis.
B. dorsalis was collected and reared as previously described . Briefly, the flies were reared in plastic cages at 27 ± 1 °C, 70 ± 5% relative humidity, and a photoperiod of 14:10 (L:D) h.
Variable amounts of yeast extract (Oxoid Ltd., Basingstoke, Hampshire, England) per 100 g of adult food (yeast:sucrose) were used in this study, as well as a sugar-only diet (yeast:sucrose = 0:100, 0% yeast), and normal diet (yeast:sucrose = 25:75, 25% yeast). Yeast mainly comprised 62.5% protein, as well as water soluble B-complex vitamins, and sodium chloride, and had a pH of 7.0 ± 0.2 (0.5% solution), and temperature of 25 °C. The B. dorsalis were given ad libitum access to all diets, assuming that sucrose in each of the diets is not limiting, and that any response to SD is primarily mediated through variations in yeast content.
Insect samples and RNA extraction
B. dorsalis were reared at a density of 15 males and 15 females per jar and randomly assigned to either the SD or ND diet. Thirty-day-old female flies from each diet group were collected. Three replicates were performed for each diet. RNA was extracted using the RNeasy plus Micro Kit (Qiagen GmbH, Hilden, Germany), following the manufacturer’s instructions. RNA was quantified by measuring the absorbance at a wavelength of 260 nm using a NanoVue UV-Vis spectrophotometer (GE Healthcare Bio-Science, Uppsala, Sweden). The purity of all RNA samples was assessed at an absorbance ratio of OD260/280 and OD260/230, and the integrity of RNA was confirmed by 1% agarose gel electrophoresis.
cDNA library construction and sequencing
Briefly, oligo (dT) magnetic beads were used to select mRNAs with a poly(A) tail, and the target RNA was obtained after purification. First-strand cDNA was synthesized using random hexamer-primers from purified poly(A) mRNA. Second-strand cDNA was synthesized using buffer, dNTPs, RNaseH, and DNA polymerase I. Short fragments were purified using a QiaQuick PCR extraction kit. These fragments were washed with ethidium bromide buffer for end reparation poly(A) addition and then ligated to sequencing adapters. Suitable fragments, as judged by agarose gel electrophoresis, were selected as templates for PCR amplification. The resulting cDNA library was sequenced on a next-generation sequencing platforms (complete genomics) using paired-end technology in a single run.
Analysis of RNA-seq data
Raw reads produced in this study were cleaned by discarding reads with adapters and reads in which unknown bases comprised more than 10%. Low quality reads (a percentage of low quality bases over 50% in a read) were also excluded from further analysis. HISAT was used to map clean reads to the B. dorsalis genome reference and Bowtie2  to the gene reference using the default parameters, respectively. Gene expression levels in terms of transcripts were quantified by RSEM (RNA-Seq by Expectation Maximization) and FPKM (Fragment Per Kilobase of exon model per million mapped reads) method . The FPKM between the biological replications was analyzed by Pearson correlation. The Pearson coefficient of gene expression in different replications was more than 0.85, indicating consistency between the replicates. When the value of either sample FPKM was zero, 0.01 was used to instead of 0 to calculate the fold change. According to the correlation results, the NOISeq method was selected to analyze differential expression between the SD and ND diets . We used the absolute value of log2 ratio ≥ 1 and probability ≥0.8 as the threshold to judge significant differences in gene expression .
Gene Ontology (GO) functional analysis provides GO functional classification annotation of DEGs as well as GO functional enrichment analysis. The annotation terms formed the GO ontology and were obtained from Blast2GO . This method first maps all DEGs to GO terms in the database (www.geneontology.org), calculating gene numbers for every term, then uses a hypergeometric test to find significantly enriched GO terms in the input list of DEGs based on ‘GO::TermFinder’ (http://www.yeastgenome.org/help/analyze/go-term-finder). The p-value for the hypothesis test was calculated using the following formula:
where N is the number of all genes with a GO annotation; n is the number of DEGs in N; M is the number of all genes that were annotated to certain GO terms; m is the number of DEGs in M. The calculated p-value was subjected to Bonferroni correction, taking the corrected p-value <0.05 as a threshold. GO terms fulfilling this condition were defined as significantly enriched GO terms in DEGs.
Pathway analysis helps in further elucidating the distinct functions of gene regulation and identify the significant pathways of the DEGs according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, which is a major public database on pathway enrichment analysis . The formula for the p-value is similar to that of the GO analysis. Here, N is the number of all genes with KEGG annotation, n is the number of DEGs in N, M is the number of all genes annotated to specific pathways, and m is the number of DEGs in M.
Effect of SD diet on ovary development
For the SD and ND diets, five jars (100 mm high × 70 mm diameter) were each set up for 10 females and 10 males, making a total of 10 jars (n = 50 females per diet). Food and water were both supplied from two small vessels and placed in each jar. Food was replaced every 7 days, and water was replaced daily.
To maintain a 1:1 sex ratio in each jar, dead males were replaced with age-matched males until the female adult died. For ovary development, 30-day-old female flies from each jar were dissected, and the maximum diameter of each entire ovary was measured. The images were captured with a Leica M205A stereomicroscope (Leica Microsystems, Wetzlar, Germany).
Starvation and desiccation treatment
For the SD and ND diets, five jars (100 mm high × 70 mm diameter) were each set up for 10 females and 10 males, making a total of 10 jars (n = 50 females per diet). Diet and water were both supplied from two small vessels and placed in each jar. To maintain a 1:1 sex ratio in each jar, dead males were replaced with age-matched males. When a female died, a male was also removed. For starvation diet, the female flies that emerged in each jar at day were transferred into a new jar containing a small vessel plugged with cotton to prevent desiccation. After transferring the flies, the number of dead flies was recorded each day until all were dead. For desiccation diet, the female flies that emerged in each jar at day 30 were transferred into a new jar containing a disc of dry filter paper. Desiccation jars were observed every 12 h after transferring the flies until all were dead.
Effect of SD diet on fly lifespan
For the SD and ND diets, we prepared 10 jars (100 mm high × 70 mm diameter) of 5 females and 5 males, making a total of 20 jars (n = 50 females per diet). Flies were maintained at a 1:1 sex ratio. When a female died, an opposite-sex individual was removed, and when a male died, it was replaced with an age-matched individual of the same sex. The flies were supplied with food and water ad libitum using the food and water-vessel system as previously described.
Protein extraction was conducted using a commercially available assay kit (Nanjing Jiancheng Bioengineering Institute, Jiangsu, China). The treated females were first homogenized in 0.9% saline at a ratio of 1:9 (Wflies:Vsaline). The crude homogenates were centrifuged at 4 °C and 10,000 g for 10 min. The supernatants were collected and stored on ice until assayed. The protein concentration was determined using the Bradford (1976)  method, with bovine serum albumin as standard.
Enzyme activity, T-AOC, and lipid peroxidation (LPO) assay
Assay kits for SOD (Item number: A001–3), POD (Item number: A084–1), CAT (Item number: A007–2), T-AOC (Item number: A015–1), and MDA (Item number: A003–1) were used in this study (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Absorbances were read in a microplate spectrophotometer (XMark™, Bio-Rad, Hercules, CA, USA). Specifically, SOD activity was determined at a wavelength of 550 nm using xanthine and xanthine oxidase systems. One unit of SOD activity was defined as the amount of enzyme required to cause 50% inhibition of xanthine oxidase in 1 mL of enzyme extraction of 1 mg protein. POD activity was determined at a wavelength of 420 nm by catalyzing the oxidation in the presence of H2O2 as a substrate. One unit of POD activity was defined as the amount that catalyzes 1 μg of substrate per min per mg protein. CAT activity was calculated by measuring the decrease in absorbance at a wavelength of 240 nm due to H2O2 decomposition. One unit of CAT activity was defined as the amount that decomposes 1 μmol of H2O2 per sec per mg protein. GST activity was determined with 1-chloro-2,4-dinitrobenzene (CDNB, Shanghai Chemicals, Shanghai, China) and reduced GSH (Sigma-Aldrich, St. Louis, MO, USA) as substrates according to Habig et al. (1974) , with slight modifications. Briefly, a 96-well microplate with 100 μL CDNB (1% ethanol (V/V) included) and 100 μL GSH in 0.9% saline in each well was incubated for 20 min at 37 °C. Subsequently, 100 μL of each enzyme solution was added to individual wells to a final concentration of 0.2 mM CDNB and 2.0 mM GSH. Absorbance was measured continuously at a wavelength of 340 nm and 37 °C for 5 min. Absorbance changes per min were converted into nmol CDNB conjugated/min/mg protein using the extinction coefficient of the resulting 2,4-dinitrophenyl-glutathione: ε340 nm = 9.6 mM−1 cm−1 .
T-AOC was measured based on the generation of the Fe2+-o-phenanthroline complex, as the overall reducing agents in the sample supernatant that reduced Fe3+ to Fe2+ in the sample supernatant. The absorbance of the Fe2+-o-phenanthroline complex was measured at a wavelength of 520 nm. One unit of T-AOC was defined as the amount necessary to increase absorbance by 0.01 per min per mg protein. LPO levels were indirectly determined by using MDA levels. MDA production was evaluated by reacting with thiobarbituric acid to yield a red species with a maximum absorbance at 532 nm. The MDA concentration was expressed as nmol of MDA produced per mg protein.
Data were analyzed using SPSS 16.0 software (SPSS Inc., Chicago, IL, USA, 2008) and are presented as the mean ± SE. Significant differences between treatments were analyzed by using an independent Student’s t-test (for comparison of two means) (* P < 0.05, ** P < 0.01, *** P < 0.001).
Evaluation of RNA-seq data
In the present study, samples from B. dorsalis were sequenced using the RNA-seq technology. More than 99.98% of the clean data was obtained from each sample for pair-end sequencing, and the length of sequence reads was 50 bp (Additional file 1: Table S1). To evaluate the quality of the sequencing data, several aspects of each B. dorsalis sample were subjected to strict quality control and the average genome mapping ratio was 83.22% (Additional file 2: Table S2). The FPKM density had a similar pattern in each sample of B. dorsalis, indicating that the transcriptome analysis of each diet was highly reproducible (Additional file 3: Figure S1).
Analysis of gene expression profiling
The gene expression profiles of the flies fed on an SD or ND were analyzed at day 30 after emergence of B. dorsalis. A total of 2600 (1656 downregulated, 944 upregulated) DEGs were generated between the SD and ND groups, and the magnitude of these changes were defined as abs (log2(Y/X) ≥ 1 and divergence probability ≥0.8) (Additional file 4: Figure S2 and Table S3).
For GO analysis, DEGs were divided into three ontologies: biological process, cellular component, and molecular function, including 51 annotations (Additional file 5: Figure S3A). The majority of the enriched and upregulated genes were assigned into the following GO categories (Additional file 6: Table S4): ‘generation of precursor metabolites and energy’ (10.7%, P = 1.15e-23), ‘energy derivation by oxidation of organic compounds’ (9.8%, P = 5.54e-23), ‘ATP metabolic process’ (2.9%, P = 9.23e-08), ‘glycosyl compound biosynthetic process’ (2.9%, P = 2.68e-07), ‘somatic muscle development’ (1.6%, P = 0.0088), and ‘aerobic respiration’ (2.5%, P = 0.035). Most of the downregulated genes were included in the ‘chromosome organization’ (10.1%, P = 2.71e-20), ‘cell cycle’ (17%, P = 6.48e-20), ‘organelle organization’ (18.4%, P = 8.01e-07), ‘negative regulation of metabolic process’ (6.3%, P = 9.10e-07), ‘reproduction’ (16%, P = 4.99e-06), and ‘female gamete generation’ (9.1%, P = 6.93e-05) categories (Additional file 7: Table S5).
For KEGG analysis, the DEGs were enriched into 32 KEGG pathways, including the categories “cellular process”, “environmental information processing”, “genetic information processing”, “metabolism”, and “organismal system” (Additional file 5: Figure S3B). The upregulated genes were assigned to the following KEGG pathways: namely, ‘oxidative phosphorylation’ (9.9%, P = 3.81e-41) (Fig. 1), ‘metabolic pathways’ (23.8%, P = 2.41e-10), ‘calcium signaling pathway’ (3.7%, P = 2.86e-06), ‘MAPK signaling pathway’ (4.3%, P = 2.82e-05), ‘folate biosynthesis’ (1.35%, P = 0.000292), ‘biosynthesis of amino acids’ (2.1%, P = 0.00158), ‘alanine, aspartate and glutamate metabolism’ (1.2%, P = 0.00281), ‘drug metabolism cytochrome P450’ (2.0%, P = 0.0056), ‘cAMP signaling pathway’ (3.1%, P = 0.0068), and ‘cysteine and methionine metabolism’ (1.2%, P = 0.0327) (Table 1). The downregulated genes were enriched in the ‘cell cycle’ (4.0%, P = 9.25e-19), ‘insulin signaling pathway’ (1.4%, P = 0.015) (Fig. 2a), ‘the target of rapamycin (TOR) signaling pathway’ (0.41%, P = 0.0042) (Additional file 8: Figure S4), ‘insect hormone biosynthesis (including juvenile hormone (JH) biosynthesis)’ (0.90%, P = 0.036) (Fig. 2b), and ‘p53 signaling pathway’ (0.98%, P-values = 0.012) categories (Table 2).
The ovary size of the female adults reared with the ND was significantly greater than those of females reared with SD (P < 0.001) (Fig. 3a). The average ovary size of the ND group was 2.04 ± 0.01 mm, whereas the average ovary size for the SD group was 0.61 ± 0.003 mm (Fig. 3b).
Starvation and desiccation treatments
Starvation resistance was significantly affected by the quality of the diets (Fig. 4a-b), and B. dorsalis fed on an SD had longer survival time than those on the ND, and the mean survival time was 4.82 ± 0.49 and 2.21 ± 0.06 days, respectively (P = 0.006). In terms of the desiccation treatment, the mean survival time was 55.60 ± 2.53 and 43.57 ± 2.53 h under the SD and ND diets, respectively, which suggested that SD significantly increased the desiccation resistance (P = 0.01) (Fig. 4c, d).
No significant effect on lifespan was observed using different diets (P = 0.27) (Fig. 5a). The mean lifespan was 37.05 ± 2.25 days and 40.64 ± 2.24 days under the SD and ND diets, respectively (Fig. 5b). These results suggested that SD repressed ovary development without extending the lifespan of B. dorsalis.
The antioxidant enzyme activities of B. dorsalis fed on the SD and ND for 30 days significantly differed (Fig. 6a-d). Flies fed on the SD showed lower SOD (P = 0.007), GST (P = 0.002), POD (P = 0.003), and CAT (P = 0.02) activities than those on ND. Compared to the ND diet, B. dorsalis exhibited lower T-AOC when fed on the SD diet (P < 0.001) (Fig. 6e). Flies fed on the SD had higher MDA concentrations (P < 0.001), suggesting that they had higher lipid peroxidation levels (Fig. 6f).
In the present study, we used RNA-seq to identify genes that were regulated by nutrient signals in B. dorsalis. Similar to Drosophila , we observed that B. dorsalis elicited a strong response to SD as revealed by changes in gene expression levels (Additional file 4: Figure S2), and a large number of DEGs were significantly enriched in metabolic pathways (Additional file 5: Figure S3). These changes may have led to the alterations in the storage and metabolism of fats and carbohydrates [21, 35], as well as induction of autophagy , which in turn resulted in the growth and developmental arrest of B. dorsalis, as well as maintenance of homeostasis under poor nutrient conditions . We also observed an inhibition of ‘cell growth and death’ pathways (Table 2), suggesting that an SD prevent cell growth and cell cycle progression in B. dorsalis . Three pathways that are associated with transcriptional activities were downregulated under an SD regimen, including: ‘RNA degradation’, ‘spliceosome’ and ‘mRNA surveillance’ (Table 2), which could affect the fidelity and quality of mRNA molecules . Calcium signaling pathway was activated in response to the SD (Table 1), and a previous study has reported that calcium and calmodulin are key components of signal transduction pathways that may be involved in stress responses to various environmental stressors . Similar to a previous report that showed that the nutritional environment of insects influences the immune status of D. melanogaster, the ‘immune system’ category was also significantly enriched with 28 upregulated genes under the SD regime (Table 1), which suggests that the chronic nutritional stress affects adult immunity by influencing plastic allocation of resources . Another interesting finding is that genes involved in the folate biosynthesis pathway were upregulated (Table 1); these have been implicated in the aging process in various organisms [41, 42]. In particular, folate is a precursor to methionine, which is then converted into s-adenosyl methionine (SAM), a cell-wide donor of methyl groups. It is through this pathway that folate levels have been associated with DNA methylation levels and aging . Protein starvation disrupts amino acid and fatty acid metabolism in Drosophila . Similarly, the ‘biosynthesis of the amino acids alanine’, ‘aspartate and glutamate metabolism’, and ‘Cysteine and methionine metabolism’ pathways were also significantly enriched in our study. The complete lack of yeast in the diet would certainly disrupt amino acid and protein biosynthesis and metabolism, and the flies responded to this stress by increasing their rate for amino acids metabolism.
Additionally, the current study found that 66 upregulated genes were significantly enriched in the ‘oxidative phosphorylation’ pathway, which have been reported to play a central role in eukaryotic metabolism to provide the energy (ATP) that is required for insect survival . The overexpression of genes that are involved in energy production in B. dorsalis might reflect the costs of survival in the chronic SD conditions. Interestingly, as one of the major pathways associated with lifespan, the oxidative phosphorylation pathway is also activated in D. melanogaster under protein restricted diet . Furthermore, 112 (16.0%) downregulated genes were enriched in the ‘reproductive process’ GO terms (Additional file 7: Table S5). Here, the upregulation of pathways involved in energy production and the downregulation of pathways involved in reproduction suggest a trade-off between reproduction versus somatic maintenance (stress resistance) and lifespan . Similarly, in most insects reproduction trades off with somatic maintenance and lifespan exists, and it has been documented at the phenotypic, physiological, or quantitative genetic level in numerous species and represents a major constraint upon physiology and the evolution of life histories [47, 48]. At the physiological level, the possible explanation of this trade-off is that the energetically costly process of reproduction ‘competes’ with the energy demands of somatic maintenance and survival . At the same time, Drosophila and other solitary insects have showed that this trade-off is achieved by an endocrine network that integrates insulin signaling pathway, JH, and the yolk precursor vitellogenin (Vg) . Although, the extended lifespan is often associated with reduced fecundity, recent studies have also shown that fecundity and lifespan have different nutritional optima, and that there is no concomitant trade-off, indicating a controversial area .
As reported, the negative effect of InR and insulin receptor substrate (IRS) RNAi-mediated silencing on reproductive events has been observed in several insect species, thereby suggesting the involvement of the insulin pathway in the control of reproductive system development [52, 53]. Two recent studies have also reported that the nutritional environment affects the hormone production in insects, and knocking down the expression of genes involved in JH biosynthesis (JH acid methyltransferase) blocks Vg synthesis and disrupts ovarian development [54, 55]. A similar negative effect on JH biosynthesis and adult Vg has also been obtained by means of InR RNAi in the penultimate and last instar nymphs of the Blattella germanica . Insulin signaling pathway affects JH production, which is the central regulator of fecundity, thereby regulating the synthesis of Vg . As precursor proteins of vitellin, the major yolk egg protein, Vg is directly involved in reproduction. Downregulation of this network, for example in response to low nutrient availability or other environmental changes (temperature or photoperiod), promotes maintenance and survival at the expense of reproduction .
In the present study, we found that SD downregulates the nutritional sensor of the TOR (Additional file 8: Figure S4), insulin (Fig. 2a), and JH biosynthesis (Fig. 2b) signaling pathways in B. dorsalis. The TOR, which is a serine/threonine kinase, is highly conserved in most eukaryotes and plays a key role in the transduction of nutritional signals. The direct inhibition of TOR signaling pathway has been shown to extend the lifespan of yeast , flies , and mammals . Here, among the down-regulated genes, Vg1 (gene ID: 105,232,170; diverge probability = 0.99) and Vg2 (gene ID: 105,222,970; diverge probability = 0.99) had more than 100-fold lower expression under the SD diet in B. dorsalis (Additional file 9: Table S3). Therefore, these results suggested that insulin signaling pathway asserts its nutritionally linked influence on Vg expression of B. dorsalis by contributing to the control of biosynthesis and secretion of JH under the SD regimen .
As expected, ovary development of B. dorsalis under the SD regimen was completely blocked, and no mature eggs were observed (Fig. 3), which implied that dietary protein plays an important role in the regulation of ovarian development in B. dorsalis . Similar to D. melanogaster, the desiccation resistance has a positive relationship with starvation resistance, and the increasing levels of starvation resistance were also more resistant to desiccation . High sugar diets are often associated with higher triglyceride levels and therefore, have better starvation resistance . In the present study, SD induced an increase in starvation and desiccation resistance in B. dorsalis, thereby suggesting a trade-off between stress resistance and reproduction under the SD diet. According to the several aging theories, the extended lifespan is often correlated with increased resistance against various stressors, and several studies have shown genetic correlations among resistance to several stresses and longevity, thereby strongly suggesting that genetic basis of lifespan and stress resistance overlap [65,66,67].
However, our study found that the downregulated Insulin-JH-Vg signaling network, decreased reproduction and increased stress resistance were not associated with the lifespan extension under the SD regimen in B. dorsalis. Antioxidant activity and the accumulation of oxidative damage caused by ROS are closely related to lifespan in various organisms . ROS are formed by the incomplete reduction of oxygen, primarily in the mitochondria during aerobic metabolism, and antioxidant systems were usually able to keep this oxidative stress under control . When ROS production exceeds antioxidant defenses, cells enter a pro-oxidant state called oxidative stress and damage occurs. When this damage is not repaired, damages accumulate in cells, which in turn results in further oxidative damage, thereby causing aging and ultimately, death [69, 70]. Moreover, a study involving Teleogryllus commodus has shown that dietary manipulation have significant effects on oxidative damage and antioxidant protection . Intriguingly, our RNA-seq data showed that the antioxidant gene expressions of SOD (gene ID: 105,232,170; diverge probability = 0.90), GST (gene ID: 105,229,686; diverge probability = 0.96), and CAT (gene ID: 105,222,970; diverge probability = 0.87) were downregulated under the SD diet (Additional file 9: Table S3). Meanwhile, in agreement with our RNA-seq data, the present study showed that the SD decreased antioxidant enzyme activities (of SOD, GST, POD, and CAT) and antioxidant capacity (T-AOC) in B. dorsalis. As a major oxidation product of peroxidized polyunsaturated fatty acids, MDA can be used as a biological marker for oxidative damage. Here, we also found that SD increased the MDA concentrations of B. dorsalis compared to ND regime (Fig. 6). These results indicated that antioxidants of B. dorsalis cannot keep up with the increase in oxidative damage that was induced by SD diets. Hence, we believe that the flies did not extend the lifespan amidst oxidative damage accumulation and under the unnatural and extreme diets (SD). Similarly, the relationship between sugar-only diet and lifespan have also been studied in Anastrepha ludens and Drosophila, which showed less fecundity without lifespan extension under SD regimen [11, 72].
Our data show that SD leads to widespread changes in the expression levels of genes involved in energy production, reproduction, as well as the nutritional sensor of the TOR, insulin, and JH biosynthesis signaling pathways in B. dorsalis. The SD also significantly affected the phenotypes and antioxidant responses, including an increase in stress resistance, disruption of ovary development, decrease in antioxidant enzyme activities, and accumulation of oxidative damage. Furthermore, yeast provides protein, micronutrients, and salt, although which component or combination of components plays the major role remains unclear. We anticipate conducting another study that will focus on elucidating this mechanism.
Cyclic adenosine 3′,5′-monophosphate
Different expression genes
Fragment per kilobase of exon model per million mapped reads
Hierarchical indexing for spliced alignment of transcripts
Insulin-like growth factors
Insulin receptor substrate
Kyoto encyclopedia of genes and genomes
Mitogen-activated protein kinase
Reactive oxygen species
RNA-Seq by expectation maximization
Total antioxidant capacity
Target of rapamycin
Okamoto N, Yamanaka N. Nutrition-dependent control of insect development by insulin-like peptides. Curr Opin Insect Sci. 2015;11:21–30.
Bruce KD, Hoxha S, Carvalho GB, Yamada R, Wang HD, Karayan P, et al. High carbohydrate-low protein consumption maximizes Drosophila lifespan. Exp Gerontol. 2013;48:1129–35.
Lee KP, Simpson SJ, Clissold FJ, Brooks R, Ballard JWO, Taylor PW, et al. Lifespan and reproduction in Drosophila: new insights from nutritional geometry. P Natl Acad of Sci USA. 2008;105:2498–503.
Lee KP. Dietary protein:carbohydrate balance is a critical modulator of lifespan and reproduction in Drosophila melanogaster: a test using a chemically defined diet. J Insect Physiol. 2015;75:12–9.
Fanson BG, Taylor PW. Additive and interactive effects of nutrient classes on longevity, reproduction, and diet consumption in the Queensland fruit fly (Bactrocera tryoni). J Insect Physiol. 2012;58:327–34.
Liedo A, Orozco D, Cruz-López L, Quintero JL, Becerra-Pérez C, del Refugio HM, et al. Effect of post-teneral diets on the performance of sterile Anastrepha ludens and Anastrepha obliqua fruit flies. J Appl Entomol. 2013;137:49–60.
Holliday R. Food, reproduction and longevity: is the extended lifespan of calorie-restricted animals an evolutionary adaptation? BioEssays. 1989;10:125–7.
Simpson SJ, Raubenheimer D. Macronutrient balance and lifespan. Aging. 2009;10:875–80.
Britton JS, Edgar BA. Environmental control of the cell cycle in Drosophila: nutrition activates mitotic and endoreplicative cells by distinct mechanisms. Development. 1998;125:2149–58.
Zera AJ, Harshman LG. The physiology of life history trade-offs in animals. Annu Rev Ecol S. 2003;32:95–126.
Carey JR, Harshman LG, Liedo P, Müller HG, Wang JL, Zhang Z. Longevity fertility trade-offs in the tephritid fruit fly, Anastrepha ludens, across dietary restriction gradients. Aging Cell. 2008;7:470–7.
Ellers J, Ruhe B, Visser B. Discriminating between energetic content and dietary composition as an explanation for dietary restriction effects. J Insect Physiol. 2011;57:1670–6.
Flatt T. Survival costs of reproduction in Drosophila. Exp Gerontol. 2011;46:369–75.
Zinke I, Schütz CS, Katzenberger JD, Bauer M, Pankratz MJ. Nutrient control of gene expression in Drosophila: microarray analysis of starvation and sugar-dependent response. EMBO J. 2002;21:6162–73.
Magnanou E, Malenke JR, Dearing MD. Expression of biotransformation genes in woodrat (Neotoma) herbivores on novel and ancestral diets: identification of candidate genes responsible for dietary shifts. Mol Ecol. 2009;18:2401–14.
Salmon AB, Richardson A, Perez VI. Update on the oxidative stress theory of aging: does oxidative stress play a role in aging or healthy aging? Free Radic Biol Med. 2010;48:642–55.
Sun XP, Komatsu T, Lim J, Laslo M, Yolitz J, Wang C, et al. Nutrient-dependent requirement for SOD1 in lifespan extension by protein restriction in Drosophila melanogaster. Aging Cell. 2012;11:783–93.
Sanz A. Mitochondrial reactive oxygen species: do they extend or shorten animal lifespan? BBA-Bioenergetics. 2016;1857:1116–26.
Koizumi A, Weindruch R, Walford RL. Influences of dietary restriction and age on liver enzyme activities and lipid peroxidation in mice. J Nutr. 1987;117:361–7.
Kabil H, Partridge L, Harshman LG. Superoxide dismutase activities in long-lived Drosophila melanogaster females: chico 1 genotypes and dietary dilution. Biogerontology. 2007;8:201–8.
Britton JS, Lockwood WK, Li L, Cohen SM, Edgar BA. Drosophila’s insulin/ PI3-kinase pathway coordinates cellular metabolism with nutritional conditions. Dev Cell. 2002;2:239–49.
Oldham S, Hafen E. Insulin/IGF and target of rapamycin signaling: a TOR de force in growth control. Trends Cell Biol. 2003;13:79–85.
Clarke AR, Armstrong KF, Carmichael AE, Milne JR, Raghu S, Roderick GK, et al. Invasive phytophagous pests arising through a recent tropical evolutionary raditation: the Bactrocera dorsalis complex of fruit flies. Annu Rev Entomol. 2005;50:293–319.
Jin T, Zeng L, Lin YY, Lu YY, Liang GW. Insecticide resistance of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), in mainland China. Pest Manag Sci. 2011;67:370–6.
Carsten LD, Watts T, Markow TA. Gene expression patterns accompanying a dietary shift in Drosophila melanogaster. Mol Ecol. 2005;14:3203–8.
Shen GM, Dou W, Huang Y, Jiang XZ, Smagghe G, Wang JJ. In silico cloning and annotation of genes involved in the digestion, detoxification and RNA interference mechanism in the midgut of Bactrocera dorsalis [Hendel (Diptera: Tephritidae)]. Insect Mol Biol. 2013;22:354–65.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:25–34.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12:323.
Tarazona S, Garcíaalcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21:2213–23.
Zhao YQ, Li J, Liu HJ, Xi Y, Xue M, Liu W, et al. Dynamic transcriptome profiles of skeletal muscle tissue across 11 developmental stages for both Tong cheng and Yorkshire pigs. BMC Genomics. 2015;16:377.
Conesa A, Götz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Re. 2007;35:182–5.
Bradford MA. A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein dye binding. Anal Biochem. 1976;72:248–54.
Habig WH, Pabst MJ, Jakoby WB. Glutathione S-transferase. J Biol Chem. 1974;249:7130–9.
Teleman AA, Chen YW, Cohen SM. 4E-BP functions as a metabolic brake used under stress conditions but not during normal growth. Genes Dev. 2005;19:1844–8.
Scott RC, Schuldiner O, Neufeld TP. Role and regulation of starvation induced autophagy in the Drosophila fat body. Dev Cell. 2004;7:167–78.
Scott RC, Juhasz G, Neufeld TP. Direct induction of autophagy by Atg1 inhibits cell growth and induces apoptotic cell death. Curr Biol. 2007;17:1–11.
Li L, Edgar BA, Grewal SS. Nutritional control of gene expression in Drosophila larvae via TOR, Myc and a novel cisregulatory element. BMC Cell Biol. 2010;11:7.
Ding F, Gil MP, Franklin M, Ferreira J, Tatar M, Helfand SL, et al. Transcriptional response to dietary restriction in Drosophila melanogaster. J Insect Physiol. 2014;69:101–6.
Fellous S, Lazzaro BP. Larval food quality affects adult (but not larval) immune gene expression independent of effects on general condition. Mol Ecol. 2010;19:1462–8.
Cabreiro F, Au C, Leung KY, Vergara-Irigaray N, Cochemé HM, Noori T, et al. Metformin retards aging in C. elegans by altering microbial folate and methionine metabolism. Cell. 2013;153:228–39.
Virk B, Correia G, Dixon DP, Feyst I, Jia J, Oberleitner N, et al. Excessive folate synthesis limits lifespan in the C. elegans: E. coli aging model. BMC Biol. 2012;10:67.
Kim KC, Friso S, Choi SW. DNA methylation, an epigenetic mechanism connecting folate to healthy embryonic development and aging. J Nutr Biochem. 2009;20:917–26.
Wilson DF, Vinogradov SA. Mitochondrial cytochrome c oxidase: mechanism of action and role in regulating oxidative phosphorylation. J Appl Physiol. 2014;117:1431–9.
Whitaker R, Gil MP, Ding F, Tatar M, Helfand SL, Neretti N. Dietary switch reveals fast coordinated gene expression changes in Drosophila melanogaster. Aging. 2014;6:355–68.
Parthasarathy R, Palli SR. Molecular analysis of nutritional and hormonal regulation of female reproduction in the red flour beetle, Tribolium castaneum. Insect Biochem Mol Biol. 2011;41:294–305.
Harshman L, Zera A. The cost of reproduction: the devil in the details. Trends Ecol Evol. 2007;22:80–6.
Flatt T, Heyland A. Mechanisms of life history evolution-the genetics and physiology of life history traits and trade-offs. Mechanisms of life history evolution. Oxford University Press; 2011.
Hansen M, Flatt T, Aguilaniu H. Reproduction, fat metabolism, and life span: what is the connection? Cell Metab. 2013;17:10–9.
Rodrigues AM, Flatt T. Endocrine uncoupling of the trade-off between reproduction and somatic maintenance in eusocial insects. Curr Opin Insect Sci. 2016;16:1–8.
Adler MI, Cassidy EJ, Fricke C, Bonduriansky R. The lifespan-reproduction trade-off under dietary restriction is sex-specific and context-dependent. Exp Gerontol. 2013;48:539–48.
Brown MR, Clark KD, Gulia M, Zhao ZW, Garczynski SF, Crim JW, et al. An insulin-like peptide regulates egg maturation and metabolism in the mosquito Aedes aegypti. Proc Natl Acad Sci USA. 2008;105:5716–21.
Abrisqueta M, Süren-Castillo S, Maestro JL. Insulin receptor-mediated nutritional signalling regulates juvenile hormone biosynthesis and vitellogenin production in the German cockroach. Insect Biochem Mol Biol. 2014;49:14–23.
Parthasarathy R, Sun Z, Bai H, Palli SR. Juvenile hormone regulation of vitellogenin synthesis in the red flour beetle, Tribolium castaneum. Insect Biochem Mol Biol. 2010;40:405–14.
Tatar M, Kopelman A, Epstein D, Tu MP, Yi CM, Garofalo RS. A mutant Drosophila insulin receptor homolog that extends life-span and impairs neuroendocrine function. Science. 2001;292:107–10.
Toivonen JM, Partridge L. Endocrine regulation of aging and reproduction in Drosophila. Mol Cell Endocrinolo. 2009;299:39–50.
Flatt T, Amdam GV, Kirkwood TBL, Omholt SW. Life-history evolution and the polyphenic regulation of somatic maintenance and survival. Q Rev Biol. 2013;88:185–218.
Kaeberlein M, Powers RW, Steffen KK. Regulation of yeast replicative life span by TOR and Sch9 in response to nutrients. Science. 2005;310:1193–6.
Kapahi P, Zid BM, Harper T. Regulation of lifespan in Drosophila by modulation of genes in the TOR signaling pathway. Curr Biol. 2004;14:885–90.
Harrison DE, Strong R, Sharp ZD. Rapamycin fed late in life extends lifespan in genetically heterogeneous mice. Nature. 2009;460:392–5.
Smykal V, Raikhel AS. Nutritional control of insect reproduction. Curr Opin Insect Sci. 2015;11:31–8.
Kai L, Xia C, Liu WT, Zhang XY, Chen MX, Qiang Z. Nutritional signaling regulates vitellogenin synthesis and egg development through juvenile hormone in Nilaparvata lugens. Int J Mol Sci. 2016;17:269.
Hoffmann AA, Hallas R, Sinclair C, Mitrovski P. Levels of variation in stress resistance in Drosophila among strains, local populations, and geographic regions: patterns for desiccation, starvation, cold resistance, and associated traits. Evolution. 2001;55:1621–30.
Burger JM, Hwangbo DS, Corby-Harris V, Promislow DE. The functional costs and benefits of dietary restriction in Drosophila. Aging Cell. 2006;6:63–71.
Vermeulen CJ, Loeschcke V. Longevity and the stress response in Drosophila. Exp Gerontol. 2000;42:153–9.
Bubliy OA, Loeschcke V. Correlated responses to selection for stress resistance and longevity in a laboratory population of Drosophila melanogaster. J Evolution Biol. 2005;18:789–803.
Mockett RJ, Orr WC, Rahmandar JJ, Sohal BH, Sohal RS. Antioxidant status and stress resistance in long- and short-lived lines of Drosophila melanogaster. Exp Gerontol. 2001;36:441–63.
Barja G. Mitochondrial oxygen consumption and reactive oxygen species production are independently modulated: Implications for aging studies. Rejuv Res. 2007;10:215–24.
Finkel T, Holbrook NJ. Oxidants, oxidative stress and the biology of ageing. Nature. 2000;408:239–47.
Harman D. A theory based on free radical and radiation chemistry. J Gerontol. 1956;11:298–300.
Archer CR, Hempenstall S, Royle NJ, Selman C, Willis S, Rapkin J, et al. Testing the effects of dl-Alpha-Tocopherol supplementation on oxidative damage, total antioxidant protection and the sex-specific responses of reproductive effort and lifespan to dietary manipulation in Australian field Crickets (Teleogryllus commodus). Antioxidants. 2015;4:768–92.
Deshpande SA, Yamada R, Mak CM, Hunter B, Soto OA, Hoxha S, et al. Acidic food pH increases palatability and consumption and extends Drosophila lifespan. J Nutr. 2015;145:2789–96.
This study was supported in part by the National Key Research and Development Program (2016YFC1200600), the earmarked fund for Modern Agro-industry (Citrus) Technology Research System (CARS-27), and the Fundamental Research Funds for the Central Universities, China (2362015xk04).
Availability of data and materials
RNA-Seq raw sequences generated using complete genomics, and have been submitted to SRA at NCBI under the Accession No. SUB2028322 and BioSample accession: SAMN05929037. Bactrocera dorsalis genome sequence data: NCBI Assembly: ASM78921v2. The phenotypic data are accessible at the Dryad repository: doi:10.5061/dryad.v73p.
Ethics approval and consent to participate
No special permits were required in the collections of oriental fruit flies.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Data output quality and mapping rates for the examined samples of Bactrocera dorsalis. (DOCX 17 kb)
The items of quality control for each sample of Bactrocera dorsalis. (DOCX 17 kb)
Histogram distribution of gene expression levels of each sample. X-axis is FPKM value (the coordinate has been changed by logarithm for better view). Y-axis is gene number of corresponding FPKM. A: ND-1, ND-2 and ND-3; B: SD-1, SD-2 and SD-3. (JPEG 565 kb)
Analysis of DEGs between the two diets. The DEGs were defined as abs (log2(Y/X) ≥ 1 and diverge probability ≥0.8). (JPEG 621 kb)
GO and KEGG classification on DEGs for sugar-only diet and normal diet. (A) GO terms for DEGs. (B) KEGG pathway for DEGs. The Bonferroni correction method was used for multiple hypothesis test correction and FDR-corrected P < 0.05 as cut-off. (JPEG 1050 kb)
Selected GO terms significantly enriched for upregulated genes in SD versus ND. (DOCX 19 kb)
Selected GO terms significantly enriched for downregulated genes in SD versus ND. (DOCX 20 kb)
Downregulated genes involved in the TOR signaling pathway in response to the sugar-only diet. The Bonferroni correction method was used for multiple hypothesis test correction and FDR-corrected P < 0.05 as cut-off. (JPEG 412 kb)
Detailed information on the DEGs of chronic SD and ND diets in Bactrocera dorsalis. (XLSX 556 kb)
About this article
Cite this article
Chen, EH., Hou, QL., Wei, DD. et al. Phenotypes, antioxidant responses, and gene expression changes accompanying a sugar-only diet in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae). BMC Evol Biol 17, 194 (2017). https://doi.org/10.1186/s12862-017-1045-5