Expression analysis of immune genes in goldfish (Carassius auratus) infected with the monogenean parasite Gyrodactylus kobayashii

Transcriptome Profile Analyses of Head Kidney in Roach (Rutilus rutilus), Common Bream (Abramis brama) and Their Hybrids: Does Infection by Monogenean Parasites in Freshwater Fish Reveal Differences in Fish Vigour among Parental Species and Their Hybrids?

Hybrid generations usually face either a heterosis advantage or a breakdown, that can be expressed by the level of parasite infection in hybrid hosts. Hybrids are less infected by parasites than parental species (especially F1 generations) or more infected than parental species (especially post-F1 generations). We performed the experiment with blood-feeding gill parasite Paradiplozoon homoion (Monogenea) infecting leuciscid species, Abramis brama and Rutilus rutilus, their F1 generation and two backcross generations. Backcross generations tended to be more parasitized than parental lines and the F1 generation. The number of differentially expressed genes (DEGs) was lower in F1 hybrids and higher in backcross hybrids when compared to each of the parental lines. The main groups of DEGs were shared among lines; however, A. brama and R. rutilus differed in some of the top gene ontology (GO) terms. DEG analyses revealed the role of heme binding and erythrocyte differentiation after infection by blood-feeding P. homoion. Two backcross generations shared some of the top GO terms, representing mostly downregulated genes associated with P. homoion infection. KEGG analysis revealed the importance of disease-associated pathways; the majority of them were shared by two backcross generations. Our study revealed the most pronounced DEGs associated with blood-feeding monogeneans in backcross hybrids, potentially (but not exclusively) explainable by hybrid breakdown. The lower DEGs reported in F1 hybrids being less parasitized than backcross hybrids is in line with the hybrid advantage.

Impact of high temperature, CO2 and parasitic infection on inflammation, immunodepression and programmed cell death in Colossoma macropomum at the transcriptional level

The production of tambaqui Colossoma macropomum has recently reached a milestone, being considered the main native species produced in South American continental waters. Despite the importance of this fish, its immunity is poorly understood, and global warming could pose severe risks to its health as increasing water temperature leads to an increase in the incidence of parasitic diseases. In an experimental context based on the high-emission scenario of the 5th Intergovernmental Panel on Climate Change (IPCC) report, we evaluated the synergistic effect of exposure to the extreme climate change scenario (RCP8.5) during two exposure periods (7 and 30 days) and two levels of parasitism by monogeneans (low and high). The goal was to understand how the tambaqui immune system will react to this challenge. To achieve this goal, we analyzed the expression of nine immunity-related genes (jak3, stat3, il-10, socs1, casp1, il-1β, tp53, bcl2, and hif-1α) in the spleen. Our main findings showed downregulation in the jak3/stat3 pathway, genes related to the control of inflammation and apoptosis, in addition to upregulation of proinflammatory genes and those related to pyroptosis during the first 7 days of exposure to the extreme climate scenario, also indicating a stage of immunodepression in these animals. After 30 days of exposure, all genes tended to return to similar levels in the current scenario, possibly due to the decrease in parasite load caused by chronic exposure to the extreme scenario. Our data strongly suggest that the increase in parasitism intensity caused by the extreme climate change scenario is responsible for disturbances in the host's immune system. However, more studies are needed to clarify this poorly understood cascade of events.

A bloody interaction: plasma proteomics reveals gilthead sea bream (Sparus aurata) impairment caused by Sparicotyle chrysophrii

Background

Sparicotylosis is an enzootic parasitic disease that is well established across the Mediterranean Sea. It is caused by the polyopisthocotylean monogenean Sparicotyle chrysophrii and affects the gills of gilthead sea bream (GSB; Sparus aurata). Current disease management, mitigation and treatment strategies are limited against sparicotylosis. To successfully develop more efficient therapeutic strategies against this disease, understanding which molecular mechanisms and metabolic pathways are altered in the host is critical. This study aims to elucidate how S. chrysophrii infection modulates the plasma proteome of GSB and to identify the main altered biological processes involved.

Methods

Experimental infections were conducted in a recirculating aquaculture system (RAS) in which naïve recipient GSB ([R]; 70 g; n = 50) were exposed to effluent water from S. chrysophrii-infected GSB (98 g; n = 50). An additional tank containing unexposed naïve fish (control [C]; 70 g; n = 50) was maintained in parallel, but with the open water flow disconnected from the RAS. Haematological and infection parameters from sampled C and R fish were recorded for 10 weeks. Plasma samples from R fish were categorised into three different groups according to their infection intensity, which was based on the number of worms fish⁻¹: low (L: 1–50), medium (51–100) and high (H: > 100). Five plasma samples from each category and five C samples were selected and subjected to a SWATH-MS proteome analysis. Additional assays on haemoglobin, cholesterol and the lytic activity of the alternative complement pathway were performed to validate the proteome analysis findings.

Results

The discriminant analysis of plasma protein abundance revealed a clear separation into three groups (H, M/L and C). A pathway analysis was performed with the differentially quantified proteins, indicating that the parasitic infection mainly affected pathways related to haemostasis, the immune system and lipid metabolism and transport. Twenty-two proteins were significantly correlated with infection intensity, highlighting the importance of apolipoproteins, globins and complement component 3. Validation assays of blood and plasma (haemoglobin, cholesterol and lytic activity of alternative complement pathway) confirmed these correlations.

Conclusions

Sparicotylosis profoundly alters the haemostasis, the innate immune system and the lipid metabolism and transport in GSB. This study gives a crucial global overview of the pathogenesis of sparicotylosis and highlights new targets for further research.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05441-1.

Skin transcriptomic analysis and immune-related gene expression of golden pompano (Trachinotus ovatus) after Amyloodinium ocellatum infection

Amyloodiniosis is a severe disease of marine and brackish water fish caused by Amyloodinium ocellatum. Golden pompano (Trachinotus ovatus) is often repeatedly infected by A. ocellatum, leading to extensive mortality. However, little is known about the immune response mechanisms of the T. ovatus following reinfection with A. ocellatum. In this study, an extensive analysis at the transcriptome level of T. ovatus skin was carried out at 24 h post-infection by A. ocellatum. During the transcriptomic analysis, 1367 differentially expressed genes (DEGs) in the skin of T. ovatus under A. ocellatum infection and control conditions were obtained. In Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotated analyses, the DEGs were significantly enriched in the immune-related pathways. To better understand the immune-related gene expression dynamics, a quantitative reverse transcription-polymerase chain reaction (RT-qPCR) was used to assess the primary and secondary infection groups of T. ovatus at different stages (3 h, 12 h, 24 h, 48 h and, 72 h post-infection) of infection with A.ocellatum. The results showed that innate immunity-related genes [interleukin (IL-8), chemokine ligand 3 (CCL3), toll-like receptor 7 (TLR7), and G-type lysosome (LZM g)] and adaptive immunity-related gene [major histocompatibility complex (MHC) alpha antigen I and MHC alpha antigen II] expression levels in the primary and secondary infection groups were significantly increased compared to the control group. The expression of MHC I and MHC II was more rapidly upregulated in the secondary infection group compared with the primary infection group after A.ocellatum infection. However, no significant differences of A.ocellatum load were observed in primary and secondary infection groups. In addition, the serum of the primary infection group had significantly higher concentrations of triglyceride (TG), higher alanine transaminase (ALT), aspartate transaminase (AST), and lactate dehydrogenase (LDH) activities than the control group. This study contributes to understanding the defense mechanisms in fish skin against ectoparasite infection.

Population growth of Gyrodactylus kobayashii in goldfish (Carassius auratus) associated with host density

Host density is a key regulatory factor in parasite transmission. The goldfish (Carassius auratus)-Gyrodactylus kobayashii model was used to investigate effects of host density on population growth of gyrodactylids. A donor fish infected by five gravid gyrodactylids was mixed with 11 parasite-free goldfish at five host densities. There was a significant positive correlation between host density and mean abundance of G. kobayashii throughout the 58-day experiment. During early infection (days 15–24), mean abundance in medium high (0.5 fish L−1) and high host density groups (1 and 2 fish L−1) was significantly higher than that in the low host density groups (0.125 and 0.25 fish L−1). At high host density, prevalence increased more rapidly, and the peak prevalence was higher. Fitting of an exponential growth model showed that the population growth rate of the parasite increased with host density. A hypothesis was proposed that higher host density contributed to increased reinfection of detached gyrodactylids. A reinfection experiment was designed to test this hypothesis. Both mean abundance and prevalence at a host density of 1 fish L−1 were significantly higher than those at 0.25 fish L−1 on days 1 and 3, which suggested that more reinfections of G. kobyashii occurred at the higher host density. Density-dependent transmission during the early infection was an important contributor of population growth of G. kobayashii, as well as density-dependent reinfection of the detached gyrodactylids.

Transmission dynamics of ectoparasitic gyrodactylids (Platyhelminthes, Monogenea): An integrative review

Parasite transmission is the ability of pathogens to move between hosts. As a key component of the interaction between hosts and parasites, it has crucial implications for the fitness of both. Here, we review the transmission dynamics of Gyrodactylus species, which are monogenean ectoparasites of teleost fishes and a prominent model for studies of parasite transmission. Particularly, we focus on the most studied host–parasite system within this genus: guppies, Poecilia reticulata, and G. turnbulli/G. bullatarudis. Through an integrative literature examination, we identify the main variables affecting Gyrodactylus spread between hosts, and the potential factors that enhance their transmission. Previous research indicates that Gyrodactylids spread when their current conditions are unsuitable. Transmission depends on abiotic factors like temperature, and biotic variables such as gyrodactylid biology, host heterogeneity, and their interaction. Variation in the degree of social contact between hosts and sexes might also result in distinct dynamics. Our review highlights a lack of mathematical models that could help predict the dynamics of gyrodactylids, and there is also a bias to study only a few species. Future research may usefully focus on how gyrodactylid reproductive traits and host heterogeneity promote transmission and should incorporate the feedbacks between host behaviour and parasite transmission.

Functional immunogenetic variation, rather than local adaptation, predicts ectoparasite infection intensity in a model fish species

Natural host populations differ in their susceptibility to infection by parasites, and these intrapopulation differences are still an incompletely understood component of host‐parasite dynamics. In this study, we used controlled infection experiments with wild‐caught guppies (Poecilia reticulata) and their ectoparasite Gyrodactylus turnbulli to investigate the roles of local adaptation and host genetic composition (immunogenetic and neutral) in explaining differences in susceptibility to infection. We found differences between our four study host populations that were consistent between two parasite source populations, with no indication of local adaptation by either host or parasite at two tested spatial scales. Greater values of host population genetic variability metrics broadly aligned with lower population mean infection intensity, with the best alignments associated with major histocompatibility complex (MHC) “supertypes”. Controlling for intrapopulation differences and potential inbreeding variance, we found a significant negative relationship between individual‐level functional MHC variability and infection: fish carrying more MHC supertypes experienced infections of lower severity, with limited evidence for supertype‐specific effects. We conclude that population‐level differences in host infection susceptibility probably reflect variation in parasite selective pressure and/or host evolutionary potential, underpinned by functional immunogenetic variation.

Molecular communication between the monogenea and fish immune system

Monogeneans parasitise mainly the outer structures of fish, such as the gills, fins, and skin, that is, tissues covered with a mucous layer. While attached by sclerotised structures to host's surface, monogeneans feed on its blood or epidermal cells and mucus. Besides being a rich source of nutrients, these tissues also contain humoral immune factors and immune cells, which are ready to launch defence mechanisms against the tegument or gastrointestinal tract of these invaders. The exploitation of hosts' resources by the Monogenea must, therefore, be accompanied by suppressive and immunomodulatory mechanisms which protect the parasites against attacks by host immune system. Elimination of hosts' cytotoxic molecules and evasion of host immune response is often mediated by proteins secreted by the parasites. The aim of this review is to summarise existing knowledge on fish immune responses against monogeneans. Results gleaned from experimental infections illustrate the various interactions between parasites and the innate and adaptive immune system of the fish. The involvement of monogenean molecules (mainly inhibitors of peptidases) in molecular communication with host immune system is discussed.

Mucosal immune response of Nile tilapia Oreochromis niloticus during Gyrodactylus cichlidarum infection

Monogenean Gyrodactylus cichlidarum can cause severe mortality of Nile tilapia (Oreochromis niloticus) fry. To date, reports about mucosal immunity of O. niloticus against this parasite have been rare. In order to explore the mucosal immunity of Nile tilapia against G. cichlidarum infection, the expressions of six adaptive immune-related genes and the contents of specific immunoglobulin IgM and IgT in the skin-associated lymphoid tissues (SALT) were dynamically analyzed after primary and secondary infections. The abundances of G. cichlidarum on the hosts after secondary infection were lower than those after primary parasite infection, which implied that hosts could initiate immune protection against G. cichlidarum reinfection to some degree. The transcription levels of TCR-β and CD4 genes in the skin tissue were significantly up-regulated after primary G. cichlidarum infection, while genes pIgR and IgT were only detected with significant up-regulations during secondary infection. With the exception of pIgR, which had remarkably higher expression in the fish with low parasite loads, all other genes studied tended to have higher mRNA level in the fish with higher parasite loads. The specific IgM content in the skin mucus increased significantly on the 2nd day after the primary exposure, higher than the corresponding value during the secondary exposure, and had significantly positive correlation with the parasite loads during the first parasite infection. These results manifested that acquired immune responses in the SALT of Nile tilapia participated in the resistance against G. cichlidarum infection, underscoring the involvement of mucosal immunity in fish against monogenean infection, and suggesting potential prophylactic treatment of gyrodactylid disease of tilapia.

Transcriptome analysis of goldfish (Carassius auratus) in response to Gyrodactylus kobayashii infection

Gyrodactylid monogeneans are widespread parasites of teleost fishes, and infection with these parasites results in high host morbidity and mortality in aquaculture. To comprehensively elucidate the immune mechanisms against Gyrodactylus kobayashii, the transcriptome profiles of goldfish (Carassius auratus) skin after challenge with G. kobayashii were first investigated using next-generation sequencing. Approximately 21 million clean reads per library were obtained, and the average percentage of these clean reads mapped to the reference genome was 82.25%. A total of 556 differentially expressed genes (DEGs), including 344 upregulated and 212 downregulated genes, were identified, and 380 DEGs were successfully annotated and assigned to 95 signaling pathways in Kyoto Encyclopedia of Genes and Genomes (KEGG). In addition, 14 pathways associated with immune response were identified mainly including mTOR signaling pathway, cytokine-cytokine receptor interaction, intestinal immune network for IgA production, toll-like receptor signaling pathway, and phagosome. Twelve genes were selected and validated using qRT-PCR. A similar trend of these genes between RNA-Seq and qRT-PCR was observed, indicating that RNA-Seq data was reliable. Besides, the ALP activity and NO content in serum were significantly higher in the infected goldfish compared with the non-infected goldfish. In summary, this study provides better understandings of immune defense mechanisms of goldfish against G. kobayashii, which will support future molecular research on gyrodactylids and facilitate the prevention and treatment of gyrodactylosis in aquaculture.

Unveiling the effect of dietary essential oils supplementation in Sparus aurata gills and its efficiency against the infestation by Sparicotyle chrysophrii

A microencapsulated feed additive composed by garlic, carvacrol and thymol essential oils (EOs) was evaluated regarding its protective effect in gills parasitized by Sparicotyle chrysophrii in Sparus aurata. A nutritional trial (65 days) followed by a cohabitation challenge with parasitized fish (39 days) were performed. Transcriptomic analysis by microarrays of gills of fish fed the EOs diet showed an up-regulation of genes related to biogenesis, vesicular transport and exocytosis, leukocyte-mediated immunity, oxidation–reduction and overall metabolism processes. The functional network obtained indicates a tissue-specific pro-inflammatory immune response arbitrated by degranulating acidophilic granulocytes, sustained by antioxidant and anti-inflammatory responses. The histochemical study of gills also showed an increase of carboxylate glycoproteins containing sialic acid in mucous and epithelial cells of fish fed the EOs diet, suggesting a mucosal defence mechanism through the modulation of mucin secretions. The outcomes of the in vivo challenge supported the transcriptomic results obtained from the nutritional trial, where a significant reduction of 78% in the abundance of S. chrysophrii total parasitation and a decrease in the prevalence of most parasitic developmental stages evaluated were observed in fish fed the EOs diet. These results suggest that the microencapsulation of garlic, carvacrol and thymol EOs could be considered an effective natural dietary strategy with antiparasitic properties against the ectoparasite S. chrysophrii.

RNA-Seq analysis of the guppy immune response against Gyrodactylus bullatarudis infection

Gyrodactylids are ubiquitous ectoparasites of teleost fish, but our understanding of the host immune response against them is fragmentary. Here, we used RNA-Seq to investigate genes involved in the primary response to infection with Gyrodactylus bullatarudis on the skin of guppies, Poecilia reticulata, an important evolutionary model, but also one of the most common fish in the global ornamental trade. Analysis of differentially expressed genes identified several immune-related categories, including IL-17 signalling pathway and Th17 cell differentiation, cytokine-cytokine receptor interaction, chemokine signalling pathway, NOD-like receptor signalling pathway, natural killer cell-mediated cytotoxicity and pathways involved in antigen recognition, processing and presentation. Components of both the innate and the adaptive immune responses play a role in response to gyrodactylid infection. Genes involved in IL-17/Th17 response were particularly enriched among differentially expressed genes, suggesting a significant role for this pathway in fish responses to ectoparasites. Our results revealed a sizable list of genes potentially involved in the teleost-gyrodactylid immune response.

Extreme climate scenario and parasitism affect the Amazonian fish Colossoma macropomum

Global warming caused by greenhouse gases accumulation, in particular carbon dioxide, is the major current environmental challenge, as it will affect all life forms over the next decades. Aquaculture, a fast growing food production sector, is already facing the effects of global warming. The fish immune system is expected to be especially affected by increased temperature and carbon dioxide, mainly when associated to infectious diseases outbreaks. Here, we experimentally investigated the associated effects of an extreme climate scenario and two levels of monogenean parasitism on the hematological and immunological response of an important Amazon fish for continental aquaculture: Colossoma macropomum. Individuals of C. macropomum with low and high levels of parasitism were exposed to current and extreme climate scenarios (4.5 °C and 900 ppm CO₂ above current levels). We characterized their hematological profile using classical methods, and their immune-related gills gene expression (HSP70, IL-1β and IL-10) using quantitative real-time polymerase chain reaction (qPCR). After 7 days of exposure, we observed that exposure to extreme climate scenario caused rapid increase of parasitism intensity and likely acute inflammation, indicated by the higher expression of HSP70 and IL-1β. The IL-10 gene was downregulated in both groups exposed to extreme climate scenario, contrasting with animals exposed to current scenario. Thus, in the current scenario, the parasitized tambaqui showed a response to the tissue damage, which was not observed in the animals exposed to the extreme scenario.

Immune response of teleost fish to helminth parasite infection

Abstract Fish immune systems have become the subject of several studies due to the growing development of aquaculture and fisheries, and the demand for healthy produce for human consumption. Among the parasites responsible for diseases in fish farming, helminths stand out because they cause infections in farmed fish and decrease food conversion, zootechnical performance and meat quality. In the present review, the components that participate in the innate and adaptive immune responses of teleost fish that have so far been described are presented in order to summarize the defenses that these hosts have recourse to, in combating different groups of helminth parasites.

Heterogeneity in helminth infections: factors influencing aggregation in a simple host-parasite system

The almost universally-occurring aggregated distributions of helminth burdens in host populations have major significance for parasite population ecology and evolutionary biology, but the mechanisms generating heterogeneity remain poorly understood. For the direct life cycle monogenean Discocotyle sagittata infecting rainbow trout, Oncorhynchus mykiss, variables potentially influencing aggregation can be analysed individually. This study was based at a fish farm where every host individual becomes infected by D. sagittata during each annual transmission period. Worm burdens were examined in one trout population maintained in isolation for 9 years, exposed to self-contained transmission. After this year-on-year recruitment, prevalence was 100% with intensities 10-2628, mean 576, worms per host. Parasite distribution, amongst hosts with the same age and environmental experience, was highly aggregated with variance to mean ratio 834 and negative binomial parameter, k, 0.64. The most heavily infected 20% of fish carried around 80% of the total adult parasite population. Aggregation develops within the first weeks post-infection; hosts typically carried intensities of successive age-specific cohorts that were consistent for that individual, such that heavily-infected individuals carried high numbers of all parasite age classes. Results suggest that host factors alone, operating post-infection, are sufficient to generate strongly overdispersed parasite distributions, rather than heterogeneity in exposure and initial invasion.

Acting locally - affecting globally: RNA sequencing of gilthead sea bream with a mild Sparicotyle chrysophrii infection reveals effects on apoptosis, immune and hypoxia related genes

BACKGROUND

Monogenean flatworms are the main fish ectoparasites inflicting serious economic losses in aquaculture. The polyopisthocotylean Sparicotyle chrysophrii parasitizes the gills of gilthead sea bream (GSB, Sparus aurata) causing anaemia, lamellae fusion and sloughing of epithelial cells, with the consequent hypoxia, emaciation, lethargy and mortality. Currently no preventive or curative measures against this disease exist and therefore information on the host-parasite interaction is crucial to find mitigation solutions for sparicotylosis. The knowledge about gene regulation in monogenean-host models mostly comes from freshwater monopysthocotyleans and almost nothing is known about polyopisthocotyleans. The current study aims to decipher the host response at local (gills) and systemic (spleen, liver) levels in farmed GSB with a mild natural S. chrysophrii infection by transcriptomic analysis.

RESULTS

Using Illumina RNA sequencing and transcriptomic analysis, a total of 2581 differentially expressed transcripts were identified in infected fish when compared to uninfected controls. Gill tissues in contact with the parasite (P gills) displayed regulation of fewer genes (700) than gill portions not in contact with the parasite (NP gills) (1235), most likely due to a local silencing effect of the parasite. The systemic reaction in the spleen was much higher than that at the parasite attachment site (local) (1240), and higher than in liver (334). NP gills displayed a strong enrichment of genes mainly related to immune response and apoptosis. Processes such as apoptosis, inflammation and cell proliferation dominated gills, whereas inhibition of apoptosis, autophagy, platelet activation, signalling and aggregation, and inflammasome were observed in spleen. Proteasome markers were increased in all tissues, whereas hypoxia-related genes were down-regulated in gills and spleen.

CONCLUSIONS

Contrasting forces seem to be acting at local and systemic levels. The splenic down-regulation could be part of a hypometabolic response, to counteract the hypoxia induced by the parasite damage to the gills and to concentrate the energy on defence and repair responses. Alternatively, it can be also interpreted as the often observed action of helminths to modify host immunity in its own interest. These results provide the first toolkit for future studies towards understanding and management of this parasitosis.

Cytokine gene expression profiles in goldfish (Carassius auratus) during Gyrodactylus kobayashii infection

Monogeneans of the genus Gyrodactylus are well-known pathogens causing huge mortalities in wild and cultured fish. Cytokine expression is one of most important host defense mechanisms against parasite infections. In this study, the expression pattern of the key pro-inflammatory (IL-1β, IL-8, IFN-γ, TNF-α, IL-12 and iNOS) and anti-inflammatory cytokine genes (IL-10, TGFβ and IL-4) of Gyrodactylus kobayashii infected goldfish (Carassius auratus) were determined by real-time quantitative PCR analysis. Our results showed that G. kobayashii infection caused increased expression of the pro-inflammatory cytokines including IFN-γ, TNF-α and iNOS in all detected tissues throughout the infection period. Among these genes, iNOS has the highest transcript level accompanied with increased nitric oxide (NO) concentration in the serum of all infected goldfish. The mRNA level of IL-1β in the liver, spleen and head kidney was significantly up-regulated during the early stage of infection (days 2-8). While high expression level of IL-8 and IL-12 was observed during the elimination phase of infection (days 10-14). As for anti-inflammatory cytokines, the expression profiles of IL-10 were distinct from those of TGF-β and IL-4. Specifically, the mRNA level of IL-10 did not increase in the spleen and head kidney during the early stage of infection, while increased expression of TGF-β and IL-4 were likewise seen. Besides, all infected fish had significantly higher complement C3 but lower IgM levels than the non-infected fish. The results provide insights into the interaction between gyrodactylids and the fish host, and indicate that systemic cytokine responses are critical for controlling parasite infection in fish.