Our understanding of the vertebrate immune system is dominated by a few model organisms such as mice. This use of a few model systems is reasonable if major features of the immune systems evolve slowly and are conserved across most vertebrates, but may be problematic if there is substantial macroevolutionary change in immune responses. Here, we present a test of the macroevolutionary stability, across 14 species of ray-finned fishes, of the transcriptomic response to a standardized immune challenge. Intraperitoneal injection of an immune adjuvant (alum) induces a fibrosis response in nearly all jawed fishes, which in some species contributes to anti-helminth protection. Despite this conserved phenotypic response, the underlying transcriptomic response is highly inconsistent across species. Although many gene orthogroups exhibit differential expression between saline versus alum-injected fish in at least one species, few orthogroups exhibit consistent differential expression across species. This result suggests that although the phenotypic response to alum (fibrosis) is highly conserved, the underlying gene regulatory architecture is very flexible and cannot readily be extrapolated from any one species to fishes (or vertebrates) more broadly. The vertebrate immune response is remarkably changeable over macroevolutionary time, requiring a diversity of model organisms to describe effectively.

Department of Biology Texas State University San Marcos TX 78666 USA

Department of Ecology and Evolutionary Biology University of Connecticut Storrs CT 06269 USA

Institute of Vertebrate Biology The Czech Academy of Sciences Brno 603 65 Czech Republic

Před aktualizací

bioRxiv. 2024 Dec 17:2024.12.13.628417. doi: 10.1101/2024.12.13.628417. PubMed

Zobrazit více v PubMed

Agius  C, Roberts  RJ. Melano-macrophage centres and their role in fish pathology. J Fish Dis. 2003:26:499–509. 10.1046/j.1365-2761.2003.00485.x. PubMed DOI

Alfaro  ME, Bolnick  DI, Wainwright  PC. Evolutionary consequences of many-to-one mapping of jaw morphology to mechanics in labrid fishes. Am Nat. 2005:165:E140–E154. 10.1086/429564. PubMed DOI

Bernhagen  J  et al.  MIF is a noncognate ligand of CXC chemokine receptors in inflammatory and atherogenic cell recruitment. Nat Med. 2007:13:587–596. 10.1038/nm1567. PubMed DOI

Bjornson-Hooper  ZB  et al.  A comprehensive atlas of immunological differences between humans, mice, and non-human primates. Front Immunol. 2022:13:867015. 10.3389/fimmu.2022.867015. PubMed DOI PMC

Bolnick  DI  et al.  The dominance of coinfecting parasites' indirect genetic effects on host traits. Am Nat. 2024:204:482–500. 10.1086/732256. PubMed DOI PMC

Boyle  EA, Li  YI, Pritchard  JK. An expanded view of complex traits: from polygenic to omnigenic. Cell. 2017:169:1177–1186. 10.1016/j.cell.2017.05.038. PubMed DOI PMC

Brownstein  CD  et al.  Synergistic innovations enabled the radiation of anglerfishes in the deep open ocean. Curr Biol. 2024:34:2541–2550.e4. 10.1016/j.cub.2024.04.066. PubMed DOI

Bucchini  F  et al.  TRAPID 2.0: a web application for taxonomic and functional analysis of de novo transcriptomes. Nucleic Acids Res. 2021:49:e101. 10.1093/nar/gkab565. PubMed DOI PMC

Buchfink  B, Xie  C, Huson  DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015:12:59–60. 10.1038/nmeth.3176. PubMed DOI

Bushmanova  E, Antipov  D, Lapidus  A, Prjibelski  AD. rnaSPAdes: a de novo transcriptome assembler and its application to RNA-Seq data. Gigascience. 2019:8:giz100. 10.1093/gigascience/giz100. PubMed DOI PMC

Chen  F  et al.  Phylogenetic comparative analysis of single-cell transcriptomes reveals constrained accumulation of gene expression heterogeneity during clonal expansion. Mol Biol Evol. 2023:40:msad113. 10.1093/molbev/msad113. PubMed DOI PMC

Colquitt  BM, Merullo  DP, Konopka  G, Roberts  TF, Brainard  MS. Cellular transcriptomics reveals evolutionary identities of songbird vocal circuits. Science. 2021:371:eabd9704. 10.1126/science.abd9704. PubMed DOI PMC

De Lisle  SP, Bolnick  DI. Male and female reproductive fitness costs of an immune response in natural populations. Evolution. 2021:75:2509–2523. 10.1111/evo.14266. PubMed DOI PMC

de Vienne  DM, Giraud  T, Martin  OC. A congruence index for testing topological similarity between trees. Bioinformatics. 2007:23:3119–3124. 10.1093/bioinformatics/btm500. PubMed DOI

Domínguez Conde  C  et al.  Cross-tissue immune cell analysis reveals tissue-specific features in humans. Science. 2022:376:eabl5197. 10.1126/science.abl5197. PubMed DOI PMC

Dubin  A, Jørgensen  TE, Moum  T, Johansen  SD, Jakt  LM. Complete loss of the MHC II pathway in an anglerfish, Lophius piscatorius. Biol Lett. 2019:15:20190594. 10.1098/rsbl.2019.0594. PubMed DOI PMC

Emms  DM, Kelly  S. STRIDE: species tree root inference from gene duplication events. Mol Biol Evol. 2017:34:3267–3278. 10.1093/molbev/msx259. PubMed DOI PMC

Emms  DM, Kelly  S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019:20:238. 10.1186/s13059-019-1832-y. PubMed DOI PMC

Fortier  AL, Pritchard  JK. The primate Major Histocompatibility Complex as a case study of gene family evolution. eLife. 2025:14:RP103545. 10.7554/eLife.103545.3 . PubMed DOI PMC

Fuess  LE  et al.  Between-population differences in constitutive and infection-induced gene expression in threespine stickleback. Mol Ecol. 2021:30:6791–6805. 10.1111/mec.16197. PubMed DOI PMC

Gilbertson  SE, Weinmann  AS. Conservation and divergence in gene regulation between mouse and human immune cells deserves equal emphasis. Trends Immunol. 2021:42:1077–1087. 10.1016/j.it.2021.10.007. PubMed DOI PMC

Greenway  R  et al.  Integrative analyses of convergent adaptation in sympatric extremophile fishes. Curr Biol. 2024:34:4968–4982.e7. 10.1016/j.cub.2024.09.027. PubMed DOI

Gros  P, Casanova  J-L. Reconciling mouse and human immunology at the altar of genetics. Annu Rev Immunol. 2023:41:39–71. 10.1146/annurev-immunol-101721-065201. PubMed DOI

Haas  BJ  et al.  De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013:8:1494–1512. 10.1038/nprot.2013.084. PubMed DOI PMC

Hadfield  JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J Stat Softw. 2010:33:1–22. 10.18637/jss.v033.i02. DOI

Henderson  N, Rieder  F, Wynn  T. Fibrosis: from mechanisms to medicines. Nature. 2020:587:555–566. 10.1038/s41586-020-2938-9. PubMed DOI PMC

Hou  X  et al.  Transcriptome analysis reveals the neuroactive receptor genes response to Streptococcus agalactiae infection in tilapia, Oreochromis niloticus. Fish Shellfish Immunol. 2023:141:109090. 10.1016/j.fsi.2023.109090. PubMed DOI

Huang  Y  et al.  Transcriptome profiling of immune tissues reveals habitat-specific gene expression between lake and river sticklebacks. Mol Ecol. 2016:25:943–958. 10.1111/mec.13520. PubMed DOI PMC

Hughes  LC  et al.  Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data. Proc Natl Acad Sci U S A. 2018:115:6249–6254. 10.1073/pnas.1719358115. PubMed DOI PMC

Hund  AK  et al.  Population-level variation in parasite resistance due to differences in immune initiation and rate of response. Evol Lett. 2022:6:162–177. 10.1002/evl3.274. PubMed DOI PMC

Kool  M, Fierens  K, Lambrecht  BN. Alum adjuvant: some of the tricks of the oldest adjuvant. J Med Microbiol. 2012:61:927–934. 10.1099/jmm.0.038943-0. PubMed DOI

Kuo  C-L  et al.  APOE e4 genotype predicts severe COVID-19 in the UK biobank community cohort. J Gerontol A Biol Sci Med Sci. 2020:75:2231–2232. 10.1093/gerona/glaa131. PubMed DOI PMC

Liu  J, Yuan  X, Fan  C, Ma  G. Application of the zebrafish model in human viral research. Virus Res. 2024:341:199327. 10.1016/j.virusres.2024.199327. PubMed DOI PMC

Lohman  BK, Steinel  NC, Weber  JN, Bolnick  DI. Gene expression contributes to the recent evolution of host resistance in a model host parasite system. Front Immunol. 2017:8:1071. 10.3389/fimmu.2017.01071. PubMed DOI PMC

Lohman  BK, Weber  JN, Bolnick  DI. Evaluation of TagSeq, a reliable low-cost alternative for RNAseq. Mol Ecol Resour. 2016:16:1315–1321. 10.1111/1755-0998.12529. PubMed DOI

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

Malmstrøm  M  et al.  Evolution of the immune system influences speciation rates in teleost fishes. Nat Genet. 2016:48:1204–1210. 10.1038/ng.3645. PubMed DOI

Matthews  DG, Maciejewski  MF, Wong  GA, Lauder  GV, Bolnick  DI. Locomotor effects of a fibrosis-based immune response in stickleback fish. J Exp Biol. 2023:226:jeb246684. 10.1242/jeb.246684. PubMed DOI

Meyer  E, Aglyamova  GV, Matz  MV. Profiling gene expression responses of coral larvae (Acropora millepora) to elevated temperature and settlement inducers using a novel RNA-Seq procedure. Mol Ecol. 2011:20:3599–3616. 10.1111/j.1365-294X.2011.05205.x. PubMed DOI

O’Connor  EA, Cornwallis  CK, Hasselquist  D, Nilsson  J-Å, Westerdahl  H. The evolution of immunity in relation to colonization and migration. Nat Ecol Evol. 2018:2:841–849. 10.1038/s41559-018-0509-3. PubMed DOI

Ostendorf  BN  et al.  Common human genetic variants of APOE impact murine COVID-19 mortality. Nature. 2022:611:346–351. 10.1038/s41586-022-05344-2. PubMed DOI PMC

Pasquier  J  et al.  Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database. BMC Genomics. 2016:17:368. 10.1186/s12864-016-2709-z. PubMed DOI PMC

Patro  R, Duggal  G, Love  MI, Irizarry  RA, Kingsford  C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017:14:417–419. 10.1038/nmeth.4197. PubMed DOI PMC

R Core Team . R: a language and environment for statistical computing. 2024. https://www.R-project.org/

Schiffman  JS, Ralph  PL. System drift and speciation. Evolution. 2022:76:236–251. 10.1111/evo.14356. PubMed DOI PMC

Singh  KK  et al.  Human SAP18 mediates assembly of a splicing regulatory multiprotein complex via its ubiquitin-like fold. RNA. 2010:16:2442–2454. 10.1261/rna.2304410. PubMed DOI PMC

Soto  E  et al.  Expression of immune markers of white sturgeon (Acipenser transmontanus) during Veronaea botryosa infection at different temperatures. Comp Biochem Physiol Part D Genomics Proteomics. 2022:41:100950. 10.1016/j.cbd.2021.100950. PubMed DOI

Star  B  et al.  The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011:477:207–210. 10.1038/nature10342. PubMed DOI PMC

Sun  T  et al.  Integrative mRNA-miRNA interaction analysis associated with the immune response in the head kidney of rainbow trout (Oncorhynchus mykiss) after infectious hematopoietic necrosis virus infection. Fish Shellfish Immunol. 2023:142:109140. 10.1016/j.fsi.2023.109140. PubMed DOI

Swann  JB, Holland  SJ, Petersen  M, Pietsch  TW, Boehm  T. The immunogenetics of sexual parasitism. Science. 2020:369:1608–1615. 10.1126/science.aaz9445. PubMed DOI

True  JR, Haag  ES. Developmental system drift and flexibility in evolutionary trajectories. Evol Dev. 2001:3:109–119. 10.1046/j.1525-142x.2001.003002109.x. PubMed DOI

Vrtílek  M, Bolnick  DI. Macroevolutionary foundations of a recently evolved innate immune defense. Evolution. 2021:75:2600–2612. 10.1111/evo.14316. PubMed DOI PMC

Wainwright  PC, Alfaro  ME, Bolnick  DI, Hulsey  CD. Many-to-one mapping of form to function: a general principle in organismal design?  Integr Comp Biol. 2005:45:256–262. 10.1093/icb/45.2.256. PubMed DOI

Weber  JN  et al.  Evolutionary gain and loss of a pathological immune response to parasitism. Science. 2022:377:1206–1211. 10.1126/science.abo3411. PubMed DOI PMC

Wegner  KM, Reusch  TBH, Kalbe  M. Multiple parasites are driving major histocompatibility complex polymorphism in the wild: multiple parasites drive MHC polymorphism. J Evol Biol. 2003:16:224–232. 10.1046/j.1420-9101.2003.00519.x. PubMed DOI

Wiley  EO, Lieberman  BS. Phylogenetics: theory and practice of phylogenetic systematics. 2nd ed. Wiley-Blackwell; 2011.

Wuitchik  DM  et al.  Seasonal temperature, the lunar cycle and diurnal rhythms interact in a combinatorial manner to modulate genomic responses to the environment in a reef-building coral. Mol Ecol. 2019:28:3629–3641. 10.1111/mec.15173. PubMed DOI PMC

Wynn  TA, Ramalingam  TR. Mechanisms of fibrosis: therapeutic translation for fibrotic disease. Nat Med. 2012:18:1028–1040. 10.1038/nm.2807. PubMed DOI PMC

Zapata  A, Diez  B, Cejalvo  T, Gutiérrez-de Frías  C, Cortés  A. Ontogeny of the immune system of fish. Fish Shellfish Immunol. 2006:20:126–136. 10.1016/j.fsi.2004.09.005. PubMed DOI