- 1Wuxi Fisheries College, Nanjing Agricultural University, Wuxi, China
- 2Key Laboratory of Integrated Rice-Fish Farming Ecology, Ministry of Agriculture and Rural Affairs, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi, China
Red claw crayfish (Cherax quadricarinatus) is an economic freshwater shrimp with great commercial potential. However, the suitable stocking density of C. quadricarinatus is still unclear in integrated rice–crayfish farming system. Thus, this study aimed to investigate the effects of stocking density on growth performance, mitophagy, endocytosis and metabolism of C. quadricarinatus. The C. quadricarinatus was reared at low density (LD, 35.73 g/m2), middle density (MD, 71.46 g/m2) and high density (HD, 107.19 g/m2) in an integrated rice–crayfish farming system. After 90 days of farming, the growth performance of C. quadricarinatus significantly decreased in the MD and HD groups relative to that in the LD group. The HD treatment caused oxidative stress and lipid peroxidation at the end of the experiment in hepatopancreas. Transcriptome analysis showed that there were 1,531 DEGs (differently expressed genes) between the LD group and HD group, including 1,028 upregulated genes and 503 downregulated genes. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis indicated that the DEGs were significantly enriched in endocytosis and mitophagy pathways. Meanwhile, four lipid metabolism pathways, including biosynthesis of unsaturated fatty acids, fatty acid biosynthesis, glycerolipid metabolism and glycerophospholipid metabolism, exhibited an upregulated tendency in the HD group. In conclusion, our data showed that when the stocking density reached up to 207.15 g/m2 in HD group, the growth performance of C. quadricarinatus was significantly inhibited in this system. Meanwhile, the data indicated that C. quadricarinatus may respond to the stressful condition via activating antioxidant defense system, endocytosis, mitophagy and metabolism-related pathways in hepatopancreas.
Introduction
In aquaculture practice, semi-intensive and intensive aquaculture modes are commonly selected by farmers to obtain higher economic benefits, where the farmed animals maintain a high stocking density (Roy et al., 2022). As one of the critical factors, stocking density influences the water quality, growth, survival, behavior, and health of farmed aquatic animals (Ani et al., 2022; Hossain et al., 2022). Studies on stocking density have been conducted in many aquatic species, such as turbot (Scophthalmus maximus) (Liu B et al., 2019), rainbow trout (Oncorhynchus mykiss) (Zahedi et al., 2019), Pacific white shrimp (Litopenaeus vannamei) (Deng et al., 2019) and Siberian sturgeon (Acipenser baerii) (Aidos et al., 2018). A variety of adverse effects are caused in aquatic animals when they are stocked in a high density. It has been reported that the growth performance and survival rate of aquatic animals are decreased under a high stocking density due to the intensification of interspecific competition and/or the deterioration of water quality (Trenzado et al., 2006; Yuan et al., 2018). High stocking density, as a stressor, may disturb the balance of homeostasis in farmed aquatic animals, which further negatively impacts on immune function and physiological status, even increases the susceptibility to diseases (Lin et al., 2015; Ellison et al., 2018; Liu B et al., 2019; Ellison et al., 2020). An increasing stocking density induces oxidative stress leading to high catalase (CAT) activity and malonaldehyde (MDA) content in channel catfish (Ictalurus punctatus) (Refaey et al., 2022). In addition, some studies demonstrated that high stocking density could decrease digestive enzyme activity (Liu et al., 2018) and change body composition (Aidos et al., 2018). Hence, an appropriate stocking density is vital in aquaculture, which may take into consideration both welfare of animal and economic benefits of farming.
Transcriptomic analysis is considered as a reliable tool and has been widely used to assess the molecular mechanism of abnormal physiological change caused by adverse stimulus in aquatic animals. Under high stocking density, transcriptomic analysis reveals adverse alterations in muscle quality and immune function of grass carp (Ctenopharyngodon idellus) (Zhao et al., 2019). According to transcriptomic analysis, Ellison et al. (2018) demonstrated that rearing density significantly impacts susceptibility of Nile tilapia (Oreochromis niloticus) to the oomycete Saprolegnia parasitica. In large yellow croaker (Larimichthys crocea), transcriptomic analysis showed that short crowding stress can induce an immune response, but long-term high stocking density may suppress the immunity (Sun et al., 2017). It should be noted that transcriptomic analysis can provide more molecular function information, which contributes to understanding the mechanism of rearing density in aquatic animals.
In China, there are multiple aquaculture models including pond farming, lake/reservoir farming, indoor recirculating aquaculture and integrated rice-aquatic animals farming. Among these, integrated rice-aquatic animals farming has expanded rapidly in the last 10 years, and the farming area is 2.56 million hectares and aquatic production reaches up to 3.25 million tons in 2020 (National Fesheries Technology Extension Center, 2021). It has been considered as a sustainable strategy that improves the utilization of land and water resources, provides food for human, and alleviates environmental pollution resulted from agricultural production (Bashir et al., 2020). In the rice-aquatic animals co-culture system, the excrement of aquatics animals as fertilizer can be utilized by rice to meet the growth requirement for nitrogen and phosphorus, while weeds, insects and plankton can be eaten by aquatic animals as food, which reduces the input of commercial diet (Vromant et al., 2001; Tsuruta et al., 2010). It has been reported that rice-aquatic animals co-culture system has higher ecosystem service value compared with rice monoculture (Liu et al., 2020). In farming practice, various aquatic animals, such as common carp (Cyprinus carpio) (Liu et al., 2022), Chinese mitten crab (Eriocheir sinensis) (Jiang et al., 2021), red swamp crayfish (Procambarus clarkia) (Si et al., 2018) and Chinese soft-shelled turtle (Trionyx sinensis) (Wu et al., 2021), have been co-cultured with rice. It is worth noting that existing research majorly focuses on the assessment of ecological and economic values, microbial diversity and soil nutrition in the integrated rice-aquatic animals farming system (Si et al., 2018; Liu et al., 2020), but the suitable stocking density of aquatic animals and the effects of stoking density on physiological function are rarely evaluated.
C. quadricarinatus, also known as red claw crayfish, has considerable potential for commercial culture due to high growth rate and well tolerance to stressful conditions (Naranjo-Páramo et al., 2004). There are two farming models for C. quadricarinatus including pond monoculture and rice-crayfish co-culture in China. The suitable stocking density of C. quadricarinatus in the pond monoculture has been reported, and the effect of high stocking density on growth performance was evaluated (Pinto and Rouse, 1996; Naranjo-Páramo et al., 2021). However, in rice-crayfish co-culture system, the suitable stocking density has not been reported. In addition, it is unclear that how high density affects the physiological, biochemical and molecular variations of C. quadricarinatus. Therefore, in this study, we set different rearing densities in the integrated rice-crayfish system, and then compared the growth performance and biochemical parameters among different stocking densities after 90 days farming. Further, we evaluated the underlying mechanism of stress induced by high stocking density in C. quadricarinatus via transcriptomic analysis.
Materials and methods
Animals, experimental design and sampling
C. quadricarinatus used in the study were purchased from Zhejiang Freshwater Fisheries Research Institute (Huzhou, China). The experiment was carried out at the farm of Freshwater Fisheries Research Center (Jingjiang, China). The integrated rice-crayfish farming system consists of a rice field (360 m2) and a canal refuge (0.8 m in depth, 40 m2) (Supplementary Figure S1). In the system, rice seeding (Nangeng 5,055) was transplanted in the middle of June 2021, and the rice was harvested in early November. The management of rice fields was based on local agricultural practice experience. Basic fertilizer (16% nitrogen, 8% phosphorus and 16% potassium) was applied before rice transplantation, and no fertilizer was used after C. quadricarinatus farming.
In the rice-crayfish system, we reared three densities of C. quadricarinatus: low density (LD, 35.73 g/m2), middle density (MD, 71.46 g/m2) and high density (HD, 107.19 g/m2). Each density included three repetitions. The average initial weight was 14.29 ± 1.05 g/crayfish. The experiment lasted 90 days from 22 July 2021. During the experiment, the crayfish were fed on a commercial feed (crude protein ≥30%, crude fat ≥3%, crude fiber ≤8%, crude ash ≤18%, total phosphorus ≥1%, and calcium 1%–3.5%) once every day. The daily feed ration was adjusted according to crayfish weight (1%–2% of weight). During the trial period, the water quality maintains within a reasonable range, such as 0.19–0.48 mg/L total ammonia nitrogen, 0.01–0.04 mg/L nitrite, 23.4–32.1°C temperature and 3.21–5.08 mg/L dissolved oxygen. The amount of fed diet and mortality of crayfish were recorded.
To assess the growth performance, the body weight was measured every 30 days via randomly catching 20% of the crayfish in each group. After 90 days of farming, 45 individuals in each group were sampled at random for biochemical and transcriptome analyses. In biochemical analysis, the hepatopancreas of five individuals was pooled into a sample, while 15 individuals’ hepatopancreas were mixed into a sample to sequence transcriptome. All samples were stored in liquid nitrogen temporarily, and then stored for a long time under −80°C. The use of crayfish in this study was approved by the Freshwater Fisheries Research Centre (FFRC), Wuxi, China. All experiment operations were performed according to the requirement of animal welfare.
Growth performance
The growth performance was assessed by calculating specific growth rate (SGR, % day−1), weight gain ratio (WGR, %) and survival rate (SR, %) under three different densities. The calculation formulas are presented as follows:
Where, W1 and W2 are initial weight (g) and final weight (g), n is the days of the feeding trial.
Determination of oxidative stress parameters
Oxidative stress parameters including total antioxidant capacity (T-AOC), malondialdehyde (MDA), superoxide dismutase (SOD), catalase (CAT), glutathione (GSH) and glutathione peroxidase (GPx) were determined in hepatopancreas. The determination of all parameters was conducted according to the method described by the manufacturers. In addition, total protein (TP) was detected to calculate the levels of oxidative stress parameters. The commercial kits for TP, T-AOC, SOD and GSH were provided by Nanjing Jiancheng Bioengineering Institute (Nanjing, China). The kits of CAT and MDA were purchased by Beyotime Biotechnology (Shanghai, China). The kit of GPx was ordered from Grace Biotechnology Co., Ltd (Suzhou, China).
RNA extraction, cDNA library construction, and transcriptome sequencing
The hepatopancreas from the LD and HD groups was used to transcriptome analysis. Total RNA of hepatopancreas was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, United States). RNA concentration and integrality was measured using Qubit2.0 Fluorometer (Life Technologies, CA, United States) and Agilent 2,100 Bioanalyzer, respectively (Agilent Technologies, CA, United States). The RNA was used to synthesize cDNA and then construct DNA libraries by PCR amplification according to a standard procedure. The libraries were sequenced using Illumina Novaseq6000 (Gene Denovo Biotechnology Co., Guangzhou, China). The raw data of transcriptome sequences have been submitted to NCBI Sequence Read Archive (SRA) database (NO. PRJNA884003).
To ensure data quality, low-quality reads in the raw data were filtered by fastp (version 0.18.0) (Chen et al., 2018). The value of Q20 (the base quality score ≥20), Q30 (the base quality score ≥30), and GC (GC content in clean reads) of clean reads were counted. The filtered reads were assembled by Trinity (version 2.8.4) (Grabherr et al., 2011), and the integrity of assemble was evaluated by Benchmarking universal Single-Copy Orthologs (BUSCO). Assembled unigenes were annotated via nr (Non-Redundant Protein Sequence Database), SwissProt, KEGG (Kyoto Encyclopedia of Genes and Genomes) and COG/KOG (Clusters of Orthologous Groups of proteins) databases. The unigene expression was calculated and normalized to RPKM (Reads Per Kilobase of transcript per Million mapped reads) by RSEM (Li and Dewey, 2011). Principal component analysis (PCA) was used to evaluate the relationship of samples.
The differently expressed genes (DEGs) was selected using DESeq2 (version 1.20.0) (Love et al., 2014) with the threshold value: false discovery rate (FDR) < 0.05 and fold change (FC) ≥2. All DEGs were annotated to GO (Gene Ontology) and KEGG databases to conduct enrichment analysis. Furthermore, gene set enrichment analysis (GSEA) was used to assess the differences of genes in important KEGG pathways (Subramanian et al., 2005).
Validation of Differently expressed genes by quantitative real-time PCR
To verify the reliability of RNA-seq, the key genes related to endocytosis and mitophagy in hepatopancreas were selected and detected via quantitative real-time PCR (qPCR). Total RNA of hepatopancreas from the LD and HD groups was isolated using TRIzol reagent (Takara Biomedical Technology Co., Ltd, Beijing, China) according to the manufacturer’s instructions. The isolated RNA was used to synthesize cDNA using a commercial kit (Takara, RR047A). The TB Green Premix EX Taq ™ kit (Takara, RR820A) was used to amplify sequence of target gene. The β-actin was selected as a reference gene, and the relative expression of target gene was calculated by the 2−ΔΔCq method (Livak and Schmittgen, 2001). The specific primers are shown in Supplementary Table S1.
Statistical analysis
SPSS (version 25.0) software was used to conduct statistical analysis in this study. All values were presented as mean ± SE (standard error). The data of growth performance and antioxidant parameters were analyzed by one-way ANOVA with LSD post-hoc test, and the relative levels of mRNA between the LD group and the HD group were analyzed by t-test. The normal distribution and homogeneity of the variances was assessed by Shapiro-Wilk test and Levene test, respectively. Pearson test was used for correlation analysis between the qPCR data and RNA-seq data. Differences were considered to be significant if the p-value < 0.05.
Results
Growth performance parameters
The variation of body weight, SR, SGR and WGR of C. quadricarinatus are presented in Figure 1. The average body weight of C. quadricarinatus exhibited an upward trend in different groups as the farming time increased. After 90 days of farming, the average body weight was lower in the MD and HD groups than that in the LD group (p < 0.05). Meanwhile, the SGR and WGR significantly decreased in the MD and HD groups relative to that in the LD group (p < 0.05). However, the SR showed similar change among different groups (p > 0.05).
FIGURE 1. Mean body weight (A), specific growth rate (B), weight gain ratio (C) and survival rate (D) of C. quadricarinatus under different densities in an integrated rice–crayfish farming system. Values are presented as means ± SE. Means with different superscripts denote significant differences (p < 0.05). LD, low stocking density; MD, middle stocking density; HD, high stocking density.
Antioxidant parameters
At the end of the trial, significant differences in antioxidant ability of C. quadricarinatus were observed among different stocking densities (Figure 2). The activities of SOD, CAT and GPx and the contents of MDA and GSH were higher in the HD group than those in the LD group (p < 0.05) after 90 days of farming, whereas the CAT activity was lower in the HD group than that in the LD group (p < 0.05). In addition, the levels of GSH and Gpx were also significantly increased in the MD group compared with the LD group (p < 0.05; Figures 2C,D).
FIGURE 2. Antioxidative parameters in hepatopancreas of C. quadricarinatus under different densities in an integrated rice–shrimp farming system after 90 days. Values are presented as means ± SE (n = 9). (A) SOD; (B) CAT; (C) GSH; (D) GPx, (E) T-AOC; (F) MDA. Means with different superscripts denote significant differences (p < 0.05). LD, low stocking density; MD, middle stocking density; HD, high stocking density.
Transcriptome sequencing and analysis of differently expressed genes
In order to better understand the adverse effects of high stocking density, we analyzed the genes expression profile in hepatopancreas via RNA-seq. After filtering of raw reads, 41,402,944 (99.64%)—47,175,044 (99.66%) clean reads were obtained, which were assembled into unigenes for further analysis. The results of base quality score are as follow: the values of Q20, Q30, and GC were 97.75%–98.05%, 93.58%–94.37% and 44.71%–48.25%, respectively. The mapped ratio was greater than 92.54% (Table 1).
The result of PCA showed that the LD group and HD group were distinctly separated, and samples had higher correlation within group (Figure 3A). In hepatopancreas, a total of 1,531 DEGs were identified between the LD group and HD group, including 1,028 upregulated genes and 503 downregulated genes (Figures 3B,C).
FIGURE 3. Differently expressed genes (DEGs) in the hepatopancreas of C. quadricarinatus between the LD and HD group. (A) The correlation of samples between the LD and HD groups. (B) The number of DEGs in the hepatopancreas between the LD and HD groups. (C) Volcano plot of DEGs between the LD and HD groups.
GO enrichment analysis of differently expressed genes
All DEGs were subjected to GO enrichment analysis to identify significantly enriched biological process, molecular function and cellular component (Figure 4). In the biological process category, the DEGs were enriched mainly in organonitrogen compound catabolic process (p.adjust = 0.003), proteolysis involved in cellular protein catabolic process (p.adjust = 0.032), modification-dependent protein catabolic process (p.adjust = 0.040), modification-dependent macromolecule catabolic process (p.adjust = 0.040) and cellular protein catabolic process (p.adjust = 0.040) (Figure 4B). In the molecular function category, binding and catalytic activity involved the most DEGs (Figure 4A), and lipid transporter activity (p.adjust = 0.090) and cobalt ion transmembrane transporter activity (p.adjust = 0.090) were the top two GO terms (Figure 4C). In the cellular component category, the DEGs mainly participated in cell, cell part and organelle (Figure 4A), and RNA polymerase III transcription factor complex (p.adjust = 0.002) and transcription factor TFIIIC complex (p.adjust = 0.003) were the top two GO terms (Figure 4D).
FIGURE 4. GO enrichment analysis in the hepatopancreas of C. quadricarinatus under two stocking densities (LD-vs-HD). (A) GO enrichment terms of DEGs in three ontologies. (B) The top 8 GO terms in the biological process. (C) Top 5 GO terms in the molecular function. (D) The top five GO terms in the cellular components.
KEGG enrichment analysis of differently expressed genes
To understand biological function and key signaling pathways, the DEGs were enriched in KEGG database (Figure 5). The result indicated that annotated DEGs were mapped to 124 specific KEGG pathways, and metabolism was the most abundant KEGG A class (Figure 5A). The top 10 enriched KEGG pathways were presented in Figure 5B. It was worth nothing that the DEGs principally enriched in immune function (e.g., endocytosis, mitophagy-animal, autophagy-other eukaryotes and phagosome) and metabolism function (e.g., taurine and hypotaurine metabolism and biosynthesis of unsaturated fatty acids). In addition, the FoxO signaling pathway and TGF-bata signaling pathway were also affected by different stocking densities.
FIGURE 5. KEGG enrichment analysis in the hepatopancreas of C. quadricarinatus under two stocking densities (LD-vs-HD). (A) The result of DEGs enrichment in the KEGG. (B) The top 10 enriched KEGG pathways.
Change of mitophagy pathway
At the end of the trial, 13 genes in mitophagy pathway were significantly changed, including 12 upregulated genes and 1 downregulated gene (p < 0.05; Figure 6A). The qPCR data showed that 7 key genes regulated mitophagy pathway, including casein kinase II subunit alpha (csnk2a), sequestosome-1 (p62), Bcl-2 nineteen kilodalton interacting protein 3 (bnip3), autophagy-related gene 9A (atg9), TBC1 domain family member 15-like isoform X2 (tbc1d15), microtubule-associated proteins 1A/1B light chain 3A (lc3a) and FUN14 domain-containing protein 1 (fundc1) were also significantly upregulated in the HD group (p < 0.05; Figure 6B), which was significantly consistent with the RNA-seq data (r = 0.872, p = 0.010; Figure 6C).
FIGURE 6. Changes of mitophagy pathway in the hepatopancreas of C. quadricarinatus between the LD group and the HD group. (A) DEGs enriched in mitophagy pathway. (B) Mitophagy-related gene expression measured by qPCR, values are presented as means ± SE (n = 3) (C) The correlation between the qPCR data and RNA-seq data.
Change of endocytosis pathway
After 90 days of farming, the endocytosis pathway in the hepatopancreas of C. quadricarinatus was significantly changed. KEGG enrichment analysis showed that a total of 32 DEGs were enriched in the endocytosis pathway, including 30 upregulated genes and 2 downregulated genes (Figure 7A). Further, the mRNA levels of 11 key genes involved in the endocytosis pathway, including heat shock protein 70 (hsp70), clathrin light chain (clta), ras-related protein Rab-5C (rab5c), rab7a, rab11a, suppressor protein of bem1/bed5 double mutants (vps4), sorting nexin-12 (snx12), actin-related protein 2/3 complex subunit 4 (arpc4), adaptor protein 2 complex subunit mu (ap-2), dynamin superfamily protein (dnm) and ras-like GTP-binding protein (roh1), were also upregulated in the HD group compared with the LD group (p < 0.05; Figure 7B). Meanwhile, correlation analysis indicated that the RNA-seq data were significantly consistent with the qPCR data (r = 0.852, p = 0.001; Figure 7C).
FIGURE 7. Changes of endocytosis pathway in the hepatopancreas of C. quadricarinatus between the LD group and the HD group. (A) DEGs enriched in endocytosis pathway. (B) Endocytosis-related gene expression measured by qPCR, values are presented as means ± SE (n = 3). (C) The correlation between the qPCR data and RNA-seq data.
Change of lipid metabolism pathways
The result of GSEA indicated that the lipid metabolism function was significantly influenced by stocking density (Figure 8). Under high stocking density, four lipid metabolism-related pathways, including biosynthesis of unsaturated fatty acids (KO01040), fatty acid biosynthesis (KO00061), glycerolipid metabolism (KO00561) and glycerophospholipid metabolism (KO00564) were more likely to be upregulated.
FIGURE 8. Changes in lipid metabolism pathways identified using GSEA in the hepatopancreas of C. quadricarinatus under two stocking densities (LD-vs-HD). (A) Biosynthesis of unsaturated fatty acids. (B) Fatty acid biosynthesis. (C) Glycerolipid metabolism. (D) Glycerophospholipid metabolism. The |Normalized ES| > 1, nominal p-value < 0.05 and FDR <0.25 in each gene set were set as threshold values for statistical significance.
Discussion
Effect of stocking density on growth performance
Growth performance was vital parameter for assessing effects of stocking density on farmed fish or shrimp. Under excessive stocking density, growth performance of aquatic animal is inhibited, which may reduce economic benefits of farming. In this study, C. quadricarinatus in the HD group exhibited a lower growth performance than that in the LD group as evidenced by decreased body weight, SGR and WGR. The growth performance was inhibited when the stocking density reached to 207.15 g/m2 in the rice-crayfish system after 90 days of farming. The similar results have also reported in other shrimp species. For example, the growth rate of white shrimp (Penaeus vannamei) reared at high density (180 shrimp/m2) was significantly lower than that at low density (90 shrimp/m2) at the 30th week (Araneda et al., 2008). The black tiger shrimp (Penaeus monodon) under low density (400 shrimp/m3) had a higher average daily weight gain, SGR, and final biomass than that at middle density (450 shrimp/m3) and high density (500 shrimp/m3), while the FCR (food conversion ratio) was lower at low density after 127 days (AftabUddin et al., 2020). There are many reasons for the poor growth performance at high density, such as the availability of food, water quality, habitable space and physiological status (Riar et al., 2021). A study on P. monodon has demonstrated that negative effects of stocking density on growth performance may be attributed to deterioration of water quality (Nga et al., 2005). Irwin et al. (1999) suggested that the turbot reared at high density had lower growth rate due to the increase of social interactions. In this work, the water quality parameters were similar in the three groups and met the C. quadricarinatus farming standard (Pan et al., 2017) because of the purification of rice fields. Thus, the inhibition of growth performance was more likely due to the interspecific competition. In the integrated rice–crayfish farming system, the shallow water and rice cropping limited the survival space, which may lead to intensification of interspecific competition under high stocking density. A considerable amount of energy was consumed to cope with intensive social interaction, resulting in the reduction of energy used for growth and metabolism activity (Liu et al., 2016; Wang et al., 2022b).
Effect of stocking density on antioxidative status
The antioxidant defense system is crucial to the stress response mechanism (Jiang et al., 2014). It has been proven that the activities of antioxidant enzymes can be regarded as indicators of oxidative stress (Wang et al., 2013). Various antioxidant enzymes and non-enzymatic antioxidants can effectively remove excess reactive oxygen species (ROS) and protect against oxidative damage (Zhang et al., 2022). However, when the stress exceeds the tolerance, the scavenging capacity of the antioxidant defense system is decreased (Winston, 1991; Mathew et al., 2007). In aquatic animals, high stocking density can cause oxidative stress and change antioxidative status, which is related to species and farming conditions (Liu et al., 2017; Dorothy et al., 2021; Campa-Córdova et al., 2022). The antioxidant parameters, such as SOD, CAT, GPx, GSH and T-AOC, were inhibited by high stocking density in largemouth bass (Micropterus salmoides) (Kommaly et al., 2021) and tiger puffer (Takifugu rubripes) (Zhang et al., 2019), which may cause an oxidative damage. However, the study on L. vannamei reared in biofloc system showed that the activities of antioxidant enzymes (SOD, CAT, and GPX) were increased with increasing stocking density (Dorothy et al., 2021). Similarly, the data in olive flounder (Paralichthys olivaceus) indicated that high stocking density significantly increased the expression and activities of SOD and CAT (Choi et al., 2019). In this study, the levels of SOD, GPx, GSH and T-AOC increased in the HD group after 90 days farming, indicating antioxidant defense system was activated to copy with the adverse stimulation caused by high stocking density. In addition, the activity of CAT in the HD group was significantly lower than that in the LD group, indicating chronic oxidative stress may depress the activity of CAT. This was probably because the CAT consumption was greater than its synthesis under chronic oxidative stress caused by high stocking density.
Under oxidative stress, the excess ROS attacks lipids, inducing lipid peroxidation (Yin et al., 2011). MDA is a final product of lipid peroxidation (Wang et al., 2022a). Hence, it is generally recognized that the content of MDA can reflect the level of oxidative damage in organisms (Han et al., 2022). It has been reported that the content of MDA significantly increased with the increasing stocking density in many aquatic animals, such as M. salmoides (Jia et al., 2022), S. maximu (Jia et al., 2016), and blunt snout bream (Megalobrama amblycephala) (Wang et al., 2018). In this study, the significant increase of MDA content in the HD group also proved that long-term high stocking density caused oxidative stress.
Effect of stocking density on mitophagy
The liver is a highly dynamic metabolic organ and a major site of protein synthesis, lipid metabolism and detoxification (Qian et al., 2021). These metabolic processes require a high amount of energy. Mitochondria is the primary energy-generating organelle and its function readily deteriorates under stress condition. To maintain mitochondrial quality, cells can invoke a mitochondria-specific form of degradative process, named mitophagy, to remove damaged and dysfunctional mitochondria (Saito et al., 2021). Mitophagy shares some key regulatory proteins, such as LC3 and p62, with macroautophagy, and it is also regulated by specific proteins including PTEN-induced putative kinase 1 (PINK1), Bnip3, and Fundc1 (Urbina-Varela et al., 2020). The mitophagy is susceptible to external stress, such as hypoxia and nutrient deprivation, which, in turn, significantly influences mitochondrial function, cell survival and energy homeostasis (Ke, 2020). It has been reported that mitophagy is triggered to scavenge impaired mitochondria and reduce ROS level under oxidative stress (Fan et al., 2019a). Thus, mitophagy is considered as a protective response against oxidative damage induced by adverse stimuli (Fan et al., 2019b; Garza-Lombo et al., 2020). In aquatic animals, lots of exogenous stimuli, such as starvation, hypoxia, and bacterial or viral infections, have been proven to cause the activation of autophagy. For example, copper exposure upregulated the mRNA levels of LC3a, PINK1, Parkin, and induced mitophagy in L. crocea (Pan et al., 2020). Dietary methionine deficiency led to induction of mitophagy via PINK1/PARKIN axis in the liver of rainbow trout (Oncorhynchus mykiss) (Seite et al., 2018). In this study, the KEGG analysis showed that the mitophagy pathway was activated under high stocking density. It is worth noting that the key genes including lc3a, p62, bnip3, and fundc1 were significantly upregulated in HD group. We speculated that high stocking density is a stressor that could induce ROS production and further trigger mitophagy. The activation of mitophagy is an adaptive response which may alleviate oxidative damage (Tian et al., 2022). Similar data were also reported in previous study, where mitophagy was activated to maintain cellular homeostasis in loach fin cells under oxidative stress induced by doxycycline exposure (Shan et al., 2022).
Effect of stocking density on endocytosis
Endocytosis is an essential and highly dynamic biological process that is responsible for the internalization of transmembrane receptor ligand complexes, lipids and pathogens (Schroeder and McNiven, 2014). It plays an important role in maintaining cellular homeostasis and interacting with environments (Mellman, 1996), which is mainly manifested in controlling the composition of plasma membrane (Doherty and McMahon, 2009), participating in immune regulation (Lv et al., 2021), absorption of nutrients and other signal transduction of physiological activities (Leborgne-Castel and Luu, 2009). It has been reported that environmental stress-induced adaptive programs can regulate bulk endocytic flux and alter endocytosis (Shin et al., 2021). Increased endocytosis may strengthen cellular resilience via elevating nutrient intake under stress condition (Sébastien, 2021). In addition, the alteration of endocytosis during stress may counteract deleterious effects through switching on or amplifying specific cellular pathways (Lopez-Hernandez et al., 2020). Previous studies also suggested that the significant change of endocytosis pathway could be induced by adverse environmental conditions in aquatic animals, such as nitrite (Yu et al., 2019), salinity (Qin et al., 2020), copper exposure (Xing et al., 2021) and cadmium exposure (Zhang et al., 2021). Vega et al. (2010) suggested the upregulated endocytosis played an important role in the recovery from damage under stress conditions. In the present study, we found that high stocking density induced the upregulation of endocytosis in hepatopancreas of C. quadricarinatus, which probably was a result of oxidative stress. The increased endocytosis may be an adaptive response which enhanced nutrients uptake to maintain the cellular homeostasis.
Effect of stocking density on metabolism
Dynamic change in metabolic function is a major mechanism responded to external stress in aquatic animals. High stocking density as a stressor has been confirmed to induce a change in amino acid carbohydrate and triglyceride metabolism in liver of Patagonian blennie (Eleginops maclovinus) (Oyarzún et al., 2020). Under a high stocking density, glucose metabolic enzymes were activated, which enhanced energy production to resist the environment stimuli in abalone (Haliotis discus hannai) (Gao et al., 2018). A transcriptomic analysis of M. salmoides indicated that high stocking density caused abnormal lipid metabolism (Jia et al., 2022). In line with previous studies, our data also showed that the metabolic function of hepatopancreas, such as organonitrogen compound catabolic process, proteolysis involved in cellular protein catabolic process, modification-dependent protein catabolic process and cellular protein catabolic process, was significantly influenced by high stocking density. We suspected the changes of metabolic function may produce more energy to adapt to the stress condition, meaning that less energy and matter were used for growth under high stocking density.
It has been reported that hepatic lipid metabolism is highly susceptible to adverse stress (Guo et al., 2022). In aquatic animals, high stocking density-induced changes of lipid metabolism has been widely reported. The lipid contents were reduced in liver of gilthead seabream (Sparus aurata) and piabanha (Brycon insignis) under high stocking density, which may reflect a higher lipid utilization to cope with the stress (Montero et al., 2001; Tolussi et al., 2010). Multiple omics analysis showed that high density caused abnormal lipid metabolism in grass carp (Ctenopharyngodon Idella) and lenok (Brachymystax lenok) (Liu Y et al., 2019; He et al., 2021). Similarly, our data exhibited altered lipid metabolism in hepatopancreas of C. quadricarinatus under high stocking density. Meanwhile, the four key pathways associated with lipid metabolism, including fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, glycerolipid metabolism and glycerophospholipid metabolism, displayed an upregulated tendency in the HD group, which may be an adaptive response to stressful condition. The results were consistent with previous studies. For example, activated biosynthesis of unsaturated fatty acids can maintain the membrane fluidity, regulate physiological state and provide energy under starvation and hypoxia stress (Jiao et al., 2020; Ma et al., 2021); Glycerolipid metabolism had an active effect on response to thermal stress in the liver of S. maximus (Zhao et al., 2021); yeast (Saccharomyces cerevisiae) maintained the membrane homeostasis by activating glycerophospholipid metabolism (Xia et al., 2019).
Conculsion
This study showed that C. quadricarinatus is suitable for growing in integrated rice-crayfish farming system according to the growth performance and survival rate. However, the growth performance of C. quadricarinatus was significantly inhibited when the stocking density reached up to 207.15 g/m2 (the stocking density of the HD group on 90th day). Meanwhile, the high stocking density caused oxidative stress after 90 days of farming. In order to cope with the adverse change of physiological state, the endocytosis, autophagy and lipid metabolism pathways were activated in the hepatopancreas of C. quadricarinatus, which may maintain cellular homeostasis, strengthen cellular resilience and provide energy. In addition, the activation of these pathways consumed a considerable amount of energy, resulting in the reduction of energy used for growth activity, which may be a potential mechanism to explain the inhibition of growth under high stocking density. In summary, our study provided a reference for optimizing the stocking density of C. quadricarinatus in an integrated rice-crayfish farming system.
Data availability statement
The datasets presented in this study can be found in online repositories. The name of the repository and accession number can be found below: NCBI; PRJNA884003.
Ethics statement
The animal study was reviewed and approved by Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences.
Author contributions
YD conceptualization, writing—original draft preparation, and software; RJ methodology, investigation validation, formal analysis, and resources; YH and WD data cu-ration, writing—review, and editing; JZ and BL visualization, supervision, project administration, and funding acquisition. All authors read and approved the final version of the manuscript.
Funding
This study was supported by Central Public-interest Scientific Institution Basal Research Fund, CAFS (2020TD60), earmarked fund for CARS (CARS-45–20), and National Key R&D Program of China (2019YFD0900305).
Conflicts of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.1040712/full#supplementary-material
References
AftabUddin S., Siddique M. A. M., Sein A., Dey P. K., Rashed-Un-Nabi M., Haque M. A. (2020). First use of biofloc technology for Penaeus monodon culture in Bangladesh: Effects of stocking density on growth performance of shrimp, water quality and bacterial growth. Aquac. Rep. 18, 100518. doi:10.1016/j.aqrep.2020.100518
Aidos L., Vasconi M., Abbate F., Valente L. M. P., Lanfranchi M., Giancamillo A. D. (2018). Effects of stocking density on reared Siberian sturgeon ( Acipenser baerii ) larval growth, muscle development and fatty acids composition in a recirculating aquaculture system. Aquac. Res. 50 (2), 588–598. doi:10.1111/are.13936
Ani J. S., Manyala J. O., Masese F. O., Fitzsimmons K. (2022). Effect of stocking density on growth performance of monosex Nile Tilapia (Oreochromis niloticus) in the aquaponic system integrated with lettuce (Lactuca sativa). Aquac. Fish. 7 (3), 328–335. doi:10.1016/j.aaf.2021.03.002
Araneda M., Pérez E. P., Gasca-Leyva E. (2008). White shrimp Penaeus vannamei culture in freshwater at three densities: Condition state based on length and weight. Aquaculture 283 (1-4), 13–18. doi:10.1016/j.aquaculture.2008.06.030
Bashir M. A., Liu J., Geng Y., Wang H., Pan J., Zhang D., et al. (2020). Co-culture of rice and aquatic animals: An integrated system to achieve production and environmental sustainability. J. Clean. Prod. 249 (C), 119310. doi:10.1016/j.jclepro.2019.119310
Campa-Córdova Á. I., Angulo C., Zarain-Herzberg M., Pacheco-Marges R., Ascencio F., Guzmán-Murillo M. A., et al. (2022). Stressing stocking density and rearing time effect on whiteleg shrimp (Penaeus vannamei) reared intensively in floating cages. Lat. Am. J. Aquat. Res. 50 (2), 158–167. doi:10.3856/vol50-issue2-fulltext-2814
Chen S., Zhou Y., Chen Y., Gu J. (2018). Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34 (17), i884–i890. doi:10.1093/bioinformatics/bty560
Choi C. Y., Choi J. Y., Choi Y. J., Kim B.-S., Kim J.-W. (2019). Effects of green wavelength light on antioxidant and non-specific immune responses of the olive flounder Paralichthys olivaceus maintained at different stocking densities. Aquac. Eng. 84, 23–28. doi:10.1016/j.aquaeng.2018.11.004
Deng Y., Xu X., Yin X., Lu H., Chen G., Yu J., et al. (2019). Effect of stock density on the microbial community in biofloc water and Pacific white shrimp (Litopenaeus vannamei) gut microbiota. Appl. Microbiol. Biotechnol. 103 (10), 4241–4252. doi:10.1007/s00253-019-09773-4
Doherty G. J., McMahon H. T. (2009). Mechanisms of endocytosis. Annu. Rev. Biochem. 78, 857–902. doi:10.1146/annurev.biochem.78.081307.110540
Dorothy M. S., Harikrishna V., Arun S., Appidi R., Babitha R. A. M. (2021). Growth, body composition and antioxidant status of Litopenaeus vannamei juveniles reared at different stocking densities in the biofloc system using inland saline groundwater. Aquac. Res. 52 (12), 6299–6307. doi:10.1111/are.15493
Ellison A. R., Uren Webster T. M., Rey O., Garcia de Leaniz C., Consuegra S., Orozco-terWengel P., et al. (2018). Transcriptomic response to parasite infection in Nile tilapia (Oreochromis niloticus) depends on rearing density. BMC Genomics 19 (1), 723. doi:10.1186/s12864-018-5098-7
Ellison A. R., Uren Webster T. M., Rodriguez-Barreto D., Leaniz C. G., Consuegra S., Orozco-terWengel P., et al. (2020). Comparative transcriptomics reveal conserved impacts of rearing density on immune response of two important aquaculture species. Fish. Shellfish Immunol. 104, 192–201. doi:10.1016/j.fsi.2020.05.043
Fan P., Xie X. H., Chen C. H., Peng X., Zhang P., Yang C., et al. (2019a). Molecular regulation mechanisms and interactions between reactive oxygen species and mitophagy. DNA Cell Biol. 38 (1), 10–22. doi:10.1089/dna.2018.4348
Fan P., Yu X. Y., Xie X. H., Chen C. H., Zhang P., Yang C., et al. (2019b). Mitophagy is a protective response against oxidative damage in bone marrow mesenchymal stem cells. Life Sci. 229, 36–45. doi:10.1016/j.lfs.2019.05.027
Gao X. L., Li X., Shi C., Wu F. C., Song C. B., Liu Y. (2018). Effects of stocking density on growth, metabolism, and energy budget of Haliotis discus hannai Ino. Aquaculture 483, 84–95. doi:10.1016/j.aquaculture.2017.09.045
Garza-Lombo C., Pappa A., Panayiotidis M. I., Franco R. (2020). Redox homeostasis, oxidative stress and mitophagy. Mitochondrion 51, 105–117. doi:10.1016/j.mito.2020.01.002
Grabherr M. G., Haas B. J., Yassour M., Levin J. Z., Thompson D. A., Amit I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29 (7), 644–652. doi:10.1038/nbt.1883
Guo W., Gao B., Zhang X., Ren Q., Xie D., Liang J., et al. (2022). Distinct responses from triglyceride and cholesterol metabolism in common carp (Cyprinus carpio) upon environmental cadmium exposure. Aquat. Toxicol. 249, 106239. doi:10.1016/j.aquatox.2022.106239
Han K. K., Zunaina E., Hanizasurana H., Che-Badariah A. A., Che-Maraina C. H. (2022). Comparison of catalase, glutathione peroxidase and malondialdehyde levels in tears among diabetic patients with and without diabetic retinopathy. J. Diabetes Metab. Disord. 21 (1), 681–688. doi:10.1007/s40200-022-01030-2
He Y., Yu H., Zhao H., Zhu H., Zhang Q., Wang A., et al. (2021). Transcriptomic analysis to elucidate the effects of high stocking density on grass carp (Ctenopharyngodon idella). BMC Genomics 22 (1), 620. doi:10.1186/s12864-021-07924-4
Hossain M. A., Hossain M. A., Haque M. A., Mondol M. M. R., Harun-Ur-Rashid M., Das S. K. (2022). Determination of suitable stocking density for good aquaculture practice-based carp fattening in ponds under drought-prone areas of Bangladesh. Aquaculture 547, 737485. doi:10.1016/j.aquaculture.2021.737485
Irwin S., Halloran J. O., FitzGerald R. D. (1999). Stocking density, growth and growth variation in juvenile turbot, Scophthalmus maximus (Rafinesque). Aquaculture 178 (1), 77–88. doi:10.1016/S0044-8486(99)00122-2
Jia R., Liu B. L., Han C., Huang B., Lei J,L. (2016). Influence of stocking density on growth performance, antioxidant status, and physiological response of juvenile turbot, Scophthalmus maximu, reared in land-based recirculating aquaculture system. J. World Aquac. Soc. 47 (4), 587–599. doi:10.1111/jwas.12295
Jia R., Wang L., Hou Y., Feng W., Li B., Zhu J. (2022). Effects of stocking density on the growth performance, physiological parameters, redox status and lipid metabolism of Micropterus salmoides in integrated rice-Fish farming systems. Antioxidants (Basel) 11 (7), 1215. doi:10.3390/antiox11071215
Jiang Q., Zhang W., Tan H., Pan D., Yang Y., Ren Q., et al. (2014). Analysis of gene expression changes, caused by exposure to nitrite, in metabolic and antioxidant enzymes in the red claw crayfish, Cherax quadricarinatus. Ecotoxicol. Environ. Saf. 104, 423–428. doi:10.1016/j.ecoenv.2013.10.031
Jiang X., Ma H., Zhao Q. L., Yang J., Xin C. Y., Chen B. (2021). Bacterial communities in paddy soil and ditch sediment under rice-crab co-culture system. Amb. Express 11 (1), 163. doi:10.1186/s13568-021-01323-4
Jiao S., Nie M., Song H., Xu D., You F. (2020). Physiological responses to cold and starvation stresses in the liver of yellow drum (Nibea albiflora) revealed by LC-MS metabolomics. Sci. Total Environ. 715, 136940. doi:10.1016/j.scitotenv.2020.136940
Ke P. Y. (2020). Mitophagy in the pathogenesis of liver diseases. Cells 9 (4), 831. doi:10.3390/cells9040831
Kommaly O., Marina P., Mrkonjić F. M., Gavrilović A., Liang X., Liu L., et al. (2021). High stocking density alters growth performance, blood biochemical profiles, and hepatic antioxidative capacity in gibel carp (Carassius gibelio). Fish. Physiol. Biochem. 47 (2), 203–212. doi:10.1007/s10695-020-00905-6
Leborgne-Castel N., Luu D. T. (2009). Regulation of endocytosis by external stimuli in plant cells. Plant Biosyst. - Int. J. Deal. all Aspects Plant Biol. 143 (3), 630–635. doi:10.1080/11263500903233375
Li B., Dewey C. N. (2011). RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 12 (1), 323. doi:10.1186/1471-2105-12-323
Lin Y. C., Chen J. C., Chen Y. Y., Yeh S. T., Chen L. L., Huang C. L., et al. (2015). Crowding of white shrimp Litopenaeus vananmei depresses their immunity to and resistance against Vibrio alginolyticus and white spot syndrome virus. Fish. Shellfish Immunol. 45 (1), 104–111. doi:10.1016/j.fsi.2015.02.012
Liu B., Fei F., Li X., Wang X., Huang B. (2019). Effects of stocking density on stress response, innate immune parameters, and welfare of turbot (Scophthalmus maximus). Aquac. Int. 27 (6), 1599–1612. doi:10.1007/s10499-019-00413-2
Liu D., Tang R., Xie J., Tian J., Shi R., Zhang K. (2020). Valuation of ecosystem services of rice–fish coculture systems in Ruyuan County, China. Ecosyst. Serv. 41 (C), 101054. doi:10.1016/j.ecoser.2019.101054
Liu G., Ye Z., Liu D., Zhu S. (2018). Influence of stocking density on growth, digestive enzyme activity, immunity, and cortisol levels of subadult half-smoothtongue sole Cynoglossus semilaevis in a recirculating aquaculture system. N. Am. J. Aquac. 80 (3), 286–293. doi:10.1002/naaq.10030
Liu G., Zhu S., Liu D., Guo X., Ye Z. (2017). Effects of stocking density of the white shrimp Litopenaeus vannamei (Boone) on immunities, antioxidant status, and resistance against Vibrio harveyi in a biofloc system. Fish. Shellfish Immunol. 67, 19–26. doi:10.1016/j.fsi.2017.05.038
Liu Q., Hou Z., Wen H., Li J., He F., Wang J., et al. (2016). Effect of stocking density on water quality and (growth, body composition and plasma cortisol content) performance of pen-reared rainbow trout (Oncorhynchus mykiss). J. Ocean. Univ. China 15 (4), 667–675. doi:10.1007/s11802-016-2956-2
Liu Y., Liu H. B., Wu W., Yin J. S., Mou Z. B., Hao F. H. (2019). Effects of stocking density on growth performance and metabolism of juvenile Lenok (Brachymystax lenok). Aquaculture 504, 107–113. doi:10.1016/j.aquaculture.2019.01.058
Liu Y. T., Huang S. Y., Li L. J., Zhao T. Z., Li H. Y., Su Z. F., et al. (2022). Comparative study on biological indexes and meat nutritional value of Cyprinus carpio under Earth pond reared mode and rice field reared mode. Chin. Agric. Sci. Bull. 38 (4), 159–164.
Livak K. J., Schmittgen T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262
Lopez-Hernandez T., Haucke V., Maritzen T. (2020). Endocytosis in the adaptation to cellular stress. Cell Stress 4 (10), 230–247. doi:10.15698/cst2020.10.232
Love M. I., Huber W., Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550. doi:10.1186/s13059-014-0550-8
Lv Z. M., Yue Z. X., Shao Y., Li C. H., Zhao X. L., Guo M. (2021). Bax functions as coelomocyte apoptosis regulator in the sea cucumber Apostichopus japonicus. Dev. Comp. Immunol. 118, 103490. doi:10.1016/j.dci.2019.103490
Ma J. L., Qiang J., Tao Y. F., Bao J. W., Zhu H. J., Li L. G., et al. (2021). Multi-omics analysis reveals the glycolipid metabolism response mechanism in the liver of genetically improved farmed Tilapia (GIFT, Oreochromis niloticus) under hypoxia stress. BMC Genomics 22 (1), 105. doi:10.1186/s12864-021-07410-x
Mathew S., Kumar K. A., Anandan R., Viswanathan Nair P. G., Devadasan K. (2007). Changes in tissue defence system in white spot syndrome virus (WSSV) infected Penaeus monodon. Comp. Biochem. Physiol. C. Toxicol. Pharmacol. 145 (3), 315–320. doi:10.1016/j.cbpc.2007.01.001
Mellman I. (1996). Endocytosis and molecular sorting. Annu. Rev. Cell Dev. Biol. 12, 575–625. doi:10.1146/annurev.cellbio.12.1.575
Montero D., Robaina L. E., Socorro J., Vergara J. M., Tort L., Izquierdo M. S. (2001). Alteration of liver and muscle fatty acid composition in gilthead seabream (Sparus aurata) juveniles held at high stocking density and fed an essential fatty acid deficient diet. Fish. Physiol. Biochem. 24 (1), 63–72. doi:10.1023/a:1011145426543
Naranjo-Páramo J., Hernández-Llamas A., Vargas-Mendieta M., Villarreal-Colmenares H. (2021). Stochastic dynamic model analysis of the effect of stocking density on the monosex production of male redclaw crayfish Cherax quadricarinatus reared in commercial gravel-lined ponds. Aquaculture 535, 736351. doi:10.1016/j.aquaculture.2021.736351
Naranjo-Páramo J., Hernandez-Llamas A., Villarreal H. (2004). Effect of stocking density on growth, survival and yield of juvenile redclaw crayfish Cherax quadricarinatus (Decapoda: Parastacidae) in gravel-lined commercial nursery ponds. Aquaculture 242 (1-4), 197–206. doi:10.1016/j.aquaculture.2004.05.017
National Fesheries Technology ExtensionCenter (2021). Industrial development report on integrated rice-fish farming in China during “13th Five-Year”. National Fesheries Technology Extension Center: Beijing.
Nga B. T., Lürling M., Peeters E. T. H. M., Roijackers R., Scheffer M., Nghia T. T. (2005). Chemical and physical effects of crowding on growth and survival of Penaeus monodon Fabricius post-larvae. Aquaculture 246 (1-4), 455–465. doi:10.1016/j.aquaculture.2005.02.026
Oyarzún R., Paredes R., Saravia J., Morera F. J., Muñoz J. L. P., Ruiz-Jarabo I., et al. (2020). Stocking density affects the growth performance, intermediary metabolism, osmoregulation, and response to stress in Patagonian blennie Eleginops maclovinus. Aquaculture 515 (C), 734565. doi:10.1016/j.aquaculture.2019.734565
Pan X. B., Zhang X. X., Lu Y. P., Zheng P. H., Wang A. L., Xian J. A. (2017). Acute toxicity of ammonia and nitrite on juvenile and subadult redclaw crayfish (Cherax quadricarinatus). J. Biosaf. 26 (4), 316–322. doi:10.3969/j.issn.2095-1787.2017.04.011
Pan Y., Ai C. X., Zeng L., Liu C., Li W. C. (2020). Modulation of copper-induced antioxidant defense, Cu transport, and mitophagy by hypoxia in the large yellow croaker (Larimichthys crocea). Fish. Physiol. Biochem. 46 (3), 997–1010. doi:10.1007/s10695-020-00765-0
Pinto G. F., Rouse D. B. (1996). Growth and survival of the Australian red claw crayfish Cherax quadricarinatus at three densities in earthen ponds. J. World Aquac. Soc. 27 (2), 187–193. doi:10.1111/j.1749-7345.1996.tb00268.x
Qian H., Chao X. J., Williams J., Fulte S., Li T. G., Yang L., et al. (2021). Autophagy in liver diseases: A review. Mol. Asp. Med. 82, 100973. doi:10.1016/j.mam.2021.100973
Qin H., Zhu Z. X., Lin H. R., Xia J. H., Jia Y. X. (2020). Exploring candidate genes in a major QTL region associated with salinity tolerance in the skin of Nile tilapia based on transcriptomic analysis. Aquaculture 526, 735380. doi:10.1016/j.aquaculture.2020.735380
Refaey M. M., Li D., Tian X., Onxayvieng K., Tang R. (2022). Physiological responses of channel catfish (Ictalurus punctatus) reared at different stocking densities in a recirculating aquaculture system. Aquaculture 557, 738329. doi:10.1016/j.aquaculture.2022.738329
Riar G. S., Raushon N. A., Paul S. K. (2021). Effect of stocking density on growth performance and the survival of golden mahseer, Tor putitora (Hamilton) fry. Asian J. Fish. Aquatic Res. 14 (5), 47–54. doi:10.9734/ajfar/2021/v14i530308
Roy S. M., Machavaram R., Moulick S., Mukherjee C. K. (2022). Economic feasibility study of aerators in aquaculture using life cycle costing (LCC) approach. J. Environ. Manage. 302, 114037. doi:10.1016/j.jenvman.2021.114037
Saito T., Hamano K., Sadoshima J. (2021). Molecular mechanisms and clinical implications of multiple forms of mitophagy in the heart. Cardiovasc. Res. 117 (14), 2730–2741. doi:10.1093/cvr/cvaa340
Schroeder B., McNiven M. A. (2014). Importance of endocytic pathways in liver function and disease. Compr. Physiol. 4 (4), 1403–1417. doi:10.1002/cphy.c140001
Sébastien L. (2021). Endocytosis and stress: From mechanisms to cellular physiology. Biol. Cell 113 (11), 439–440. doi:10.1111/boc.202100072
Seite S., Mourier A., Camougrand N., Salin B., Figueiredo-Silva A. C., Fontagne-Dicharry S., et al. (2018). Dietary methionine deficiency affects oxidative status, mitochondrial integrity and mitophagy in the liver of rainbow trout (Oncorhynchus mykiss). Sci. Rep. 8 (1), 10151. doi:10.1038/s41598-018-28559-8
Shan J., Xiaoqian D., Xia L., Yu W., Zhilong Z., Zhihui S., et al. (2022). Oxidative stress, autophagy, and apoptosis induced by doxycycline in loach fin cells in vitro. Sci. Total Environ. 839, 156379. doi:10.1016/j.scitotenv.2022.156379
Shin E. Y., Soung N. K., Schwartz M. A., Kim E. G. (2021). Altered endocytosis in cellular senescence. Ageing Res. Rev. 68, 101332. doi:10.1016/j.arr.2021.101332
Si G. H., Yuan J. F., Xu X. Y., Zhao S. J., Peng C. L., Wu J. S., et al. (2018). Effects of an integrated rice-crayfish farming system on soil organic carbon, enzyme activity, and microbial diversity in waterlogged paddy soil. Acta Ecol. Sin. 38 (1), 29–35. doi:10.1016/j.chnaes.2018.01.005
Subramanian A., Tamayo P., Mootha V. K., Mukherjee S., Ebert B. L., Gillette M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Sun P., Bao P. B., Tang B. J. (2017). Transcriptome analysis and discovery of genes involved in immune pathways in large yellow croaker (Larimichthys crocea) under high stocking density stress. Fish. Shellfish Immunol. 68, 332–340. doi:10.1016/j.fsi.2017.07.013
Tian L., Li N., Li K., Tan Y. Z., Han J., Lin B. C., et al. (2022). Ambient ozone exposure induces ROS related-mitophagy and pyroptosis via NLRP3 inflammasome activation in rat lung cells. Ecotoxicol. Environ. Saf. 240, 113663. doi:10.1016/j.ecoenv.2022.113663
Tolussi C. E., Hilsdorf A. W. S., Caneppele D., Moreira R. G. (2010). The effects of stocking density in physiological parameters and growth of the endangered teleost species piabanha, Brycon insignis (Steindachner, 1877). Aquaculture 310 (1), 221–228. doi:10.1016/j.aquaculture.2010.10.007
Trenzado C. E., Morales A. E., Manuel D. L. H. (2006). Physiological effects of crowding in rainbow trout, Oncorhynchus mykiss, selected for low and high stress responsiveness. Aquaculture 258 (1), 583–593. doi:10.1016/j.aquaculture.2006.03.045
Tsuruta T., Yamaguchi M., Abe S. I., Iguchi K. I. (2010). Effect of fish in rice-fish culture on the rice yield. Fish. Sci. 77 (1), 95–106. doi:10.1007/s12562-010-0299-2
Urbina-Varela R., Castillo N., Videla L. A., Del Campo A., Rodrigo U., Castillo N., et al. (2020). Impact of mitophagy and mitochondrial unfolded protein response as new adaptive mechanisms underlying old pathologies: sarcopenia and non-alcoholic fatty liver disease. Int. J. Mol. Sci. 21 (20), 7704. doi:10.3390/ijms21207704
Vega V. L., Charles W., De Maio A. (2010). A new feature of the stress response: increase in endocytosis mediated by Hsp70. Cell Stress Chaperones 15 (5), 517–527. doi:10.1007/s12192-009-0165-2
Vromant N., Chau N. T. H., Ollevier F. (2001). The effect of rice-seeding rate and fish stocking on the floodwater ecology of the trench of a concurrent, direct-seeded rice-fish system. Hydrobiologia 457, 105–117. doi:10.1023/A:1012218725126
Wang X., Xing H. J., Jiang Y., Wu H. D., Sun G., Xu Q. Y., et al. (2013). Accumulation, histopathological effects and response of biochemical markers in the spleens and head kidneys of common carp exposed to atrazine and chlorpyrifos. Food Chem. Toxicol. 62, 148–158. doi:10.1016/j.fct.2013.08.044
Wang X., Yao Q., Lei X. Y., Zhang D. M., Wang S., Wan J. W., et al. (2022a). Effects of different stocking densities on the growth performance and antioxidant capacity of Chinese mitten crab (Eriocheir sinensis) in rice crab culture system. Aquac. Int. 30 (2), 883–898. doi:10.1007/s10499-021-00829-9
Wang X., Zheng R. X., Zhang D. M., Lei X. Y., Wang S., Wan J. W., et al. (2022b). Effects of different stocking densities on growth performance, nutritional quality and economic benefit of juvenile female Chinese mitten crab (Eriocheir sinensis) in rice-crab culture systems. Aquaculture 553, 738111. doi:10.1016/j.aquaculture.2022.738111
Wang Y. W., Zhu J., Ge X. P., Sun S. M., Su Y. L., Li B., et al. (2018). Effects of stocking density on the growth performance, digestive enzyme activities, antioxidant resistance, and intestinal microflora of blunt snout bream ( Megalobrama amblycephala) juveniles. Aquac. Res. 50 (1), 236–246. doi:10.1111/are.13889
Winston G. W. (1991). Oxidants and antioxidants in aquatic animals. Comp. Biochem. Physiol. C Comp. Pharmacol. Toxicol. 100 (1-2), 173–176. doi:10.1016/0742-8413(91)90148-M
Wu B. L., Chen J., Huang L., Zhang Y., Fang T., He J. X. (2021). Dynamics of soil fertility and microbial community response to stocking density in rice-turtle co-culture. Aquac. Rep. 20, 100765. doi:10.1016/j.aqrep.2021.100765
Xia Z. C., Zhou X. L., Li J. Y., Li L., Ma Y., Wu Y., et al. (2019). Multiple-omics techniques reveal the role of glycerophospholipid metabolic pathway in the response of saccharomyces cerevisiae against hypoxic stress. Front. Microbiol. 10, 1398. doi:10.3389/fmicb.2019.01398
Xing K. F., Liu Y. J., Yan C. C., Zhou Y. Z., Sun Y. Y., Su N. K., et al. (2021). Transcriptome analysis of Neocaridina denticulate sinensis under copper exposure. Gene 764, 145098. doi:10.1016/j.gene.2020.145098
Yin H. Y., Xu L. B., Porter N. A. (2011). Free radical lipid peroxidation: mechanisms and analysis. Chem. Rev. 111 (10), 5944–5972. doi:10.1021/cr200084z
Yu J. L., Ji X. S., Wang X. P., Li T. M., Wang H., Zeng Q. F. (2019). Identification and characterization of differentially expressed genes in hepatopancreas of oriental river prawn Macrobrachium nipponense under nitrite stress. Fish. Shellfish Immunol. 87, 144–154. doi:10.1016/j.fsi.2018.12.075
Yuan Q., Qian J., Ren Y., Zhang T. L., Li Z. J., Liu J. S. (2018). Effects of stocking density and water temperature on survival and growth of the juvenile Chinese mitten crab, Eriocheir sinensis, reared under laboratory conditions. Aquaculture 495, 631–636. doi:10.1016/j.aquaculture.2018.06.029
Zahedi S., Akbarzadeh A., Mehrzad J., Noori A., Harsij M. (2019). Effect of stocking density on growth performance, plasma biochemistry and muscle gene expression in rainbow trout (Oncorhynchus mykiss). Aquaculture 498, 271–278. doi:10.1016/j.aquaculture.2018.07.044
Zhang L., Qu K. M., Zang P., Wang L., Li W. D., Zhu J. X. (2019). Effects of the stocking density on stress response and antioxidant status of Takifugu rubripes in recirculating aquaculture systems. Fish. Mod. 46 (4), 14–23. doi:10.3969/j.issn.1007-9580.2019.04.003
Zhang Y., Chen C. F., Shen W. L., Chen J. Q., Wu X. W., Lin Z. H. (2022). Long-term effects of three compound probiotics on water quality, growth performances, microbiota distributions and resistance to Aeromonas veronii in crucian carp Carassius auratus gibelio. Fish. Shellfish Immunol. 120, 233–241. doi:10.1016/j.fsi.2021.11.036
Zhang Y., Chen C., Shen W., Chen J., Wu X., Lin Z. (2021). Comparative transcriptome analysis reveals the biological mechanism of selective cadmium enrichment in Tegillarca granosa. Aquac. Rep. 21, 100960. doi:10.1016/j.aqrep.2021.100960
Zhao H. H., Soufan O., Xia J. G., Tang R., Li L., Li D. P. (2019). Transcriptome and physiological analysis reveal alterations in muscle metabolisms and immune responses of grass carp (Ctenopharyngodon idellus) cultured at different stocking densities. Aquaculture 503, 186–197. doi:10.1016/j.aquaculture.2019.01.003
Zhao T. T., Ma A. J., Huang Z. H., Liu Z. F., Sun Z. B., Zhu C. Y., et al. (2021). Transcriptome analysis reveals that high temperatures alter modes of lipid metabolism in juvenile turbot (Scophthalmus maximus) liver. Comp. Biochem. Physiol. Part D. Genomics Proteomics 40, 100887. doi:10.1016/j.cbd.2021.100887
Keywords: Cherax quadricarinatus, stocking density, transcriptome analysis, growth performance, oxidative stress, metabolism
Citation: Dong Y, Jia R, Hou Y, Diao W, Li B and Zhu J (2022) Effects of stocking density on the growth performance, mitophagy, endocytosis and metabolism of Cherax quadricarinatus in integrated rice–crayfish farming systems. Front. Physiol. 13:1040712. doi: 10.3389/fphys.2022.1040712
Received: 09 September 2022; Accepted: 16 November 2022;
Published: 28 November 2022.
Edited by:
Youji Wang, Shanghai Ocean University, ChinaCopyright © 2022 Dong, Jia, Hou, Diao, Li and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bing Li, bGliQGZmcmMuY24=; Jian Zhu, emh1akBmZnJjLmNu
†These authors have contributed equally to this work