Costimulatory receptors in the channel catfish: CD28 family members and their ligands.

The CD28-B7 interaction is required to deliver a second signal necessary for T-cell activation. Additional membrane receptors of the CD28 and B7 families are also involved in immune checkpoints that positively or negatively regulate leukocyte activation, in particular T lymphocytes. BTLA is an inhibitory receptor that belongs to a third receptor family. Fish orthologs exist only for some of these genes, and the potential interactions between the corresponding ligands remain mostly unclear. In this work, we focused on the channel catfish (Ictalurus punctatus), a long-standing model for fish immunology, to analyze these co-stimulatory and co-inhibitory receptors. We identified one copy of cd28, ctla4, cd80/86, b7h1/dc, b7h3, b7h4, b7h5, two btla, and four b7h7 genes. Catfish CD28 contains the highly conserved mammalian cytoplasmic motif for PI3K and GRB2 recruitment, however this motif is absent in cyprinids. Fish CTLA4 share a C-terminal putative GRB2-binding site but lacks the mammalian PI3K/GRB2-binding motif. While critical V-domain residues for human CD80 or CD86 binding to CD28/CTLA4 show low conservation in fish CD80/86, C-domain residues are highly conserved, underscoring their significance. Catfish B7H1/DC had a long intracytoplasmic domain with a P-loop-NTPase domain that is absent in mammalian sequences, while the lack of NLS motif in fish B7H4 suggests this protein may not regulate cell growth when expressed intracellularly. Finally, there is a notable expansion of fish B7H7s, which likely play diverse roles in leukocyte regulation. Overall, our work contributes to a better understanding of fish leukocyte co-stimulatory and co-inhibitory receptors.

Identification of a pharyngeal mucosal lymphoid organ in zebrafish and other teleosts: Tonsils in fish?

The constant exposure of the fish branchial cavity to aquatic pathogens causes local mucosal immune responses to be extremely important for their survival. Here, we used a marker for T lymphocytes/natural killer (NK) cells (ZAP70) and advanced imaging techniques to investigate the lymphoid architecture of the zebrafish branchial cavity. We identified a sub-pharyngeal lymphoid organ, which we tentatively named "Nemausean lymphoid organ" (NELO). NELO is enriched in T/NK cells, plasma/B cells, and antigen-presenting cells embedded in a network of reticulated epithelial cells. The presence of activated T cells and lymphocyte proliferation, but not V(D)J recombination or hematopoiesis, suggests that NELO is a secondary lymphoid organ. In response to infection, NELO displays structural changes including the formation of T/NK cell clusters. NELO and gill lymphoid tissues form a cohesive unit within a large mucosal lymphoid network. Collectively, we reveal an unreported mucosal lymphoid organ reminiscent of mammalian tonsils that evolved in multiple teleost fish families.

Cold-blooded vertebrates evolved organized germinal center-like structures.

Germinal centers (GCs) or analogous secondary lymphoid microstructures (SLMs) are thought to have evolved in endothermic species. However, living representatives of their ectothermic ancestors can mount potent secondary antibody responses upon infection or immunization, despite the apparent lack of SLMs in these cold-blooded vertebrates. How and where adaptive immune responses are induced in ectothermic species in the absence of GCs or analogous SLMs remain poorly understood. Here, we infected a teleost fish (trout) with the parasite Ichthyophthirius multifiliis (Ich) and identified the formation of large aggregates of highly proliferating IgM+ B cells and CD4+ T cells, contiguous to splenic melanomacrophage centers (MMCs). Most of these MMC-associated lymphoid aggregates (M-LAs) contained numerous antigen (Ag)-specific B cells. Analysis of the IgM heavy chain CDR3 repertoire of microdissected splenic M-LAs and non-M-LA areas revealed that the most frequent B cell clones induced after Ich infection were highly shared only within the M-LAs of infected animals. These M-LAs represented highly polyclonal SLMs in which Ag-specific B cell clonal expansion occurred. M-LA-associated B cells expressed high levels of activation-induced cytidine deaminase and underwent significant apoptosis, and somatic hypermutation of Igµ genes occurred prevalently in these cells. Our findings demonstrate that ectotherms evolved organized SLMs with GC-like roles. Moreover, our results also point to primordially conserved mechanisms by which M-LAs and mammalian polyclonal GCs develop and function.

Evolution of T cell receptor beta loci in salmonids.

T-cell mediated immunity relies on a vast array of antigen specific T cell receptors (TR). Characterizing the structure of TR loci is essential to study the diversity and composition of T cell responses in vertebrate species. The lack of good-quality genome assemblies, and the difficulty to perform a reliably mapping of multiple highly similar TR sequences, have hindered the study of these loci in non-model organisms. High-quality genome assemblies are now available for the two main genera of Salmonids, Salmo and Oncorhynchus. We present here a full description and annotation of the TRB loci located on chromosomes 19 and 25 of rainbow trout (Oncorhynchus mykiss). To get insight about variations of the structure and composition of TRB locus across salmonids, we compared rainbow trout TRB loci with other salmonid species and confirmed that the basic structure of salmonid TRB locus is a double set of two TRBV-D-J-C loci in opposite orientation on two different chromosomes. Our data shed light on the evolution of TRB loci in Salmonids after their whole genome duplication (WGD). We established a coherent nomenclature of salmonid TRB loci based on comprehensive annotation. Our work provides a fundamental basis for monitoring salmonid T cell responses by TRB repertoire sequencing.

Conserved and divergent arms of the antiviral response in the duplicated genomes of salmonid fishes.

Antiviral innate immunity is orchestrated by the interferon system, which appeared in ancestors of jawed vertebrates. Interferon upregulation induces hundreds of interferon-stimulated-genes (ISGs) with effector or regulatory functions. Here we investigated the evolutionary diversification of ISG responses through comparison of two salmonid fishes, accounting for the impact of sequential whole genome duplications ancestral to teleosts and salmonids. We analysed the transcriptomic response of the IFN pathway in the head kidney of rainbow trout and Atlantic salmon, which separated 25-30 Mya. We identified a large set of ISGs conserved in both species and cross-referenced them with zebrafish and human ISGs. In contrast, around one-third of salmonid ISG lacked orthologs in human, mouse, chicken or frog, and often between rainbow trout and Atlantic salmon, revealing a fast-evolving, lineage-specific arm of the antiviral response. This study also provides a key resource for in-depth functional analysis of ISGs in salmonids of commercial significance.

Differentiation and traffic of IgM+ B cells between focal dark spots in skeletal muscle of Atlantic salmon, lymphoid and adipose tissues.

Focal dark spots (DS) in farmed Atlantic salmon fillets contain a significant number of B cells as revealed by the high abundance of immunoglobulin (Ig) transcripts in transcriptome data. The immune response in DS remains unknown while they represent a major problem in commercial aquaculture. Here, we characterized the diversity and clonal composition of B cells in DS. Sixteen gene markers of immune cells and antigen presentation were analyzed with RT-qPCR. All genes expression showed a positive correlation with DS area and intensity. The flatter the DS, the higher the expression of cd28, csfr, ctla, igt, and sigm, the lower expression of cd83 and btla, and the larger the cumulative frequency within DS. The expression of most of the analyzed immune genes, including three Ig types and markers of B cells was lower in DS than in the lymphatic organs, head kidney and spleen, but significantly higher compared to skeletal muscle. High levels of ctla4 and cd28 in DS might indicate the recruitment of T cells. Sequencing of IgM repertoire (Ig-seq) assessed migration of B cells by co-occurrence of identical CDR3 sequences in different tissues. The combination of gene expression and Ig-seq revealed the presence of several stages of B cell differentiation in DS. B cells at the earliest stage, with high ratio of membrane to secretory IgM (migm and sigm), showed minor Ig repertoire overlap with other tissues. Further differentiation stage (increased sigm to migm ratio and high expression of pax5 and cd79) was associated with active movement of B cells from DS towards lymphatic organs and visceral fat. Traffic and expression of immune genes decreased at later stages. These B cells could be involved in a response directed against viruses, pathogenic or opportunistic bacteria in DS. Seven of eight fish were positive for salmon alphavirus, and levels were higher in DS than in unstained muscle. PCR with universal primers to the 16S rRNA gene did not detect bacteria in DS. Although the evolution of DS most likely implies local exposure to antigens, neither this nor previous studies have found a necessary association between DS and pathogens or self-antigens.

Double-stranded RNA-dependent protein kinase (PKR) in antiviral defence in fish and mammals.

The interferon-inducible double-stranded RNA-dependent protein kinase (PKR) is one of the key antiviral arms of the innate immune system. Upon binding of viral double stranded RNA, a viral Pattern Associated Molecular Pattern (PAMP), PKR gets activated and phosphorylates the eukaryotic translation initiation factor 2α (eIF2α) resulting in a protein shut-down that limits viral replication. Since its discovery in the mid-seventies, PKR has been shown to be involved in multiple important cellular processes including apoptosis, proinflammatory and innate immune responses. Viral subversion mechanisms of PKR underline its importance in the antiviral response of the host. PKR activation pathways and its mechanisms of action were previously identified and characterised mostly in mammalian models. However, fish Pkr and fish-specific paralogue Z-DNA-dependent protein kinase (Pkz) also play key role in antiviral defence. This review gives an update on the current knowledge on fish Pkr/Pkz, their conditions of activation and their implication in the immune responses to viruses, in comparison to their mammalian counterparts.

Phylogeny and expression of tetraspanin CD9 paralogues in rainbow trout (Oncorhynchus mykiss).

CD9 is a member of the tetraspanin family, which is characterised by a unique domain structure and conserved motifs. In mammals, CD9 is found in tetraspanin-enriched microdomains (TEMs) on the surface of virtually every cell type. CD9 has a wide variety of roles, including functions within the immune system. Here we show the first in-depth analysis of the cd9 gene family in salmonids, showing that this gene has expanded to six paralogues in three groups (cd9a, cd9b, cd9c) through whole genome duplication events. We suggest that through genome duplications, cd9 has undergone subfunctionalisation in the paralogues and that cd9c1 and cd9c2 in particular are involved in antiviral responses in salmonid fish. We show that these paralogues are significantly upregulated in parallel to classic interferon-stimulated genes (ISGs) active in the antiviral response. Expression analysis of cd9 may therefore become an interesting target to assess teleost responses to viruses.

Interplay between a bacterial pathogen and its host in rainbow trout isogenic lines with contrasted susceptibility to cold water disease.

Infectious diseases are a major constraint on aquaculture. Genetic lines with different susceptibilities to diseases are useful models to identify resistance mechanisms to pathogens and to improve prophylaxis. Bacterial cold-water disease (BCWD) caused by Flavobacterium psychrophilum represents a major threat for freshwater salmonid farming worldwide. A collection of rainbow trout (Oncorhynchus mykiss) isogenic lines was previously produced from a French domestic population. Here, we compared BCWD resistance phenotypes using a subset of isogenic lines chosen for their contrasted susceptibilities to F. psychrophilum. We applied individual monitoring to document the infection process, including time-course quantification of bacteremia and innate immune response. Strikingly, BCWD resistance was correlated with a lower bacterial growth rate in blood. Several immune genes were expressed at higher levels in resistant fish regardless of infection: the Type II arginase (arg2), a marker for M2 macrophages involved in anti-inflammatory responses and tissue repair, and two Toll-like receptors (tlr2/tlr7), responsible for pathogen detection and inflammatory responses. This study highlights the importance of innate and intrinsic defense mechanisms in determining the outcome of F. psychrophilum infections, and illustrates that non-lethal time-course blood sampling for individual monitoring of bacteremia is a powerful tool to resolve within-host pathogen behavior in bacterial fish diseases.

AKT activity orchestrates marginal zone B cell development in mice and humans.

The signals controlling marginal zone (MZ) and follicular (FO) B cell development remain incompletely understood. Here, we show that AKT orchestrates MZ B cell formation in mice and humans. Genetic models that increase AKT signaling in B cells or abolish its impact on FoxO transcription factors highlight the AKT-FoxO axis as an on-off switch for MZ B cell formation in mice. In humans, splenic immunoglobulin (Ig) D+CD27+ B cells, proposed as an MZ B cell equivalent, display higher AKT signaling than naive IgD+CD27- and memory IgD-CD27+ B cells and develop in an AKT-dependent manner from their precursors in vitro, underlining the conservation of this developmental pathway. Consistently, CD148 is identified as a receptor indicative of the level of AKT signaling in B cells, expressed at a higher level in MZ B cells than FO B cells in mice as well as humans.

Toll-like receptors control the accumulation of neutrophils in lymph nodes that expand CD4+ T cells during experimental autoimmune encephalomyelitis.

Toll-like receptors (TLR) control the activation of dendritic cells that prime CD4+ T cells in draining lymph nodes, where these T cells then undergo massive clonal expansion. The mechanisms controlling this clonal T cell expansion are poorly defined. Using the CD4+ T cell-mediated disease experimental autoimmune encephalomyelitis (EAE), we show here that this process is markedly suppressed when TLR9 signaling is increased, without noticeably affecting the transcriptome of primed T cells, indicating a purely quantitative effect on CD4+ T cell expansion. Addressing the underpinning mechanisms revealed that CD4+ T cell expansion was preceded and depended on the accumulation of neutrophils in lymph nodes a few days after immunization. Underlying the importance of this immune regulation pathway, blocking neutrophil accumulation in lymph nodes by treating mice with a TLR9 agonist inhibited EAE progression in mice with defects in regulatory T cells or regulatory B cells, which otherwise developed a severe chronic disease. Collectively, this study demonstrates the key role of neutrophils in the quantitative regulation of antigen-specific CD4+ T cell expansion in lymph nodes, and the counter-regulatory role of TLR signaling in this process.

Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida.

The liver is a multitasking organ with essential functions for vertebrate health spanning metabolism and immunity. In contrast to mammals, our understanding of liver cellular heterogeneity and its role in regulating immunological status remains poorly defined in fishes. Addressing this knowledge gap, we generated a transcriptomic atlas of 47,432 nuclei isolated from the liver of Atlantic salmon (Salmo salar L.) contrasting control fish with those challenged with a pathogenic strain of Aeromonas salmonicida, a problematic bacterial pathogen in global aquaculture. We identified the major liver cell types and their sub-populations, revealing poor conservation of many hepatic cell marker genes utilized in mammals, while identifying novel heterogeneity within the hepatocyte, lymphoid, and myeloid lineages. This included polyploid hepatocytes, multiple T cell populations including γδ T cells, and candidate populations of monocytes/macrophages and dendritic cells. A dominant hepatocyte population radically remodeled its transcriptome following infection to activate the acute phase response and other defense functions, while repressing routine functions such as metabolism. These defense-specialized hepatocytes showed strong activation of genes controlling protein synthesis and secretion, presumably to support the release of acute phase proteins into circulation. The infection response further involved up-regulation of numerous genes in an immune-cell specific manner, reflecting functions in pathogen recognition and killing, antigen presentation, phagocytosis, regulation of inflammation, B cell differentiation and T cell activation. Overall, this study greatly enhances our understanding of the multifaceted role played by liver immune and non-immune cells in host defense and metabolic remodeling following infection and provides many novel cell-specific marker genes to empower future studies of this organ in fishes.

Molecular and Cellular Analysis of the Repair of Zebrafish Optic Tectum Meninges Following Laser Injury.

We studied cell recruitment following optic tectum (OT) injury in zebrafish (Danio rerio), which has a remarkable ability to regenerate many of its organs, including the brain. The OT is the largest dorsal layered structure in the zebrafish brain. In juveniles, it is an ideal structure for imaging and dissection. We investigated the recruited cells within the juvenile OT during regeneration in a Pdgfrß-Gal4:UAS-EGFP line in which pericytes, vascular, circulating, and meningeal cells are labeled, together with neurons and progenitors. We first performed high-resolution confocal microscopy and single-cell RNA-sequencing (scRNAseq) on EGFP-positive cells. We then tested three types of injury with very different outcomes (needle (mean depth in the OT of 200 µm); deep-laser (depth: 100 to 200 µm depth); surface-laser (depth: 0 to 100 µm)). Laser had the additional advantage of better mimicking of ischemic cerebral accidents. No massive recruitment of EGFP-positive cells was observed following laser injury deep in the OT. This type of injury does not perturb the meninx/brain-blood barrier (BBB). We also performed laser injuries at the surface of the OT, which in contrast create a breach in the meninges. Surprisingly, one day after such injury, we observed the migration to the injury site of various EGFP-positive cell types at the surface of the OT. The migrating cells included midline roof cells, which activated the PI3K-AKT pathway; fibroblast-like cells expressing numerous collagen genes and most prominently in 3D imaging; and a large number of arachnoid cells that probably migrate to the injury site through the activation of cilia motility genes, most likely being direct targets of the FOXJ1a gene. This study, combining high-content imaging and scRNAseq in physiological and pathological conditions, sheds light on meninges repair mechanisms in zebrafish that probably also operate in mammalian meninges.

Clonotypic IgH Response against Systemic Viral infection in Pronephros and Spleen of a Teleost Fish.

Upon infection, B lymphocytes develop clonal responses. In teleost fish, which lack lymph nodes, the kinetics and location of B cell responses remain poorly characterized. Fish pronephros is the site of B cell differentiation and the main niche for persistence of plasma cells. In this study, we undertook the analysis of the rainbow trout IgHµ repertoire in this critical tissue for humoral adaptive immunity after primary immunization and boost with a rhabdovirus, the viral hemorrhagic septicemia virus (VHSV). We used a barcoded 5' RACE-cDNA sequencing approach to characterize modifications of the IgHµ repertoire, including VH usage in expressed V(D)J rearrangements, clonal diversity, and clonotype sharing between individual fish and treatments. In the pronephros, our approach quantified the clonotype frequency across the whole IgH repertoire (i.e., with all VH), measuring the frequency of Ag-responding clonotypes. Viral infection led to extensive modifications of the pronephros B cell repertoire, implicating several VH subgroups after primary infection. In contrast, only modest changes in repertoire persisted 5 mo later, including VHSV-specific public expansions. The IgM public response implicating IgHV1-18 and JH5, previously described in spleen, was confirmed in pronephros in all infected fish, strongly correlated to the response. However, the distribution of top clonotypes showed that pronephros and spleen B cells constitute distinct compartments with different IgH repertoires. Unexpectedly, after boost, the frequency of anti-VHSV clonotypes decreased both in pronephros and spleen, raising questions about B cell circulation. A better monitoring of B cell response kinetics in lymphoid tissues will be an essential step to understand B memory and plasmocyte formation mechanisms in fish.

Antigen receptor-engineered Tregs inhibit CNS autoimmunity in cell therapy using nonredundant immune mechanisms in mice.

CD4+ FOXP3+ Tregs are currently explored to develop cell therapies against immune-mediated disorders, with an increasing focus on antigen receptor-engineered Tregs. Deciphering their mode of action is necessary to identify the strengths and limits of this approach. Here, we addressed this issue in an autoimmune disease of the CNS, EAE. Following disease induction, autoreactive Tregs upregulated LAG-3 and CTLA-4 in LNs, while IL-10 and amphiregulin (AREG) were increased in CNS Tregs. Using genetic approaches, we demonstrated that IL-10, CTLA-4, and LAG-3 were nonredundantly required for the protective function of antigen receptor-engineered Tregs against EAE in cell therapy whereas AREG was dispensable. Treg-derived IL-10 and CTLA-4 were both required to suppress acute autoreactive CD4+ T-cell activation, which correlated with disease control. These molecules also affected the accumulation in the recipients of engineered Tregs themselves, underlying complex roles for these molecules. Noteworthy, despite the persistence of the transferred Tregs and their protective effect, autoreactive T cells eventually accumulated in the spleen of treated mice. In conclusion, this study highlights the remarkable power of antigen receptor-engineered Tregs to appropriately provide multiple suppressive factors nonredundantly necessary to prevent autoimmune attacks.

Exploring Zebrafish Larvae as a COVID-19 Model: Probable Abortive SARS-CoV-2 Replication in the Swim Bladder.

Animal models are essential to understanding COVID-19 pathophysiology and for preclinical assessment of drugs and other therapeutic or prophylactic interventions. We explored the small, cheap, and transparent zebrafish larva as a potential host for SARS-CoV-2. Bath exposure, as well as microinjection in the coelom, pericardium, brain ventricle, or bloodstream, resulted in a rapid decrease of SARS-CoV-2 RNA in wild-type larvae. However, when the virus was inoculated in the swim bladder, viral RNA stabilized after 24 h. By immunohistochemistry, epithelial cells containing SARS-CoV-2 nucleoprotein were observed in the swim bladder wall. Our data suggest an abortive infection of the swim bladder. In some animals, several variants of concern were also tested with no evidence of increased infectivity in our model. Low infectivity of SARS-CoV-2 in zebrafish larvae was not due to the host type I interferon response, as comparable viral loads were detected in type I interferon-deficient animals. A mosaic overexpression of human ACE2 was not sufficient to increase SARS-CoV-2 infectivity in zebrafish embryos or in fish cells in vitro. In conclusion, wild-type zebrafish larvae appear mostly non-permissive to SARS-CoV-2, except in the swim bladder, an aerial organ sharing similarities with the mammalian lung.

Interferons and interferon receptors in the channel catfish, Ictalurus punctatus.

In this work, we describe the complete repertoire of channel catfish, Ictalurus punctatus, IFNs and IFN receptor genes. Based on multiple genomic and transcriptomic resources we identified 16 type I IFN genes, which represent the six type I IFN subgroups previously defined in salmonids (a-f.) No representatives of subgroup h previously only found in percomorphs were identified. An expansion in copy numbers of subgroup d IFN genes was of particular interest, as this has not been reported in other fish species to date. Furthermore, we confirmed the presence of two type II ifn genes encoding orthologs of IFNγ and the teleost-specific IFNγRel. Six homologs of IFN type I receptor genes were found in an array that shows conserved synteny with human chromosome 21. Three homologs of type II IFN receptor genes were also identified. These type I and type II receptor sequences are compatible with the dual type I IFN receptors, and the potentially more complex type II IFN receptors described in teleosts. Our data provide a comprehensive resource for future studies of channel catfish innate antiviral immunity.

Intranasal delivery of SARS-CoV-2 spike protein is sufficient to cause olfactory damage, inflammation and olfactory dysfunction in zebrafish.

Anosmia, loss of smell, is a prevalent symptom of SARS-CoV-2 infection. Anosmia may be explained by several mechanisms driven by infection of non-neuronal cells and damage in the nasal epithelium rather than direct infection of olfactory sensory neurons (OSNs). Previously, we showed that viral proteins are sufficient to cause neuroimmune responses in the teleost olfactory organ (OO). We hypothesize that SARS-CoV-2 spike (S) protein is sufficient to cause olfactory damage and olfactory dysfunction. Using an adult zebrafish model, we report that intranasally delivered SARS-CoV-2 S RBD mostly binds to the non-sensory epithelium of the olfactory organ and causes severe olfactory histopathology characterized by loss of cilia, hemorrhages and edema. Electrophysiological recordings reveal impaired olfactory function to both food and bile odorants in animals treated intranasally with SARS-CoV-2 S RBD. However, no loss of behavioral preference for food was detected in SARS-CoV-2 S RBD treated fish. Single cell RNA-Seq of the adult zebrafish olfactory organ indicated widespread loss of olfactory receptor expression and inflammatory responses in sustentacular, endothelial, and myeloid cell clusters along with reduced numbers of Tregs. Combined, our results demonstrate that intranasal SARS-CoV-2 S RBD is sufficient to cause structural and functional damage to the zebrafish olfactory system. These findings may have implications for intranasally delivered vaccines against SARS-CoV-2.

Recurrent expansions of B30.2-associated immune receptor families in fish.

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.