Next Article in Journal
Identification of Potential Hub Genes Related to Aflatoxin B1, Liver Fibrosis and Hepatocellular Carcinoma via Integrated Bioinformatics Analysis
Next Article in Special Issue
Bisphenol A (BPA) Directly Activates the G Protein-Coupled Estrogen Receptor 1 and Triggers the Metabolic Disruption in the Gonadal Tissue of Apostichopus japonicus
Previous Article in Journal
Spatially and Temporally Confined Response of Gastrointestinal Antibiotic Resistance Gene Levels to Sulfadiazine and Extracellular Antibiotic Resistance Gene Exposure in Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of a High-Fat Diet on Intestinal and Gonadal Metabolism in Female and Male Sea Cucumber Apostichopus japonicus

1
CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
2
Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China
3
CAS Engineering Laboratory for Marine Ranching, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
4
University of Chinese Academy of Sciences, Beijing 100049, China
5
Shandong Province Key Laboratory of Experimental Marine Biology, Qingdao 266071, China
6
Institute of Science and Environment, University of Saint Joseph, Macao SAR 999078, China
7
The Innovation of Seed Design, Chinese Academy of Sciences, Wuhan 430072, China
8
Zhongke Tonghe (Shandong) Ocean Technology Co., Ltd., Dongying 257200, China
*
Author to whom correspondence should be addressed.
Biology 2023, 12(2), 212; https://doi.org/10.3390/biology12020212
Submission received: 19 December 2022 / Revised: 24 January 2023 / Accepted: 27 January 2023 / Published: 29 January 2023
(This article belongs to the Special Issue Biology and Ecology of Sea Cucumbers)

Abstract

:

Simple Summary

The sea cucumber Apostichopus japonicus is an important aquatic invertebrate in the aquaculture industry. It has high nutritional and medicinal value and is popular with consumers, in particular in Asia. The high demand for the species has led to declines in wild populations and A. japonicus is listed as an endangered species by IUCN. Aquaculture of sea cucumber is important to safeguard the conservation of the species in the wild, and its successful reproduction in captivity is a prerequisite for industrial production. It is known that parental nutrient reserves, in particular of lipids, are directly related to the reproductive performance of sea cucumbers. In this study, we investigated the impact of a high-fat diet on the intestinal and gonadal metabolism of male and female A. japonicus. Diet lipid content significantly affected the intestinal microbiome and metabolite profile, and also interfered with gonadal metabolites. Dietary lipid supplementation enhanced the dominance of Proteobacteria in the intestines. Other microbes in the intestines responded to a high-fat diet, which may contribute to maintain intestinal homeostasis. The physiological and metabolic responses of sea cucumbers to a high-fat diet showed sex differences. Interestingly, a correlation between intestinal and gonadal tissue levels of polyunsaturated fatty acid was found. The results of this study indicate that the lipid content in diets can be differentially adjusted for male and female sea cucumbers to improve nutrition and promote reproduction. This data also contributes to a better understanding of the reproductive biology and sex differences of sea cucumbers.

Abstract

Parental nutrient reserves are directly related to reproductive performance in sea cucumbers. This study focused on the lipid requirements of male and female sea cucumbers Apostichopus japonicus during the reproductive stage and analyzed their physiological responses to a high-fat diet (HFD). The intestinal lipid metabolites and microbiome profile changed significantly in animals fed with the HFD, as given by an upregulation of metabolites related to lipid metabolism and an increase in the predominance of Proteobacteria in the microbiome, respectively. The metabolic responses of male and female sea cucumbers to the HFD differed, which in turn could have triggered sex-related differences in the intestinal microbiome. These results suggest that the lipid content in diets can be differentially adjusted for male and female sea cucumbers to improve nutrition and promote reproduction. This data contributes to a better understanding of the reproductive biology and sex differences of sea cucumbers.

1. Introduction

Optimizing breeding technology and reproductive performance is crucial for the large-scale sustainable production of aquatic animals [1]. The nutrient reserves of parental animals have a significant role in determining their reproductive performance, with the quantity and quality of these reserves being associated with fecundity and gamete quality [2]. Thus, the energy source required for embryonic and larval development of most aquatic animals is dependent on the parental nutrient reserves [3]. Nutrient reserves are related to dietary intake and to the proportion of different nutrients in the diet, with essential dietary nutrients, such as lipids, directly affecting gonadal development in aquatic animals [4,5]. Lipids have complex functions in organisms, not only in energy storage but also in cell signaling and as membrane components [6]. Dietary lipids are an important source of nutrition, which can not only conserve dietary protein but also improve feed efficiency [7]. In addition, the amount of essential fatty acids significantly affects the reproductive performance of fish, while linolenic acid can affect the content of polyunsaturated fatty acids (PUFA) related to reproduction [8].
The growth and feeding efficiency of aquatic animals can be improved by appropriately increasing dietary lipids [9,10]. When a high-fat diet (HFD) contains a large amount of unsaturated fatty acids, the feed is easily oxidized during storage [11]. However, excessive intake of lipids can affect metabolic processes and significantly interfere with the intestinal microbial community structure [12,13]. In male mice, HFDs were found to alter the seminal fluid and gut microbiome [14]. Aquatic animals are also affected by HFDs, which might impair intestinal health by damaging the intestinal structure and the microbiota [15,16]. The microbiome is thought to influence metabolism and immune responses, which can affect both energy homeostasis and health [17]. Also, metabolic homeostasis was found to be affected by sex in most terrestrial animals, and sex differences in intestinal physiological metabolism might be an important factor leading to other sex-related differences in animals [18,19]. As an example, the physiological responses of zebrafish to HFDs, including changes in intestinal microbiota and fat deposition, were shown to be sex-specific [20,21,22]. Therefore, dietary lipids are important for both animal reproduction and intestinal health.
Sea cucumber Apostichopus japonicus is an important economic aquaculture product in China, with its production area increasing, and associated aquaculture technologies improving, over recent years [23]. The body wall of A. japonicus contains more proteins than lipids and is an excellent source of proteins, monounsaturated fatty acids, and n-3 polyunsaturated fatty acids [24]. The amount of eicosapentaenoic acid (EPA) in the body wall of juvenile sea cucumbers increases with dietary lipid levels [25], but the effects of the latter on reproduction have not been studied. Therefore, our research aimed to investigate the sex differences in lipid metabolism of sea cucumbers to determine whether males and females have different lipid nutritional requirements. For this, we manipulated the lipid content ingested by male and female sea cucumbers and analyzed the physiological changes in the intestines and gonads. Changes in intestinal microorganisms under the influence of an HFD were analyzed by 16s rDNA sequencing, and lipid metabolic pathways and related metabolites were studied using metabolomics approaches. Overall, the study aimed to investigate if the nutrition and reproduction of sea cucumbers could be improved by adjusting the fat content in the diet differently for males and females.

2. Materials and Methods

2.1. Experimental Design

Sea cucumbers A. japonicus were purchased from Zhongke Tonghe (Shandong) Ocean Technology Co., Ltd. (Dongying, China). Animals were allowed to acclimatize to stock ponds for 1 week before the start of the experiment, after which the wet weight of each individual was recorded. Sixty sea cucumbers were equally divided between two groups. The control group (C group) was fed with purchased compound feed mixed with three times the mass of sea mud. The HFD group (H group) was fed compound feed mixed with linseed oil (17:3 compound feed: linseed oil), and with three times the mass of sea mud. The ponds were continuously aerated, and the sea cucumbers were fed three times a day for 40 days. At the end of the experiment, all sea cucumbers were anesthetized with ice-cold water and their wet weights were recorded. They were then sacrificed, and their gonads, intestines, and other tissues removed and weighed. Subsequently, the gonads were checked under a microscope to determine the sex and reproductive stage, and samples from animals in the stage of gonadal proliferation were frozen in liquid nitrogen and stored at −80 °C for further analysis. Among them, 8 female and 8 male sea cucumbers were selected in H group and C group for subsequent experiments. The female sea cucumber intestines and gonads in the C group were marked as ICF and GCF, and in the H group were marked as IHF and GHF, respectively. The male sea cucumber groups were designated as ICM, GCM, IHM, and GHM, respectively.

2.2. Crude Fat and Fatty Acid Profile in Feed

2.2.1. Detection of Crude Fat

To determine the composition of the feed used in the study, five replicate feed samples for each of the two groups were mixed, and the fat content of the mixed samples determined. Samples were weighed and placed in a tube to which 8 mL hydrochloric acid was added. The tubes were placed in a 70–80 °C water bath until the samples were completely digested. The tubes were then removed, and 10 mL ethanol was added and mixed. Once the tubes had cooled, the solution was transferred into a measuring cylinder with a plug, and 25 mL anhydrous ether was added, followed by oscillation. The solution was then left to stand for a few minutes and the supernatant was removed and put into a constant weight flask to which 5 mL anhydrous ether was added for repeated extraction. Anhydrous ether was recovered from the solution and the supernatant was dried in a water bath, followed by drying at 100 ± 5 °C to a constant weight. The crude fat content in the feed was calculated using Equation (1):
X = m 2 m 1 m × 100
where, m1 is the original weight of the flask, m2 is the weight of the flask after final drying, and m is the weight of the tested feed sample.

2.2.2. Analysis of Fatty Acid Profile

Feed samples from the two groups were weighed, placed in flasks, and successively mixed with 100 mg pyrogallic acid, several zeolites, 2 mL 95% ethanol, and 10 mL hydrochloric acid. Flasks were hydrolyzed for 40 min in a 70–80 °C water bath, and 10 mL 95% ethanol was added to flasks after they had cooled. Next, the hydrolysate was transferred to a separating funnel in which a 50 mL mixture of ether and petroleum ether had been added. The separating funnel was shaken for 5 min and left to stand for 10 min. The ether layer was collected into a 250 mL flask, dried in a water bath, and placed in an oven for 2 h at 100 ± 5 °C. Next, 2 mL 2% NaOH methanol solution was added to the flask for 30 min at 85 °C in a water bath, after which 3 mL 14% boron trifluoride methanol solution was added for another 30 min. Once the flask had cooled, 1 mL normal hexane was added; the flask was shaken for 2 min and left to stand for 1 h. After stratification, 100 μL supernatant was added to 1 mL normal hexane, and the solution was filtered through a 0.45 μm membrane. Finally, samples were analyzed by using a gas chromatograph mass spectrometer (Trace1310 ISQ, Thermo Fisher Scientific, Shanghai, China).

2.3. Intestinal Microbiota Analysis

2.3.1. DNA Extraction, PCR Amplification of 16S rDNA Genes, and Sequencing

DNA was extracted from 32 samples of sea cucumber using a QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). DNA purity was assessed on a 1% agarose gel and its concentration was measured using a NanoDrop 2100 device (Thermo Fisher Scientific). PCR amplification and sequencing of 16S rDNA was performed by Omicsmart (Guangdong, China). 16S rDNA was amplified by PCR using 341f341F(CCTACGGGNGGCWGCAG) and 806R(GGACTACHVGGGTATCTAAT) primers targeting the V4 region of the bacterial 16S rRNA gene. After purification of the PCR products, the sequencing library was constructed and sequenced using the Ion S5TMXL platform.

2.3.2. Bacterial Community Composition, Alpha Diversity Analysis, and Function Prediction

The bacterial composition of the samples was analyzed with Uparse software, and the Effective Tags of all samples were clustered. The stacked bar chart of community composition was generated using the R project ggplot2 package (version 2.2.1). The Shannon Index was calculated in QIIME (version 1.9.1). Alpha index comparisons between groups were performed using Welch’s t-test. KEGG pathway analysis was performed using PICRUSt. Analysis of functional differences between groups was performed using the Welch t-test in R project Vegan package (version 2.5.3). Pearson correlation coefficient was used to analyze correlations between different microbiota and pathways.

2.4. Intestinal and Gonadal Metabolite Analysis

2.4.1. Metabolite Extraction and LC-MS/MS Analysis

Tissue samples were added to tubes containing 1 mL extraction solvent and rotated for 30 s. The samples were then homogenized at 45 Hz for 4 min and treated for 5 min in an ice-water bath, with the procedure repeated three times. Next, the samples were incubated for 1 h at −20 °C and centrifuged for 15 min at 12,000 rpm and 4 °C. The supernatant was removed for UHPLC-QE Orbitrap/MS analysis. LC-MS/MS analysis was performed by Gene Denovo Biotechnology Co. (Guangzhou, China) using a UHPLC system (1290, Agilent Technologies, Santa Clara, CA, USA) with a UPLC HSS T3 column (2.1 mm × 100 mm, 1.8 μm) and Q Exactive Orbitrap MS (Thermo Fisher Scientific). The QE mass spectrometer registered MS/MS spectra on an information-dependent basis (IDA). Acquisition software (Xcalibur 4.0.27, Thermo) continuously evaluated full scan measures of the MS data as collected and recorded the acquisition of MS/MS spectra. MS raw data were converted to mzML format using ProteoWizard, and then processed using the R package XCMS (version 3.2) for retention time alignment, peak detection, and peak matching. After data processing using the MS/MS database, peak annotation was performed using OSI-SMMS (version 1.0, Dalian Chem Data Solution Information Technology Co., Ltd.).

2.4.2. Differential Metabolite and Pathway Enrichment Analysis

Differential metabolites were analyzed by orthogonal partial least squares discriminant analysis. The threshold of variable importance in projection (VIP) was set to 1, and univariate analysis was performed in combination with the p-value of the t-test to analyze the differences in metabolites among the tissue samples. Metabolites with p-values < 0.05 and VIP ≥ 1 were considered to be differentially expressed between the experimental groups. The enrichment pathway analysis of differential metabolites was performed through KEGG databases, and metabolites related to lipid metabolism were screened by searching the KEGG B class. Pearson correlation coefficient was used to analyze correlations between differential metabolites and fatty acids in the tissue samples.

3. Results

3.1. Detection of Crude Fat and Fatty Acid Profile in Feeds

The proportion of crude fat content in the feed of the C group and H group was 2.6% and 6.5%, respectively (Table 1). Given the addition of linseed oil to the HFD, the C18:0 of the diet in the H group was more than five times that of the diet in the C group.

3.2. Analysis of Intestinal Microbiota in Sea Cucumbers

3.2.1. Species Composition

The 16s rDNA sequencing results of the intestinal microbiota of the sea cucumbers showed that the top ten most-abundant phyla were Proteobacteria, Firmicutes, Bacteroidetes, Epsilonbacteraeota, Actinobacteria, Fusobacteria, Planctomycetes, Cyanobacteria, Patescibacteria, and Chloroflexi (Figure 1A). The top ten most abundant genera were Pseudophaeobacter, Staphylococcus, Lutibacter, Arcobacter, Vibrio, Bacillus, Psychromonas, Salinicoccus, Propionigenium, and Oceanobacillus (Figure 1A). Welch’s t-test results showed that the intestinal microbiome in the H group was significantly different from that of the C group, regardless of sex. Proteobacteria were significantly upregulated in the H group (Figure 1B,C), mainly as a result of the increased number of Pseudophaeobacter in this group. Regarding sex differences, they were not significant in the C group, whereas variation in Protobacteria underlined differences between males and females in the H group (Figure 1D). Among these, Pseudophaeobacter was more abundant in HF than in HM, whereas there were fewer Methylobacterium in HF than in HM.

3.2.2. Alpha Diversity and Functional Prediction Analysis

Alpha diversity analysis can reflect the richness and evenness of microorganisms in a group. The Shannon index showed that species diversity was high in both the CM and CF groups, but there was no significant difference between the two groups (Figure 2A). However, there were significant differences in intestinal microbial diversity between the sexes in the H group. The microbial diversity in male and female sea cucumbers also decreased significantly after HFD treatment. PICRUSt, used for functional prediction, revealed 33 enriched pathways in four groups, mainly related to metabolism. The top ten pathways were carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins, metabolism of other amino acids, metabolism of terpenoids and polyketides, xenobiotics biodegradation and metabolism, lipid metabolism, energy metabolism, replication and repair, and protein folding, sorting, and degradation (Figure 2B). The results of the Pearson correlation coefficient analysis showed that there was a significant negative correlation between the Proteobacteria and several pathways (Figure 2C), including transcription, signaling molecules and interaction, immune system, digestive system, and infectious diseases. In addition, Pseudophaeobacter was significantly negatively correlated with the pathways of replication and repair, protein folding, sorting and degradation, translation, transcription, signaling molecules and interaction, immune system, and infectious diseases (Figure 2C).

3.3. Analysis of Intestinal and Gonadal Metabolites in Sea Cucumbers

3.3.1. Differences in Metabolites and Enrichment Pathways

The results of the differential analysis of A. japonicus intestinal and gonadal metabolites are shown in Figure 3A. The HFD had effects on the composition of intestinal and gonadal metabolites in sea cucumbers, with 45 upregulated and 100 downregulated metabolites in IHF versus ICF; 72 upregulated and 65 downregulated metabolites in IHM versus ICM; 23 upregulated and 66 downregulated metabolites in GHF versus GCF; and 36 upregulated and 44 downregulated metabolites in GHM versus GCM.
Sex also affected the intestinal and gonadal metabolites of A. japonicus, with 55 upregulated and 6 downregulated metabolites in ICF versus ICM; 13 upregulated and 12 downregulated metabolites in IHF versus IHM; 56 upregulated and 46 downregulated metabolites in GCF versus GCM; and 50 upregulated and 66 downregulated metabolites in GHF versus GHM.
KEGG enrichment analysis of differential metabolites among each group showed that the top-20 enriched pathways were mainly related to metabolism (Figure 3B). Further analysis showed that most of these were lipid metabolism pathways, including fatty acid metabolism (KO00062, KO00071, KO01212), phospholipid metabolism (KO00564), and PUFA metabolism (KO00590, KO00591, KO00592). In addition, a few metabolites were enriched in organismal systems (KO04723) and human diseases (KO05231).

3.3.2. Metabolites Related to Lipid Metabolism Pathways and Their Correlation with Diet

In fatty acid metabolism, changes were mainly seen in the activity of coenzymes (e.g., acetyl-CoA) involved in fatty acid metabolism, elongation, and decomposition pathways. Metabolites, such as choline, lecithin (PC), and lysophosphatidylcholine (LysoPC), varied significantly in each group, and were significantly enriched in the phospholipid metabolic pathway. PUFAs, such as α-linolenic acid, linoleic acid, arachidonic acid (ARA), and EPA, were also significantly altered (Figure 4A). Pearson correlation analysis (Figure 4B) showed that choline, cholesterol, Cer [d18:1/24:1(15Z)], 2-acetyl-1-alkyl-sn-glycero-3-phosphocholine, PE [18:3(9Z,12Z,15Z)/P-18:0], oleic acid, stearidonic acid, ARA, docosahexaenoic acid, cholesterol sulfate, PS (16:0/16:0), 11 types of PC, and 7 types of LysoPC were significantly correlated with dietary fatty acids. In addition, metabolites in gonads associated with dietary fatty acids included oleic acid, alpha-linolenic acid, cholesterol sulfate, and eight types of PC. Given the proportion of C14:0, C16:0, and C18:0 in the experimental diet, if these metabolites were positively correlated with C14:0 and C16:0, they were negatively correlated with C18:0, and vice versa.

4. Discussion

Dietary lipids impact the gonadal development of aquatic animals and the level of essential fatty acids in their diet is directly related to reproductive performance [4,8]. However, excessive intake of dietary lipids can significantly affect lipid metabolism and intestinal microorganisms [26]. The intestinal microbiota of aquatic animals is fluid and more sensitive to the influence of diet, compared with terrestrial animals, and lipid levels, lipid sources, and PUFA in dietary lipids can all affect the gut microbiota [27]. In the current study, the intestinal microbial community structure of the sea cucumber changed significantly after being fed an HFD. Notoriously, Proteobacteria, the most abundant group in all treatments, were significantly upregulated in both male and female sea cucumbers in the H group. These results are in agreement with previous studies in A. japonicus that showed that Proteobacteria is the main phylum in the intestines [28,29,30]. They also agree with studies in crustaceans where Proteobacteria in intestines were the dominant phylum and their abundance was positively correlated with dietary lipid content in Pacific white shrimp (Litopenaeus vannamei) and swimming crabs (Portunus trituberculatus) [31,32]. Taken together with the results of the current study, this suggests that the dominance of Proteobacteria in the intestines is common across different taxa and enhanced under a high lipid diet. Proteobacteria play an important role in the intestines by producing lipopolysaccharides that can trigger inflammation [33]. There was a significant negative correlation between the abundance of Proteobacteria and various pathways, including transcription, signaling molecules and interaction, immune system, digestive system, and infectious diseases, which might indicate that an HFD could affect intestinal functions in sea cucumber, such as digestive and immune function. Analysis at the genera level showed that the enrichment in Proteobacteria was mainly a consequence of a relative increase in Pseudophaeobacter under an HFD. In marine medaka (Oryzias melastigma), it was found that the stress response, glycerol metabolism, and metabolism of aromatic compounds were positively correlated with the abundance of Pseudophaeobacter [34]. It has also been reported that Pseudophaeobacter is rarely identified as a harmful microorganism and rather could serve as a potential antagonist of pathogenic microorganisms [35]. Accordingly, Pseudophaeobacter abundance was significantly negatively correlated with the pathways of replication and repair, protein folding, sorting, and degradation; translation, transcription, signaling molecules and interaction, immune system, and infectious diseases. This suggests that an increased abundance of Pseudophaeobacter in the intestines may have a protective effect against pathogens in sea cucumbers, although further studies would be needed to confirm this hypothesis. Some metabolites of host microbes can promote healthy intestinal function, and dietary components, such as lipids, can affect the intestinal microbiota and host fermentation products [36]. Thus, by changing the original intestinal microbiome structure, it is possible that an HFD contributes to maintaining intestinal homeostasis. Studies have also shown that the intestinal microbiome in mammals is sex-specific [37,38]. However, sex-based differences in the intestinal microbiome have rarely been reported in aquatic animals [39]. A comparison of the intestinal microbiota between male and female sea cucumbers showed no significant difference in community structure in control groups. However, when sea cucumbers were fed an HFD, the intestinal microbes of both sexes changed significantly but not similarly, such that the relative abundance of Pseudophaeobacter was higher in females than in males and the reverse was true for Methylobacterium. Methylobacterium can use single carbon substances, such as methanol and other methylated compounds, as substrates to provide carbon and energy via the serine cycle [40]. According to the changes observed for the intestinal microbiome, it can be confirmed that an HFD affects the intestinal physiology, that is, the colonization environment of intestinal microbes. It is possible that diet had a different impact in male and female physiological metabolism which, in turn changed intestinal microbes differently.
Dietary lipids have been shown to affect specific intestinal metabolites [36]. The current results corroborated these findings but also showed that an HFD affected gonadal metabolites. These metabolites were mainly involved in metabolic pathways, with a few related to disease pathways. Dietary lipid supplementation led to obvious changes in lipid metabolism, including fatty acid, phospholipid, and PUFA metabolism. To analyze the correlation of dietary lipids with intestinal and gonadal metabolites, metabolites related to lipid metabolism pathways were compared to the fatty acids in diets. The results showed that there were 29 and 11 metabolites related to lipid metabolism in the intestines and gonads respectively, that were significantly correlated with dietary fatty acids. In other words, these metabolites were clearly affected by dietary lipids. Moreover, the correlation between these metabolites and dietary lipids was related to the contents of C14:0, C16:0, and C18:0 in the diet. Choline is a vitamin that is a precursor of PC, which is a component of the cell membrane, and is often used as a feed additive because of this important physiological role [41,42]. Studies showed that dietary choline supplementation reduced fat deposition, whereas choline deficiency led to stunted growth, impaired lipid metabolism, and ruptured cell membranes [43]. When the lipid content in the diet was higher, the antioxidant capacity of aquatic animals decreased, which might be why an HFD reduced the choline content in the intestines and gonads of A. japonicus [44]. In addition, choline has an important role in ARA metabolism, and an HFD reduced not only the choline content, but also the ARA content in the intestines [45]. Another important PFUA, DHA, was also reduced in response to the HFD. According to the correlation analysis with dietary fatty acids, DHA and ARA were significantly negatively correlated with C18:0 in the diet. Although the content of C18:0 in the H group was approximately five times higher than that in the C group, there was a significant negative correlation between intestinal stearidonic acid and C18:0 in the diet. A positive correlation between stearidonic acid, ARA, and dietary DHA was found in blunt snout bream (Megalobrama amblycephala), while synergism between these metabolites was also found in the intestines of sea cucumber [46]. In the H group, the main addition was flaxseed oil, of which alpha-linolenic acid was the main component; unsurprisingly, the HFD increased alpha-linolenic acid levels in tissues in both males and females. A significant correlation between an HFD and alpha-linolenic acid was found in the gonads, which not only confirmed the reliability of dietary lipid effects in this experiment, but also suggests that a diet supplemented with linolenic acid affected gonadal metabolites. In rainbow trout studies, females were able to synthesize EPA and DHA from dietary linolenic acid [47]. In tilapia, increased dietary 18:3n-3 caused a corresponding increase in long-chain n-3 PUFAs in muscle lipids [48]. However, in the current study, only alpha-linolenic acid was increased, and DHA, EPA, and ARA in the intestines were significantly reduced. It was also reported that EPA in the body wall generally increased with dietary lipid levels in juvenile A. japonicus [25]. Therefore, it is possible that PUFAs are transferred from the intestines to other tissues in sea cucumbers [49]. Most metabolites were significantly influenced by diet, but sex also had an important role. In the current study, dietary lipids significantly affected DHA, EPA, and ARA in the intestines of sea cucumber, whereas sex had a significant impact on EPA and DHA in the gonads and EPA in the intestines. Meanwhile, DHA levels in the H group were found to be correlated in the gonads and intestines, which supports the previous hypothesis about the possible transfer of PUFA from intestines to gonads in sea cucumber [50]. In zebrafish, sex differences have been found not only in the intestinal microbiome, but also in fat deposition, pointing to sex differences in physiological metabolic responses to HFD [20,21,22]. In the intestines of sea cucumbers, an HFD affected the sphingolipid metabolic pathway and showed gender-based differences, whereas, in the gonads, triglyceride metabolism and biosynthesis of unsaturated fatty acid pathway were also impacted. Therefore, there were no sex differences in the intestines and gonads of the C group, whereas metabolites showed sex differences in the H group, suggesting that sea cucumbers show sex-based differences in physiological metabolic responses to an HFD.

5. Conclusions

This study aimed to determine whether male and female sea cucumbers A. japonicus have different nutritional requirements for lipids by exploring sex differences in lipid metabolism. It was found that an HFD significantly affected the intestinal microbiome and metabolite profile, and also interfered with gonadal metabolites. Dietary lipid supplementation enhanced the dominance of Proteobacteria in the intestines. Other microbes showed changes in the intestines, which may contribute to maintaining intestinal homeostasis. Physiological and metabolic responses of sea cucumbers to an HFD showed sex differences, which has probably led to differences in the intestinal microbiome of male and female sea cucumbers. However, it is also possible that the changes in physiological metabolism are caused by the direct influence of diet on intestinal microbes. Interestingly, an intestinal-gonadal tissue correlation was found in PUFA levels. These results are of relevance to the reproductive biology of sea cucumbers as they show that an HFD impacts metabolism, the intestinal microbiome and reproductive physiology. The study also demonstrates that an HFD has different impacts in the physiology of male and female sea cucumbers. Overall, the data may contribute to improving the nutrition and reproductive performance of sea cucumbers, highlighting the need to account for sex differences when considering the introduction of dietary changes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology12020212/s1, File S1: 16s rDNA SEQ raw data; Table S1: All Metabolites raw data.

Author Contributions

Conceptualization, L.Z. and H.Y.; methodology, S.Z. and X.R.; validation, S.Z.; formal analysis, S.Z.; investigation, S.Z. and J.X.; resources, J.X.; data curation, S.Z.; writing—original draft preparation, S.Z. and H.Y.; writing—review and editing, S.Z., X.R., L.Z., D.G., H.Y. and J.X.; visualization, S.Z.; supervision, H.Y. and L.Z.; project administration, L.Z. and D.G.; funding acquisition, L.Z. and D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 42061160365; the Macao Science and Technology Development Fund, grant number FDCT 0001/2020/AFJ.

Institutional Review Board Statement

This study was approved by the Experimental Animal Ethics Committee of the Institute of Oceanology, Chinese Academy of Sciences (approval date: 1 November 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Izquierdo, M.S.; Fernandez-Palacios, H.; Tacon, A.G.J. Effect of broodstock nutrition on reproductive performance of fish. Aquaculture 2001, 197, 25–42. [Google Scholar] [CrossRef]
  2. Schreck, C.B.; Contreras-Sanchez, W.; Fitzpatrick, M.S. Effects of stress on fish reproduction, gamete quality, and progeny. In Reproductive Biotechnology in Finfish Aquaculture; Elsevier: Amsterdam, The Netherlands, 2001; pp. 3–24. [Google Scholar]
  3. Huang, X.X.; Yin, Y.Q.; Shi, Z.H.; Li, W.W.; Zhou, H.Q.; Lv, W.Q. Lipid content and fatty acid composition in wild-caught silver pomfret (Pampus argenteus) broodstocks: Effects on gonad development. Aquaculture 2010, 310, 192–199. [Google Scholar] [CrossRef]
  4. Duray, M.; Kohno, H.; Pascual, F. The effect of lipid-enriched broodstock diets on spawning and on egg and larval quality of hatchery-bred rabbitfish (Siganus guttatus). Philipp. Sci. 1994, 31, 42–57. [Google Scholar]
  5. Rowe, D.K.; Thorpe, J.E.; Shanks, A.M. Role of fat stores in the maturation of male Atlantic salmon (Salmo salar) parr. Can. J. Fish. Aquat. Sci. 1991, 48, 405–413. [Google Scholar] [CrossRef]
  6. Sunshine, H.; Iruela-Arispe, M.L. Membrane lipids and cell signaling. Curr. Opin. Lipidol. 2017, 28, 408. [Google Scholar] [CrossRef]
  7. Pohlenz, C.; Gatlin III, D.M. Interrelationships between fish nutrition and health. Aquaculture 2014, 431, 111–117. [Google Scholar] [CrossRef]
  8. Watanabe, T.; Arakawa, T.; Kitajima, C.; Fujita, S. Effect of nutritional quality of broodstock diets on reproduction of red sea bream. Nippon Suisan Gakkaishi 1984, 50, 495–501. [Google Scholar] [CrossRef] [Green Version]
  9. Chou, B.S.; Shiau, S.Y. Optimal dietary lipid level for growth of juvenile hybrid tilapia, Oreochromis niloticus × Oreochromis aureus. Aquaculture 1996, 143, 185–195. [Google Scholar] [CrossRef]
  10. Zhang, Y.F.; Sun, Z.Z.; Wang, A.L.; Ye, C.X.; Zhu, X. Effects of dietary protein and lipid levels on growth, body and plasma biochemical composition and selective gene expression in liver of hybrid snakehead (Channa maculata♀× Channa argus♂) fingerlings. Aquaculture 2017, 468, 1–9. [Google Scholar] [CrossRef]
  11. Ma, Q.; Li, L.Y.; Le, J.Y.; Lu, D.L.; Qiao, F.; Zhang, M.L.; Du, Z.Y.; Li, D.L. Dietary microencapsulated oil improves immune function and intestinal health in Nile tilapia fed with high-fat diet. Aquaculture 2018, 496, 19–29. [Google Scholar] [CrossRef]
  12. He, A.Y.; Ning, L.J.; Chen, L.Q.; Chen, Y.L.; Xing, Q.; Li, J.M.; Qiao, F.; Li, D.L.; Zhang, M.L.; Du, Z.Y. Systemic adaptation of lipid metabolism in response to low- and high-fat diet in Nile tilapia (Oreochromis niloticus). Physiol. Rep. 2015, 3, e12485. [Google Scholar] [CrossRef] [PubMed]
  13. Peng, M.; Xue, J.J.; Hu, Y.; Wen, C.G.; Hu, B.Q.; Jian, S.Q.; Liang, L.F.; Yang, G. Disturbance in the homeostasis of intestinal microbiota by a high-fat diet in the rice field eel (Monopterus albus). Aquaculture 2019, 502, 347–355. [Google Scholar] [CrossRef]
  14. Javurek, A.B.; Spollen, W.G.; Johnson, S.A.; Bivens, N.J.; Bromert, K.H.; Givan, S.A.; Rosenfeld, C.S. Consumption of a high-fat diet alters the seminal fluid and gut microbiomes in male mice. Reprod. Fertil. Dev. 2017, 29, 1602–1612. [Google Scholar] [CrossRef] [PubMed]
  15. Falcinelli, S.; Rodiles, A.; Hatef, A.; Picchietti, S.; Cossignani, L.; Merrifield, D.L.; Unniappan, S.; Carnevali, O. Dietary lipid content reorganizes gut microbiota and probiotic L. rhamnosus attenuates obesity and enhances catabolic hormonal milieu in zebrafish. Sci. Rep. 2017, 7, 5512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Tian, Y.; Wang, H.L.; Yuan, F.H.; Li, N.; Huang, Q.; He, L.; Wang, L.M.; Liu, Z.G. Perilla oil has similar protective effects of fish oil on high-fat diet-induced nonalcoholic fatty liver disease and gut dysbiosis. BioMed Res. Int. 2016, 2016, 9462571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Butt, R.L.; Volkoff, H. Gut microbiota and energy homeostasis in fish. Front. Endocrinol. 2019, 10, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Mauvais-Jarvis, F.; Arnold, A.P.; Reue, K. A guide for the design of pre-clinical studies on sex differences in metabolism. Cell Metab. 2017, 25, 1216–1230. [Google Scholar] [CrossRef] [Green Version]
  19. Hudry, B.; Khadayate, S.; Miguel-Aliaga, I. The sexual identity of adult intestinal stem cells controls organ size and plasticity. Nature 2016, 530, 344–348. [Google Scholar] [CrossRef] [Green Version]
  20. Fowler, L.A.; Powers, A.D.; Williams, M.B.; Davis, J.L.; Barry, R.J.; D’Abramo, L.R.; Watts, S.A. The effects of dietary saturated fat source on weight gain and adiposity are influenced by both sex and total dietary lipid intake in zebrafish. PLoS ONE 2021, 16, e0257914. [Google Scholar] [CrossRef]
  21. Erick, N.B.; Crisantema, H.; Raúl, L.H.; Alejandra, G.G.; Bruno, G.G. Overfeeding a high-fat diet promotes sex-specific alterations on the gut microbiota of the zebrafish (Danio rerio). Zebrafish 2019, 16, 268–279. [Google Scholar]
  22. Meguro, S.; Hosoi, S.; Hasumura, T. High-fat diet impairs cognitive function of zebrafish. Sci. Rep. 2019, 9, 17063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Fisheries Administration Bureau of the Ministry of Agriculture and Rural Affairs. 2022 China Fishery Statistical Yearbook; China Agriculture Press: Beijing, China, 2022. [Google Scholar]
  24. Gao, F.; Xu, Q.; Yang, H.S. Seasonal biochemical changes in composition of body wall tissues of sea cucumber Apostichopus japonicus. Chin. J. Oceanol. Limnol. 2011, 29, 252–260. [Google Scholar] [CrossRef]
  25. Liao, M.L.; Ren, T.J.; Chen, W.; Han, Y.Z.; Liu, C.M.; Jiang, Z.Q.; Wang, F.Q. Effects of dietary lipid level on growth performance, body composition and digestive enzymes activity of juvenile sea cucumber, Apostichopus japonicus. Aquac. Res. 2017, 48, 92–101. [Google Scholar] [CrossRef]
  26. Schoeler, M.; Caesar, R. Dietary lipids, gut microbiota and lipid metabolism. Rev. Endocr. Metab. Disord. 2019, 20, 461–472. [Google Scholar] [CrossRef] [Green Version]
  27. Ringø, E.Z.Z.V.; Zhou, Z.; Vecino, J.L.G.; Wadsworth, S.; Romero, J.; Krogdahl, Å.; Olsen, R.; Dimitroglou, A.; Foey, A.; Davies, S. Effect of dietary components on the gut microbiota of aquatic animals. A never-ending story? Aquac. Nutr. 2016, 22, 219–282. [Google Scholar] [CrossRef] [Green Version]
  28. Gao, F.; Li, F.H.; Tan, J.; Yan, J.P.; Sun, H.L. Bacterial community composition in the gut content and ambient sediment of sea cucumber Apostichopus japonicus revealed by 16S rRNA gene pyrosequencing. PLoS ONE 2014, 9, e100092. [Google Scholar] [CrossRef]
  29. Sha, Y.J.; Liu, M.; Wang, B.J.; Jiang, K.Y.; Sun, G.Q.; Wang, L. Gut bacterial diversity of farmed sea cucumbers Apostichopus japonicus with different growth rates. Microbiology 2016, 85, 109–115. [Google Scholar] [CrossRef]
  30. Wang, L.; Zhao, X.W.; Xu, H.C.; Bao, X.Y.; Liu, X.J.; Chang, Y.Q.; Ding, J. Characterization of the bacterial community in different parts of the gut of sea cucumber (Apostichopus japonicus) and its variation during gut regeneration. Aquac. Res. 2018, 49, 1987–1996. [Google Scholar] [CrossRef]
  31. Zhang, M.L.; Sun, Y.H.; Chen, K.; Yu, N.; Zhou, Z.G.; Chen, L.Q.; Du, Z.Y.; Li, E. Characterization of the intestinal microbiota in Pacific white shrimp, Litopenaeus vannamei, fed diets with different lipid sources. Aquaculture 2014, 434, 449–455. [Google Scholar] [CrossRef]
  32. Sun, P.; Jin, M.; Ding, L.Y.; Lu, Y.; Ma, H.N.; Yuan, Y.; Zhou, Q.C. Dietary lipid levels could improve growth and intestinal microbiota of juvenile swimming crab, Portunus trituberculatus. Aquaculture 2018, 490, 208–216. [Google Scholar] [CrossRef]
  33. Mukhopadhya, I.; Hansen, R.; El-Omar, E.M.; Hold, G.L. IBD—What role do Proteobacteria play? Nat. Rev. Gastroenterol. Hepatol. 2012, 9, 219–230. [Google Scholar] [CrossRef] [PubMed]
  34. Feng, S.B.; Zeng, Y.H.; Cai, Z.H.; Wu, J.J.; Chan, L.L.; Zhu, J.M.; Zhou, J. Polystyrene microplastics alter the intestinal microbiota function and the hepatic metabolism status in marine medaka (Oryzias melastigma). Sci. Total Environ. 2021, 759, 143558. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, K.; Liu, L.; He, Y.Y.; Qu, C.F.; Miao, J.L. Effects of Dietary Supplementation with κ-Selenocarrageenan on the Selenium Accumulation and Intestinal Microbiota of the Sea Cucumbers Apostichopus japonicus. Biol. Trace Elem. Res. 2021, 199, 2753–2763. [Google Scholar] [CrossRef] [PubMed]
  36. Shen, W.; Gaskins, H.R.; McIntosh, M.K. Influence of dietary fat on intestinal microbes, inflammation, barrier function and metabolic outcomes. J. Nutr. Biochem. 2014, 25, 270–280. [Google Scholar] [CrossRef]
  37. Shastri, P.; McCarville, J.; Kalmokoff, M.; Brooks, S.P.J.; Green-Johnson, J.M. Sex differences in gut fermentation and immune parameters in rats fed an oligofructose-supplemented diet. Biol. Sex Differ. 2015, 6, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. de la Cuesta-Zuluaga, J.; Kelley, S.T.; Chen, Y.; Escobar, J.S.; Mueller, N.T.; Ley, R.E.; McDonald, D.; Huang, S.; Swafford, A.D.; Knight, R. Age-and sex-dependent patterns of gut microbial diversity in human adults. Msystems 2019, 4, e00261-19. [Google Scholar] [CrossRef] [Green Version]
  39. Liu, Y.K.; Yao, Y.Y.; Li, H.; Qiao, F.; Wu, J.L.; Du, Z.Y.; Zhang, M.L. Influence of endogenous and exogenous estrogenic endocrine on intestinal microbiota in zebrafish. PLoS ONE 2016, 11, e0163895. [Google Scholar] [CrossRef] [Green Version]
  40. Green, P.N.; Ardley, J.K. Review of the genus Methylobacterium and closely related organisms: A proposal that some Methylobacterium species be reclassified into a new genus, Methylorubrum gen. nov. Int. J. Syst. Evol. Microbiol. 2018, 68, 2727–2748. [Google Scholar] [CrossRef]
  41. Qin, D.G.; Dong, X.H.; Tan, B.P.; Yang, Q.H.; Chi, S.Y.; Liu, H.Y.; Zhang, S. Effects of dietary choline on growth performance, lipid deposition and hepatic lipid transport of grouper (Epinephelus coioides). Aquac. Nutr. 2017, 23, 453–459. [Google Scholar] [CrossRef]
  42. Adjoumani, J.J.Y.; Abasubong, K.P.; Phiri, F.; Xu, C.; Liu, W.; Zhang, D. Effect of dietary betaine and choline association on lipid metabolism in blunt snout bream fed a high-fat diet. Aquac. Nutr. 2019, 25, 1017–1027. [Google Scholar] [CrossRef]
  43. Michael, F.R.; Teshima, S.; Koshio, S.; Ishikawa, M.; Uyan, O.; Ren, T. Effect of two choline sources on the performance of postlarval Marsupenaeus japonicus bate. Aquac. Nutr. 2007, 13, 59–64. [Google Scholar] [CrossRef]
  44. Lu, K.L.; Xu, W.N.; Li, J.Y.; Li, X.F.; Huang, G.Q.; Liu, W.B. Alterations of liver histology and blood biochemistry in blunt snout bream Megalobrama amblycephala fed high-fat diets. Fish. Sci. 2013, 79, 661–671. [Google Scholar] [CrossRef]
  45. Huang, M.X.; Lin, H.X.; Xu, C.; Yu, Q.R.; Wang, X.D.; Qin, J.G.; Chen, L.Q.; Han, F.L.; Li, E. Growth, metabolite, antioxidative capacity, transcriptome, and the metabolome response to dietary choline chloride in Pacific white shrimp Litopenaeus vannamei. Animals 2020, 10, 2246. [Google Scholar] [CrossRef]
  46. Wang, C.C.; Jiang, G.Z.; Cao, X.F.; Dai, Y.J.; Huang, Y.Y.; Li, X.F.; Liu, W.B. Effects of dietary docosahexaenoic acid on growth performance, fatty acid profile and lipogenesis of blunt snout bream (Megalobrama amblycephala). Aquac. Nutr. 2020, 26, 502–515. [Google Scholar] [CrossRef]
  47. Lazzarotto, V.; Corraze, G.; Leprevost, A.; Quillet, E.; Dupont-Nivet, M.; Médale, F. Three-year breeding cycle of rainbow trout (Oncorhynchus mykiss) fed a plant-based diet, totally free of marine resources: Consequences for reproduction, fatty acid composition and progeny survival. PLoS ONE 2015, 10, e0117609. [Google Scholar] [CrossRef] [Green Version]
  48. Karapanagiotidis, I.T.; Bell, M.V.; Little, D.C.; Yakupitiyage, A. Replacement of dietary fish oils by alpha-linolenic acid-rich oils lowers omega 3 content in tilapia flesh. Lipids 2007, 42, 547–559. [Google Scholar] [CrossRef]
  49. Ghazali, A.; Azra, M.N.; Noordin, N.M.; Abol-Munafi, A.B.; Ikhwanuddin, M. Ovarian morphological development and fatty acids profile of mud crab (Scylla olivacea) fed with various diets. Aquaculture 2017, 468, 45–52. [Google Scholar] [CrossRef]
  50. Zhang, S.; Ru, X.; Zhang, L.; Gonçalves, D.; Yang, H. Examination of sex-related differences in intestinal and gonadal lipid metabolism in the sea cucumber Apostichopus japonicus. Aquaculture 2022, 562, 738787. [Google Scholar] [CrossRef]
Figure 1. Species composition analysis of intestinal microbiota in sea cucumbers [(i): phyla; (ii): genus]. (A) The top ten most-abundant microbiota in each experimental group; (B) Microbiota that were significantly different between CF and HF. (C) Microbiota that were significantly different between CM and HM. (D) Microbiota that were significantly different between HM and HF.
Figure 1. Species composition analysis of intestinal microbiota in sea cucumbers [(i): phyla; (ii): genus]. (A) The top ten most-abundant microbiota in each experimental group; (B) Microbiota that were significantly different between CF and HF. (C) Microbiota that were significantly different between CM and HM. (D) Microbiota that were significantly different between HM and HF.
Biology 12 00212 g001
Figure 2. Alpha diversity and function analysis of intestinal microbiota in sea cucumbers. (A) Alpha diversity by Shannon index. (B) The top ten most enriched microbiota pathways. (C) Functional correlation analysis of microbiota that were significantly different between experimental groups [(i): phyla; (ii): genus].
Figure 2. Alpha diversity and function analysis of intestinal microbiota in sea cucumbers. (A) Alpha diversity by Shannon index. (B) The top ten most enriched microbiota pathways. (C) Functional correlation analysis of microbiota that were significantly different between experimental groups [(i): phyla; (ii): genus].
Biology 12 00212 g002
Figure 3. Differential intestinal and gonadal metabolites and enrichment pathways in sea cucumbers. (A) Comparison of differential metabolites between experimental groups. (B) The top-20 most-enriched differential metabolite pathways between the experimental groups.
Figure 3. Differential intestinal and gonadal metabolites and enrichment pathways in sea cucumbers. (A) Comparison of differential metabolites between experimental groups. (B) The top-20 most-enriched differential metabolite pathways between the experimental groups.
Biology 12 00212 g003
Figure 4. Intestinal and gonadal metabolites related to lipid metabolism pathways in sea cucumbers and their correlation with fatty acids in experimental diets. (A) Heatmap of metabolites related to lipid metabolism pathways. (B) Lipid metabolites significantly associated with fatty acids in experimental diets [(i): intestines; (ii): gonads].
Figure 4. Intestinal and gonadal metabolites related to lipid metabolism pathways in sea cucumbers and their correlation with fatty acids in experimental diets. (A) Heatmap of metabolites related to lipid metabolism pathways. (B) Lipid metabolites significantly associated with fatty acids in experimental diets [(i): intestines; (ii): gonads].
Biology 12 00212 g004
Table 1. Levels of crude fat and fatty acids in experimental diets.
Table 1. Levels of crude fat and fatty acids in experimental diets.
Content (g/100 g)Experimental Group
CH
Crude fat2.66.5
C14:00.0095N/D
C16:00.09420.0816
C18:00.01030.0583
N/D (not detected): fatty acids were not detected when the content was <0.0033 (g/100 g).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, S.; Ru, X.; Zhang, L.; Gonçalves, D.; Yang, H.; Xu, J. Effects of a High-Fat Diet on Intestinal and Gonadal Metabolism in Female and Male Sea Cucumber Apostichopus japonicus. Biology 2023, 12, 212. https://doi.org/10.3390/biology12020212

AMA Style

Zhang S, Ru X, Zhang L, Gonçalves D, Yang H, Xu J. Effects of a High-Fat Diet on Intestinal and Gonadal Metabolism in Female and Male Sea Cucumber Apostichopus japonicus. Biology. 2023; 12(2):212. https://doi.org/10.3390/biology12020212

Chicago/Turabian Style

Zhang, Shuangyan, Xiaoshang Ru, Libin Zhang, David Gonçalves, Hongsheng Yang, and Jialei Xu. 2023. "Effects of a High-Fat Diet on Intestinal and Gonadal Metabolism in Female and Male Sea Cucumber Apostichopus japonicus" Biology 12, no. 2: 212. https://doi.org/10.3390/biology12020212

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop