Root-level loss of immunoglobulin and B-cell immune genes in clingfishes

Abstract

Immunoglobulin genes are a central component of jawed-vertebrate adaptive immunity. A previous study showed that the blunt-snouted clingfish Gouania willdenowi lacks immunoglobulin genes and T-cell receptor gamma/delta loci, while retaining T-cell receptor alpha/beta genes, MHC genes, and RAG1 /RAG2. Here I extend that observation to the family Gobiesocidae using all seven chromosome-level Gobiesocidae genome assemblies currently available. Manual tblastn and synteny-guided searches found no convincing immunoglobulin heavy-chain or light-chain loci in G. willdenowi, Gouania pigra, Gobiesox punctulatus, Apletodon dentatus, Lepadogaster candolii, Lepadogaster purpurea, or Diplecogaster bimaculata. Thus, the absence of antibody genes is best interpreted as a root-level character of clingfishes. The latest seven-species screen of 40 additional immune-associated genes shifts the broader interpretation in the same direction: the B-cell/adaptive core genes CD79A, CD79B, CIITA, TNFRSF13B, and TNFSF13B lack strong tblastn support in all sampled Gobiesocidae, and 37 of the 40 tested targets show an all-zero binary pattern at the presence threshold. Only IL21R.1, TYROBP, and TNFRSF11A show strong hits in one or more species. I therefore interpret the principal immune-gene erosion as occurring at or near the Gobiesocidae root rather than as a recent Gouania-specific process, while keeping weak, paralog-sensitive, and patchy loci provisional. RAG2 comparisons show a shared Gobiesocidae PHD-domain C-to-S replacement in the zinc-binding motif, with apparently intact RAG2 coding sequence. A family-wide TRG/TRD screen did not recover TRGV V segments or accepted TRDC constant-region exons, but it did detect TRGC-like constant exons in several genomes. These TRGC-like sequences are probably not canonical TRG constant exons without further validation, so I treat the gamma/delta system as eroded or rearranged rather than as a complete root-level loss equivalent to the Ig loss. The RAG2 variant provides a plausible molecular context for antigen-receptor remodeling, but it is not evidence that RAG genes are pseudogenized, because TCR alpha/beta, MHC genes, and RAG1 /RAG2 are retained. Gobiesocidae are therefore best described as a vertebrate family with ancestral loss of canonical immunoglobulin genes and associated root-level erosion of B-cell and immune-related genes, not as a lineage lacking adaptive immunity in its entirety. HighlightsO_LISeven chromosome-level Gobiesocidae genomes lack convincing canonical IgH and IgL loci. C_LIO_LIThe strongest non-Ig losses map to the B-cell/adaptive core: CD79A, CD79B, CIITA, TNFRSF13B, and TNFSF13B. C_LIO_LITCR alpha/beta, MHC genes, and RAG1 /RAG2 are retained, so Gobiesocidae should not be described as lacking adaptive immunity in full. C_LIO_LIA shared Gobiesocidae RAG2 PHD-domain C-to-S variant provides candidate molecular context for antigen-receptor remodeling. C_LI