DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

Bivalves Present the Largest and Most Diversified Repertoire of Toll-Like Receptors in the Animal Kingdom, Suggesting Broad-Spectrum Pathogen Recognition in Marine Waters.

Toll-like receptors (TLRs) are the most widespread class of membrane-bound innate immune receptors, responsible of specific pathogen recognition and production of immune effectors through the activation of intracellular signaling cascades. The repertoire of TLRs was analyzed in 85 metazoans, enriched on molluscan species, an underrepresented phylum in previous studies. Following an ancient evolutionary origin, suggested by the presence of TLR genes in Anthozoa (Cnidaria), these receptors underwent multiple independent gene family expansions, the most significant of which occurred in bivalve molluscs. Marine mussels (Mytilus spp.) had the largest TLR repertoire in the animal kingdom, with evidence of several lineage-specific expanded TLR subfamilies with different degrees of orthology conservation within bivalves. Phylogenetic analyses revealed that bivalve TLR repertoires were more diversified than their counterparts in deuterostomes or ecdysozoans. The complex evolutionary history of TLRs, characterized by lineage-specific expansions and losses, along with episodic positive selection acting on the extracellular recognition domains, suggests that functional diversification might be a leading evolutionary force. We analyzed a comprehensive transcriptomic data set from Mytilus galloprovincialis and built transcriptomic correlation clusters with the TLRs expressed in gills and in hemocytes. The implication of specific TLRs in different immune pathways was evidenced, as well as their specific modulation in response to different biotic and abiotic stimuli. We propose that, in a similar fashion to the remarkable functional specialization of vertebrate TLRs, the expansion of the TLR gene family in bivalves attends to a functional specification motivated by the biological particularities of these organisms and their living environment.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Multifunctional Cell Regulation Activities of the Mussel Lectin SeviL: Induction of Macrophage Polarization toward the M1 Functional Phenotype.

SeviL, a galactoside-binding lectin previously isolated from the mussel Mytilisepta virgata, was demonstrated to trigger apoptosis in HeLa ovarian cancer cells. Here, we show that this lectin can promote the polarization of macrophage cell lines toward an M1 functional phenotype at low concentrations. The administration of SeviL to monocyte and basophil cell lines reduced their growth in a dose-dependent manner. However, low lectin concentrations induced proliferation in the RAW264.7 macrophage cell line, which was supported by the significant up-regulation of TOM22, a component of the mitochondrial outer membrane. Furthermore, the morphology of lectin-treated macrophage cells markedly changed, shifting from a spherical to an elongated shape. The ability of SeviL to induce the polarization of RAW264.7 cells to M1 macrophages at low concentrations is supported by the secretion of proinflammatory cytokines and chemokines, as well as by the enhancement in the expression of IL-6- and TNF-α-encoding mRNAs, both of which encode inflammatory molecular markers. Moreover, we also observed a number of accessory molecular alterations, such as the activation of MAP kinases and the JAK/STAT pathway and the phosphorylation of platelet-derived growth factor receptor-α, which altogether support the functional reprogramming of RAW264.7 following SeviL treatment. These results indicate that this mussel ß-trefoil lectin has a concentration-dependent multifunctional role in regulating cell proliferation, phenotype, and death in macrophages, suggesting its possible involvement in regulating hemocyte activity in vivo.

A high-quality reference genome for the critically endangered Aeolian wall lizard, Podarcis raffonei.

The Aeolian wall lizard, Podarcis raffonei, is an endangered species endemic to the Aeolian archipelago, Italy, where it is present only in 3 tiny islets and a narrow promontory of a larger island. Because of the extremely limited area of occupancy, severe population fragmentation and observed decline, it has been classified as Critically Endangered by the International Union for the Conservation of Nature (IUCN). Using Pacific Biosciences (PacBio) High Fidelity (HiFi) long-read sequencing, Bionano optical mapping and Arima chromatin conformation capture sequencing (Hi-C), we produced a high-quality, chromosome-scale reference genome for the Aeolian wall lizard, including Z and W sexual chromosomes. The final assembly spans 1.51 Gb across 28 scaffolds with a contig N50 of 61.4 Mb, a scaffold N50 of 93.6 Mb, and a BUSCO completeness score of 97.3%. This genome constitutes a valuable resource for the species to guide potential conservation efforts and more generally for the squamate reptiles that are underrepresented in terms of available high-quality genomic resources.

Taxonomic Distribution and Molecular Evolution of Mytilectins.

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a ß-trefoil fold. Mytilectins represent a recently described subgroup of ß-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the ß-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

Biochemical and Functional Characterization of the Three Zebrafish Transglutaminases 2.

Transglutaminase 2 (TG2) is a multifunctional protein widely distributed in various tissues and involved in many physiological and pathological processes. However, its actual role in biological processes is often controversial as TG2 shows different effects in these processes depending on its localization, cell type, or experimental conditions. We characterized the enzymatic and functional properties of TG2 proteins expressed in Danio rerio (zebrafish) to provide the basis for using this established animal model as a reliable tool to characterize TG2 functions in vivo. We confirmed the existence of three genes orthologous to human TG2 (zTGs2) in the zebrafish genome and their expression and function during embryonic development. We produced and purified the zTGs2s as recombinant proteins and showed that, like the human enzyme, zTGs2 catalyzes a Ca2+ dependent transamidation reaction that can be inhibited with TG2-specific inhibitors. In a cell model of human fibroblasts, we also demonstrated that zTGs2 can mediate RGD-independent cell adhesion in the extracellular environment. Finally, we transfected and selected zTGs2-overexpressing HEK293 cells and demonstrated that intracellular zTGs2 plays a very comparable protective/damaging role in the apoptotic process, as hTG2. Overall, our results suggest that zTGs2 proteins behave very similarly to the human ortholog and pave the way for future in vivo studies of TG2 functions in zebrafish.

Comparative analysis of novel and common reference genes in adult tissues of the mussel Mytilus galloprovincialis.

BACKGROUND: Real-time quantitative PCR is a widely used method for gene expression analyses in various organisms. Its accuracy mainly relies on the correct selection of reference genes. Any experimental plan involving real-time PCR needs to evaluate the characteristics of the samples to be examined and the relative stability of reference genes. Most studies in mollusks rely on reference genes commonly used in vertebrates. RESULTS: In this study, we focused on the transcriptome of the bivalve mollusk Mytilus galloprovincialis in physiological state to identify suitable reference genes in several adult tissues. Candidate genes with highly stable expression across 51 RNA-seq datasets from multiple tissues were selected through genome-wide bioinformatics analysis. This approach led to the identification of three genes (Rpl14, Rpl32 and Rpl34), whose suitability was evaluated together with 7 other reference genes commonly reported in literature (Act, Cyp-A, Ef1α, Gapdh, 18S, 28S and Rps4). The stability analyses performed with geNorm, NormFinder and Bestkeeper identified specific either single or pairs of genes suitable as references for gene expression analyses in specific tissues and revealed the Act/Cyp-A pair as the most appropriate to analyze gene expression across different tissues. CONCLUSION: Mytilus galloprovincialis is a model system increasingly used in ecotoxicology and molecular studies. Our transcriptome-wide approach represents the first comprehensive investigation aimed at the identification of suitable reference genes for expression studies in this species.

First Insights into the Repertoire of Secretory Lectins in Rotifers.

Due to their high biodiversity and adaptation to a mutable and challenging environment, aquatic lophotrochozoan animals are regarded as a virtually unlimited source of bioactive molecules. Among these, lectins, i.e., proteins with remarkable carbohydrate-recognition properties involved in immunity, reproduction, self/nonself recognition and several other biological processes, are particularly attractive targets for biotechnological research. To date, lectin research in the Lophotrochozoa has been restricted to the most widespread phyla, which are the usual targets of comparative immunology studies, such as Mollusca and Annelida. Here we provide the first overview of the repertoire of the secretory lectin-like molecules encoded by the genomes of six target rotifer species: Brachionus calyciflorus, Brachionus plicatilis, Proales similis (class Monogononta), Adineta ricciae, Didymodactylos carnosus and Rotaria sordida (class Bdelloidea). Overall, while rotifer secretory lectins display a high molecular diversity and belong to nine different structural classes, their total number is significantly lower than for other groups of lophotrochozoans, with no evidence of lineage-specific expansion events. Considering the high evolutionary divergence between rotifers and the other major sister phyla, their widespread distribution in aquatic environments and the ease of their collection and rearing in laboratory conditions, these organisms may represent interesting targets for glycobiological studies, which may allow the identification of novel carbohydrate-binding proteins with peculiar biological properties.

Coming New Age of Marine Glycomics: The Fundamental, Medical, and Ecological Aspects.

This Special Issue "Marine Glycomics" (https://www.mdpi.com/journal/marinedrugs/special_issues/Marine_Glycomics, accessed on 12 September 2022) provided new approachesand information on bioactive compounds, such as glycans and lectins from marine animals,seaweeds, and microorganisms for the application of clinical therapy and elucidation of thephysiological functions of marine organisms [...].

A Survey on the Distribution of Ovothiol and ovoA Gene Expression in Different Tissues and Cells: A Comparative Analysis in Sea Urchins and Mussels.

Ovothiols are histidine-derived thiols produced by a variety of marine invertebrates, protists and bacteria. These compounds, which are among the strongest natural antioxidants, are involved in controlling the cellular redox balance due to their redox exchange with glutathione. Although ovothiols were initially reported as protective agents against environmental stressors, new evidence suggests that they can also act as pheromones and participate in fundamental biological processes such as embryogenesis. To get further insight into the biological roles of ovothiols, we compared ovothiol biosynthesis in the sea urchin Paracentrotus lividus and in the mussel Mytilus galloprovincialis, the two species that represent the richest sources of these compounds among marine invertebrates. Ovothiol content was measured in different tissues and in the immune cells from both species and the expression levels of ovoA, the gene responsible for ovothiol biosynthesis, was inferred from publicly available transcriptomes. A comparative analysis of ovothiol biosynthesis in the two species allowed the identification of the tissues and cells synthesizing the metabolite and highlighted analogies and differences between sea urchins and mussels. By improving our knowledge on the biological roles of ovothiols and pointing out the existence of sustainable natural sources for their isolation, this study provides the basis for future biotechnological investigations on these valuable compounds.

Morphological analysis of erythrocytes of an Antarctic teleost under heat stress: Bias of the stabling effect.

The stenothermal Antarctic fish that live in the coastal waters of the Terra Nova Bay (Ross Sea) are rarely exposed to temperatures above zero during the year. We tested whether a slight temperature rise of 1.5 °C affects a sensitive biomarker such as erythrocytes morphology in sections of blood pellets of a small demersal notothen. The erythrocytes' shape descriptors showed significant or highly significant differences temporally from the capture of fish to the conclusion of the experiment. Surprisingly, the erythrocyte's morphology did not show significant differences between the two experimental conditions, returning similar results in control fish stabled at -0.9 °C and in the fish treated at +0.6 °C, although the values of the shape descriptors were often lower in the latter. This study demonstrates the critical issues of comparative physiology in the study of extremely sensitive organisms, such as the fish of the High Antarctic Zone. Moreover, the stabling effect inside the aquarium facilities appears to significantly obscure the effects of the experimental heat treatment.

Transcriptional Contribution of Transposable Elements in Relation to Salinity Conditions in Teleosts and Silencing Mechanisms Involved.

Fish are an interesting taxon comprising species adapted to a wide range of environments. In this work, we analyzed the transcriptional contribution of transposable elements (TEs) in the gill transcriptomes of three fish species exposed to different salinity conditions. We considered the giant marbled eel Anguilla marmorata and the chum salmon Oncorhynchus keta, both diadromous, and the marine medaka Oryzias melastigma, an euryhaline organism sensu stricto. Our analyses revealed an interesting activity of TEs in the case of juvenile eels, commonly adapted to salty water, when exposed to brackish and freshwater conditions. Moreover, the expression assessment of genes involved in TE silencing mechanisms (six in heterochromatin formation, fourteen known to be part of the nucleosome remodeling deacetylase (NuRD) complex, and four of the Argonaute subfamily) unveiled that they are active. Finally, our results evidenced for the first time a krüppel-associated box (KRAB)-like domain specific to actinopterygians that, together with TRIM33, might allow the functioning of NuRD complex also in fish species. The possible interaction between these two proteins was supported by structural prediction analyses.

The evolutionary puzzle solution for the origins of the partial loss of the Cτ2 exon in notothenioid fishes.

Cryonotothenioidea is the main group of fishes that thrive in the extremely cold Antarctic environment, thanks to the acquisition of peculiar morphological, physiological and molecular adaptations. We have previously disclosed that IgM, the main immunoglobulin isotype in teleosts, display typical cold-adapted features. Recently, we have analyzed the gene encoding the heavy chain constant region (CH) of the IgT isotype from the Antarctic teleost Trematomus bernacchii (family Nototheniidae), characterized by the near-complete deletion of the CH2 domain. Here, we aimed to track the loss of the CH2 domain along notothenioid phylogeny and to identify its ancestral origins. To this end, we obtained the IgT gene sequences from several species belonging to the Antarctic families Nototheniidae, Bathydraconidae and Artedidraconidae. All species display a CH2 remnant of variable size, encoded by a short Cτ2 exon, which retains functional splicing sites and therefore is included in the mature transcript. We also considered representative species from the three non-Antarctic families: Eleginopsioidea (Eleginops maclovinus), Pseudaphritioidea (Pseudaphritis urvillii) and Bovichtidae (Bovichtus diacanthus and Cottoperca gobio). Even though only E. maclovinus, the sister taxa of Cryonotothenioidea, shared the partial loss of Cτ2, the other non-Antarctic notothenioid species displayed early molecular signatures of this event. These results shed light on the evolutionary path that underlies the origins of this remarkable gene structural modification.

Comparative Genomics Reveals 13 Different Isoforms of Mytimycins (A-M) in Mytilus galloprovincialis.

Mytimycins are cysteine-rich antimicrobial peptides that show antifungal properties. These peptides are part of the immune network that constitutes the defense system of the Mediterranean mussel (Mytilus galloprovincialis). The immune system of mussels has been increasingly studied in the last decade due to its great efficiency, since these molluscs, particularly resistant to adverse conditions and pathogens, are present all over the world, being considered as an invasive species. The recent sequencing of the mussel genome has greatly simplified the genetic study of some of its immune genes. In the present work, we describe a total of 106 different mytimycin variants in 16 individual mussel genomes. The 13 highly supported mytimycin clusters (A-M) identified with phylogenetic inference were found to be subject to the presence/absence variation, a widespread phenomenon in mussels. We also identified a block of conserved residues evolving under purifying selection, which may indicate the "functional core" of the mature peptide, and a conserved set of 10 invariable plus 6 accessory cysteines which constitute a plastic disulfide array. Finally, we extended the taxonomic range of distribution of mytimycins among Mytilida, identifying novel sequences in M. coruscus, M. californianus, P. viridis, L. fortunei, M. philippinarum, M. modiolus, and P. purpuratus.

Transposable Elements and Teleost Migratory Behaviour.

Transposable elements (TEs) represent a considerable fraction of eukaryotic genomes, thereby contributing to genome size, chromosomal rearrangements, and to the generation of new coding genes or regulatory elements. An increasing number of works have reported a link between the genomic abundance of TEs and the adaptation to specific environmental conditions. Diadromy represents a fascinating feature of fish, protagonists of migratory routes between marine and freshwater for reproduction. In this work, we investigated the genomes of 24 fish species, including 15 teleosts with a migratory behaviour. The expected higher relative abundance of DNA transposons in ray-finned fish compared with the other fish groups was not confirmed by the analysis of the dataset considered. The relative contribution of different TE types in migratory ray-finned species did not show clear differences between oceanodromous and potamodromous fish. On the contrary, a remarkable relationship between migratory behaviour and the quantitative difference reported for short interspersed nuclear (retro)elements (SINEs) emerged from the comparison between anadromous and catadromous species, independently from their phylogenetic position. This aspect is likely due to the substantial environmental changes faced by diadromous species during their migratory routes.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency.

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

Identification, molecular characterization and functional analysis of interleukin (IL)-2 and IL-2like (IL-2L) cytokines in sea bass (Dicentrarchus labrax L.).

In mammals, interleukin (IL)-2, initially known as a T-cell grow factor, is an immunomodulatory cytokine involved in the proliferation of T cells upon antigen activation. In bony fish, some IL-2 orthologs have been identified, but, recently, an additional IL-2like (IL-2L) gene has been found. In this paper, we report the presence of these two divergent IL-2 isoforms in sea bass (Dicentrarchus labrax L.). Genomic analyses revealed that they originated from a gene duplication event, as happened in most percomorphs. These two IL-2 paralogs show differences in the amino acid sequence and in the exon 4 size, and these features could be an indication that they bind preferentially to different specific IL-2 receptors. Sea bass IL-2 paralogs are highly expressed in gut and spleen, which are tissues and organs involved in fish T cell immune functions, and the two cytokines could be up-regulated by both PHA stimulation and vaccination with a bacterial vaccine, with IL-2L being more inducible. To investigate the functional activities of sea bass IL-2 and IL-2L we produced the corresponding recombinant molecules in E. coli and used them to in vitro stimulate HK and spleen leukocytes. IL-2L is able to up-regulate the expression of markers related to different T cell subsets (Th1, Th2 and Th17) and to Treg cells in HK, whereas it has little effect in spleen. IL-2 is not active on these markers in HK, but shows an effect on Th1 markers in spleen. Finally, the stimulation with recombinant IL-2 and IL-2L is also able to induce in vitro proliferation of HK- and spleen-derived leukocytes. In conclusion, we have demonstrated that sea bass possess two IL-2 paralogs that likely have an important role in regulating T cell development in this species and that show distinct bioactivities.

Identification of an IgD/IgT chimera in the European sea bass (Dicentrarchus labrax L.).

Three classes of immunoglobulins have been identified in Teleosts: IgM, IgT/Z and IgD. They are fundamental for fish immune responses and, therefore, their functional activities are heavily investigated. In this paper, we describe the identification of a new IgD/IgT chimera in sea bass (Dicentrarchus labrax) from a gills transcriptome. This transcript joined the first six constant domains of the IgD chain with the two terminal constant domains of IgT, generating a long in-frame coding sequence with a junction between the canonical δ6 exon splicing donor site and the τ3 exon splicing acceptor site. Studies performed on genomic DNA confirmed the presence of the sequence and identifies and intronic region of 656 bp within this joining region. The basal expression of the IgD/IgT chimera was investigated both in silico and in vivo: high level of expression was found in gills, gut and head kidney. Moreover, IgD/IgT transcripts were up-regulated after in vitro stimulation of sea bass HK leukocytes with LPS. The IgD/IgT chimera was found also in two congener species, Morone saxatilis and Morone chrysops. It is not possible to have a precise idea on the evolutionary scenario that lead to the appearance of this sequence due to the lack of genomic information, but we could speculate that an ancestral duplication of the entire IgH locus was followed by the chimerization of Cδ/Cτ in one of the two loci. Finally, the IgD/IgT high basal expression in tissues and organs fundamental for sea bass immune response and its modulation after LPS stimulation provide a very preliminary indication that this unusual Ig variant could have a functional activity.

Differential Effects of Dietary Supplementation of Krill Meal, Soybean Meal, Butyrate, and Bactocell® on the Gene Expression of Atlantic Salmon Head Kidney.

The head kidney is a key organ that plays a fundamental role in the regulation of the fish immune response and in the maintenance of endocrine homeostasis. Previous studies indicate that the supplementation of exogenous dietary components, such as krill meal (KM), soybean meal (SM), Bactocell® (BA), and butyrate (BU), can have a significant effect on the immune function of the head kidney. The aim of this study was to investigate the differential effect of these four dietary ingredients on the transcriptional profiles of the head kidney of the Atlantic salmon. This study revealed that just a small number of genes were responsive to the feeding regime after a long-term (12 weeks) treatment, and evidenced that the most significant alterations, both in terms of the number of affected genes and magnitude of changes in gene expression, were detectable in the BU- and KM-fed groups compared with controls, while the SM diet had a nearly negligible effect, and BA had no significant effects at all. Most of the differentially expressed genes were involved in the immune response and, in line with data previously obtained from pyloric caeca, major components of the complement system were significantly affected. These alterations were accompanied by an increase in the density of melanomacrophage centers in the KM- and SM-fed group and their reduction in the BU-fed group. While three types of dietary supplements (BU, KM, and SM) were able to produce a significant modulation of some molecular players of the immune system, the butyrate-rich diet was revealed as the one with the most relevant immune-stimulating properties in the head kidney. These preliminary results suggest that further investigations should be aimed towards the elucidation of the potential beneficial effects of butyrate and krill meal supplementation on farmed salmon health and growth performance.