Evolutionary history of the T cell receptor complex as revealed by small-spotted catshark (Scyliorhinus canicula)

CD4 and LAG-3 from sharks to humans: related molecules with motifs for opposing functions

CD4 and LAG-3 are related molecules that are receptors for MHC class II molecules. Their major functional differences are situated in their cytoplasmic tails, in which CD4 has an activation motif and LAG-3 an inhibitory motif. Here, we identify shark LAG-3 and show that a previously identified shark CD4-like gene has a genomic location, expression pattern, and motifs similar to CD4 in other vertebrates. In nurse shark (Ginglymostoma cirratum) and cloudy catshark (Scyliorhinus torazame), the highest CD4 expression was consistently found in the thymus whereas such was not the case for LAG-3. Throughout jawed vertebrates, the CD4 cytoplasmic tail possesses a Cx(C/H) motif for binding kinase LCK, and the LAG-3 cytoplasmic tail possesses (F/Y)xxL(D/E) including the previously determined FxxL inhibitory motif resembling an immunoreceptor tyrosine-based inhibition motif (ITIM). On the other hand, the acidic end of the mammalian LAG-3 cytoplasmic tail, which is believed to have an inhibitory function as well, was acquired later in evolution. The present study also identified CD4-1, CD4-2, and LAG-3 in the primitive ray-finned fishes bichirs, sturgeons, and gars, and experimentally determined these sequences for sterlet sturgeon (Acipenser ruthenus). Therefore, with CD4-1 and CD4-2 already known in teleosts (modern ray-finned fish), these two CD4 lineages have now been found within all major clades of ray-finned fish. Although different from each other, the cytoplasmic tails of ray-finned fish CD4-1 and chondrichthyan CD4 not only contain the Cx(C/H) motif but also an additional highly conserved motif which we expect to confer a function. Thus, although restricted to some species and gene copies, in evolution both CD4 and LAG-3 molecules appear to have acquired functional motifs besides their canonical Cx(C/H) and ITIM-like motifs, respectively. The presence of CD4 and LAG-3 molecules with seemingly opposing functions from the level of sharks, the oldest living vertebrates with a human-like adaptive immune system, underlines their importance for the jawed vertebrate immune system. It also emphasizes the general need of the immune system to always find a balance, leading to trade-offs, between activating and inhibiting processes.

Harnessing CD3 diversity to optimize CAR T cells

Current US Food and Drug Administration-approved chimeric antigen receptor (CAR) T cells harbor the T cell receptor (TCR)-derived ζ chain as an intracellular activation domain in addition to costimulatory domains. The functionality in a CAR format of the other chains of the TCR complex, namely CD3δ, CD3ε and CD3γ, instead of ζ, remains unknown. In the present study, we have systematically engineered new CD3 CARs, each containing only one of the CD3 intracellular domains. We found that CARs containing CD3δ, CD3ε or CD3γ cytoplasmic tails outperformed the conventional ζ CAR T cells in vivo. Transcriptomic and proteomic analysis revealed differences in activation potential, metabolism and stimulation-induced T cell dysfunctionality that mechanistically explain the enhanced anti-tumor performance. Furthermore, dimerization of the CARs improved their overall functionality. Using these CARs as minimalistic and synthetic surrogate TCRs, we have identified the phosphatase SHP-1 as a new interaction partner of CD3δ that binds the CD3δ-ITAM on phosphorylation of its C-terminal tyrosine. SHP-1 attenuates and restrains activation signals and might thus prevent exhaustion and dysfunction. These new insights into T cell activation could promote the rational redesign of synthetic antigen receptors to improve cancer immunotherapy.

Organized B cell sites in cartilaginous fishes reveal the evolutionary foundation of germinal centers

The absence of germinal centers (GCs) in cartilaginous fishes lies at odds with data showing that nurse sharks can produce robust antigen-specific responses and affinity mature their B cell repertoires. To investigate this apparent incongruity, we performed RNA sequencing on single nuclei, allowing us to characterize the cell types present in the nurse shark spleen, and RNAscope to provide in situ cellular resolution of key marker gene expression following immunization with R-phycoerythrin (PE). We tracked PE to the splenic follicles where it co-localizes with CXCR5high centrocyte-like B cells and a population of putative T follicular helper (Tfh) cells, surrounded by a peripheral ring of Ki67+ AID+ CXCR4+ centroblast-like B cells. Further, we reveal selection of mutations in B cell clones dissected from these follicles. We propose that the B cell sites identified here represent the evolutionary foundation of GCs, dating back to the jawed vertebrate ancestor.

Origin and evolutionary malleability of T cell receptor α diversity

Lymphocytes of vertebrate adaptive immune systems acquired the capability to assemble, from split genes in the germline, billions of functional antigen receptors^1-3. These receptors show specificity; unlike the broadly tuned receptors of the innate system, antibodies (Ig) expressed by B cells, for instance, can accurately distinguish between the two enantiomers of organic acids⁴, whereas T cell receptors (TCRs) reliably recognize single amino acid replacements in their peptide antigens⁵. In developing lymphocytes, antigen receptor genes are assembled from a comparatively small set of germline-encoded genetic elements in a process referred to as V(D)J recombination^6,7. Potential self-reactivity of some antigen receptors arising from the quasi-random somatic diversification is suppressed by several robust control mechanisms^8-12. For decades, scientists have puzzled over the evolutionary origin of somatically diversifying antigen receptors^13-16. It has remained unclear how, at the inception of this mechanism, immunologically beneficial expanded receptor diversity was traded against the emerging risk of destructive self-recognition. Here we explore the hypothesis that in early vertebrates, sequence microhomologies marking the ends of recombining elements became the crucial targets of selection determining the outcome of non-homologous end joining-based repair of DNA double-strand breaks generated during RAG-mediated recombination. We find that, across the main clades of jawed vertebrates, TCRα repertoire diversity is best explained by species-specific extents of such sequence microhomologies. Thus, selection of germline sequence composition of rearranging elements emerges as a major factor determining the degree of diversity of somatically generated antigen receptors.

IL-4/13 expressing CD3γ/δ+ T cells regulate mucosal immunity in response to Flavobacterium columnare infection in grass carp

Interleukin (IL) 4 and 13 are signature cytokines orchestrating Th2 immune response. Teleost fish have two homologs, termed IL-4/13A and IL-4/13B, and have been functionally characterized. However, what cells express IL-4/13A and IL-4/13B has not been investigated in fish. In this work, the recombinant IL-4/13A and IL-4/13B proteins of grass carp (Ctenopharyngodon idella) were produced in the Escherichia coli (E. coli) cells and purified. Monoclonal antibodies (mAbs) against the recombinant CiIL-4/13A and CiIL-4/13B proteins were prepared and characterized. Western blotting analysis showed that the CiIL-4/13A and CiIL-4/13B mAbs could specifically recognize the recombinant proteins expressed in the E. coli cells and HEK293T cells and did not cross-react with each other. Confocal microscopy revealed that the CiIL-4/13A⁺ and CiIL-4/13B⁺ cells were present in the gills, intestine and spleen and could be upregulated in fish infected with Flavobacterium columnare (F. columnare). Interestingly, the cells expressing CiIL-4/13A and CiIL-4/13B were mostly CD3γ/δ⁺ cells. The CD3γ/δ⁺/IL-4/13A⁺ and CD3γ/δ⁺/IL-4/13B⁺ cells were significantly upregulated in the gill filaments and the intestinal mucosa after F. columnare infection. Our results imply that the CD3γ/δ⁺/IL-4/13A⁺ and CD3γ/δ⁺/IL-4/13B⁺ cells are important for homeostasis and the regulation of mucosal immunity.

Pre-Hatching Ontogenetic Changes of Morphological Characters of Small-Spotted Catshark (Scyliorhinus canicula)

The small-spotted catshark, Scyliorhinus canicula, provides an optimal model organism to include chondrichthyans in studies comparing morphology or physiology through vertebrate evolution. In particular, for investigations with ontogenetic aspects, there are only a limited number of alternative taxa. Therefore, a detailed staging system is a prerequisite to allowing comparison between different studies. This study supplements information on the latest stages of the established system by Ballard and colleagues in 1993 and complements the respective staging system by including the latest pre-hatching stages. During this phase, some significant ontogenetic shifts happen, e.g., reduction of external gill filament length and complete flattening of rostral angle until Size Class 6, change in the ratio of pre- to post-vent length, and establishment of body pigmentation in Size Classes 7 and 8. All these shifts finally transform the embryo into a hatchling prepared for living outside the eggshell. This study provides a framework allowing comparison of investigations on pre-hatchings of the small-spotted catshark.

Analysis of shark NCR3 family genes reveals primordial features of vertebrate NKp30

Natural killer (NK) cells play major roles in innate immunity against viruses and cancer. Natural killer receptors (NKR) expressed by NK cells recognize foreign- or self-ligands on infected and transformed cells as well as healthy cells. NKR genes are the most rapidly evolving loci in vertebrates, and it is generally difficult to detect orthologues in different taxa. The unique exception is NKp30, an activating NKR in mammals that binds to the self-ligand B7H6. The NKp30-encoding gene, NCR3, has been found in most vertebrates including sharks, the oldest vertebrates with human-type adaptive immunity. NCR3 has a special, non-rearranging VJ-type immunoglobulin superfamily (IgSF) domain that predates the emergence of the rearranging antigen receptors. Herein we show that NCR3 loci are linked to the shark major histocompatibility complex (MHC), proving NCR3's primordial association with the MHC. We identified eight subtypes of differentially expressed highly divergent shark NCR3 family genes. Using in situ hybridization, we detected one subtype, NS344823, to be expressed by predominantly single cells outside of splenic B cell zones. The expression by non-B cells was also confirmed by PCR in peripheral blood lymphocytes. Surprisingly, high expression of NS344823 was detected in the thymic cortex, demonstrating NS344823 expression in developing T cells. Finally, we show for the first time that shark T cells are found as single cells or in small clusters in the splenic red pulp, also unassociated with the large B cell follicles we previously identified.

Molecular and cellular characterization of European sea bass CD3ε+ T lymphocytes and their modulation by microalgal feed supplementation

The CD3 coreceptor is a master T cell surface marker, and genes encoding CD3ζ, γδ, and ε chains have been reported in several teleost fish. Here, a complete cDNA sequence of CD3ɛ chain was identified from a sea bass (Dicentrarchus labrax L.) gill transcriptome. Its basal expression was quantified in both lymphoid and non-lymphoid organs of sea bass juveniles with real-time qPCR analysis. After either in vitro stimulation of head kidney leukocytes with the T-cell mitogen phytohaemagglutinin or in vivo stimulation with an orally administered Vibrio anguillarum vaccine, CD3ε expression levels increased in head kidney leukocytes, confirming that CD3ε T cells may play important roles in fish systemic protection against pathogens. Further, three peptides were designed on the CD3ɛ cytoplasmic tail region and employed as immunogens for antibody production in rabbit. One antiserum so obtained, named RACD3/1, immunostained a band of the expected size in a western blot of a sea bass thymocyte lysate. The distribution of CD3ε⁺ lymphocyte population in the lymphoid organs and mucosal tissues was addressed in healthy fish by IHC. In decreasing percentage order, CD3ε⁺ lymphocytes were detected by flow cytometry in thymus, peripheral blood leukocytes, gills, head kidney, gut, and spleen. Finally, a significant in vivo enhancement of CD3ε⁺ T intestinal lymphocytes was found in fish fed on diets in which 100% fish meal was replaced by the microalgae Nannochloropsis sp. biomass. These results indicate that CD3ε⁺ T cells are involved in nutritional immune responses.

Characterization of CD3γ/δ+ cells in grass carp (Ctenopharyngodon idella)

CD3 is an essential component of the TCR-CD3 complex which plays a key role in adaptive immunity. Non-mammalian CD3 complex consists of CD3γ/δ, CD3ε and CD3ζ subunits. In this study, homologues of CD3γ/δ and CD3ε (termed CiCD3γ/δ and CiCD3ε) have been identified in grass carp (Ctenopharyngodon idella). Like their counterparts from other vertebrates, the CiCD3γ/δ and CiCD3ε are clustered in the same locus in the genome and encode proteins which are structurally conserved, comprising a signal peptide, an extracellular domain, a transmembrane domain and a cytoplasmic tail containing two ITAM motifs. Sequence analyses identified two novel conserved motifs in the cytoplasmic tail of CiCD3γ/δ and CiCD3ε, one is composed of an arginine and lysine motif (RK or RR) at the C terminus of CiCD3γ/δ and a proline rich domain (PxxPxP/Q) located at the N terminus of ITAM motifs of CiCD3ε. Both genes were highly expressed at the mRNA level in the spleen and gills of healthy fish and could be modulated by infection of Flavobacterium columnare and grass carp reovirus. A monoclonal antibody against the CiCD3γ/δ (GC38T) was produced and showed good reactivity with the native molecule in Western blotting analysis and flow cytometry. The CiCD3γ/δ⁺ cells were analysed in the primary leucocytes, accounting for 5.5% of lymphocytes isolated from spleen, 4.5% from head kidney and 2.8% from peripheral blood. The CiCD3γ/δ⁺ cells were localized in the gills and head kidney by fluorescent confocal microscopy.

Multiple Signaling Roles of CD3ε and Its Application in CAR-T Cell Therapy

A T cell receptor (TCR) mediates antigen-induced signaling through its associated CD3ε, δ, γ, and ζ, but the contributions of different CD3 chains remain elusive. Using quantitative mass spectrometry, we simultaneously quantitated the phosphorylation of the immunoreceptor tyrosine-based activation motif (ITAM) of all CD3 chains upon TCR stimulation. A subpopulation of CD3ε ITAMs was mono-phosphorylated, owing to Lck kinase selectivity, and specifically recruited the inhibitory Csk kinase to attenuate TCR signaling, suggesting that TCR is a self-restrained signaling machinery containing both activating and inhibitory motifs. Moreover, we found that incorporation of the CD3ε cytoplasmic domain into a second-generation chimeric antigen receptor (CAR) improved antitumor activity of CAR-T cells. Mechanistically, the Csk-recruiting ITAM of CD3ε reduced CAR-T cytokine production whereas the basic residue rich sequence (BRS) of CD3ε promoted CAR-T persistence via p85 recruitment. Collectively, CD3ε is a built-in multifunctional signal tuner, and increasing CD3 diversity represents a strategy to design next-generation CAR.

Adverse effects of plastic ingestion on the Mediterranean small-spotted catshark (Scyliorhinus canicula)

Plastics are widely diffused in the oceans and their ingestion by marine organisms is raising concern for potentially adverse effects. The risk of harmful interactions with marine plastic pollution depends on the biology of the species as well as the distribution and abundance of the different plastic types. The aim of this study was to assess the occurrence of plastic ingestion by the small-spotted catshark (Scyliorhinus canicula), one of the most abundant elasmobranchs in the Mediterranean Sea. The expression levels of genes indicative of total immune system function were analyzed to gather preliminary data for further investigation of any potential correlations between plastic presence and immune activation. One hundred catsharks were collected during the Spring 2018 in two geographic locations in the southern region of the central Mediterranean Sea: 1) near Mazara del Vallo, SW Sicily and 2) near Lampedusa island, Italy's southernmost. Standard measurements were recorded for each specimen and its organs and sex was determined. The gastrointestinal tract (GIT) was preserved for plastic detection and identification. Where present, plastics (macro- and micro-) were characterized in terms of size, shape and polymer typology through microscopy and μ-Raman spectroscopy. Spleen from a subset of thirty samples was preserved for RNA extraction, then used to quantify by real time PCR the transcripts of T cell receptor beta (TCRB), T cell receptor delta (TCRD) and IgM genes. The results indicated that ingestion of plastic is widespread, with microplastics (MP, from 1 μm to <1 mm) abundantly present in nearly all samples and macroplasticplastic (MaP, > 1 cm) in approximately 18% of the specimens collected. A significant increase in the expression of TCRB, TCRD and IgM was observed in the spleen of MaP + specimens from Mazara del Vallo waters, in parallel with 67% increase in liver weight. While the presence of MP alone is not enough to induce a strong activation of the immunity, some type of plastics falling into the MaP category may be more toxic than others and crucial in the activation of the immune response. The results of this study represent a first evidence that plastic pollution represents an emerging threat to S. canicula, the Mediterranean food web and human consumers.

Teleost cytotoxic T cells

Cell-mediated cytotoxicity is one of the major mechanisms by which vertebrates control intracellular pathogens. Two cell types are the main players in this immune response, natural killer (NK) cells and cytotoxic T lymphocytes (CTL). While NK cells recognize altered target cells in a relatively unspecific manner CTLs use their T cell receptor to identify pathogen-specific peptides that are presented by major histocompatibility (MHC) class I molecules on the surface of infected cells. However, several other signals are needed to regulate cell-mediated cytotoxicity involving a complex network of cytokine- and ligand-receptor interactions. Since the first description of MHC class I molecules in teleosts during the early 90s of the last century a remarkable amount of information on teleost immune responses has been published. The corresponding studies describe teleost cells and molecules that are involved in CTL responses of higher vertebrates. These studies are backed by functional investigations on the killing activity of CTLs in a few teleost species. The present knowledge on teleost CTLs still leaves considerable room for further investigations on the mechanisms by which CTLs act. Nevertheless the information on teleost CTLs and their regulation might already be useful for the control of fish diseases by designing efficient vaccines against such diseases where CTL responses are known to be decisive for the elimination of the corresponding pathogen. This review summarizes the present knowledge on CTL regulation and functions in teleosts. In a special chapter, the role of CTLs in vaccination is discussed.

A Comparison of the Innate and Adaptive Immune Systems in Cartilaginous Fish, Ray-Finned Fish, and Lobe-Finned Fish

The immune system is composed of two subsystems-the innate immune system and the adaptive immune system. The innate immune system is the first to respond to pathogens and does not retain memory of previous responses. Innate immune responses are evolutionarily older than adaptive responses and elements of innate immunity can be found in all multicellular organisms. If a pathogen persists, the adaptive immune system will engage the pathogen with specificity and memory. Several components of the adaptive system including immunoglobulins (Igs), T cell receptors (TCR), and major histocompatibility complex (MHC), are assumed to have arisen in the first jawed vertebrates-the Gnathostomata. This review will discuss and compare components of both the innate and adaptive immune systems in Gnathostomes, particularly in Chondrichthyes (cartilaginous fish) and in Osteichthyes [bony fish: the Actinopterygii (ray-finned fish) and the Sarcopterygii (lobe-finned fish)]. While many elements of both the innate and adaptive immune systems are conserved within these species and with higher level vertebrates, some elements have marked differences. Components of the innate immune system covered here include physical barriers, such as the skin and gastrointestinal tract, cellular components, such as pattern recognition receptors and immune cells including macrophages and neutrophils, and humoral components, such as the complement system. Components of the adaptive system covered include the fundamental cells and molecules of adaptive immunity: B lymphocytes (B cells), T lymphocytes (T cells), immunoglobulins (Igs), and major histocompatibility complex (MHC). Comparative studies in fish such as those discussed here are essential for developing a comprehensive understanding of the evolution of the immune system.

Discovery of All Three Types in Cartilaginous Fishes Enables Phylogenetic Resolution of the Origins and Evolution of Interferons

Interferons orchestrate host antiviral responses in jawed vertebrates. They are categorized into three classes; IFN1 and IFN3 are the primary antiviral cytokine lineages, while IFN2 responds to a broader variety of pathogens. The evolutionary relationships within and between these three classes have proven difficult to resolve. Here, we reassess interferon evolution, considering key phylogenetic pitfalls including taxon sampling, alignment quality, model adequacy, and outgroup choice. We reveal that cartilaginous fishes, and hence the jawed vertebrate ancestor, possess(ed) orthologs of all three interferon classes. We show that IFN3 groups sister to IFN1, resolve the origins of the human IFN3 lineages, and find that intronless IFN3s emerged at least three times. IFN2 genes are highly conserved, except for IFN-γ-rel, which we confirm resulted from a teleost-specific duplication. Our analyses show that IFN1 phylogeny is highly sensitive to phylogenetic error. By accounting for this, we describe a new backbone IFN1 phylogeny that implies several IFN1 genes existed in the jawed vertebrate ancestor. One of these is represented by the intronless IFN1s of tetrapods, including mammalian-like repertoires of reptile IFN1s and a subset of amphibian IFN1s, in addition to newly-identified intron-containing shark IFN1 genes. IFN-f, previously only found in teleosts, likely represents another ancestral jawed vertebrate IFN1 family member, suggesting the current classification of fish IFN1s into two groups based on the number of cysteines may need revision. The providence of the remaining fish IFN1s and the coelacanth IFN1s proved difficult to resolve, but they may also be ancestral jawed vertebrate IFN1 lineages. Finally, a large group of amphibian-specific IFN1s falls sister to all other IFN1s and was likely also present in the jawed vertebrate ancestor. Our results verify that intronless IFN1s have evolved multiple times in amphibians and indicate that no one-to-one orthology exists between mammal and reptile IFN1s. Our data also imply that diversification of the multiple IFN1s present in the jawed vertebrate ancestor has occurred through a rapid birth-death process, consistent with functional maintenance over a 450-million-year host-pathogen arms race. In summary, this study reveals a new model of interferon evolution important to our understanding of jawed vertebrate antiviral immunity.

Molecular characterization and expression analysis of T cell receptor (TCR) γ and δ genes in dojo loach (Misgurnus anguillicaudatus) in response to bacterial, parasitic and fungal challenge

In mammalian, T-cell receptors (TCRs) play a key role in recognizing the presented antigen from external to protect organisms against environmental pathogens. To understand the potential roles of TCRγ and TCRδ in dojo loach (Misgurnus anguillicaudatus), Ma-TCRγ and Ma-TCRδ cDNAs were cloned and their gene expression profiles were investigated after bacterial, parasitic and fungal challenge. The open reading frame (ORF) of Ma-TCRγ and Ma-TCRδ cDNAs contained 948 and 867 bp, encoding 316 and 288 amino acid residues, respectively. Structurally, Ma-TCRγ and Ma-TCRδ were consisted of a signal peptide, a variable region, a constant region (IgC), a connecting peptide (CPS), a transmembrane region (TM) and a cytoplasmic domain (CYT), which were similar to those of other vertebrates. Multiple sequence alignment and phylogenetic analysis showed Ma-TCRγ and Ma-TCRδ were closely related to fish of Cyprinidae family. Ma-TCRγ and Ma-TCRδ were widely expressed in all tested organs/tissues, as the highest expressions of Ma-TCRγ and Ma-TCRδ were detected in kidney and gill, respectively. In addition, three infection models of dojo loach with bacteria (F. columnare G₄), parasite (Ichthyophthirius multifiliis) and fungus (Saprolegnia sp.) were constructed. The morphological changes of gills and skin after challenged with F. columnare G₄ and Ichthyophthirius multifiliis were investigated. Compared to F. columnare G₄ infection, mRNA expression of both TCRγ and TCRδ showed higher sensitivity in classical immune organs (kidney and spleen) and mucosal tissues (skin and gill) after challenge with Ichthyophthirius multifiliis and Saprolegnia sp. Our results first indicated that TCRγ and TCRδ of dojo loach might function differently in response to challenge with different pathogens.

Phylotranscriptomics suggests the jawed vertebrate ancestor could generate diverse helper and regulatory T cell subsets

BACKGROUND

The cartilaginous fishes diverged from other jawed vertebrates ~ 450 million years ago (mya). Despite this key evolutionary position, the only high-quality cartilaginous fish genome available is for the elephant shark (Callorhinchus milii), a chimaera whose ancestors split from the elasmobranch lineage ~ 420 mya. Initial analysis of this resource led to proposals that key components of the cartilaginous fish adaptive immune system, most notably their array of T cell subsets, was primitive compared to mammals. This proposal is at odds with the robust, antigen-specific antibody responses reported in elasmobranchs following immunization. To explore this discrepancy, we generated a multi-tissue transcriptome for small-spotted catshark (Scyliorhinus canicula), a tractable elasmobranch model for functional studies. We searched this, and other newly available sequence datasets, for CD4+ T cell subset-defining genes, aiming to confirm the presence or absence of each subset in cartilaginous fishes.

RESULTS

We generated a new transcriptome based on a normalised, multi-tissue RNA pool, aiming to maximise representation of tissue-specific and lowly expressed genes. We utilized multiple transcriptomic datasets and assembly variants in phylogenetic reconstructions to unambiguously identify several T cell subset-specific molecules in cartilaginous fishes for the first time, including interleukins, interleukin receptors, and key transcription factors. Our results reveal the inability of standard phylogenetic reconstruction approaches to capture the site-specific evolutionary processes of fast-evolving immune genes but show that site-heterogeneous mixture models can adequately do so.

CONCLUSIONS

Our analyses reveal that cartilaginous fishes are capable of producing a range of CD4+ T cell subsets comparable to that of mammals. Further, that the key molecules required for the differentiation and functioning of these subsets existed in the jawed vertebrate ancestor. Additionally, we highlight the importance of considering phylogenetic diversity and, where possible, utilizing multiple datasets for individual species, to accurately infer gene presence or absence at higher taxonomic levels.