Search

NF-κB signalling is largely controlled by the family of ‘inhibitors of NF-κB’ (IκB). The relevant databases indicate that the genome of rainbow trout contains multiple gene copies coding for iκbα (nfkbia), iκbε (nfkbie), iκbδ (nkfbid), iκbζ (nfkbiz), and bcl3, but it lacks iκbβ (nfkbib) and iκbη (ankrd42). Strikingly, three nfkbia paralogs are apparently present in salmonid fish, two of which share a high sequence identity, while the third putative nfkbia gene is significantly less like its two paralogs. This particular nfkbia gene product, iκbα, clusters with the human IκBβ in a phylogenetic analysis, while the other two iκbα proteins from trout associate with their human IκBα counterpart. The transcript concentrations were significantly higher for the structurally more closely related nfkbia paralogs than for the structurally less similar paralog, suggesting that iκbβ probably has not been lost from the salmonid genomes but has been incorrectly designated as iκbα. In the present study, two gene variants coding for iκbα (nfkbia) and iκbε (nfkbie) were prominently expressed in the immune tissues and, particularly, in a cell fraction enriched with granulocytes, monocytes/macrophages, and dendritic cells from the head kidney of rainbow trout. Stimulation of salmonid CHSE-214 cells with zymosan significantly upregulated the iκbα-encoding gene while elevating the copy numbers of the inflammatory markers interleukin-1-beta and interleukin-8. Overexpression of iκbα and iκbε in CHSE-214 cells dose-dependently quenched both the basal and stimulated activity of an NF-κB promoter suggesting their involvement in immune-regulatory processes. This study provides the first functional data on iκbε—versus the well-researched iκbα factor—in a non-mammalian model species.

There are still numerous difficulties in the successful farming of pikeperch in the anthropogenic environment of various aquaculture systems, especially during early developmental steps in the hatchery. To investigate the physiological processes involved on the molecular level, we determined the basal expression patterns of 21 genes involved in stress and immune responses and early ontogenesis of pikeperch between 0 and 175 days post hatch (dph). Their transcription patterns most likely reflect the challenges of growth and feed conversion. The gene coding for apolipoprotein A (APOE) was strongly expressed at 0 dph, indicating its importance for yolk sac utilization. Genes encoding bone morphogenetic proteins 4 and 7 (BMP4,BMP7), creatine kinase M (CKM), and SRY-box transcription factor 9 (SOX9) were highly abundant during the peak phases of morphological changes and acclimatization processes at 4–18 dph. The high expression of genes coding for peroxisome proliferator-activated receptors alpha and delta (PPARA,PPARD) at 121 and 175 dph, respectively, suggests their importance during this strong growth phase of juvenile stages. As an alternative experimental model to replace further in vivo investigations of ontogenetically important processes, we initiated the first approach towards a long-lasting primary cell culture from whole pikeperch embryos. The present study provides a set of possible biomarkers to support the monitoring of pikeperch farming and provides a first basis for the establishment of a suitable cell model of this emerging aquaculture species.

Two structurally similar NF-κB-inhibitor-interacting Ras-like proteins (NKIRAS) regulate the activity of the transcription factor NF-κB and thereby control several early immune mechanisms in mammals. We identified the orthologous sequences of NKIRAS1 and NKIRAS2 from the rainbow trout Oncorhynchus mykiss. The level of sequence identity was similarly high (≥68%) between the two and in comparison to their mammalian counterparts. Strikingly, NKIRAS2 was present as four transcript variants. These variants differed only in length and in the nucleotide composition of their 5’ termini and were most likely generated by splicing along unconventional splice sites. The shortest NKIRAS2 variant was most strongly expressed in a lymphocyte-enriched population, while NKIRAS1 was most strongly expressed in cells of myeloid origin. Fluorescent-labelled NKIRAS1 and NKIRAS2 proteins from rainbow trout were detected in close association with the p65 subunit of NF-κB in the nucleus and cytoplasm of CHSE-214 cells. Subsequent reporter-gene experiments revealed that NKIRAS1 and a longer NKIRAS2 variant in rainbow trout decreased the level of activated NF-κB, while the two shortest NKIRAS2 variants increased the NF-κB activity. In addition, the overexpression of the shortest NKIRAS2 variant in CHSE-214 cells induced a stronger transcription of the genes encoding the pro-inflammatory cytokines TNF, CXCL8, and IL1B compared to non-transfected control cells. This is the first characterisation of NKIRAS orthologues in bony fish and provides additional information to the as yet underexplored inhibition pathways of NF-κB in lower vertebrates.

The functional spectrum of the teleostean head kidney covers haematopoietic, immune and endocrine signalling pathways with physiological effects that are likely to conflict if activated at the same time. An in vivo experiment on the salmonid fish maraena whitefish (Coregonus maraena) revealed that the head kidney shows a remarkably strong response after injection of Aeromonas salmonicida within 48 h. In order to investigate the potential influence of endocrine signalling on the initiation of immune responses, we established a primary culture of head-kidney cells of maraena whitefish. For the characterisation of this model system, we used flow cytometry complemented with an extensive panel of immunological/haematological and stress-physiological/neuroendocrinological qPCR assays. More than one third of the cells expressed the characteristic signature of myeloid cells, while more than half of the cells expressed those genes typical for lymphocytes and monocytes. In parallel, we quantified the expression of genes encoding endocrine receptors and identified ADRA2D as by far the most highly expressed adrenergic-receptor gene in head-kidney cells. The stimulation of the head-kidney cells with toll-like receptor ligands induced the expression of typical immune genes (IL1B, CXCL8, TNF, SAA) after only 1 h. The incubation with the stress hormones cortisol, adrenaline and noradrenaline also had an immune-activating effect, though less pronounced. However, cortisol had the strongest suppressive effect on the stimulation-induced immune response, while adrenaline exerted a comparably weaker effect and noradrenaline was almost ineffective. Moreover, we found that cortisol reduced the expression of genes coding for adrenergic and some glucocorticoid receptors, while noradrenaline increased it. In conclusion, the primary head-kidney cells of maraena whitefish reflect the immunological and neuroendocrinological diversity of the entire organ. This in vitro system allowed thus identifying the correlative changes between the activities of hormones and immune factors in salmonid fish in order to contribute to a better understanding of the regulation circuit between stress and immune defence.

Four ‘protein inhibitors of activated STAT’ (PIAS) control STAT-dependent and NF-κB-dependent immune signalling in humans. The genome of rainbow trout (Oncorhynchus mykiss) contains eight pias genes, which encode at least 14 different pias transcripts that are differentially expressed in a tissue- and cell-specific manner. Pias1a2 was the most strongly expressed variant among the analysed pias genes in most tissues, while pias4a2 was commonly low or absent. Since the knock-out of Pias factors in salmonid CHSE cells using CRISPR/Cas9 technology failed, three structurally different Pias protein variants were selected for overexpression studies in CHSE-214 cells. All three factors quenched the basal activity of an NF-κB promoter in a dose-dependent fashion, while the activity of an Mx promoter remained unaffected. Nevertheless, all three overexpressed Pias variants from trout strongly reduced the transcript level of the antiviral Stat-dependent mx gene in ifnγ-expressing CHSE-214 cells. Unlike mx, the overexpressed Pias factors modulated the transcript levels of NF-κB-dependent immune genes (mainly il6, il10, ifna3, and stat4) in ifnγ-expressing CHSE-214 cells in different ways. This dissimilar modulation of expression may result from the physical cooperation of the Pias proteins from trout with differential sets of interacting factors bound to distinct nuclear structures, as reflected by the differential nuclear localisation of trout Pias factors. In conclusion, this study provides evidence for the multiplication of pias genes and their sub-functionalisation during salmonid evolution.

Water temperature is doubtlessly one of the most important factors affecting the well-being of fish. For pikeperch Sander lucioperca an ambient water temperature of >20 °C has been suggested to achieve optimal growth in aquaculture farms, although natural habitats of pikeperch seldom reach such temperatures. We used multiplex qPCR arrays to assess the impact of temperature between 15 °C and 25 °C based on the expression of 38 potential biomarker genes characterising heat shock and hypoxia response, immune activation, energy supply and development in the liver and gills of pikeperch. In addition, liquid chromatography/mass spectrometry was used to assess the individual levels of glucocorticoids. Our data revealed that pikeperch do not develop severe responses to ambient temperatures between 15 °C and 25 °C, although specific parameters indicated a phased restriction of optimal conditions, above all the heat-shock genes SERPINH1, HSP90AA1 and HSPA8b; the acute-phase genes HP, C3 and SAA; and the transcription factor-encoding gene HMX1. In contrast to these in vivo observations, a temperature rise up to 37 °C led to an acute response in WF2-Sander cells characterised by the pronounced upregulation of HSP90AA1, HSP1A1, SERPINH1, HSPB1, HSPA8b and HSPB8. This report paves the ground for future assessment studies on the imbalanced homeostasis of pikeperch by means of molecular tools.

Single-cell RNA-sequencing (scRNA-seq) provides high-resolution insights into complex tissues. Cardiac tissue, however, poses a major challenge due to the delicate isolation process and the large size of mature cardiomyocytes. Regardless of the experimental technique, captured cells are often impaired and some capture sites may contain multiple or no cells at all. All this refers to “low quality” potentially leading to data misinterpretation. Common standard quality control parameters involve the number of detected genes, transcripts per cell, and the fraction of transcripts from mitochondrial genes. While cutoffs for transcripts and genes per cell are usually user-defined for each experiment or individually calculated, a fixed threshold of 5% mitochondrial transcripts is standard and often set as default in scRNA-seq software. However, this parameter is highly dependent on the tissue type. In the heart, mitochondrial transcripts comprise almost 30% of total mRNA due to high energy demands. Here, we demonstrate that a 5%-threshold not only causes an unacceptable exclusion of cardiomyocytes but also introduces a bias that particularly discriminates pacemaker cells. This effect is apparent for our in vitro generated induced-sinoatrial-bodies (iSABs; highly enriched physiologically functional pacemaker cells), and also evident in a public data set of cells isolated from embryonal murine sinoatrial node tissue (Goodyer William et al. in Circ Res 125:379–397, 2019). Taken together, we recommend omitting this filtering parameter for scRNA-seq in cardiovascular applications whenever possible. Supplementary Information The online version contains supplementary material available at 10.1007/s00018-021-03916-5.

Pikeperch (Sander lucioperca) is a fish species with growing economic significance in the aquaculture industry. However, successful positioning of pikeperch in large-scale aquaculture requires advances in our understanding of its genome organization. In this study, an ultra-high density linkage map for pikeperch comprising 24 linkage groups and 1,023,625 single nucleotide polymorphisms markers was constructed after genotyping whole-genome sequencing data from 11 broodstock and 363 progeny, belonging to 6 full-sib families. The sex-specific linkage maps spanned a total of 2985.16 cM in females and 2540.47 cM in males with an average inter-marker distance of 0.0030 and 0.0026 cM, respectively. The sex-averaged map spanned a total of 2725.53 cM with an average inter-marker distance of 0.0028 cM. Furthermore, the sex-averaged map was used for improving the contiguity and accuracy of the current pikeperch genome assembly. Based on 723,360 markers, 706 contigs were anchored and oriented into 24 pseudomolecules, covering a total of 896.48 Mb and accounting for 99.47% of the assembled genome size. The overall contiguity of the assembly improved with a scaffold N50 length of 41.06 Mb. Finally, an updated annotation of protein-coding genes and repetitive elements of the enhanced genome assembly is provided at NCBI.

Cardiovascular diseases represent a major health concern worldwide with few therapy options for ischemic injuries due to the limited regeneration potential of affected cardiomyocytes. Innovative cell replacement approaches could facilitate efficient regenerative therapy. However, despite extensive attempts to expand primary human cells in vitro, present technological limitations and the lack of human donors have so far prevented their broad clinical use. Cell xenotransplantation might provide an ethically acceptable unlimited source for cell replacement therapies and bridge the gap between waiting recipients and available donors. Pigs are considered the most suitable candidates as a source for xenogeneic cells and tissues due to their anatomical and physiological similarities with humans. The potential of porcine cells in the field of stem cell-based therapy and regenerative medicine is under intensive investigation. This review outlines the current progress and highlights the most promising approaches in xenogeneic cell therapy with a focus on the cardiovascular system.

Simple Summary Inadequate oxygen saturation, or hypoxia, belongs to one of the critical stress factors in intensive aquaculture. Exposure of fish to low oxygen levels over prolonged periods substantially affects their well-being and immune competence, resulting in increased disease susceptibility and consequent economic losses. In this interdisciplinary research, we aimed to provide a deeper understanding of the effect of chronic low oxygen saturation on pikeperch farmed in recirculating aquaculture systems. The obtained data offer unprecedented insights into the changes in the immunocompetence of studied fish and suggest high robustness of this new aquaculture species to the stress factors of intensive aquaculture. Abstract Inadequate oxygen saturation can induce stress responses in fish and further affect their immunity. Pikeperch, recently introduced in intensive aquaculture, is suggested to be reared at nearly 100% DO (dissolved oxygen), yet this recommendation can be compromised by several factors including the water temperature, stocking densities or low circulation. Herein, we aimed to investigate the effect of low oxygen saturation of 40% DO (±3.2 mg/L) over 28 days on pikeperch farmed in recirculating aquaculture systems. The obtained data suggest that—although the standard blood and health parameters did not reveal any significant differences at any timepoint—the flow cytometric analysis identified a slightly decreased proportion of lymphocytes in the HK (head kidney) of fish exposed to hypoxia. This has been complemented by marginally downregulated expression of investigated immune and stress genes in HK and liver (including FTH1, HIF1A and NR3C1). Additionally, in the model of acute peritoneal inflammation induced with inactivated Aeromonas hydrophila, we observed a striking dichotomy in the sensitivity to the low DO between innate and adaptive immunity. Thus, while the mobilization of myeloid cells from HK to blood, spleen and peritoneal cavity, underlined by changes in the expression of key proinflammatory cytokines (including MPO, IL1B and TNF) was not influenced by the low DO, hypoxia impaired the influx of lymphocytes to the peritoneal niche in the later phases of the immune reaction. Taken together, our data suggest high robustness of pikeperch towards the low oxygen saturation and further encourage its introduction to the intensive aquaculture systems.

A rapid decline in temperature poses a major challenge for poikilothermic fish, as their entire metabolism depends on ambient temperature. The gene expression of rainbow trout Oncorhynchus mykiss having undergone such a cold shock (0◦C) was compared to a control (5◦C) in a microarray and quantitative real-time PCR based study. The tissues of gill, kidney and liver were examined. The most differently expressed genes were found in liver, many of them contributing to the network 'cellular compromise, cellular growth and proliferation'.However, the number of genes found to be regulated at 0◦Cwas surprisingly low. Instead of classical genes involved in temperature shock, the three genes encoding fibroblast growth factor 1 (fgf1), growth arrest and DNA-damageinducible, alpha (gadd45a) and sclerostin domain-containing protein 1 (sostdc1) were upregulated in the liver upon cold shock in two different rainbow trout strains, suggesting that these genes may be considered as general biomarkers for cold shock in rainbow trout.

Stocking density is considered as a key factor determining the productivity of fish aquaculture systems. The transcriptomic response to crowding stress is, however, still poorly investigated. We aimed to identify potential biomarker genes via microarray analyses to gain insight into molecular pathways modulated by density-induced stress in farmed rainbow trout Oncorhynchus mykiss. Transcriptome profiling in the liver, kidney, and gills was complemented by analysis of classical plasma parameters and behavioural observations. Individuals of two trout strains were exposed for nine days to definite stocking densities, ~ 10 kg/m3 (low density); ~ 30 kg/m3 (moderate); ~ 60 kg/m3 (elevated); and ~ 120 kg/m3 (high). Only the PYGMA gene, encoding the key enzyme of glycogenolysis, was upregulated in gills and kidney, while about 90 genes were differentially regulated in the liver of trout exposed to high density compared with the moderate-density group. Pathway enrichment analyses predicted the upregulation of HIF1A signalling in liver, contributing to glucose homeostasis during stress conditions. This is consistent with the increased plasma glucose levels detected in crowding-stressed trout. The mTOR and PI3K/AKT signalling pathways were impaired in liver at higher stocking densities indicating stress at the cellular level. Contrary to our expectations, protein ubiquitination and cholesterol biosynthesis pathways were found to be underexpressed. Cholesterol is a precursor of the stress hormone cortisol, whose plasma concentrations were indeed within the normal physiological range in all treatment groups. In addition, we did not observe any aggressive interactions, which are known to be triggered through cortisol release. Taken together, we identified a relatively modest number of rainbow trout-specific crowding markers suggesting that farmed rainbow trout have gained enhanced stress-tolerant phenotypes through decades of selective breeding. Statement of relevance The comprehensive recording of animal-based biomarkers is becoming increasingly important for the further development of the aquaculture industry with importance for fish well-being and consumer opinion. Previous studies on crowding stress in rainbow trout have used single or multi-gene approaches to assess the relevance of putative biomarkers, but transcriptomic approaches to define crowding stress-induced pathways in salmonids are overdue. The present study compares the transcriptomic profiles of rainbow trout exposed to different stocking densities and correlates these data sets with classical plasma parameters and behavioural observations.

In stock enhancement and sea-ranching procedures, the adipose fin of hundreds of millions of salmonids is removed for marking purposes annually. However, recent studies proved the significance of the adipose fin as a flow sensor and attraction feature. In the present study, we profiled the specific expression of 20 neuron- and glial cell-marker genes in the adipose fin and seven other tissues (including dorsal and pectoral fin, brain, skin, muscle, head kidney, and liver) of the salmonid species rainbow trout Oncorhynchus mykiss and maraena whitefish Coregonus maraena. Moreover, we measured the transcript abundance of genes coding for 15 mechanoreceptive channel proteins from a variety of mechanoreceptors known in vertebrates. The overall expression patterns indicate the presence of the entire repertoire of neurons, glial cells and receptor proteins on the RNA level. This quantification suggests that the adipose fin contains considerable amounts of small nerve fibers with unmyelinated or slightly myelinated axons and most likely mechanoreceptive potential. The findings are consistent for both rainbow trout and maraena whitefish and support a previous hypothesis about the innervation and potential flow sensory function of the adipose fin. Moreover, our data suggest that the resection of the adipose fin has a stronger impact on the welfare of salmonid fish than previously assumed.

Analyses on the cellular level are indispensable to expand our understanding of complex tissues like the mammalian heart. Single-nucleus sequencing (snRNA-seq) allows for the exploration of cellular composition and cell features without major hurdles of single-cell sequencing. We used snRNA-seq to investigate for the first time an entire adult mammalian heart. Single-nucleus quantification and clustering led to an accurate representation of cell types, revealing 24 distinct clusters with endothelial cells (28.8%), fibroblasts (25.3%), and cardiomyocytes (22.8%) constituting the major cell populations. An additional RNA velocity analysis allowed us to study transcription kinetics and was utilized to visualize the transitions between mature and nascent cellular states of the cell types. We identified subgroups of cardiomyocytes with distinct marker profiles. For example, the expression of Hand2os1 distinguished immature cardiomyocytes from differentiated cardiomyocyte populations. Moreover, we found a cell population that comprises endothelial markers as well as markers clearly related to cardiomyocyte function. Our velocity data support the idea that this population is in a trans-differentiation process from an endothelial cell-like phenotype towards a cardiomyocyte-like phenotype. In summary, we present the first report of sequencing an entire adult mammalian heart, providing realistic cell-type distributions combined with RNA velocity kinetics hinting at interrelations.

The interleukin-1-receptor-associated kinase 3 (IRAK3) is known in mammals as a negative feedback regulator of NF-κB-mediated innate-immune mechanisms. Our RNA-seq experiments revealed a prototypic 1920-nt sequence encoding irak3 from rainbow trout (Oncorhynchus mykiss), as well as 20 variants that vary in length and nucleotide composition. Based on the DNA-sequence information from two closely related irak3 genes from rainbow trout and an irak3-sequence fragment from Atlantic salmon retrieved from public databases, we elucidated the underlying genetic causes for this striking irak3 diversity. Infecting rainbow trout with a lethal dose of Aeromonas salmonicida enhanced the expression of all variants in the liver, head kidney, and peripheral blood leucocytes. We analyzed the functional impact of the full-length factor and selected structural variants by overexpressing them in mammalian HEK-293 cells. The full-length factor enhanced the basal activity of NF-κB, but did not dampen the TLR2-signaling-induced levels of NF-κB activation. Increasing the basal NF-κB-activity through Irak3 apparently does not involve its C-terminal domain. However, more severely truncated factors had only a minor impact on the activity of NF-κB. The TLR2-mediated stimulation did not alter the spatial distribution of Irak3 inside the cells. In salmonid CHSE-214 cells, we observed that the Irak3-splice variant that prominently expresses the C-terminal domain significantly quenched the stimulation-dependent production of interleukin-1β and interleukin-8, but not the production of other immune regulators. We conclude that the different gene and splice variants of Irak3 from trout play distinct roles in the activation of immune-regulatory mechanisms.

The pikeperch (Sander lucioperca) is a fresh and brackish water Percid fish natively inhabiting the northern hemisphere. This species is emerging as a promising candidate for intensive aquaculture production in Europe. Specific traits like cannibalism, growth rate and meat quality require genomics based understanding, for an optimal husbandry and domestication process. Still, the aquaculture community is lacking an annotated genome sequence to facilitate genome-wide studies on pikeperch. Here, we report the first highly contiguous draft genome assembly of Sander lucioperca. In total, 413 and 66 giga base pairs of DNA sequencing raw data were generated with the Illumina platform and PacBio Sequel System, respectively. The PacBio data were assembled into a final assembly size of ~900 Mb covering 89% of the 1,014 Mb estimated genome size. The draft genome consisted of 1966 contigs ordered into 1,313 scaffolds. The contig and scaffold N50 lengths are 3.0 Mb and 4.9 Mb, respectively. The identified repetitive structures accounted for 39% of the genome. We utilized homologies to other ray-finned fishes, and ab initio gene prediction methods to predict 21,249 protein-coding genes in the Sander lucioperca genome, of which 88% were functionally annotated by either sequence homology or protein domains and signatures search. The assembled genome spans 97.6% and 96.3% of Vertebrate and Actinopterygii single-copy orthologs, respectively. The outstanding mapping rate (99.9%) of genomic PE-reads on the assembly suggests an accurate and nearly complete genome reconstruction. This draft genome sequence is the first genomic resource for this promising aquaculture species. It will provide an impetus for genomic-based breeding studies targeting phenotypic and performance traits of captive pikeperch.

The creatine/phosphocreatine system is the principal energy buffer in mammals, but is scarcely documented in fish. We measured the gene expression of major enzymes of this system, glycine amidinotransferase (GATM), guanidinoacetate N-methyltransferase (GAMT) and muscle-type creatine kinase (CKM) in kidney, liver, and muscle tissues of fish and mammals. CKM was expressed strongly in the muscles of all examined species. In contrast, GATM and GAMT were strongly expressed in the muscle tissue of fish, but not of mammals. This indicates that creatine synthesis and usage are spatially separated in mammals, but not in fish, which is supported by RNA-Seq data of 25 species. Differences in amino acid metabolism along with methionine adenosyltransferase gene expression in muscle from fishes but not mammals further support a central metabolic role of muscle in fish, and hence different organization of the creatine/phosphocreatine biosynthesis system in higher and lower vertebrates.

The mammalian toll-like receptor 2 (TLR2) is a dominant receptor for the recognition of Gram-positive bacteria. Its structure and functional properties were unknown in salmonid fish. In RT-PCR and RACE experiments, we obtained the full-length cDNA sequence encoding Tlr2 from rainbow trout (Oncorhynchus mykiss) as well as a copy of an unspliced nonsense message from a highly segmented gene. The primary structure of the encoded receptor complies with the domain structure and ligand-binding sites known from mammals and other fish species and sorts well into the evolutionary tree of teleostean Tlr2s. We retrieved a gene version encoding the receptor on a single exon (tlr2a) and also a partial sequence of a second gene variant being segmented into multiple exons (tlr2b). Surprisingly, the abundances of both transcript variants accounted only for ∼10% of all Tlr2-encoding transcripts in various tissues and cell types of healthy fish. This suggests the expression of several distinct tlr2 gene variants in rainbow trout. We expressed tlr2a in HEK-293 cells, but were unable to demonstrate its functionality through NF-κB activation. Neither synthetic lipopeptides known to stimulate mammalian TLR2 nor different bacterial challenges induced OmTLR2-mediated NF-κB activation, not in HEK-293 or in salmon CHSE-214 cells. Positive demonstration of TLR2-MYD88 interaction excluded that its functional impairment caused the failure of NF-κB activation. We discuss impaired heterodimerization with a necessary Tlr partner as one from among several alternatives to explain the dysfunction of Tlr2a in the interspecies reconstitution system of TLR signaling.

Selective breeding can significantly improve the establishment of sustainable and profitable aquaculture fish farming. For rainbow trout (Oncorhynchus mykiss), one of the main aquaculture coldwater species in Europe, a variety of selected hatchery strains are commercially available. In this study, we investigated the genetic variation between the local Born strain, selected for survival, and the commercially available Silver Steelhead strain, selected for growth. We sequenced the transcriptome of six tissues (gills, head kidney, heart, liver, spleen, and white muscle) from eight healthy individuals per strain, using RNA-seq technology to identify strain-specific gene-expression patterns and single nucleotide polymorphisms (SNPs). In total, 1760 annotated genes were differentially expressed across all tissues. Pathway analysis assigned them to different gene networks. We also identified a set of SNPs, which are heterozygous for one of the two breeding strains: 1229 of which represent polymorphisms over all tissues and individuals. Our data indicate a strong genetic differentiation between Born and Silver Steelhead trout, despite the relatively short time of evolutionary separation of the two breeding strains. The results most likely reflect their specifically adapted genotypes and might contribute to the understanding of differences regarding their robustness toward high stress and pathogenic challenge described in former studies.

Apoptosis is an integral part of homeostasis and supports multiple physiological processes such as development and immune defense, thereby directly targeting damaged or unwanted cells without affecting neighbor cells. In the present study, we characterized the apoptotic key factors caspase-3, −7, and − 8 as well as regulator protein TPT1 (translationally-controlled tumor protein 1) from rainbow trout (Oncorhynchus mykiss). We identified multiple single-nucleotide changes in their coding sequences and showed that the CASP3 gene is present in at least three variants. Caspase genes were clustered to their orthologs in bony fish and human by using evolutionary analysis. Expression profiling in seven tissues of unchallenged adult fish revealed predominant transcript levels in the head kidney (CASP3, 7, and 8) or brain (TPT1). Further, we analyzed the expression of a more comprehensive panel of 16 trout genes encoding pro- and anti-apoptotic factors and associated proteins during development and upon stress exposure (in vitro temperature and staurosporine treatment). Previously published transcriptome data suggested that the induction of apoptotic processes is mirrored on the transcript level, but this could not be confirmed by the present gene-profiling study. Yet on the protein level, treatment of trout cell line RT-gill-W1 with 1 μM staurosporine for up to 120 min led to a significant increase of CASP3/7 activity. Moreover, a meta-analysis on published data showed that stress-related expression could only be detected sporadically for apoptotic key factors. In conclusion, there seems to be no reliable pattern or marker representing the stress-related induction of apoptosis in salmonids.

Siglecs (sialic acid-binding immunoglobulin-type lectins) are a family of immune regulatory receptors predominantly found on the cells of the hematopoietic system. A V-set Ig-like domain mediates the recognition of different sialylated glycoconjugates, which can lead to the activation or inhibition of the immune response, depending on the involved Siglecs. Siglecs are categorized into two subgroups: one including all CD33-related Siglecs and the other consisting of Siglec-1 (Sialoadhesin), Siglec-2 (CD22), Siglec-4 (myelin-associated glycoprotein, MAG) and Siglec-15. In contrast to the members of the CD33-related Siglecs, which share ∼50–99% sequence identity, Siglecs of the other subgroup show quite low homology (approximately 25–30% sequence identity). Based on the published sequences and functions of Siglecs, we performed phylogenetic analyses and sequence alignments to reveal the conservation of Siglecs throughout evolution. Therefore, we focused on the presence of Siglecs in different classes of vertebrates (fishes, amphibians, birds, reptiles and mammals), offering a bridge between the presence of different Siglecs and the biological situations of the selected animals.

The innate immune response involves a concerted network of induced gene products, preformed immune effectors, biochemical signalling cascades and specialised cells. However, the multifaceted activation of these defensive measures can derail or overshoot and, if left unchecked, overwhelm the host. A plenty of regulatory devices therefore mediate the fragile equilibrium between pathogen defence and pathophysiological manifestations. Over the past decade in particular, an almost complete set of teleostean sequences orthologous to mammalian immunoregulatory factors has been identified in various fish species, which prove the remarkable conservation of innate immune-control concepts among vertebrates. This review will present the current knowledge on more than 50 teleostean regulatory factors (plus additional fish-specific paralogs) that are of paramount importance for controlling the clotting cascade, the complement system, pattern-recognition pathways and cytokine-signalling networks. A special focus lies on those immunoregulatory features that have emerged as potential biomarker genes in transcriptome-wide research studies. Moreover, we report on the latest progress in elucidating control elements that act directly with immune-gene-encoding nucleic acids, such as transcription factors, hormone receptors and micro- and long noncoding RNAs. Investigations into the function of teleostean inhibitory factors are still mainly based on gene-expression profiling or overexpression studies. However, in support of structural and in-vitro analyses, evidence from in-vivo trials is also available and revealed many biochemical details on piscine immune regulation. The presence of multiple gene copies in fish adds a degree of complexity, as it is so far hardly understood if they might play distinct roles during inflammation. The present review addresses this and other open questions that should be tackled by fish immunologists in future.

A key common feature of all but three known mammalian genera is the strict seven cervical vertebrae blueprint, suggesting the involvement of strong conserving selection forces during mammalian radiation. This is further supported by reports indicating that children with cervical ribs die before they reach reproductive age. Hypotheses were put up, associating cervical ribs (homeotic transformations) to embryonal cancer (e.g., neuroblastoma) or ascribing the constraint in cervical vertebral count to the development of the mammalian diaphragm. Here, we describe a spontaneous mutation c.196A > G in the Bos taurus T gene (also known as brachyury) associated with a cervical vertebral homeotic transformation that violates the fundamental mammalian cervical blueprint, but does not preclude reproduction of the affected individual. Genome-wide mapping, haplotype tracking within a large pedigree, resequencing of target genome regions, and bioinformatic analyses unambiguously confirmed the mutant c.196G allele as causal for this previously unknown defect termed vertebral and spinal dysplasia (VSD) by providing evidence for the mutation event. The nonsynonymous VSD mutation is located within the highly conserved T box of the T gene, which plays a fundamental role in eumetazoan body organization and vertebral development. To our knowledge, VSD is the first unequivocally approved spontaneous mutation decreasing cervical vertebrae number in a large mammal. The spontaneous VSD mutation in the bovine T gene is the first in vivo evidence for the hypothesis that the T protein is directly involved in the maintenance of the mammalian seven-cervical vertebra blueprint. It therefore furthers our knowledge of the T-protein function and early mammalian notochord development.

Toll-like receptors (TLRs) are known to detect a defined spectrum of microbial structures. However, the knowledge about the specificity of teleost Tlr factors for distinct pathogens is limited so far. We measured baseline expression profiles of 18tlrgenes and associated signaling factors in four immune-relevant tissues of rainbow troutOncorhynchus mykiss. Intraperitoneal injection of a lethal dose ofAeromonas salmonicidasubsp.salmonicidainduced highly increased levels of cytokine mRNAs during a 72-hour postinfection (hpi) period. In contrast, only the fish-specifictlr22a2and the downstream factorirak1featured clearly increased transcript levels, while the mRNA concentrations of many othertlrgenes decreased. Flow cytometry quantified cell trafficking after infection indicating a dramatic influx of myeloid cells into the peritoneum and a belated low level immigration of lymphoid cells. T and B lymphocytes were differentiated with RT-qPCR revealing that B lymphocytes emigrated from and T lymphocytes immigrated into head kidney. In conclusion, no specific TLR can be singled out as a dominant receptor forA. salmonicida. The recruitment of cellular factors of innate immunity rather than induced expression of pathogen receptors is hence of key importance for mounting a first immune defense against invadingA. salmonicida.

Salmon lice (Lepeophtheirus salmonis) are parasitic copepods, living mainly on Atlantic salmon and leading to large economical losses in aquaculture every year. Due to the emergence of resistances to several drugs, alternative treatments are developed, including treatment with hydrogen peroxide, freshwater or thermal treatment. The present study gives a first overview of the thermotolerance and stress response of salmon lice. Sea lice nauplii acclimated to 10 °C can survive heat shocks up to 30 °C and are capable of hardening by a sublethal heat shock. We searched in the genome for heat shock protein (HSP) encoding genes and tested their inducibility after heat shock, changes in salinity and treatment with hydrogen peroxide, employing microfluidic qPCRs. We assessed 38 candidate genes, belonging to the small HSP, HSP40, HSP70 and HSP90 families. Nine of these genes showed strong induction after a non-lethal heat shock. In contrast, only three and two of these genes were induced after changes in salinity and incubation in hydrogen peroxide, respectively. This work provides the basis for further work on the stress response on the economically important parasite L. salmonis.

We describe an emerging initiative - the ‘Functional Annotation of All Salmonid Genomes’ (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance. © The Author(s). 2017 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/)

Interleukin-8 (IL8) is an immediate-early chemokine that has been well characterized in several fish species. Ten IL8 gene variants have already been described in rainbow trout, but none of their promoters has structurally been defined or functionally characterized in teleost fish. To uncover key factors regulating IL8 expression, we intended to functionally characterize an IL8 promoter from rainbow trout. Incidentally, we isolated a novel IL8 gene variant (IL8-G). It is structurally highly similar to the other trout IL8 gene variants and its mRNA concentration increased significantly in secondary lymphoid tissues after infecting healthy fish with Aeromonas salmonicida. The proximal promoter sequence of the IL8-G-encoding gene features in close proximity two consensus elements for CEBP attachment. The proximal site overlaps with a NF-κB-binding site. Cotransfection of an IL8-G promoter-driven reporter gene together with vectors expressing various mammalian CEBP or NF-κB factors revealed in human HEK-293 cells that CEBPA and NF-κB p50, but not NF-κB p65 activate this promoter. The stimulatory effect of NF-κB p50 is likely conveyed by synergizing with CEBPA. Deletion or mutation of either the distal or the proximal CEBP-binding site, respectively, caused a significant decrease in IL8-G promoter activation. We confirmed the significance of the CEBPA factor for IL8-G expression by comparing the stimulatory capacity of the trout CEBPA and -B factors, thereby reducing the evolutionary distance in the inter-species expression assays. Similar promoter induction potential and intracellular localization of the mammalian and teleostean CEBPA and -B factors suggests their functional conservation throughout evolution.

Adverse life circumstances evoke a common “conserved transcriptional response to adversity” (CTRA) in mammalian leukocytes. To investigate whether this pattern is preserved in lower vertebrates, maraena whitefish (Coregonus maraena) were exposed for nine days to different stocking densities: ~ 10 kg/m³ (low density), ~ 33 kg/m³ (moderate), ~ 60 kg/m³ (elevated), and ~ 100 kg/m³ (high). Transcriptome profiling in the liver and kidney of individuals from each group suggested that crowding conditions activate stress-related signaling and effector pathways. Remarkably, about one quarter of the genes differentially expressed under crowding conditions were involved in the activation of immune pathways such as acute-phase response and interleukin/TNF signaling attended by the simultaneous reduction of antiviral potency. Network analysis confirmed the complex interdigitation of immune- and stress-relevant pathways with interleukin-1 playing a central role. Antibody-based techniques revealed remarkable changes in the blood composition of whitefish and demonstrated the correlation between increasing stocking densities and elevated number of myeloid cells together with the increased phagocytic activity of peripheral blood leukocytes. In line with current studies in mammals, we conclude that crowding stress triggers in whitefish hallmarks of a CTRA, indicating that the stress-induced molecular mechanisms regulating the immune responses not only are conserved within mammals but were established earlier in evolution.

We describe an emerging initiative - the ‘Functional Analysis of All Salmonid Genomes’ (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, through to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

The mammalian interleukin 1 receptor-like 1 receptor (IL1RL1), commonly known as ST2, is thought to downregulate TLR signalling by sequestering the signalling adapter MYD88 (myeloid differentiation primary response protein 88). ST2 sequences are known in several fish species, but none of them have functionally been examined. We characterised ST2 from rainbow trout (Oncorhynchus mykiss) and the structure of its encoding gene. The primary sequence of ST2 is only weakly conserved from fish to human. However, the amino acid sequences forming the interfaces for ST2 and MYD88 interaction are well conserved throughout evolution. High similarity of the gene segmentation unambiguously proves the common ancestry of fish and mammalian ST2. Trout ST2 and trout MYD88 genes were constitutively expressed in embryonic, larval and adult trout. In vivo infection with Aeromonas salmonicida did not modulate the mRNA levels of both factors. Overexpressing trout ST2 in the mammalian HEK-293 reconstitution system of TLR2 signalling quenched the Escherichia coli–induced activation of NF-κB and SAA promoters in a dose-dependent fashion. The expression of GFP-tagged trout ST2 in human HEK-293 or trout CHSE-214 cells surprisingly revealed that (i) ST2 localised abundantly at the nuclear membrane rather than at the cell membrane and (ii) the coexpression of both ST2 and MYD88 allowed the translocation of trout MYD88 from cytoplasm to nucleus, as assessed using confocal microscopy and Western blotting. Hence, we validated that trout ST2 is a dampener of TLR signalling and interacts with MYD88. The spatial distribution of these factors raises questions about how this repressive mechanism functions.