Immunologic Profiling of the Atlantic Salmon Gill by Single Nuclei Transcriptomics

Hypersalinity tolerance of mummichogs (Fundulus heteroclitus): A branchial transcriptomic analysis

Along the east coast of North America, mummichogs (Fundulus heteroclitus) are subjected to a broad range of salinities in their nearshore habitats. However, there is a paucity of information regarding the molecular and cellular processes that mummichogs (and other highly osmotolerant fishes) engage to survive environmental salinities greater than seawater (SW). To reveal branchial processes underlying their extraordinarily broad salinity tolerance, we performed an RNA-Seq analysis to identify differentially expressed genes (DEGs) in mummichogs residing in 3, 35, and 105 ppt conditions. We identified a series of DEGs previously associated with both freshwater (FW)- and SW-type ionocytes; however, the heightened expression of anoctamin 1a, a Ca2+-activated Cl- channel, in 35 and 105 ppt indicates that an undescribed Cl--secretion pathway may operate within the SW-type ionocytes of mummichogs. Concerning FW-adaptive branchial processes, we identified claudin 5a as a gene whose product may limit the diffusive loss of ions between cellular tight junctions. Further, in response to hypersaline conditions, we identified DEGs linked with myo-inositol synthesis and kinase signaling. This study provides new molecular targets for future physiological investigations that promise to reveal the mechanistic bases for how mummichogs and other euryhaline species tolerate hypersaline conditions.

Cell proliferation and regeneration in the gill

Seminal studies from the early 20th century defined the structural changes associated with development and regeneration of the gills in goldfish at the gross morphological and cellular levels using standard techniques of light and electron microscopy. More recently, investigations using cell lineage tracing, molecular biology, immunohistochemistry and single-cell RNA-sequencing have pushed the field forward and have begun to reveal the cellular and molecular processes that orchestrate cell proliferation and regeneration in the gills. The gill is a multifunctional organ that mediates an array of important physiological functions, including respiration, ion regulation and excretion of waste products. It is comprised of unique cell types, such as pavement cells, ionocytes, chemoreceptors and undifferentiated stem or progenitor cells that regulate growth and replenish cell populations. The gills develop from the embryonic endoderm and are rich in cell types derived from the neural crest. The gills have the capacity to remodel themselves in response to environmental change, such as in the case of ionocytes, chemoreceptors and the interlamellar cell mass, and can completely regenerate gill filaments and lamellae. Both processes of remodeling and regeneration invariably involve cell proliferation. Although gill regeneration has been reported in only a limited number of fish species, the process appears to have many similarities to regeneration of other organs in fish and amphibians. The present article reviews the studies that have described gill development and growth, and that demonstrate a suite of genes, transcription factors and other proteins involved in cell proliferation and regeneration in the gills.

Tanycytes from a bird's eye view: gene expression profiling of the tanycytic region under different seasonal states in the Svalbard ptarmigan

In mammals and birds, tanycytes are known to regulate thyroid hormone conversion, and this process is central to the control of seasonal reproduction. In mammals, this cell type is also implicated in retinoic acid signalling, neurogenesis, and nutritional gatekeeping, all of which have been linked to hypothalamic regulation of energy metabolism. Less is known about these potential wider roles of tanycytes in birds. To address this gap, we combined LASER capture microdissection and transcriptomics to profile the tanycytic region in male Svalbard ptarmigan, a High Arctic species with photoperiod-dependent seasonal rhythms in reproductive activation and body mass. Short photoperiod (SP) adapted birds were transferred to constant light (LL) to trigger breeding and body mass loss. After five months under LL, the development of photorefractoriness led to spontaneous re-emergence of the winter phenotype, marked by the termination of breeding and gain in body mass. The transfer from SP to LL initiated gene expression changes in both thyroid hormone and retinoic acid pathways, as described in seasonal mammals. Furthermore, transcriptomic signatures of cell differentiation and migration were observed. Comparison to data from Siberian hamsters demonstrated that a photoperiod-dependent re-organisation of the hypothalamic tanycytic region is likely a conserved feature. Conversely, the spontaneous development of photorefractoriness showed a surprisingly small number of genes that reverted in expression level, despite reversal of the reproductive and metabolic phenotype. Our data suggest general conservation of tanycyte biology between photoperiodic birds and mammals and raise questions about the mechanistic origins of the photorefractory state.

Unveiling osmoregulation and immunological adaptations in Eleutheronema tetradactylum gills through high-throughput single-cell transcriptome sequencing

The fourfinger threadfin fish (Eleutheronema tetradactylum) is an economically significant species renowned for its ability to adapt to varying salinity environments, with gills serving as their primary organs for osmoregulation and immune defense. Previous studies focused on tissue and morphological levels, whereas ignored the cellular heterogeneity and the crucial gene information related to core cell subsets within E. tetradactylum gills. In this study, we utilized high-throughput single-cell RNA sequencing (scRNA-seq) to analyze the gills of E. tetradactylum, characterizing 16 distinct cell types and identifying unique gene markers and enriched functions associated within each cell type. Additionally, we subdivided ionocyte cells into four distinct subpopulations for the first time in E. tetradactylum gills. By employing weighted gene co-expression network analysis (WGCNA), we further investigated the cellular heterogeneity and specific response mechanisms to salinity fluctuant. Our findings revealed the intricate osmoregulation and immune functions of gill cells, highlighting their crucial roles in maintaining homeostasis and adapting to fluctuating salinity levels. This comprehensive cell-type atlas provides valuable insights into the species adaptive strategies, contributing to the conservation and management of this commercially significant fish as well as other euryhaline species.

Impact of freshwater rearing history on Atlantic salmon gill response to viral stimulation post seawater transfer

Land-based recirculating aquaculture systems (RAS) have risen in prevalence in recent years for Atlantic salmon production, enabling intensive production which allows increased growth and environmental control, but also having the potential for reducing water use and eutrophication. The Atlantic salmon has an anadromous life history with juvenile stages in freshwater (FW) and on-growing in seawater (SW), enabled by a transformational process known as smoltification. The timing of smoltification and transfer of smolts from FW to SW is critical under commercial production with high mortalities during this period. The impact of FW rearing system on immune function following seawater transfer (SWT) is not well understood. In this study parr were raised in either RAS or a traditional open-LOCH system until smolting and then transferred to a common marine environment. Two-weeks post-SWT fish were immune stimulated with a viral mimic (poly I:C) for 24 h to assess the ability to mount an antiviral immune response, assessed by whole transcriptome analysis of gill tissue, an important immune organ in fish. We show that unstimulated smolts reared in the LOCH had higher immune gene expression than those reared in RAS as determined by functional analysis. However, following stimulation, smolts reared in the RAS mounted a greater magnitude of response with a suite of immune genes displaying higher fold induction of transcription compared to LOCH reared smolts. We suggest RAS smolts have a lower steady state immune-associated transcriptome likely due to an unvarying environment, in terms of environmental factors and lack of exposure to pathogens, which shows a compensatory mechanism following stimulation allowing immune 'catch-up' with those reared in the LOCH. Alternatively, the RAS fish are experiencing an excessive response to the immune stimulation.

Salinity change evokes stress and immune responses in Atlantic salmon with microalgae showing limited potential for dietary mitigation

Smoltification was found to impact both immune and stress responses of farmed Atlantic salmon (Salmo salar), but little is known about how salinity change affects salmon months after completed smoltification. Here, we examined (1) the effect of salinity change from brackish water to seawater on the stress and immune responses in Atlantic salmon and (2) evaluated if functional diets enriched with microalgae can mitigate stress- and immune-related changes. Groups of Atlantic salmon were fed for 8 weeks with different microalgae-enriched diets in brackish water and were then transferred into seawater. Samples of the head kidney, gill, liver and plasma were taken before seawater transfer (SWT), 20 h after SWT, and 2 weeks after SWT for gene-expression analysis, plasma biochemistry and protein quantification. The salmon showed full osmoregulatory ability upon transfer to seawater reflected by high nkaα1b levels in the gill and tight plasma ion regulation. In the gill, one-third of 44 investigated genes were reduced at either 20 h or 2 weeks in seawater, including genes involved in cytokine signaling (il1b) and antiviral defense (isg15, rsad2, ifit5). In contrast, an acute response after 20 h in SW was apparent in the head kidney reflected by increased plasma stress indicators and induced expression of genes involved in acute-phase response (drtp1), antimicrobial defense (camp) and stress response (hspa5). However, after 2 weeks in seawater, the expression of antiviral genes (isg15, rsad2, znfx1) was reduced in the head kidney. Few genes (camp, clra, c1ql2) in the gill were downregulated by a diet with 8% inclusion of Athrospira platensis. The results of the present study indicate that salinity change months after smoltification evokes molecular stress- and immune responses in Atlantic salmon. However, microalgae-enriched functional diets seem to have only limited potential to mitigate the related changes.

Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.

Physiological and immune profiling of tilapia Oreochromis niloticus gills by high-throughput single-cell transcriptome sequencing

The physiological and immune functions of fish gills are largely recognized, but their following functional heterogeneity at the single cell scale has been rarely reported. Here, we performed single cell RNA sequencing (scRNA-seq) on the gills of tilapia fish Oreochromis niloticus. We identified a total of 12 cell populations and analyzed their functional heterogeneity. To investigate the physiological function of O. niloticus gills, expression patterns of genes encoding ion transporters were selected from the identified H+-ATPase-rich cells (HR cells), Na+/K+-ATPase-rich cells (NaR cells), and pavement cells. Specific enrichment of ca4a, slc9a1a, and LOC100692482 in the HR cells of O. niloticus gills explained their functions in acid-base regulation. Genes encoding Ca2+ transporters, including atp2b1, LOC100696627, and LOC 100706765, were specifically expressed in the NaR cells. Pavement cells were presumably the main sites responsible for ammonia and urea transports in O. niloticus gills with specific enrichment of Rhbg and LOC100693008, respectively. The expression patterns of the four immune cell subtypes varied greatly, with B cells being enriched with the most immune-related GO terms. KEGG enrichment analysis showed that MAPK signaling pathway was the most enriched pathway among the four types of immune cells in O. niloticus gills. Our results are important in understanding the physiological and immune responses of fish gills at the cellular resolution.

A single-cell atlas of West African lungfish respiratory system reveals evolutionary adaptations to terrestrialization

The six species of lungfish possess both lungs and gills and are the closest extant relatives of tetrapods. Here, we report a single-cell transcriptome atlas of the West African lungfish (Protopterus annectens). This species manifests the most extreme form of terrestrialization, a life history strategy to survive dry periods that can last for years, characterized by dormancy and reversible adaptive changes of the gills and lungs. Our atlas highlights the cell type diversity of the West African lungfish, including gene expression consistent with phenotype changes of terrestrialization. Comparison with terrestrial tetrapods and ray-finned fishes reveals broad homology between the swim bladder and lung cell types as well as shared and idiosyncratic changes of the external gills of the West African lungfish and the internal gills of Atlantic salmon. The single-cell atlas presented here provides a valuable resource for further exploration of the respiratory system evolution in vertebrates and the diversity of lungfish terrestrialization.

Microbiomes in the context of developing sustainable intensified aquaculture

With an ever-growing human population, the need for sustainable production of nutritional food sources has never been greater. Aquaculture is a key industry engaged in active development to increase production in line with this need while remaining sustainable in terms of environmental impact and promoting good welfare and health in farmed species. Microbiomes fundamentally underpin animal health, being a key part of their digestive, metabolic and defense systems, in the latter case protecting against opportunistic pathogens in the environment. The potential to manipulate the microbiome to the advantage of enhancing health, welfare and production is an intriguing prospect that has gained considerable traction in recent years. In this review we first set out what is known about the role of the microbiome in aquaculture production systems across the phylogenetic spectrum of cultured animals, from invertebrates to finfish. With a view to reducing environmental footprint and tightening biological and physical control, investment in "closed" aquaculture systems is on the rise, but little is known about how the microbial systems of these closed systems affect the health of cultured organisms. Through comparisons of the microbiomes and their dynamics across phylogenetically distinct animals and different aquaculture systems, we focus on microbial communities in terms of their functionality in order to identify what features within these microbiomes need to be harnessed for optimizing healthy intensified production in support of a sustainable future for aquaculture.

Understanding host response to infectious salmon anaemia virus in an Atlantic salmon cell line using single-cell RNA sequencing

BACKGROUND

Infectious Salmon Anaemia Virus (ISAV) is an Orthomixovirus that represents a large problem for salmonid aquaculture worldwide. Current prevention and treatment methods are only partially effective. Genetic selection and genome engineering have the potential to develop ISAV resistant salmon stocks. Both strategies can benefit from an improved understanding of the genomic regulation of ISAV pathogenesis. Here, we used single-cell RNA sequencing of an Atlantic salmon cell line to provide the first high dimensional insight into the transcriptional landscape that underpins host-virus interaction during early ISAV infection.

RESULTS

Salmon head kidney (SHK-1) cells were single-cell RNA sequenced at 24, 48 and 96 h post-ISAV challenge. At 24 h post infection, cells showed expression signatures consistent with viral entry, with genes such as PI3K, FAK or JNK being upregulated relative to uninfected cells. At 48 and 96 h, infected cells showed a clear anti-viral response, characterised by the expression of IFNA2 or IRF2. Uninfected bystander cells at 48 and 96 h also showed clear transcriptional differences, potentially suggesting paracrine signalling from infected cells. These bystander cells expressed pathways such as mRNA sensing, RNA degradation, ubiquitination or proteasome; and up-regulation of mitochondrial ribosome genes also seemed to play a role in the host response to the infection. Correlation between viral and host genes revealed novel genes potentially key for this fish-virus interaction.

CONCLUSIONS

This study has increased our understanding of the cellular response of Atlantic salmon during ISAV infection and revealed host-virus interactions at the cellular level. Our results highlight various potential key genes in this host-virus interaction, which can be manipulated in future functional studies to increase the resistance of Atlantic salmon to ISAV.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Differential Expression of miRNAs and Their Predicted Target Genes Indicates That Gene Expression in Atlantic Salmon Gill Is Post-Transcriptionally Regulated by miRNAs in the Parr-Smolt Transformation and Adaptation to Sea Water

Smoltification (parr-smolt transformation) is a complex developmental process consisting of developmental changes that lead to remodeling of the Atlantic salmon gill. Here, the expression changes of miRNAs and mRNAs were studied by small-RNA sequencing and microarray analysis, respectively, to identify miRNAs and their predicted targets associated with smoltification and subsequent sea water adaptation (SWA). In total, 18 guide miRNAs were identified as differentially expressed (gDE miRNAs). Hierarchical clustering analysis of expression changes divided these into one cluster of 13 gDE miRNAs with decreasing expression during smoltification and SWA that included the miRNA-146, miRNA-30 and miRNA-7132 families. Another smaller cluster that showed increasing expression consisted of miR-101a-3p, miR-193b-5p, miR-499a-5p, miR-727a-3p and miR-8159-5p. The gDE miRNAs were predicted to target 747 of the genes (DE mRNAs), showing expression changes in the microarray analysis. The predicted targets included genes encoding NKA-subunits, aquaporin-subunits, cystic fibrosis transmembrane conductance regulator and the solute carrier family. Furthermore, the predicted target genes were enriched in biological processes associated with smoltification and SWA (e.g., immune system, reactive oxygen species, stress response and extracellular matrix organization). Collectively, the results indicate that remodeling of the gill involves the post-transcriptional regulation of gene expression by the characterized gDE miRNAs.

Immune responses to prebiotics in farmed salmonid fish: How transcriptomic approaches help interpret responses

Within aquaculture, prebiotics are composed of complex carbohydrate molecules that cannot be digested by the fish directly but are metabolised by the microbial communities within the host gut, with the desire that "healthy" bacterial species are promoted with subsequently improved performance of the fish, there are likely some direct responses of intestinal cells to these dietary components. The sources and processing of prebiotics, which fall under the overarching theme of "functional feeds" are highly varied between species and types of prebiotics administered. How these feeds exert their effect, and the host responses are hard to determine, but new technologies and the development of high-throughput technologies (omics) are enabling the mechanisms and methods of action to be further understood. The recent advances in the availability of 'omics' technologies with the transition from single gene assays to microarray and RNA-seq in fish health have enabled novel functional ingredients to be analysed. This review will focus on recent studies on targeted gene expression and 'omics' technologies to characterize immune responses. Comparisons between the immunomodulatory effect of different prebiotics have been made and specific examples of how transcriptomics techniques have been used to identify immune responses to prebiotics are given.

Single-cell transcriptomic analysis of neuroepithelial cells and other cell types of the gills of zebrafish (Danio rerio) exposed to hypoxia

The fish gill is a multifunctional organ involved in numerous physiological processes, such as gas exchange and sensing of hypoxia by respiratory chemoreceptors, called neuroepithelial cells (NECs). Many studies have focused on zebrafish (Danio rerio) to investigate the structure, function and development of the gills, yet the transcriptomic profile of most gill cells remains obscure. We present the results of a comprehensive transcriptomic analysis of the gills of zebrafish using single-cell RNA sequencing (scRNA‐seq). Gill cells from ETvmat2:EGFP zebrafish were individually labelled before scRNA‐seq library construction using 10× Genomics Chromium technology. 12,819 cells were sequenced with an average depth of over 27,000 reads per cell. We identified a median of 485 genes per cell and 16 cell clusters, including NECs, neurons, pavement cells, endothelial cells and mitochondrion-rich cells. The identity of NECs was confirmed by expression of slc18a2, encoding the vesicular monoamine transporter, Vmat2. Highly differentially-expressed genes in NECs included tph1a, encoding tryptophan hydroxylase, sv2 (synaptic vesicle protein), and proteins implicated in O₂ sensing (ndufa4l2a, cox8al and epas1a). In addition, NECs and neurons expressed genes encoding transmembrane receptors for serotonergic, cholinergic or dopaminergic neurotransmission. Differential expression analysis showed a clear shift in the transcriptome of NECs following 14 days of acclimation to hypoxia. NECs in the hypoxia group showed high expression of genes involved in cell cycle control and proliferation. The present article provides a complete cell atlas for the zebrafish gill and serves as a platform for future studies investigating the molecular biology and physiology of this organ.

Novel Insight Into Nutritional Regulation in Enhancement of Immune Status and Mediation of Inflammation Dynamics Integrated Study In Vivo and In Vitro of Teleost Grass Carp (Ctenopharyngodon idella): Administration of Threonine

This study aims to investigate the effects of threonine (Thr) on immunoregulation in vivo and in vitro of teleost grass carp (Ctenopharyngodon idella). Juveniles (9.53 ± 0.02 g) were reared for 8 weeks with respective Thr diet (3.99, 7.70, 10.72, 14.10, 17.96, and 21.66 g/kg) and then challenged with Aeromonas hydrophila for in vivo study. Macrophages isolated from head kidney were treated in vitro for 48 h with L-Thr (0, 0.5, 1.0, 2.0, 4.0, and 8.0 mM) after 6 h of lipopolysaccharide induction. The results showed that, compared with Thr deficiency (3.99 g/kg), the optimal dietary Thr (14.10g/kg) affected the immunocyte activation in the head kidney (HK) and spleen (SP) by downregulating the mRNA expressions of MHC-II and upregulating CD4 (not CD8), and it mediated the innate immune by enhancing the activities of lysozyme (LZ), acid phosphatase content of complement 3 (C3) and C4, increasing the mRNA abundances of hepcidin, liver expressed antimicrobial peptide-2A (LEAP-2A), LEAP-2B, β-defensin1, downregulating tumor necrosis factor α (TNF-α), IL-6, IL-1β, IL-12p35, IL-12p40, IL-17AF1, and IL-17D partly by attenuating RORγ1 transcriptional factor and nuclear factor kappa B p65 (NF-κBp65) signaling cascades [IKKβ/IκBα/NF-κBp65] and upregulating transforming growth factor β1 (TGF-β1), IL-4/13A, -4/13B, IL-10, and IL-22 partly by GATA-3. Besides these, the optimal dietary Thr regulated the adaptive immune by upregulating the mRNAs of immunoglobulin M (IgM) and IgZ (not IgD). Moreover, 2 mM Thr downregulated in vitro the mRNA abundances of colony stimulating factor-1, inducible nitric oxide synthase, mannose receptor 1, matrix metalloproteinase2 (MMP-2), and MMP-9 significantly (P < 0.05), indicating that Thr could attenuate the M1-type macrophages’ activation. Moreover, L-Thr downregulated the mRNA transcripts of TNF-α, IL-6, and IL-1β associated with impairing the SOCS1/STAT1 signaling and upregulated IL-10 and TGF-β1 partly by accentuating the SOCS3/STAT3 pathway. The above-mentioned observations suggested that Thr improved the immune status in the immune organs of fish by enhancing the immune defense and mediating the inflammation process. Finally, based on the immune indices of LZ activity in HK and C3 content in SP, the optimal Thr for immune enhancement in juvenile grass carp (9.53–53.43 g) was determined to be 15.70 g/kg diet (4.85 g/100 g protein) and 14.49 g/kg diet (4.47 g/100 g protein), respectively.

Gill Transcriptomic Responses to Toxin-producing Alga Prymnesium parvum in Rainbow Trout

The gill of teleost fish is a multifunctional organ involved in many physiological processes, including protection of the mucosal gill surface against pathogens and other environmental antigens by the gill-associated lymphoid tissue (GIALT). Climate change associated phenomena, such as increasing frequency and magnitude of harmful algal blooms (HABs) put extra strain on gill function, contributing to enhanced fish mortality and fish kills. However, the molecular basis of the HAB-induced gill injury remains largely unknown due to the lack of high-throughput transcriptomic studies performed on teleost fish in laboratory conditions. We used juvenile rainbow trout (Oncorhynchus mykiss) to investigate the transcriptomic responses of the gill tissue to two (high and low) sublethal densities of the toxin-producing alga Prymnesium parvum, in relation to non-exposed control fish. The exposure time to P. parvum (4–5 h) was sufficient to identify three different phenotypic responses among the exposed fish, enabling us to focus on the common gill transcriptomic responses to P. parvum that were independent of dose and phenotype. The inspection of common differentially expressed genes (DEGs), canonical pathways, upstream regulators and downstream effects pointed towards P. parvum-induced inflammatory response and gill inflammation driven by alterations of Acute Phase Response Signalling, IL-6 Signalling, IL-10 Signalling, Role of PKR in Interferon Induction and Antiviral Response, IL-8 Signalling and IL-17 Signalling pathways. While we could not determine if the inferred gill inflammation was progressing or resolving, our study clearly suggests that P. parvum blooms may contribute to the serious gill disorders in fish. By providing insights into the gill transcriptomic responses to toxin-producing P. parvum in teleost fish, our research opens new avenues for investigating how to monitor and mitigate toxicity of HABs before they become lethal.