Transcriptome analysis of genes and pathways associated with metabolism in Scylla paramamosain under different light intensities during indoor overwintering
BMC Genomics volume 21, Article number: 775 (2020)
Scylla paramamosain is one of the commercially crucial marine crustaceans belonging to the genus Scylla, which is commonly distributed along the coasts of China, Vietnam, and Japan. Genomic and transcriptomic data are scarce for the mud crab. Light intensity is one of the ecological factors that affect S. paramamosain during indoor overwintering. To understand the energy metabolism mechanism adapted to light intensity, we analyzed the transcriptome of S. paramamosain hepatopancreas in response to different light intensities (0, 1.43, 40.31 μmol·m− 2·s− 1).
A total of 5052 differentially expressed genes were identified in low light group (LL group, 3104 genes were up-regulated and 1948 genes were down-regulated). A total of 7403 differentially expressed genes were identified in high light group (HL group, 5262 genes were up-regulated and 2141 genes were down-regulated). S. paramamosain adapts to different light intensity environments through the regulation of amino acids, fatty acids, carbon and energy metabolism. Different light intensities had a strong impact on the energy generation of S. paramamosain by influencing oxygen consumption rate, aerobic respiration, glycolysis/gluconeogenesis pathway, the citrate cycle (TCA cycle) and fatty acid degradation.
Low light is more conducive to the survival of S. paramamosain, which needs to produce and consume relatively less energy to sustain physiological activities. In contrast, S. paramamosain produced more energy to adapt to the pressure of high light intensities. The findings of the study add to the knowledge of regulatory mechanisms related to S. paramamosain metabolism under different light intensities.
There are many different stages in the life history of aquatic animals. Overwintering is an important life stage for aquatic animals living in some areas to adapt to cold winter conditions . The process of overwintering determines the quality of parents and is a key factor in seedling development in the coming year . There are many factors affecting the environment of aquatic animals during overwintering, such as salinity, temperature, food availability, and light.
In the aquatic environment, absorption and reflection of incident light is easily changed due to the presence of plankton, suspended particulate matter and soluble organic substances , which influences the light intensity of water. Therefore, light is one of the most variable water quality factors that directly or indirectly affects behavior [4, 5], survival and growth [6, 7], skin color [8, 9], digestion and immunity , fatty acid composition , and oxidative stress  in aquatic animals. Metabolism is an extremely important activity and has been studied extensively. Nalbach et al. demonstrated that the activity of the neuron-secreting cells in the eye stalk and the thoracic brain of crabs are closely related to light intensity . Increased activity of neurosecretory cells under high light intensity promoted the activity of hormone-promoting hormones and metabolism-related hormones. Most shrimp and larvae react positively to light in the larval stage and a negative reaction to light in the adult stage . Appropriate light conditions can improve the feeding rate of shrimp and crab larvae and promote growth and metamorphosis . This phenomenon is particularly evident in the early developmental stages of shrimp and crab larvae.
S. paramamosain is an important marine crab species in China. This species has fast growth characteristics, large sized individuals, highly adaptable, delicious meat and high nutritional value . S. paramamosain has played an important role in wild fisheries and aquaculture over the past few decades in China . In 2018, the yield of S. paramamosain increased to 157,712 tons, accounting for 53.68% of the total marine crab aquaculture industry in China . Although the culture of S. paramamosain has grown continuously for years, there are still very low yields, only 15–60 kg/667 m2 a year, which is unable to meet ever increasing market demand [19,20,21]. It is of urgent importance to improve the artificial breeding and breeding technology of S. paramamosain by improving the overwintering environment . Therefore, understanding light intensity conditions required during the overwintering period are the focus of this study.
In recent years, with the development of molecular sequencing technology, non-parametric transcriptome sequencing technology has gradually become the main means of studying molecular mechanisms in biological research . Transcriptome sequencing is the study of all mRNAs that a particular cell can transcribe under a certain functional state. For species without reference genomes, splicing small fragments out of Unigene (non repetitive sequence genes, a transcriptome database of the National Center for Biotechnology Information) and constructing reference sequences for subsequent analysis is an effective means of studying the molecular mechanisms and regulatory networks of non-parametric species. Analysis by Ding et al. identified key metabolic changes in the liver of the large yellow croaker (Larimichthys crocea) in response to acute hypoxia by transcriptome . Wang et al. revealed the potential influencing mechanism of dietary astaxanthin on growth and metabolism in Litopenaeus vannamei by transcriptome . Research carried out by Zou et al. described the change of aroma volatile metabolism in tomato (Solanum lycopersicum) fruit under different storage temperatures and 1-MCP treatment by transcriptome . While there has been a lot of research on S. paramamosain, no studies have reported the applications of the transcriptome technique for investigating molecular mechanisms associated with metabolism under different light intensities. There has been no research to date on the molecular mechanism of light intensity on S. paramamosain metabolism.
In this study, the effects of different light intensities on S. paramamosain were investigated using transcriptome analysis, in particular the metabolic changes of gene expression under different light intensities. This study will assist in explaining the potential molecular mechanisms of S. paramamosain under different light intensities. These data may help to provide the necessary theoretical basis and technical support for breeding S. paramamosain.
Differentially expressed genes (DEGs) in the hepatopancreas of S. paramamosain
A total of nine samples were used in this analysis. The Q30 of raw data for each sample was distributed from 93.00 to 94.09%, the effective data was distributed from 6.97 to 8.06 G, and the average GC content was 45.90% (Additional file 1: Table S3). Unigene 54,537 pieces were spliced with a total length of 60,890,086 bp and an average length of 1116 bp (Additional file 1: Table S2). As the genome sequencing of S. paramamosain has not yet been elucidated, after reads filtering, Trinity was used to perform de novo assembly with clean reads. Tgicl was used on cluster transcripts to remove abundance and obtain Unigenes.
The reads were compared to Unigene and the alignment was 90.63–93.10% (Additional file 1: Table S3). Stochasticity evaluation of sequencing data showed that the distribution of reads in the various parts of the gene was relatively uniform, indicating that the randomness of the break was good. The RPKM values of the individual genes determined under each 5% increment were compared to the expression levels of the final corresponding genes. The results differed by less than 15%, and therefore, sequencing data met quantitative requirements. The distribution of FPKM mean values for Unigenes ranged from 11.50 to 13.33 (Additional file 1: Fig. S1). Inter-sample correlation tests confirmed that the study was reliable and the sample selection was reasonable. Principal component analysis showed that each group of samples was distributed in different regions, and the same group of samples was concentrated spatially. Cluster analysis also showed that similar samples grouped together (Fig. 1).
Finally, high-quality transcripts were obtained and used as reference sequences. For database annotations of Unigenes, 16,049 (29.43%) genes were annotated into the NR library, 11,876 (21.78%) genes were annotated into the Swissprot library, 7881 (14.45%) genes were annotated into the KEGG library, 10,586 (19.41%) were annotated into the KOG library, 13,351 (24.48%) genes were annotated into the eggNOG library, 11,204 (20.54%) genes were annotated into the GO library, and 284 (0.52%) genes were annotated into the Pfam library (Additional file 1: Table S4). For functional annotation results, there were 54,287 SSRs predicted, including 25,440 Unigenes containing SSRs, 13,056 Unigenes containing more than one SSR and 11,842 composite SSRs (Additional file 1: Table S5). A total of 29,221 CDS sequences were predicted, of which 16,111 were predicted by the database comparison method and 13,110 were predicted by ESTScan. A total of 1672 Unigenes were annotated to the transcription factor database and distributed across 63 families.
To compare the difference between the experimental group and control group in terms of the change in global proteomes in the hepatopancreas, volcano plots were created of P-value (−log10 P-value) versus the log2fold change for each gene analyzed. Accordingly, relative to the control group, the LL and HL group exhibited a change of a certain proportion of proteins (Fig. 2). A total of 5052 differentially expressed Unigenes were obtained between the LL group and the control, including 3104 up-regulated Unigenes and 1948 down-regulated Unigenes, while 7403 differentially expressed Unigenes were obtained between the HL group and the control, including 5262 up-regulated Unigenes and 2141 down-regulated Unigenes (p value< 0.05 & |log2FC| > 1) (Fig. 3a). Interestingly, 2482 differentially expressed Unigenes were intersected (Fig. 3b).
It is well known that carbon metabolism plays a key role in metabolic regulation. An important pathway in carbon metabolism is oxidative phosphorylation (ko00190), in which all significantly different genes are down-regulated. In the LL group, the three most highly down-regulated genes were NADH-ubiquinone oxidoreductase chain 3 (ID: P18930; 0.21-fold), NADH-ubiquinone oxidoreductase chain 2 (ID: P34848; 0.24-fold), and V-type proton ATPase subunit F (ID: Q1HQK8; 0.31-fold) (Table 1). In the HL group, the three most highly down-regulated genes were NADH-ubiquinone oxidoreductase chain 2 (ID: P34848; 0.10-fold), cytochrome c oxidase subunit 1 (ID: P00399; 0.17-fold), and NADH-ubiquinone oxidoreductase chain 4 (ID: Q34048; 0.17-fold) (Table 1). Two important pathways in energy metabolism are glycolysis / gluconeogenesis (ko00010) and the citrate cycle (TCA cycle, ko00020). In this process, only pyruvate dehydrogenase E1 component (ID: P49432) was significantly up-regulated. The fold change was 3.26 in the LL group and 2.64 in the HL group. In the LL group, the three most highly down-regulated genes were glyceraldehyde-3-phosphate dehydrogenase (ID: P56649; 0.25-fold), hexokinase type 2 (ID: Q9NFT7; 0.28-fold) and enolase (ID: P56252; 0.34-fold). In the HL group, the three most highly down-regulated genes were triosephosphate isomerase B (ID: Q90XG0 0.33-fold), aconitate hydratase (ID: Q99KI0; 0.37-fold), and glyceraldehyde-3-phosphate dehydrogenase (ID: P56649; 0.38-fold) (Table 1). Glycine, serine and threonine metabolism (ko00260) and beta-Alanine metabolism (ko00410) are important pathways in amino acid metabolism. 2-amino-3-ketobutyrate coenzyme A ligase (ID: O75600; LL: 2.13-fold, HL: 2.49-fold) and cystathionine gamma-lyase (ID: Q8VCN5; LL: 4.08-fold, HL: 16.10-fold) were up-regulated in both pathways. Betaine--homocysteine S-methyltransferase 1 (ID: Q5XGM3) was down-regulated in the LL group (0.53-fold) and up-regulated in the HL group (7.34-fold). In the LL group, the three most highly down-regulated genes were cytosolic non-specific dipeptidase (ID: Q96KP4; 0.13-fold), peroxisomal sarcosine oxidase (ID: Q29RU9; 0.16-fold), and aldehyde dehydrogenase (ID: P81178; 0.25-fold). In the HL group, the three most highly down-regulated genes were serine--pyruvate aminotransferase (ID: P41689; 0.19-fold), cytosolic non-specific dipeptidase (ID: Q96KP4; 0.21-fold) and D-3-phosphoglycerate dehydrogenase (ID: A5GFY8; 0.30-fold) (Table 1). Two important pathways in lipid metabolism are fatty acid elongation (ko00062) and fatty acid degradation (ko00071). Lysosomal thioesterase PPT2 (ID: O70489; LL: 3.47-fold, HL: 3.12-fold) and lysosomal thioesterase PPT2-B (ID: Q6GNY7; LL: 1.97-fold, HL: 3.00-fold) were up-regulated in both the LL and HL groups. In the LL group, fatty aldehyde dehydrogenase (ID: P47740) was the most down-regulated gene (0.31-fold), followed by alcohol dehydrogenase class-3 (ID: P79896, 0.38-fold) and hydroxyacyl-coenzyme A dehydrogenase (ID: Q16836; 0.44-fold) (Table 1). In the HL group, the three most down-regulated genes were hydroxyacyl-coenzyme A dehydrogenase (ID: Q16836; 0.35-fold), alcohol dehydrogenase class-3 (ID: P79896; 0.40-fold) and enoyl-CoA delta isomerase 2 (ID: Q9WUR2; 0.41-fold) (Table 1).
Gene Ontology (GO) terms could be divided into three ontologies: cellular components, molecular function, and biological processes. In the LL group, 49 GO terms were assigned to the up-regulated group and 54 to the down-regulated group (Fig. 4a). In the HL group, 48 GO terms were assigned to the up-regulated group and 52 to the down-regulated group (Fig. 4b). A total of 23 processes were identified in the biological process category, comprising three biological processes (cellular process, metabolic process and biological regulation) as the most strongly affected in the hepatopancreas of S. paramamosain under different light intensities (Fig. 4). In the cellular component category, cells, cell parts, and membranes were most involved (Fig. 4). In the molecular function category, two items comprising binding and catalytic activities were most involved (Fig. 4).
Metabolic related DEG pathway analysis
KEGG pathway classification (Fig. 5) was performed according to the KEGG database website. All DEGs were graded into five categories according to their biological function, including cellular processes (LL, 150; HL, 175), environmental information processing (LL, 119; HL, 121), genetic information processing (LL, 198; HL, 309), metabolism (LL, 426; HL, 519) and organismal systems (LL, 227; HL, 247) (Fig. 5).
The pathways related to metabolism were subdivided into twelve subsets (level 2): xenobiotic biodegradation and metabolism (LL, 23; HL, 25), nucleotide metabolism (LL, 18; HL, 26), metabolism of terpenoids and polyketides (LL, 5; HL, 7), metabolism of other amino acids (LL, 24; HL, 26), metabolism of cofactors and vitamins (LL, 37; HL, 34), lipid metabolism (LL, 63; HL, 66), glycan biosynthesis and metabolism (LL, 27; HL, 36), global and overview maps (LL, 42; HL, 57), energy metabolism (LL, 55; HL, 99), carbohydrate metabolism (LL, 72; HL, 75), biosynthesis of secondary metabolism (LL, 4; HL, 3) and amino acid metabolism (LL, 56; HL, 65) (Fig. 5a&b and Additional file 1: Table S6).
Functional enrichment was also performed on DEGs according to the above KEGG pathway classification. The top 20 pathways showed that there were 15 and 12 pathways directly related to metabolism in the LL and HL groups, respectively. The most common metabolic pathways were biosynthesis of amino acids (ko01230), degradation of aromatic compounds (ko01220), carbon metabolism (ko01200), methane metabolism (ko00680), oxidative phosphorylation (ko00190) and glycolysis/Gluconeogenesis (ko00010) (Fig. 5c&d). The pathways related to metabolism in the LL group were pentose and glucuronate interconversions (ko00040), drug metabolism-cytochrome P450 (ko00982), starch and sucrose metabolism (ko00500), metabolism of xenobiotics by cytochrome P450 (ko00980), porphyrin and chlorophyll metabolism (ko00860), retinol metabolism (ko00830), glycine, serine and threonine metabolism (ko00260), steroid hormone biosynthesis (ko00140) and ascorbate and aldarate metabolism (ko00053) (Fig. 5c). The pathways related to metabolism in the HL group were caprolactam degradation (ko00930), carbon fixation in photosynthetic organisms (ko00710), aminobenzoate degradation (ko00627), geraniol degradation (ko00281), fatty acid elongation (ko00062) and the citrate cycle (TCA cycle) (ko00020) (Fig. 5d).
Seven important metabolic pathways were selected for further analysis: Glycine, serine and threonine metabolism (ko00260), beta-Alanine metabolism (ko00410), fatty acid elongation (ko00062), fatty acid degradation (ko00071), oxidative phosphorylation (ko00190), glycolysis / gluconeogenesis (ko00010) and the TCA cycle (ko00020) (Table 1). We focused on these pathways and the differential genes they contain.
Validity of DEGs in transcriptomic data
To validate the transcriptome results, qRT-PCR was used to check the transcript levels in six identified DEGs (PPT2, ODPB, KBL, ATP5H, G3P, and SPYA). As indicated by qRT-PCR, in the six genes selected in a random manner, the majority of the genes exhibited a similar transcription level expression to the transcriptome results in different light intensities, with the exception of G3P in the HL group and SPYA in the LL group (Fig. 6). The qRT-PCR results showed a high level of consistency with the transcriptome data obtained from transcriptome, indicating that the transcriptome results were reliable.
Light is an important and variable factor in both natural and farmed ponds . S. paramamosain living in water is easily impacted by light intensity and the impacts of different light intensities on aquatic crustaceans have received some attention in recent years. One of the main functions of the hepatopancreas is metabolism [28, 29]. This study focuses on using transcriptome technology to investigate metabolism gene expression changes in the hepatopancreas of S. paramamosain in different light intensities. Four functional categories were selected for evaluating the dataset, mainly referring to carbohydrate metabolism, energy metabolism, amino acid metabolism, as well as lipid metabolism. The following sections will focus on discussing the biological relevance exhibited by the above-mentioned seven metabolic pathways.
Effect of light intensity on energy metabolism of S. paramamosain
Among the energy metabolism pathways that were enriched, oxidative phosphorylation (ko00190) was the most affected. Five types of genes (NADH dehydrogenase, NADH-ubiquinone oxidoreductase, ATP synthase, cytochrome c oxidase, inorganic pyrophosphatase) were significantly down-regulated, and the fold change of the HL group was significantly higher than that of the LL group.
Mitochondria are the driving force for all living things to carry out life activities, and are the main sites for intracellular oxidative phosphorylation and formation of ATP which drives cellular functions . Oxidative paosphorylation is a coupling reaction between the energy released by the body, ATP, and inorganic phosphate synthesis through the respiratory chain . NADH-ubiquinone oxidoreductase and cytochrome c oxidase (COX) participate in the transfer of electrons to oxygen, and thus, are considered the key components of the respiratory chain . Previous studies have shown that light affects oxygen consumption in the respiratory chain of crustaceans. Crayfish (Procambarus clarkia and Procambarus digueti) have a reduced oxygen consumption rate during long light periods . The oxygen consumption rate of Chinese prawn and Orconectes nais was not affected by the light cycle . The larvae of Macrobrachium rosenbergii have a high oxygen consumption rate in the dark, while lobster larvae and Chinese prawn exhibit high oxygen consumption under light conditions . The oxygen consumption rate of P. vannamei was not affected by light intensity . The average oxygen consumption rate of Penaeus chinensis in different shades of light was different under dark and light conditions. In this study, the change of NADH-ubiquinone oxidoreductase and COX genes showed that different light intensities affected the respiratory process of S. paramamosain. Significant down-regulation of NADH-ubiquinone oxidoreductase and COX genes indicated the reduced transfer of electrons to oxygen and a reduced oxygen consumption rate under light conditions. In addition, the oxygen consumption rate under high light intensity was lower than that under low light intensity.
The reaction catalyzed by the inorganic pyrophosphatase produces inorganic pyrophosphate, which relies on several important nucleotide triphosphates and is necessary for the synthesis of RNA, DNA, proteins, and lipids. While providing Pi for biomolecules, the inorganic pyrophosphatase can synthesize ATP, a terminal product of cellular energy metabolism . NADH dehydrogenase which is produced by reduced coenzyme NADH is involved in many biochemical metabolic processes, such as glycolysis, the Krebs cycle, and beta oxidation . This coenzyme contains electrons with a high electrode potential, which release a lot of energy when oxidized . The results of the current study showed that inorganic pyrophosphatase genes in both experimental groups were significantly down-regulated, indicating that the amount of ATP synthesis was significantly reduced under light conditions. This conclusion can also be illustrated from the negative regulation of NADH dehydrogenase complex and ATP synthase subunit enrichment under light conditions in the current study.
Effect of light intensity on carbohydrate metabolism of S. paramamosain
The hepatopancreas is the main organ of carbohydrate metabolism in crustaceans, and facilitates the maintenance of blood glucose levels through four metabolic compounds and the metabolic pathways they participate in; glycosylation, glycogenolysis, glycolysis and gluconeogenesis [40, 41]. According to enrichment analysis of the KEGG pathway, two important pathways related to carbohydrate metabolism were found: glycolysis/gluconeogenesis (ko00010) and the the citrate cycle (TCA cycle) (ko00020). Many of the genes in these pathways changed significantly, indicating that carbohydrate metabolism was of great significance for S. paramamosain to adapt to different light conditions.
Because the hepatopancreas stores a limited amount of glucose or glycogen, once this has been used up the body needs to synthesize glucose through gluconeogenesis to maintain normal energy requirements. During the process of gluconeogenesis, glucone-6-phosphatase catalyzes the conversion of non-sugar substances (such as lactic acid, pyruvate, propionic acid, glycerol, etc.) into glucose or glycogen [42, 43]. Triosephosphate isomerase, an important enzyme of the glycolysis/gluconeogenesis pathway, interconverts dihydroxyacetone phosphate and glyceraldehyde-3-phosphate dehydrogenase (G3P) to prevent the accumulation of cell-toxic dihydroxyacetone phosphate . Under light conditions, the glycolysis /gluconeogenesis-related enzymes, such as triosephosphate isomerase, G3P, L-lactate dehydrogenase, acetyl-coenzyme A (acetyl-coA) and glucone-6-phosphatase were significantly down-regulated, indicating that the synthesis pathway of glycolysis /gluconeogenesis was inhibited to some extent. It can be seen from the oxidative phosphorylation pathway of energy metabolism that aerobic respiration is more vigorous under light conditions. This process produces carbon dioxide and water, releases energy and synthesizes large amounts of ATP by completely oxidizing and breaking down organic matter (usually decomposing glucose). This may be the main reason for inhibiting the glycolysis/gluconeogenesis pathway. In this study, pyruvate kinase gene expression was up-regulated by 9.7-fold in the HL group, which promoted the production of pyruvate which can be converted into acetyl-coA and enter the fatty acid biosynthesis pathway .
The citrate cycle (TCA cycle) is the third stage of catabolism, in which acetyl-coA is completely oxidized and decomposed into carbon dioxide . The related genes, except for the pyruvate dehydrogenase (PDH) gene, were significantly down-regulated. PDH is one of the key enzymes in the central metabolic system. When pyruvate enters the TCA cycle, it first forms acetyl-coA under the action of PDH in the mitochondrial inner membrane, and then enters the TCA cycle . The PDH gene was significantly up-regulated by 2.6-fold in the HL group and 3.2-fold in the LL group. However, there was no significant difference in acetyl-coA gene between the HL and LL groups, indicating that not enough pyruvate entered the TCA cycle to produce acetyl-coA. From the experimental results, the main genes that were down-regulated in this cycle were catalytic enzymes genes, such as malate dehydrogenase, isocitrate dehydrogenase, ATP-citrate synthase, aconitate hydratase, and fumarate hydratase. This indicates that light conditions limit the TCA cycle of the organism and the specific mechanism behind this requires further study.
Effect of light intensity on amino acid metabolism in S. paramamosain
In terms of amino acid metabolism, two genes were significantly up-regulated in the two groups: 2-amino-3-ketobutyrate coenzyme A ligase and cystathionine gamma-lyase. Betaine--homocysteine S-methyltransferase 1 and sarcosine dehydrogenase were significantly up-regulated in the HL group, but significantly down-regulated in the LL group. From the results of this study, it can be seen that the regulation patterns of amino acids under different light intensities have similarities and differences. Cystathionine is a metabolic intermediate of sulfur-containing amino acids, which is involved in the conversion between cysteine and methionine. Cystathionine gamma-lyase acts on cystathionine so that cystathionine eliminates pyruvate and NH3 and produces homocysteine . The fold change of GCL in the HL group (16.10-fold change) was significantly higher than that in the LL group (4.08-fold change), indicating that the HL group produced more homocysteine. Homocysteine has been widely studied in human diseases, but it has not been studied in crustaceans. Sarcosine usually exists in the form of creatine phosphate, and organisms mainly rely on ATP to provide energy under stress . However, ATP reserves in the body are very small, so continuous synthesis is needed, and creatine phosphate can promote the synthesis of ATP. Creatine dehydrogenase is an enzyme that catalyzes chemical reactions, and its elevation significantly reduces the level of creatine. The sarcosine content in high light conditions was significantly higher than that in low-light conditions, which may be related to the fact that S. paramamosain requires more energy under high light conditions.
Some studies have shown that some glycosylated amino acids (such as alanine, glycine and serine) are decomposed in the hepatopancreas, and the carbon skeleton can be converted into metabolites such as pyruvate, ketoglutaric acid, succinyl coA, fumaric acid and oxaloacetic acid through catabolism under stress conditions, thus, affecting amino acid metabolism . However, under opposite environmental conditions, these metabolites are restricted from entering the TCA cycle to produce ATP for energy or glucose through the gluconeogenic pathway. This conclusion was confirmed by a significant decrease of the gluconeogenic pathway. The combination of amino acid catabolism and gluconeogenic pathways may be a mechanism for adaptation to light intensity. However, the specific role of some amino acids needs further study to fully understand amino acid dynamics in different light intensities.
Effect of light intensity on lipid metabolism in S. paramamosain
Lipids are the most important substrates for maintaining tissue function and are also important energy stores, providing energy for organisms . In this study, KEGG analysis screened two key metabolic pathways related to the synthesis of fatty acids including fatty acid elongation and fatty acid degradation. A total of 14 DEGs were found to be involved in these two metabolic pathways, indicating that light intensity significantly affects lipids metabolism.
Elongation of fatty acids is an important biological process in the biosynthesis of polyunsaturated fatty acids [51, 52]. Polyunsaturated fatty acids play an important role in maintaining the relative fluidity of cell membranes and ensuring the normal physiological function of cells. In this study, four genes were found to be significantly correlated with the content of polyunsaturated fatty acids: lysosomal thioesterase PPT2, lysosomal thioesterase PPT2-B, 3-ketoacyl-CoA thiolase and hydroxyacyl-coA. The lysosomal thiesterase gene plays an important role in fatty acid elongation, but the specific mechanism has not been reported and needs further study. Thiesterase plays an important catalytic role in the synthesis and degradation of fatty acids. Thiesterases can be divided into two groups according to their different functions, in which 3-ketoacyl-CoA thiolase plays an important role in fatty acid decomposition . The significant down-regulation of the 3-ketoacyl-CoA thiolase gene indicates that the solution rate of fatty acids decreases and more fatty acids might be accumulated under light conditions. In addition, different genes in fatty acid degradation pathways differed between high-light and low-light conditions, but were significantly down-regulated in both. This result indicates that fatty acid degradation was somewhat inhibited in the presence of light, but the specific mechanisms of these actions requires further study.
In this study, we used RNA-sequencing to construct gene expression profiles related to metabolism mechanisms in S. paramamosain. To our knowledge, this is the first transcriptome study of these mechanisms in S. paramamosain. The results of this study revealed that glycine, serine and threonine metabolism, beta-Alanine metabolism, fatty acid elongation, fatty acid degradation, oxidative phosphorylation, glycolysis/gluconeogenesis, and the TCA cycle are involved in metabolism mechanisms in S. paramamosain. In conclusion, energy metabolism, amino acid metabolism, carbohydrate metabolism, and lipid metabolism are important means of regulation to cope with different light intensities. These results contribute to the existing knowledge of energy metabolism mechanisms in aquatic animals when adapting to different light intensities. These data may help to understand the molecular mechanisms of metabolism under different light intensities in S. paramamosain from a new perspective. We also suggested that low light intensity could be used for indoor overwintering.
Source, breeding and grouping of Scylla paramamosain
S. paramamosain individuals were collected from the sea near Xiangshan City (121 57 14 E, 29 28 29 N) and reared in crab apartment systems. The crabs were acclimatized in the carb apartments for 10 days prior to experimentation. And the crabs (290 ± 40 g) were randomly placed into nine groups, which were placed in nine crab apartments. Three groups of crab apartments comprised an experimental group: HL group (high light illumination with a photon flux density of about 40.31 μmol·m− 2·s− 1), LL group (low light illumination with a photon flux density of about 1.43 μmol·m− 2·s− 1), and the control group (full darkness). The light sources of each experimental group were white LED lights, simulating natural light with full spectrum. Peak light intensity was measured using a spectroradiometer. The light cycle was controlled by a timer (10 L: 14D), which turned the lights on at 07:00 every day and turned them off at 17:00. During the experiment period, the crabs were fed once daily (at 17:00) with Ruditapes philippinarums. The food debris and feces were cleaned at the outfall under the tanks at 07:00 the next day. The content of moisture, crude protein, crude fat and ash in whole organs (wet basis) of Ruditapes philippinarums were 77.21, 16.40, 1.46 and 3.01%, respectively . Apart from light intensities, all other daily feeding and management operations were the same for all experimental groups during the trial.
The experiment was conducted in December 2018–March 2019 for a period of 4 months. After the experiment, the hepatopancreases of three crabs were collected and pooled to generate a single pooled sample. A total of 9 hepatopancreases (3 pooled samples) were obtained from each treatment group and used for further analyses. All samples were stored at − 80 °C.
RNA extraction, library construction and sequencing
Total RNA was extracted with Trizol Reagent (Invitrogen, Shanghai, China) from 100 mg of hepatopancreas tissue following the manufacturer’s instructions. The sample was detected in 1.5% agarose gel electrophoresis. All RNA samples were of high quality (OD260/280 = 2.10–2.15, OD260/230 ≥ 2.0).
After digesting the DNA with DNase, eukaryotic mRNA was enriched using magnetic beads with Oligo (dT). A disruption reagent was added to break the mRNA into short fragments. Using the interrupted mRNA as a template, single-strand cDNA was synthesized using a six-base random primer, and then a two-stranded reaction system was used to synthesize double-stranded cDNA, which was purified using the kit. The purified double-stranded cDNA was subjected to end repair, a tail was added and the sequencing linker was ligated. The fragment size was then selected, and finally PCR amplification was performed. The constructed library was qualified with an Agilent 2100 Bioanalyzer and sequenced using an Illumina HiSeq X Ten sequencer to generate 150 bp double-ended data.
Gene annotation and classification
After removing adaptor sequences, ambiguous N nucleotides (with a ratio of N > 10%) and low-quality sequences (with a quality score < 5), the remaining clean reads were assembled using Trinity software as described for de novo transcriptome assembly without a reference genome, and the longest copies of redundant transcripts were regarded as unigenes. Using diamond software, Unigenes were compared to the NR (NCBI non-redundant protein sequences, http://www.ncbi.nlm.nih.gov/), KOG (clusters of euKaryotic Orthologous Groups, http://www.ncbi.nlm.nih.gov/COG/), GO (the Gene Ontology, http://www.geneontology.org/), Swiss-Prot (http://www.ebi.ac.uk/uniprot/), and KEGG (Kyoto Encyclopedia of Genes and Genomes, http://www.genome.jp/kegg/) databases, and PFMER software was used to compare the Pfam database for unigene functional analysis.
Differential gene identification, enrichment, and pathway analysis
Gene expression levels were positively correlated with their abundance. In transcriptome sequencing analysis, Unigene expression levels could be estimated by counting the sequencing reads of Unigenes. The reads count was directly proportional to the true expression level of the gene and also positively correlated with the length of the gene and the depth of sequencing. The uniplexed Unigene database was used to identify the abundance of each Unigene by sequence similarity alignment. Bowtie2 software was used to get the number of Unigene reads in each sample, and Express software was applied to calculate the Unigene expression FPKM value. The number of Unigene counts for each sample was normalized using DESeq software to calculate the difference multiple. The negative binomial distribution test was used to test the significance of the difference of read numbers. Finally, the differentially expressed Unigenes were screened based on the difference multiple and the difference significance test results. The threshold for screening differences was p value < 0.05 and a fold change > 2.
The differential expression of Unigenes was analyzed by GO enrichment and its function was described. The number of differential mRNAs included in each GO term was counted and the significance of the differential Unigene enrichment in each GO entry was calculated using the hypergeometric distribution test. The pathway analysis of the differential Unigenes was performed based on the KEGG database, and the significance of the differential Unigene enrichment in each pathway entry was calculated by the hypergeometric distribution test.
Quantitative real-time PCR (qRT-PCR) analysis
Trizol reagent (Takara Bio, Otsu, Japan) was used to extract the total RNA from the hepatopancreas samples following the manufacturer’s instructions. Electrophoresis using 1.0% agarose gel helped to assess the quality of RNA. After the purification of total RNA, PrimeScript®RT Reagent Kit With gDNA Eraser (Takara Bio Inc., Shiga, Japan) was used to synthesize first-strand cDNA following the manufacturer’s instructions. ddH2O was used to dilute the cDNA obtained to 1:10, which was used as a template for RT-PCR. Table 1 lists the primer information. SYBR®Premix Ex Taq™ (QIGAN) was used for conducting the qRT-PCR amplifications with a 20-mL reaction mixture on the ABI7500 Real-Time PCR Detection System (Applied Biosystems). All of these experiments were conducted in triplicate. The comparative threshold cycle method (2-△△CT) was used to analyze the expression levels of target genes, taking β-actin as the reference gene. Each experiment involved three biological replicates. All the data underwent tests to meet the requirements for ANOVA, and a one-way ANOVA together with Tukey’s HSD test were carried out to analyze the data. All data were represented by the mean ± standard deviation. SPSS statistical software 22.0 was used for statistical analyses. P values less than 0.05 were considered to be statistically significant.
Availability of data and materials
Raw Illumina sequences were deposited in the National Center for Biotechnology Information (NCBI) and our SRA records will be accessible with the following link after the indicated release date: https://www.ncbi.nlm.nih.gov/sra/PRJNA640880, SRA accession: PRJNA640880; Temporary Submission ID: SUB7645760.
Kyoto Encyclopedia of Genes and Genomes
Low light intensity
High light intensity
- S. paramamosain :
Quantitative real-time PCR
Lysosomal thioesterase PPT2
Pyruvate dehydrogenase E1 component subunit beta
2-amino-3-ketobutyrate coenzyme A ligase
ATP synthase subunit d
Fasola M, Biddau CL. Foraging habits of crab plovers dromas ardeola overwintering on the Kenya coast. Colon Waterbirds. 1996;19:207–13.
Zhen JT, Shao YF. Experimental study on thermal environment in over-wintering greenhouse for Scylla serrate. J Zhejiang Ocean Univ. 2001;20:205–8.
Peirson SN, Halford S, Foster RG. The evolution of irradiance detection: melanopsin and the non-visual opsins. Philos Trans R Soc. 2009;364:2849–65.
Valdimarsson SK, Metcalfe NB. Is the level of aggression and dispersion in territorial fish dependent on light intensity? Anim Behav. 2001;61(6):1143–9.
Kestemont P, Jourdan S, Houbart M, Me’lard C, Paspatis M, Fontaine P, Cuvier A, Kentouri M, Baras E. Size heterogeneity, cannibalism and competition in cultured predatory fish larvae: biotic and abiotic influences. Aquaculture. 2003;227:333–56.
Puvanendran V, Brown JA. Effect of light intensity on the foraging and growth of Atlantic cod larvae: inter population difference? Mar Ecol Prog Se. 2000;167:207–14.
Trippel EA, Neil SRE. Effects of photoperiod and light intensity on growth and activity of juvenile haddock (Melanogrammus aeglefinus). Aquaculture. 2003;217:633–45.
Han D, Xie S, Lei W, Lei W, Zhu XM, Yang YX. Effect of light intensity on growth, survival and skin color of juvenile Chinese longsnout catfish (Leiocassis longirostris Gunther). Aquaculture. 2005;248:299–306.
Booth MA, Warner-Smith RJ, Allan GL, Glencross BD. Effects of dietary astaxanthin source and light manipulation on the skin color of Australian snapper Pagrus auratus (bloch & schneider, 1801). Aquac Res. 2015;35(5):458–64.
Tian HY, Zhang DD, Xu C, Wang F, Liu WB. Effects of light intensity on growth, immune responses, antioxidant capability and disease resistance of juvenile blunt snout bream Megalobrama amblycephala. Fish Shellfish Immunol. 2015;47(2):674–80.
Wang T, Cheng Y, Liu Z, Yan S, Long X. Effects of light intensity on growth, immune response, plasma cortisol and fatty acid composition of juvenile Epinephelus coioides reared in artificial seawater. Aquaculture. 2013;414-415(2):135–9.
Honryo T, Kurata M, Okada T, Ishibashi Y. Effects of night-time light intensity on the survival rate and stress responses in juvenile Pacific bluefin tuna Thunnus orientalis (Temminck and Schlegel). Aquac Res. 2013;44(7):1058–65.
Nalbach HO, Nalbach G, Forzin L. Visual control of eye-stalk orientation in crabs: vertical optokinetics, visual fixation of the horizon, and eye design. J Comp Physiol A. 1989;165(5):577–87.
Bai AS. Photic effects on embryonation and phototactic responses by the larvae of Argulus siamensis. Proc Anim Sci. 1981;90(5):513–7.
Gao XL, Zhang M, Li X, Shi C, Song CB, Liu Y. Effects of led light quality on the growth, metabolism, and energy budgets of Haliotis discus. Aquaculture. 2016;453:31–9.
Wang H, Tang L, Wei H, Mu CK, Wang CL. “Butter crab”: an environment-induced phenotypic variation of Scylla paramamosain with special nutrition and flavor. Aquac Res. 2019;50(2):541–9.
Zhou JM, Li N, Wang H, Wang CL, Mu CK, Shi C, Liu L, Li RH, Ye YF, Song WW. Effects of salinity on growth, nutrient composition, fatty acid composition and energy metabolism of Scylla paramamosain during indoor overwintering. Aquac Res. 2020;51:1834–43.
Fishery Bureau, Ministry of Agriculture, China. China fisheries statistical yearbook 2018. Beijing: Chinese Agriculture Express; 2019.
Ma HY, Li SJ, Feng NN, Ma CY, Wang W, Chen W, Ma LB. First genetic linkage map for the mud crab (Scylla paramamosain) constructed using microsatellite and AFLP markers. Genet Mol Res. 2016;15(2):gmr.15026929.
Shi X, Waiho K, Li XC, Ikhwanuddin M, Miao GD, Lin F, Zhang Y, Li SK, Zheng HP, Liu WH, Aweya JJ, Azmie G, Baylon J, Quinitio ET, Ma HY. Female-specific SNP markers provide insights into a WZ/ZZ sex determination system for mud crab Scylla paramamosain, S tranquebarica and S. serrata with a rapid method for genetic sex identification. BMC Genomics. 2018;19:981.
Wang H, Tang L, Wei H, Lu JK, Mu CK, Wang CL. Transcriptomic analysis of adaptive mechanisms in response to sudden salinity drops in the mud crab, Scylla paramamosain. BMC Genomics. 2018;19:421.
Shi X, Lu JX, Wu QY, Khor W, Jude JA, Hanafiah F, Zhang YL, Li SK, Zheng HP, Lin F, You CH, Mhd I, Ma HY. Comparative analysis of growth performance between female and male mud crab Scylla paramamosain crablets: evidences from a four-month successive growth experiment. Aquaculture. 2019;505:351–62.
Farlora R, Araya-Garay J, Gallardo-Escárate C. Discovery of sex-related genes through high-throughput transcriptome sequencing from the salmon louse Caligus rogercresseyi. Mar Genomics. 2014;15(6):85–93.
Ding J, Liu C, Luo SY, Zhang Y, Zhu J. Transcriptome and physiology analysis identify key metabolic changes in the liver of the large yellow croaker (Larimichthys crocea) in response to acute hypoxia. Ecotoxicol Environ Saf. 2019;189:109957.
Wang YC, Wang BJ, Liu M, Jiang KJ, Wang MQ, Wang L. Comparative transcriptome analysis reveals the potential influencing mechanism of dietary astaxanthin on growth and metabolism in Litopenaeus vannamei. Aquaculture Rep. 2020;16:100259.
Zou J, Chen J, Tang N, Gao YQ, Hong MS, Wei W, Cao HH, Jian W, Li N, Deng W, Li ZG. Transcriptome analysis of aroma volatile metabolism change in tomato (Solanum lycopersicum) fruit under different storage temperatures and 1-mcp treatment. Postharvest Biol Technol. 2018;135:57–67.
Unwin MJ, Poortenaar CW, Rowe DK, Boustead NC, Porter MJR. Seasonal profiles in growth, energy reserves, gonad development, and plasma steroids in age 1+ cultured Chinook salmon (Oncorhynchus tshawytscha) females. N Z J Mar Freshwater Res. 2004;38(1):29–41.
Xie SW, Tian LX, Li YM, Zhou W, Zeng SL, Yang HJ, Liu YJ. Effect of proline supplementation on anti-oxidative capacity, immune response and stress tolerance of juvenile Pacific white shrimp, Litopenaeus vannamei. Aquaculture. 2015;448:105–11.
Xuan RJ, Wu H, Lin CD, Ma DD, Li YJ, Xu TA, Wang L. Oxygen consumption and metabolic responses of freshwater crab Sinopotamon henanense to acute and sub-chronic cadmium exposure. Ecotoxicol Environ Saf. 2013;89:29–35.
Horst M, Azem A, Schatz G, Glick BS. What is the driving force for protein import into mitochondria? Biochim Biophys Acta. 1997;1318(1–2):71–8.
Gonzalvez F, Pariselli F, Dupaigne P, Budihardjo I, Lutter M, Antonsson B, et al. tBid interaction with cardiolipin primarily orchestrates mitochondrial dysfunctions and subsequently activates bax and bak. Cell Death Differ. 2005;12(6):614–26.
Hatefi Y. The mitochondrial electron transport and oxidative phosphorylation system. Annu Rev Biochem. 1985;54:1015–69.
Castanoncervantes O, BosquesTistler T, FernandezRiveraRio L, Fanjul ML, Prieto J. Effect of variation in photoperiod and light intensity on oxygen consumption, lactate concentration and behavior in crayfish Procambarus clarkii and Procambarus digueti. Comp Biochem Physiol A Mol Integr Physiol. 1998;119(1):263.
Armitage KB, Wall TJ. The effects of body size, starvation and temperature acclimation on oxygen consumption of the crayfish Orconectes nais. Comp Biochem Physiol Part A Physiol. 1982;73(1):63–8.
Withyachumnarnkul B, Poolsanguan B, Poolsanguan W. Continuous darkness stimulates body growth of the juvenile giant freshwater prawn, Macrobrachiumrosenbergiide Man. Chronobiol Int. 1990;7(2):93–7.
Villarreal H, Hinojosa, p., Naranjo, J. Effect of temperature and salinity on the oxygen consumption of laboratory produced Penaeus vannamei postlarvae. Comp Biochem Physiol Part A Physiol. 1994;108(2–3):331–6.
Kennedy H, Haack T, Hartill V, Matakovi L, Baumgartner ER, Potter H, et al. Sudden cardiac death due to deficiency of the mitochondrial inorganic pyrophosphatase PPA2. Am J Hum Genet. 2016;99(3):675–82.
Georgakoudi I, Quinn KP. Optical imaging using endogenous contrast to assess metabolic state. Annu Rev Biomed Eng. 2012;14(1):351–67.
Nantapong N, Otofuji A, Miguit CT, Adachi O, Toyama H, Matsushita K. Electron transfer ability from NADH to Menaquinone and from NADPH to oxygen of type II NADH dehydrogenase of Corynebacterium glutamicum. Biosci Biotechnol Biochem. 2005;69(1):149–59.
Petersen MC, Vatner DF, Shulman GI. Regulation of hepatic glucose metabolism in health and disease. Nat Rev Endocrinol. 2017;13(10):572–87.
Roach PJ, Depaoli-Roach AA, Hurley TD, Tagliabracci VS. Glycogen and its metabolism: some new developments and old themes. Biochem J. 2012;441:763–87.
Jeong Hae C, Jung PM, Kook Whan K, Yoon Ha C, Hee PS, Won Gun A, Ung Suk Y, Jaehun C. Molecular mechanism of hypoxia-mediated hepatic gluconeogenesis by transcriptional regulation. FEBS Lett. 2005;579:2795–801.
Cota-Ruiz K, Peregrino-Uriarte AB, Felix-Portillo M, Martínez-Quintana JA, Yepiz-Plascencia G. Expression of fructose 1, 6-bisphosphatase and phosphofructokinase is induced in hepatopancreas of the white shrimp Litopenaeus vannamei by hypoxia. Mar Environ Res. 2015;106:1–9.
Wierenga RK, Kapetaniou EG, Venkatesan R. Triosephosphate isomerase: a highly evolved biocatalyst. Cell Mol Life Sci. 2010;67:3961–82.
Allen AE, Moustafa A, Montsant A, Eckert A, Kroth PG, Bowler C. Evolution and functional diversification of fructose bisphosphate aldolase genes in photosynthetic marine diatoms. Mol Biol Evol. 2012;29:367–79.
Blencke HM, Reif I, Commichau FM, Detsch C. Jrg Stülke. Regulation of citB expression in bacillus subtilis: integration of multiple metabolic signals in the citrate pool and by the general nitrogen regulatory system. Arch Microbiol. 2006;185(2):136–46.
Linn TC, Pettit FH, Hucho F, Reed LJ. Alpha-keto acid dehydrogenase complexes. XI. Comparative studies of regulatory properties of the pyruvate dehydrogenase complexes from kidney, heart, and liver mitochondria. Proc Natl Acad Sci U S A. 1969;64:227–34.
Ngo PTH, Kim JK, Kim H, Jung J, Ahn YJ, Kim JG, Lee BM, Kang HW, Kang LW. Expression, crystallization and preliminary x-ray crystallographic analysis of XometC, a cystathionine γ-lyase-like protein from Xanthomonas oryzae pv. Oryzae. Acta Crystallogr. 2008;64:750–8.
Wheatley AM, McLoughlin JV. The effect of haloperidol, spiperone and dantrolene on the concentrations of creatine phosphate, ATP and lactate in brain and skeletal muscle of the rat during halothane anesthesia. Res Commun Chem Pathol Pharmacol. 1991;73(3):259–68.
Jean-Michel W. Metabolic fuels: regulating fluxes to select mix. J Exp Biol. 2011;214:286–94.
Warude D, Joshi K, Harsulkar A. Polyunsaturated fatty acids: biotechnology. Crit Rev Biotechnol. 2006;26:83–93.
Khadake R, Khonde V, Mhaske V, Ranjekar P, Harsulkar A. Functional and bioinformatic characterisation of sequence variants of Fad3 gene from flax. J Sci Food Agric. 2011;91:2689–96.
Huang LM, Lai CP, Chen LF, Chan MT, Shaw JF. Arabidopsis SFAR4 is a novel GDSL-type esterase involved in fatty acid degradation and glucose tolerance. Bot Stud. 2015;56:33.
Yang JL, Li LD, Huang K, Li ZG, Wang Q, Wang XF, Yang HL, Chen PJ. Analysis and evaluation on nutrients in whole viscera of Ruditapes philippinarum. Chin Fishery Qual Stand. 2014;4(2):26–31.
Aside from funding support, we also thank OE biotech Co., Ltd. (Shanghai, China) for sequencing consultation and support.
This work was supported by the Major Agricultural Technology Collaborative Extension Program of Zhejiang province (2020XTTGSC03), Basic Public Welfare Research Program of Zhejiang Province (LY20D060001), Major Sci & Tech Special Project of Zhejiang Province (2016C02055–8), Ministry of Agriculture of China & China Agriculture Research System (CARS-48), the K. C. Wong Magna Fund in Ningbo University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
The animal subjects used in the present study are crabs, which are invertebrates and are exempt from this requirement. S. paramamosain is not an endangered or protected species. All animal work has been conducted according to the relevant national and international guidelines. No specific permissions are required to work with invertebrates in China. Similarly, no specific permissions are required for the collection of S. paramamosain from sample sites because they were not collected from protected areas of land.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The gene-specific primers used in this study, Table S2. Result statistics of splicing, Table S3. Summary for transcriptome sequencing data generated from nine cDNA libraries, Table S4. Summary of annotated genes in different databases, Table S5. Statistical table of SSRs, Table S6. Sub categorization of the metabolism group, Figure S1. Distribution of FPKM value in each sample.
About this article
Cite this article
Li, N., Zhou, J., Wang, H. et al. Transcriptome analysis of genes and pathways associated with metabolism in Scylla paramamosain under different light intensities during indoor overwintering. BMC Genomics 21, 775 (2020). https://doi.org/10.1186/s12864-020-07190-w