Origin and evolution of the vertebrate leukocyte receptors: the lesson from tunicates

Understanding vertebrate immunity through comparative immunology

Evolutionary immunology has entered a new era. Classical studies, using just a handful of model animal species, combined with clinical observations, provided an outline of how innate and adaptive immunity work together to ensure tissue homeostasis and to coordinate the fight against infections. However, revolutionary advances in cellular and molecular biology, genomics and methods of genetic modification now offer unprecedented opportunities. They provide immunologists with the possibility to consider, at unprecedented scale, the impact of the astounding phenotypic diversity of vertebrates on immune system function. This Perspective is intended to highlight some of the many interesting, but largely unexplored, biological phenomena that are related to immune function among the roughly 60,000 existing vertebrate species. Importantly, hypotheses arising from such wide-ranging comparative studies can be tested in representative and genetically tractable species. The emerging general principles and the discovery of their evolutionarily selected variations may inspire the future development of novel therapeutic strategies for human immune disorders.

Genomic view of the origins of cell-mediated immunity

NKp30 is an activating natural killer cell receptor (NKR) with a single-exon variable (VJ)-type immunoglobulin superfamily (IgSF) domain. Such VJ-IgSF domains predate the emergence of the antigen receptors (immunoglobulin and T cell receptor), which possess the same domain but undergo gene rearrangement. NCR3, the gene encoding NKp30, is present in jawed vertebrates from sharks to mammals; thus, unlike most NKR that are highly divergent among vertebrate taxa, NKp30 is uniquely conserved. We previously hypothesized that an ancestral NCR3 gene was encoded in the proto-major histocompatibility complex (MHC), the region where many immune-related genes have accumulated. Herein, we searched in silico databases to identify NCR3 paralogues and examined their genomic locations. We found a paralogue, NCR3H, in many vertebrates but was lost in mammals. Additionally, we identified a set of voltage-gated sodium channel beta (SCNB) genes as NCR3-distantly-related genes. Like NCR3, both NCR3H and SCNB proteins contain a single VJ-IgSF domain followed by a transmembrane region. These genes map to MHC paralogous regions, originally described in an invertebrate, along with genes encoding cell adhesion molecules involved in NK cell recognition networks. Other genes having no obvious relationship to immunity also map to these paralogous regions. These gene complexes were traced to several invertebrates, suggesting that the foundation of these cellular networks emerged before the genome-wide duplications in early gnathostome history. Here, we propose that this ancestral region was involved in cell-mediated immunity prior to the emergence of adaptive immunity and that NCR3 piggybacked onto this primordial complex, heralding the emergence of vertebrate NK cell/T cells.

An Ancient MHC-Linked Gene Encodes a Nonrearranging Shark Antibody, UrIg, Convergent with IgG

Gnathostome adaptive immunity is defined by the Ag receptors, Igs and TCRs, and the MHC. Cartilaginous fish are the oldest vertebrates with these adaptive hallmarks. We and others have unearthed nonrearranging Ag receptor-like genes in several vertebrates, some of which are encoded in the MHC or in MHC paralogous regions. One of these genes, named UrIg, was detected in the class III region of the shark MHC that encodes a protein with typical V and C domains such as those found in conventional Igs and TCRs. As no transmembrane region was detected in gene models or cDNAs, the protein does not appear to act as a receptor. Unlike some other shark Ig genes, the UrIg V region shows no evidence of RAG-mediated rearrangement, and thus it is likely related to other V genes that predated the invasion of the RAG transposon. The UrIg gene is present in all elasmobranchs and evolves conservatively, unlike Igs and TCRs. Also, unlike Ig/TCR, the gene is not expressed in secondary lymphoid tissues, but mainly in the liver. Recombinant forms of the molecule form disulfide-linked homodimers, which is the form also detected in many shark tissues by Western blotting. mAbs specific for UrIg identify the protein in the extracellular matrix of several shark tissues by immunohistochemistry. We propose that UrIg is related to the V gene invaded by the RAG transposon, consistent with the speculation of emergence of Ig/TCR within the MHC or proto-MHC.

Cytokine Receptor Diversity in the Lamprey Predicts the Minimal Essential Cytokine Networks of Vertebrates

The vertebrate adaptive immune systems (Agnatha and Gnathostomata) use sets of T and B lymphocyte lineages that somatically generate highly diverse repertoires of Ag-specific receptors and Abs. In Gnathostomata, cytokine networks regulate the activation of lymphoid and myeloid cells, whereas little is known about these components in Agnathans. Most gnathostome cytokines are four-helix bundle cytokines with poorly conserved primary sequences. In contrast, sequence conservation across bilaterians has been observed for cognate cytokine receptor chains, allowing their structural classification into two classes, and for downstream JAK/STAT signaling mediators. With conserved numbers among Gnathostomata, human cytokine receptor chains (comprising 34 class I and 12 class II) are able to interact with 28 class I helical cytokines (including most ILs) and 16 class II cytokines (including all IFNs), respectively. Hypothesizing that the arsenal of cytokine receptors and transducers may reflect homologous cytokine networks, we analyzed the lamprey genome and transcriptome to identify genes and transcripts for 23 class I and five class II cytokine receptors alongside one JAK signal mediator and four STAT transcription factors. On the basis of deduction of their respective orthologs, we predict that these receptors may interact with 16 class I and 3 class II helical cytokines (including IL-4, IL-6, IL-7, IL-12, IL-10, IFN-γ, and thymic stromal lymphoprotein homologs). On the basis of their respective activities in mammals, this analysis suggests the existence of lamprey cytokine networks that may regulate myeloid and lymphoid cell differentiation, including potential Th1/Th2 polarization. The predicted networks thus appear remarkably homologous to those of Gnathostomata, albeit reduced to essential functions.

Genomic features of humoral immunity support tolerance model in Egyptian rousette bats

Bats asymptomatically harbor many viruses that can cause severe human diseases. The Egyptian rousette bat (ERB) is the only known reservoir for Marburgviruses and Sosuga virus, making it an exceptional animal model to study antiviral mechanisms in an asymptomatic host. With this goal in mind, we constructed and annotated the immunoglobulin heavy chain locus, finding an expansion on immunoglobulin variable genes associated with protective human antibodies to different viruses. We also annotated two functional and distinct immunoglobulin epsilon genes and four distinctive functional immunoglobulin gamma genes. We described the Fc receptor repertoire in ERBs, including features that may affect activation potential, and discovered the lack of evolutionary conserved short pentraxins. These findings reinforce the hypothesis that a differential threshold of regulation and/or absence of key immune mediators may promote tolerance and decrease inflammation in ERBs.

Retention of duplicated ITAM-containing transmembrane signaling subunits in the tetraploid amphibian species Xenopus laevis

The ITAM-bearing transmembrane signaling subunits (TSS) are indispensable components of activating leukocyte receptor complexes. The TSS-encoding genes map to paralogous chromosomal regions, which are thought to arise from ancient genome tetraploidization(s). To assess a possible role of tetraploidization in the TSS evolution, we studied TSS and other functionally linked genes in the amphibian species Xenopus laevis whose genome was duplicated about 40 MYR ago. We found that X. laevis has retained a duplicated set of sixteen TSS genes, all except one being transcribed. Furthermore, duplicated TCRα loci and genes encoding TSS-coupling protein kinases have also been retained. No clear evidence for functional divergence of the TSS paralogs was obtained from gene expression and sequence analyses. We suggest that the main factor of maintenance of duplicated TSS genes in X. laevis was a protein dosage effect and that this effect might have facilitated the TSS set expansion in early vertebrates.

Immune-directed support of rich microbial communities in the gut has ancient roots

The animal gut serves as a primary location for the complex host-microbe interplay that is essential for homeostasis and may also reflect the types of ancient selective pressures that spawned the emergence of immunity in metazoans. In this review, we present a phylogenetic survey of gut host-microbe interactions and suggest that host defense systems arose not only to protect tissue directly from pathogenic attack but also to actively support growth of specific communities of mutualists. This functional dichotomy resulted in the evolution of immune systems much more tuned for harmonious existence with microbes than previously thought, existing as dynamic but primarily cooperative entities in the present day. We further present the protochordate Ciona intestinalis as a promising model for studying gut host-bacterial dialogue. The taxonomic position, gut physiology and experimental tractability of Ciona offer unique advantages in dissecting host-microbe interplay and can complement studies in other model systems.

Fc receptors for immunoglobulins and their appearance during vertebrate evolution

Receptors interacting with the constant domain of immunoglobulins (Igs) have a number of important functions in vertebrates. They facilitate phagocytosis by opsonization, are key components in antibody-dependent cellular cytotoxicity as well as activating cells to release granules. In mammals, four major types of classical Fc receptors (FcRs) for IgG have been identified, one high-affinity receptor for IgE, one for both IgM and IgA, one for IgM and one for IgA. All of these receptors are related in structure and all of them, except the IgA receptor, are found in primates on chromosome 1, indicating that they originate from a common ancestor by successive gene duplications. The number of Ig isotypes has increased gradually during vertebrate evolution and this increase has likely been accompanied by a similar increase in isotype-specific receptors. To test this hypothesis we have performed a detailed bioinformatics analysis of a panel of vertebrate genomes. The first components to appear are the poly-Ig receptors (PIGRs), receptors similar to the classic FcRs in mammals, so called FcRL receptors, and the FcR γ chain. These molecules are not found in cartilagous fish and may first appear within bony fishes, indicating a major step in Fc receptor evolution at the appearance of bony fish. In contrast, the receptor for IgA is only found in placental mammals, indicating a relatively late appearance. The IgM and IgA/M receptors are first observed in the monotremes, exemplified by the platypus, indicating an appearance during early mammalian evolution. Clearly identifiable classical receptors for IgG and IgE are found only in marsupials and placental mammals, but closely related receptors are found in the platypus, indicating a second major step in Fc receptor evolution during early mammalian evolution, involving the appearance of classical IgG and IgE receptors from FcRL molecules and IgM and IgA/M receptors from PIGR.

Chicken CRTAM binds nectin-like 2 ligand and is upregulated on CD8+ αβ and γδ T lymphocytes with different kinetics

During a search for immunomodulatory receptors in the chicken genome, we identified a previously cloned chicken sequence as CRTAM homologue by its overall identity and several conserved sequence features. For further characterization, we generated a CRTAM specific mab. No staining was detectable in freshly isolated cell preparations from thymus, bursa, caecal tonsils, spleen, blood and intestine. Activation of splenocytes with recombinant IL-2 increased rapid CRTAM expression within a 2 h period on about 30% of the cells. These CRTAM⁺ cells were identified as CD8⁺ γδ T lymphocytes. In contrast, CRTAM expression could not be stimulated on PBL with IL-2, even within a 48 h stimulation period. As a second means of activation, T cell receptor (TCR) crosslinking using an anti-αβ-TCR induced CRTAM on both PBL and splenocytes. While CRTAM expression was again rapidly upregulated on splenocytes within 2 h, it took 48 h to reach maximum levels of CRTAM expression in PBL. Strikingly, albeit the stimulation of splenocytes was performed with anti-αβ-TCR, CRTAM expression after 2 h was mainly restricted to CD8⁺ γδ T lymphocytes, however, the longer anti-TCR stimulation of peripheral blood lymphocytes (PBL) resulted in CRTAM expression on αβ T lymphocytes. In order to characterize the potential ligand we cloned and expressed chicken Necl-2, a member of the nectin and nectin-like family which is highly homologous to its mammalian counterpart. Three independent assays including a reporter assay, staining with a CRTAM-Ig fusion protein and a cell conjugate assay confirmed the interaction of CRTAM with Necl-2 which could also be blocked by a soluble CRTAM-Ig fusion protein or a CRTAM specific mab. These results suggest that chicken CRTAM represents an early activation antigen on CD8⁺ T cells which binds to Necl-2 and is upregulated with distinct kinetics on αβ versus γδ T lymphocytes.

A Basal chordate model for studies of gut microbial immune interactions

Complex symbiotic interactions at the surface of host epithelia govern most encounters between host and microbe. The epithelium of the gut is a physiologically ancient structure that is comprised of a single layer of cells and is thought to possess fully developed immunological capabilities. Ciona intestinalis (sea squirt), which is a descendant of the last common ancestor of all vertebrates, is a potentially valuable model for studying barrier defenses and gut microbial immune interactions. A variety of innate immunological phenomena have been well characterized in Ciona, of which many are active in the gut tissues. Interactions with gut microbiota likely involve surface epithelium, secreted immune molecules including variable region-containing chitin-binding proteins, and hemocytes from a densely populated laminar tissue space. The microbial composition of representative gut luminal contents has been characterized by molecular screening and a potentially relevant, reproducible, dysbiosis can be induced via starvation. The dialog between host and microbe in the gut can be investigated in Ciona against the background of a competent innate immune system and in the absence of the integral elements and processes that are characteristic of vertebrate adaptive immunity.

Evolution of the B7 family: co-evolution of B7H6 and NKp30, identification of a new B7 family member, B7H7, and of B7's historical relationship with the MHC

The B7 family of genes is essential in the regulation of the adaptive immune system. Most B7 family members contain both variable (V)- and constant (C)-type domains of the immunoglobulin superfamily (IgSF). Through in silico screening of the Xenopus genome and subsequent phylogenetic analysis, we found novel genes belonging to the B7 family, one of which is the recently discovered B7H6. Humans and rats have a single B7H6 gene; however, many B7H6 genes were detected in a single large cluster in the Xenopus genome. The B7H6 expression patterns also varied in a species-specific manner. Human B7H6 binds to the activating natural killer receptor, NKp30. While the NKp30 gene is single-copy and maps to the MHC in most vertebrates, many Xenopus NKp30 genes were found in a cluster on a separate chromosome that does not harbor the MHC. Indeed, in all species so far analyzed from sharks to mammals, the number of NKp30 and B7H6 genes correlates well, suggestive of receptor-ligand co-evolution. Furthermore, we identified a Xenopus-specific B7 homolog (B7HXen) and revealed its close linkage to B2M, which we have demonstrated previously to have been originally encoded in the MHC. Thus, our study provides further proof that the B7 precursor was included in the proto MHC. Additionally, the comparative analysis revealed a new B7 family member, B7H7, which was previously designated in the literature as an unknown gene, HHLA2.

Identification of a chicken CLEC-2 homologue, an activating C-type lectin expressed by thrombocytes

Receptors on natural killer (NK) cells are classified as C-type lectins or as Ig-like molecules, and many of them are encoded by two genomic clusters designated natural killer gene complex (NKC) and leukocyte receptor complex, respectively. Here, we describe the analysis of an NKC-encoded chicken C-type lectin, previously annotated as homologue to CD94 and NKG2 and thus designated chicken CD94/NKG2. To further elucidate its potential function on NK cells, we produced a specific mab by immunizing with stably transfected HEK293 cells expressing this lectin. Staining of various chicken tissues revealed minimal reactivity with bursal, or thymus cells. In peripheral blood mononuclear cell and spleen, however, the mab reacted with virtually all thrombocytes, whereas most NK cells in organs such as embryonic spleen, lung and intestine were found to be negative. These findings indicate that the gene may not resemble CD94/NKG2, but rather a CLEC-2 homologue, a claim further supported by sequence features such as an additional extracellular cysteine residue and the presence of a cytoplasmic motif known as a hem immunoreceptor tyrosine-based activation motif, found in C-type lectins such as Dectin-1, CLEC-2, but not CD94/NKG2. The biochemical analyses demonstrated that CLEC-2 is present on the cell surface as heavily glycosylated homodimer, which upon mab crosslinking induced thrombocyte activation, as measured by CD107 expression. These analyses reveal that the chicken NKC may not encode NK cell receptor genes, in particular not CD94 or NKG2 genes, and identifies a chicken CLEC-2 homologue.

A role for variable region-containing chitin-binding proteins (VCBPs) in host gut-bacteria interactions

A number of different classes of molecules function as structural matrices for effecting innate and adaptive immunity. The most extensively characterized mediators of adaptive immunity are the immunoglobulins and T-cell antigen receptors found in jawed vertebrates. In both classes of molecules, unique receptor specificity is effected through somatic variation in the variable (V) structural domain. V region-containing chitin-binding proteins (VCBPs) consist of two tandem Ig V domains as well as a chitin-binding domain. VCBPs are encoded at four loci (i.e., VCBPA-VCBPD) in Ciona, a urochordate, and are expressed by distinct epithelial cells of the stomach and intestine, as well as by granular amoebocytes present in the lamina propria of the gut and in circulating blood. VCBPs are secreted into the gut lumen, and direct binding to bacterial surfaces can be detected by immunogold analysis. Affinity-purified native and recombinant VCBP-C, as well as a construct consisting only of the tandem V domains, enhance bacterial phagocytosis by granular amoebocytes in vitro. Various aspects of VCBP expression and function suggest an early origin for the key elements that are central to the dialogue between the immune system of the host and gut microflora.

The chicken leukocyte receptor cluster

Receptors of the immunoglobulin-like superfamily are critically involved in virtually every aspect of immune responses. One large chromosomal area encoding such immunoregulatory receptors is the leukocyte receptor cluster. Here we review various aspects of the chicken Ig-like receptor (CHIR) family, located on microchromosome 31, an orthologous position to the mammalian leukocyte receptor cluster. The CHIR family has been massively expanded with over hundred CHIR genes that are further distinguished into activating, inhibitory and bifunctional receptors. Comparisons of various features such as amino acid motifs, genomic structure, expression and associated adaptor molecules reveal the homology of CHIR to both the killer Ig-like and the leukocyte Ig-like receptor families, with most pronounced correlation of certain CHIR to the NK cell receptor KIR2DL4. To date the CHIR ligands remain largely obscure with the exception of CHIR-AB1 that binds to chicken IgY. Detailed analyses of CHIR-AB1, its crystal structure, the interaction to IgY and functional capabilities allow us to draw conclusions regarding Fc receptor phylogeny and function.

The invention of lymphocytes

Lamprey and hagfish are surviving representatives of the most ancient vertebrates. They possess adaptive immune systems based on a vast, somatically diversified repertoire of lymphocyte-bound antigen receptors. Despite these similarities to antibody and T cell receptors (TCR) of later vertebrates, the variable lymphocyte receptors (VLR) are not related to the immunoglobulin (Ig)-superfamily of genes; and instead of V(D)J recombination VLR are somatically assembled by a gene conversion process. However, recent studies have revealed two lamprey lymphocyte subsets so closely resembling B cells and T cells that separate lymphocyte lineages must have already existed in the ancestral vertebrate, before Ig/TCR emergence. VLR and Ig/TCR arose independently, but the convergent evolution they display actually reflects their selection in cells with specialized functions.

Hydractinia allodeterminant alr1 resides in an immunoglobulin superfamily-like gene complex

Allorecognition, the ability to discriminate between self and nonself, is ubiquitous among colonial metazoans and widespread among aclonal taxa. Genetic models for the study of allorecognition have been developed in the jawed vertebrates, invertebrate chordate Botryllus, and cnidarian Hydractinia. In Botryllus, two genes contribute to the histocompatibility response, FuHC and fester. In the cnidarian Hydractinia, one of the two known allorecognition loci, alr2, has been isolated, and a second linked locus, alr1, has been mapped to the same chromosomal region, called the allorecognition complex (ARC). Here we isolate alr1 by positional cloning and report it to encode a transmembrane receptor protein with two hypervariable extracellular regions similar to immunoglobulin (Ig)-like domains. Variation in the extracellular domain largely predicts fusibility within and between laboratory strains and wild-type isolates. alr1 was found embedded in a family of immunoglobulin superfamily (IgSF)-like genes, thus establishing that the ARC histocompatibility complex is an invertebrate IgSF-like gene complex.

The red jungle fowl leukocyte receptor complex contains a large, highly diverse number of chicken immunoglobulin-like receptor (CHIR) genes

The chicken Ig-like receptor (CHIR) gene family is located on microchromosome 31, the orthologous region to the mammalian leukocyte receptor complex. CHIR are equally related to the mammalian killer Ig-like receptors and leukocyte Ig-like transcripts, but they occur in a much higher number and diversity. The chicken microchromosome 31 has been neglected in the genome sequence analysis. Here, we provide a first analysis of this region. For this purpose bacterial artificial chromosome (BAC) sequences originating from a single inbred red jungle fowl that served as basis for the chicken genome project were screened for the presence of CHIR sequences and eight BACs were identified as major CHIR containing regions. Since the sequences of these BACs that were available in the database were not complete, sequence gaps were further closed by novel data from the chicken genome project. The entire sequence was aligned into 26 contigs covering 875kbp that contained 84 functional CHIR and 46 CHIR pseudogenes that were hampered by different reasons such as premature stop codons. The 84 functional CHIR were further categorized into 35 activating (CHIRA), 26 inhibitory (CHIRB) and 23 bifunctional (CHIRAB) genes. A detailed comparison of the annotated sequence taking also into account the previously published CHIR BAC sequence originating from an Lohman selected leghorn chicken revealed that the CHIR locus seems to be a very active region with a high degree of gene reorganization that resembles a constant birth and death evolution. The present report provides a framework for the future completion of the CHIR locus. It further suggests that the entire microchromosome 31 may resemble a locus of extraordinary genomic diversity that is beneficial for the development of a large CHIR repertoire, but that has therefore lost all other genes, where such a diversification would be fatal.

Fish 'n' TRIMs

A novel diversified multigene family of tripartite-motif (TRIM) intracellular receptors with putative antiviral activity has been identified in teleost fish and published in BMC Biology. The history of these receptors involves ancient linkage to paralogs of the major histocompatibility complex, and the family has invertebrate precursors.