Research

Exploring single cell microbial protein as a sustainable fishmeal alternative in yellowtail kingfish (Seriola lalandi) diets: impacts on health and gut microbiome

1,2 ,1,2 ,1 ,1,2,3 ,4,5 ,4 ,6

Abstract

Background

With the global expansion of aquaculture and the increasing demand for fish meal, identifying appropriate and sustainable alternative protein sources for aquafeeds has become essential. Single-cell protein (SCP), derived from methanotrophic bacteria, presents a promising alternative by converting methane into protein, potentially addressing both the need for alternative protein sources and reducing industrial greenhouse gas emissions. This study aimed to evaluate the effects of different levels of SCP inclusion (0%, 25%, 50%, and 75% fish meal replacement) on the health, gene expression, and gut microbiome of yellowtail kingfish (YTK,Seriola lalandi) following a 35-day growth trial.

Results

The study found that SCP inclusion at the highest level of fishmeal replacement (75%) induced a mild inflammatory response in the hindgut of the fish. However, micromorphological assessments of the hindgut, serum biochemistry, and gene expression analyses revealed no significant detrimental effects from SCP replacement. Notably, there were indications of improved lipid digestibility with SCP. Furthermore, SCP inclusion significantly enhanced microbial richness and altered the composition of the gut microbiome, introducing beneficial bacterial taxa that may contribute to improved gut health and resilience.

Conclusions

This study highlights SCP as a viable and sustainable alternative to fish meal in YTK diets. The findings suggest that SCP can be included in YTK diets without adverse health effects at moderate levels and may even offer benefits in terms of lipid digestibility and gut microbiome diversity. These results contribute to the advancement of more sustainable aquaculture practices.

Introduction

Seriola lalandi, commonly known as yellowtail kingfish (YTK), is a temperate, carnivorous, marine species, that is popular in the global sashimi market [1]. It is well-suited to aquaculture as it commands a high market price and is fast-growing [2]. YTK culture is an established industry in Australia and Japan [3, 4] and emerging in New Zealand, South Africa, Europe, and the Americas [5, 6]. Currently, commercial YTK diets require large amounts of fishmeal to optimise the health and growth of the fish, with most commercial diets for the species incorporating approximately 40% fishmeal to meet nutritional needs [7].

The rapid growth of the aquaculture industry and an increasing emphasis on environmental sustainability has encouraged researchers to investigate alternative sources of protein that can completely or partially replace fishmeal from the diets of cultured fish [8,9,10,11,12,13,14]. Previously studied alternative sources of protein include plants and legumes, animal by-products, algae, insects and single-cell proteins (SCP) [10, 15,16,17,18]. Of these, SCP derived from methanotrophic bacteria are of particular interest due to their ability to convert methane into protein. Methanotrophic bacteria utilize methane as a carbon and energy source and ammonia as a nitrogen source [19]. With increasing global concerns around industrial greenhouse gas emissions, the ability to turn waste methane into a product has both economic and environmental benefits [20].

Recent investigations have highlighted the effectiveness of methanotrophic bacterial meal, a promising alternative to traditional fishmeal, in the diets of various aquaculture species. Studies have shown positive outcomes for growth and digestibility in species such as Pacific white shrimp (Litopenaeus vannamei), Atlantic halibut (Hippoglossus hippoglossus), Atlantic salmon (Salmo salar), red sea bream (Pagrus major), rainbow trout (Oncorhynchus mykiss), and Japanese yellowtail (Seriola quinqueradiata) [9, 13, 14, 21,22,23,24,25]. In a previous study [7], we explored SCP inclusion at four dietary levels in YTK diets of 10%, 20% and 30% equating to 0%, 25%, 50%, and 75% fish meal replacement from a control diet containing 40% fish meal. In this 35-day trial, growth rates and feed conversion ratio (FCR) of fish were significantly improved in the diets that replaced 25% of fishmeal with SCP compared to the control. A consistent FCR, protein retention ability, and digestibility coefficients were obtained even with higher SCP replacement levels, despite reduced feed consumption and subsequently growth rates, indicating a potential for higher SCP replacements in YTK diets if palatability issues can be overcome.

While previous studies on methanotrophic bacterial meals have focused primarily on growth metrics, there has been limited exploration into the broader health impacts of dietary SCP inclusion, particularly concerning gene expression and microbiome composition in YTK. Emerging research suggests that dietary interventions can significantly influence gut microbiota, which in turn affects nutrient absorption, immune function, and overall fish health [26, 27]. Additionally, alterations in gene expression related to immune responses, metabolism, and stress have been observed in fish subjected to varying dietary compositions, including those containing alternative protein sources [28]. The YTK microbiome, which plays a crucial role in digestion and pathogen resistance, may be particularly sensitive to changes in diet composition, such as the replacement of fishmeal with SCP [29, 30].

This study aimed to examine the effects of SCP dietary replacement levels on a range of health and immune parameters in YTK, including changes in gene expression and microbiome composition. We hypothesised that there would be no adverse health impacts in fish fed for 35 d with SCP replacing 25%, 50%, and 75% of fishmeal, compared to a control diet containing 40% fishmeal.

Methods

Experimental design

This study used samples from fish in the palatability and growth experiment reported in Pilmer er al. [7]. Detailed experimental protocols are described therein. Briefly, 312 juvenile YTK, averaging 169.3 ± 0.6 g in weight, were individually tagged with RFID pit tags, then evenly and randomly distributed among 24 tanks, comprising 8 treatments (four levels of SCP inclusion, either with or without dietary additives, such as liquid tuna hydrolysate and garlic powder) with 3 replicate tanks per treatment. The SCP used was derived from methanotrophic bacteria (StringBio Pty Ltd., Bengaluru, India) provided by Ridley AgriProducts Pty Ltd. (Melbourne, Australia). All diets were formulated to be isonitrogenous (50% crude protein) and isolipidic (19% crude fat). The four levels of SCP inclusions were 0%, 10%, 20%, and 30% of the diet, which corresponded to 0%, 25%, 50%, and 75% fish meal replacement, respectively. The formulated diets were processed into 5 mm pellets using a commercial extrusion system with vacuum oil coating. During the 35-day growth trial, the fish were fed the experimental diets twice daily until satiety, the tanks were supplied with oxygenated seawater at a flow rate of 6 L/min, with ambient temperature maintained at 20 ± 1 °C. At the end of the trial, individual body weight was recorded for all fish and three fish per tank were then sampled for health measurements.

Morphological and biochemical parameters

Whole blood was collected from the caudal vein using EDTA (2% solution; Sigma) coated needles and stored in 1-mL foetal tubes (K3E EDTA Mini collect®, Greiner Bio-One, Austria). Further samples of whole blood were collected without EDTA from the same fish and stored in 2-mL Eppendorf tubes for serum collection. The latter blood samples were allowed to clot at room temperature (~ 20 °C) for 1 h and then stored at 4 °C overnight before being spun for 10 min at 3,000 × g at 4 °C. Blood serum was removed, divided into triplicate 200 μL aliquots and frozen before analyses of biochemical parameters and serum antioxidant capacity. The biochemical parameters of cholesterol, alanine aminotransferase (ALT), glutamate dehydrogenase (GLDH), lipase, urea, total protein and triglyceride were measured on an ‘AU 680’ analyser. Fast, slow and total antioxidant capacity were measured using an e-BQC lab (BioQuoChem, Spain).

Samples of the liver, hindgut and pyloric caeca were dissected and fixed in 10% neutral buffered formalin for traditional histological examination. Histology sections were stained in H&E and PAS-Alcian Blue. Hindgut sections were used for the determination of lamina propria area as a percentage of villus area (% lamina propria) and muscular thickness. Further sections of the hindgut were dissected and immediately frozen in liquid nitrogen for analysis of myeloperoxidase (MPO) activity as a marker of neutrophil accumulation in the hindgut. The method used for MPO analysis was the same as described in Woolley et al. [31]. Briefly, frozen tissue was ground using a pre-chilled mortar and pestle and homogenised in a 0.5% hexadecyltrimethylammonium bromide (HTAB) solution prepared in 50 mmol/L phosphate-buffered saline (PBS). The homogenized samples were centrifuged, and the supernatant was collected. To determine MPO activity, a fluorescence-based assay was conducted: samples were diluted, and triplicates were loaded into a black 384-well microplate along with a working solution of 10 µmol/L 2-[6-(4-aminophenoxy)-3-oxo-3H-xanthen-9-yl] benzoic acid (APF) and 10 µmol/L hydrogen peroxide in PBS. Fluorescence was measured every minute for 30 min at an excitation wavelength of 485 nm and emission between 515–530 nm.

Gene expression

Twelve genes involved in immune function and growth were targeted for analysis of expression in different tissues (Table 1). Total RNA was isolated from 30 mg samples of intestinal, liver, and brain tissues, all preserved in RNAlater. Samples from all tissues, were only analysed from fish in the SCP0%, SCP25%, and SCP75% treatment groups due to the lack of significant differences observed in other health parameters. These tissues were processed using a TissueRuptor (Qiagen, Melbourne, Victoria, Australia) in 600 µL of RLT buffer for thorough disruption and homogenization, followed by RNA purification with the RNeasy Mini Kit (Qiagen, Melbourne, Victoria, Australia), according to the manufacturer’s guidelines. The concentration of the extracted RNA was quantified using NanoDrop spectrophotometry (ND-2000), with detailed concentration data presented in Table S1. Subsequently, complementary DNA (cDNA) was synthesized from 2 µg of the total RNA via a reverse transcriptase reaction facilitated by the Tetro cDNA Synthesis Kit (Bioline, Australia), which includes a crucial step for the elimination of any potential genomic DNA contamination. Primers and probes specific to the genes under investigation (Table 1) were designed based on species-specific mRNA sequences utilizing the Roche Assay Design platform, as outlined in Table S2.

Table 1 Genes were analysed in the different tissues sampled for gene expression

Quantitative PCR (qPCR) assays were executed in duplicate on a Rotor-gene 6000 Thermocycler (Qiagen, Germany), employing 20 µL reaction mixtures that comprised the FastStart Universal Probe Master Mix (Rox; Roche Diagnostics GmbH), Universal Probe Library Probe (Roche), primers, molecular biology-grade water, and synthesised cDNA. Each assay included a non-template control for background amplification. Cycling conditions involved an initial uracil N-glycosylase pre-treatment, followed by a denaturation step and 35 amplification cycles, each with specific temperature and time settings.

For normalisation purposes, the stability of potential reference genes, including those coding for 18S rRNA, β-actin, and elongation factor 1 alpha, was evaluated, with the 18S rRNA gene being selected due to its consistent expression levels across samples. Data were normalised, log2 transformed and analysed using REST 2008 2.0.7 software and the GeneX tool for efficiency-corrected expression ratio calculations.

Microbiome

Microbial DNA was isolated from the hindgut of experimental fish using the QIAamp 96 DNA QIACube HT kit (Cat. #51331, Qiagen, Germany), as per the manufacturer’s protocol. Intestinal tissue was removed before DNA extraction to minimize eukaryotic contamination. Whole tissue samples were preserved in RNAlater immediately after dissection. For microbial DNA extraction, 20 mg was digested with proteinase at 56 °C for an 1 h, followed by DNA extraction. The quantity of extracted DNA was determined through the Qubit dsDNA High Sensitivity assay (Cat. #Q32851, ThermoFisher Scientific) and then preserved at −20 °C until further processing. Prior to library construction, DNA samples were standardised to 10 ng/μL and re-measured using the Qubit dsDNA HS assay.

The V3–V4 region of the 16S rRNA gene, a highly conserved genetic marker widely used for bacterial identification and community profiling, was selectively amplified from the purified DNA using the Quick-16S Plus NGS Library Prep kit. The 16S rRNA gene was chosen because it allows for efficient characterisation of bacterial diversity and composition in complex microbiomes, such as the fish gut, making it an appropriate target for understanding microbial changes in response to dietary SCP. Amplification employed a single-step PCR with specific V3–V4 16S rRNA and Illumina i7 and i5 index primers, under defined thermocycling conditions. Post-amplification, products were electrophoretically assessed and, if necessary, subjected to additional amplification cycles. Purification of the correctly sized products (~ 600 bp) was achieved using MagBeads, with subsequent quantification and pooling following the kit’s guidelines. Paired-end sequencing was executed on the MiSeq platform at Genomics WA, Australia.

The ensuing Illumina reads underwent initial processing in qiime2, encompassing quality trimming, denoising, chimera filtering, and amplicon sequence variant (ASV) picking via the q2-dada2 plugin [32, 33]. This step included specific trimming parameters to eliminate low-quality sequences. Multiple sequence alignment and phylogenetic tree generation were conducted using Mafft and Fasttree, respectively [34, 35]. ASVs were classified at various taxonomic levels using a pre-trained naive Bayes classifier and the consensus blast method against the SILVA 138 database [36]. Taxonomy-based filtering was applied to remove sequences related to mitochondria and chloroplasts. For diversity and community composition analyses, samples were normalized to the lowest sequencing depth, and samples with fewer than 1,000 high-quality reads were excluded.

Data analysis

The results were subjected to a comprehensive statistical analysis to ascertain the differences between treatments across various parameters, including growth immune and biochemical parameters, serum antioxidant capacity, and histomorphological measurements. Two-way ANOVA was used to examine the impacts of SCP levels and feed additives on immune and biochemical parameters, serum antioxidant capacity, and histomorphological measurements, using tank mean values (n = 3) as the unit of observation. When significant differences were identified, a post hoc Tukey’s test was performed to determine where these differences occurred. Proportional data were arcsine transformed and residual plots were examined to verify the assumption of normality in all analyses to account for the risk of type I errors due to multiple comparisons, the Benjamini–Hochberg procedure was employed, with significance threshold set at P < 0.05. All statistical analyses were performed using JMP (Version 15.0.0).

Microbial diversity analyses were conducted using the phyloseq, ape, and vegan (adonis) packages in R [37, 38]. Diversity within fish samples (alpha diversity) was estimated by ASV richness (number of ASVs) and the Shannon-Weiner diversity index that considers both ASV richness and evenness of abundances. Diversity indices were compared between treatment groups using a non-parametric two-way ANOVA (Scheirer–Ray–Hare test). Differences in microbial community composition between fish samples (beta diversity) were assessed using unweighted and weighted UniFrac distance metrics. The effect of SCP replacement level and feed additives on beta diversity was tested by a permutational analysis of variance (PERMANOVA, 1,000 permutations) applied to the distance matrices. Venn diagrams, generated with the MicroEco R package [39], illustrated the shared and unique ASVs among treatments, while taxa abundances were analysed and visualized using the phyloseq R package, focusing on taxa with over 1% read abundance. Differential abundances of higher taxa (genus and phylum-levels) among groups were identified using the MicrobiomeMarker R package [40], with ambiguous taxa categorized as “unclassified” for clarity in composition analysis.

Results

Histology

The lamina propria to villi ratio (LP/V) was significantly influenced by SCP replacement level (F = 9.248, P = 0.0064), but not by additives or the interaction of additives and replacement level (Table 2). The LP/V ratio increased with increasing SCP levels, being greatest in fish fed the SCP75% diet with additives (8.85% ± 0.21%) and lowest those fed the SCP25% diet without additives (6.66% ± 0.45%). No significant effects of SCP level, additives, or their interaction were observed on any of the other hindgut morphology measurements (Table 2).

Table 2 Histomorphology measurements in yellowtail kingfish after 35-day growth trial

Biochemical parameters

There were significant effects of SCP replacement level on serum lipase, cholesterol, and triglycerides, with higher levels of lipase, but lower levels of cholesterol, and triglycerides as SCP replacement level increased (Table 3). There were no significant effects of dietary additives or an interaction between additives and SCP replacement level for any of these parameters. Serum urea levels were also influenced by SCP replacement level, but with a significant interaction with additives; the increase in urea levels with SCP replacement level was greater when feed contained additives (Table 3). All other biochemical parameters and measures of total antioxidant capacity showed no significant effect of SCP replacement level, dietary additives, or their interaction.

Table 3 Biochemical analysis and total antioxidant capacity measured in serum using e-BQC in yellowtail kingfish after 35-day growth trial

Myeloperoxidase (MPO) activity

MPO in the hindgut of the fish was significantly affected by SCP replacement level (F = 4.121, P = 0.024), with MPO activity decreasing as SCP levels increased (Fig. 1). There was no significant effect of dietary additives on MPO or an interaction of additives and SCP replacement level.

Fig. 1
figure 1

Myeloperoxidase activity in yellowtail kingfish hindgut after 35-day growth trial. Values are mean, with SE bars (n = 3). Different letters indicate significant differences between FM replacement levels, as determined by post hoc Tukey’s test

Gene expression

There were no significant effects of SCP replacement level or dietary additives, and no significant interaction between these terms, on gene expression levels in the intestine (Table 4) or brain (Table 5).

Table 4 Gene expressionm in the gut of yellowtail kingfish after 35-day growth trial
Table 5 Gene expression in the brain of yellowtail kingfish after 35-day growth trial

In the liver, however, there was a significant interaction effect between additives and SCP replacement level on the expression of glutathione peroxidase (gpx; F = 13.14, P = 0.003), with expression increasing with SCP replacement levels when there were no additives in the diet, but not when additives were present (Table 6). Interleukin-1 (itl1) expression in the liver showed significant effects of SCP replacement level (F = 121.34, P ≤ 0.0001), additives (F = 17.41, P ≤ 0.001), and their interaction (F = 10.13, P = 0.007). Expression increased when additives were added at 0% replacement level, but not at 25% and 75% replacement levels, and increased with increasing replacement level, although this was more marked with dietary additives (Table 6). There were no significant effects of additive inclusion, SCP replacement level, or their interaction on the expression of superoxide dismutase (sod), catalase (cat), and hepcidin (hep).

Table 6 Gene expression in the liver of yellowtail kingfish after 35-day growth trial

Microbiome

The analysis yielded over 1.4 million reads, with an average of approximately 21,928 high-quality, non-chimeric reads per sample, equating to 612 ASVs. After excluding singletons, a mean ASV frequency of over 20,000 per sample was obtained. After removing sequences related to mitochondria and chloroplasts, 12 unique phyla, 44 orders, and 65 genera were identified.

Rarefaction curve analysis indicated adequate sampling depth for over 90% of the samples, with a notable exception in fish fed the SCP50% diet, where some samples exhibited elevated lines at the lowest sampling depth, suggesting potential unique ASVs left in the community representing more species richness (observed and Chao1) (Fig. S1). Similarly, samples with dietary additives also displayed upward trending rarefaction lines (Fig. S2).

There were significant effects of both SCP replacement level (P < 0.001) and dietary additives (P = 0.016) on ASV richness, with no significant interaction. Richness was greater in fish that received dietary additives, and at intermediate (25% and 50%) replacement levels (Fig. 2A). In contrast to ASV richness, there were no significant effects of SCP replacement level or dietary additives on the Shannon-Weiner diversity index (Fig. 2B).

Fig. 2
figure 2

Alpha–beta diversity analysis of gut microbiota in yellowtail kingfish after 35-day growth trial. A and B Diversity of gut bacteria in terms of species richness and Shannon-Weiner index. Boxplots with different superscript letters at the top indicate significantly different mean values between additives (with and without) within a replacement level. C Number of shared and unique ASVs in fish-fed diets with and without additives (top), and at different SCP replacement levels (bottom). D and E PCoA plots representing beta-diversity for unweighted and weighted UniFrac distance metric

Analysis of unshared ASVs revealed that 40%–46% of ASVs were unique to samples with or without additives, with only 14% shared between the two groups. The SCP50% replacement level contributed to approximately 40% of the unshared ASVs, followed by SCP 25% (21%), SCP 75% (14.5%), and SCP 0% (10%) (Fig. 2C).

Both unweighted (R2 = 0.118, P = 0.00002) and weighted (R2 = 0.118, P = 0.042) UniFrac analyses demonstrated significant differences in microbial community composition among SCP replacement levels (P < 0.001), with no effects of dietary additives or the interaction between SCP replacement level and additives (Fig. 2D and E).

Analysis of the relative abundance of bacterial phyla revealed that Proteobacteria-dominated bacterial communities were observed across all samples, regardless of SCP replacement levels and additives. In the 0% and 25% SCP diets, Proteobacteria and Firmicutes together comprised approximately 75% of the gut microbial communities. In contrast, fish from the 50% and 75% replacement levels exhibited a notable increase in Cyanobacteria abundance, accounting for 24% and 28% of the microbial reads, respectively, alongside a consistent presence of Proteobacteria at 21%. Actinobacteriota showed a high abundance (24%) in the 50% replacement level with additives, whereas other treatments displayed uniformly low Actinobacteriota abundance (5%–15%) without additives. Firmicutes abundance was lowest (5%) in the 50% replacement level with additives but ranged from 10%–18% in all other samples (Fig. 3A, Fig. S3).

Fig. 3
figure 3

Relative abundance of bacteria at phylum (A) and genus (B) level in the gut of yellowtail kingfish after 35-day growth trial. Phyla and genera with more than 1% relative abundance were considered for plotting

At the genus level (Fig. 3B), Vibrio was more abundant (14% to 26%) in fish from the 0% and 25% replacement diet, while its presence was minimal (< 1%) in fish fed the higher replacement diets (SCP50% and SCP75%). Bradyrhizobium was consistently present across all samples, making it a ‘core’ component of the YTK gut microbiota. Sphingomonas, absent at 0% replacement level, was found in 5% to 12% of fish fed with higher SCP replacement diets. Notably, fish fed the 75% SCP diet had more opportunistic pathogens such as Staphylococcus and Corynebacterium. Additionally, 26% of bacteria in fish fed the 75% SCP replacement diet with additives remained unclassified, indicating a potential for undiscovered microbial diversity.

Fish in the 50% SCP replacement diet had a more diverse gut microbiota compared to fish in other replacement levels. Differential abundance analysis highlighted a significantly higher abundance of Vibrio and Cutibacterium at 0% replacement, whereas Sphingomonas and Candidatus obscuribacter were more abundant in fish from the 50% SCP group (Fig. 4). Fish from the 75% SCP group showed significantly higher abundance of Acinetobacter, Corynebacterium, Staphylococcus, and several unclassified genera (Fig. 4). No specific genera were enriched in fish from the 25% SCP group, underscoring the nuanced impact of SCP levels on community composition.

Fig. 4
figure 4

Differential abundance (DA) of gut microbiota in yellowtail kingfish after 35-day growth trial. An LDA score of 2.0 and above following the Kruskal–Wallis, with a Bonferroni adjustment P-value of < 0.05, was considered for DA analysis

Discussion

As global aquaculture industries continue to grow, sustainable alternatives to fishmeal are continuously being explored. Single-cell microbial protein has shown promise as an alternative with its ability to capture and utilise greenhouse gases. Our previous research demonstrated that single-cell microbial protein (SCP), capable of capturing and utilizing greenhouse gases, serves as a promising alternative in the diet of YTK. In that study, fish fed a diet with 25% SCP replacement consumed less feed than those on the control diet but achieved comparable growth rates, leading to an improved feed conversion ratio (FCR). Higher levels of SCP replacement resulted in decreased feed consumption and growth rate, but FCR, protein retention efficiency, and apparent digestibility coefficients were maintained [7]. Building on these findings, the current study found no major negative health impacts, and increased gut microbial richness when juvenile YTK were fed diets substituting fishmeal with SCP at various levels.

Health

Over the course of 35 d, this study observed no major adverse health effects from any of the diets, although there was some evidence for mild inflammation at the highest replacement level, as indicated by the lamina propria to villi area ratio, decreased myeloperoxidase (MPO) activity, and upregulation of some immune-related genes in the hindgut. The lamina propria to villi area ratio is crucial as it reflects the structural integrity and functionality of the gut, with a balanced ratio indicating effective nutrient absorption and a healthy gut environment [41]. While the ratio was elevated at increased SCP replacement levels, it still fell within a healthy range for juvenile YTK [41]. Similar results were found in previous studies on Atlantic salmon [42], barramundi [31], and black seabream [43], where no morphological damage to the gut was observed with the use of SCP in these diets. Importantly, although this study was 35 d in duration, it is sufficient to detect early signs of gut inflammation, as previous studies have shown that acute enteritis can manifest between 34–42 d [44, 45]. This suggests that while no major adverse health impacts were observed, the timeframe of this study was adequate to capture early indicators of gut health issues.

MPO activity has been used to indicate neutrophil infiltration and subsequent inflammation within the lamina propria of YTK, with a decrease in MPO levels indicative of inflammation [46]. MPO activity is a key indicator of neutrophil infiltration and inflammation within the lamina propria in YTK. Stone et al. [44] demonstrated that a reduction in MPO levels was associated with a decrease in inflammation, highlighting its utility as a marker for assessing inflammatory responses in the gut. In the current study, fish fed with the highest replacement level of SCP had significantly less MPO activity compared to fish fed with no FM replacement, suggesting that although the lamina propria area was significantly higher in the diets containing SCP, the reduced MPO levels indicate that this was not indicative of enteritis.

Serum biochemistry is routinely used in diagnostics to determine the health of aquaculture fish species. This study found significant differences in the lipase, cholesterol, and triglycerides in fish fed with SCP, but no significant changes in any of the other measured parameters. The significant rise in lipase levels in fish fed increasing levels of SCP and at the same time lower levels of triglycerides are similar to that found in barramundi fed the same SCP [31]. The reduction in serum cholesterol is likely linked to reducing dietary cholesterol as fishmeal inclusion decreased but could also be the result of improved lipid digestibility, consistent with the reduced triglycerides. Increases in urea are documented in both Atlantic salmon and barramundi fed SCP diets [22, 31] and these increased levels have been attributed to higher nucleic acids in the bacterial meal compared to fishmeal. Despite changes in these biochemical parameters with SCP inclusion, all biochemical indicators still fell within healthy ranges for YTK [47].

Although no differential gene expression was found in the intestine or brain of fish on the diets containing SCP, there were differences in expression levels of several genes related to oxidative stress and immune response in the liver. Specifically, glutathione peroxidase (gpx) showed a significant interaction between additives and SCP replacement, with the highest expression in fish fed the SCP75% diet without additives. This shows that high levels of SCP replacement can influence antioxidant enzyme activity and suggests that fish fed this diet were dealing with oxidative stress, consistent with the indications of inflammation in the hindgut. Interleukin-1 (itl1) expression was also significantly affected by SCP replacement and additives, with the highest expression in fish fed the SCP75% diet. Itl1 was included as it is a pro-inflammatory cytokine that plays a crucial role in initiating immune responses to stress or injury, including dietary stress. Previous studies in YTK have shown that inflammation in the gut can occur when fishmeal is substituted with terrestrial alternatives, such as soy[44]. This suggests that SCP replacement may play a role in modulating immune responses, possibly contributing to the mild inflammatory effects observed in histomorphological analyses. In contrast, the expression of superoxide dismutase (sod), catalase (cat), and hepcidin (hep) did not show significant changes with SCP replacement, indicating that these aspects of the oxidative stress response and immune function were not markedly affected by the dietary treatments [48].

Microbiome

Microbiome analysis of YTK gut found an increase in bacterial ASV richness and a change in ASV composition with SCP replacement in the diet, and this was most marked at replacement levels of 25% and 50%. This finding underscores the potential of a moderate level of SCP replacement to introduce beneficial new bacterial taxa to the gut microbiome, which could positively affect the host’s nutrient metabolism and immune function. Similar to how Aguilera et al. [49] observed shifts in microbial communities in response to fish growth stages, SCP inclusion may introduce beneficial new taxa that impact nutrient metabolism and immune responses in yellowtail kingfish. This diversification is particularly encouraging, given that research on cage-cultured YTK by Legrand et al. [50] highlighted a diminished microbial diversity in fish afflicted with enteritis, underscoring the potential health benefits of a richer gut microbiome. Despite the increased species richness, the lack of significant difference in the Shannon-Weiner diversity index suggests that the overall distribution of these species remains relatively uniform, which may indicate a resilient core microbiome that adapts to changes in SCP replacement levels while maintaining a balanced gut ecosystem.

The dominance of Proteobacteria across all samples is consistent with previous studies on aquatic species, indicating a stable presence of this phylum in the gut environment [51]. This aligns with the findings of Aguilera et al. [49], where Proteobacteria, particularly genera like Pseudomonas, were predominant in juvenile yellowtail kingfish. The consistent detection of Bradyrhizobium as a ‘core’ member of the gut microbiota in all samples underscores its likely contribution to the gut ecosystem’s stability and functionality in YTK. The association of Bradyrhizobium with lupin meal in the diet suggests that dietary components can have a direct influence on gut microbial composition. [52]. Despite these consistencies, the beta-diversity analysis revealed significant impacts of SCP replacement on the gut microbiota, with both presence-absence and relative abundance of bacterial taxa being affected, which could have implications for the host’s health. Notably, the reduction in Vibrio spp. in diets with higher SCP replacement is particularly encouraging, as Vibrio spp. are often associated with pathogenic outbreaks in aquaculture [53]. This reduction could signify a shift towards a microbial community composition that supports fish health and reduces disease risk. As highlighted in Aguilera et al. [49], Vibrio was one of the dominant bacterial genera in the gut of juvenile yellowtail kingfish. The reduction of this genus with SCP inclusion in our study may signify a shift toward a healthier gut microbiota. However, the increase in cyanobacterial abundance with higher SCP levels may have implications for gut health and nutrient absorption, warranting further investigation. In addition, the occurrence of opportunistic pathogens such as Staphylococcus and Corynebacterium in fish fed the SCP 75% diet raises concerns that elevated SCP levels might disturb the microbial equilibrium, potentially fostering conditions favourable for pathogen colonization. This aligns with the findings of Aguilera et al. [49], who observed Staphylococcus as a predominant genus in the early growth stages of yellowtail kingfish, suggesting that high SCP inclusion might create conditions conducive to the growth of opportunistic pathogens. Romarheim et al. [54] found that large and/or water-insoluble components found in M. capsulatus are important for maintaining the normal intestinal homeostasis in fish, further supporting SCP’s potential in gut health management. The benefits of increased microbiome diversity, such as enhanced nutrient metabolism, improved immune function, and protection against pathogens, underscore the importance of balancing SCP inclusion to optimize gut health. Greater microbial diversity helps create a stable gut environment that supports nutrient absorption and resilience to dietary changes [55]. These findings emphasize the nuanced influence of SCP replacement on the gut microbial ecosystem and highlight the complex balance between promoting beneficial microbial diversity and managing the risk of opportunistic pathogens.

Nevertheless, the differential abundance of bacterial genera across different SCP replacement levels and additives suggests that dietary manipulation can be a powerful tool in shaping the gut microbiome. The increased diversity observed in the SCP25% and SCP50% diets could be beneficial for nutrient absorption and immune modulation. Stone et al. [47] found that disease symptoms such as enteritis are associated with decreased microbial richness, diversity, and evenness, alongside a dominance of certain bacterial taxa in YTK. This parallels findings from a study in gilthead sea bream (Sparus aurata) that demonstrated increased gut microbiome diversity with 20% SCP replacements [56]. Notably, previous studies in YTK, which ran for 30–33 d, found significant shifts in microbial composition and function over similar durations to this study [49, 57]. This justifies the trial length, demonstrating that similar timeframes are adequate to capture meaningful changes in microbiome composition and associated functional outcomes. However, this benefit is juxtaposed with the potential risk of opportunistic pathogens at higher replacement levels, and balancing these factors is crucial for optimizing gut health and the overall well-being of the host species.

Conclusion

This study employed a comprehensive approach, encompassing micromorphological assessments, serum biochemistry, and gene expression analyses, to investigate the effects of SCP on fish health. Notably, the microbiome analysis sheds light on SCP’s role in shaping gut microbial diversity and composition, which are pivotal for the overall health and growth of YTK. As the quest for sustainable feed alternatives in aquaculture continues, this research marks a significant step toward diminishing the dependence on traditional fishmeal, highlighting SCP’s potential as a promising option. This is particularly relevant when considering optimal replacement levels and the palatability of SCP. Despite the lack of observed clinical health impacts within the duration of this study, the complexity of fish nutrition and health warrants further, more detailed investigations. Future research should focus on unravelling the extensive physiological, immunological, and microbial dynamics associated with prolonged SCP usage, aiming for a comprehensive understanding of its implications for fish health and growth in the context of sustainable aquaculture.

Data Availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. The raw sequence data in fastq format have been deposited to National Centre for Biotechnology Information (NCBI) and currently available under the BioProject accession number PRJNA726745.

Abbreviations

  • ALT:: Alanine aminotransferase
  • APF:: 2-[6-(4-Aminophenoxy)-3-oxo-3H-xanthen-9-yl] benzoic acid
  • ASV:: Amplicon sequence variant
  • cat:: Catalase
  • cDNA:: Complementary DNA
  • CTAB:: Cetyltrimethylammonium bromide
  • EDTA:: Ethylenediaminetetraacetic acid
  • FCR:: Feed conversion ratio
  • FM:: Fish meal
  • GLDH:: Glutamate dehydrogenase
  • gpx:: Glutathione peroxidase
  • H&E:: Hematoxylin and eosin
  • hep:: Hepcidin
  • itl1:: Interleukin-1
  • MPO:: Myeloperoxidase
  • PAS:: Periodic acid-Schiff
  • PBS:: Phosphate buffered saline
  • PERMANOVA:: Permutational multivariate analysis of variance
  • qPCR:: Quantitative polymerase chain reaction
  • REST:: Relative expression software tool
  • RFID:: Radio frequency identification
  • RNA:: Ribonucleic acid
  • rRNA:: Ribosomal RNA
  • SCP:: Single-cell protein
  • sod:: Superoxide dismutase
  • YTK:: Yellowtail kingfish

References

  1. 1.Whatmore P, Nguyen NH, Miller A, Lamont R, Powell D, D’Antignana T, et al. Genetic parameters for economically important traits in yellowtail kingfishSeriola lalandi. Aquaculture. 2013;400–401.(2013)77–84.: 77.
  2. 2.Tanner JE, Fernandes M. Environmental effects of yellowtail kingfish aquaculture in south australia. Aquac Environ Interact. 2012;1.(2012)155–65.: 155.
  3. 3.Booth MA, Allan GL, Pirozzi I. Estimation of digestible protein and energy requirements of yellowtail kingfishSeriola lalandiusing a factorial approach. Aquaculture. 2010;307.(2010)247–59.: 247.
  4. 4.Gupta N, Partridge GJ, Buller NB, Pilmer L, Currie A, Lymbery AJ. Genetic diversity of Australian isolates ofPhotobacterium damselaesubsp.damselaeis associated with virulence to yellowtail kingfish (Seriola lalandi). Aquaculture. 2021;538.(2021)736552.: 736552.
  5. 5.Horstmann P, Maas RM, de Boer XV, Staessen TWO, Kokou F, Schrama JW. Faecal waste characteristics of yellowtail kingfish (Seriola lalandi) fed with pelleted and natural feed. Anim Feed Sci Technol. 2023;299.(2023)115625.: 115625.
  6. 6.van der Heijden PGM, Rohaan AC. Kweek van Yellowtail Kingfish in recirculatiesystemen. Aquacultuur. 2019;34.(2019)35–8.: 35.
  7. 7.Pilmer LW, Woolley LD, Lymbery AJ, Salini M, Partridge GJ. Using dietary additives to improve palatability of diets containing single-cell protein from methanotrophic bacteria in yellowtail kingfish (Seriola lalandi) diets. Aquac Res. 2022.https.(2022)1111/are.15986.
  8. 8.Henry M, Gasco L, Piccolo G, Fountoulaki E. Review on the use of insects in the diet of farmed fish.(2015)Past and future.Anim Feed Sci Technol.: 1.
  9. 9.Zamani A, Khajavi M, Nazarpak MH, Gisbert E. Evaluation of a bacterial single-cell protein in compound diets for rainbow trout (Oncorhynchus mykiss) fry as an alternative protein source. Animals (Basel). 2020;10(9).(2020)org/10.3390/ani10091676.: 1676.
  10. 10.Jintasataporn O, Chumkam S, Triwutanon S, LeBlanc A, Sawanboonchun J. Effects of a single cell protein (Methylococcus capsulatus, Bath) in Pacific white shrimp (Penaeus vannamei) diet on growth performance, survival rate and resistance toVibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis diseas. Front Mar Sci. 2021;8.(2021)1643.: 1643.
  11. 11.Belghit I, Liland NS, Gjesdal P, Biancarosa I, Menchetti E, Li Y, et al. Black soldier fly larvae meal can replace fish meal in diets of sea-water phase Atlantic salmon (Salmo salar). Aquaculture. 2019;503.(2019)609–19.: 609.
  12. 12.Zhang Q, Liang H, Longshaw M, Wang J, Ge X, Zhu J, et al. Effects of replacing fishmeal with methanotroph (Methylococcus capsulatus, Bath) bacteria meal (FeedKind®) on growth and intestinal health status of juvenile largemouth bass (Micropterus salmoides). Fish Shellfish Immunol. 2022;122.(2022)298–305.: 298.
  13. 13.Biswas A, Takakuwa F, Yamada S, Matsuda A, Saville RM, LeBlanc A, et al. Methanotroph (Methylococcus capsulatus, Bath) bacteria meal as an alternative protein source for Japanese yellowtail. Seriola quinqueradiata Aquaculture. 2020;529.(2020)735700.: 735700.
  14. 14.Chen Y, Chi S, Zhang S, Dong X, Yang Q, Liu H, et al. Replacement of fish meal with methanotroph (Methylococcus capsulatus, Bath) bacteria meal in the diets of Pacific white shrimp (Litopenaeus vannamei). Aquaculture. 2021;541.(2021)736801.: 736801.
  15. 15.Hua K, Cobcroft JM, Cole A, Condon K, Jerry DR, Mangott A, et al. The future of aquatic protein.(2019)implications for protein sources in aquaculture diets.One Earth.: 316.
  16. 16.Gasco L, Gai F, Maricchiolo G, Genovese L, Ragonese S, Bottari T, et al. Fishmeal alternative protein sources for aquaculture feeds. Feeds for the aquaculture sector, Chemistry of Foods. Cham: Springer; 2018. p. 1–28.
  17. 17.Sillman J, Nygren L, Kahiluoto H, Ruuskanen V, Tamminen A, Bajamundi C, et al. Bacterial protein for food and feed generated via renewable energy and direct air capture of CO2: Can it reduce land and water use? Glob Food Sec. 2019;22:25–32.
  18. 18.Wan AHL, Davies SJ, Soler-Vila A, Fitzgerald R, Johnson MP. Macroalgae as a sustainable aquafeed ingredient. Rev Aquac. 2019;11.(2019)458–92.: 458.
  19. 19.Glencross B, Huyben D, Schrama J. The application of single-cell ingredients in aquaculture feeds—A review. Fishes. 2020;5.(2020)22.: 22.
  20. 20.Strong PJ, Xie S, Clarke WP. Methane as a resource.(2015)can the methanotrophs add value? Environ Sci Technol.: 4001.
  21. 21.Chen Y, Chi S, Zhang S, Dong X, Yang Q, Liu H, et al. Evaluation of methanotroph (Methylococcus capsulatus, Bath) bacteria meal on body composition, lipid metabolism, protein synthesis and muscle metabolites of Pacific white shrimp (Litopenaeus vannamei). Aquaculture. 2022;547.(2022)737517.: 737517.
  22. 22.Aas TS, Hatlen B, Grisdale-Helland B, Terjesen BF, Penn M, Bakke-McKellep AM, et al. Feed intake, growth and nutrient utilization in Atlantic halibut (Hippoglossus hippoglossus) fed diets containing a bacterial protein meal. Aquac Res. 2007;38.(2007)351–60.: 351.
  23. 23.Storebakken T, Baeverfjord G, Skrede A, Olli JJ, Berge GM. Bacterial protein grown on natural gas in diets for Atlantic salmon,Salmo salar, in freshwater. Aquaculture. 2004;241.(2004)413–25.: 413.
  24. 24.Biswas A, Takakuwa F, Yamada S, Furukawa A, Shiratori M, Koshiishi T, et al. Partial replacement of fishmeal by a protein-rich product from methanotrophic bacteria in juvenile red sea bream (Pagrus major). Aquacult Nutr. 2021.https.(2021)1111/anu.13398.
  25. 25.Rajesh M, Kamalam BS, Sharma P, Verma VC, Pandey A, Dubey MK, et al. Evaluation of a novel methanotroph bacteria meal grown on natural gas as fish meal substitute in rainbow trout. Oncorhynchus mykiss Aquac Res. 2022;53.(2022)2159–74.: 2159.
  26. 26.Kleiveland CR, Hult LTO, Spetalen S, Kaldhusdal M, Christofferesen TE, Bengtsson O, et al. The noncommensal bacteriummethylococcus capsulatus(Bath) ameliorates dextran sulfate (sodium salt)-induced ulcerative colitis by influencing mechanisms essential for maintenance of the colonic barrier function. Appl Environ Microbiol. 2013;79.(2013)48.: 48.
  27. 27.Øverland M, Tauson AH, Shearer K, Skrede A. Evaluation of methane-utilising bacteria products as feed ingredients for monogastric animals. Arch Anim Nutr. 2010;64(3).(2010)org/10.1080/17450391003691534.: 171.
  28. 28.Martin SAM, Król E. Nutrigenomics and immune function in fish.(2017)new insights from omics technologies.Dev Comp Immunol.: 86.
  29. 29.Wang M, Fan Z, Zhang Z, Yi M, Liu Z, Ke X, et al. Effects of diet on the gut microbial communities of Nile tilapia (Oreochromis niloticus) across their different life stages. Front Mar Sci. 2022;9.(2022)926132.: 926132.
  30. 30.Blanco AM, Calo J, Soengas JL. The gut-brain axis in vertebrates.(2021)Implications for food intake regulation.J Experiment Biol.
  31. 31.Woolley L, Chaklader MR, Pilmer L, Stephens F, Wingate C, Salini M, et al. Gas to protein.(2023)Microbial single cell protein is an alternative to fishmeal in aquaculture.Sci Total Environ.: 160141.
  32. 32.Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8).(2019)852–7.: 852.
  33. 33.Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2.(2016)High-resolution sample inference from Illumina amplicon data.Nat Methods.: 581.
  34. 34.Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7.(2013)Improvements in performance and usability.Mol Biol Evol.: 772.
  35. 35.Price MN, Dehal PS, Arkin AP. FastTree 2 – Approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5.(2010)e9490.
  36. 36.Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project.(2013)improved data processing and web-based tools.Nucleic Acids Res.
  37. 37.McMurdie PJ, Holmes S. phyloseq.(2013)An R package for reproducible interactive analysis and graphics of microbiome census data.PLoS ONE.
  38. 38.Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14.(2003)927–30.: 927.
  39. 39.Liu C, Cui Y, Li X, Yao M. microeco.(2021)an R package for data mining in microbial community ecology.FEMS Microbiol Ecol.: 255.
  40. 40.Cao Y, Dong Q, Wang D, Zhang P, Liu Y, Niu C. microbiomeMarker.(2022)an R/Bioconductor package for microbiome marker identification and visualization.Bioinformatics.: 4027.
  41. 41.Candebat CL, Stephens F, Booth MA, Fernando F, Lopata A, Pirozzi I. Adequate levels of dietary sulphur amino acids impart improved liver and gut health in juvenile yellowtail kingfish (Seriola lalandi). Br J Nutr. 2023;129(8).(2023)org/10.1017/S0007114522002458.: 1289.
  42. 42.Berge G, Baeverfjord G, Skrede A, Storebakken T. Bacterial protein grown on natural gas as protein source in diets for Atlantic salmon,Salmo salar, in saltwater. Aquaculture. 2005;244.(2005)233–40.: 233.
  43. 43.Xu B, Liu Y, Chen K, Wang L, Sagada G, Tegomo AF, et al. Evaluation of methanotroph (Methylococcus capsulatus, Bath) bacteria meal (FeedKind®) as an alternative protein source for juvenile black sea bream. Acanthopagrus schlegelii Front Mar Sci. 2021;8.(2021)1702.: 1702.
  44. 44.Stone DAJ, Bellgrove EJ, Forder REA, Howarth GS, Bansemer MS. Inducing subacute enteritis in yellowtail kingfishSeriola lalandi.(2018)the effect of dietary inclusion of soybean meal and grape seed extract on hindgut morphology and inflammation.N Am J Aquac.: 59.
  45. 45.Bansemer MS, Forder REA, Howarth GS, Suitor GM, Bowyer J, Stone DAJ. The effect of dietary soybean meal and soy protein concentrate on the intestinal mucus layer and development of subacute enteritis in Yellowtail Kingfish (Seriola lalandi) at suboptimal water temperature. Aquac Nutr. 2015;21.(2015)1111/anu.12160.: 300.
  46. 46.Faith M, Sukumaran A, Pulimood AB, Jacob M. How reliable an indicator of inflammation is myeloperoxidase activity? Clin Chim Acta. 2008;396.(2008)23–5.: 23.
  47. 47.Stone DAJ, Booth MA, Clarke SM. Growing a profitable, innovative and collaborative Australian yellowtail kingfish aquaculture industry bringing ‘white’ fish to the market. Fisheries Research and Development Corporation, South Australian Research and Development Institute and New South Wales Department of Primary Industries; 2019.
  48. 48.Deng J, Kang B, Tao L, Rong H, Zhang X. Effects of dietary cholesterol on antioxidant capacity, non-specific immune response, and resistance toAeromonas hydrophilain rainbow trout (Oncorhynchus mykiss) fed soybean meal-based diets. Fish Shellfish Immunol. 2013;34.(2013)324–31.: 324.
  49. 49.Aguilera E, Yany G, Romero J. Cultivable intestinal microbiota of yellowtail juveniles (Seriola lalandi) in an aquaculture system. Lat Am J Aquat Res. 2013;41(3).(2013)395–403.: 395.
  50. 50.Legrand TPRA, Wynne JW, Weyrich LS, Oxley APA. Investigating both mucosal immunity and microbiota in response to gut enteritis in yellowtail kingfish. Microorganisms. 2020;8(9).(2020)1267.: 1267.
  51. 51.Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, et al. Environmental and ecological factors that shape the gut bacterial communities of fish.(2012)A meta-analysis.Mol Ecol.: 3363.
  52. 52.Songwattana P, Chaintreuil C, Wongdee J, Teulet A, Mbaye M, Piromyou P, et al. Identification of type III effectors modulating the symbiotic properties ofBradyrhizobium vignaestrain ORS3257 with variousVignaspecies. Sci Rep. 2021;11.(2021)org/10.1038/s41598-021-84205-w.: 4874.
  53. 53.Sanches-Fernandes GMM, Sá-Correia I, Costa R. Vibriosis Outbreaks in Aquaculture.(2022)Addressing environmental and public health concerns and preventive therapies using gilthead seabream farming as a model system.Front Microbiol.: 904815.
  54. 54.Romarheim OH, Landsverk T, Mydland LT, Skrede A, Øverland M. Cell wall fractions fromMethylococcus capsulatusprevent soybean meal-induced enteritis in Atlantic salmon (Salmo salar). Aquaculture. 2013;402–403.(2013)13–8.: 13.
  55. 55.Tarnecki AM, Burgos FA, Ray CL, Arias CR. Fish intestinal microbiome.(2017)diversity and symbiosis unravelled by metagenomics.J Appl Microbiol.: 2.
  56. 56.Marchi A, Bonaldo A, Scicchitano D, Candela M, De Marco A, Falciglia S, et al. Feeding gilthead sea bream with increasing dietary bacterial single cell protein level.(2023)Implication on growth, plasma biochemistry, gut histology, and gut microbiota.Aquaculture.: 739132.
  57. 57.Horlick J, Booth MA, Tetu SG. Alternative dietary protein and water temperature influence the skin and gut microbial communities of yellowtail kingfish (Seriola lalandi). PeerJ. 2020;8.(2020)7717/peerj.8705.

Acknowledgements

We would like to thank the Department of Primary Industries and Regional Development (DPIRD) staff for supplying the fish and facilities for this project, Ridley Agriproducts and Specialty feeds for helping develop the diets.

Funding

Fisheries Research Development Corporation, Grant/Award Number: 2017-030.

Ethics Declaration

Ethics approval and consent to participate

All experiments were conducted under the Department of Primary Industries and Regional Development (DPIRD) Animal Ethics Committee approval (permit number: AEC 20-3-07).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Rights and Permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ . The Creative Commons Public Domain Dedication waiver ( http://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Reprints and permissions