# Differential effects of silencing crustacean hyperglycemic hormone gene expression on the metabolic profiles of the muscle and hepatopancreas in the crayfish Procambarus clarkii

Wenfeng Li, Kuo-Hsun Chiu, Yi-Chun Tien, Shih-Fu Shih, Li-Jane Shih, Chien-Hsun Lee, Jean-Yves Toullec, Chi-Ying Lee

<figure><img src="image_1.png" /></figure>

► To cite this version:

Wenfeng Li, Kuo-Hsun Chiu, Yi-Chun Tien, Shih-Fu Shih, Li-Jane Shih, et al.. Differential effects of silencing crustacean hyperglycemic hormone gene expression on the metabolic profiles of the muscle and hepatopancreas in the crayfish Procambarus clarkii. PLoS ONE, 2017, 12 (2), pp.e0172557. 10.1371/journal.pone.0172557.s004 . hal-01474112

HAL Id: hal-01474112

https://hal.sorbonne-universite.fr/hal-01474112

Submitted on 22 Feb 2017

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Differential effects of silencing crustacean hyperglycemic hormone gene expression on the metabolic profiles of the muscle and hepatopancreas in the crayfish *Procambarus clarkii*

Wenfeng Li¹, Kuo-Hsun Chiu², Yi-Chun Tien¹, Shih-Fu Tsai¹, Li-Jane Shih³, Chien-Hsun Lee¹, Jean-Yves Toullec⁴,⁵, Chi-Ying Lee¹*

1 Department of Biology, National Changhua University of Education, Changhua, Taiwan, 2 Department of Aquaculture, National Kaohsiung Marine University, Kaohsiung, Taiwan, 3 Department of Medical Laboratory, Taoyuan Armed Forces General Hospital, Taoyuan, Taiwan, 4 Sorbonne Universités, UPMC Université Paris 06, UMR 7144 CNRS, Equipe ABICE, Station Biologique de Roscoff, Roscoff, France, 5 CNRS, UMR 7144, Adaptation et Diversité en Milieu Marin, Station Biologique de Roscoff, Roscoff, France

* bicylee@cc.ncue.edu.tw

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Li W, Chiu K-H, Tien Y-C, Tsai S-F, Shih L-J, Lee C-H, et al. (2017) Differential effects of silencing crustacean hyperglycemic hormone gene expression on the metabolic profiles of the muscle and hepatopancreas in the crayfish *Procambarus clarkii*. PLoS ONE 12(2): e0172557. doi:10.1371/journal.pone.0172557

**Editor:** Emanuele Buratti, International Centre for Genetic Engineering and Biotechnology, ITALY

**Received:** December 12, 2016

**Accepted:** February 6, 2017

**Published:** February 16, 2017

**Copyright:** © 2017 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The present study was supported financially by Ministry of Science and Technology, Taiwan (https://www.most.gov.tw/) through grants (102-2311-B-018-002 and 104-2311-B-018-002) to C-YL. WL was supported by a postdoctoral fellowship (104-2811-B-018-001) from Ministry of Science and Technology, Taiwan. The funders had

## Abstract

In order to functionally characterize the metabolic roles of crustacean hyperglycemic hormone (CHH), gene expression of CHH in the crayfish (*Procambarus clarkii*) was knocked down by *in vivo* injection of CHH double-stranded RNA (dsRNA), followed by metabolomic analysis of 2 CHH target tissues (the muscle and hepatopancreas) using nuclear magnetic resonance spectroscopy. Compared to the levels in untreated and saline-injected (SAI) animals, levels of CHH transcript, but not those of molt-inhibiting hormone (a CHH-family peptide), in the eyestalk ganglia of CHH dsRNA-injected (DSI) animals were significantly decreased at 24, 48, and 72 hour post injection (hpi), with concomitant changes in levels of CHH peptide in the sinus gland (a neurohemal organ) and hemolymph. Green fluorescence protein (GFP) dsRNA failed to affect levels of CHH transcript in the eyestalk ganglia of GFP DSI animals. Number of metabolites whose levels were significantly changed by CHH dsRNA was 149 and 181 in the muscle and 24 and 12 in the hepatopancreas, at 24 and 48 hpi, respectively. Principal component analysis of these metabolites show that metabolic effects of silencing CHH gene expression were more pronounced in the muscle (with the cluster of CHH DSI group clearly being separated from that of SAI group at 24 hpi) than in the hepatopancreas. Moreover, pathway analysis of the metabolites closely related to carbohydrate and energy metabolism indicate that, for CHH DSI animals at 24 hpi, metabolic profile of the muscle was characterized by reduced synthesis of NAD⁺ and adenine ribonucleotides, diminished levels of ATP, lower rate of utilization of carbohydrates through glycolysis, and a partially rescued TCA cycle, whereas that of the hepatopancreas by unaffected levels of ATP, lower rate of utilization of carbohydrates, and increased levels of ketone bodies. The combined results of metabolic changes in response to silenced CHH gene
no role in study design, data collection and
analysis, decision to publish, or preparation of the
manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

expression reveal that metabolic functions of CHH on the muscle and hepatopancreas are
more diverse than previously thought and are differential between the two tissues.

**Introduction**

Crustacean eyestalks contain a neuroendocrine tissue, the X-organ/sinus gland complex, which synthesizes and secretes several neurohormones that are implicated in regulating a wide variety of physiological functions, including metabolism, reproduction, growth, ionic balance, and color changes. The wide array of biological processes that are under the hormonal control of the eyestalk factors highlights the importance of the XO/SG complex in crustacean physiology [1-3]. Among these neuropeptides, crustacean hyperglycemic hormone (CHH) is the prototypic member of the CHH family, which also includes molt-inhibiting hormone (MIH), vitellogenesis-inhibiting hormone (VIH), and mandibular organ-inhibiting hormone (MOIH), and insect ion transport peptide (ITP) [4, 5]. CHH, a 72- or 73-residue polypeptide, is characterized by 3 disulfide bridges formed by 6 highly conserved cysteine residues and N- and C-terminally blocked, with the C-terminal amide being critical for its hyperglycemic activity [6-10]. Molecular characterizations of CHH precursor indicated that the precursor protein consists of a signal peptide, a CHH precursor-related peptide, and a mature CHH peptide [11].

Functionally, the best characterized roles for CHH is its effects on carbohydrate metabolism [12] and CHH is considered a stress hormone that elicits stress-induced hyperglycemia [13-17]. The majority of the data concerned the effects of CHH on glycogen metabolism. Thus, early *in vivo* and *in vitro* studies have shown that CHH preparations (*i.e.*, eyestalk or sinus gland extracts) rapidly mobilized glycogen stores in the hepatopancreas and muscle of *Orconectes limosus* and *Uca pugilator*, resulting in decreased levels of tissue glycogen or increased levels of glucose release from the tissues [18-22], and that activity of glycogen synthase was inhibited by CHH preparations (*i.e.*, sinus gland extracts) and increased by eyestalk ablation [20]. A more recent study showed that transcript levels of glycogen synthase were decreased and those of glycogen phosphorylase increased in hepatopancreatic tissues of *Marsupenaeus japonicus* treated *in vitro* with recombinant CHH [23]. Although most early studies employed tissue extracts (eyestalk ganglia or sinus glands) rather than purified CHH, the combined data are consistent with the notion that CHH mobilizes tissue glycogen store through regulating activity and amount of enzymes in glycogen metabolism. Specific CHH binding sites have been characterized in multiple tissues (including the hepatopancreas and muscle) [24-27].

Surprisingly, with the exception of regulation of glycogen metabolism as just mentioned, evidence for other metabolic effects of CHH was rather limited. Again, relevant data were those obtained from early studies. Thus, it has been proposed that CHH stimulates glycolysis, based on observations that hypoxia increased glucose and lactate levels in the hemolymph of intact animals (*Carcinus maenas* and *Orconectes limosus*), but not in eyestalk-ablated animals; the hypoxia-increased hemolymph levels of glucose and lactate were however observed in eyestalk-ablated animals receiving injection of CHH [28]. The authors considered that these data indicated that CHH stimulates glycolytic flux by increasing the availability of glycolysis substrate (glucose) through its stimulatory effect on glycogen mobilization, and hence leading to elevated levels of lactate [28, 29].

During the past 2 decades or so, much more research effort was directed towards molecular and biochemical characterization of CHH-family peptides and their encoding gene in both crustacean and non-crustacean species [2, 3, 4, 5, 11, 30, 31, 32], while the regulatory function
of its prototypic member is still poorly characterized. Thus, the aim of the present study was to elucidate in a more comprehensive manner the regulatory effects of CHH on the metabolism of its target tissues. A double-stranded RNA (dsRNA)-based gene silencing method was developed and validated for specific and effective silencing of CHH gene expression in the eyestalk ganglia, with concomitant decreases in CHH levels in the sinus gland and hemolymph. Subsequently, metabolites in the tissues (the hepatopancreas and muscle) harvested from saline-injected and CHH dsRNA-injected animals were identified and quantified using nuclear magnetic resonance spectroscopy. The resulting data clearly indicate that CHH has more diverse regulatory effects on tissue metabolism than previously known and exerts its metabolic effects in tissue-specific and complementary manners.

# Materials and methods

## Experimental animals and procedures

Experimental protocols involving animals employed by the present study were approved by the review committee of National Changhua University of Education (Permit number: NCUE-103320001), in full accordance with the recommendations (Guidelines for Management and Use of Experimental Animals) set by the Council of Agriculture, Taiwan. Animals (*P. clarkii*) used in the present study were trapped and supplied by local fisherman from Xihu River, Miaoli County, Taiwan (24.25° N, 120.45° E), and reared in the laboratory in water tanks with temperature held at 24 ± 1°C and a 12L/12D photoperiod regime [16]. Animals used in the present study were all in intermolt stage [33] with body weights ranging from 5 to 6 grams.

Experimental procedures described below were always carried out in the morning (9-11 a.m.) in order to reduce possible interference due to circadian changes in CHH and hemolymph glucose levels [34, 35]. Animals were separated into 3 groups (6 animals per group): the untreated group, saline-injected (SAI) group, and dsRNA-injected (DSI) group.

Animals in the SAI and DSI groups were intramuscularly injected using an insulin syringe coupled to a 29-gauge needle (Terumo Medical Corp.) with crayfish saline (50µl/animal, 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES)-buffered Van Harreveld saline, pH = 7.4, [16]) and saline containing CHH dsRNA or green fluorescence protein (GFP) dsRNA (50µl/animal, 30µg dsRNA/g body weight), respectively. Hemolymph, eyestalk ganglia, muscle, and hepatopancreas were collected from animals of the untreated group and from animals of the SAI and DSI groups at 24, 48, and 72 h post injection (hpi). Briefly, hemolymph was withdrawn using a syringe coupled to a 26G syringe and processed according to Yeh et al. [36] for quantification of CHH and glucose levels (see below "CHH sandwich enzyme-linked immunosorbent assay (ELISA) and glucose assay"). Sinus glands were completely removed from the eyestalk ganglia, homogenized in HEPES-buffered saline, and centrifuged at 16,000g for 30 min at 4°C; the supernatants were stored at -80°C until used for quantification of CHH levels. The rest of the eyestalk ganglia were used for extraction of total RNA (Purescript RNA Isolation Kit, Gentra), which was then treated with RQ1 RNase-free DNase (Promega) according to the supplier's procedures. For cDNA synthesis, a reverse transcription reaction was performed as previously described [37] and the synthesized cDNA was kept at -20°C until used for quantitative real-time polymerase chain reaction (qPCR) (see below "Semi-quantitative real-time polymerase chain reaction"). Dissected muscle and hepatopancreas were boiled using a hot plate (100°C, 5 min, Corning) and frozen in liquid nitrogen until used for nuclear magnetic resonance (NMR) analysis (see below "Nuclear magnetic resonance analysis of metabolites").
## Production of CHH and GFP double-strained RNAs (dsRNAs)

Double-stranded RNAs were produced using an *in vitro* transcription reaction driven by T7 promoter and T7 RNA Polymerase (T7 RiboMAX™ Express RNAi System, Promega).

For the production of CHH dsRNA, sense and antisense CHH DNA templates each with a T7 promoter sequence were separately amplified by PCR with primers containing T7 promoter sequence. Primer pairs used were: T7-CHH1-352-357-F and CHH1-472-492-R (for sense strand) or CHH1-352-357-F and T7-CHH1-472-492-R (for antisense strand) (S1 Table). A 25-µl PCR reaction contained 5µl 5X Colorless GoTaq Flexi Buffer (Promage), 2.5µl MgCl₂ (25 mM), 0.5µl dNTP (10mM), 0.5µl forward primer (10mM), 0.5µl reverse primer (10mM), 0.13µl GoTaq DNA polymerase (5U/µl) and 1µl cDNA (eyestalk cDNA synthesized as mentioned in "Experimental animals and procedures"). PCR products of expected size were purified using a commercial kit (PCR Advanced™, VIOGENE) and auto-sequenced to confirm their sequences. Using the CHH sense strand and antisense strand as DNA templates, complementary single-stranded RNA was produced *in vitro* using T7 RiboMAX™ Express RNAi System (Promega) according to the supplier's instruction. Resulting sense and antisense RNAs were mixed, heated at 70°C for 10min, and annealed to each other by cooling the mixture to 25°C at a rate of -15°C/hr. The dsRNA mixture was then treated by RQ1 RNase-Free DNase and RNase A Solution (Promega) to remove unwanted DNA and single-stranded RNA, precipitated with isopropanol, air-fried, and re-suspended in nuclease-free water following the supplier's instruction. The purified dsRNA samples were spectrophotometrically quantified (NanoDrop1000, Thermo), confirmed for integrity using non-denaturing 2% agarose gel electrophoresis (S1 Fig), and kept at -80°C before being used for *in vivo* injection.

Sense and antisense GFP DNA templates each with a T7 promoter sequence were separately amplified by PCR. Briefly, a plasmid containing a GFP insert (L4417, Addgene) was amplified using a primer pair: GFP-F-58-78-P144 and GFP-R-181-201-P144 (S1 Table). The resulting amplified GFP product was purified from the reaction mixture and then used as the template for a second PCR using the PCR primers containing T7 promoter sequence (S1 Table): F-T7-GFP-P140-ds & GFP-R-181-201-P144 (for GFP sense strand) or GFP-F-58-78-P144 & R-t7-gfp-p140-ds (for GFP antisense strand). Subsequently, using the GFP sense strand and antisense strand as DNA templates for *in vitro* transcription reaction, complementary single-stranded RNA was produced and annealed, and the resulting GFP dsRNA mixture was treated and processed as mentioned above for CHH dsRNA.

## CHH sandwich enzyme-linked immunosorbent assay (ELISA) and glucose assay

CHH peptide levels in sinus gland and hemolymph samples were quantified using an established sandwich ELISA with 2 anti-CHH antibodies; specificity of the assay for the crayfish CHH has been demonstrated and reported [38, 39]. Glucose levels in the hemolymph samples were quantified using an established glucose assay (GAGO20-1KT, Sigma). Protocols and reagents for CHH ELISA and glucose assay have been described previously [16, 38, 39].

## Semi-quantitative real-time polymerase chain reaction (qPCR)

A semi-quantitative real-time PCR was used to assess the effects of CHH or GFP dsRNA on the expression of the target gene (*CHH*) and a related gene *MIH* (molt-inhibiting hormone gene, a CHH family-peptide gene) that serves as an off-target gene. Primers used were CHHe2F and SG-R-400-427 for CHH gene [40] and MIH-F and MIH-R for MIH gene (S1 Table). 18S rRNA (18S ribosomal RNA gene) was used as the reference gene [40], and the
primers used for its amplification were 18S-F-1280 (forward) and 18S-R-1360 (reverse) (S1 Table). Briefly, eyestalk ganglia cDNA templates were amplified in a 10-µl reaction using a commercially available kit (LightCycler FastStart DNA Master SYBR Green I, Roche) according to manufacturer's protocol. Reactions were performed under the following conditions: an initial denaturation (10 min, 95°C), 40 cycles of denaturation (10 s, 95°C), annealing (7 s, 60°C), and extension (14 s, 72°C); the rate of temperature change was 20°C/s.

Amplification efficiency and specificity were validated according to previously described methods [37]. For each of the 3 amplicons (*CHH*, *MIH*, and *18S rRNA*), PCR conditions were optimized to ensure an amplification efficiency greater than 90%. To ensure amplification specificity, melting curve analysis was performed from 65°C to 94°C at the rate of 0.1°C/s after each PCR assay. For each of the 3 amplicons, melting curve analyses showed only a single peak indicating the amplification specificity of the reaction.

The comparative threshold cycle method [41] was used to determine the transcript levels. Threshold cycle number for each PCR was determined using the Second Derivative Maximum algorithm that identifies the first turning point of the fluorescence curve (Light-Cycler software, v. 3.5). Transcript levels of *CHH* or *MIH* were normalized to those of *18S rRNA* and the normalized levels at each time points after injection were expressed relative to the levels of untreated animals (arbitrarily designated as the calibrator), according to previously described equations [37].

¹H Nuclear magnetic resonance (NMR) analysis of metabolites

Tissue metabolites were extracted using an extraction procedure according to Teng et al. [42] with modifications. In brief, tissues (the muscle and hepatopancreas) collected from the animals in different treatment groups (see "Experimental animals and procedures") were frozen in liquid nitrogen and ground in ceramic mortar using a pestle. Ground tissue powders were added to a centrifuge vial containing cold 50% methanol (250 µl per 0.1g tissue powders), and homogenized at 4°C for 15 min using a tissue homogenizer (Bullet Blender®, Next Advance). Resulting homogenate was collected and an aliquot of the homogenate (250 µl) was transferred to a new vial containing an equal volume of 50% methanol, followed by adding 0.5ml of chloroform to the mixture. The mixture was further homogenized for 20 min at 4°C, centrifuged at 12,000g for 15 min at 4°C. The resulting aqueous phase was collected. Extraction solvents were completely removed using a vacuum concentrator (Savant SpeedVac, Thermo Scientific). Samples were then reconstituted in deuterium oxide ($D_2O$, Merck) containing 1 mM sodium salt of 3-(trimethylsilyl) propionate-2,2,3,3-d₄ (TSP) for NMR analysis.

All NMR spectra were acquired at 25°C on a Varian Inova 500 MHz NMR spectrometer (Varian, Inc.). The residual water resonance was pre-saturated using the PRESAT sequence. Typically, 128 transients were summed over a spectral width of 8 kHz into 32,000 data points per sample, and the summed free induction decays (FIDs) were Fourier transformed with exponential line-broadening of 0.3 Hz [43].

Statistical and principal component analyses of metabolites

The spectra were phased, baseline-corrected, and referenced to TSP (at 0.0 ppm) as well as quantified and qualified using a commercial software package (Chenomx NMR suite 4.6, Chenomx, Inc., Edmonton, Alberta, Canada) with the Chenomx 500-MHz (pH 6-8) library. The spectra data were pre-processed using normalization and scaling to remove possible bias arising due to sample handling and sample variability. Normalization (by sum) was performed in order to minimize possible differences in concentration between samples. In addition, the
spectral region from 5.0 to 4.7 ppm containing the residual water resonance was excluded, and
all subsequent analysis used centred scaling.

After data pre-processing, a univariate statistical analysis (Student's t-test) was used to ver-
ify the significance of difference in metabolite levels between the CHH dsRNA-injected group
and the saline-injected group. p<0.05 was considered to indicate a statistically significant dif-
ference. These significantly changed metabolites were further analyzed by heuristic methods of
dimension reduction (principal component analysis [PCA]) using commercially available soft-
ware (SIMCA-P 11.5; Umetrics, Umeå, Sweden), and all subsequent analysis used centred scal-
ing. The number of PCA components was calculated using an auto-fit model in SIMCA-P
(Umetrics).

Statistical analysis of transcript, peptide and glucose levels

Data for CHH transcript and peptide levels and glucose levels were analyzed by one-way analy-
sis of variance followed by Tukey's HSD test or by Student's t-test using computer software
(SPSS Manager, SPSS Inc.) and are reported as the mean ± standard deviation.

**Results**

Effects of dsRNA treatments on transcript levels in the eyestalk ganglia

Transcript levels in the eyestalk ganglia were quantified by a relative real-time PCR. Animals of the CHH dsRNA-injected (DSI) group exhibited a progressive decrease in CHH transcript levels over a period of 72 h post injection (hpi). Thus, compared to the those in untreated animals, CHH transcript levels in CHH DSI animals were significantly lower, 52.0 ± 15.8%, 18.5 ± 27.0%, and 4.0 ± 2.4% of the untreated value at 24, 48, and 72 hpi, respectively (Fig 1A). On the contrary, CHH transcript levels in the saline-injected (SAI) animals were relatively stable from 24 to 72 hpi, not significantly different from those in the untreated animals at any time point examined (82.4 ± 26.4%, 88.5 ± 4.7%, and 70.9 ± 21.1% of the untreated value at 24, 48, and 72 hpi, respectively), but were higher than those in the CHH DSI animals at each time point (with statistical significance being recorded at 48 and 72 hpi) (Fig 1A). These results clearly indicate that expression of CHH gene in the eyestalk ganglia was significantly knocked down by CHH dsRNA treatment.

Specificity of the CHH dsRNA treatment was tested by examining its effect on the tran-
script levels of a gene encoding another CHH family-peptide, molt-inhibiting hormone
(MIH). Compare to the untreated levels, MIH transcript levels in the CHH DSI animals were
not significantly different, 95.3 ± 32.8%, 117.8 ± 37.5%, and 86.3 ± 10.2% of the untreated
value at 24, 48, and 72 hpi, respectively (Fig 1B). Similarly, MIH transcript levels in the SAI
animals were not significantly different from the untreated levels (95.7 ± 9.7%, 77.7 ± 40.7%,
and 143.4 ± 67.0% of the untreated value at 24, 48, and 72 hpi, respectively), nor were they sig-
nificantly different from those in the CHH DSI group at any time point (Fig 1B). Specificity of
the silencing effect of CHH dsRNA treatment on CHH was additionally tested by another
dsRNA with a sequence unrelated to CHH-green fluorescence protein (GFP) dsRNA. CHH
transcript levels in the GFP DSI or SAI animals were not significantly from those in the
untreated animals (DSI: 72.0 ± 15.8%, 77.0 ± 23.0%, and 62.0 ± 17.0% of the untreated value;
SAI: 82.4 ± 26.4%, 88.5 ± 4.7%, and 70.9 ± 21.1% of the untreated value, at 24, 48, and 72 hpi,
respectively), nor were they significantly different from each other at any time point examined
(Fig 1C).
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Effects of dsRNA treatment on transcript levels in the eyestalk ganglia of the crayfish *Procambarus clarkii*. Animals were separated into 3 groups: the untreated, saline-injected, and dsRNA-injected groups. Using a relative real-time qPCR, CHH (A, C) and MIH (B) transcript levels were quantified for the samples of the eyestalk ganglia harvested from the untreated animals (blue bar), and the saline-injected (crimson bar) and dsRNA-injected animals (green bar) at 24, 48, and 72 h post injection (hpi) of CHH dsRNA (A, B), GFP dsRNA (C), or saline (A, B, C). Transcript levels are normalized to a reference gene (18S rRNA) and expressed relative to the untreated levels. **, and *** represent significantly different from the untreated or the saline-injected values at the level of 0.01 and 0.005, respectively (n = 5 for each data).</figcaption></figure>

doi:10.1371/journal.pone.0172557.g001

Effects of CHH dsRNA on the levels of CHH peptide in the sinus gland and hemolymph

CHH peptide levels in the sinus gland and hemolymph were measured by a sandwich ELISA. CHH DSI animals exhibited a progressive decrease in CHH peptide levels in the sinus gland over a period of 72 h, with a temporal profile similar to that of changes in transcript levels. Compared to the untreated levels (3.1 ± 0.7 x10³ pmole/mg protein), CHH levels in CHH DSI animals were lower at 24 hpi (2.3 ± 0.7 x10³ pmole/mg protein) and significantly lower at 48 and 72 hpi (1.5 ± 0.5 x10³ and 0.7 ± 0.4 x10³ pmole/mg protein, respectively) (Fig 2A). On the contrary, CHH peptide levels in the SAI animals remained relatively stable from 24 to 72 hpi (3.8 ± 1.5, 3.5 ± 1.1, 2.9 ± 0.6 x10³ pmole/mg protein, respectively), not significantly different from those in the untreated animals at any time point examined, but were significantly higher than those in the CHH DSI animals at 48 and 72 hpi (Fig 2A).

Compared to those of the untreated animals (39.5 ± 7.0 fmole/ml), CHH peptide levels in the hemolymph of the CHH DSI animals at 24 hpi (12.0 ± 1.6 fmole/ml) were significantly decreased, but were not significantly changed at 48 and 72 hpi (17.8 ± 4.4 and 34.7 ± 6.9 fmole/ml, respectively) (Fig 2B). CHH peptide levels in the SAI animals remained relatively
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Effects of CHH dsRNA treatment on CHH peptide levels in the crayfish *Procambarus clarkii*. Animals were separated into 3 groups: the untreated, saline-injected, and dsRNA-injected groups. Using a CHH-specific sandwich ELISA, CHH peptide levels were quantified for the samples of the sinus gland (A) and hemolymph (B) harvested from the untreated animals (blue bar), and the saline-injected animals (crimson bar) and CHH dsRNA-injected animals (green bar) at 24, 48, and 72 h post injection (hpi) of dsRNA or saline. *, **, and *** represent significantly different from the untreated or saline-injected values at the level of 0.05, 0.01 and 0.005, respectively (n = 5 for each data).</figcaption></figure>

doi:10.1371/journal.pone.0172557.g002

stable from 24 to 72 hpi (42.9 ± 9.1, 41.7 ± 7.4, 34.4 ± 10.6 fmole/ml, respectively), not significantly different from those in the untreated animals at any time point examined, but significantly higher than those in the CHH DSI animals at 24 and 48 hpi (Fig 2B).

## Effects of CHH dsRNA on the levels of hemolymph glucose

There was no significant change in hemolymph glucose levels for DSI or SAI animals from 24 to 72 hpi (DSI: 21.9 ± 9.2, 24.4 ± 9.4, and 17.9 ± 4.8 mg/dL, respectively; SAI: 24.5 ± 7.5, 24.9 ± 14.1, and 21.1 ± 1.9 mg/dL, respectively), nor were the glucose levels significantly different between the 2 groups of animals at each time point (Fig 3). Hemolymph glucose levels for DSI or SAI animals at any time point examined were not significantly different from the untreated levels (Fig 3).

## Effects of CHH dsRNA on the levels of tissue metabolites

Tissue metabolites detected, quantified, and found significantly different in metabolite levels between SAI and CHH DSI groups were listed in supplementary tables (S2 and S3 Tables), which includes 149 and 181 metabolites from the muscle at 24 and 48 hpi, respectively, and 24 and 12 metabolites from the hepatopancreas at 24 and 48 hpi, respectively. Levels of these metabolites in the untreated animals are also given for reference and they are not significantly different from those of the SAI group (S2 and S3 Tables). Principal component analysis (PCA) of the metabolites in the muscle data show a clear separation between SAI versus CHH dsRNA groups at 24 hpi (Fig 4A; combined contribution of the first two PCs explained 83.3% of the total variation), while the separation became less distinguishable at 48 hpi (Fig 4B; combined contribution of the first two PCs explained 89.8% of the total variation). PCA for the metabolites in the hepatopancreas, at both 24 and 48 hpi, show that the cluster of CHH DSI group was somewhat separable from that of SAI group but the separations were less clear-cut than those in the muscle (Fig 4C and 4D; combined contribution of the first two PCs explained 79.7% and 90.2% of the total variation for 24 and 48 hpi, respectively).

Levels of the metabolites that are closely related to carbohydrate and energy metabolism in the tissues of DSI and SAI animals are presented and statistically compared in Tables 1
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Effects of CHH dsRNA treatment on hemolymph glucose levels in the crayfish *Procambarus clarkii*. Animals were separated into 3 groups: the untreated, saline-injected, and dsRNA-injected groups. Using a glucose assay, glucose levels were quantified for the hemolymph samples harvested from the untreated animals (blue bar), and the saline-injected animals (crimson bar) and CHH dsRNA-injected animals (green bar) at 24, 48, and 72 h post injection (hpi) of dsRNA or saline (n = 5 for each data).</figcaption></figure>

doi:10.1371/journal.pone.0172557.g003

and 2 and described below. In the muscle, 24 h after CHH dsRNA treatment, when hemolymph CHH levels were significantly decreased compared to those of the untreated and the SAI groups, levels of the reduced form of nicotinamide adenine dinucleotide (NADH), as well as those of its oxidized form (NAD⁺), were significantly reduced, with ratio of NADH/NAD⁺ being significantly increased. Levels of 2 precursors to NAD⁺ (quinolinate, deamino-NAD⁺) were also significantly decreased (Table 1). In addition, levels of ATP and ADP were significantly decreased, with the ratio of ATP/ADP remained relatively constant; likewise, levels of AMP were significantly decreased (Table 1). Regarding glycolytic metabolites, levels of 2-phosphoglycerate, pyruvate, and lactate were significantly decreased (Table 1). Levels of 2 monosaccharides (fructose and galactose) and 2 disaccharides (sucrose and lactose) were also significantly decreased, while the decrease in levels of glucose was not significant. Tri-carboxylic acid (TCA) cycle was also affected, with levels of oxaloacetate, citrate, and isocitrate significantly increased, but those of malate significantly decreased. Moreover, levels of 2 ketone bodies (β-hydroxybutyrate and acetoacetate) were also significantly reduced.

Metabolic profile in the muscle of CHH DSI group at 48 hpi was similar to that at 24 hpi (Table 1). Levels of the majority of the metabolites were similarly significantly changed at 48 hpi (but in general with weaker levels of statistical significance), except those of several metabolites (pyruvate, lactose, and oxaloacetate) that became not statistically significantly changed (compared to those of SAI group) (Table 1).

Effects of CHH dsRNA treatment on the hepatopancreas were relatively mild and short-lasting (Table 2), compared to those on the muscle (Table 1). At 24 hpi, levels of ATP and
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

**Fig 4.** Principal component analysis plot of tissue samples of the crayfish *Procambarus clarkii* showing differential metabolic profiles between saline-injected and CHH dsRNA-injected groups. A, B: the muscle at 24 and 48 hpi, respectively; C, D: the hepatopancreas at 24 and 48 hpi, respectively. Triangle: CHH dsRNA injected group; Square: saline-injected group. Each symbol represents the metabolome of a single sample along two principal components.

doi:10.1371/journal.pone.0172557.g004

ADP, as well as the ratio of ATP/ADP, remained relatively constant; so did levels of NADH and NAD+, and the ratio of NADH/NAD+ (Table 2). Levels of fructose, galactose, sucrose, and lactose were slightly increased (although not significantly) (Table 2). Levels of 2-phosphoglycerate was increased (not significantly) and those of pyruvate, acetaldehyde, and acetate significantly decreased (Table 2). As for the TCA cycle intermediates, only levels of oxaloacetate were increased (although not significantly), while those of the rest of intermediates remained relatively constant (Table 2). Levels of the 2 ketone bodies were either remained relatively constant (acetoacetate) or slightly increased (β-hydroxybutyrate).

Metabolic profile in the hepatopancreas of CHH DSI animals at 48 hpi is essentially not different from that of SAI animals. None of the metabolites was significantly different in levels between CHH DSI and SAI animals (Table 2), including pyruvate, acetaldehyde, and acetate that were significantly decreased at 24 hpi.

## Discussion

To evaluate metabolic effects of crustacean hyperglycemic hormone on its targets, a double-stranded RNA-based gene silencing method was developed and validated. Resulting data indicate that CHH dsRNA injected *in vivo* effectively silenced CHH expression in the crayfish *P*.
<table><caption>Table 2. Effects of CHH dsRNA treatment on levels of metabolites related to energy and carbohydrate metabolism in the hepatopancreas of the crayfish *Procambarus clarkii*.</caption><thead><tr><th></th><th></th><th></th><th>2-Phosphoglycerate</th><th>Acetaldehyde</th><th>Acetate</th><th>Acetoacetate</th><th>ADP</th><th>AMP</th><th>ATP</th><th>ATP/ADP</th><th>Citrate</th></tr></thead><tbody><tr><td rowspan="6">24 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>1.23</td><td>0.00131</td><td>0.00498</td><td>0.0026</td><td>0.0024</td><td>0.0137</td><td>0.0022</td><td>0.95707</td><td>0.0213</td></tr><tr><td>Std</td><td>1.06405</td><td>0.0005</td><td>0.0007</td><td>0.00053</td><td>0.00211</td><td>0.01575</td><td>0.00195</td><td>0.06961</td><td>0.00681</td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>1.84987</td><td>0.00079</td><td>0.00412</td><td>0.00247</td><td>0.00214</td><td>0.01213</td><td>0.00198</td><td>0.95029</td><td>0.02479</td></tr><tr><td>Std</td><td>0.01851</td><td>0.00023</td><td>0.00059</td><td>0.00076</td><td>0.00122</td><td>0.00754</td><td>0.0011</td><td>0.06176</td><td>0.00697</td></tr><tr><td colspan="2">p value</td><td>0.28186</td><td>0.02157</td><td>0.02756</td><td>0.72378</td><td>0.78838</td><td>0.8154</td><td>0.80171</td><td>0.87844</td><td>0.56125</td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>150.4</td><td>60.4</td><td>82.7</td><td>94.9</td><td>89.2</td><td>88.5</td><td>89.9</td><td>99.3</td><td>116.4</td></tr><tr><td rowspan="6">48 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>1.83896</td><td>0.00102</td><td>0.00458</td><td>0.00331</td><td>0.00255</td><td>0.01344</td><td>0.00239</td><td>0.94588</td><td>0.02019</td></tr><tr><td>Std</td><td>0.04985</td><td>0.00054</td><td>0.00038</td><td>0.00078</td><td>0.0016</td><td>0.0076</td><td>0.00146</td><td>0.04946</td><td>0.00148</td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>1.86875</td><td>0.00099</td><td>0.00499</td><td>0.00474</td><td>0.00216</td><td>0.01174</td><td>0.00204</td><td>0.96204</td><td>0.01844</td></tr><tr><td>Std</td><td>0.01007</td><td>0.00057</td><td>0.00142</td><td>0.00384</td><td>0.00111</td><td>0.00732</td><td>0.00101</td><td>0.05158</td><td>0.00278</td></tr><tr><td colspan="2">p value</td><td>0.25824</td><td>0.9018</td><td>0.49009</td><td>0.37256</td><td>0.67888</td><td>0.72922</td><td>0.69084</td><td>0.6246</td><td>0.18055</td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>101.6</td><td>96.2</td><td>108.9</td><td>143.2</td><td>84.4</td><td>87.3</td><td>85.4</td><td>101.7</td><td>91.3</td></tr><tr><td rowspan="6">24 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>1.52608</td><td>0.00015</td><td>0.69862</td><td>1.83434</td><td>0.38443</td><td>0.00318</td><td>1.90983</td><td>0.01394</td><td>0.00905</td></tr><tr><td>Std</td><td>1.31517</td><td>0.00005</td><td>0.60071</td><td>1.57217</td><td>0.50044</td><td>0.00271</td><td>1.64139</td><td>0.00258</td><td>0.00578</td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>2.26947</td><td>0.00013</td><td>1.01698</td><td>2.77756</td><td>0.37728</td><td>0.0026</td><td>2.81394</td><td>0.01538</td><td>0.00719</td></tr><tr><td>Std</td><td>0.00634</td><td>0.00003</td><td>0.00627</td><td>0.03435</td><td>0.46156</td><td>0.00143</td><td>0.01891</td><td>0.00617</td><td>0.00301</td></tr><tr><td colspan="2">p value</td><td>0.29468</td><td>0.43559</td><td>0.20367</td><td>0.15667</td><td>0.98354</td><td>0.64652</td><td>0.30586</td><td>0.63267</td><td>0.50686</td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>148.7</td><td>87.8</td><td>145.6</td><td>151.4</td><td>98.1</td><td>81.8</td><td>147.3</td><td>110.4</td><td>79.4</td></tr><tr><td rowspan="6">48 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>1.51985</td><td>0.00016</td><td>0.86221</td><td>2.29038</td><td>0.38168</td><td>0.00376</td><td>2.17322</td><td>0.01445</td><td>0.00676</td></tr><tr><td>Std</td><td>1.28494</td><td>0.00009</td><td>0.40964</td><td>0.63201</td><td>0.46777</td><td>0.00193</td><td>1.10769</td><td>0.0018</td><td>0.00265</td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>2.28056</td><td>0.00016</td><td>1.02053</td><td>2.80172</td><td>0.36094</td><td>0.00306</td><td>2.51827</td><td>0.01367</td><td>0.00777</td></tr><tr><td>Std</td><td>0.00641</td><td>0.00007</td><td>0.00376</td><td>0.0251</td><td>0.45123</td><td>0.00101</td><td>0.77801</td><td>0.0023</td><td>0.00283</td></tr><tr><td colspan="2">p value</td><td>0.41299</td><td>0.91008</td><td>0.38726</td><td>0.2504</td><td>0.9525</td><td>0.49503</td><td>0.76795</td><td>0.52557</td><td>0.6201</td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>150.1</td><td>96.6</td><td>118.4</td><td>122.3</td><td>94.6</td><td>81.4</td><td>115.9</td><td>94.6</td><td>115</td></tr><tr><td rowspan="6">24 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>0.03327</td><td>3.9597</td><td>0.63987</td><td>0.00085</td><td>0.00448</td><td>2.53967</td><td>0.31075</td><td>0.01883</td><td></td></tr><tr><td>Std</td><td>0.01823</td><td>0.68821</td><td>0.54205</td><td>0.0002</td><td>0.00089</td><td>2.19471</td><td>0.52572</td><td>0.00446</td><td></td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>0.02699</td><td>3.93722</td><td>0.95191</td><td>0.00062</td><td>0.00436</td><td>3.87046</td><td>0.24418</td><td>0.02559</td><td></td></tr><tr><td>Std</td><td>0.00965</td><td>0.57783</td><td>0.00922</td><td>0.00017</td><td>0.00149</td><td>0.01908</td><td>0.43899</td><td>0.0091</td><td></td></tr><tr><td colspan="2">p value</td><td>0.48237</td><td>0.96332</td><td>0.1717</td><td>0.03826</td><td>0.86702</td><td>0.35288</td><td>0.83539</td><td>0.14069</td><td></td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>81.1</td><td>99.4</td><td>148.8</td><td>72.6</td><td>97.3</td><td>152.4</td><td>78.6</td><td>135.9</td><td></td></tr><tr><td rowspan="6">48 hpi</td><td rowspan="2">Saline-injected group (SAI)</td><td>Mean</td><td>0.02525</td><td>3.72475</td><td>0.94828</td><td>0.00073</td><td>0.00458</td><td>2.49687</td><td>0.23749</td><td>0.01872</td><td></td></tr><tr><td>Std</td><td>0.00995</td><td>0.61928</td><td>0.01216</td><td>0.00035</td><td>0.00085</td><td>2.13694</td><td>0.45763</td><td>0.00048</td><td></td></tr><tr><td rowspan="2">CHH dsRNA-injected group (DSI)</td><td>Mean</td><td>0.0275</td><td>3.67413</td><td>0.95547</td><td>0.00076</td><td>0.00483</td><td>3.92309</td><td>0.27263</td><td>0.01908</td><td></td></tr><tr><td>Std</td><td>0.00623</td><td>0.52807</td><td>0.00803</td><td>0.00041</td><td>0.00197</td><td>0.02581</td><td>0.45606</td><td>0.00434</td><td></td></tr><tr><td colspan="2">p value</td><td>0.71723</td><td>0.90518</td><td>0.2498</td><td>0.87741</td><td>0.76946</td><td>0.36712</td><td>0.90617</td><td>0.83544</td><td></td></tr><tr><td colspan="2">% SAI<sup>a</sup></td><td>108.9</td><td>98.6</td><td>100.8</td><td>104.6</td><td>105.4</td><td>157.1</td><td>114.8</td><td>101.9</td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0172557.t002

clarkii at the levels of CHH transcript expressed in the eyestalk ganglia, as well as those of CHH protein in the sinus gland and hemolymph, at least from 24 to 48 h after dsRNA treatment. Silencing specificity was demonstrated by a lack of effect of the CHH dsRNA on the expression of a gene encoding another CHH-family peptide-MIH, and additionally confirmed by a treatment of GFP dsRNA, which did not affect CHH gene expression.

Metabolic data reported in the present study clearly show that silencing of CHH gene expression significantly altered metabolic profiles of 2 CHH target tissues—the muscle and hepatopancreas, with the effects being metabolically manifested in tissue-specific but complementary manners (see below and Fig 5). PCA score plots indicate that the impacts were more

<footnote>a % SAI = (CHH DSI value/SAI value) x 100.</footnote>

<footnote>n = 6 for each data.</footnote>
<figure><img src="image_12.png" /><figcaption>**Fig 5.** Metabolic profiles of the muscle and hepatopancreas of the crayfish *Procambarus clarkii* 24 h after CHH dsRNA treatment. Red rectangle: significantly decreased from saline-injected levels; blue rectangle: significantly increased from saline-injected levels; dark gray rectangle: not significantly changed; light gray rectangle: not detected. Data obtained at 24 hpi (as listed in Tables 1 and 2) were used for constructing the metabolic networks relating to carbohydrate and energy metabolism. Dotted lines indicate multiple metabolic steps are involved that are not individually specified.</figcaption></figure>

doi:10.1371/journal.pone.0172557.g005

pronounced, especially at 24 h after dsRNA treatment, for the muscle than for the hepatopan-
creas. Importantly, the observed metabolic changes suggest a much wider array of effects of
CHH than previously realized [12].

Both tissues were more significantly impacted by CHH dsRNA treatment at 24 hpi, with the metabolic status at 48 hpi being reversing towards the saline-treated levels. Hence metabolic profiles for each tissue are discussed below based on changes at 24 h post CHH dsRNA treatment (see Fig 5). Central to these metabolic changes in the muscle are decreased levels of NAD⁺ and adenine ribonucleotides (AMP and ADP) (Table 1). It is well-established that glycolytic flux is dependent on a continuous supply of NAD⁺ [44]. Hence the observed inhibition of glycolysis (as evidenced by significantly lower levels of 2-phosphoglycerate, pyruvate, and lactate) (Fig 5) is at least in part due to low NAD⁺ availability. Consequently, glucose, instead of entering glycolysis, was presumably diverted for glycogen synthesis either intramuscularly or, perhaps more likely, in the hepatopancreas, after it was transported to this tissue (Fig 5). The suggested transport of glucose from the muscle to the hepatopancreas would therefore explain why hemolymph glucose levels in the CHH DSI animals, compared to those in the SAI animals, were not decreased in response to significantly lower levels of hemolymph CHH. Regarding other carbohydrates, significantly lower levels of sucrose, lactose, fructose, and galactose would most likely reflect contribution of these sugars to glycogen synthesis via glucose (Fig 5). Moreover, TCA cycle, where NAD⁺ is reduced to NADH in 3 of the reactions,
was apparently inhibited (Fig 5). The inhibition was likely also due to low availability of NAD+, which inhibits the three oxidation steps for which NAD+ is the cofactor (Lehninger et al., 1993); the reactions of TCA cycle that carry isocitrate through to oxaloacetate were inhibited, with levels of all the detected intermediates being decreased (Fig 5). Thus, diminished levels of NAD+ (and high NADH/NAD+ ratio), observed when CHH gene expression was silenced, is responsible at least partly for reducing TCA cycle flux. The flux was, however, partially rescued by (1) the 2 ketone bodies, β-hydroxybutyrate and acetoacetate (presumably derived from the hepatopancreas; see below) through their conversion to acetyl CoA, which in turn enters the cycle at the point of citrate formation and (2) by an anaplerotic reaction that feeds into the TCA cycle at the point of succinyl-CoA (from α-ketobutyrate that was derived from serine and methionine) (Fig 5). Finally, another prominent feature observed when CHH gene was silenced is severely diminished levels of ATP, which is closely related to the metabolic changes mentioned above. Thus, significantly lower levels of NADH, due mainly to a negatively impacted TCA cycle flux, is one of the factors that limited production of ATP through oxidative phosphorylation. Another factor that may lead to significantly lower levels of ATP would be low availability of its mono- and di-phosphate precursors-AMP and ADP (Table 1).

Thus, metabolic profile of the muscle under the condition of silenced CHH expression was characterized by significantly reduced synthesis of NAD⁺ and adenine ribonucleotides, lower rate of utilization of carbohydrates (glucose and other sugars) through glycolysis, and a negatively impacted (but partially rescued) TCA cycle, and greatly diminished supply of ATP.

Several earlier studies had suggested that CHH stimulates glycolytic flux [28, 29], although such a proposition has never been directly tested. Our observation that glycolytic flux was significantly negatively impacted by silencing expression of CHH gene provides the first piece of evidence indicating that CHH indeed exerts stimulatory effects on glycolysis. As discussed above, metabolic effects of CHH in the muscle are not limited to glycolysis, and regulation of synthesis of NAD⁺ and adenine ribonucleotides is central to the metabolic roles played by CHH. In this regard, it is interesting to note that, 5-φsophoribosyl 1-φrophosphae (PRPP) synthase is a tightly regulated enzyme whose product (PRPP) is a precursor for the synthesis of not only ribonucleotides but also for the synthesis of NAD⁺ [45, 46]. Thus, it is entirely possible that CHH could have diverse metabolic effects by regulating both biosynthetic pathways (NAD⁺ and adenine ribonucleotides) via affecting the activity of certain critical enzymes such as PRPP synthase.

In light of the metabolic effects of CHH suggested by the present study, it is relevant to discuss these data with findings reported by earlier studies [40, 47]. First, it has been shown that white spot syndrome virus (WSSV) activated CHH gene expression and enhanced release of CHH from the sinus glands in the crayfish *Procambarus clarkii*, resulting in an acute onset of significant increase, commencing as early as 3 h post infection, in hemolymph CHH levels that lasted for at least 2 days [40]. In other line of studies, an invertebrate Warburg effect, characterized by up-regulation of several metabolic pathways, including glycolysis, the pentose phosphate pathway, ribonucleotide biosynthesis, glutaminolysis and amino acid biosynthesis, was observed in the hosts (*Litopenaeus vannamei*) at the viral genome replication stage (12 h after being infected by WSSV) [47]. Taken together, these results and those presented in this study strongly indicate that CHH released by WSSV infection is responsible for inducing at least certain aspects of the WSSV-induced Warburg effect, including enhanced glycolysis, increased adenine ribonucleotide synthesis.

Compared to the muscle, the hepatopancreas exhibited less drastic changes in response to silencing of CHH gene expression. Both ratios of NADH/NAD+ and ATP/ADP were not significantly affected, nor were the levels of each metabolite. Accumulation of fructose, galactose, sucrose, and lactose probably reflected an increased uptake by the hepatopancreas of these
sugars, which were likely to be converted into glycogen via glucose (Fig 5). Glycolysis in the hepatopancreas was apparently also inhibited, as levels of pyruvate and lactate were decreased (Fig 5). Given the observation that levels of pyruvate and lactate were lower, increases (although not significant) in levels of 2-phosphoglycerate (a common intermediate for both glycolysis and gluconeogenesis) and glucose (Fig 5) suggest that rate of gluconeogenesis was elevated. Moreover, levels of β-hydroxybutyrate (a ketone body) were increased (Fig 5), an indication of excess acetyl CoA exceeding the capacity of TCA cycle, and presumably reflected enhanced β-oxidation of fatty acids under low glucose oxidation through glycolysis [48]. In addition, the 2 ketone bodies (β-hydroxybutyrate and acetoacetate) were most likely to be transported to the muscle, where they were converted into acetyl CoA and subsequently enter TCA cycle via formation of citrate (Fig 5).

Thus, metabolic status of the hepatopancreas under the condition of silenced CHH expression, in many ways different from that of the muscle, was characterized by unaffected supply of ATP, lower rate of utilization of carbohydrates (inhibited glycolysis and probably stimulated glycogen synthesis), and increased synthesis of ketone bodies (that were likely to be transported to the muscle) presumably via enhanced β-oxidation of fatty acids (Fig 5).

The metabolic responses of the muscle and hepatopancreas to silencing of CHH expression were obviously different but complementary. It has been showed that eyestalk ablation significantly decreased transcript levels of glycogen phosphorylase (GP) and increased those of glycogen synthase (GS) in the hepatopancreas, but not in muscle, in the kuruma prawn *Marsupenaeus japonicas* [23]. Conversely, transcript levels of GP were significantly increased and those of GS decreased, after hepatopancreatic tissues were treated *ex vivo* with recombinant CHH [23]. Thus, these data not only support our suggestion that with silenced CHH expression glycogen synthesis would be increased in the hepatopancreas, but also indicate that muscle and hepatopancreas were differentially regulated in terms of glycogen metabolism by CHH, highlighting the point we made that CHH regulates its target tissues in a tissue-specific manner. On the other hand, the differential responses after CHH gene silencing were complementary in that glucose was suggested to be transported from the muscle to the hepatopancreas for glycogen synthesis and that acetyl CoA synthesized via β-oxidation of fatty acids in the hepatopancreas was used for production of ketone bodies to be transported to the muscle (Fig 5).

In summary, a double-stranded RNA-based method was developed and validated for effective and specific silencing of CHH gene. When applied *in vivo* to the experimental animals (*Procambarus clarkii*), it allowed us to clearly reveal the differential and complementary effects of silencing CHH gene expression on two CHH target tissues. These data imply that CHH has a range of effects on these tissues much more diverse than previously known. Finally, but not least importantly, combined data of this and other previous studies [40, 47] suggest that the WSSC-released CHH is at least partly responsible for inducing the invertebrate Warburg effect that was manifested in infected hosts and thought to be beneficial for viral replication.

**Supporting information**

**S1 Fig. Confirmation of integrity of the double-stranded RNA.** CHH (A) or GFP (B) dsRNA produced using *in vitro* transcription reactions were separated by 2% agarose electrophoresis. M: RNA markers. Positions of 100-bp and 200-bp makers are labeled.
(TIF)

**S1 Table. Sequences of primers used in the present study.**
(PDF)
S2 Table. List of metabolites in the muscle of the crayfish *Procambarus clarkii* whose levels were significantly different between the saline-injected and CHH dsRNA-injected animals at 24 and 48h post injection (hpi). Data are expressed as mean ± standard deviation (n = 6 for each data). <sup>a</sup> Levels of the metabolites in the untreated animals are also included for reference.

(PDF)

S3 Table. List of metabolites in the hepatopancreas of the crayfish *Procambarus clarkii* whose levels were significantly different between the saline-injected and CHH dsRNA-injected animals at 24 and 48h post injection (hpi). Data are expressed as mean ± standard deviation (n = 6 for each data). <sup>a</sup> Levels of the metabolites in the untreated animals are also included for reference.

(PDF)

**Acknowledgments**

The authors thank Nuclear Magnetic Resonance Core Facility at Department of Chemistry, National Sun Yat-sen University for providing metabolomic analysis and related technical support.

**Author Contributions**

**Conceptualization:** C-YL.

**Data curation:** WL K-HC Y-CT S-FT C-HL C-YL.

**Formal analysis:** WL Y-CT S-FT L-JS C-HL C-YL.

**Funding acquisition:** C-YL WL.

**Investigation:** Y-CT S-FT L-JS.

**Methodology:** C-YL K-HC.

**Project administration:** C-YL.

**Resources:** C-YL K-HC.

**Supervision:** C-YL.

**Validation:** WL K-HC C-YL.

**Visualization:** WL C-YL.

**Writing – original draft:** C-YL WL.

**Writing – review & editing:** J-YT C-YL WL.

**References**

1. Cooke IM, Sullivan RE. Hormones and neurosecretion. In: Atwood HL, Sandeman DC, editors. New York: Academic Press. The Biology of Crustacea. 1982; 3:206–290.

2. Keller R. Crustacean neuropeptides: structures, functions, and comparative aspects. Experientia. 1992; 48:439–448. PMID: 1601108

3. Soyez D. Occurrence and diversity of neuropeptides from the crustacean hyperglycemic hormone family in arthropods. Ann N Y Acad Sci. 1997; 814:319–323. PMID: 9160986

4. Montagné N, Desdevises Y, Soyez D, Toullec J-Y. Molecular evolution of the crustacean hyperglycemic hormone family in ecdysozoans. BMC Evo Biol. 2010; 10:62.
5. Webster SG, Keller R, Dircksen H. The CHH-superfamily of multifunctional peptide hormones controlling crustacean metabolism, osmoregulation, moulting, and reproduction. *Gen Comp Endocrinol.* 2012; 175:217-233. doi: 10.1016/j.ygcen.2011.11.035 PMID: 22146796

6. Marco HG, Brandt W, Stoeva S, Voelter W, Gäde G. Primary structures of a second hyperglycemic peptide and of two truncated forms in the spiny lobster, *Jasus lalandii*. Peptides. 2000; 21:19-27. PMID: 10704715

7. Katayama H, Ohira T, Aida K, Nagasawa H. Significance of a carboxyl-terminal amide moiety in the folding and biological activity of crustacean hyperglycemic hormone. Peptides. 2002; 23:1537-1546. PMID: 12217413

8. Mosco A, Edomi P, Guarnaccia C, Lorenzon S, Pongor S, Ferrero EA, et al. Functional aspects of CHH C-terminal amidation in crayfish species. *Regul Pept.* 2008; 147:88-95. doi: 10.1016/j.regpep.2008.01.005 PMID: 18281112

9. Chang C-C, Tsai T-W, Hsiao N-W, Chang C-Y, Lin C-L, Watson RD, et al. Structural and functional comparisons and production of recombinant crustacean hyperglycemic hormone (CHH) and CHH-like peptides from the mud crab *Scylla olivacea*. *Gen Comp Endocrinol.* 2010; 167:68-76. doi: 10.1016/j.ygcen.2010.02.013 PMID: 20171218

10. Liu C-J, Huang S-S, Toullec J-Y, Chang C-Y, Chen Y-R, Huang W-S, et al. Functional assessment of residues in the amino- and carboxyl-termini of crustacean hyperglycemic hormone (CHH) in the mud crab *Scylla olivacea* using point-mutated peptides. *PLoS ONE*. 2015; 10(8):e0134983. doi: 10.1371/journal.pone.0134983 PMID: 26261986

11. De Kleijn DP, Van Herp F. Molecular biology of neurohormone precursors in the eyestalk of Crustacea. *Comp Biochem Physiol.* 1995; 112:573-579.

12. Santos EA, Keller R. Crustacean hyperglycemic hormone (CHH) and the regulation of carbohydrate metabolism: current perspectives. *Comp Biochem Physiol A.* 1993; 106:405-411.

13. Webster SG. Measurement of crustacean hyperglycemic hormone levels in the edible crab *Cancer pagurus* during emersion stress. *J Exp Bio.* 1996; 199:1579-1585.

14. Chang ES, Keller R, Chang SA. Quantification of crustacean hyperglycemic hormone by ELISA in hemolymph of the lobster, *Homarus americanus*, following various stresses. *Gen Comp Endocrinol.* 1998; 111:359-366. doi: 10.1006/gcen.1998.7120 PMID: 9707481

15. Chang ES, Chang SA, Beltz BS, Kravitz EA. Crustacean hyperglycemic hormone in the lobster nervous system: localization and release from cells in the subesophageal ganglion and thoracic second roots. *J Comp Neurol.* 1999; 414:50-56. PMID: 10494077

16. Zou H-S, Juan C-C, Chen S-C, Wang H-Y, Lee C-Y. Dopaminergic regulation of crustacean hyperglycemic hormone and glucose levels in the hemolymph of the crayfish *Procambarus clarkii*. *J Exp Zool.* 2003; 298:44-52.

17. Lorenzon S, Edomi P, Giulianini PG, Mettulio R, Ferrero EA. Variation of crustacean hyperglycemic hormone (CHH) level in the eyestalk and haemolymph of the shrimp *Palaemon elegans* following stress. *J Exp Biol.* 2004; 207:4205-4213. doi: 10.1242/jeb.01264 PMID: 15531641

18. Keller R, Andrew EM. The site of action of the crustacean hyperglycemic hormone. *Gen Comp Endocrinol.* 1973; 20:572-578. PMID: 4715238

19. Parvathy D. Endocrine regulation of carbohydrate metabolism during the moult cycle in crustaceans I. Effect of eyestalk removal in *Ocypoda platytarsis*. *Mar Biol.* 1972; 14:58.

20. Sedlmeier D, Keller R. The mode of action of the crustacean neurosecretory hyperglycemic hormone. I. Involvement of cyclic nucleotides. *Gen Comp Endocrinol.* 1981; 45 (1):82-90. PMID: 6269955

21. Sedlmeier D. The role of hepatopancreatic glycogen in the action of the crustacean hyperglycemic hormone (CHH). *Comp Biochem Physiol A.* 1987; 87:423-425.

22. Sedlmeier D. Studies on glycogen synthesis and breakdown in isolated hepatopancreas of *Orconectes limoae* (Crustacea, Decapoda). *Biochem Soc Trans.* 1981; 9:240.

23. Nagai C, Nagata S, Nagasawa H. Effects of crustacean hyperglycemic hormone (CHH) on the tran-script expression of carbohydrate metabolism-related enzyme genes in the kuruma prawn *Marsupe-naeus japonicus*. *Gen Comp Endocrinol.* 2011; 172:293-304. doi: 10.1016/j.ygcen.2011.03.019 PMID: 21447337

24. Kummer G, Keller R. High-affinity binding of crustacean hyperglycemic hormone (CHH) to hepatopan-creatic plasma membranes of the crab *Carcinus maenas* and the crayfish *Orconectes limosus*. Pept-tides. 1993; 14(1): 103-108. PMID: 8382808

25. Webster SG. High-affinity binding of putative molt-inhibiting horome (MIH) and crustacean hyperglyce-mic hormone (CHH) to membrane-bound receptors on the Y-organ of the shore crab *Carcinus maenus*. *Proc R Soc Lond Ser B Biol Sci.* 1993; 251:53-59.
26. Chung JS, Webster SG. Binding sites of crustacean hyperglycemic hormone and its second messengers on gills and hindgut of the green shore crab, *Carcinus maenas*: a possible osmoregulatory role. *Gen Comp Endocrinol.* 2006; 147:206–213. doi: 10.1016/j.ygcen.2006.01.002 PMID: 16504184

27. Katayama H, Chung JS. The specific binding sites of eyestalk- and pericardial organ-crustacean hyperglycaemic hormones (CHHs) in multiple tissues of the blue crab, *Callinectes sapidus*. *J Exp Biol.* 2009; 212:542–549. doi: 10.1242/jeb.022889 PMID: 19181902

28. Santos EA, Keller R. Effect of exposure to atmospheric air on blood glucose and lactate concentrations in two crustacean species: A role of the crustacean hyperglycemic hormone (CHH). *Comp Biochem Physiol A.* 1993; 106:343–347.

29. Santos EA, Keller R. Regulation of circulating levels of the crustacean hyperglycemic hormone: evidence for a dual feedback control system. *J Comp Physiol B.* 1993; 163:374–379.

30. Christie AE. Neuropeptide discovery in Ixodoidea: An *in silico* investigation using publicly accessible expressed sequence tags. *Gen Comp Endocrinol.* 2008; 157:174–185. doi: 10.1016/j.ygcen.2008.03.027 PMID: 18495123

31. Christie AE, Nolan DH, Garcia ZA, McCoole MD, Harmon SM, Congdon-Jones B, et al. Bioinformatic prediction of arthropod/nematode-like peptides in non-arthropod, non-nematode members of the Ecdysozoa. *Gen Comp Endocrinol.* 2011; 170:480–486. doi: 10.1016/j.ygcen.2010.11.002 PMID: 21074533

32. Lee C-Y, Tsai K-W, Tsai W-S, Huang S-Y, Chen Y-J. Crustacean hyperglycemic hormone: structural variants, physiological function, and cellular mode of action. *J Mar Sci Tech.* 2014; 22:75–81.

33. Nakatsuji T, Keino H, Tamura K, Yoshimura S, Kawakami T, Aimoto S, et al. Changes in the amounts of the molt-inhibiting hormone in sinus glands during the molt cycle of the american crayfish, *Procambarus clarkii*. *Zool Sci.* 2000; 17:1129–1136. doi: 10.2108/zzj.17.1129 PMID: 18522469

34. Gorgels-Kallen JL, Voorter CEM. The secretory dynamics of the CHH-producing cell group in the eye-stalk of the crayfish, *Astacus leptodactylus*, in the course of the day/night cycle. *Cell Tissue Res.* 1985; 241:361.

35. Kallen JL, Abrahamse SL, Van Herp F. Circadian rhythmicity of the crustacean hyperglycemic hormone (CHH) in the hemolymph of the crayfish. *Biol Bull.* 1990; 179:351–357.

36. Yeh F-C, Wu S-H, Lai C-Y, Lee C-Y. Demonstration of nitric oxide synthase activity in crustacean hemocytes and anti-microbial activity of hemocyte-derived nitric oxide. *Comp Biochem Physiol B.* 2006; 144:11–17. doi: 10.1016/j.cbpb.2006.01.007 PMID: 16516518

37. Chen H-Y, Watson RD, Chen J-C, Lee C-Y. Molecular characterization and gene expression pattern of two putative molt-inhibiting hormones from *Litopenaeus vannamei*. *Gen Comp Endocrinol.* 2007; 151:72–81. doi: 10.1016/j.ygcen.2006.11.016 PMID: 17222842

38. Wu H-J, Tsai W-S, Huang S-Y, Chen Y-J, Chen Y-H, Hsieh Y-R, et al. Identification of crustacean hyperglycemic hormone (CHH) and CHH-like (CHH-L) peptides in the crayfish *Procambarus clarkii* and localization of functionally important regions of CHH. *Zool Stud.* 2012; 51:288–297.

39. Wu S-H, Chen Y-J, Huang S-Y, Tsai W-S, Wu H-J, Hsu T-T, et al. Demonstration of expression of a neuropeptide-encoding gene in crustacean hemocytes. *Comp Biochem Physiol A.* 2012; 161A:463–468.

40. Lin L-J, Chen Y-J, Chang Y-S, Lee C-Y. Neuroendocrine responses of a crustacean host to viral infection: effects of infection of white spot syndrome virus on the expression and release of crustacean hyperglycemic hormone in the crayfish *Procambarus clarkii*. *Comp Biochem Physiol A.* 2013; 164:327–332.

41. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2^ΔΔCt^ Method. *Methods.* 2001; 25:402–408. doi: 10.1006/meth.2001.1262 PMID: 11846609

42. Teng Q, Ekman DR, Huang W, Collette TW. Impacts of 17-ethynylestradiol exposure on metabolite profiles of zebrafish (*Danio rerio*) liver cells. *Aquat Toxicol.* 2013; 130–131:184–191. doi: 10.1016/j.aquatox.2013.01.011 PMID: 23416411

43. Chiu K-H, Ding S, Chen Y-W, Lee C-H, Mok H-K. A NMR-based metabolomic approach for differentiation of hagfish dental and somatic skeletal muscles. *Fish Physiol Biochem.* 2011; 37:701–707. doi: 10.1007/s10695-011-9470-7 PMID: 21327496

44. Canto C, Menzies KJ, Auwerx J. NAD(+) metabolism and the control of energy homeostasis: a balancing act between mitochondria and the nucleus. *Cell Metab.* 2015; 22:31–53. doi: 10.1016/j.cmet.2015.05.023 PMID: 26118927

45. Lehninger AL, Nelson DL, Cox M. Principles of Biochemistry, Worth Publishers Inc., U.S.A. 2nd edition. 1993.

46. Chiarugi A, Dölle C, Felici R, Ziegler M. The NAD metabolome—a key determinant of cancer cell biology. *Nat Rev Cancer.* 2012; 12:741–752. doi: 10.1038/nrc3340 PMID: 23018234
47. Su M-A, Huang Y-T, Chen I-T, Lee D-Y, Hsieh Y-C, Li C-Y, et al. An invertebrate Warburg effect: a shrimp virus achieves successful replication by altering the host metabolome via the PI3K-Akt-mTOR pathway. PLoS Pathog. 2014; 10(6): e1004196. doi: 10.1371/journal.ppat.1004196 PMID: 24945378

48. Moran LA, Scrimgeour KG, Horton HR, Ochs RS, Rawn JD. Biochemistry. 1994. 2nd edition, Englewood Cliffs, New Jersey: Neil Patterson Publishers/Prentice-Hall, Inc.