Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The transcriptome sequencing and functional analysis of eyestalk ganglions in Chinese mitten crab (Eriocheir sinensis) treated with different photoperiods

  • Yang-yang Pang ,

    Contributed equally to this work with: Yang-yang Pang, Cong Zhang

    Roles Writing – original draft

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

  • Cong Zhang ,

    Contributed equally to this work with: Yang-yang Pang, Cong Zhang

    Roles Conceptualization

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

  • Min-jie Xu,

    Roles Conceptualization

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

  • Gen-yong Huang,

    Roles Conceptualization

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

  • Yong-xu Cheng ,

    Roles Conceptualization

    zny119@163.com (XZY); Chengyongxucrab@163.com (YXC)

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

  • Xiao-zhen Yang

    Roles Writing – review & editing

    zny119@163.com (XZY); Chengyongxucrab@163.com (YXC)

    Affiliations National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China, Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture, Shanghai Ocean University, Shanghai, China, Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China

Abstract

Photoperiod plays an important role in individual growth, development, and metabolism in crustaceans. The growth and reproduction of crabs are closely related to the photoperiod. However, as of yet, there are still no transcriptomic reports of eyestalk ganglions treated under different photoperiods in the Chinese mitten crab (Eriocheir sinensis), which is a benthonic crab with high commercial value in Asia. In this study, we collected the eyestalk ganglions of crabs that were reared under different photoperiods, including a control group (L: D = 12 h: 12 h, named CC), a constant light group (L: D = 24 h: 0 h, named LL) and a constant darkness group (L: D = 0 h: 24 h, named DD). RNA sequencing was performed on these tissues in order to examine the effects of different photoperiods. The total numbers of clean reads from the CC, LL and DD groups were 48,772,584 bp, 53,943,281 bp and 53,815,178 bp, respectively. After de novo assembly, 161,380 unigenes were obtained and were matched with different databases. The DEGs were significantly enriched in phototransduction and energy metabolism pathways. Results from RT-qPCR showed that TRP channel protein (TRP) in the phototransduction pathway had a significantly higher level of expression in LL and DD groups than in the CC group. We found that the downregulation of the pyruvate dehydrogenase complex (PDC) gene and the upregulation phosphoenolpyruvate carboxykinase (PPC) gene were involved in energy metabolism processes in LL or DD. In addition, we also found that the upregulation of the expression level of the genes Gαq, pyruvate kinase (PK), NADH peroxidase (NADH) and ATPase is involved in phototransduction and energy metabolism. These results may shed some light on the molecular mechanism underlying the effect of photoperiod in physiological activity of E. sinensis.

Introduction

Photoperiod is important to individual growth, development, reproduction and energy metabolism [13]. Researchers have found that the survival was high when combined with longer photoperiod, and weight gain and specific growth rate were higher at shorter photoperiods of crayfish (Astacus leptodactylus) [4]. There were also low hatching rates and population numbers when calanoid copepods (Acartia grani) were kept under constant light conditions [5]. The photoperiod also affects molting and influences the sexual maturation of crustaceans [69]. For example, a long photoperiod can accelerate the rate of metamorphosis in the lobster (Panulirus japonicus)[10]. The effects of photoperiod on energy metabolism have also been reported in crustaceans. For example, the digestive enzymatic activity of the general proteases trypsin and chymotrypsin were affected by the photoperiod in the prawn (Macrobrachium tenellum) [7]. Similarly, reactive oxygen species and lipid peroxidation are regulated by the glutathione system, which is influenced by photoperiods in crayfish (Procambarus clarkii)[11]. In addition, changes in energy metabolism due to the photoperiod have been reported in adult krill (Euphausia superba) and Arctic copepod (Calanus glacialis) [12, 13]. However, the study of the effects of photoperiod on E. sinensis, which is an economically important aquaculture freshwater species in China, is still sparse.

Eyestalks are an important phototransduction organ that can receive light signals through photoreceptors [1416]. Eyestalk ganglions, located in eyestalk, play an important role in phototransduction, energy metabolism, and endocrine regulation. A previous study showed that damaged eyestalk ganglion can lead to a destruction of the phototransduction pathway [17]. Moreover, some studies have found that eyestalks are involved in hatchability of embryos and gonadal development, such as in the crab (Dyspanopeus sayi) [18] and the mud crab (Scylla paramamosain)[19]. It was previously known that levels of the crustacean hyperglycemic hormone (CHH) in eyestalks are important for energy metabolism via the regulation glucose levels [20, 21], and these levels were controlled by photoperiods in the crayfish (P. clarkii) [22]. In our previous studies, we found that the eyestalks of E. sinensis express the DA2 receptor which participates in light adaptation during the dark hours [23]. However, there still lack information about the influence of photoperiods on the gene expression of eyestalk ganglion in E. sinensis.

The TRP channel protein (TRP) is a transient receptor potential cation channel with very diverse permeation and gating properties and participates in sensory and motile regulation process [24]. It is a vital gene in the phototransduction pathway [25]. Previous studies have shown that photoperiods can affect the TRP expression level in mice, TRP is involved in clock gene oscillations and energy balance control [26]. TRP can also change the osmolarity and fluid flow of individual cells in different environments [27]. For example, TRP can induce the increased levels of Ca2+, which further reduces the photoreceptor sensitivity of crayfish [28]. Additionally, both the pyruvate dehydrogenase complex (PDC) gene and the phosphoenolpyruvate carboxykinase (PPC) gene are related to energy metabolism [29, 30], and their activity can be influenced by photoperiods in mammals [31, 32].

The current study, the first to our knowledge, analyzes eyestalks ganglion transcriptomes of E. sinensis in all bright (LL), all dark (DD) and normal light (CC) conditions using Illumina HiSeqTM2000 technology. We attempt to identify important pathways and genes involved in photoperiod regulation. This study also provides a basis for further research.

Materials and methods

Experimental animals and sampling

Chinese mitten crabs (E. sinensis) were obtained from the Shuxin crab base in Chongming Island (121°30′~121°40′ E, 31°34′~31°37′ N), Shanghai (China) with body weights of 17.18 ± 2.2 g. They were stocked in clear glass aquaria (Length: Width: Height = 130: 60: 40 cm, Water depth = 25 cm) for 1 week, with a circulating system containing thoroughly aerated freshwater and UV-treated PVC tubes as shelter. The crabs were fed a basal diet once daily at 20:00 and were kept under a natural photoperiod of 12L: 12D. During the experiment, crabs were kept in a temperature of 19~22°C, a pH of 7.6~7.8, a dissolved oxygen concentration of at least 6.0 mg/L.

This study investigated whether different photoperiods affected the nervous system of E. sinensis. Crabs were exposed to one of three photoperiods including control (L: D = 12 h: 12 h, named CC), constant light (L: D = 24 h: 0 h, named LL) or constant darkness (L: D = 0 h: 24 h, named DD). Within the tank (Length: Width: Height = 37: 24.5: 11 cm) 10 crabs were randomly placed as a group, a light intensity of 100 lux was used, and treatments continued for 7 d. Then, crabs were frozen on ice, and the eyestalks were harvested and stored at -80°C until RNA isolation.

RNA extracted, read alignment and RNA-seq analysis

Total RNA was extracted using RNAisoTMPlus (TaKaRa) according to manufacturer’s instructions. Qualified total RNA was further purified with a RNeasy micro kit (Cat#74004, QIAGEN, GmBH, Germany) and RNase-Free DNase Set (Cat#79254, QIAGEN, GmBH, Germany). The RNA concentration and purity were determined using a Nanodrop2000, and RNA quantity was measured by denaturing formaldehyde agarose gel electrophoresis to examine integrity and measured with an Agilent Bioanalyzer 2100 (Agilent Technologies). Only high-quality eyestalk RNA samples were used for cDNA synthesis, so in the end, there were two qualified samples in each group for cDNA synthesis.

The RNA integrity numbers (RINs) of samples was required to be more than 7.0 when RNA-seq transcriptome libraries were prepared. Poly-T oligo-linked magnetic beads were used to purify PolyA mRNA from total RNA, and the intact mRNA was broken into fragments with bead washing buffer and metal bath. Aforementioned mRNAs were used as templates to synthesize first-strand complementary DNA (cDNA). Then, second-strand cDNA was synthesized using Resuspension Buffer, EtOH and Buffer. Next, cDNA, was end-repaired, a base was added to the 3’ end, and the cDNA was amplified with PCR. Finally, the constructed cDNA libraries were sequenced with Illumina HiSeqTM2000.

De novo assembly and annotation of the transcriptome

Raw reads from Illumina HiSeqTM2000 may contain sequencing primers and low-quality sequence, which can affect analytical quality. Therefore, the raw reads were cleaned through three steps: (a) linger sequences were discarded; (b) Q < 20 (Q = -10log10E) bases were removed; and (c) reads length shorter than 25 bp also discarded. Then, the clean reads were used for de novo assembly with Trinity (Minimum contig length> = 400 bp).

For annotation, the assembled final unigenes with screening condition (E-value<1e-5) were annotated using the NCBI protein nonredundant (Nr) and UniProt database. The top five unigenes that were compared with the CDD database by rpstblastn (http://www.biomedcentral.com/content/supplementary/1471-2105-13-42-s1/Cloud-BioLinux-Package-Documentation/docs/rpstblastn.html) were annotated using COG (Cluster of Orthologous Groups of proteins) classification. The unigenes were also classified with the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases.

Analysis of differentially expressed genes

The clean reads were mapped per sample to the corresponding gene for read counts. The expression levels of genes were determined by using Express software and the FPKM (fragments per kilobase of exon model per million mapped reads) method. Analysis of differential expression levels across samples was performed in edgeR, and significantly differentially expressed genes (DEGs) were calculated (q-value <= 0.05, Fold-change>=2). Significantly enriched terms were obtained by mapping DEGs to the corresponding GO term. Similarly, significantly enriched KEGG pathways were obtained.

Quantitative real-time PCR verification

Three DEGs were randomly chosen for quantitative real-time PCR (RT-qPCR) to validate the accuracy of the RNA-seq results. The total RNA which was returned after sequencing (three samples in each group) was used as template to synthesize first-strand cDNA with PrimeScript RT Master Mix (Cat. No. RR036A, TaKaRa). Primers were designed based on the DEGs with Primer Premier 5.0 (Table 1).

thumbnail
Table 1. Primers were selected for RNA-seq validation by RT-qPCR.

Three genes were named that TRP channel protein (TRP), pyruvate dehydrogenase complex (PDC), Phosphoenolpyruvate carboxykinase (PPC).

https://doi.org/10.1371/journal.pone.0210414.t001

The reaction system of RT-qPCR contained 5 μl of 2×SYBR Premix Ex TaqTM (TaKaRa, Japan), 1 μl of diluted first-strand cDNA, 3.4 μl of PCR-grade water, 0.2 μl of ROX Reference Dye II, and 0.2 μl of each primer [33]. The mixtures were run under the following conditions: 95°C for 30 s; followed by 40 cycles of 95°C for 5 s, and 60°C for 34 s; in addition, 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s using the ABI 7500 Real-Time PCR System (Life Technology, USA). We used 18S ribosome RNA (18S) as a reference gene. The DEGs expression levels were calculated using the 2-ΔΔCt method [23].

Statistical analysis

The data were expressed as mean ± S.D. values, and One-way analysis of variance was used for comparisons among the groups. A P value <0.05 indicated statistically significant difference.

Results

Transcriptome sequencing and reads assembly

We obtained 51,786,884 bp raw reads from CC, 58,383,804 bp raw reads from LL, and 59,026,446bp raw reads from DD. Clean reads were selected by excluding reads that did not conform to the requirements. Thus, the clean reads of CC, LL and DD were 48,772,584 bp, 53,943,281bp and 53,815,178 bp, respectively, (Table 2) and they were used for de novo assembly. Using this, we obtained 330,686 unigenes that contained the 161,380 unigenes that were 200–400 bp, 70,473 unigenes that were 400-600bp and 22,260 unigenes that were more than 2000 bp in length (Fig 1). Unigenes were used for further functional analysis.

thumbnail
Table 2. Summary of the RNA-Seq reads production after quality trimming and adapter clipping of CC, DD and LL groups.

https://doi.org/10.1371/journal.pone.0210414.t002

Unigenes functional annotation

COG classification was used to analyze unigenes function, and 184,664 unigenes were classified into 25 COG clusters, which were mainly enriched in signal transduction mechanisms (14.65%), general function prediction only (11.93%), cytoskeleton (10.16%), RNA processing and modification (9.11%) and transcription (6.86%) (Fig 2).

A total of 74,804 unigenes were obtained by GO annotation, which were assigned to three groups: molecular function, cellular component and biological process, and further divided into 57 categories. Most of them were concentrated to metabolic process, organelle and catalytic activity, etc. (Fig 3).

Finally, we analyzed the signaling pathways of unigenes with a KEGG KAAS online pathway analysis tool (http://www.kegg.jp/blastkoala/), a total of 3412 unigenes were mapped to five processes: metabolism (48.1%), environmental information processing (18.4%), genetic information processing (13.6%), cellular processes (13.1%), and organismal systems (6.8%). The maximum concentration of unigenes in metabolism was in the global and overview maps pathways, followed by carbohydrate metabolism, amino acid metabolism and lipid metabolism pathway. The pathways which contained the most unigenes of the other four processes were signal transduction, translation, transport and catabolism, and endocrine system, respectively (Fig 4).

DEGs analysis

To analyze the DEGs among the different samples, we found that the number of up-regulated DEGs between the LL vs. CC was 6,008, and down-regulated DEGs were 2,064. The up-regulated DEGs in DD vs. CC were 4,704 and down-regulated DEGs were 2,210. Between the LL vs. DD, the up-regulated and down-regulated DEGs were 2,817 and 4,341, respectively (Fig 5).

thumbnail
Fig 5. Numbers of DEGs in LL vs. CC, LL vs. DD, and DD vs. CC.

https://doi.org/10.1371/journal.pone.0210414.g005

The significant DEGs were assigned to 196, 180 and 212 KEGG pathways for LL vs. CC, LL vs. DD and DD vs. CC, respectively. The significantly enriched pathways (q-value<0.05) are shown in Table 3. We found that the phototransduction (Fig 6) and glycolysis / gluconeogenesis (Fig 7) pathways were only included in the DD vs. CC comparison; oxidative phosphorylation (Fig 8) pathway was included only in the LL vs. CC comparison, but no significantly enriched pathway was found in the LL vs. DD comparison. We also found that the TRP has a significant upregulation in the phototransduction pathway. The PDC and PPC genes were found in glycolysis/ gluconeogenesis pathway, but PDC was significantly down-regulated and PPC was significantly up-regulated. Next, we obtained the significantly enriched GO terms from the DEGs with a p-value<0.05. There we had 23, 16 and 42 significantly enriched terms in LL vs. CC, LL vs. DD and DD vs. CC, respectively. The terms of ‘metabolic process’, ‘cell’, and ‘catalytic activity’ were the most highly enriched in the three groups (S1, S2 and S3 Figs). Some of the DEGs clustered in ‘developmental process’ and ‘growth’ in the LL vs. CC and DD vs. CC comparisons (S1 and S3 Figs).

thumbnail
Table 3. Results of KEGG Pathway of Chinese mitten crab kept in LL vs. CC, LL vs. DD and DD vs. CC.

https://doi.org/10.1371/journal.pone.0210414.t003

RT-qPCR identified DEGs expression

Three DEGs were selected to verify the RNA-Seq results. The results of RT-qPCR confirmed that TRP was significantly up-regulated in the constant light and constant darkness groups when compared with the control group (Fig 9A). Additionally, PDC was significantly down-regulated in the LL and DD groups (Fig 9B). We also found that PPC had no expression in CC group, but significant expression levels in the LL and DD groups, confirming the results from the RNA-Seq (Fig 9C).

thumbnail
Fig 9. RT-qPCR analysis of 3 selected DEGs.

Relative expression levels are shown for (A) TRP; (B) PDC; (C) PPC. Asterisks above the bars indicates significant differences (P < 0.05) between control (CC) and LL or DD group. ——indicates that PPC had no expression in CC group.

https://doi.org/10.1371/journal.pone.0210414.g009

Discussion

This is the first report on the transcriptome of eyestalks ganglion in E. sinensis that were obtained with Illumina HiSeqTM2000 after exposure to different photoperiods. From the results we obtained, unigenes found to be primarily enriched in signal transduction mechanisms and metabolism with COG, GO and KEGG annotation analysis These results indicate that photoperiod conditions can also affect signal transduction and energy metabolism in E. sinensis as in other animals, such as E. superba, C. glacialis and P. clarkii [12, 13, 22].

Further analysis of DEGs, we found that the phototransduction pathway was significantly enriched in DEGs (Fig 6). In the phototransduction pathway, the gene TRP has a significant upregulation (Fig 6A) in LL vs. CC and in DD vs. CC. Researchers had previously confirmed that TRP was involved in phototransduction [28], energy balance control [26], and sensory functions [24]. In mammalian, the intrinsically photosensitive ganglion cells (ipRGCs) and melanopsin can be photoactivated, and then triggers G protein signaling pathway which can activate phospholipase C (PLC) and subsequent opening TRP [34]. As same time, TRP can be activated after light stimulus in Drosophila [35]. These studies indicated that TRP channel involved the light signal transduction in both vertebrates and invertebrates. So, when crabs were treated with different photoperiod, the expression level of TRP gene may be affected. However, the mechanism of how photoperiods affect TRP level is not clear. Some studies have proven that photoperiods can influence the Ca2+ concentration [36], and thereby influence the growth and survival of cells [37]. In crayfish, Ca2+ is related to photoreceptor sensitivity has been confirmed [15]. And TRP can be activated by Ca2+, to further participate in the process of phototransduction [38, 39]. In this study, we found that the expression level of TRP was changed by constant light or constant darkness. The above statement may provide an idea to explain how photoperiods affect TRP level. But the effects of TRP are on survival and growth in E. sinensis needs further research.

Additionally, we also found that guanine nucleotide–binding proteins (G proteins) gene had significant expression from the phototransduction pathway. G proteins which are made up of Gαq, Gβe and Gγe can be activated by the light rhodopsin [4042]. As mentioned earlier, G protein can stimulate PLC activity and further mediated phototransduction[34, 42]. In fact, light signal transduction is a complex process, which can be regulated by a variety of substance interactions including G proteins, cGMP, and K+ exchanger [43]. In this study, the Gαq gene is significantly upregulated in the phototransduction pathway when crabs were placed in constant darkness seven days (Fig 6). Given this, we can predict that a condition of constant darkness may improve the Gαq expression in eyestalks and therefore influence the phototransduction pathway.

As mentioned previously, photoperiod can affect energy metabolism. We also found that the glycogenolysis and glycolysis pathways in E. sinensis eyestalk ganglion were significantly enriched for DEGs, and the PDC and PPC genes were involved in these pathways (Fig 7). PDC is an abbreviation for pyruvate dehydrogenase complex, as controlled by the pyruvate dehydrogenase kinase (PDK). It catalyses the oxidative decarboxylation of pyruvate to the acetyl coenzymeA (acetyl-CoA) which produced from the glycogenolysis and glycolysis pathways [44]. And acetyl-CoA plays an important role in the tricarboxylic acid (TCA) cycle which produce chemical energy in the form of ATP [45]. Therefore, PDC is indispensable in the metabolic process. On the other hand, researchers have found PDC E1 subunit deficiency in zebrafish visual mutants, and visual defect can affect energy metabolism [46]. Moreover, PDC subunits participate in severe neurological dysfunctions, including optokinetic response [47]. Meanwhile, in this study, we also found that photoperiods can affect PDC gene expression of E. sinensis eyestalk ganglion. Therefore, this study suggests that photoperiod possibly regulates the PDC activity by eyestalk ganglion, and then affects the energy metabolism process.

PPC is the only enzyme that can catalyze the conversion of oxaloacetate to the intermediate phosphoenolpyruvate in the glycogenolysis/glycolysis progress [48]. From our RNA-Seq results, we found the upregulation of PPC gene expression in the glycogenolysis/glycolysis KEGG pathway (Fig 8). In mammalian, hepatic glucose production is controlled by the gluconeogenic enzyme activity in the liver, PPC [49]. And PPC can be regulated by the forkhead transcription factor (Foxo1) which is a member of a highly-conserved DNA binding motif protein family and plays an important role in insulin signaling transduction [50]. There was reported that Foxo 1 mRNA level can be affected on different photoperiod and subsequent regulating PPC level in rats [32]. In a short, photoperiods regulated the PPC expression level which was possibly controlled by Foxo 1 expression. But, there was no research on the effect of photoperiod on PPC activity in crustaceans. This study provides new evidence for the effects of photoperiod on PPC gene expression in E. sinensis. In addition, we found that the pyruvate kinase (PK) gene was up-regulated in the glycogenolysis/glycolysis pathway. PK can catalyze phosphate group transfer and phosphoenolpyruvate conversion to adenosine diphosphate [51, 52]. PK not only takes part in carbohydrate metabolism [5355] but is also involved in immunological reaction progress in crustaceans [56, 57]. A study on Drosophila suzukii found that photoperiod could regulate the PK levels, and then mediate the adult reproductive diapause [58]. Currently, there are few reports about the effects of different photoperiods on PK activity in crustaceans [59]. From our discoveries however, we found that different photoperiods regulate the expression of PK, and PK as an important enzyme, participates in energy metabolism. Therefore, photoperiod could affect energy metabolism progress though PK in crabs.

In the transcriptome, the oxidative phosphorylation pathway was significantly enriched in DEGs. The NADH peroxidase (NADH) and ATPase were included in this pathway (Fig 9). NADH, as a part of the pyridine nucleotide, can decompose hydrogen peroxide [60], which covers the cellular metabolism process [61]. For example, a study showed that temperature reduction could drop the oxidation matrix volume of NADH, and then decrease the oxidative phosphorylation capacity of mitochondria [62]. Another study also showed that radiation affected the oxidative phosphorylation process by oxidizing NADH to NAD+ [63]. However, NADH is important to mitochondrial metabolism in animals. There have been few reports about the effects of photoperiod on NADH activity [64]. In this study, constant light and darkness had upregulated the NADH peroxidase gene, which indicates that photoperiods affected metabolism activity of crabs. There was also a significant upregulation of ATPase gene activity in the LL and DD groups. ATPase is generated in the mitochondria and participates in oxidative phosphorylation in order to produce ATP [65, 66]. It is well known that ATP decomposition releases energy to promote cell growth, cell division and other vital movements [67]. Recently, research has shown that photoperiod significantly influences the Na+/K+-ATPase activity of gills of animals [68, 69]. The vacuolar H+- ATPase (V-ATPase) daily rhythm was also controlled by photoperiod in insects [70]. Together, these studies indicate that photoperiod may affect energy metabolism through affecting the activity of ATP enzyme.

In conclusion, photoperiod is important for animals, and man-made changes in photoperiods are bound to influence the physiological activity of animals. Among these physiological activities, the phototransduction and metabolism processes are significantly enriched in DEGs. Though RT-qPCR, we verified the upregulation of TRP in the phototransduction pathway and the downregulation of PDC and upregulation PPC in energy metabolism processes. Additionally, we also found that there is significantly increased gene expression of cGMP, PK, NADH peroxidase gene and ATPase genes in different pathways.

Supporting information

S1 Fig. GO annotation of DEGs in DD vs. CC.

https://doi.org/10.1371/journal.pone.0210414.s001

(ZIP)

S2 Fig. GO annotation of DEGs in DD vs. LL.

https://doi.org/10.1371/journal.pone.0210414.s002

(ZIP)

S3 Fig. GO annotation of DEGs in LL vs. CC.

https://doi.org/10.1371/journal.pone.0210414.s003

(ZIP)

References

  1. 1. Hamasaki K, Ogiso Y, Dan S, Kitada S. Survival, development and growth of larvae of the coconut crab, Birgus latro, cultured under different photoperiod conditions. Aquaculture Research. 2016;47(8):2506–17.
  2. 2. Mendoza-Vargas L, Báez-Saldaña A, Alvarado R, Fuentes-Pardo B, Flores-Soto E, Solís-Chagoyán H. Circadian rhythm in melatonin release as a mechanism to reinforce the temporal organization of the circadian system in crayfish. Invert Neurosci. 2017;17(2):6. pmid:28540583.
  3. 3. Yufera M, Perera E, Mata-Sotres JA, Calduch-Giner J, Martinez-Rodriguez G, Perez-Sanchez J. The circadian transcriptome of marine fish (Sparus aurata) larvae reveals highly synchronized biological processes at the whole organism level. Sci Rep. 2017;7(1):12943. pmid:29021622
  4. 4. Farhadi A, Jensen MA. Effects of photoperiod and stocking density on survival, growth and physiological responses of narrow clawed crayfish (Astacus leptodactylus). Aquaculture Research. 2016;47(8):2518–27.
  5. 5. Nogueira N, Sumares B, Andrade CAP, Afonso A. The effects of temperature and photoperiod on egg hatching success, egg production and population growth of the calanoid copepod, Acartia grani (Calanoida: Acartiidae). Aquaculture Research. 2018;49(1):93–103.
  6. 6. Fitzgibbon QP, Battaglene SC. Effect of photoperiod on the culture of early-stage phyllosoma and metamorphosis of spiny lobster (Sagmariasus verreauxi). Aquaculture. 2012;368–369:48–54.
  7. 7. Espinosa-Chaurand D, Vega-Villasante F, Carrillo-Farnés O, Nolasco-Soria H. Effect of circadian rhythm, photoperiod, and molt cycle on digestive enzymatic activity of Macrobrachium tenellum juveniles. Aquaculture. 2017;479:225–32.
  8. 8. Drezon‐Tepler M. Moult activity-related spawning success in the Mediterranean deep-water shrimp Aristeus antennatus (Decapoda:Dendrobranchiata): MARINE ECOLOGY PROGRESS SERIES; 1995.
  9. 9. Wadsworth T, Carriman A, Gutierrez AA, Moffatt C, Fuse M. Ecdysis behaviors and circadian rhythm of ecdysis in the stick insect, Carausius morosus. J Insect Physiol. 2014;71:68–77. pmid:25450561.
  10. 10. Matsuda H, Abe F, Tanaka S. Effect of photoperiod on metamorphosis from phyllosoma larvae to puerulus postlarvae in the Japanese spiny lobster Panulirus japonicus. Aquaculture. 2012;326–329:136–40.
  11. 11. Fanjul-Moles ML, Prieto-Sagredo J, López DS, Bartolo-Orozco R, Cruz-Rosas H. Crayfish Procambarus clarkii retina and nervous system exhibit antioxidant circadian rhythms coupled with metabolic and luminous daily cycles. Photochemistry & Photobiology. 2009;85(1):78–87. pmid:18643905.
  12. 12. Groeneveld J, Johst K, Kawaguchi S, Meyer B, Teschke M, Grimm V. How biological clocks and changing environmental conditions determine local population growth and species distribution in Antarctic krill (Euphausia superba): a conceptual model. Ecological Modelling. 2015;303:78–86.
  13. 13. Morata N, Søreide JE. Effect of light and food on the metabolism of the Arctic copepod Calanus glacialis. Polar Biology. 2013;38(1):67–73.
  14. 14. Picones A AH. Electrophysiological properties of crayfish retinal photoreceptors. Journal of Experimental Biology. 1990;150(1):111. pmid:24003478
  15. 15. Barriga-Montoya C, de la OMA, Fuentes-Pardo B, Gomez-Lagunas F. Desensitization and recovery of crayfish photoreceptors. Dependency on circadian time, and pigment-dispersing hormone. Comp Biochem Physiol A Mol Integr Physiol. 2017;203:297–303. pmid:27783925
  16. 16. Hui M, Song C, Liu Y, Li C, Cui Z. Exploring the molecular basis of adaptive evolution in hydrothermal vent crab Austinograea alayseae by transcriptome analysis. PLoS One. 2017;12(5):e0178417. pmid:28552991
  17. 17. Muto A, Orger MB, Wehman AM, Smear MC, Kay JN, Page-McCaw PS, et al. Forward genetic analysis of visual behavior in zebrafish. PLoS Genet. 2005;1(5):e66. pmid:16311625
  18. 18. Jr FR, Sanchez KG, Riley PP. Entrainment of the circadian rhythm in egg hatching of the crab Dyspanopeus sayi by chemical cues from ovigerous females. Biological Bulletin. 2016;230(1):15. pmid:26896174
  19. 19. Huang X, Ye H, Huang H, Yu K, Huang Y. Two beta-pigment-dispersing hormone (beta-PDH) isoforms in the mud crab, Scylla paramamosain: implication for regulation of ovarian maturation and a photoperiod-related daily rhythmicity. Anim Reprod Sci. 2014;150(3–4):139–47. pmid:25262380.
  20. 20. Loredo-Ranjel R, Fanjul-Moles ML, Escamilla-Chimal EG. Crustacean hyperglycemic hormone is synthesized in the eyestalk and brain of the crayfish Procambarus clarkii. Plos One. 2017;12(4):e0175046. pmid:28369112.
  21. 21. Hoelters L, O'Grady JF, Webster SG, Wilcockson DC. Characterization, localization and temporal expression of crustacean hyperglycemic hormone (CHH) in behaviorally rhythmic peracarid crustaceans, Eurydice pulchra (Leach) and Talitrus saltator (Montagu). General & Comparative Endocrinology. 2016;237:43–52. pmid:27468954.
  22. 22. Nelson-Mora J, Prieto-Sagredo J, Loredo-Ranjel R, Fanjul-Moles ML. Putative pacemakers in the eyestalk and brain of the crayfish Procambarus clarkii show circadian oscillations in levels of mRNA for crustacean hyperglycemic hormone. PLoS One. 2013;8(12):e83937. pmid:24391849.
  23. 23. Yang X, Huang G, Xu MJ, Zhang C, Cheng Y, Yang Z. Cloning and functional characterization of the DA2 receptor gene in Chinese mitten crab (Eriocheir sinensis). Plos One. 2018;13(3):e0193999. pmid:29554147.
  24. 24. Peng G, Shi X, Kadowaki T. Evolution of TRP channels inferred by their classification in diverse animal species. Molecular Phylogenetics & Evolution. 2015;84:145–57. pmid:24981559.
  25. 25. Bellono NW, Oancea E. UV light phototransduction depolarizes human melanocytes. Channels. 2013;7(4):243–8. pmid:23764911.
  26. 26. Moraes MN, Mezzalira N, Menaker M, Guler A, Castrucci AML. TRPV1 participates in the activation of clock molecular machinery in the brown adipose tissue in response to light-dark cycle. Biochimica Et Biophysica Acta. 2016;1864(2):324–35. pmid:27864077.
  27. 27. Nilius B, Owsianik G. The transient receptor potential family of ion channels. Genome Biology. 2011;12(3):218. pmid:21401968.
  28. 28. Glantz RM, Miller CS, Nässel DR. Tachykinin-related peptide and GABA-mediated presynaptic inhibition of crayfish photoreceptors. Journal of Neuroscience. 2000;20(5):1780–90. pmid:10684879.
  29. 29. Patil YN, Marden B, Brand MD, Hand SC. Metabolic downregulation and inhibition of carbohydrate catabolism during diapause in embryos of Artemia franciscana. Physiological & Biochemical Zoology Pbz. 2013;86(1):106–18. pmid:23303325.
  30. 30. Li S, Liu H, Bai J, Zhu X. Transcriptome assembly and identification of genes and SNPs associated with growth traits in largemouth bass (Micropterus salmoides). Genetica. 2017;145(2):175–87. pmid:28204905.
  31. 31. Heldmaier G, Klingenspor M, Werneyer M, Lampi BJ, Brooks SP, Storey KB. Metabolic adjustments during daily torpor in the Djungarian hamster. American Journal of Physiology. 1999;276(5 Pt 1):E896. pmid:10329984.
  32. 32. Xie X, Zhao B, Huang L, Shen Q, Ma L, Chen Y, et al. Effects of altered photoperiod on circadian clock and lipid metabolism in rats. Chronobiology International. 2017;34(8):1094–104. pmid:28708429.
  33. 33. Wei B, Yang Z, Wang J, Chen A, Shi Q, Cheng Y. Effects of dietary lipids on the hepatopancreas transcriptome of Chinese mitten crab (Eriocheir sinensis). Plos One. 2017;12(7):e0182087. pmid:28753670.
  34. 34. Hankins MW, Peirson SN, Foster RG. Melanopsin: an exciting photopigment. Trends Neurosci. 2008;31(1):27–36. pmid:18054803.
  35. 35. Tian Y, Hu W, Tong HW, Han JH. Phototransduction in Drosophila. Sci China-Life Sci. 2012;55(1):27–34. pmid:22314488.
  36. 36. Molnar T, Barabas P, Birnbaumer L, Punzo C, Kefalov V, Križaj D. Store-operated channels regulate intracellular calcium in mammalian rods. Journal of Physiology. 2012;590(15):3465–81. pmid:22674725.
  37. 37. Makino CL, Wen XH, Michaud N, Peshenko IV, Pawlyk B, Brush RS, et al. Effects of low AIPL1 expression on phototransduction in rods. Invest Ophthalmol Vis Sci. 2006;47(5):2185–94. pmid:16639031
  38. 38. Duan J, Li Z, Li J, Santa-Cruz A, Sanchez-Martinez S, Zhang J, et al. Structure of full-length human TRPM4. Proc Natl Acad Sci U S A. 2018;115(10):201722038. pmid:29463718.
  39. 39. Minke B. The History of the Drosophila TRP Channel: The birth of a new channel superfamily. Journal of Neurogenetics. 2010;24(4):216. pmid:21067449.
  40. 40. Waldo GL, Ricks TK, Hicks SN, Cheever ML, Kawano T, Tsuboi K, et al. Kinetic scaffolding mediated by a phospholipase C–β and Gq signaling complex. Science. 2010;330(6006):974. pmid:20966218.
  41. 41. Wang T, Montell C. Phototransduction and retinal degeneration in Drosophila. Pflugers Arch. 2007;454(5):821–47. pmid:17487503
  42. 42. And JM, Mayeenuddin LH. Purification, G protein activation, and partial amino acid sequence of a novel phospholipase C from squid photoreceptors. Biochemistry. 1998;37(25):9064. pmid:9636052.
  43. 43. Wensel TG. Signal transducing membrane complexes of photoreceptor outer segments. Vision Research. 2008;48(20):2052–61. pmid:18456304.
  44. 44. Suagee JK, Corl BA, Crisman MV, Wearn JG, McCutcheon LJ, Geor RJ. De novo fatty acid synthesis and NADPH generation in equine adipose and liver tissue. Comp Biochem Physiol B-Biochem Mol Biol. 2010;155(3):322–6. pmid:19962447.
  45. 45. Trivedi AK, Malik S, Rani S, Kumar V. Adaptation of oxidative phosphorylation to photoperiod-induced seasonal metabolic states in migratory songbirds. Comp Biochem Physiol A-Mol Integr Physiol. 2015;184:34–40. pmid:25636903.
  46. 46. Maurer CM, Schönthaler HB, Mueller KP, Neuhauss SC. Distinct retinal deficits in a zebrafish pyruvate dehydrogenase-deficient mutant. Journal of Neuroscience the Official Journal of the Society for Neuroscience. 2010;30(36):11962–72. pmid:20826660.
  47. 47. Taylor MR, Dowling JE. A zebrafish model for pyruvate dehydrogenase deficiency: Rescue of neurological dysfunction and embryonic lethality using a ketogenic diet. Proceedings of the National Academy of Sciences of the United States of America. 2004;101(13):4584–9. pmid:15070761.
  48. 48. Guo L, Zhang F, Zhang C, Hu G, Gao C, Chen X, et al. Enhancement of malate production through engineering of the periplasmic rTCA pathway in Escherichia coli. Biotechnology & Bioengineering. 2018;115(6):1571–80. pmid:29476618.
  49. 49. Altomonte J, Richter A, Harbaran S, Suriawinata J, Nakae J, Thung SN, et al. Inhibition of Foxo1 function is associated with improved fasting glycemia in diabetic mice. Am J Physiol-Endocrinol Metab. 2003;285(4):E718–E28. pmid:12783775.
  50. 50. O-Sullivan I, Zhang WW, Wasserman DH, Liew CW, Liu J, Paik J, et al. FoxO1 integrates direct and indirect effects of insulin on hepatic glucose production and glucose utilization. Nat Commun. 2015;6:12. pmid:26206703.
  51. 51. Biasibetti-Brendler H, Schmitz F, Pierozan P, Zanotto BS, Prezzi CA, de Andrade RB, et al. Hypoxanthine induces neuroenergetic impairment and cell death in striatum of young adult wistar rats. Molecular Neurobiology. 2017;55(5):4098–106. pmid:28593435.
  52. 52. Lee CH, Wu CC, Tyan YC, Yu WT, Huang ES, Yu HS. Identification of pyruvate kinase as a novel allergen in whiteleg shrimp (Litopenaeus vannamei) by specific-IgE present in patients with shrimp allergy. Food Chemistry. 2018;258:359–65. pmid:29655746.
  53. 53. Wang Z, Ren X, Gao B, Yu X, Zhang X, Zhang J, et al. Comparison of carbohydrate metabolism key enzymes in different generations of growth-selected Portunus trituberculatus families. Aquaculture. 2017;477:6–14. .
  54. 54. Ding ZL, Kong YQ, Li JF, Cao F, Zhang YX, Du ZY, et al. Growth and metabolic responses of juvenile Macrobrachium nipponense to different dietary carbohydrate levels. Aquaculture Nutrition. 2017;23(5):1136–44. .
  55. 55. Wang X, Li E, Qin JG, Wang S, Chen X, Cai Y, et al. Growth, body composition, and ammonia tolerance of juvenile white shrimp Litopenaeus vannamei fed diets containing different carbohydrate levels at low salinity. Journal of Shellfish Research. 1943;33(2):511–7.
  56. 56. Olguín‐León P, Enríquez‐Espinoza T, Mendoza‐Cano F, Encinas‐García T, Sánchez‐Paz A. Transcriptional profile of pyruvate kinase and pancreatic lipase encoding mRNAs of the Pacific whiteleg shrimp Penaeus vannamei during PstDV‐1 infection. Aquaculture Research. 2017;48(11):5587–94. .
  57. 57. Liu Y, Xu Y, Li S, Xu X, Gao Q, Yuan M, et al. Identification of proteome, antigen protein and antigen membrane protein from Spiroplasma eriocheiris. Letters in Applied Microbiology. 2017;65(5):395–402. pmid:28763106.
  58. 58. Zhai Y, Lin Q, Zhang J, Zhang F, Zheng L, Yu Y. Adult reproductive diapause in Drosophila suzukii females. Journal of Pest Science. 2016;89(3):679–88.
  59. 59. Gao XL, Zhang M, Tian HQ, Zheng JM, Li X, Song CB, et al. Effect of LED light quality on respiratory metabolism and activities of related enzymes of Haliotis discus hannai. Aquaculture. 2016;452:52–61.
  60. 60. Keirsse-Haquin J, Picaud T, Bordes L, Gracia AGD, Desbois A. Modulation of the flavin–protein interactions in NADH peroxidase and mercuric ion reductase: a resonance Raman study. European Biophysics Journal Ebj. 2017;47:205–23. pmid:28889232.
  61. 61. Moreno A, Kuzmiakglancy S, Rd RJ, Kay MW. Enzyme-dependent fluorescence recovery of NADH after photobleaching to assess dehydrogenase activity of isolated perfused hearts. Scientific Reports. 2017;7:45744. pmid:28361886
  62. 62. Brown JC, Gerson AR, Staples JF. Mitochondrial metabolism during daily torpor in the dwarf Siberian hamster: role of active regulated changes and passive thermal effects. Am J Physiol Regul Integr Comp Physiol. 2007;293(5):R1833. pmid:17804585.
  63. 63. Le M, Mcneill FE, Seymour CB, Rusin A, Diamond K, Rainbow AJ, et al. Modulation of oxidative phosphorylation (OXPHOS) by radiation- induced biophotons. Environmental Research. 2018;163:80–7. pmid:29427954.
  64. 64. Jia ZH, Wang MQ, Wang XD, Wang LL, Qiu LM, Song LS. Transcriptome sequencing reveals the involvement of reactive oxygen species in the hematopoiesis from Chinese mitten crab Eriocheir sinensis. Dev Comp Immunol. 2018;82:94–103. pmid:29307815.
  65. 65. Wang HW, Zhao WP, Liu J, Tan PP, Tian WS, Zhou BH. ATP5J and ATP5H proactive expression correlates with cardiomyocyte mitochondrial dysfunction induced by fluoride. Biological Trace Element Research. 2017;180(1):1–7. .
  66. 66. Bernardi P, Di LF, Fogolari F, Lippe G. From ATP to PTP and Back: A dual function for the mitochondrial ATP synthase. Circulation Research. 2015;116(11):1850–62. pmid:25999424.
  67. 67. Hardie DG. Keeping the home fires burning†: AMP-activated protein kinase. Journal of the Royal Society Interface. 2018;15(138):20170774. pmid:29343628.
  68. 68. Dutra BK, Fernandes FA, Oliveira GT. Carbofuran-induced alterations in biochemical composition, lipoperoxidation, and Na+/K+ATPase activity of Hyalella pleoacuta and Hyalella curvispina in bioassays. Comp Biochem Physiol C Toxicol Pharmacol. 2008;147(2):179–88. pmid:17936079.
  69. 69. Ban M. Rearing conditions to develop seawater tolerance in underyearling sockeye salmon smolt. Fisheries Science. 2006;72(1):128–35.
  70. 70. Bebas P, Cymborowski B, Giebultowicz JM. Circadian rhythm of acidification in insect vas deferens regulated by rhythmic expression of vacuolar H(+)-ATPase. Journal of Experimental Biology. 2002;205(1):37–44. pmid:11818410.