INTRODUCTION: Parasite-mediated selection is considered one of the potential mechanisms contributing to the coexistence of asexual-sexual complexes. Gibel carp (Carassius gibelio), an invasive fish species in Europe, often forms populations composed of gynogenetic and sexual specimens. METHODS: The experimental infection was induced in gynogenetic and sexual gibel carp using eye-fluke Diplostomum pseudospathaceum (Trematoda), and the transcriptome profile of the spleen as a major immune organ in fish was analyzed to reveal the differentially expressed immunity-associated genes related to D. pseudospathaceum infection differing between gynogenetic and sexual gibel carp. RESULTS: High parasite infection was found in gynogenetic fish when compared to genetically diverse sexuals. Although metacercariae of D. pseudospathaceum are situated in an immune-privileged organ, our results show that eye trematodes may induce a host immune response. We found differential gene expression induced by eye-fluke infection, with various impacts on gynogenetic and sexual hosts, documenting for the majority of DEGs upregulation in sexuals, and downregulation in asexuals. Differences in gene regulation between gynogenetic and sexual gibel carp were evidenced in many immunity-associated genes. GO analyses revealed the importance of genes assigned to the GO terms: immune function, the Notch signaling pathway, MAP kinase tyrosine/threonine/phosphatase activity, and chemokine receptor activity. KEGG analyses revealed the importance of the genes involved in 12 immunity-associated pathways - specifically, FoxO signaling, adipocytokine signaling, TGF-beta signaling, apoptosis, Notch signaling, C-type lectin receptor signaling, efferocytosis, intestinal immune network for IgA production, insulin signaling, virion - human immunodeficiency virus, Toll-like receptor signaling, and phosphatidylinositol signaling system. DISCUSSION: Our study indicates the limited potential of asexual fish to cope with higher parasite infection (likely a loss of capacity to induce an effective immune response) and highlights the important role of molecular mechanisms associated with immunity for the coexistence of gynogenetic and sexual gibel carp, potentially contributing to its invasiveness.

Central European Institute of Technology Masaryk University Brno Czechia

Department of Botany and Zoology Faculty of Science Masaryk University Brno Czechia

Institute of Vertebrate Biology of the Czech Academy of Sciences Brno Czechia

Laboratory of Non Mendelian Evolution Institute of Animal Physiology and Genetics of the Czech Academy of Science Liběchov Czechia

University of South Bohemia in České Budějovice Faculty of Fisheries and Protection of Waters South Bohemian Research Center of Aquaculture and Biodiversity of Hydrocenoses Vodňany Czechia

See more in PubMed

Lusková V, Lusk S, Halačka K, Vetešník L. DOI

Fuad MMH, Vetešník L, Šimková A. Is gynogenetic reproduction in gibel carp ( DOI

Tobler M, Schlupp I. Parasites in sexual and asexual mollies ( PubMed DOI PMC

Neaves WB, Baumann P. Unisexual reproduction among vertebrates. Trends Genet. (2011) 27:81–8. doi:  10.1016/j.tig.2010.12.002 PubMed DOI

Lusková V, Halačka K, Vetešník L, Lusk S. Changes of ploidy and sexuality status of “

Przybył A, Przybylski M, Spóz A, Juchno D, Szabelska A, Kowalewska K, et al. Sex, size and ploidy ratios of DOI

Paschos I, Nathanailides C, Tsoumani M, Perdikaris C, Gouva E, Leonardos I. Intra and inter-specific mating options for gynogenetic reproduction of

Barbuti R, Mautner S, Carnevale G, Milazzo P, Rama A, Sturmbauer C. Population dynamics with a mixed type of sexual and asexual reproduction in a fluctuating environment. BMC Evol Biol. (2012) 12:49. doi:  10.1186/1471-2148-12-49 PubMed DOI PMC

Mee JA, Rowe L. A comparison of parasite loads on asexual and sexual DOI

Mee JA, Taylor EB. The cybrid invasion: widespread postglacial dispersal by DOI

Lamatsch DK, Schmid M, Schartl M. A somatic mosaic of the gynogenetic Amazon molly. J Fish Biol. (2002) 60:1417–22. doi:  10.1111/j.1095-8649.2002.tb02436.x DOI

Lamatsch DK, Lampert KP, Fischer P, Geiger M, Schlupp I, Schartl M. Diploid Amazon mollies ( DOI

Tobler M, Schlupp I. Differential susceptibility to food stress in neonates of sexual and asexual mollies ( DOI

Schlupp I. Mate choice and the amazon molly: how sexuality and unisexuality can coexist. J Hered. (2010) 101:S55–61. doi:  10.1093/jhered/esq015 PubMed DOI

Janko K, Eisner J. Sperm-dependent parthenogens delay the spatial expansion of their sexual hosts. J Theor Biol. (2009) 261:431–40. doi:  10.1016/j.jtbi.2009.08.012 PubMed DOI

Choleva L, Apostolou A, Ráb P, Janko K. Making it on their own: sperm-dependent hybrid fishes ( PubMed DOI PMC

Hakoyama H, Nishimura T, Matsubara N, Iguchi K. Difference in parasite load and nonspecific immune reaction between sexual and gynogenetic forms of DOI

Hakoyama H, Iwasa Y. Coexistence of a sexual and an unisexual form stabilized by parasites. J Theor Biol. (2004) 226:186–94. doi:  10.1016/j.jtbi.2003.08.012 PubMed DOI

Šimková A, Košař M, Vetešník L, Vyskočilová M. MHC genes and parasitism in PubMed DOI PMC

Šimková A, Hyršl P, Halačka K, Vetešník L. Physiological and condition-related traits in the gynogenetic-sexual PubMed DOI PMC

Vetešník L, Lusk S, Halačka K, Spurný P. Morphometric characteristics and growth of

Vetešník L, Halačka K, Šimková A. The effect of ploidy and temporal changes in the biochemical profile of gibel carp ( PubMed DOI

Xie J, Wen JJ, Chen B, Gui JF. Differential gene expression in fully-grown oocytes between gynogenetic and gonochoristic crucian carps. Gene. (2001) 271:109–16. doi:  10.1016/s0378-1119(01)00491-7 PubMed DOI

Pakosta T, Vetešník L, Šimková A. A long temporal study of parasitism in asexualsexual populations of PubMed DOI PMC

Hamilton WD. Sex versus non-sex versus parasite. Oikos. (1980) 35:282–90. doi:  10.2307/3544435 DOI

Hamilton WD, Axelrod R, Tanese R. Sexual reproduction as an adaptation to resist parasites (a review). Proc Natl Acad Sci U S A. (1990) 87:3566–73. doi:  10.1073/pnas.87.9.3566 PubMed DOI PMC

Lively CM, Craddock C, Vrijenhoek RC. Red Queen hypothesis supported by parasitism in sexual and clonal fish. Nature. (1990) 344:864–67. doi:  10.1038/344864a0 DOI

Weeks SC. A reevaluation of the Red Queen model for the maintenance of sex in a clonal-sexual fish complex (Poeciliidae: DOI

Scharsack JP, Kalbe M. Differences in susceptibility and immune responses of three-spined sticklebacks ( PubMed DOI PMC

Niewiadomska K. Verification of the life-cycles of DOI

Frey RA, Barrett LM, Parkin L, Blakeley B, Alund M, Byford G, et al. Eye flukes ( PubMed DOI

Hakalahti T, Karvonen A, Valtonen ET. Climate warming and disease risks in temperate regions - PubMed DOI

Palmieri JR, Cali A, Heckmann RA. Experimental biological-control of eye fluke, PubMed DOI

Michálková V, Ondračková M. Experimental evidence for parasite-induced over-winter mortality in juvenile. Rhodeus amarus. J Fish Biol. (2014) 84:1377–88. doi:  10.1111/jfb.12363 PubMed DOI

Crowden AE, Broom DM. Effects of the eyefluke, DOI

Mikheev VN, Pasternak AF, Taskinen J, Valtonen ET. Parasite-induced aggression and impaired contest ability in a fish host. Parasit Vectors. (2010) 3:17. doi:  10.1186/1756-3305-3-17 PubMed DOI PMC

Seppälä O, Karvonen A, Valtonen ET. Parasite-induced change in host behaviour and susceptibility to predation in an eye fluke - fish interaction. Anim Behav. (2004) 68:257–63. doi:  10.1016/j.anbehav.2003.10.021 DOI

Seppälä O, Karvonen A, Valtonen ET. Manipulation of fish host by eye flukes in relation to cataract formation and parasite infectivity. Anim Behav. (2005) 70:889–94. doi:  10.1016/j.anbehav.2005.01.020 DOI

Mikheev VN, Pasternak AF. Structure of aggressive behavior in underyearlings of the rainbow trout DOI

Nezhybová V, Reichard M, Methling C, Ondračková M. Limited impacts of chronic eye fluke infection on the reproductive success of a fish host. Biol J Linn Soc. (2020) 129:334–46. doi:  10.1093/biolinnean/blz189 DOI

Jakobsson S, Brick O, Kullberg C. Escalated fighting behavior incurs increased predation risk. Anim Behav. (1995) 49:235–39. doi:  10.1016/0003-3472(95)80172-3 DOI

McDougall PT, Kramer DL. Short-term behavioral consequences of territory relocation in a Caribbean damselfish, DOI

Buchmann K. Antiparasitic immune mechanisms in teleost fish: a two-edged sword? Bull Eur Ass Fish Pathol. (2000) 20:48–59.

Rauch G, Kalbe M, Reusch TBH. One day is enough: rapid and specific host–parasite interactions in a stickleback-trematode system. Biol Lett. (2006) 2:382–84. doi:  10.1098/rsbl.2006.0462 PubMed DOI PMC

Haase D, Rieger JK, Witten A, Stoll M, Bornberg-Bauer E, Kalbe M, et al. Specific gene expression responses to parasite genotypes reveal redundancy of innate immunity in vertebrates. PLoS One. (2014) 9:e108001. doi:  10.1371/journal.pone.0108001 PubMed DOI PMC

Whyte SK, Chappell LH, Secombes CJ. Cytotoxic reactions of rainbow trout, DOI

Woo PTK. Immunological responses of fish to parasitic organisms. Ann Rev Fish Dis. (1992) 2:339–66. doi:  10.1016/0169-4758(87)90178-5 DOI

Haase D, Rieger JK, Witten A, Stoll M, Bornberg-Bauer E, Kalbe M, et al. Immunity comes first: The effect of parasite genotypes on adaptive immunity and immunization in three-spined sticklebacks. Dev Comp Immunol. (2016) 54:137e144. doi:  10.1016/j.dci.2015.09.008 PubMed DOI

Karvonen A, Hudson PJ, Seppälä O, Valtonen ET. Transmission dynamics of a trematode parasite: exposure, acquired resistance and parasite aggregation. Parasitol Res. (2004) 92:183–88. doi:  10.1007/s00436-003-1035-y PubMed DOI

Georgieva S, Soldánová M, Pérez-del-Olmo A, Dangel RD, Sitko J, Sures B, et al. Molecular prospecting for European PubMed DOI

Bush AO, Lafferty KD, Lotz JM, Shostaket AW. Parasitology meets ecology on its own terms: Margolis et al. revisited. J Parasitol. (1997) 83:575–83. doi:  10.2307/3284227 PubMed DOI

Andrews S. FastQC: A Quality Control Tool for High Throughput Sequence Data. (2010). Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

Bushnell B. BBMap: A fast, accurate, splice-aware aligner. Lawrence Berkeley National Laboratory. LBNL Report #: LBNL-7065E. (2014). Available at: https://escholarship.org/uc/item/1h3515gn

Dobin A, Carrie A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. (2013) 29:15–21. doi:  10.1093/bioinformatics/bts635 PubMed DOI PMC

Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. (2014) 30:923–30. doi:  10.1093/bioinformatics/btt656 PubMed DOI

R Core Team . A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; (2023). Available at: https://www.R-project.org/.

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genom Biol. (2014) 15:550. doi:  10.1186/s13059-014-0550-8 PubMed DOI PMC

Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. (1995) 57:289–300. doi:  10.1111/j.2517-6161.1995.tb02031.x DOI

Korotkevich G, Sukhov V, Budin N, Shpak B, Artyomov MN, Sergushichev A. Fast gene set enrichment analysis. bioRxiv. (2021). doi:  10.1101/060012 DOI

Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. (2011) 6:7. doi:  10.1371/journal.pone.0021800 PubMed DOI PMC

Durinck S, Spellman P, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. (2009) 4:1184–91. doi:  10.1038/nprot.2009.97 PubMed DOI PMC

Yu G, Wang L, Han Y, He Q. ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi:  10.1089/omi.2011.0118 PubMed DOI PMC

Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. (2001) 29:e45. doi:  10.1093/nar/29.9.e45 PubMed DOI PMC

Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. (2002) 3:1–12. doi:  10.1186/gb-2002-3-7-research0034 PubMed DOI PMC

Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real - time quantitative PCR and the 2–ΔΔCT method. Methods. (2001) 25:402–8. doi:  10.1006/meth.2001.1262 PubMed DOI

Kalbe M, Kurtz J. Local differences in immunocompetence reflect resistance of sticklebacks against the eye fluke PubMed DOI

Wegner KM, Kalbe M, Reusch TBH. Innate versus adaptive immunity in sticklebacks: evidence for trade-offs from a selection experiment. Evol Ecol. (2007) 21:473–83. doi:  10.1007/s10682-006-9129-5 DOI

Noreikiene K, Ozerov M, Ahmad F, Kõiv T, Kahar S, Gross R, et al. Humic- acid-driven escape from eye parasites revealed by RNA-seq and target-specific metabarcoding. Parasit Vectors. (2020) 13:433. doi:  10.1186/s13071-020-04306-9 PubMed DOI PMC

Carbon S, Ireland A, Mungall CJ, Shu SQ, Marshall B, Lewis S. AmiGO: online access to ontology and annotation data. Bioinformatics. (2009) 25:288–9. doi:  10.1093/bioinformatics/btn615 PubMed DOI PMC

Neagu M, Constantin C. Signal transduction in immune cells and protein kinases. Adv Exp Med Biol. (2021) 1275:133–49. doi:  10.1007/978-3-030-49844-3_5 PubMed DOI

He L, Zhao L, Li Q, Huang L, Qin Y, Zhuang Z, et al. PubMed DOI

Gong D, Xu L, Li W, Shang R, Chen J, Hu F, et al. Comparative analysis of liver transcriptomes associated with hypoxia tolerance in the gynogenetic blunt snout bream. Aquaculture. (2020) 523:735163. doi:  10.1016/j.aquaculture.2020.735163 DOI

Sauerwein H, Bendixen E, Restelli L, Ceciliani F. The adipose tissue in farm animals: a proteomic approach. Curr Protein Pept Sci. (2014) 15:146–55. doi:  10.2174/1389203715666140221123105 PubMed DOI

Sun J, Bian CC, Ji SH, Luo XL, Ji H. Greater potency of adipocytes compared with preadipocytes under lipopolysaccharide exposure in grass carp PubMed DOI

Pignatelli J, Castro R, Granja AG, Abós B, González L, Jensen LB, et al. Immunological characterization of the teleost adipose tissue and its modulation in response to viral infection and fat-content in the diet. PLoS One. (2014) 9:e110920. doi:  10.1371/journal.pone.0110920 PubMed DOI PMC

Wang YD, Rajanbabu V, Chen JY. Transcriptome analysis of medaka following epinecidin-1 and TH1–5 treatment of NNV infection. Fish Shellfish Immunol. (2015) 42:121–31. doi:  10.1016/j.fsi.2014.10.040 PubMed DOI

Chen J, Cai B, Tian C, Jiang D, Shi H, Huang Y, et al. RNA sequencing (RNA-Seq) analysis reveals liver lipid metabolism divergent adaptive response to low- and high-salinity stress in spotted scat ( PubMed DOI PMC

Jiang JL, Xu J, Ye L, Sun ML, Jiang ZQ, Mao MG. Identification of differentially expressed genes in gills of tiger puffer ( PubMed DOI

Li C, Wang L. Molecular characterization, expression and functional analysis of TGFβ1-b in crucian carp ( PubMed DOI

Gahr SA, Weber GM, Rexroad CE. Identification and expression of Smads associated with TGF-b/activin/nodal signaling pathways in the rainbow trout ( PubMed DOI

Zhou S, Liu Y, Dong J, Yang Q, Xu N, Yang Y, et al. Transcriptome analysis of goldfish ( PubMed DOI

O’Connor MN, Salles II, Cvejic A, Watkins NA, Walker A, Garner SF, et al. Functional genomics in zebrafish permits rapid characterization of novel platelet membrane proteins. Blood. (2009) 113:4754–62. doi:  10.1182/blood-2008-06-162693 PubMed DOI PMC

Harel T, Levy-Lahad E, Daana M, Mechoulam H, Horowitz-Cederboim S, Gur M, et al. Homozygous stop-gain variant in LRRC32, encoding a TGFβ receptor, associated with cleft palate, proliferative retinopathy, and developmental delay. Eur J Hum Genet. (2019) 27:1315–9. doi:  10.1038/s41431-019-0380-y PubMed DOI PMC

Lehmkuhl P, Gentz M, Garcia de Otezya AC, Grimbacher B, Schulze-Koops H, Skapenko A. Dysregulated immunity in PID patients with low GARP expression on Tregs due to mutations in LRRC32. Cell Mol Immunol. (2021) 18:1677–91. doi:  10.1038/s41423-021-00701-z PubMed DOI PMC

Qu S, Liu Z, Wang B. Down-regulation of Gremlin1 inhibits inflammatory response and vascular permeability in chronic idiopathic urticaria through suppression of TGF-β signaling pathway. Gene. (2020) 756:144916. doi:  10.1016/j.gene.2020.144916 PubMed DOI

Xie Z, Zhou G, Zhang M, Han J, Wang Y, Li X, et al. Recent developments on BMPs and their antagonists in inflammatory bowel diseases. Cell Death Discovery. (2023) 9:210. doi:  10.1038/s41420-023-01520-z PubMed DOI PMC

Lü AJ, Hu XC, Wang Y, Zhu AH, Shen LL, Tian J, et al. Skin immune response in the zebrafish, PubMed DOI

Du K, Lu F, Xie C, Ding H, Shen Y, Gao Y, et al. PubMed DOI PMC

Ma F, Zhao L, Ma R, Wang J, Du L. FoxO signaling and mitochondria-related apoptosis pathways mediate tsinling lenok trout ( PubMed DOI

Xiang Y, Jia P, Liu W, Yi M, Jia K. Comparative transcriptome analysis reveals the role of p53 signalling pathway during red-spotted grouper nervous necrosis virus infection in PubMed DOI PMC

Karaji N, Sattentau QJ. Efferocytosis of pathogen-infected cells. Front Immunol. (2017) 22:1863. doi:  10.3389/fimmu.2017.01863 PubMed DOI PMC

Jiang X, Sun J, Li C, Hu X, Ge Y, Li B, et al. Molecular cloning and sequence characterization of common carp ( DOI

Podok P, Wang H, Xu L, Xu D, Lu L. Characterization of myeloid-specific peroxidase, keratin 8, and dual specificity phosphatase 1 as innate immune genes involved in the resistance of crucian carp ( PubMed DOI

Shi X, Gao F, Zhao X, Pei C, Zhu L, Zhang J, et al. Role of HIF in fish inflammation. Fish Shellfish Immunol. (2023) 143:109222. doi:  10.1016/j.fsi.2023.109222 PubMed DOI

Kaczorek E, Szarek J, Mikiewicz M, Terech-Majewska E, Schulz P, Małaczewska J, et al. Effect of feed supplementation with kynurenic acid on the morphology of the liver, kidney and gills in rainbow trout ( PubMed DOI

Rojo I, de Ilárduya ÓM, Estonba A, Pardo MA. Innate immune gene expression in individual zebrafish after PubMed DOI

Liu J, Yan Y, Yan J, Wang J, Wei J, Xiao J, et al. Multi-omics analysis revealed crucial genes and pathways associated with black carp antiviral innate immunity. Fish Shellfish Immunol. (2020) 106:724–32. doi:  10.1016/j.fsi.2020.08.047 PubMed DOI

Lu WJ, Zhou L, Gao FX, Zhou YL, Li Z, Zhang XJ, et al. Dynamic and differential expression of duplicated cxcr4/cxcl12 genes facilitates antiviral response in hexaploid gibel carp. Front Immunol. (2020) 11:2176. doi:  10.3389/fimmu.2020.02176 PubMed DOI PMC

Lin CY, Chen YM, Hsu HH, Shiu CT, Kuo HC, Chen TY. Grouper ( PubMed DOI

Aquilino C, Castro R, Fischer U, Tafall C. Transcriptomic responses in rainbow trout gills upon infection with viral hemorrhagic septicemia virus (VHSV). Dev Comp Immunol. (2014) 44:12–20. doi:  10.1016/j.dci.2013.11.006 PubMed DOI

Zaballos A, Gutiérrez J, Varona R, Ardavín C, Márquez G. Cutting edge: Identification of the orphan chemokine receptor GPR-9–6 as CCR9, the receptor for the chemokine TECK. J Immunol. (1999) 162:5671–5. doi:  10.4049/jimmunol.162.10.5671 PubMed DOI

Liu F, Su B, Fu Q, Shang M, Gao C, Tan F, et al. Identification, characterization and expression analysis of TLR5 in the mucosal tissues of turbot ( PubMed DOI

Xue X, Caballero-Solares A, Hall JR, Umasuthan N, Kumar S, Jakob E, et al. Transcriptome profiling of atlantic salmon ( PubMed DOI PMC

Zhan F, Li Y, Shi F, Lu Z, Yang M, Li Q, et al. Characterization analysis of TLR5a and TLR5b immune response after different bacterial infection in grass carp ( PubMed DOI

Cordero H, Mauro M, Cuesta A, Cammarata M, Esteban MA. PubMed DOI

Huang P, Cai J, Yu D, Tang J, Lu Y, Wu Z, et al. An IL-6 gene in humphead snapper ( PubMed DOI

Huang L, Zhao L, Liu W, Xu X, Su Y, Qin Y, et al. Dual RNA-seq unveils PubMed DOI PMC

Ling X, Dong W, Zhang Y, Qian X, Zhang W, He W, et al. Comparative transcriptomics and histopathological analysis of crucian carp infection by atypical PubMed DOI

Valle A, Leiro JM, Pereiro P, Figueras A, Novoa B, Dirks RPH, et al. Interactions between the parasite PubMed DOI PMC

King KC, Seppälä O, Neiman M. Is more better? Polyploidy and parasite resistance. Biol Lett. (2012) 8:598–600. doi:  10.1098/rsbl.2011.1152 PubMed DOI PMC

Guo Z, Liliom K, Fischer DJ, Bathurst IC, Tomei LD, Kiefer MC, et al. Molecular cloning of a high-affinity receptor for the growth factor- like lipid mediator lysophosphatidic acid from PubMed DOI PMC