Next Article in Journal
MicroRNAs Associated with Metformin Treatment in the Diabetes Prevention Program
Next Article in Special Issue
Identification of Marker Genes in Infectious Diseases from ScRNA-seq Data Using Interpretable Machine Learning
Previous Article in Journal
Transmission of Oyster Mushroom Spherical Virus to Progeny via Basidiospores and Horizontally to a New Host Pleurotus floridanus
Previous Article in Special Issue
Applicability of a Chemiluminescence Immunoassay to Screen Postmortem Bile Specimens and Its Agreement with Confirmation Analysis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 11
10.3390/ijms25115683
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function

by
Keren Maor-Landaw
1,
Margarita Smirnov
2 and
Tamar Lotan
1,*
1
Marine Biology Department, The Leon H. Charney School of Marine Sciences, University of Haifa, Mt. Carmel, Haifa 3103301, Israel
2
Central Fish Health Laboratory, Department of Fisheries and Aquaculture, Ministry of Agriculture and Rural Development, Nir David 1080300, Israel
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(11), 5683; https://doi.org/10.3390/ijms25115683
Submission received: 25 April 2024 / Revised: 19 May 2024 / Accepted: 22 May 2024 / Published: 23 May 2024
(This article belongs to the Special Issue Targeted Therapy for Immune Diseases)
Browse Figures
Versions Notes

Abstract

:
Myxozoa, a unique group of obligate endoparasites within the phylum Cnidaria, can cause emerging diseases in wild and cultured fish populations. Recently, the myxozoan Myxobolus bejeranoi has been identified as a prevalent pathogen infecting the gills of cultured hybrid tilapia, leading to systemic immune suppression and considerable mortality. Here, we employed a proteomic approach to examine the impact of M. bejeranoi infection on fish gills, focusing on the structure of the granulomata, or cyst, formed around the proliferating parasite to prevent its spread to surrounding tissue. Enrichment analysis showed increased immune response and oxidative stress in infected gill tissue, most markedly in the cyst’s wall. The intense immune reaction included a consortium of endopeptidase inhibitors, potentially combating the myxozoan arsenal of secreted proteases. Analysis of the cyst’s proteome and histology staining indicated that keratin intermediate filaments contribute to its structural rigidity. Moreover, we uncovered skin-specific proteins, including a grainyhead-like transcription factor and a teleost-specific S100 calcium-binding protein that may play a role in epithelial morphogenesis and cysts formation. These findings deepen our understanding of the proteomic elements that grant the cyst its distinctive nature at the critical interface between the fish host and myxozoan parasite.

1. Introduction

Tilapia, belonging to the Cichlidae family, ranks as the second most widely cultivated fish globally [1] and comprises 60% of the total freshwater fish production in Israel [2]. Intensive aquaculture systems of earthen ponds housing all-male hybrids of Oreochromis niloticus (Nile tilapia) females and O. aureus (Jordan/blue tilapia) males have become widespread worldwide [1,3,4], including in Israel [3,5]. Hybrid tilapia have gained popularity due to their rapid growth, adaptability to various environmental conditions, stress resistance, and efficient reproduction in captivity [1]. However, in Israeli fish ponds, hybrid tilapia are severely affected by the myxozoan parasite M. bejeranoi, which infects these fish at a prevalence exceeding 80% [6].
Myxozoa is a vast group of microscopic obligate endoparasites within the phylum Cnidaria [7]. Myxozoans impact both wild and farmed fish populations, causing diseases such as whirling disease and proliferative kidney disease [8]. The complex myxozoan life cycle involves two hosts; a vertebrate, mostly fish, and an invertebrate, mostly a worm [9,10]. Transmission between hosts is achieved by two distinct types of waterborne spores termed actinospores and myxospores [11,12].
Typically, the first contact of myxozoans with fish hosts is through mucosal surfaces, where they might be challenged by various cell types [13,14]. Teleost fish are armed with both innate and adaptive immunity [15,16,17], and they share a large repertoire of immune system cells and molecules with mammals [18,19]. After passing mucosal and epithelial barriers, the myxozoan parasite travels, typically through the bloodstream, to its specific target tissue where it proliferates [20,21,22]. Once the parasite is at the target tissue, the host activates immune mechanisms including immunoactivating and immunosuppressive cytokines [18,23,24]. Another defensive strategy is the formation of a granuloma, or cyst, around the myxozoan plasmodia. The cyst wall encapsulates the parasite by connective and epithelioid tissue layers originating from the host [25], thereby isolating it and preventing its dispersal to surrounding tissues [18,26].
M. bejeranoi sporulation site is the gill filament base [27], where it spills its infective sporoplasm. Myxospores are formed in fish within plasmodia, where multicellular proliferative stages occur. The plasmodia undergo intense cell differentiation to generate sporogonic cell stages, culminating in the production of mature spores [28]. M. bejeranoi is a highly efficient parasite that is capable of rapid proliferation within its host. Upon infection, M. bejeranoi cells express genes, such as histones, that facilitate rapid cell divisions and proper supply of the required energy [29,30]. Further, similar to other parasites [29,31,32], it employs tactics, such as an increase in calreticulin and an arsenal of proteases, to invade host tissues and eventually shut down its immune system [6].
In a recent transcriptomic study, we characterized the temporal progression of M. bejeranoi infection and the immune response of hybrid tilapia. While in the gills, the immune system was triggered at a severe infection stage, in the immuno-organs, the head, kidney, and spleen, an immune suppression was observed [6]. Our recent investigations of the interactions between tilapia fish and M. bejeranoi have shown that the infection can occur at very early life stages [33] and that the systemic immune suppression induced by the parasite renders the host vulnerable to other opportunistic pathogens [6]. This might be detrimental to tilapia health and, thereby, have a high economic impact on commercial fish farms. Therefore, we investigated the effects of M. bejeranoi infection on fish gills using a proteomic approach, reflecting the tissue’s functional status and its dynamic responses to infection, with a particular focus on the cyst walls. Our findings demonstrate a significant impact on the structure of the gills, along with increased oxidative stress and a local up-regulation of the immune response. Additionally, within the cysts, a robust immune reaction in the form of endopeptidase inhibitor activity was evident. Furthermore, we revealed that the structure of the cysts is composed of specific proteins that support a network of keratin intermediate filaments, which provide the cysts with their rigid structure. These findings shed light on fish defense mechanisms as well as on the critical physical interface between the host and its parasite.

2. Results

To gain insight into the effect of M. bejeranoi infection on the gills of hybrid tilapia, we conducted proteomic profiling of three sample types: gill tissue from infected and control uninfected fish and cysts dissected from infected fish. Our analysis identified 2446 protein groups, 1659 of which were shared between the three sample types (Figure 1, Supplementary Table S1).
To identify proteomic changes in gill tissue following infection with M. bejeranoi, four biological replicates from control and infected gills were analyzed. Results showed 324 proteins that were significantly differentially expressed between infected and healthy gills, of which 212 were up-regulated and 112 were down-regulated (Welch’s t-test, p < 0.05) (Supplementary Table S1). To the list of up-regulated proteins, we added 102 proteins that were expressed exclusively in infected tissue, whereas to the list of down-regulated proteins, we added 387 proteins that were found only in control samples (Figure 1).
Enrichment analysis of up-regulated proteins in infected vs. control gills showed changes in proteins associated with oxidative stress and glutathione metabolism (Figure 2, Supplementary Table S2), among which superoxide dismutase 2 and glutathione S-transferase scored high fold change values (3.42 and 3.13, respectively) (Figure 3, Supplementary Table S1). Other up-regulated clusters contained actin and cytoskeleton and closely related proteins associated with the GO term of tight junction, RNA splicing and transport, and translation (Figure 2, Supplementary Table S2). Following infection with M. bejeranoi, proteins related to the fish immune system were increased, most prominently Scavenger Receptor Cysteine-Rich (SRCR) domain-containing protein (m130/cd163) (fold change 6.12) and pentraxin 3 (fold change 4.38) (Figure 3, Supplementary Table S1). The immune response included a notable group of endopeptidase inhibitors, some belonging to the serpin superfamily (Figure 2, Supplementary Table S2).
Enrichment analysis of down-regulated proteins revealed that those related to cellular respiration formed a distinct cluster apart from proteins responsible for oxygen binding (globins) (Figure 2, Supplementary Table S2). Two of the latter category, hemoglobin subunit beta-1 and hemoglobin subunit alpha, had remarkably low fold change values of −68.05 and −35.25, respectively (Figure 3). The next two most extremely down-regulated proteins were two collagen types, alpha-1(I) chain-like (collagen α1) and type I alpha 2 chain (Col1α2), with fold change values of −22.69 and −20.69, respectively. Other down-regulated processes were lysosome, phagosome, and protein transport with a high enrichment score of over 9 (Figure 2).
To gain insights into the cellular processes occurring at the fish–myxozoan interface, we analyzed the proteomic profile of cysts from infected fish. We sought to identify factors that contribute to the structural integrity of the rigid cysts, isolating their contents from the surrounding fish tissues (Figure 4A). To achieve this goal, we performed proteomic profiling on manually dissected cysts and compared them with the proteomes of gills from healthy, unexposed fish. Two of the analyzed four biological replicates of cyst samples were partial, and therefore, only the other two were used in the analysis. Among the identified 2075 fish protein groups, 179 groups were exclusively expressed in cysts, whereas 337 groups were exclusively expressed in healthy gills (Figure 1). Since gills and cysts are distinct tissues, we employed iBAQ values for the proteomic comparison. These values represent the relative abundance of proteins in each sample [34]. We found 367 protein groups that were significantly up-regulated and 584 groups that were significantly down-regulated in the cyst compared to control gills (Supplementary Table S1). A similar comparison between cysts and control gills was performed using STRING enrichment analysis. Several GO terms that were found for infected gills were also observed in the cyst analysis (Figure 4B, Supplementary Table S3), including oxidoreductase activity and a larger group of endopeptidase inhibitors. Table 1 details the up-regulated proteins within the prominent group of endopeptidase inhibitors, as well as their expression in cysts versus control gills. Similar to infected gills, cysts were also enriched in oxidative-stress-related proteins, with prominent representatives such as glutathione S-transferase and thioredoxin, which had among the highest fold change values (12.12 and 11.13, respectively) (Supplementary Table S1). A distinct cluster of endocytosis and glycolysis, along with proteasome, and catabolic processes, with nearly a 10-fold enrichment, was evident in up-regulated cysts compared to the control. Other highly expressed cyst proteins were the calcium-binding EF-hand domain-containing proteins, ictacalcin and calmodulin 2, with fold change values of 19.63 and 14.27, respectively (Supplementary Table S1). Notably, keratin 98, which is an intermediate filament protein, had a high fold change value of 9.36 in cysts. The pattern of the enriched down-regulated categories in cysts resembled the one observed in infected gills, with oxidative phosphorylation, extracellular matrix, and protein folding. The leading enriched categories were collagen-containing extracellular matrix and focal adhesion (Figure 4C, Supplementary Table S3).
Next, we mined our proteomic data to identify candidates for structural proteins that could reinforce the barrier formed by the cyst tissue. We therefore focused on differentially expressed proteins in cysts compared to control gill samples, as well as proteins exclusive to cysts. Two proteins involved with cell–cell adhesion were found in cyst proteome and were missing from control samples: cadherin 2 and claudin I, of adherens junction and tight junction, respectively (Supplementary Table S1). We identified three keratin proteins: keratin 98, si:dkey-222n6.2, and keratin, type I cytoskeletal 50 kDa. To visualize keratin localization within the cysts, we applied three different keratin stains, Hematoxylin and Eosin, Ayoub-Shklar, and Dane-Herman stain, all of which are known to effectively stain keratin [35,36]. Our results demonstrate that the cysts’ walls are indeed composed of keratin (Figure 5A–C). Additionally, alongside the identified keratins, we identified two proteins specifically expressed in cysts: sciellin, a precursor of cornified envelope [37] and grainyhead-like transcription factor 1, which plays a role in epidermal barrier structure [38,39,40]. To gain insight into their functions within cysts, we proceeded to explore the protein interactions of these keratins. For this purpose, we generated a protein interaction network of the top 2000 expressed proteins identified within the cysts. A derivative network of keratins and sciellin ‘first-neighbors’ was further drafted and is shown in Figure 5D. The three keratins that were activated in cysts were connected to other cyst-expressed keratins, heat shock proteins, nucleoporins, and several ubiquitins. Sciellin was related to envoplakin, periplakin, and plakoglobin, which are characteristic of cornified epithelia [41] (Figure 5D). Ran binding protein 2 (RanBP2) appeared to be a central node in the network, exhibiting multiple connections. Taken together, the identified protein network indicated a rigid cornified structure of the cyst’s wall.

3. Discussion

A common pathological reaction of fish hosts to myxozoan infection is the formation of a granulomatous complex. This cellular barrier encapsulates and isolates the parasite, thereby preventing its spreading to surrounding tissues [18,26]. During the interaction between hybrid tilapia and M. bejeranoi, cysts are formed around developing plasmodia in the striated muscle of the gill filament, where infection takes place [27]. While the origin of the cyst tissue is debated [42], a previous study has demonstrated that the cyst walls surrounding the myxozoan are of fish origin [25].
Confined within the cyst walls, a chronic inflammatory reaction involves a complex and diverse array of leukocytes, macrophages, mast cells, and epithelioid cells [18,43]. Upon infection with M. bejeranoi, the entire gill tissue exhibits an immune response; however, the intensity of this is higher within the cysts. The SRCR domain-containing protein CD163, which we identified as the most highly expressed immune response protein, exhibited increased expression in cysts. Little is known about SRCR-domain containing proteins in teleosts [44,45]; nonetheless, the essential role of CD163 in the immune system has been demonstrated in humans [46,47]. CD163 is expressed on the surface of macrophages, where its binding to haptoglobin/hemoglobin complexes stimulates macrophage digestion. This process is essential for eliminating free hemoglobin, preventing oxidative stress and cells injury. Additionally, CD163 binding triggers an anti-inflammatory reaction and was thus suggested as a biomarker for inflammation [46,47,48]. Free hemoglobin, the product of hemolysis, can be induced by various pathogens [49], which was recently reported as a feeding behavior of the myxozoan Sphaerospora molnari [50,51]. Upon infection, M. bejeranoi deploys an array of activated proteases, including cathepsin D and cathepsin L [29], which may act to degrade hemoglobin, as in other parasites [32,52]. This aligns with our findings, indicating decrease in globin proteins and an increase in endocytosis, the process in which macrophages engulf their targets [47].
Myxozoan parasites exploit their host and migrate through its tissues using a battery of secreted proteolytic enzymes [29,32,52,53,54,55]. Thus, a vital component in the host immune arsenal is protease inhibitors [56]. This host–parasite interaction culminates at the cyst walls, further stressing the importance of the cyst in the host strategy to contain the pathogen. Nevertheless, this cellular barrier does not offer absolute protection against the devastating effects of the parasite. Our results reveal an impact on the entire gill tissue manifested, for instance, in pronounced oxidative stress. Moreover, the effects of a progressing infection are systemic, impacting the hematopoietic organs head, kidney, and spleen, leading to immune suppression and high mortality [6].
As sporogenesis progresses, M. bejeranoi cells continue to proliferate inside the cysts and the number of oxygen-consuming cells increases, leading to hypoxic conditions [29]. Nevertheless, we did not observe signs of angiogenesis in the host connective tissue, such as in the case of Myxobolus pendula, where the host vascular network was suggested to provide the parasite with nutrients and oxygen [42]. This may be related to the changes we reported in the fish cellular respiration.
One of our main questions was related to the protein elements that contribute to the integrity and rigidity of the cyst. Histopathological studies of granulomata forming in fish following exposure to various pathogens, including myxozoans, describe layers of fibroblasts intermingled with collagen fibers and densely packed epithelioid cells [25,42,43,57,58,59,60]. Our proteomic results did not confirm the histological evidence for the presence of collagen. Rather, we observed a substantial down-regulation of collagen in both infected tissue and cysts. This was accompanied by alterations in extracellular matrix organization, either due to the overall stress response or to enable cyst expansion.
Conversely, the elevated levels of keratins we observed align with myxozoan studies that employed anti-keratin antibodies for immunostaining in epidermal and epithelioid cell layers [25,43]. Keratins are a vast family of intermediate filaments that can aggregate into bundles of up to 12 nm in diameter [61], allowing them to endure mechanical forces from neighboring cells and uphold the integrity of the cytoskeletal [62]. Further, there are indications that complexes such as those consisting of cadherin-2 and claudin, which were up-regulated in cysts, play a crucial role in this structure of mechanical stability [63]. The stringent cysts that form in M. bejeranoi-infected fish physically press on the local nerves and muscle tissue of the gill filament base [6].
Sciellin, which was expressed only in cysts, is a protein regulator of keratinocytes differentiation and is found in cornified tissue [64,65]. Evoplakin and periplakin, that are apparent in the network, are recruited during this differentiation as well and are associated with keratin filament scaffold formation [41]. The transcription factor grainyhead like 1, which was specific to the cysts, plays a role in epithelial morphogenesis and wound repair [38,39,40]. In gilthead seabream, grainyhead transcription factor 1 along with keratin, are considered markers for skin regeneration [39]. Additionally, keratins are known to be regulated by post-translational modifications, including sumoylation [66]. Possibly RanBP2, which has a SUMO1 E3-like activity [67] and was highly connected to the keratin network in the cyst, plays a role in mediating the assembly of newly formed keratin bundles [68].
Another highly expressed protein in cyst walls was ictacalcin, also known as S100I [69]. This new member in the S100 calcium-binding protein family is exclusively found in teleosts [69]. Ictacalcin is considered to be skin-specific in fish, sometimes related to keratinocytes [70,71]. Its levels were up-regulated in fish infected with amoeba [72] and sea lice [73] and were down-regulated after vaccination with attenuated bacteria [74]. However, S100 proteins are poorly understood in fish. In mammals, they play roles in calcium signaling and inflammatory processes [69,75] and are shown to associate with dental and skin cysts [76,77,78]. Nevertheless, the function of ictacalcin in cysts of hybrid tilapia following M. bejeranoi infection remains unclear.
To conclude, our proteomic results indicate that fish cells construct the cyst walls, reinforced with keratin intermediate filaments, producing a rigid structure. At this cellular barrier between the rapidly proliferating myxozoan and its host, the interplay between the fish defense mechanisms and the parasite’s efforts to evade and overcome them is at its highest. However, although the parasite is confined within the cyst, this barrier is often breached, and the infection affects the entire gill tissue, leading to cellular respiratory and oxidative stress. Eventually, the impact will reach the immune organs, resulting in an immune-suppressed fish.

4. Materials and Methods

4.1. Gill Isolation

Exposure of naïve hybrid tilapia (O. niloticus × O. aureus) to M. bejeranoi-infected fish pond was conducted on August 2021, as we previously reported [29]. Briefly, healthy hybrid tilapia fish with a mean weight of 2.98 g were exposed to the pond water for 24 h using confined cages of ~100 L. The mean water temperature, which was recorded constantly during the experiment using a temperature data logger (HOBO), was 30.25 °C. The fish were then relocated to indoor tanks at the Central Fish Health Laboratory in Nir David. Tanks had a flow-through system with dechlorinated tap water at a temperature of ~25 °C. Fish were sampled at 10- and 20-days post-exposure. Following euthanasia using 1 mL/L 2-phenoxyethanol, gill tissue was dissected from control and infected fish, snap-frozen in liquid nitrogen, and kept at −80 °C until used.

4.2. Cyst Isolation

Intact cysts were dissected from infected gills under a Zeiss Discovery V8 binocular (Zeiss, Jena, Germany), using sterilized syringe needles. About 20 isolated cysts from different fish were pooled together into an Eppendorf tube in four biological replicates, snap-frozen in liquid nitrogen and kept at −80 °C until protein extraction.

4.3. DNA and Protein Extractions

Four frozen gill lamellae from control or infected fish were dissected under a Zeiss Discovery V8 binocular. Half of the lamellae were transferred into TRIzol reagent (Thermo Scientific, Waltham, MA, USA) for DNA extraction, while the remaining half were placed into a protein extraction buffer containing 0.5% deoxycholate, 10 mM Tris (pH 7.0), and 2 mM EDTA. DNA was extracted according to TRIzol reagent manufacturer’s instructions and as previously reported [6]. Briefly, the tissue was lysed in TRIzol reagent, and the DNA phase was separated using chloroform. The concentration of DNA was measured using a NanoDrop 2000c spectrophotometer (Thermo Scientific, Waltham, MA, USA).The presence of M. bejeranoi infection in the fish gills was determined by qPCR using specific primers targeted to amplify the M. bejeranoi small subunit ribosomal RNA gene (SSU rDNA), as previously reported [6]. Protein samples from gills and isolated cysts were homogenized in protein extraction buffer using 3 mm glass beads (CS Chemicals Ltd., Ahmedabad, India) in a TissueLyser II (Qiagen, Hilden, Germany) for 3 min at 30 Hz. Samples were then centrifuged, and the protein supernatants were kept at −80 °C until use. For proteolysis, the samples were reduced with 3 mM DTT (54 °C for 45 min), modified with 9 mM iodoacetamide in 400 mM ammonium bicarbonate (in the dark at room temperature for 30 min) and digested with modified trypsin (Promega) at a 1:50 enzyme-to-substrate ratio overnight at 37 °C. A second trypsinization was performed for 4 h. The resulting peptides were acidified with 1% formic acid and precipitated, whereas the peptides in the supernatant were desalted using C18 Stage Tips, dried and re-suspended in 0.1% formic acid.

4.4. Mass Spectrometry Analysis

The peptides were resolved by reverse-phase chromatography on 0.075 × 180 mm fused silica capillaries (J&W) packed with Reprosil reversed phase material (Dr Maisch GmbH, Germany). The peptides were eluted with the following concentrations of acetonitrile with 0.1% of formic acid in water at a flow rate of 0.15 μL/min: a linear 180 min gradient of 5 to 28%, followed by a 15 min gradient of 28 to 95% and 25 min at 95% acetonitrile. Mass spectrometry was performed by Q Executive HFX mass spectrometer (Thermo) in a positive mode (m/z 350–1200, resolution 120,000 for MS1 and 15,000 for MS2) using repetitively full MS scan. This was followed by collision-induced dissociation (HCD, at 27 normalized collision energy) of the 30 most dominant ions (>1 charges) selected from the first MS scan. The AGC settings were 3 × 106 for the full MS and 1 × 105 for the MS/MS scans. A dynamic exclusion list was enabled with an exclusion duration of 20 s. The mass spectrometry proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository [79] with the dataset identifier PXD051329.

4.5. Gills and Cyst Analysis

The mass spectrometry data were analyzed using the MaxQuant software 2.1.3.0 [80] for peak picking and identification. We used the Andromeda search engine against the O. niloticus databases (Ensembl release 110; https://ftp.ensembl.org/pub/release-110/fasta/oreochromis_niloticus (accessed on 12 September 2023)), with mass tolerance of 6 ppm for the precursor masses and the fragment ions. Oxidation on methionine and protein N-terminus acetylation were accepted as variable modifications, and carbamidomethyl on cysteine was accepted as static modifications. Minimal peptide length was set to seven amino acids, and a maximum of two miscleavages were allowed. The data were quantified by label-free analysis using the same software. False discovery rates (FDRs) for peptide and protein levels were filtered to 1% using the target-decoy strategy. Protein tables were filtered to eliminate the identifications from the reverse database and common contaminants. Protein groups were filtered according to the number of unique peptides in each group, with a minimum of one unique peptide in all replicates. Downstream proteomic analyses were performed using Perseus 2.0.11 [81,82]. LFQ intensity values were utilized for the comparison of infected versus control gill samples using Welch’s t-test in Perseus 2.0.11 [81,82]. Percent iBAQ values were calculated for cyst and control gill protein groups and further used for Welch’s t-test in Perseus 2.0.11 [81,82]. STRING v12.0 [83] was utilized for enrichment analysis and to produce protein interaction networks of significantly and uniquely expressed protein groups (p < 0.05) between infected and control gills and between cysts and control gills. An enrichment score was calculated for each group of GO/KEGG/STRING cluster terms, which was defined as the minus log of the geometric mean of all p-values of the GO categories within the group [84]. A STRING protein interaction network of cysts-expressed proteins was constructed using the top 2000 expressed proteins identified within the cysts. A derivative network of keratins and sciellin ‘first-neighbors’ was further drafted and was graphically edited in Cytoscape [85]. A Venn diagram representing the filtered protein groups was drawn using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/, accessed on 12 September 2023).

4.6. Keratin Staining

Gills dissected from infected gills were fixed in 10% neutral-buffered formalin, dehydrated in a graded ethanol series, and embedded in 2-hydroxyethyl methacrylate. Then, 3 μm thick sections were generated using a Leica RM 2245 microtome (Leica Biosystems, Nussloch, Germany). Sections were stained with haematoxylin and eosin, Ayoub-Shklar, and Dane-Herman methods [35,36].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25115683/s1.

Author Contributions

Conceptualization, K.M.-L. and T.L.; methodology, K.M.-L., M.S. and T.L.; formal analysis, K.M.-L. and T.L.; writing—original draft preparation, K.M.-L. and T.L.; writing—review and editing, K.M.-L., M.S. and T.L.; supervision, T.L.; funding acquisition, T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Israel Science Foundation, grant number 297/19. K.M.-L. was funded by the University of Haifa post-doc fellowship.

Institutional Review Board Statement

The animal study protocol was approved by the Institu-tional Ethics Committee of University of Haifa (permit 652/19).

Data Availability Statement

The data supporting the findings of this study are presented in the main text and its additional files. The mass spectrometry proteomic data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD051329.

Acknowledgments

We thank the Bioinformatics Facility at the University of Haifa and the Smoler Proteomics Center at the Technion for their assistance in the bioinformatics analysis and proteomic profiling.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. El-Sayed, A.-F.M. (Ed.) Tilapia Culture, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2020; ISBN 9780128165096. [Google Scholar]
  2. Hulata, G. Tilapias—Biology and Culture; Ma’arechet Publishing House: Kibbutz Dalia, Israel, 2014. [Google Scholar]
  3. Milstein, A.; Zoran, M.; Kochba, M.; Avnimelech, Y. Effect of different management practices on water quality of intensive tilapia culture systems in Israel. Aquac. Int. 2001, 9, 133–152. [Google Scholar] [CrossRef]
  4. Eknath, A.E.; Hulata, G. Use and exchange of genetic resources of Nile tilapia (Oreochromis niloticus). Rev. Aquac. 2009, 1, 197–213. [Google Scholar] [CrossRef]
  5. Miao, W.; Wang, W. Trends of aquaculture production and trade: Carp, tilapia, and shrimp. Asian Fish. Sci. 2020, 33, 1–10. [Google Scholar] [CrossRef]
  6. Maor-Landaw, K.; Smirnov, M.; Brekhman, V.; Ofek-Lalzar, M.; Yahav, T.; Lotan, T. Infection by the Parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) suppresses the immune system of hybrid tilapia. Microorganisms 2022, 10, 1893. [Google Scholar] [CrossRef]
  7. Atkinson, S.D.; Bartholomew, J.L.; Lotan, T. Myxozoans: Ancient metazoan parasites find a home in phylum Cnidaria. Zoology 2018, 129, 66–68. [Google Scholar] [CrossRef] [PubMed]
  8. Okamura, B.; Gruhl, A.; Bartholomew, J.L. An Introduction to Myxozoan Evolution, Ecology and Development; Springer International Publishing: Berlin/Heidelberg, Germany, 2015; pp. 1–441. [Google Scholar]
  9. Eszterbauer, E.; Atkinson, S.; Diamant, A.; Morris, D.; El-Matbouli, M.; Hartikainen, H. Myxozoan Life Cycles: Practical approaches and insights. In Myxozoan Evolution, Ecology and Development; Okamura, B., Gruhl, A.B.J., Eds.; Springer: Cham, Switzerland, 2015; pp. 175–198. [Google Scholar]
  10. Holzer, A.S.; Bartošová-Sojková, P.; Born-Torrijos, A.; Alena, L.; Hartigan, A.; Fiala, I. The joint evolution of the Myxozoa and their alternate hosts: A cnidarian recipe for success and vast biodiversity. Mol. Ecol. 2018, 27, 1651–1666. [Google Scholar] [CrossRef]
  11. Wolf, K.; Markiw, M.E. Biology contravenes taxonomy in the Myxozoa: New discoveries show alternation of invertebrate and vertebrate hosts. Science 1984, 225, 1449–1452. [Google Scholar] [CrossRef]
  12. Kent, M.; Margolis, L.; Corliss, J. The demise of a class of protists: Taxonomic and nomenclatural revisions proposed for the protist phylum Myxozoa Grasse, 1970. Can. J. Zool. 1994, 72, 932–937. [Google Scholar] [CrossRef]
  13. Estensoro, I.; Jung-Schroers, V.; Álvarez-Pellitero, P.; Streinhagen, D.; Sitjà-Bobadilla, A. Effects of Enteromyxum leei (Myxozoa) infection on gilthead sea bream (Sparus aurata) (Teleostei) intestinal mucus: Glycoprotein profile and bacterial adhesion. Parasitol. Res. 2013, 112, 567–576. [Google Scholar] [CrossRef]
  14. Piazzon, M.C.; Galindo-Villegas, J.; Pereiro, P.; Estensoro, I.; Calduch-Giner, J.A.; Gomez-Casado, E.; Novoa, B.; Mulero, V.; Sitjà-Bobadilla, A.; Pérez-Sánchez, J. Differential modulation of IgT and IgM upon parasitic, bacterial, viral, and dietary challenges in a perciform Fish. Front. Immunol. 2016, 7, 236169. [Google Scholar] [CrossRef]
  15. Holzer, A.S.; Piazzon, M.C.; Barrett, D.; Bartholomew, J.L.; Sitjà-Bobadilla, A. To react or not to react: The dilemma of fish immune systems facing myxozoan infections. Front. Immunol. 2021, 12, 734238. [Google Scholar] [CrossRef] [PubMed]
  16. Tort, L.; Balasch, J.C.; Mackenzie, S. Fish immune system. A crossroads between innate and adaptive responses. Inmunología 2003, 22, 277–286. [Google Scholar]
  17. Zhu, L.; Nie, L.; Zhu, G.; Xiang, L.; Shao, J. Advances in research of fish immune-relevant genes: A comparative overview of innate and adaptive immunity in teleosts. Dev. Comp. Immunol. 2013, 39, 39–62. [Google Scholar] [CrossRef]
  18. Sitjà-Bobadilla, A.; Schmidt-Posthaus, H.; Wahli, T.; Holland, J.W.; Secombes, C.J. Fish immune responses to Myxozoa. In Myxozoan Evolution, Ecology and Development; Okamura, B., Gruhl, A., Bartholome, J.L., Eds.; Springer: Cham, Switzerland, 2015; pp. 253–280. ISBN 9783319147536. [Google Scholar]
  19. Bjørgen, H.; Koppang, E.O. Anatomy of teleost fish immune structures and organs. In Principles of Fish Immunology, from Host Cells and Molecules to Host Protection; Buchmann, K., Secombes, C.J., Eds.; Springer: Cham, Switzerland, 2021; pp. 1–30. [Google Scholar]
  20. Sitjà-Bobadilla, A.; Calduch-Giner, J.; Saera-Vila, A.; Palenzuela, O.; Álvarez-Pellitero, P.; Pérez-Sánchez, J. Chronic exposure to the parasite Enteromyxum leei (Myxozoa: Myxosporea) modulates the immune response and the expression of growth, redox and immune relevant genes in gilthead sea bream, Sparus aurata L. Fish Shellfish Immunol. 2008, 24, 610–619. [Google Scholar] [CrossRef] [PubMed]
  21. Sitjà-Bobadilla, A.; Redondo, M.J.; Bermúdez, R.; Palenzuela, O.; Ferreiro, I.; Riaza, A.; Quiroga, I.; Nieto, J.M.; Alvarez-Pellitero, P. Innate and adaptive immune responses of turbot, Scophthalmus maximus (L.), following experimental infection with Enteromyxum scophthalmi (Myxosporea: Myxozoa). Fish Shellfish Immunol. 2006, 21, 485–500. [Google Scholar] [CrossRef] [PubMed]
  22. Barrett, D.E.; Estensoro, I.; Sitj, A.; Bartholomew, J.L. Intestinal transcriptomic and histologic profiling reveals tissue repair mechanisms underlying resistance to the parasite Ceratonova shasta. Pathogens 2021, 10, 1179. [Google Scholar] [CrossRef] [PubMed]
  23. Cuesta, A.; Salinas, I.; Rodríguez, A.; Muñoz, P.; Sitjà-Bobadilla, A.; Álvarez-Pellitero, P.; Meseguer, J.; Esteban, M.Á. Cell-mediated cytotoxicity is the main innate immune mechanism involved in the cellular defence of gilthead seabream (Teleostei: Sparidae) against Enteromyxum leei (Myxozoa). Parasite Immunol. 2006, 28, 657–665. [Google Scholar] [CrossRef] [PubMed]
  24. Zou, J.; Secombes, C.J. The function of fish cytokines. Biology 2016, 5, 23. [Google Scholar] [CrossRef] [PubMed]
  25. Koehler, A.; Romans, P.; Desser, S.; Ringuette, M. Encapsulation of Myxobolus pendula (Myxosporidia) by epithelioid cells of its cyprinid host Semotilus atromaculatus. J. Parasitol 2004, 90, 1401–1405. [Google Scholar] [CrossRef]
  26. Balovet, G.; Baudin Laurencin, F. Granulomatous nodules in fish: An experimental assessment in rainbow trout, Salmo gairdneri Richardson, and turbot, Scophthalmus maximus (L.). J. Fish Dis. 1986, 9, 417–429. [Google Scholar] [CrossRef]
  27. Lövy, A.; Smirnov, M.; Brekhman, V.; Ofek, T.; Lotan, T. Morphological and molecular characterization of a novel myxosporean parasite Myxobolus bejeranoi n. sp. (Cnidaria: Myxosporea) from hybrid tilapia in Israel. Parasitol. Res. 2018, 117, 491–499. [Google Scholar] [CrossRef] [PubMed]
  28. Kallert, D.M.; Grabner, D.; Yokoyama, H.; El-Matbouli, M.; Eszterbauer, E. Transmission of Myxozoans to vertebrate hosts. In Myxozoan Evolution, Ecology and Development; Okamura, B., Gruhl, A., Bartholomew, J.L., Eds.; Springer: Cham, Switzerland, 2015; pp. 235–251. [Google Scholar]
  29. Maor-Landaw, K.; Avidor, I.; Rostowsky, N.; Salti, B.; Smirnov, M.; Ofek-Lalzar, M.; Levin, L.; Brekhman, V.; Lotan, T. The molecular mechanisms employed by the parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) from invasion through sporulation for successful proliferation in its fish host. Int. J. Mol. Sci. 2023, 24, 12824. [Google Scholar] [CrossRef] [PubMed]
  30. Miao, J.; Fan, Q.; Cui, L.; Li, J.; Li, J.; Cui, L. The malaria parasite Plasmodium falciparum histones: Organization, expression, and acetylation. Gene 2006, 369, 53–65. [Google Scholar] [CrossRef] [PubMed]
  31. Ramírez-Toloza, G.; Aguilar-Guzmán, L.; Valck, C.; Ferreira, V.P.; Ferreira, A. The Interactions of parasite calreticulin with initial complement components: Consequences in immunity and virulence. Front. Immunol. 2020, 11, 556854. [Google Scholar] [CrossRef]
  32. Alama-Bermejo, G.; Bartošová-Sojková, P.; Atkinson, S.D.; Holzer, A.S.; Bartholomew, J.L. Proteases as therapeutic targets against the parasitic cnidarian Ceratonova shasta: Characterization of molecules key to parasite virulence in Salmonid hosts. Front. Cell Infect. Microbiol. 2022, 11, 804864. [Google Scholar] [CrossRef] [PubMed]
  33. Maor-Landaw, K.; Avidor, I.; Salti, B.; Smirnov, M.; Brekhman, V.; Lotan, T. The myxozoan parasite Myxobolus bejeranoi (Cnidaria: Myxozoa) infection dynamics and host specificity in hybrid tilapia aquaculture. Parasitology 2023, 150, 524–530. [Google Scholar] [CrossRef] [PubMed]
  34. Schwanhüusser, B.; Busse, D.; Li, N.; Dittmar, G.; Schuchhardt, J.; Wolf, J.; Chen, W.; Selbach, M. Global quantification of mammalian gene expression control. Nature 2011, 473, 337–342. [Google Scholar] [CrossRef] [PubMed]
  35. Ramulu, S.; Kale, A.D.; Hallikerimath, S.; Kotrashetti, V. Comparing modified papanicolaou stain with ayoub-shklar and haematoxylin-eosin stain for demonstration of keratin in paraffin embedded tissue sections. J. Oral Maxillofac. Pathol. 2013, 17, 23–30. [Google Scholar] [CrossRef] [PubMed]
  36. Anthwal, N.; Gupta, S.; Singh, R.P.; Gupta, H. Comparison of ayoub shklar stain, Dane Herman, modified pap and routine hematoxylin and eosin stain for the keratin identification. Int. Arch. Integr. Med. 2020, 7, 1–8. [Google Scholar]
  37. Kvedar, J.C.; Manabe, M.; Phillips, S.B.; Ross, B.S.; Baden, H.P. Characterization of sciellin, a precursor to the cornified envelope of human keratinocytes. Differentiation 1992, 49, 195–204. [Google Scholar] [CrossRef]
  38. Sundararajan, V.; Pang, Q.Y.; Choolani, M.; Huang, R.Y.J. Spotlight on the granules (grainyhead-like proteins)—From an evolutionary conserved controller of epithelial trait to pioneering the chromatin landscape. Front. Mol. Biosci. 2020, 7, 213. [Google Scholar] [CrossRef] [PubMed]
  39. Ceballos-Francisco, D.; Cordero, H.; Guardiola, F.A.; Cuesta, A.; Esteban, M.Á. Healing and mucosal immunity in the skin of experimentally wounded gilthead seabream (Sparus aurata L.). Fish Shellfish Immunol. 2017, 71, 210–219. [Google Scholar] [CrossRef]
  40. Wang, S.; Samakovlis, C. Grainy Head and its target genes in epithelial morphogenesis and wound healing. Curr. Top. Dev. Biol. 2012, 98, 35–63. [Google Scholar] [PubMed]
  41. Kalinin, A.E.; Kajava, A.V.; Steinert, P.M. Epithelial barrier function: Assembly and structural features of the cornified cell envelope. BioEssays 2002, 24, 789–800. [Google Scholar] [CrossRef] [PubMed]
  42. Martyn, A.A.; Hong, H.; Ringuette, M.J.; Desser, S.S. Changes in host and parasite-derived cellular and extracellular matrix components in developing cysts of of Myxobolus pendula (Myxozoa). J. Eukaryot. Microbiol. 2002, 19, 175–182. [Google Scholar] [CrossRef] [PubMed]
  43. Polinas, M.; Padr, F.; Merella, P.; Prearo, M.; Sanna, M.A.; Marino, F.; Burrai, G.P.; Antuofermo, E. Stages of granulomatous response against histozoic Metazoan parasites in Mullets (Osteichthyes: Mugilidae). Animals 2021, 11, 1501. [Google Scholar] [CrossRef] [PubMed]
  44. Qiu, R.; Sun, B.G.; Li, J.; Liu, X.; Sun, L. Identification and characterization of a cell surface scavenger receptor cysteine-rich protein of Sciaenops ocellatus: Bacterial interaction and its dependence on the conserved structural features of the SRCR domain. Fish Shellfish Immunol. 2013, 34, 810–818. [Google Scholar] [CrossRef] [PubMed]
  45. He, J.; Liu, H.; Yang, J.; Dong, X.; Wu, C. Abundant members of Scavenger receptors family and their identification, characterization and expression against Vibrio alginolyticus infection in juvenile Larimichthys crocea. Fish Shellfish Immunol. 2016, 50, 297–309. [Google Scholar] [CrossRef] [PubMed]
  46. Law, S.K.A.; Micklem, K.J.; Shaw, J.M.; Zhang, X.-P.; Dong, Y.; Willis, A.C.; Mason, D.Y. A new macrophage differentiation antigen which is a member of the scavenger receptor superfamily. Eur. J. Immunol. 1993, 23, 2320–2325. [Google Scholar] [CrossRef]
  47. Onofre, G.; Kolácková, M.; Jankovicová, K.; Krejsek, J. Scavenger receptor CD163 and its biological functions. Acta Medica Cordoba 2009, 52, 57–61. [Google Scholar] [CrossRef]
  48. Martínez, V.G.; Moestrup, S.K.; Holmskov, U.; Mollenhauer, J.; Lozano, F. The conserved scavenger receptor cysteine-rich super family in therapy and diagnosis. Pharmacol. Rev. 2011, 63, 967–1000. [Google Scholar] [CrossRef]
  49. Orf, K.; Cunnington, A.J. Infection-related hemolysis and susceptibility to Gram-negative bacterial co-infection. Front. Microbiol. 2015, 6, 666. [Google Scholar] [CrossRef]
  50. Korytář, T.; Wiegertjes, G.F.; Zusková, E.; Tomanová, A.; Lisnerová, M.; Patra, S.; Sieranski, V.; Šíma, R.; Born-Torrijos, A.; Wentzel, A.S.; et al. The kinetics of cellular and humoral immune responses of common carp to presporogonic development of the myxozoan Sphaerospora molnari. Parasites Vectors 2019, 12, 208. [Google Scholar] [CrossRef]
  51. Korytář, T.; Chan, J.T.H.; Vancová, M.; Holzer, A.S. Blood feast: Exploring the erythrocyte-feeding behaviour of the myxozoan Sphaerospora molnari. Parasite Immunol. 2020, 42, 1–7. [Google Scholar] [CrossRef] [PubMed]
  52. Hartigan, A.; Kosakyan, A.; Pecková, H.; Eszterbauer, E.; Holzer, A.S. Transcriptome of Sphaerospora molnari (Cnidaria, Myxosporea) blood stages provides proteolytic arsenal as potential therapeutic targets against sphaerosporosis in common carp. BMC Genom. 2020, 21, 404. [Google Scholar] [CrossRef] [PubMed]
  53. Yang, Y.; Xiong, J.; Zhou, Z.; Huo, F.; Miao, W.; Ran, C.; Liu, Y.; Zhang, J.; Feng, J.; Wang, M.; et al. The genome of the myxosporean Thelohanellus kitauei shows adaptation to nutrient acquisition within its fish host. Genome Biol. Evol. 2014, 6, 3182–3198. [Google Scholar] [CrossRef]
  54. Ahmad, F.; Debes, P.V.; Pukk, L.; Kahar, S.; Hartikainen, H.; Gross, R.; Vasemägi, A. Know your enemy—Transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids. Parasitology 2021, 148, 726–739. [Google Scholar] [CrossRef] [PubMed]
  55. Faber, M.; Shaw, S.; Yoon, S.; de Paiva Alves, E.; Wang, B.; Qi, Z.; Okamura, B.; Hartikainen, H.; Secombes, C.J.; Holland, J.W. Comparative transcriptomics and host-specific parasite gene expression profiles inform on drivers of proliferative kidney disease. Sci. Rep. 2021, 11, 2149. [Google Scholar] [CrossRef]
  56. Davey, G.C.; Calduch-Giner, J.A.; Houeix, B.; Talbot, A.; Sitjà-Bobadilla, A.; Prunet, P.; Pérez-Sánchez, J.; Cairns, M.T. Molecular profiling of the gilthead sea bream (Sparus aurata L.) response to chronic exposure to the myxosporean parasite Enteromyxum leei. Mol. Immunol. 2011, 48, 2102–2112. [Google Scholar] [CrossRef]
  57. Noga, E.J.; Dykstra, M.J.; Wright, J.F. Chronic inflammatory cells with epithelial cell characteristics in teleost fishes. Vet. Pathol. 1989, 26, 429–437. [Google Scholar] [CrossRef]
  58. Sayyaf Dezfuli, B.; Fernandes, C.E.; Galindo, G.M.; Castaldelli, G.; Manera, M.; DePasquale, J.A.; Lorenzoni, M.; Bertin, S.; Giari, L. Nematode infection in liver of the fish Gymnotus inaequilabiatus (Gymnotiformes: Gymnotidae) from the Pantanal Region in Brazil: Pathobiology and inflammatory response. Parasites Vectors 2016, 9, 473. [Google Scholar] [CrossRef]
  59. Manrique, W.G.; Claudianoda, G.S.; DeCastro, M.P.; Petrillo, T.R.; Figueiredo, M.A.P.; De Andrade Belo, M.A.; Berdeal, M.I.Q.; De Moraes, J.E.R.; De Moraes, F.R. Expression of cellular components in granulomatous inflammatory response in Piaractus mesopotamicus model. PLoS ONE 2015, 10, e0121625. [Google Scholar] [CrossRef] [PubMed]
  60. Skorobrechova, E.M.; Nikishin, V.P. Structure of capsule surrounding acanthocephalans Corynosoma strumosum in paratenic hosts of three species. Parasitol. Res. 2011, 108, 467–475. [Google Scholar] [CrossRef] [PubMed]
  61. Bragulla, H.H.; Homberger, D.G. Structure and functions of keratin proteins in simple, stratified, keratinized and cornified epithelia. J. Anat. 2009, 214, 516–559. [Google Scholar] [CrossRef]
  62. Rübsam, M.; Broussard, J.A.; Wickström, S.A.; Nekrasova, O.; Green, K.J.; Niessen, C.M. Adherens junctions and desmosomes coordinate mechanics and signaling to orchestrate tissue morphogenesis and function: An evolutionary perspective. Cold Spring Harb. Perspect. Biol. 2018, 10, a029207. [Google Scholar] [CrossRef]
  63. Niessen, C.M. Tight junctions/adherens junctions: Basic structure and function. J. Investig. Dermatol. 2007, 127, 2525–2532. [Google Scholar] [CrossRef]
  64. Champliaud, M.; Burgeson, R.E.; Jin, W.; Baden, H.P.; Olson, P.F. cDNA cloning and characterization of Sciellin, a LIM domain protein of the keratinocyte cornified envelope. J. Biol. Chem. 1998, 273, 31547–31554. [Google Scholar] [CrossRef]
  65. Champliaud, M.F.; Baden, H.P.; Koch, M.; Jin, W.; Burgeson, R.E.; Viel, A. Gene characterization of sciellin (SCEL) and protein localization in vertebrate epithelia displaying barrier properties. Genomics 2000, 70, 264–268. [Google Scholar] [CrossRef] [PubMed]
  66. Snider, N.T.; Weerasinghe, S.V.W.; In, J.A.; Herrmann, H.; Omary, M.B. Keratin hypersumoylation alters filament dynamics and is a marker for human liver disease and keratin mutation. J. Biol. Chem. 2011, 286, 2273–2284. [Google Scholar] [CrossRef]
  67. Pichler, A.; Gast, A.; Seeler, J.S.; Dejean, A.; Melchior, F. The nucleoporin RanBP2 has SUMO1 E3 ligase activity. Cell 2002, 108, 109–120. [Google Scholar] [CrossRef]
  68. Alonso, A.; Greenlee, M.; Matts, J.; Kline, J.; Davis, K.J.; Miller, R.K. Emerging roles of sumoylation in the regulation of actin, microtubules, intermediate filaments, and septins. Cytoskeleton 2015, 339, 305–339. [Google Scholar] [CrossRef]
  69. Kraemer, A.M.; Saraiva, L.R.; Korsching, S.I. Structural and functional diversification in the teleost S100 family of calcium-binding proteins. BMC Evol. Biol. 2008, 8, 48. [Google Scholar] [CrossRef]
  70. Hsiao, C.D.; Ekker, M.; Tsai, H.J. Skin-specific expression of ictacalcin, a homolog of the S100 genes, during zebrafish embryogenesis. Dev. Dyn. 2003, 228, 745–750. [Google Scholar] [CrossRef]
  71. Pao, H.Y.; Wu, C.Y.; Huang, C.H.; Wen, C.M. Development, characterization and virus susceptibility of a continuous cell line from the caudal fin of marbled eel (Anguilla marmorata). J. Fish Dis. 2018, 41, 1331–1338. [Google Scholar] [CrossRef]
  72. Marcos-López, M.; Rodger, H.D.; O’Connor, I.; Braceland, M.; Burchmore, R.J.S.; Eckersall, P.D.; MacCarthy, E. A proteomic approach to assess the host response in gills of farmed Atlantic salmon Salmo salar L. affected by amoebic gill disease. Aquaculture 2017, 470, 1–10. [Google Scholar] [CrossRef]
  73. Easy, R.H.; Ross, N.W. Changes in Atlantic salmon (Salmo salar) epidermal mucus protein composition profiles following infection with sea lice (Lepeophtheirus salmonis). Comp. Biochem. Physiol. Part D Genom. Proteom. 2009, 4, 159–167. [Google Scholar] [CrossRef] [PubMed]
  74. Pridgeon, J.W.; Yeh, H.Y.; Shoemaker, C.A.; Mu, X.; Klesius, P.H. Global gene expression in channel catfish after vaccination with an attenuated Edwardsiella ictaluri. Fish Shellfish Immunol. 2012, 32, 524–533. [Google Scholar] [CrossRef]
  75. Blaufuss, P.C.; Gaylord, T.G.; Sealey, W.M.; Powell, M.S. Effects of high-soy diet on S100 gene expression in liver and intestine of rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 2019, 86, 764–771. [Google Scholar] [CrossRef] [PubMed]
  76. Kuberappa, P.H.; Bagalad, B.S.; Ananthaneni, A.; Kiresur, M.A.; Srinivas, G.V. Certainty of S100 from physiology to pathology. J. Clin. Diagn. Res. 2016, 10, ZE10–ZE15. [Google Scholar] [CrossRef]
  77. Pelc, P.; Vanmuylder, N.; Lefranc, F.; Heizmann, C.W.; Hassid, S.; Salmon, I.; Kiss, R.; Louryan, S.; Decaestecker, C. Differential expression of S100 calcium-binding proteins in epidermoid cysts, branchial cysts, craniopharyngiomas and cholesteatomas. Histopathology 2003, 42, 387–394. [Google Scholar] [CrossRef]
  78. Chang, C.-H.; Wu, Y.-C.; Wu, Y.-H.; Sun, A.; Kuo, Y.-S.; Chiang, C.-P. S100 protein-positive Langerhans cells in 80 dentigerous cysts. J. Dent. Sci. 2017, 12, 405–412. [Google Scholar] [CrossRef] [PubMed]
  79. Perez-Riverol, Y.; Bai, J.; Bandla, C.; García-Seisdedos, D.; Hewapathirana, S.; Kamatchinathan, S.; Kundu, D.J.; Prakash, A.; Frericks-Zipper, A.; Eisenacher, M.; et al. The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 2022, 50, D543–D552. [Google Scholar] [CrossRef] [PubMed]
  80. Cox, J.; Hein, M.Y.; Luber, C.A.; Paron, I.; Nagaraj, N.; Mann, M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol. Cell Proteom. 2014, 13, 2513–2526. [Google Scholar] [CrossRef] [PubMed]
  81. Tyanova, S.; Cox, J. Perseus: A Bioinformatics Platform for Integrative Analysis of Proteomics Data in Cancer Research. In Cancer Systems Biology: Methods and Protocols; von Stechow, L., Ed.; Humana Press: New York, NY, USA, 2018; Volume 1711, pp. 133–148. ISBN 9781493974931. [Google Scholar]
  82. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote) omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar] [CrossRef] [PubMed]
  83. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein—Protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed]
  84. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  85. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
Ijms 25 05683 g001
Figure 1. Venn diagram of identified protein groups in healthy control gills, infected gills, and cysts’ proteomes.
Figure 1. Venn diagram of identified protein groups in healthy control gills, infected gills, and cysts’ proteomes.
Ijms 25 05683 g001
Ijms 25 05683 g002
Figure 2. Enrichment analysis of parasite-infected fish gills versus control samples. (A,B) Up-regulated processes are presented in an enrichment score graph and protein interaction networks (STRING database). (C,D) Down-regulated processes. Node colors correspond with label in the graph. Details can be found in Supplementary Table S2.
Figure 2. Enrichment analysis of parasite-infected fish gills versus control samples. (A,B) Up-regulated processes are presented in an enrichment score graph and protein interaction networks (STRING database). (C,D) Down-regulated processes. Node colors correspond with label in the graph. Details can be found in Supplementary Table S2.
Ijms 25 05683 g002
Ijms 25 05683 g003
Figure 3. Log10 LFQ intensity values of the highest (top row) (SRCR, pentraxin 3, superoxide dismutase 2, and glutathione S-transferase-) and the lowest (bottom row) (hemoglobin subunit beta1, hemoglobin subunit alpha, collagen alpha-1, and col1α2), ranking proteins in infected gills (Inf) versus control (Cnt) analysis.
Figure 3. Log10 LFQ intensity values of the highest (top row) (SRCR, pentraxin 3, superoxide dismutase 2, and glutathione S-transferase-) and the lowest (bottom row) (hemoglobin subunit beta1, hemoglobin subunit alpha, collagen alpha-1, and col1α2), ranking proteins in infected gills (Inf) versus control (Cnt) analysis.
Ijms 25 05683 g003
Ijms 25 05683 g004
Figure 4. The proteomic profile of cysts. (A) Developed cysts (mark with arrows) are localized at the base of the gill filament. Bar 250 μc. (B,C) Enrichment analysis of cysts versus control samples. Enrichment scores, calculated from adjusted p-values, are presented for up-regulated (B) and down-regulated processes (C). Details can be found in Supplementary Table S3.
Figure 4. The proteomic profile of cysts. (A) Developed cysts (mark with arrows) are localized at the base of the gill filament. Bar 250 μc. (B,C) Enrichment analysis of cysts versus control samples. Enrichment scores, calculated from adjusted p-values, are presented for up-regulated (B) and down-regulated processes (C). Details can be found in Supplementary Table S3.
Ijms 25 05683 g004
Ijms 25 05683 g005
Figure 5. Keratin structure of the cyst’s wall and its protein network. Sections of cysts stained with (A) Haematoxylin and eosin, (B) Dane-Herman, (C) Ayoub-Shklar, showing a keratin-specific staining (mark with arrow); Bar 50 μc. (D) Protein interaction network (STRING database) of keratin and sciellin ‘first-neighbors’ that are expressed in the cysts. Intermediate filament proteins are labelled by a square node shape and node fill color is indicative of protein expression pattern.
Figure 5. Keratin structure of the cyst’s wall and its protein network. Sections of cysts stained with (A) Haematoxylin and eosin, (B) Dane-Herman, (C) Ayoub-Shklar, showing a keratin-specific staining (mark with arrow); Bar 50 μc. (D) Protein interaction network (STRING database) of keratin and sciellin ‘first-neighbors’ that are expressed in the cysts. Intermediate filament proteins are labelled by a square node shape and node fill color is indicative of protein expression pattern.
Ijms 25 05683 g005
Table 1. Endopeptidase inhibitors in cyst wall. Fold change in significant iBAQ up-regulated proteins are shown for cyst versus control. Proteins that are expressed only in cysts and absent in control samples are indicated as ‘cyst only’.
Table 1. Endopeptidase inhibitors in cyst wall. Fold change in significant iBAQ up-regulated proteins are shown for cyst versus control. Proteins that are expressed only in cysts and absent in control samples are indicated as ‘cyst only’.
Protein IDProtein NameFold Change
ENSONIP00000013827.2fetuin-B cyst only
ENSONIP00000042223.1alpha-2-HS-glycoprotein 2 cyst only
ENSONIP00000023717.1kunitz-type protease inhibitor 1 cyst only
ENSONIP00000074191.1alpha-macroglobulin receptor-binding domaincyst only
ENSONIP00000067950.1alpha-2-macroglobulin-like cyst only
ENSONIP00000005221.1alpha-2-antiplasmin cyst only
ENSONIP00000018223.2inter-alpha-trypsin inhibitor heavy chain 3 cyst only
ENSONIP00000022637.2leukocyte elastase inhibitor 7.71
ENSONIP00000025380.1serpin family A member 10 4.70
ENSONIP00000056451.1legumain 4.15
ENSONIP00000046993.1cystatin-B 2.51
ENSONIP00000043069.1cystatin fetuin-A-type domain2.15
ENSONIP00000001504.1alpha-1-antitrypsin homolog 1.45
ENSONIP00000009801.2serpin family C member 1 1.43
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maor-Landaw, K.; Smirnov, M.; Lotan, T. The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function. Int. J. Mol. Sci. 2024, 25, 5683. https://doi.org/10.3390/ijms25115683

AMA Style

Maor-Landaw K, Smirnov M, Lotan T. The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function. International Journal of Molecular Sciences. 2024; 25(11):5683. https://doi.org/10.3390/ijms25115683

Chicago/Turabian Style

Maor-Landaw, Keren, Margarita Smirnov, and Tamar Lotan. 2024. "The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function" International Journal of Molecular Sciences 25, no. 11: 5683. https://doi.org/10.3390/ijms25115683

APA Style

Maor-Landaw, K., Smirnov, M., & Lotan, T. (2024). The Tilapia Cyst Tissue Enclosing the Proliferating Myxobolus bejeranoi Parasite Exhibits Cornified Structure and Immune Barrier Function. International Journal of Molecular Sciences, 25(11), 5683. https://doi.org/10.3390/ijms25115683

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop