Chronic hyperosmotic stress inhibits renal Toll-Like Receptors expression in striped catfish (Pangasianodon hypophthalmus, Sauvage) exposed or not to bacterial infection

Effect of salinity on the physiological response and transcriptome of spotted seabass (Lateolabrax maculatus)

Salinity fluctuations significantly impact the reproduction, growth, development, as well as physiological and metabolic activities of fish. To explore the osmoregulation mechanism of aquatic organisms acclimating to salinity stress, the physiological and transcriptomic characteristics of spotted seabass (Lateolabrax maculatus) in response to varying salinity gradients were investigated. In this study, different salinity stress exerted inhibitory effects on lipase activity, while the impact on amylase activity was not statistically significant. Notably, a moderate increase in salinity (24 psu) demonstrated the potential to enhance the efficient utilization of proteins by spotted seabass. Both Na+/K+-ATPase and malondialdehyde showed a fluctuating trend of increasing and then decreasing, peaking at 72 h. Transcriptomic analysis revealed that most differentially expressed genes were involved in energy metabolism, signal transduction, the immune response, and osmoregulation. These results will provide insights into the molecular mechanisms of salinity adaptation and contribute to sustainable development of the global aquaculture industry.

A novel rare variant of CNPY3 from familial NMOSD impairs the TLR-mediated immune response

Toll-like receptors (TLRs) are a class of proteins that play essential roles in innate and adaptive immune responses. Recently, accumulating evidence has demonstrated that impairments in the TLR signalling pathway contribute to the development and progression of neuroimmune diseases, such as neuromyelitis optica spectrum disorder (NMOSD). However, the cellular and molecular mechanisms are still largely unknown. In this study, we report a novel variant, C52Y, of canopy FGF signalling regulator 3 (CNPY3) from patients with familial NMOSD and demonstrate that this variant shows a stronger interaction with GP96 and TLRs than with wild-type CNPY3. We find that C52Y has dominant negative effects on TLR4 surface expression. Importantly, the TLR4 surface expression level is decreased in RAW264.7 cells infected with the C52Y virus upon LPS stimulation. We further demonstrate that bone marrow-derived macrophages (BMDMs) from CNPY3C52Y/+ transgenic mice secrete less tumour necrosis factor (TNF) and interleukin (IL)-6 than BMDMs from wild-type mice upon stimulation with LPS. These data suggest that impairment of TLR trafficking may contribute to the development of neuroimmune disorders.

Effects of dietary supplementation of Pseudomonas aeruginosa FARP72 on the immunomodulation and resistance to Edwardsiella tarda in Pangasius pangasius

•

Pseudomonas aeruginosa FARP72 (PA) supplemented diet impacts the innate immunity.

•

Innate immune responses of P. pangasius are differentially stimulated by PA diet.

•

It improved the resistance of P. pangasius against E. tarda infection significantly.

•

IL-1β & C3 genes upregulated significantly in kidney of PA diet-fed challenged fish.

•

Highest upregulation of transferrin seen in liver of PA diet-fed challenged fish.

Edwardsiella tarda is one of the serious bacterial pathogens infecting both cultured and wild catfish urging an immediate need for effective protection strategies. This study assessed the effects of dietary supplementation of Pseudomonas aeruginosa FARP72 at 10⁸ cells/g feed (PA diet) for 30 days on the innate immunity parameters, viz., respiratory oxidative burst (ROB) activity, lysozyme, ceruloplasmin, myeloperoxidase, in-vitro nitric oxide (NO) production in addition to the expression of immune genes encoding interleukin-1β, C3 and transferrin in yellowtail catfish Pangasius pangasius and their resistance to Edwardsiella tarda challenge at a sub-lethal dose of 1.50 × 10⁷ cells/fish. A significant increase in the innate immunity parameters was noted in PA diet-fed catfish on 30 dpf compared to the control. Post E. tarda challenge, the levels of immune parameters increased significantly and peaked at 5 dpi irrespective of feeding to confer protection against E. tarda. Their levels, however, decreased on and from 10 dpi. The results on the expression of immune genes encoding interleukin-1β, C3 and transferrin in the kidney and liver tissue samples of PA diet-fed P. pangasius upon challenge with E. tarda further confirmed the ability of P. aeruginosa to stimulate primary immune organs at the gene level. The effects of feeding P. aeruginosa FARP72 on the immune functions of catfish as examined by the functional immune assays, thus, demonstrating the innate immune responses of catfish that are differentially stimulated by the PA diet. The findings of our study would help evolve management strategies to confer protection against E. tarda infection in commercial catfish aquaculture.

Salinity significantly affects intestinal microbiota and gene expression in striped catfish juveniles

In the present study, juvenile striped catfish (Pangasianodon hypophthalmus), a freshwater fish species, have been chronically exposed to a salinity gradient from freshwater to 20 psu (practical salinity unit) and were sampled at the beginning (D20) and the end (D34) of exposure. The results revealed that the intestinal microbial profile of striped catfish reared in freshwater conditions were dominated by the phyla Bacteroidetes, Firmicutes, Proteobacteria, and Verrucomicrobia. Alpha diversity measures (observed OTUs (operational taxonomic units), Shannon and Faith's PD (phylogenetic diversity)) showed a decreasing pattern as the salinities increased, except for the phylogenetic diversity at D34, which was showing an opposite trend. Furthermore, the beta diversity between groups was significantly different. Vibrio and Akkermansia genera were affected differentially with increasing salinity, the former being increased while the latter was decreased. The genus Sulfurospirillium was found predominantly in fish submitted to salinity treatments. Regarding the host response, the fish intestine likely contributed to osmoregulation by modifying the expression of osmoregulatory genes such as nka1a, nka1b, slc12a1, slc12a2, cftr, and aqp1, especially in fish exposed to 15 and 20 psu. The expression of heat shock proteins (hsp) hsp60, hsp70, and hsp90 was significantly increased in fish reared in 15 and 20 psu. On the other hand, the expression of pattern recognition receptors (PRRs) were inhibited in fish exposed to 20 psu at D20. In conclusion, the fish intestinal microbiota was significantly disrupted in salinities higher than 10 psu and these effects were proportional to the exposure time. In addition, the modifications of intestinal gene expression related to ion exchange and stressful responses may help the fish to adapt hyperosmotic environment. KEY POINTS: • It is the first study to provide detailed information on the gut microbiota of fish using the amplicon sequencing method. • Salinity environment significantly modified the intestinal microbiota of striped catfish. • Intestinal responses may help the fish adapt to hyperosmotic environment.

RNA Sequencing Analysis Reveals Divergent Adaptive Response to Hypo- and Hyper-Salinity in Greater Amberjack (Seriola dumerili) Juveniles

Salinity significantly affects physiological and metabolic activities, breeding, development, survival, and growth of marine fish. The greater amberjack (Seriola dumerili) is a fast-growing species that has immensely contributed to global aquaculture diversification. However, the tolerance, adaptation, and molecular responses of greater amberjack to salinity are unclear. This study reared greater amberjack juveniles under different salinity stresses (40, 30, 20, and 10 ppt) for 30 days to assess their tolerance, adaptation, and molecular responses to salinity. RNA sequencing analysis of gill tissue was used to identify genes and biological processes involved in greater amberjack response to salinity stress at 40, 30, and 20 ppt. Eighteen differentially expressed genes (DEGs) (nine upregulated and nine downregulated) were identified in the 40 vs. 30 ppt group. Moreover, 417 DEGs (205 up-regulated and 212 down-regulated) were identified in the 20 vs. 30 ppt group. qPCR and transcriptomic analysis indicated that salinity stress affected the expression of genes involved in steroid biosynthesis (ebp, sqle, lss, dhcr7, dhcr24, and cyp51a1), lipid metabolism (msmo1, nsdhl, ogdh, and edar), ion transporters (slc25a48, slc37a4, slc44a4, and apq4), and immune response (wnt4 and tlr5). Furthermore, KEGG pathway enrichment analysis showed that the DEGs were enriched in steroid biosynthesis, lipids metabolism, cytokine-cytokine receptor interaction, tryptophan metabolism, and insulin signaling pathway. Therefore, this study provides insights into the molecular mechanisms of marine fish adaptation to salinity.

Effects of salinity on gills' chloride cells, stress indices, and gene expression of Asian seabass (Lates calcarifer, Bloch, 1790)

A 2-week research was carried out to assess water salinity (WS) effects including 0, 15, 35, and 50‰ on osmoregulatory mechanisms and stress indices in Asian sea bass (34.4 g) juveniles. Except for fish reared at 50‰, in the other treatments, it gradually decreased to the prescribed WS during a 10-day period (- 5‰ a day). After a 10-day acclimation period, fish were reared at the prescribed WS for 2 weeks. Fish reared at 15 and 35‰ had higher chloride cell (CC) counts in the interlamellar region. The number of CC in the interlamellar region elevated with increment of WS up to 35‰, but they were pronouncedly reduced in 50‰ group. The diameter of CC in the interlamellar region was not affected by WS. The smallest nucleus diameter of CC in the interlamellar region was observed in fish reared at 15‰ (P < 0.05). The largest and the smallest amounts of serum aspartate aminotransferase were observed in fish reared at freshwater and 15‰, respectively. Fish reared at 35‰ had the highest serum sodium and potassium contents. Serum chloride content and total osmolality increased with increment of WS (P < 0.05). Serum cortisol and glucose contents gradually increased with elevation of WS up to 35‰; then, their contents remarkably decreased. The relative expression of insulin like growth factor-1 in the liver of fish reared at 35‰ was strikingly higher than that in the other groups. The relative expression of HSP70 gene in fresh water group was pronouncedly elevated compared to other treatments. The relative expression of interleukin-1β in 15 and 35‰ groups was higher than that in the other groups; however, the relative expression of lysozyme gene in the liver of fish reared at fresh water was pronouncedly lower than that in the other treatments. The results of this study suggested rearing L. calcarifer at 15‰ closer to the isosmotic point and better provide its welfare.

Gill Transcriptome Analysis Revealed the Difference in Gene Expression Between Freshwater and Seawater Acclimated Guppy (Poecilia reticulata)

Guppy (Poecilia reticulata) can adapt to a wide range of salinity changes. To investigate the gene expression changes in the guppy exposed to seawater, we characterized its gill transcriptome using RNA sequencing. Experimental fish were exposed to salinity increase from 0 to 30‰ within 4 days, while control fish were cultured in freshwater (0‰ salinity). Seven days after salinity exposure, the gills were sampled and the mortality within 2 weeks was recorded. No significant difference in the cumulative mortality at the second week was found between the two groups. Transcriptomic analysis identified 3477 differentially expressed genes (DEGs), including 1067 upregulated and 2410 downregulated genes. These DEGs were enriched in several biological processes, including ion transport, ion homeostasis, ATP biosynthetic process, metabolic process, and immune system process. Oxidative phosphorylation was the most activated pathway. DEGs involved in the pathway "endoplasmic reticulum (ER)-mediated phagocytosis," "starch and sucrose metabolism," and "steroid biosynthesis" were mainly downregulated; chemokines and interleukins involved in "cytokine-cytokine receptor interaction" were differentially expressed. The present results suggested that oxidative phosphorylation had essential roles in osmoregulation in the gills of seawater acclimated guppy, during which the decline in the expression of genes encoding V-ATPases and calreticulin had a negative effect on the phagocytosis and immune response. Besides, several metabolic processes including "starch and sucrose metabolism" and "steroid biosynthesis" were affected. This study elucidates transcriptomic changes in osmotic regulation, metabolism, and immunity in seawater acclimated guppy.

Environmental salinity influences the branchial expression of TCR pathway related genes based on transcriptome of a catadromous fish

Environmental salinity not only affects the physiological processes such as osmoregulation and hormonal control, but also changes the immune system in fishes. Studies are limited in fish on the roles of the T cell receptor (TCR)-related genes in relation to changes in environmental salinity. A large group of salinity-challenged transcripts was obtained in gills of marbled eel (Anguilla marmorata). Moreover, bioinformatic ways were used to identify the enriched TCR pathway related genes which were significantly different expressed in fresh water (FW), brackish water (BW) and seawater (SW). Meanwhile, the RT-qPCR results were validated and consistent with the RNA-seq results. TCR a, TCR b, CD45, CD28, PI3K, LCK and LAT were up-regulated when the salinity increases in BW and SW, which connected with the related signaling pathways (Ras-MAPK and PKC pathway). CD4 and Zap70 were down-regulated when the salinity increases in BW and SW, which connected with the PLC pathway. The research offers a novel viewpoint to explore the immune pathways including the TCR pathway in fish based on transcriptome.

Integrated analysis of transcriptomic and metabolomic data to evaluate responses to hypersalinity stress in the gill of the razor clam (Sinonovacula constricta)

Salinity is an important ecological factor that affects physiological metabolism, survival, and distribution of marine organisms. Despite changes in the osmolarity and composition of the cytosol during salinity shifts, marine mollusks are able to maintain their metabolic function. The razor clam (Sinonovacula constricta) survives the wide range of salinity in the intertidal zone via changes in behavior and physiology. To explore the stress responses and mechanisms of salinity tolerance in razor clams, we collected transcriptomic and metabolomic data from a control group (salinity 20‰, S20) and a salinity-stress group (salinity 35‰, S35). The transcriptome data showed that genes related to the immune system, cytoskeleton remodeling, and signal transduction pathways dominated in the S35 group to counteract hypersalinity stress in the gill. The metabolomic analysis showed that 142 metabolites were significantly different between the S35 and S20 groups and that amino acid and carbohydrate metabolism were affected by hypersalinity stress. Levels of amino acids and energy substances, such as l-proline, isoleucine, and fructose, were higher in the gill of the S35 group. The combination of transcriptomic and metabolomic data indicated that metabolism of amino acids, carbohydrates, and lipids was enhanced in the gill during adaptation to high salinity. These results clarified the complex physiological processes involved in the response to hyperosmotic stress and maintenance of metabolism in the gill of razor clams. These findings provide a reference for further study of the biological responses of euryhaline shellfish to hyperosmotic stress and a molecular basis for the search for populations with high salinity tolerance.

Identification and Expression Analysis of Long Non-coding RNA in Large Yellow Croaker (Larimichthys crocea) in Response to Cryptocaryon irritans Infection

Large-scale transcription studies have revealed numerous lncRNAs (long non-coding RNAs). lncRNAs have been proposed to participate in the regulation of a diverse range of biological processes, including transcriptional regulation. Although lncRNAs have attracted increasing attention, the studies in large yellow croaker (Larimichthys crocea) are still rare, and they lack systematic analysis. In this study, 101 RNA-seq datasets varied in ages, sexes, and tissues were retrieved from the NCBI database to generate a comprehensive catalog of large yellow croaker transcriptome database. A set of 14,599 high-confidence lncRNAs from 13,673 loci were identified and characterized. Furthermore, RNA-seq datasets obtained from the infection of C. irritans were employed to investigate the differential expression pattern of lncRNAs and analyze potential biological functions. A total of 77 differentially expressed lncRNAs targeting to 567 protein-coding genes were identified by using expression analysis. Several immune genes, including TLR5, CD2AP, and MMP9, were highlighted. With GO enrichment and KEGG pathway analysis, the immune-related terms or pathways were enriched. This study created a comprehensive dataset of lncRNAs for large yellow croaker, which would be helpful for the researches of functional roles of lncRNAs in large yellow croaker.

Gill transcriptomes reveal expression changes of genes related with immune and ion transport under salinity stress in silvery pomfret (Pampus argenteus)

Salinity is a major ecological factor in the marine environment, and extremely important for the survival, development, and growth of fish. In this study, gill transcriptomes were examined by high-throughput sequencing at three different salinities (12 ppt as low salinity, 22 ppt as control salinity, and 32 ppt as high salinity) in an importantly economical fish silvery pomfret. A total of 187 genes were differentially expressed, including 111 up-regulated and 76 down-regulated transcripts in low-salinity treatment group and 107 genes differentially expressed, including 74 up-regulated and 33 down-regulated transcripts in high-salinity treatment group compared with the control group, respectively. Some pathways including NOD-like receptor signaling pathway, cytokine-cytokine receptor interaction, Toll-like receptor pathway, cardiac muscle contraction, and vascular smooth muscle contraction were significantly enriched. qPCR analysis further confirmed that mRNA expression levels of immune (HSP90A, IL-1β, TNFα, TLR2, IP-10, MIG, CCL19, and IL-11) and ion transport-related genes (WNK2, NPY2R, CFTR, and SLC4A2) significantly changed under salinity stress. Low salinity stress caused more intensive expression changes of immune-related genes than high salinity. These results imply that salinity stress may affect immune function in addition to regulating osmotic pressure in silvery pomfret.

Interactive effects of salinity variation and exposure to ZnO nanoparticles on the innate immune system of a sentinel marine bivalve, Mytilus edulis

ZnO nanoparticles (nZnO) are released into the coastal environment from multiple sources, yet their toxicity to marine organisms is not well understood. We investigated the interactive effects of salinity (normal 15, low 5, and fluctuating 5-15) and nZnO (100 μg l-1) on innate immunity of the blue mussels Mytilus edulis from a brackish area of the Baltic Sea. Exposure to ionic Zn (100 μg l-1) was used to test whether the toxic effects of nZnO can be attributed to the potential release of Zn2+. Functional parameters and the expression of key immune-related genes were investigated in the mussels exposed to nZnO or ionic Zn under different salinity regimes for 21 days. nZnO exposures elevated hemocyte mortality, suppressed adhesion, stimulated phagocytosis, and led to an apparent increase in lysosomal volume. At salinity 15, nZnO suppressed the mRNA expression of the Toll-like receptors TLRb and c, C-lectin, and the complement system component C3q indicating impaired ability for pathogen recognition. In contrast, the mRNA levels of an antimicrobial peptide defensin increased during nZnO exposure at salinity 15. At fluctuating salinity (5-15), nZnO exposure increased expression of multiple immune-related genes in hemocytes including the complement system components C1 and C3q, and the Toll-like receptors TLRa, b and c. Low salinity (5) had strong immunosuppressive effects on the functional and molecular immune traits of M. edulis that overshadowed the effects of nZnO. The salinity-dependent modulation of immune response to nZnO cannot be attributed to the differences in the aggregation or solubility of nZnO, and likely reflects the interaction of the toxic effects of nanoparticles and physiological effects of the osmotic stress. These findings have implications for the environmental risk assessment of nanomaterials and the development of the context-specific biomarker baselines for coastal pollution monitoring.

Identification of differentially expressed genes in gills of tiger puffer (Takifugu rubripes) in response to low-salinity stress

Salinity is an important abiotic factor for aquatic organisms. In fish, changes in salinity affect physiological responses and alter the immune system. Takifugu rubripes is an important economic marine fish, and mechanisms of T. rubripes adaptation to salinity changes need to be further documented. In this study, a transcriptome sequencing technique was used to analyse genes that were differentially expressed in the T. rubripes gill after low-salinity stress for 30 d, and differential gene expression was further validated by quantitative real-time PCR (qPCR). After assembly, 385 differentially expressed genes (DEGs) were identified, including 182 upregulated genes and 203 downregulated genes. The DEGs were assigned to Gene Ontology (GO) classes with a total of 1647 functional terms. Most DEGs were assigned to biological process (984; 59.8%) followed by molecular function (445; 27.0%) and cellular component (218; 13.2%). Further KEGG analysis allocated 385 DEGs to 95 KEGG pathways. After q-value correction, 7 pathways (Glycolysis/Gluconeogenesis; Biosynthesis of amino acids; Carbon metabolism; Fructose and mannose metabolism; Pentose phosphate pathway; Metabolism of xenobiotics by cytochrome P450; and Glycine, serine and threonine metabolism) remained significant. qPCR results indicated that the transcripts of six selected genes sharply increased after 30 d of low-salinity stress. Low-salinity stress obviously increased SLC39A6, SLC5A9, NKAα1, CYP1A1, CYP1B1, and GSTA expression. In contrast, the genes encoding Aldoaa, GPI, FBP2 and GAPDH exhibited downregulation. In addition, three solute carrier (SLC) genes selected from the DEGs were further studied for differential expression patterns after low-salinity exposure, and the results showed that the SLCs were upregulated in T. rubripes after 72 h of low-salinity exposure. This investigation provides data for understanding the molecular mechanisms of fish responses to low-salinity stress and provides a reference for rationally setting salinity levels in aquaculture.

Optimum salinity for Nile tilapia (Oreochromis niloticus) growth and mRNA transcripts of ion-regulation, inflammatory, stress- and immune-related genes

We aim to study the optimum salinity concentration for Nile tilapia, through the assessment of its growth performance and the expression of its related genes (Gh and Igf-1), as well as its salinity adaptation and immune status through the assessment of the gene expression of ion-regulation genes (Na⁺/K⁺-ATPase α-1a and α-1b), stress-related genes (GST, HSP27, and HSP70), inflammatory-related genes (IL1, IL8, CC, and CXC chemokine), and immune-related genes (IgMH TLR7, MHC, and MX) at the osmoregulatory organs (gills, liver, and kidney). Based on the least mortality percentage and the physical appearance of the fish, three salt concentrations (6, 16, and 20 ppt) were chosen following a 6-month preliminary study using serial salt concentrations ranged from 6 to 36 ppt, which were obtained by rearing the fish in gradual elevated pond salinity through daily addition of 0.5 ppt saline water. The fish size was 10.2-12 cm and weight was 25.5-26.15 g. No significant differences in the fish weight gain were observed among the studied groups. The group reared at 16-ppt salt showed better performance than that of 20 ppt, as they have lower morality % and higher expression of ion-regulated gene (Na⁺/K⁺-ATPase α1-b), stress-related genes (GST, HSP27, and HSP70) of the gills and also GST, inflammatory-related genes (IL-1β and IL8), and TLR in the liver tissue. Higher expression of kidney-immune-related genes at 20-ppt salt may indicate that higher salinity predispose to fish infection and increased mortality. We concluded that 16-ppt salinity concentration is suitable for rearing O. niloticus as the fish are more adaptive to salinity condition without changes in their growth rate. Also, we indicate the use of immune stimulant feed additive to overcome the immune suppressive effect of hyper-salinity. Additionally, the survival of some fish at higher salinity concentrations (30-34 ppt) increase the chance for selection for salinity resistance in the Nile tilapia.

Identification of immune-related genes in gill cells of Japanese eels (Anguilla japonica) in adaptation to water salinity changes

The changes in ambient salinity influence ion and water homeostasis, hormones secretion, and immune response in fish gills. The physiological functions of hormones and ion transporters in the regulation of gill-osmoregulation have been widely studied, however the modulation of immune response under salinity changes is not determined. Using transcriptome sequencing, we obtained a comprehensive profile of osmo-responsive genes in gill cells of Japanese eel (Anguilla japonica). Herein, we applied bioinformatics analysis to identify the immune-related genes that were significantly higher expressed in gill pavement cells (PVCs) and mitochondrial-rich cells (MRCs) in freshwater (FW) than seawater (SW) adapted fish. We validated the data using the real-time qPCR, which showed a high correlation between the RNA-seq and real-time qPCR data. In addition, the immunohistochemistry results confirmed the changes of the expression of selected immune-related genes, including C-reactive protein (CRP) in PVCs, toll-like receptor 2 (TLR2) in MRCs and interleukin-1 receptor type 2 (IL-1R2) in both PVCs and MRCs. Collectively our results demonstrated that those immune-related genes respond to salinity changes, and might trigger related special signaling pathways and network. This study provides new insights into the impacts of ambient salinity changes on adaptive immune response in fish gill cells.

Comparative and evolutionary analysis of an adapter molecule MyD88 in invertebrate metazoans

The myeloid differentiation factor 88 (MyD88) is an essential adapter in Toll-like receptor (TLR) signalling pathways, with TLR the first pattern-recognition receptor (PRR) that was discovered in Drosophila. In the present study, a MyD88 gene was identified and characterized from a commercially important shellfish, Scapharca subcrenata, including a DEATH domain and TIR domain conserved within other molluscs. Furthermore, comparative genomic evidence revealed that MyD88 was of different lengths and contained quantitative exon and intron regions, which might be involved in specific mechanisms. To further explore the phylogenetic relationships of invertebrate metazoan MyD88, we applied MrBayes and PhyML software to construct phylogenetic trees using Bayesian and maximum likelihood approaches, respectively, which suggested that the MyD88 of Arthropoda was closely related to lower invertebrates, in contrast to morphological taxonomy. Finally, we investigated the evolutionary patterns and location of positive selection sites (PSSs) in the MyD88 gene from Arthropoda, Mollusca and Insecta using PAML software with the maximum likelihood method. The data showed that positive selection sites were detected in these groups, and partial sites were located in the TIR domain but were not found in the DEATH domain. To summarize, in this study, we report on the diversification of MyD88 in invertebrate metazoans, the specific evolutionary position of Arthropoda MyD88, and the positive selection pressures on MyD88 of Arthropoda, Mollusca and Insecta. These results are a valuable contribution to understand and clarify the evolutionary pattern of TLR/MyD88 signalling pathways in invertebrate and vertebrate taxa.