Metabolic and transcriptome responses of RNAi-mediated AMPKα knockdown in Tribolium castaneum

doi:10.21203/rs.3.rs-24958/v1

Download PDF

Research article

Metabolic and transcriptome responses of RNAi-mediated AMPKα knockdown in Tribolium castaneum

https://doi.org/10.21203/rs.3.rs-24958/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 23 Sep, 2020

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

Background

The AMP-activated protein kinase (AMPK) is an intracellular fuel sensor for lipid and glucose metabolism. In addition to the short-term regulation of metabolic enzymes by phosphorylation, AMPK may also exert long-term effects on the transcription of downstream genes through the regulation of transcription factors and coactivators. In this study, RNA interference (RNAi) was conducted to investigate the effects of knockdown of TcAMPKα on lipid and carbohydrate metabolism in the red flour beetle, Tribolium castaneum, and the transcriptome profiles of dsTcAMPKα-injected and dsEGFP-injected beetles under normal conditions were comparatively analyzed by RNA-sequencing.

Results

RNAi-mediated suppression of TcAMPKα increased whole-body triglyceride (TG) level and the ratio between glucose and trehalose, as was confirmed by in vivo treatment with the AMPK-activating compound, 5-Aminoimidazole-4-carboxamide1-β-D-ribofuranoside (AICAR). A total of 1184 differentially expressed genes (DEGs) were identified between dsTcAMPKα-injected and dsEGFP-injected beetles. These include genes involved in lipid and carbohydrate metabolism as well as insulin/insulin-like growth factor signaling (IIS). Real-time quantitative polymerase chain reaction analysis confirmed the differential expression of selected genes. Interestingly, metabolism-related transcription factors such as sterol regulatory element-binding protein 1 (SREBP1) and carbohydrate response element-binding protein (ChREBP) were also significantly upregulated in dsTcAMPKα-injected beetles.

Conclusions

AMPK plays a critical role in the regulation of beetle metabolism. The findings of DEGs involved in lipid and carbohydrate metabolism provide valuable insight into the role of AMPK signaling in the transcriptional regulation of insect metabolism.

Epigenetics & Genomics

Tribolium castaneum

AMPK

RNA interference

transcriptome

lipid metabolism

carbohydrate metabolism IIS pathway

The survival of all organisms depends on the maintenance of energy homeostasis. AMP-activated protein kinase (AMPK) is a cellular energy sensor conserved across all eukaryotic species (Garcia and Shaw 2017). As a serine/threonine protein kinase complex, AMPK consists of a catalytic subunit α and two regulatory subunits, β and γ, and is activated in response to energy stress by sensing increases in ADP/ATP and AMP/ATP ratios, which leads to the activation of ATP-generating catabolic pathways and the inhibition of ATP-consuming anabolic pathways (Hardie 2007). While nucleotide-dependent phosphorylation of Thr172 in the α subunit by liver kinase B1 (LKB1) is the principal event required for full activation of AMPK (Hawley et al. 2003; Woods et al. 2003), several studies have revealed the nucleotide-independent regulation of AMPK via the phosphorylation of Thr172 by calcium/calmodulin-dependent kinase kinase 2 (CAMKK2) (Hawley et al. 2005; Hurley et al. 2005; Woods et al. 2005).

Lipids and carbohydrates are major sources for energy storage and supply in cells. It is well known that AMPK regulates lipids and carbohydrates homeostasis via direct phosphorylation of multiple downstream effectors. Both acetyl-CoA carboxylase (ACC), the first rate-limiting enzyme in fatty acid synthesis (FAS), and 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), the rate-limiting enzyme of cholesterol biosynthesis pathway, can be inhibited by AMPK through phosphorylation (Carling et al. 1987; Winder and Hardie 1996). Subsequently, AMPK was found to promote lipids absorption and release by directly phosphorylating lipases like hormone-sensitive lipase (HSL) and adipocyte-triglyceride lipase (ATGL) (Watt et al. 2006; Narbonne and Roy 2009; Ahmadian et al. 2011). On the other hand, AMPK could regulate glycolysis via the phosphorylation and activation of 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFKFB3) and 6-phosphofructo-2-kinase (PFK2) (Marsin et al. 2000; Marsin et al. 2002). Additionally, AMPK activation was also shown to increase whole-body insulin sensitivity by phosphorylation of insulin receptor substrate-1 (IRS-1) (Jakobsen et al. 2001) and insulin receptor (InR) (Chopra et al. 2012). Given the functional attributes of AMPK in lipids and carbohydrates metabolism, AMPK is considered as an important therapeutic target for treating metabolic diseases including obesity and type 2 diabetes (Musi and Goodyear 2002).

Although the role of AMPK in the regulation of cell metabolism is well studied in mammals, related research is still limited in insects. Notably, AMPK can also contribute to energy balance via modulation of transcriptional expression of metabolic enzymes in the long term, however, its downstream transcriptional pathways remains largely elusive (Cantó and Auwerx 2010). Recently, we reported the transcriptional and post-translational activation of TcAMPKα by oxidative, heat and cold stresses in the red flour beetle, Tribolium castaneum (Jiang et al. 2019). In this study, RNAi was employed to determine the roles of TcAMPKα in lipid and carbohydrate metabolism. Comparison, annotation and classification of DEGs between dsTcAMPKα treatment and control groups were also conducted by high-throughput transcriptome sequencing to identify metabolism-related genes modulated by AMPK in T. castaneum.

Effects of TcAMPKα suppression on TG, glucose and trehalose levels

RNAi was conducted to determine the effects of TcAMPKα knockdown on TG, glucose and trehalose levels. The injection of 20-day-old larvae with dsTcAMPKα reduced transcription levels by 95.50% ± 1.86% on the sixth day after injection. TG measurement showed that the TG level in dsTcAMPKα group (9.03 ± 0.76 mmol/mgprot) was significantly increased when compared with the dsEGFP group (5.96 ± 0.78 mmol/mgprot) (Fig. 1A).

Similarly, increased glucose level was observed in the beetles injected with dsTcAMPKα (25.85 ± 6.12 µmol/g) when compared to the control beetles (8.32 ± 0.08 µmol/g) (Fig. 1B). However, the trehalose level in dsTcAMPKα group (3.73 ± 0.10 mg/g) was significantly lower than that in dsEGFP group (4.08 ± 0.05 mg/g) (Fig. 1C). These data suggested that RNAi of TcAMPKα increased TG production and the ratio between glucose and trehalose (Fig. 2A).

Effects of AICAR treatment on TG, glucose and trehalose levels

To confirm the RNAi results, the 20-day-old larvae were treated with 5-Aminoimidazole-4-carboxamide1-β-D-ribofuranoside (AICAR), and the TG, glucose and trehalose levels were measured. The results showed that the levels of TG and glucose in AICAR group were significantly decreased by 34.60% ± 5.74% and 41.89% ± 2.27%, respectively, compared with injection buffer (IB) group, whereas the trehalose level increased by 17.07% ± 4.02% in beetles treated with AICAR. These data suggested that activation of TcAMPK decrased TG production and the ratio between glucose and trehalose (Fig. 1D-F; Fig. 2B).

Transcriptome Sequence And Reads Mapping

The dsTcAMPKα and dsEGFP groups were analyzed by RNA-Seq, respectively. A mean of 23,570,938 clean reads were generated among six independent libraries (T01-T06) (Table 1). Evaluation of clean data quality showed that the GC counts were ranged from 42–45% and Q30 ratios were > 93%, indicating a high level of data quality. The alignment of clean reads to the reference genome database of T. castaneum showed that 83.50% and 77.15% reads of the dsEGFP and dsTcAMPKα groups were aligned on average, respectively (Table 2).

Table 1

Summary of the transcriptome sequencing data from the controls and dsTcAMPKα treated samples
Samples	ID	Clean Read Number	Clean Base Number	GC (%)	Q30 (%)
EGFP1	T01	21,746,430	6,523,929,000	42.79	93.39
EGFP2	T02	20,376,206	6,112,861,800	42.86	93.94
EGFP3	T03	26,557,829	7,967,348,700	43.85	92.88
dsTcAMPKα1	T04	24,192,194	7,257,658,200	44.91	93.25
dsTcAMPKα2	T05	23,215,801	6,964,740,300	45.01	92.85
dsTcAMPKα3	T06	25,337,170	7,601,151,000	45.54	92.75

Table 2

Summary of average read numbers based on the RNA-sequencing data
	dsTcAMPKα	dsEGFP
Total alignments	38,451,597	39,074,583
Reads aligned	37,340,338	37,891,202
Unique alignments	36,870,045	37,564,766
Not aligned	11,156,438	7,895,775

Correlation analysis was conducted with Pearson’s Correlation Coefficient (R value) to evaluate the biological repeatability based on expression values of each library (Schulze et al. 2012). The results showed that the Pearson’s Correlation Coefficient between the three control samples (T01-T03) was 0.81, and that of dsTcAMPKα treatment samples (T04-T06) was 0.94 (Fig. 3A), whereas a relative low R value (0.78) was found between control and treatment groups. Box plot analysis revealed that the three samples in each group had similar expression distributions of reads, while the control and treatment groups had significantly different expression distributions (Fig. 3B).

Annotation Of Expressed Unigenes

To obtain the annotation information and functional classification, several nucleotide and protein databases were used for alignment using BLAST including Non-redundant (Nr), EuKaryotic Orthologous Groups (KOG), Clusters of Orthologous Groups of proteins (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), Protein family (Pfam), Gene Ontology (GO) and Swiss-Prot databases. Of the 31,944 unique sequences, a total of 14,095 unigenes were annotated and classified into at least one database (Fig. 4A).

GO is an international standard classification system of gene function. To classify the functions of the predicted T. castaneum unigenes, the GO enrichment analysis was carried out. The results showed that a total of 9,142 unigenes retrieved the associated GO terms, and were classified into three main GO categories. Of which, the majority of the unigenes were assigned to category of molecular function (13,229), followed by the biological process (13,083) and the cellular component (6344) (Fig. 4B, Figure S1 and Table S2). It is noteworthy that the metabolic process (4,635 unigenes) was the largest subcategory in biological process category, which account for 35.43% of this category (Fig. 4B, Figure S1 and Table S2).

KOG database is based on the phylogenetic relationship of bacteria, algae and eukaryotes. Direct homologous classification of unigenes can be carried out by using KOG database. Our assembled unigenes were subjected to KOG classification to predict their possible functions. A total of 9,144 unigenes were clustered into 25 categories. Among these categories, “General function prediction only” and “Signal transduction mechanisms” are the largest two groups, the number of annotated unigenes were 1,940 and 1,302, respectively. In addition, 911 unigenes were assigned to the categories of “lipid and amino acid transport and metabolism” (Fig. 4C and Table S2).

The KEGG database provides a powerful tool to discover the pathways in which genes are involved. There were 4,644 unigenes classified into five main categories and 22 subcategories. Among them, the global metabolism (23.38%) and “Translation” (24.88%) were the most abundant subcategories (Figure S2).

Changes in gene expression profiles

To identify the effects of knock-down of TcAMPKα on global unigene expression patterns of T. castaneum, DEGs between dsTcAMPKα and dsEGFP groups were identified based on their Fragments Per Kilobase of transcript sequence per Million base pairs sequenced (FPKM) values. The number of 1184 DEGs were obtained including 349 upregulated and 835 downregulated unigenes (Figure 5A, Figure 5B and Table S3). The log₂-fold variation range of DEGs was between -6.07 and 3.75.

DEGs were classified by searching against GO and KOG databases. GO term enrichments was used to further analyze the functional orientation of DEGs in T. castaneum. GO enrichments revealed that the DEGs involved in cellular component category were enriched in cell part, cell, membrane and organelle (Figure S3). In molecular function category, most of the DEGs were dramatically enriched in binding and catalytic activity such as Retrovirus-related Pol polyprotein and ATPase inhibitor, while metabolic and cellular process were the most enriched subcategories in biological process (Figure S3).

KOG database was used to annotate DEGs with specific physiological function (Figure 6, Table S4). Although the “General function prediction only” was the most abundant group among all subcategories in KOG database, a large amount of DEGs were annotated into the metabolism classifications such as “Lipid transport and metabolism”, “Carbohydrate transport and metabolism”, “Amino acid transport and metabolism”, “Energy production and conversion” and “Secondary metabolites biosynthesis, transport and catabolism” (Table 3). Furthermore, most DEGs involved in protein translation were downregulated such as some ribosomal proteins (Table 3). Similarly, among the 42 DEGs in the post-translational modification class, 35 DEGs were down-regulated, including 3 heat shock proteins (Hsps), whereas the phosphatidylinositol 4,5-bisphosphate 3-kinase (PIK3) and InR2 involved in the signal transduction mechanisms were upregulated.

Expressions of genes involved in lipid metabolism, carbohydrate metabolism and insulin signaling

The DEGs involved in lipid and carbohydrate metabolism were identified by transcriptome sequencing analysis. The results showed that the insect adipose triacylglycerol lipase homologue, brummer, which was involved in the lipid metabolism pathway, was significantly downregulated with the expression log₂-fold change of -2.20, whereas two FAS genes (FAS1-2) and the transcription factor ChREBP, a key regulator of glucose metabolism and fat storage (Uyeda et al. 2002), were upregulated with the expression log₂-fold change from 1.02 to 1.76 (Table 4). Knock-down of TcAMPKα also caused upregulation of genes involved in IIS pathway, including PI3K, IRS1 and InR2. To confirm the reliability of the DEG data, the expression levels of these DEGs were determined using RT-qPCR (Figure 7 and Table 4). Gene expression levels validated by RT-qPCR showed the high consistency with transcriptome sequencing.

To be as inclusive as possible, less strict screening criteria was applied, and additional lipid and carbohydrate metabolism-related genes with the expression changed in transcriptome were selected for further verification using RT-qPCR. The results showed that three FAS (FAS3-5), one ACC, one glycerol-3-phosphate acyltransferases (GPAT), four trehalases (TRE), carboxylase (PC) and phosphoenolpyruvate carboxykinase (PEPCK) were significantly upregulated (Table 4). Interestingly, SREBP1, the master regulator of lipid homeostasis, and SREBP cleavage-activating protein (SCAP), a central regulator of lipogenesis that controls the activity of SREBP (Shao and Espenshade 2014) were also significantly upregulated (Table 4).

Table 3. DEGs involved in different functional categories of KOG database

Categories	Up-regulated DEGs		Down-regulated DEGs
Categories	NO.	Partial gene description	NO.	Partial gene description
Signal transduction mechanisms	26	tyrosine-protein phosphatase; phorbol ester/diacylglycerol-binding protein; lachesin； arf-GAP; citron Rho-interacting kinase; cyclic nucleotide-gated cation channel; E3 ubiquitin; sortilin; lipid transfer protein	23	neurogenic locus protein; extensin; SNF1A/AMP-activated protein kinase; cGMP-dependent phosphodiesterase; troponin C; sortilin; atrial natriuretic; sensory neuron membrane protein; regulator complex protein; tetra phosphatase
Translation, ribosomal structure and biogenesis	1	eukaryotic translation initiation factor	42	ribosomal protein; exosome complex component; ribonuclease; H/ACA ribonucleoprotein; eukaryotic translation initiation factor
Posttranslational modification, protein turnover, chaperones	7	papilin; heat shock protein 68a; pregnancy zone protein; fucosyltransferase; brain tumor protein; E3 ubiquitin;	35	mannosyltransferase; sulfotransferase; GILT-like protein; heat shock protein 23; ubiquitin; suppressor protein; collagenase; glutathione S-transferase; protein transport protein; NEDD8
Amino acid transport and metabolism	4	hydroxylase; protease; transporter; glucose dehydrogenase	30	proteinase; Carboxypeptidase; glucose dehydrogenase; trypsin
Lipid transport and metabolism	5	nose resistant to fluoxetine protein; fatty acyl-CoA reductase; apolipophorins; ATP-binding cassette; Fatty acid synthetase	17	dehydrogenase/reductase; NADPH; alpha-tocopherol; acyl-CoA-binding protein; nose resistant to fluoxetine protein; fatty acyl-CoA reductase; desaturase; monooxygenase
Energy production and conversion	1	titin isoform X4	23	cytochrome b-c1 complex; stunted; ATP synthase; V-type proton ATPase; NADH dehydrogenase; cytochrome c oxidase; cytochrome b5; acylphosphatase;
Secondary metabolites biosynthesis, transport and catabolism	5	fatty acyl-CoA reductase; cytochrome P450s; multidrug resistance-associated protein lethal	12	dehydrogenase/reductase; NADPH; laccase; fatty acyl-CoA reductase; cytochrome P450s
Carbohydrate transport and metabolism	1	mucin	6	amylase; chitinase; lactoylglutathione lyase; glycolipid transfer protein; myrosinase; peritrophic matrix protein

Table 4. DEGs encoding metabolism related proteins and transcription factors/ co-activators from T. castaneum responding to dsTcAMPKα treatment.

Unigene name	padj	Description (blast)	Length (ORF bp)	Log₂Ratio Transcriptome	Log₂Ratio qRT-PCR
Lipid metabolism
FAS1	0.038775298	fatty acid synthase	12981	1.02	0.37
FAS2	0.056982009	fatty acid synthase	6522	1.76	0.71
FAS3	0.260213539	fatty acid synthase	7152	0.79	1.05
FAS4	0.425672932	fatty acid synthase	6630	0.63	1.16
FAS5	0.170696071	fatty acid synthase	6450	0.58	1.53
ACC	0.186316184	acetyl-CoA carboxylase	7005	0.83	0.10
GPAT3	0.046646995	glycerol-3-phosphate acyltransferases	1440	0.67	1.51
Brummer	1.60E-05	triacylglycerol lipase	1635	-2.20	-0.37
Carbohydrate metabolism
TRE1-1	0.036902117	Trehalase1-1	1662	0.89	1.46
TRE1-3	0.450009964	Trehalase1-3	>507	0.56	0.58
TRE1-4	0.007152475	Trehalase1-4	1812	0.68	0.59
TRE2	0.036902117	Trehalase2	1647	0.89	1.61
Insulin signaling pathway
IRS1	0.018687502	insulin receptor substrate	2760	1.13	0.63
InR2	0.004419643	insulin-like receptor	4185	1.09	0.23
PI3K	3.37E-08	phosphatidylinositol 4,5-bisphosphate 3-kinase	3186	1.02	0.63
Transcription factor and co-activator
SCAP	0.031742692	sterol regulatory element binding protein cleavage-activating protein	3783	0.81	1.12
SREBP1	0.036908768	sterol regulatory element binding protein 1	3078	0.94	1.74
ChREBP	0.000606393	carbohydrate response-element-binding protein	510	1.18	1.94

In insects, the energy for growth and development is mainly derived from the storage and utilization of lipids and carbohydrates in specific tissues such as fat body, midgut and oenocytes (Martins et al. 2011). TG is the master form of lipids stored in fat body and plays an important role in energy storage and release (Warnakula et al. 2011). Meanwhile, glucose (monosaccharide) and trehalose (disaccharide) provide energy through glycolysis and gluconeogenesis (Reyes-DelaTorre et al. 2012). Given that AMPK activates ATP-generating pathways and inhibits energy-consuming processes under conditions of low energy status (Towler and Hardie 2007), the attenuation of AMPK signaling in insects may disrupt the energy balance in vivo. In this study, we found that the TG and glucose levels in dsTcAMPKα-injected T. castaneum were significantly increased compared with dsEGFP-injected insects, while a decreased trehalose content was observed in dsTcAMPKα-injected insects. These results were further confirmed by in vivo AICAR treatment, which resulted in the decreased TG and glucose levels and increased trehalose content. Similarly, AICAR-induced AMPK activation resulted in significantly decreased TG level in lean and obese rodents, in vivo (Bergeron et al. 2001). In human HepG2 cells, the kinase-inactive AMPKα increased lipid content and prevent the metformin from decreasing lipid accumulation (Zang et al. 2004). However, lower TG levels were observed in D. melanogaster with reduced AMPK function during fed conditions (Johnson et al. 2010). On the other hand, while activation of AMPK triggered a reduction in glucose levels in vivo in mice (Yamauchi et al. 2002), and increased hepatic glucose production was observed in AMPKα2 knockout mice (Andreelli et al. 2006), knockdown of hypothalamic AMPK activity in male Sprague-Dawley rats led to a significant suppression of glucose production (Yang et al. 2010). These results indicate the complexity of the role of AMPK in the regulation of metabolic processes.

As an energy sensor that regulates cellular metabolism, AMPK not only has acute effects on metabolic enzymes by direct phosphorylation, but also shows long term action to change the transcriptional levels of metabolic proteins and enzymes. For example, activation of AMPK in liver and adipocytes can downregulate FAS activity and inhibit lipogenesis (Madsen et al. 2015). Activation of AMPK by 4-Hydroxyderricin and xanthoangelol downregulates GPAT in 3T3‐L1 cells (Zhang et al. 2013), an enzyme necessary for triacylglycerol synthesis (Wendel et al. 2009). Similarly, in C. elegans, AAK-2 can inhibit fat synthesis under stress condition by downregulating lipid synthesis-related genes such as Δ9 fatty acid desaturases which can produce monounsaturated fatty acids to constitute TGs (Shin et al. 2011; Watts 2009). To investigate the downstream transcriptional pathways of AMPK in T. castaneum, the transcriptomes of dsTcAMPKα-injected and dsEGFP-injected larvae under normal conditions were comparatively analyzed. Global gene expression profiles of dsTcAMPKα group were distinct from dsEGFP group with 1184 DEGs. Of particular interest is the upregulation of genes involved in lipid synthesis, such as FAS, ACC and GPAT, whereas the relative expression level of brummer was dramatically downregulated. We also observed the upregulation of several carbohydrate metabolism-related genes, such as TRE, and two key enzymes in gluconeogenesis, PC and PEPCK (Jitrapakdee et al. 2008; Hanson and Reshef 1997), in the dsTcAMPKα-injected insects. The transcriptional changes of genes encoding metabolism enzymes might contribute to the increased TG and glucose levels and decreased trehalose content in dsEGFP-injected beetles.

AMPK activation has been reported to cause a reduction in transcriptional activity of several metabolism-related transcription factors. As a key transcription factor that regulates cellular lipogenesis in liver, skeletal muscle and adipose tissue, SREBP1 acts in synergy with ChREBP, which mediates the response to dietary carbohydrates (Postic and Ortega-Prieto 2019). AMPK-inhibited SREBP1 can block expression of some gluconeogenic and lipogenic genes, whereas SREBP1 overexpression can increase their transcription (Shimomura et al. 1999; Foretz et al. 1999; Kohjima et al. 2008). Activation of AMPK by metformin or an adenosine analogue suppresses the expression of SREBP1. In metformin-treated rats, hepatic expression of SREBP1 mRNAs and protein is reduced (Zhou et al. 2001). A recent study indicated that AMPK interacts with and directly phosphorylates SREBP1, suppresses SREBP1 cleavage and nuclear translocation, and represses SREBP1 target gene Fas expression in hepatocytes in response to metformin treatment, leading to reduced lipogenesis (Li et al. 2011b). On the other hand, AMPK also phosphorylates Ser568 and reduce DNA binding capacity and promote nuclear exclusion of murine ChREBP (Kawaguchi et al. 2001; Kawaguchi et al. 2002). In this study, SREBP1, ChREBP and SCAP were significantly upregulated in dsTcAMPKα-injected insects, which in turn might modulated the expression of genes involved in lipid and carbohydrate metabolism.

IIS pathway was involved in the regulation of glucose and lipid metabolism (Saltie and Kahn 2001). In addition to the regulation of lipid synthesis, studies reveal expanding roles for SREBP1 in controlling pathways for insulin resistance (Li et al. 2011b), in which the pathological process involves a series of cascades, including defective activation of IRS and PI3K (Kerouz et al. 1997; Huang et al. 2005; Taniguchi et al. 2005). Overexpression of SREBP1 decreased Irs-1 mRNA levels in a dose-dependent manner, and SREBP1 knockdown led to an upregulation of IRS-1 levels (Li et al. 2011b). Further luciferase reporter assay confirmed that Irs-1 promoter activity was repressed by SREBP1 overexpression (Li et al. 2011b). However, in liver with Nonalcoholic fatty liver disease (NAFLD), IRS-1 expression was enhanced and correlated positively with SREBP1 expression. In contrast, IRS-2 expression decreased by 50% and was not correlated with SREBP1 (Kohjima et al. 2008). In this study, we observed the upregulation of IIS-related DEGs like IRS1, InR2 and PI3K in dsTcAMPKα-injected beetles. Interestingly, it has been reported that knockdown in expression of insulin like peptide 2 (ILP2) caused a decrease in TRE mRNA levels in T. castaneum (Xu et al. 2013), suggesting the positive control of TRE transcription by IIS signaling. Further study is needed to clarify the mechanism of upregulation of IIS signaling in dsTcAMPKα-injected insects.

This study confirmed that AMPK has an important role in the regulation of beetle metabolism. Specifically, our study showed that knockdown of AMPK causes alterations in expression levels of genes involved in lipid and carbohydrate metabolism as well as IIS signaling. Such investigations will help us understand the function of AMPK in transcriptional regulation of insect metabolism.

Experimental insects

The Georgia-1 (GA-1) strain of T. castaneum was reared at 30 °C and 50% relative humidity in 5% yeasted flour under standard conditions as described previously (Haliscak and Beeman 1983; Li et al. 2011a).

Double-strand RNA Synthesis And Injection

Gene specific primers (Table S1) with T7 promoter were designed to synthesize the dsRNAs targeting nucleotides 844–1285 (442 bp) of the ORF region of the TcAMPKα using TranscriptAid™ T7 High Yield Transcription Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer's instructions. The synthesized dsRNAs were diluted in diethyl pyrocarbonate (DEPC)-treated water with a concentration of 2 µg/µL, and about 200 ng of dsRNA in 200 nL IB was injected into 20-day-old larvae using a Nanoliter 2010 injector system (WPI, Sarasota, FL, USA) under a stereomicroscope. A total of 10 unsexed insect larvae were collected as a sample for RNA extraction at six- days post-injection. The dsEGFP-injected larvae (CK group) were set as controls in all injection experiments. Three biologically independent replicates were carried out with at least 100 insects (≥ 200 mg) in each replicate.

Triglyceride (TG) Measurement

Total TG levels were determined using the liquid TG (GPO-PAP) method (Annoni et al. 1982). Briefly, each replicate with 10 injected larvae were homogenized in 270 µL of PBS (0.1 mol/L PH 7.4) and centrifuged at 2,500 rpm for 10 min. The supernatant was collected, and the TG level was analyzed by using TG Assay Kit (catalogue no. A110-1, Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer's instructions. Three independent biological replicates and four technical replicates were performed for every treatment.

Glucose And Trehalose Measurement

Total glucose level was measured using glucose content assay kit (catalogue no. BC2500, Solarbio Science & Technology, Beijing, China) according to the manufacturer's instructions. Briefly, 10 larvae on day 6 after injection of dsTcAMPKα or dsEGFP were weighed and homogenized in 0.3 mL double distilled water. The supernatant was collected after centrifugation at 8,000 g for 10 minutes, and used to quantify the level of glucose at 505 nm with Multiskan FC microplate reader (Thermo Fisher Scientific, Waltham, MA, USA). Three independent biological replicates and 4 technical replicates were performed for every treatment.

Total trehalose level was measured using Trehalose content detection kit (catalogue no. A149-1-1, Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer's instructions. Six days after injection of dsTcAMPKα or dsEGFP, 10 larvae were homogenized in 0.3 mL extraction solution and incubated at room temperature for 45 minutes. The supernatant was collected after centrifugation at 8,000 g for 10 minutes, and used to quantify the level of trehalose at 620 nm with Multiskan FC microplate reader (Thermo Fisher Scientific, Waltham, MA, USA). Three independent biological replicates and four technical replicates were performed for every treatment.

AICAR Treatment

To investigate the effects of AMPK activation on lipids and carbohydrates metabolism, 20-day-old larvae were injected with 200 nL of 1.6 mg/mL AICAR or IB as control. AMPK could be activated immediately by AICAR treatment among a couple of hours (Matsumoto et al. 2011). The total amounts of TG, glucose and glucose were determined at one- hour post-injection.

RNA Extraction, Library Construction And Sequencing

Total RNA was isolated from 10 larvae of dsTcAMPKα or dsEGFP group on the sixth day after injection using TRIzol Reagent (Invitrogen, USA), and digested by RNase-free DNase I (Takara, Dalian, China) to remove genomic DNA contaminants. To ensure the qualified samples for transcriptome sequencing, concentration and integrity of RNA samples were checked using a Nanodrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE) and an Agilent Bioanalyzer 2100 (Agilent, Palo Alto, CA), respectively. The qualified RNA samples were used for mRNA preparation and cDNA library construction. Three independent biological replicates were performed for every treatment. cDNA libraries were sequenced as described by Meng et al. (2019).

Sequence Annotation

To ensure the accuracy of sequence alignment, raw reads were cleaned by removing adapter and primer sequences, reads with ambiguous nucleotides larger than 5% and low-quality reads. Clean reads were aligned to the reference genome of T. castaneum (Reference genome version: GCF_000002335.3_Tcas5.2_genomic.fna; Reference genome site: ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/002/335/GCF_000002335.3_Tcas5.2/GCF_000002335.3_Tcas5.2_genomic.fna.gz) using HISAT2 (CCB, Johns Hopkins University, USA). Subsequently, the transcripts were subjected to blast against Nr, KOG, COG, KEGG, Pfam, GO and Swiss-Prot databases, and the unigenes were annotated. The E-value threshold was set to 10^− 5.

Differentially Expressed Genes Analysis

Transcript abundances were measured by FPKM (Trapnell et al. 2010). Differential expression analysis of dsTcAMPKα and dsEGFP groups were performed using DESeq2 software (Anders and Huber 2010). Unigene with False Discovery Rate (FDR) values < 0.05 and fold change |log₂(treatment/control) | > 1 were set as the threshold for DEGs. GO and KOG database annotations were further analyzed to understand the functions of DEGs in T. castaneum.

Reverse Transcription Quantitative Real-time Polymerase Chain Reaction (RT-qPCR)

To validate the DEGs from the RNA-sequencing, RT-qPCR reactions were performed on the Bio-Rad CFX 96 Real-time PCR system using TB Green™ Premix Ex Taq™ (Takara, Dalian, China) and gene specific primers (Table S1). The stably expressed gene encoding ribosomal protein S3 (rps3, GenBank: CB335975) was used as a reference gene (Arakane et al. 2008). PCR conditions were set as an initial incubation of 95℃ for 30 s, 40 cycles of 95 °C for 5 s and 60 °C for 30 s, and a final melting curve analysis was performed. The mRNA levels were normalized to reference gene with the 2^−ΔΔCT method (Livak and Schmittgen 2001). The means and standard errors (mean ± SE) for each time point were obtained from the average of at least three biologically independent sample sets.

Statistical analysis

Statistical analysis was performed using Graphpad Prism 6 (GraphPad Software Inc., San Diego, USA) by one-way analysis of variance, followed by Tukey’s Honestly Significant Difference test. All data are presented as the mean ± SE.

ACC

acetyl-CoA carboxylase

AICAR

5-Aminoimidazole-4-carboxamide1-β-D-ribofuranoside

ATGL

adipocyte-triglyceride lipase

AMPK

AMP-activated protein kinase

CAMKK2

calcium/calmodulin-dependent kinase kinase 2

CHREBP

carbohydrate response element-binding protein

COG

Clusters of Orthologous Groups of proteins

DEGs

differentially expressed genes

DEPC

diethyl pyrocarbonate

FAS

fatty acid synthesis

FDR

False Discovery Rate

FPKM

Fragments Per Kilobase of transcript sequence per Million base pairs sequenced

Gene Ontology

GPAT

glycerol-3-phosphate acyltransferases

HMGR

3-hydroxy-3-methylglutaryl-CoA reductase

HSL

hormone-sensitive lipase

Hsps

heat shock proteins

injection buffer

IIS

insulin/insulin-like growth factor signaling

ILP2

insulin like peptide 2

InR

insulin receptor

IRS-1

insulin receptor substrate-1

KEGG

Kyoto Encyclopedia of Genes and Genomes

KOG

EuKaryotic Orthologous Groups

LKB1

liver kinase B1

NAFLD

Nonalcoholic fatty liver disease

Non-redundant

carboxylase

PEPCK

phosphoenolpyruvate carboxykinase

Pfam

Protein family

PFK2

6-phosphofructo-2-kinase

PFKFB3

6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3

PIK3

phosphatidylinositol 4,5-bisphosphate 3-kinase

RNAi

RNA interference

SCAP

SREBP cleavage-activating protein

SREBP1

sterol regulatory element-binding protein 1

triglyceride

TRE

trehalase

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and material

All data generated or analyzed during this study are included in this article and its supplementary information files.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was supported by the National Natural Science Foundation of China under grant nos. 31572000 and 31871974.

The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Authors’ Contributions

JW conceived and designed research. HJ, NZ and CJ conducted experiments. XM, YZ and KQ analyzed the data. HJ and JW wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Ahmadian M, Abbott MJ, Tang T, Hudak CS, Kim Y, Bruss M, Hellerstein MK, Lee HY, Samuel VT, Shulman GI. Desnutrin/ATGL is regulated by AMPK and is required for a brown adipose phenotype. Cell Metab. 2011;13:739–48.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Andreelli F, Foretz M, Knauf C, Cani PD, Perrin C, Iglesias MA, Pillot B, Bado A, Tronche F, Mithieux G. Liver adenosine monophosphate-activated kinase-α2 catalytic subunit is a key target for the control of hepatic glucose production by adiponectin and leptin but not insulin. Endocrinology. 2006;147:2432–41.
Annoni G, Botasso BM, Ciaci D, Donato MF, Tripodi A. Liquid triglycerides (GPO-PAP). Medi. Diagnostic Italy. J Labor Clin Med. 1982;9:115.
Arakane Y, Specht CA, Kramer KJ, Muthukrishnan S, Beeman RW. Chitin synthases are required for survival, fecundity and egg hatch in the red flour beetle, Tribolium castaneum. Insect Biochem Mol Biol. 2008;38:959–62.
Bergeron R, Previs SF, Cline GW, Perret P, Russell RR III, Young LH, Shulman GI. Effect of 5-aminoimidazole-4-carboxamide-1-β-d-ribofuranoside infusion on in vivo glucose and lipid metabolism in lean and obese Zucker rats. Diabetes. 2001;50:1076–82.
Cantó C, Auwerx J. AMP-activated protein kinase and its downstream transcriptional pathways. Cell Mol Life Sci. 2010;67:3407–23.
Carling D, Zammit VA, Hardie DG. A common bicyclic protein kinase cascade inactivates the regulatory enzymes of fatty acid and cholesterol biosynthesis. FEBS Lett. 1987;223:217–22.
Chopra I, Li H, Wang H, Webster KA. Phosphorylation of the insulin receptor by AMP-activated protein kinase (AMPK) promotes ligand-independent activation of the insulin signalling pathway in rodent muscle. Diabetologia. 2012;55:783–94.
Foretz M, Guichard C, Ferré P, Foufelle F. Sterol regulatory element binding protein-1c is a major mediator of insulin action on the hepatic expression of glucokinase and lipogenesis-related genes. Proc. Natl. Acad. Sci. U. S. A. 1999; 96: 12737–12742.
Garcia D, Shaw RJ. AMPK: mechanisms of cellular energy sensing and restoration of metabolic balance. Mol Cell. 2017;66:789–800.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644.
Haliscak JP, Beeman RW. Status of malathion resistance in five genera of beetles infesting farm-stored corn, wheat, and oats in the United States. J Econ Entomol. 1983;76:717–22.
Hanson RW, Reshef L. Regulation of phosphoenolpyruvate carboxykinase (GTP) gene expression. Annu Rev Biochem. 1997;66:581–611.
Hardie DG. AMP-activated/SNF1 protein kinases: conserved guardians of cellular energy. Nat Rev Mol Cell Biol. 2007;8:774.
Hawley SA, Boudeau J, Reid JL, Mustard KJ, Udd L, Mäkelä TP, Alessi DR, Hardie DG. Complexes between the LKB1 tumor suppressor, STRADα/β and MO25α/β are upstream kinases in the AMP-activated protein kinase cascade. J Biol. 2003;2:28.
Hawley SA, Pan DA, Mustard KJ, Ross L, Bain J, Edelman AM, Frenguelli BG, Hardie DG. Calmodulin-dependent protein kinase kinase-β is an alternative upstream kinase for AMP-activated protein kinase. Cell Metab. 2005;2:9–19.
Huang TJ, Verkhratsky A, Fernyhough P. Insulin enhances mitochondrial inner membrane potential and increases ATP levels through phosphoinositide 3-kinase in adult sensory neurons. Mol Cell Neurosci. 2005;28:42–54.
Hurley RL, Anderson KA, Franzone JM, Kemp BE, Means AR, Witters LA. The Ca²⁺/calmodulin-dependent protein kinase kinases are AMP-activated protein kinase kinases. J Biol Chem. 2005;280:29060–6.
Jakobsen SN, Hardie DG, Morrice N, Tornqvist HE. 5′-AMP-activated protein kinase phosphorylates IRS-1 on Ser-789 in mouse C2C12 myotubes in response to 5-aminoimidazole-4-carboxamide riboside. J Biol Chem. 2001;276:46912–6.
Jiang H, Zhang N, Chen M, Meng X, Ji C, Ge H, Dong F, Miao L, Yang X, Xu X, Qian K, Wang J. Transcriptional and post-translational activation of AMPKα by oxidative, heat, and cold stresses in the red flour beetle, Tribolium castaneum. Cell Stress Chaperones. 2019;24:1079–89.
Jitrapakdee S, St Maurice M, Rayment I, Cleland WW, Wallace JC, Attwood PV. Structure, mechanism and regulation of pyruvate carboxylase. Biochem J. 2008;413:369–87.
Johnson EC, Kazgan N, Bretz CA, Forsberg LJ, Hector CE, Worthen RJ, Onyenwoke R, Brenman JE. Altered metabolism and persistent starvation behaviors caused by reduced AMPK function in Drosophila. PLoS One. 2010;5:e12799.
Kawaguchi T, Osatomi K, Yamashita H, Kabashima T, Uyeda K. Mechanism for fatty acid “sparing” effect on glucose-induced transcription regulation of carbohydrate-responsive element-binding protein by AMP-activated protein kinase. J Biol Chem. 2002;277:3829–35.
Kawaguchi T, Takenoshita M, Kabashima T, Uyeda K. Glucose and cAMP regulate the L-type pyruvate kinase gene by phosphorylation/dephosphorylation of the carbohydrate response element binding protein. Proc. Natl. Acad. Sci. U. S. A. 2001; 98: 13710–13715.
Kerouz NJ, Hörsch D, Pons S, Kahn CR. Differential regulation of insulin receptor substrates-1 and-2 (IRS-1 and IRS-2) and phosphatidylinositol 3-kinase isoforms in liver and muscle of the obese diabetic (ob/ob) mouse. J Clin Invest. 1997;100:3164–72.
Kohjima M, Higuchi N, Kato M, Kotoh K, Yoshimoto T, Fujino T, Yada M, Yada R, Harada N, Enjoji M. SREBP-1c, regulated by the insulin and AMPK signaling pathways, plays a role in nonalcoholic fatty liver disease. Int J Mol Med. 2008;21:507–11.
Li B, Beeman RW, Park Y. Functions of duplicated genes encoding CCAP receptors in the red flour beetle, Tribolium castaneum. J Insect Physiol. 2011a;57:1190–7.
Li Y, Xu S, Mihaylova MM, Zheng B, Hou X, Jiang B, Park O, Luo Z, Lefai E, Shyy YJ. AMPK phosphorylates and inhibits SREBP activity to attenuate hepatic steatosis and atherosclerosis in diet-induced insulin-resistant mice. Cell Metab. 2011b;13:376–88.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^–∆∆CT method. Methods. 2011;25:402–8.
Madsen A, Bozickovic O, Bjune JI, Mellgren G, Sagen JV. Metformin inhibits hepatocellular glucose, lipid and cholesterol biosynthetic pathways by transcriptionally suppressing steroid receptor coactivator 2 (SRC-2). Sci Rep. 2015;5:16430.
Marsin AS, Bouzin C, Bertrand L, Hue L. The stimulation of glycolysis by hypoxia in activated monocytes is mediated by AMP-activated protein kinase and inducible 6-phosphofructo-2-kinase. J Biol Chem. 2002;277:30778–83.
Marsin A, Bertrand L, Rider M, Deprez J, Beauloye C, Vincent MF, Van den Berghe G, Carling D, Hue L. Phosphorylation and activation of heart PFK-2 by AMPK has a role in the stimulation of glycolysis during ischaemia. Curr Biol. 2000;10:1247–55.
Martins GF, Serrão JE, Ramalho-OrtigãO JM, Pimenta PFP. Histochemical and ultrastructural studies of the mosquito Aedes aegypti fat body: effects of aging and diet type. Microsc Res Tech. 2011;74:1032–9.
Matsumoto Y, Sumiya E, Sugita T, Sekimizu K. An invertebrate hyperglycemic model for the identification of anti-diabetic drugs. PLoS One. 2011;6:e18292.
Musi N, Goodyear LJ. Targeting the AMP-activated protein kinase for the treatment of type 2 diabetes. Curr Drug Targets Immune Endocr Metabol Disord. 2002;2:119–27.
Narbonne P, Roy R. Caenorhabditis elegans dauers need LKB1/AMPK to ration lipid reserves and ensure long-term survival. Nature. 2009;457:210.
Postic C, Ortega-Prieto P. Carbohydrate sensing through the transcription factor ChREBP. Front Genet. 2019;10:472.
Reyes-DelaTorre A, Peña-Rangel MT, Riesgo-Escovar JR. Carbohydrate metabolism in Drosophila: reliance on the disaccharide trehalose, Carbohydrates-comprehensive studies on glycobiology and glycotechnology. 2012; IntechOpen.
Schulze SK, Kanwar R, Gölzenleuchter M, Therneau TM, Beutler AS. SERE: single-parameter quality control and sample comparison for RNA-SEq. BMC Genom. 2012;13:524.
Shao W, Espenshade PJ. Sterol regulatory element-binding protein (SREBP) cleavage regulates Golgi-to-endoplasmic reticulum recycling of SREBP cleavage-activating protein (SCAP). J Biol Chem. 2014;289:7547–57.
Shimomura I, Bashmakov Y, Horton JD. Increased levels of nuclear SREBP-1c associated with fatty livers in two mouse models of diabetes mellitus. J Biol Chem. 1999;274:30028–32.
Shin H, Lee H, Fejes AP, Baillie DL, Koo HS, Jones SJ. Gene expression profiling of oxidative stress response of C. elegans aging defective AMPK mutants using massively parallel transcriptome sequencing. BMC Res Notes. 2011;4:34.
Taniguchi CM, Ueki K, Kahn CR. Complementary roles of IRS-1 and IRS-2 in the hepatic regulation of metabolism. J Clin Invest. 2005;115:718–27.
Towler MC, Hardie DG. AMP-activated protein kinase in metabolic control and insulin signaling. Circ Res. 2007;100:328–41.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–5.
Uyeda K, Yamashita H, Kawaguchi T. Carbohydrate responsive element-binding protein (ChREBP): a key regulator of glucose metabolism and fat storage. Biochem Pharmacol. 2002;63:2075–80.
Warnakula S, Hsieh J, Adeli K, Hussain MM, Tso P, Proctor SD. New insights into how the intestine can regulate lipid homeostasis and impact vascular disease: frontiers for new pharmaceutical therapies to lower cardiovascular disease risk. Can J Cardiol. 2011;27:183–91.
Watt MJ, Dzamko N, Thomas WG, Rose-John S, Ernst M, Carling D, Kemp BE, Febbraio MA, Steinberg GR. CNTF reverses obesity-induced insulin resistance by activating skeletal muscle AMPK. Nat Med. 2006;12:541.
Watts JL. Fat synthesis and adiposity regulation in Caenorhabditis elegans. Trends Endocrinol Metab. 2009;20:58–65.
Wendel AA, Lewin TM, Coleman RA. Glycerol-3-phosphate acyltransferases: rate limiting enzymes of triacylglycerol biosynthesis. Biochim Biophys Acta Mol Cell Biol Lipids. 2009;1791:501–6.
Winder W, Hardie D. Inactivation of acetyl-CoA carboxylase and activation of AMP-activated protein kinase in muscle during exercise. Am J Physiol Endocrinol Metab. 1996;270:E299–304.
Woods A, Dickerson K, Heath R, Hong SP, Momcilovic M, Johnstone SR, Carlson M, Carling D. Ca²⁺/calmodulin-dependent protein kinase kinase-β acts upstream of AMP-activated protein kinase in mammalian cells. Cell Metab. 2005;2:21–33.
Woods A, Johnstone SR, Dickerson K, Leiper FC, Fryer LG, Neumann D, Schlattner U, Wallimann T, Carlson M, Carling D. LKB1 is the upstream kinase in the AMP-activated protein kinase cascade. Curr Biol. 2003;13:2004–8.
Xu J, Sheng Z, Palli SR. Juvenile hormone and insulin regulate trehalose homeostasis in the red flour beetle, Tribolium castaneum. PLoS Genet. 2013;9:e1003535.
Yamauchi T, Kamon J, Minokoshi Y, Ito Y, Waki H, Uchida S, Yamashita S, Noda M, Kita S, Ueki K. Adiponectin stimulates glucose utilization and fatty-acid oxidation by activating AMP-activated protein kinase. Nat Med. 2002;8:1288.
Yang CS, Lam CK, Chari M, Cheung GW, Kokorovic A, Gao S, Leclerc I, Rutter GA, Lam TK. Hypothalamic AMP-activated protein kinase regulates glucose production. Diabetes. 2010;59:2: 2435–43.
Zang M, Zuccollo A, Hou X, Nagata D, Walsh K, Herscovitz H, Brecher P, Ruderman NB, Cohen RA. AMP-activated protein kinase is required for the lipid-lowering effect of metformin in insulin-resistant human HepG2 cells. J Biol Chem. 2004;279:47898–905.
Zhang T, Sawada K, Yamamoto N, Ashida H. 4-Hydroxyderricin and xanthoangelol from Ashitaba (Angelica keiskei) suppress differentiation of preadiopocytes to adipocytes via AMPK and MAPK pathways. Mol Nutr Food Res. 2013;57:1729–40.
Zhou G, Myers R, Li Y, Chen Y, Shen X, Fenyk-Melody J, Wu M, Ventre J, Doebber T, Fujii N. Role of AMP-activated protein kinase in mechanism of metformin action. J Clin Invest. 2001;108:1167–74.

Download PDF

Journal Publication

published 23 Sep, 2020

Read the published version in BMC Genomics →

Editorial decision: Minor revision
21 Jul, 2020
Review #2 received at journal
14 Jul, 2020
Reviewer #2 agreed at journal
11 Jun, 2020
Review #1 received at journal
09 Jun, 2020
Reviewers invited by journal
18 May, 2020
Reviewer #1 agreed at journal
18 May, 2020
Editor assigned by journal
30 Apr, 2020
Submission checks completed at journal
29 Apr, 2020
Editor invited by journal
29 Apr, 2020
First submitted to journal
28 Apr, 2020

You are reading this older preprint version

Read the latest preprint version →

Metabolic and transcriptome responses of RNAi-mediated AMPKα knockdown in Tribolium castaneum

Status:

Journal Publication

Version 1

Abstract

Background

Results

Conclusions

Figures

Background

Results

Effects of TcAMPKα suppression on TG, glucose and trehalose levels

Transcriptome Sequence And Reads Mapping

Annotation Of Expressed Unigenes

Changes in gene expression profiles

Expressions of genes involved in lipid metabolism, carbohydrate metabolism and insulin signaling

Discussion

Conclusion

Methods

Experimental insects

Double-strand RNA Synthesis And Injection

Triglyceride (TG) Measurement

Glucose And Trehalose Measurement

AICAR Treatment

RNA Extraction, Library Construction And Sequencing

Sequence Annotation

Differentially Expressed Genes Analysis

Reverse Transcription Quantitative Real-time Polymerase Chain Reaction (RT-qPCR)

Statistical analysis

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1