Environmental stressors, such as temperature, are relevant factors that could generate a negative effect on several tissues in fish. A key fish species for Chilean aquaculture diversification is the red cusk-eel (Genypterus chilensis), a native fish for which knowledge on environmental stressors effects is limited. This study evaluated the effects of high-temperature stress on the liver of red cusk-eel in control (14 °C) and high-temperature (19 °C) groups using multiple approaches: determination of plasmatic hepatic enzymes (ALT, AST, and AP), oxidative damage evaluation (AP sites, lipid peroxidation, and carbonylated proteins), and RNA-seq analysis. High-temperature stress generated a significant increase in hepatic enzyme activity in plasma. In the liver, a transcriptional regulation was observed, with 1239 down-regulated and 1339 up-regulated transcripts. Additionally, high-temperature stress generated oxidative stress in the liver, with oxidative damage and transcriptional modulation of the antioxidant response. Furthermore, an unfolded protein response was observed, with several pathways enriched, as well as a heat shock response, with several heat shock proteins up regulated, suggesting candidate biomarkers (i.e., serpinh1) for thermal stress evaluation in this species. The present study shows that high-temperature stress generated a major effect on the liver of red cusk-eel, knowledge to consider for the aquaculture and fisheries of this species. Full article

In aquaculture, the Chinese soft-shelled turtle (Pelodiscus sinensis) is an economically important species with remarkable gender dimorphism in its growth patterns. However, the underlying molecular mechanisms of this phenomenon have not been elucidated well. Here, we conducted a whole-transcriptome analysis of the female and male gonads of P. sinensis. Overall, 7833 DE mRNAs, 619 DE lncRNAs, 231 DE circRNAs, and 520 DE miRNAs were identified. Some “star genes” associated with sex differentiation containing dmrt1, sox9, and foxl2 were identified. Additionally, some potential genes linked to sex differentiation, such as bmp2, ran, and sox3, were also isolated in P. sinensis. Functional analysis showed that the DE miRNAs and DE ncRNAs were enriched in the pathways related to sex differentiation, including ovarian steroidogenesis, the hippo signaling pathway, and the calcium signaling pathway. Remarkably, a lncRNA/circRNA–miRNA–mRNA interaction network was constructed, containing the key genes associated with sex differentiation, including fgf9, foxl3, and dmrta2. Collectively, we constructed a gender dimorphism profile of the female and male gonads of P. sinensis, profoundly contributing to the exploration of the major genes and potential ncRNAs involved in the sex differentiation of P. sinensis. More importantly, we highlighted the potential functions of ncRNAs for gene regulation during sex differentiation in P. sinensis as well as in other turtles. Full article

Stress-inducing husbandry and rearing conditions, bacterial infections or parasitic diseases may all lead to chronic inflammation. The immune response will then channel energy away from growth, reproduction and other important physiological processes, to fuel immune-related metabolic responses. The present study aims to unravel the mechanisms and contribute with new information on the molecular, cellular and humoral parameters of European seabass (Dicentrarchus labrax) undergoing chronic inflammation that can be used as health indicators for application in fish health management. European seabass individuals were intra-peritoneally injected with either Freund’s Incomplete Adjuvant (FIA) to induce inflammation or Hanks Balanced Salt Solution (HBSS) to serve as sham. Fish were sampled at 24 h, 7, 14 and 21 days post-injection and blood, plasma and head-kidney were collected. The results found were clear indicators of an inflamed peritoneal cavity and an ongoing systemic immune response that persisted for at least 21 days. Locally, inflammation was characterized by an intense recruitment of immune cells that was still evident 21 days after injection, thus illustrating the chronic character of the immune response. Cellular response was also noticed peripherally with leukocyte numbers rising in the blood of FIA-injected fish. Furthermore, the cellular-mediated respiratory burst peaked at 21 days post-FIA injection, suggesting that phagocytes were still actively fighting the phlogistic agent. Regarding the head-kidney molecular analysis, cxcr4 and il34 appear to be good markers of a chronic inflammation response due to their importance for pathways with high relevance in chronic inflammation settings. In addition, opioid receptor nopr seems to be a good marker of a chronic inflammation response due to its role in detecting noxious stimuli. The present study can serve as a baseline to assess long-term immune-related responses in future studies. For that, more research is nonetheless required to select more responsive and specific molecular markers. Full article

Optimal smoltification is crucial for normal development, growth, and health of farmed Atlantic salmon in seawater. Here, we characterize miRNA expression in liver to reveal whether miRNAs regulate gene expression during this developmental transition. Expression changes of miRNAs and mRNAs was studied by small-RNA sequencing and microarray analysis, respectively. This revealed 62 differentially expressed guide miRNAs (gDE-miRNAs) that could be divided into three groups with characteristic dynamic expression patterns. Three of miRNA families are known as highly expressed in liver. A rare arm shift was observed during smoltification in the Atlantic salmon-specific novel-ssa-miR-16. The gDE-miRNAs were predicted to target 2804 of the genes revealing expression changes in the microarray analysis. Enrichment analysis revealed that targets were significantly enriched in smoltification-associated biological process groups. These included lipid and cholesterol synthesis, carbohydrate metabolism, protein metabolism and protein transport, immune system genes, circadian rhythm and stress response. The results indicate that gDE-miRNAs may regulate many of the changes associated with this developmental transition in liver. The results pave the way for validation of the predicted target genes and further study of gDE-miRNA and their targets by functional assays. Full article

The allotetraploid common carp is one of the most important freshwater food fish. However, the IBs found in allotetraploid common carp increase the difficulty in fish meat processing and consumption. Although candidate genes associated with the total IB number have been identified, the SNPs associated with the numbers of the total IBs and different forms of IBs have not yet been identified, hindering the breeding of IB-reduced common carp. Herein, the numbers of different types of IBs in three common carp strains were measured. Using whole-genome resequencing and bulked segregant analysis in three pairs of IB-more and IB-less groups, we identified the consensus nonsynonymous SNPs in three strains of common carp. Screening the flanking regions of these SNPs led to the detection of other SNPs. Association study detected 21 SNPs significantly associated with the number of total IBs, epineural-IBs, and ten detailed types of IBs. We observed the joint effects of multiple SNPs on each associated IB number with an improved explained percentage of phenotypic variation. The resulting dataset provides a resource to understand the molecular mechanisms of IB development in different common carp strains. These SNPs are potential markers for future selection to generate IB-reduced common carp. Full article

The allo-tetraploid common carp, one widely cultured food fish, is able to produce poly-unsaturated fatty acids (PUFAs). The genetic markers on the PUFA contents for breeding was limited. The polymorphisms in elovl5a and elovl5b, the rate-limiting enzymes in the PUFA biosynthesis, have not been investigated yet. Herein, we identified one coding SNP (cSNP) in elovl5a associated with the content of one PUFA and two cSNPs in elovl5b with the contents of eight PUFAs. The heterozygous genotypes in these three loci were associated with higher contents than the homozygotes. Together with previously identified two associated cSNPs in fads2b, we found the joint effect of these four cSNPs in fads2b and elovl5b on the PUFA contents with the increased explained percentages of PUFA contents. The genotype combinations of more heterozygotes were associated with higher PUFA contents than the other combinations. Using ten genomic selection programs with all cSNPs in fads2b and elovl5b, we obtained the high and positive correlations between the phenotypes and the estimated breeding values of eight PUFAs. These results suggested that elovl5b might be the major gene corresponding to common carp PUFA contents compared with elovl5a. The cSNP combinations in fads2b and elovl5b and the optimal genomic selection program will be used in the future selection breeding to improve the PUFA contents of common carp. Full article

In fish, as observed in mammals, any stressful event affects the immune system to a larger or shorter extent. The neuroendocrine-immune axis is a bi-directional network of mobile compounds and their receptors that are shared between both systems (neuroendocrine and immune) and that regulate their respective responses. However, how and to what extent immunity modulates the neuroendocrine system is not yet fully elucidated. This study was carried out to understand better central gene expression response patterns in a high-valued farmed fish species to an acute peripheral inflammation, focusing on genes related to the hypothalamus-pituitary-interrenal axis and the opioid system. European seabass, Dicentrarchus labrax, were intra-peritoneally injected with either Freund’s Incomplete Adjuvant to induce a local inflammatory response or Hanks Balances Salt Solution to serve as the control. An undisturbed group was also included to take into account the effects due to handling procedures. To evaluate the outcomes of an acute immune response, fish were sampled at 4, 24, 48, and 72 h post-injection. The brain was sampled and dissected for isolation of different regions: telencephalon, optic tectum, hypothalamus, and pituitary gland. The expression of several genes related to the neuroendocrine response was measured by real-time PCR. Data were statistically analyzed by ANOVA and discriminant analyses to obtain these genes’ responsiveness for the different brain regions. Serotonergic receptors were upregulated in the telencephalon, whereas the optic tectum inhibited these transcription genes. The hypothalamus showed a somewhat delayed response in which serotonin and glucocorticoid receptors were concerned. Still, the hypothalamic corticotropin-releasing hormone played an important role in differentiating fish undergoing an inflammatory response from those not under such conditions. Opioid receptors gene expression increased in both the hypothalamus and the telencephalon, while in the optic tectum, most were downregulated. However, no changes in the pituitary gland were observed. The different brain regions under immune stimulation demonstrated clear, distinct responses regarding gene transcription rates as well as the time period needed for the effect to occur. Further, more integrative studies are required to associate functions to the evaluated genes more safely and better understand the triggering mechanisms. Full article

The heat shock protein (Hsp) superfamily has received accumulated attention because it is ubiquitous and conserved in almost all living organisms and is involved in a wide spectrum of cellular responses against diverse environmental stresses. However, our knowledge about the Hsp co-chaperon network is still limited in non-model organisms. In this study, we provided the systematic analysis of 95 Hsp genes (LmHsps) in the genome of spotted sea bass (Lateolabrax maculatus), an important aquaculture species in China that can widely adapt to diverse salinities from fresh to sea water, and moderately adapt to high alkaline water. Through in silico analysis using transcriptome and genome database, we determined the expression profiles of LmHsps in response to salinity change and alkalinity stress in L. maculatus gills. The results revealed that LmHsps were sensitive in response to alkalinity stress, and the LmHsp40-70-90 members were more actively regulated than other LmHsps and may also be coordinately interacted as co-chaperons. This was in accordance with the fact that members of LmHsp40, LmHsp70, and LmHsp90 evolved more rapidly in L. maculatus than other teleost lineages with positively selected sites detected in their functional domains. Our results revealed the diverse and cooperated regulation of LmHsps under alkaline stress, which may have arisen through the functional divergence and adaptive recruitment of the Hsp40-70-90 co-chaperons and will provide vital insights for the development of L. maculatus cultivation in alkaline water. Full article

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies. Full article

The hybrid between female channel catfish (Ictalurus punctatus) and male blue catfish (Ictalurus furcatus) is superior in feed conversion, disease resistance, carcass yield, and harvestability compared to both parental species. However, heterosis and heterobeltiosis only occur in pond culture, and channel catfish grow much faster than the other genetic types in small culture units. This environment-dependent heterosis is intriguing, but the underlying genetic mechanisms are not well understood. In this study, phenotypic characterization and transcriptomic analyses were performed in the channel catfish, blue catfish, and their reciprocal F₁s reared in tanks. The results showed that the channel catfish is superior in growth-related morphometrics, presumably due to significantly lower innate immune function, as investigated by reduced lysozyme activity and alternative complement activity. RNA-seq analysis revealed that genes involved in fatty acid metabolism/transport are significantly upregulated in channel catfish compared to blue catfish and hybrids, which also contributes to the growth phenotype. Interestingly, hybrids have a 40–80% elevation in blood glucose than the parental species, which can be explained by a phenomenon called transgressive expression (overexpression/underexpression in F₁s than the parental species). A total of 1140 transgressive genes were identified in F₁ hybrids, indicating that 8.5% of the transcriptome displayed transgressive expression. Transgressive genes upregulated in F₁s are enriched for glycan degradation function, directly related to the increase in blood glucose level. This study is the first to explore molecular mechanisms of environment-dependent heterosis/heterobeltiosis in a vertebrate species and sheds light on the regulation and evolution of heterosis vs. hybrid incompatibility. Full article

The Chinese soft-shelled turtle Pelodiscus sinensis shows obvious sexual dimorphism. The economic and nutrition value of male individuals are significantly higher than those of female individuals. Pseudo-females which are base to all-male breeding have been obtained by estrogen induction, while the gene function and molecular mechanism of sex reversal remain unclear in P. sinensis. Here, comparative transcriptome analyses of female, male, and pseudo-female gonads were performed, and 14,430 genes differentially expressed were identified in the pairwise comparison of three groups. GO and KEGG analyses were performed on the differentially expressed genes (DEGs), which mainly concentrated on steroid hormone synthesis. Furthermore, the results of gonadal transcriptome analysis revealed that 10 sex-related sox genes were differentially expressed in males vs. female, male vs. pseudo-female, and female vs. pseudo-female. Through the differential expression analysis of these 10 sox genes in mature gonads, six sox genes related to sex reversal were further screened. The molecular mechanism of the six sox genes in the embryo were analyzed during sex reversal after E₂ treatment. In mature gonads, some sox family genes, such as sox9sox12, and sox30 were highly expressed in the testis, while sox1, sox3, sox6, sox11, and sox17 were lowly expressed. In the male embryos, exogenous estrogen can activate the expression of sox3 and inhibit the expression of sox8, sox9, and sox11. In summary, sox3 may have a role in the process of sex reversal from male to pseudo-female, when sox8 and sox9 are inhibited. Sox family genes affect both female and male pathways in the process of sex reversal, which provides a new insight for the all-male breeding of the Chinese soft-shelled turtle. Full article

The growing amount of genome information and transcriptomes data available allows for a better understanding of biological processes. However, analysis of complex transcriptomic experimental designs involving different conditions, tissues, or times is relevant. This study proposes a novel approach to analyze complex data sets combining transcriptomes and miRNAs at the chromosome-level genome. Atlantic salmon smolts were transferred to seawater under two strategies: (i) fish group exposed to gradual salinity changes (GSC) and (ii) fish group exposed to a salinity shock (SS). Gills, intestine, and head kidney samples were used for total RNA extraction, followed by mRNA and small RNA illumina sequencing. Different expression patterns among the tissues and treatments were observed through a whole-genome transcriptomic approach. Chromosome regions highly expressed between experimental conditions included a great abundance of transposable elements. In addition, differential expression analysis showed a greater number of transcripts modulated in response to SS in gills and head kidney. miRNA expression analysis suggested a small number of miRNAs involved in the smoltification process. However, target analysis of these miRNAs showed a regulatory role in growth, stress response, and immunity. This study is the first to evidence the interplaying among mRNAs and miRNAs and the structural relationship at the genome level during Atlantic salmon smoltification. Full article

Intermuscular bones (IBs) are small spicule-like bones located in the myosepta of basal teleosts, which negatively affect the edibleness and economic value of fish. Blunt snout bream (Megalobrama amblycephala, with epineural and epipleural IBs) and tilapia (Oreochromis niloticus, without epineural and epipleural IBs) are two major aquaculture species and ideal models for studying the formation mechanisms of fish IBs. Here, we compared myosepta development between M. amblycephala and O. niloticus, based on histological analysis, transcriptome profiling, and expression analysis of bone-related genes. The histological results showed that dye condensation began to appear in the myosepta 20 days post hatching (dph) in M. amblycephala, and IBs could be clearly observed 50 dph in the myosepta, based on different staining methods. However, in O. niloticus, dye condensation was not observed in the myosepta from 10 to 60 dph. Differentially expressed genes (DEGs) at different developmental stages were screened by comparing the transcriptomes of M. amblycephala and O. niloticus, and KEGG analysis demonstrated that these DEGs were enriched in many bone-related pathways, such as focal adhesion, calcium, and Wnt signaling pathways. Quantitative PCR was performed to further compare the expression levels of some bone-related genes (scxa, scxb, runx2a, runx2b, bgp, sp7, col1a2, entpd5a, entpd5b, phex, alpl, and fgf23). All the tested genes (except for alpl) exhibited higher expression levels in M. amblycephala than in O. niloticus. A comparison of the dorsal and abdominal muscle tissues between the two species also revealed significant expression differences for most of the tested genes. The results suggest that scxa, scxb, runx2a, runx2b, entpd5a, col1a2, and bgp may play important roles in IB development. Our findings provide some insights into the molecular mechanisms of IB formation, as well as clues for further functional analysis of the identified genes to better understand the development of IBs. Full article

Metamorphosis is a captivating process of change during which the morphology of the larva is completely reshaped to face the new challenges of adult life. In the case of fish, this process initiated in the brain has traditionally been considered to be a critical rearing point and despite the pioneering molecular work carried out in other flatfishes, the underlying molecular basis is still relatively poorly characterized. Turbot brain transcriptome of three developmental stages (pre-metamorphic, climax of metamorphosis and post-metamorphic) were analyzed to study the gene expression dynamics throughout the metamorphic process. A total of 1570 genes were differentially expressed in the three developmental stages and we found a specific pattern of gene expression at each stage. Unexpectedly, at the climax stage of metamorphosis, we found highly expressed genes related to the immune response, while the biological pathway enrichment analysis in pre-metamorphic and post-metamorphic were related to cell differentiation and oxygen carrier activity, respectively. In addition, our results confirm the importance of thyroid stimulating hormone, increasing its expression during metamorphosis. Based on our findings, we assume that immune system activation during the climax of metamorphosis stage could be related to processes of larval tissue inflammation, resorption and replacement, as occurs in other vertebrates. Full article

Tripartite motif proteins (TRIMs), especially B30.2 domain-containing TRIMs (TRIMs-B30.2), are increasingly well known for their antiviral immune functions in mammals, while antiviral TRIMs are far from being identified in teleosts. In the present study, we identified a total of 42 CiTRIMs from the genome of grass carp, Ctenopharyngodon idella, an important cultured teleost in China, based on hmmsearch and SMART analysis. Among these CiTRIMs, the gene loci of 37 CiTRIMs were located on different chromosomes and shared gene collinearities with homologous counterparts from human and zebrafish genomes. They possessed intact conserved RBCC or RB domain assemblies at their N-termini and eight different domains, including the B30.2 domain, at their C-termini. A total of 19 TRIMs-B30.2 were identified, and most of them were clustered into a large branch of CiTRIMs in the dendrogram. Tissue expression analysis showed that 42 CiTRIMs were universally expressed in various grass carp tissues. A total of 11 significantly differentially expressed CiTRIMs were found in two sets of grass carp transcriptomes during grass carp reovirus (GCRV) infection. Three of them, including Cibtr40, CiTRIM103 and CiTRIM109, which all belonged to TRIMs-B30.2, were associated with the type I interferon response during GCRV infection by weighted network co-expression and gene expression trend analyses, suggesting their involvement in antiviral immunity. These findings may offer useful information for understanding the structure, evolution, and function of TRIMs in teleosts and provide potential antiviral immune molecule markers for grass carp. Full article

Although hemorrhagic sepsis caused by Aeromonas hydrophila infection is the dominant disease in the aquaculture of Chinese soft-shelled turtle, information on its molecular pathology is seriously limited. In this study, ninety turtles intraperitoneally injected with A. hydrophila exhibited two different phenotypes based on the pathological symptoms, referred to as active and inactive turtles. Comparative transcriptomes of liver and spleen from these two groups at 6, 24, and 72 h post-injection (hpi) were further analyzed. The results showed that cytokine–cytokine receptor interaction, PRRs mediated signaling pathway, apoptosis, and phagocytosis enriched in active and inactive turtles were significantly different. Pro-inflammatory cytokines, the TLR signaling pathway, NLR signaling pathway, and RLR signaling pathway mediating cytokine expression, and apoptosis-related genes, were significantly up-regulated in inactive turtles at the early stage (6 hpi). The significant up-regulation of phagocytosis-related genes occurred at 24 hpi in inactive turtles and relatively lagged behind those in active turtles. The anti-inflammatory cytokine, IL10, was significantly up-regulated during the tested periods (6, 24, and 72 hpi) in active turtles. These findings offer valuable information for the understanding of molecular immunopathogenesis after A. hydrophila infection, and facilitate further investigations on strategies against hemorrhagic sepsis in Chinese soft-shelled turtle T. sinensis. Full article

Edible red sea urchin (Loxechinus albus) is an endemic echinoderm species of the Chilean coasts. The worldwide demand for high-quality gonads of this species has addressed the depletion of its natural populations. Studies on this sea urchin are limited, and genomic information is almost nonexistent. Hence, generate a transcriptome is crucial information that will considerably enrich molecular data and promote future findings for the L. albus aquaculture. Here, we obtained transcriptomic data of the edible red sea urchin by Illumina platform. Total RNA was extracted from gonads, intestines, and coelomocytes of juvenile urchins, and samples were sequenced using MiSeq Illumina technology. A total of 91,119,300 paired-end reads were de novo assembled, 185,239 transcripts produced, and a reference transcriptome created with 38.8% GC content and an N50 of 1769 bp. Gene ontology analysis revealed notable differences in the expression profiles between gonads, intestines, and coelomocytes, allowing the detection of transcripts associated with specific biological processes and KEGG pathways. These data were validated using 12 candidate transcripts by real-time qPCR. This dataset will provide a valuable molecular resource for L. albus and other species of sea urchins. Full article

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms. Full article

The importance of dietary omega-6 to omega-3 (ω6:ω3) fatty acid (FA) ratios for human health has been extensively examined. However, its impact on fish physiology, and the underlying molecular mechanisms, are less well understood. This study investigated the influence of plant-based diets (12-week exposure) with varying ω6:ω3 (0.4–2.7) on the hepatic transcriptome of Atlantic salmon. Using 44 K microarray analysis, genes involved in immune and inflammatory response (lect2a, itgb5, helz2a, p43), lipid metabolism (helz2a), cell proliferation (htra1b), control of muscle and neuronal development (mef2d) and translation (eif2a, eif4b1, p43) were identified; these were differentially expressed between the two extreme ω6:ω3 dietary treatments (high ω6 vs. high ω3) at week 12. Eight out of 10 microarray-identified transcripts showed an agreement in the direction of expression fold-change between the microarray and qPCR studies. The PPARα activation-related transcript helz2a was confirmed by qPCR to be down-regulated by high ω6 diet compared with high ω3 diet. The transcript expression of two helz2 paralogues was positively correlated with ω3, and negatively with ω6 FA in both liver and muscle, thus indicating their potential as biomarkers of tissue ω6:ω3 variation. Mef2d expression in liver was suppressed in the high ω6 compared to the balanced diet (ω6:ω3 of 2.7 and 0.9, respectively) fed fish, and showed negative correlations with ω6:ω3 in both tissues. The hepatic expression of two lect2 paralogues was negatively correlated with viscerosomatic index, while htra1b correlated negatively with salmon weight gain and condition factor. Finally, p43 and eif2a were positively correlated with liver Σω3, while these transcripts and eif4b2 showed negative correlations with 18:2ω6 in the liver. This suggested that some aspects of protein synthesis were influenced by dietary ω6:ω3. In summary, this nutrigenomic study identified hepatic transcripts responsive to dietary variation in ω6:ω3, and relationships of transcript expression with tissue (liver, muscle) lipid composition and other phenotypic traits. Full article

Single-cell sequencing technologies capable of providing us with immune information from dozens to thousands of individual cells simultaneously have revolutionized the field of immunology these past years. However, to date, most of these novel technologies have not been broadly applied to non-model organisms such as teleost fish. In this study, we used the 10× Genomics single cell RNA sequencing technology and used it to analyze for the first time in teleost fish the transcriptional pattern of single B cells from peripheral blood. The analysis of the data obtained in rainbow trout revealed ten distinct cell clusters that seem to be associated with different subsets and/or maturation/differentiation stages of circulating B cells. The potential characteristics and functions of these different B cell subpopulations are discussed on the basis of their transcriptomic profile. The results obtained provide us with valuable information to understand the biology of teleost B cells and offer us a repertoire of potential markers that could be used in the future to differentiate trout B cell subsets. Full article

Hepatopancreas necrosis disease (HPND) of the Chinese mitten crab Eriocheir sinensis causes huge economic loss in China. However, the pathogenic factors and pathogenesis are still a matter of dissension. To search for potential pathogens, the hepatopancreatic flora of diseased crabs with mild symptoms, diseased crabs with severe symptoms, and crabs without visible symptoms were investigated using metatranscriptomics sequencing. The prevalence of Absidia glauca and Candidatus Synechococcus spongiarum decreased, whereas the prevalence of Spiroplasma eriocheiris increased in the hepatopancreatic flora of crabs with HPND. Homologous sequences of 34 viral species and 4 Microsporidian species were found in the crab hepatopancreas without any significant differences between crabs with and without HPND. Moreover, DEGs in the hepatopancreatic flora between crabs with severe symptoms and without visible symptoms were enriched in the ribosome, retinol metabolism, metabolism of xenobiotics by cytochrome P450, drug metabolism—cytochrome P450, biosynthesis of unsaturated fatty acids, and other glycan degradation. Moreover, the relative abundance of functions of DEDs in the hepatopancreatic flora changed with the pathogenesis process. These results suggested that imbalance of hepatopancreatic flora was associated with crab HPND. The identified DEGs were perhaps involved in the pathological mechanism of HPND; nonetheless, HPND did not occur due to virus or microsporidia infection. Full article

On-growing juveniles of gilthead sea bream were acclimated for 45 days to mild-hypoxia (M-HYP, 40–60% O₂ saturation), whereas normoxic fish (85–90% O₂ saturation) constituted two different groups, depending on if they were fed to visual satiety (control fish) or pair-fed to M-HYP fish. Following the hypoxia conditioning period, all fish were maintained in normoxia and continued to be fed until visual satiation for 3 weeks. The time course of hypoxia-induced changes was assessed by changes in blood metabolic landmarks and muscle transcriptomics before and after exhaustive exercise in a swim tunnel respirometer. In M-HYP fish, our results highlighted a higher contribution of aerobic metabolism to whole energy supply, shifting towards a higher anaerobic fitness following normoxia restoration. Despite these changes in substrate preference, M-HYP fish shared a persistent improvement in swimming performance with a higher critical speed at exercise exhaustion. The machinery of muscle contraction and protein synthesis and breakdown was also largely altered by mild-hypoxia conditioning, contributing this metabolic re-adjustment to the positive regulation of locomotion and to the catch-up growth response during the normoxia recovery period. Altogether, these results reinforce the presence of large phenotypic plasticity in gilthead sea bream, and highlights mild-hypoxia as a promising prophylactic measure to prepare these fish for predictable stressful events. Full article

Recent studies suggest that short pentraxins in fish might serve as biomarkers for not only bacterial infections, as in higher vertebrates including humans, but also for viral ones. These fish orthologs of mammalian short pentraxins are currently attracting interest because of their newly discovered antiviral activity. In the present work, the modulation of the gene expression of all zebrafish short pentraxins (CRP-like proteins, CRP1-7) was extensively analyzed by quantitative polymerase chain reaction. Initially, the tissue distribution of crp1-7 transcripts and how the transcripts varied in response to a bath infection with the spring viremia of carp virus, were determined. The expression of crp1-7 was widely distributed and generally increased after infection (mostly at 5 days post infection), except for crp1 (downregulated). Interestingly, several crp transcription levels significantly increased in skin. Further assays in mutant zebrafish of recombinant activation gene 1 (rag1) showed that all crps (except for crp2, downregulated) were already constitutively highly expressed in skin from rag1 knockouts and only increased moderately after viral infection. Similar results were obtained for most mx isoforms (a reporter gene of the interferon response), suggesting a general overcompensation of the innate immunity in the absence of the adaptive one. Full article

Chinese mitten crab (CMC) or Eriocheir sinensis is a strong osmoregulator that can keep rigorous cellular homeostasis. CMC can flourish in freshwater, as well as seawater, habitats and represents the most important species for freshwater aquaculture. Salt stress can have direct effects on several stages (e.g., reproduction, molting, growth, etc.) of the CMC life cycle. To get a better overview of the genes involved in the gills of CMC under different salinity conditions, we conducted an RNA-Seq meta-analysis on the transcriptomes of four publicly available datasets. The meta-analysis identified 405 differentially expressed transcripts (DETs), of which 40% were classified into various transporter classes, including accessory factors and primary active transporters as the major transport classes. A network analysis of the DETs revealed that adaptation to salinity is a highly regulated mechanism in which different functional modules play essential roles. To the best of our knowledge, this study is the first to conduct a transcriptome meta-analysis of gills from crab RNA-Seq datasets under salinity. Additionally, this study is also the first to focus on the differential expression of diverse transporters and channels (transportome) in CMC. Our meta-analysis opens new avenues for a better understanding of the osmoregulation mechanism and the selection of potential transporters associated with salinity change. Full article

The present study analyses the interactions between Philasterides dicentrarchi (a ciliate parasite that causes high mortalities in cultured flatfish) and the peritoneal cells of the turbot Scophthalmus maximus during an experimental infection. The transcriptomic response was evaluated in the parasites and in the fish peritoneal cells, at 1, 2 and 4 h post-infection (hpi) in turbot injected intraperitoneally (ip) with 10⁷ ciliates and at 12 and 48 hpi in turbot injected ip with 10⁵ ciliates. Numerous genes were differentially expressed (DE) in P. dicentrarchi, relative to their expression in control ciliates (0 hpi): 407 (369 were up-regulated) at 1 hpi, 769 (415 were up-regulated) at 2 hpi and 507 (119 were up-regulated) at 4 hpi. Gene ontology (GO) analysis of the DE genes showed that the most representative categories of biological processes affected at 1, 2 and 4 hpi were biosynthetic processes, catabolic processes, biogenesis, proteolysis and transmembrane transport. Twelve genes of the ABC transporter family and eight genes of the leishmanolysin family were DE at 1, 2 and 4 hpi. Most of these genes were strongly up-regulated (UR), suggesting that they are involved in P. dicentrarchi infection. A third group of UR genes included several genes related to ribosome biogenesis, DNA transcription and RNA translation. However, expression of tubulins and tubulin associated proteins, such as kinesins or dyneins, which play key roles in ciliate division and movement, was down-regulated (DR). Similarly, genes that coded for lysosomal proteins or that participate in the cell cycle mitotic control, glycolysis, the Krebs cycle and/or in the electron transport chain were also DR. The transcriptomic analysis also revealed that in contrast to many parasites, which passively evade the host immune system, P. dicentrarchi strongly stimulated turbot peritoneal cells. Many genes related to inflammation were DE in peritoneal cells at 1, 2 and 4 hpi. However, the response was much lower at 12 hpi and almost disappeared completely at 48 hpi in fish that were able to kill P. dicentrarchi during the first few hpi. The genes that were DE at 1, 2 and 4 hpi were mainly related to the apoptotic process, the immune response, the Fc-epsilon receptor signalling pathway, the innate immune response, cell adhesion, cell surface receptors, the NF-kappaB signalling pathway and the MAPK cascade. Expression of toll-like receptors 2, 5 and 13 and of several components of NF-κB, MAPK and JAK/STAT signalling pathways was UR in the turbot peritoneal cells. Genes expressing chemokines and chemokine receptors, genes involved in prostaglandin and leukotriene synthesis, prostaglandins, leukotriene receptors, proinflammatory cytokines and genes involved in apoptosis were strongly UR during the first four hours of infection. However, expression of anti-inflammatory cytokines such as Il-10 and lipoxygenases with anti-inflammatory activity (i.e., arachidonate 15-lipoxygenase) were only UR at 12 and/or 48 hpi, indicating an anti-inflammatory state in these groups of fish. In conclusion, the present study shows the regulation of several genes in P. dicentrarchi during the early stages of infection, some of which probably play important roles in this process. The infection induced a potent acute inflammatory response, and many inflammatory genes were regulated in peritoneal cells, showing that the turbot uses all the protective mechanisms it has available to prevent the entry of the parasite. Full article

The thymus is a primary lymphoid organ that plays a pivotal role in the adaptive immune system. The immunobiology of the thymus in fish is considered to be similar to that of mammals, but it is actually poorly characterized in several cultured teleost species. In particular, while investigations in human and veterinary medicine have highlighted that the thymus can be affected by different pathological conditions, little is known about its response during disease in fish. To better understand the role of the thymus under physiological and pathological conditions, we conducted a study in turbot (Scophthalmus maximus), a commercially valuable flatfish species, combining transcriptomic and histopathological analyses. The myxozoan parasite Enteromyxum scophthalmi, which represents a major challenge to turbot production, was used as a model of infection. The thymus tissues of healthy fish showed overrepresented functions related to its immunological role in T-cell development and maturation. Large differences were observed between the transcriptomes of control and severely infected fish. Evidence of inflammatory response, apoptosis modulation, and declined thymic function associated with loss of cellularity was revealed by both genomic and morphopathological analyses. This study presents the first description of the turbot thymus transcriptome and provides novel insights into the role of this organ in teleosts’ immune responses. Full article

Long noncoding RNAs (lncRNAs) are being increasingly recognised as key modulators of various biological mechanisms, including the immune response. Although investigations in teleosts are still lagging behind those conducted in mammals, current research indicates that lncRNAs play a pivotal role in the response of fish to a variety of pathogens. During the last several years, interest in lncRNAs has increased considerably, and a small but notable number of publications have reported the modulation of the lncRNA profile in some fish species after pathogen challenge. This study was the first to identify lncRNAs in the commercial species European sea bass. A total of 12,158 potential lncRNAs were detected in the head kidney and brain. We found that some lncRNAs were not common for both tissues, and these lncRNAs were located near coding genes that are primarily involved in tissue-specific processes, reflecting a degree of cellular specialisation in the synthesis of lncRNAs. Moreover, lncRNA modulation was analysed in both tissues at 24 and 72 h after infection with nodavirus. Enrichment analysis of the neighbouring coding genes of the modulated lncRNAs revealed many terms related to the immune response and viral infectivity but also related to the stress response. An integrated analysis of the lncRNAs and coding genes showed a strong correlation between the expression of the lncRNAs and their flanking coding genes. Our study represents the first systematic identification of lncRNAs in European sea bass and provides evidence regarding the involvement of these lncRNAs in the response to nodavirus. Full article