- Research article
- Open Access
The medaka novel immune-type receptor (NITR) gene clusters reveal an extraordinary degree of divergence in variable domains
BMC Evolutionary Biology volume 8, Article number: 177 (2008)
Novel immune-type receptor (NITR) genes are members of diversified multigene families that are found in bony fish and encode type I transmembrane proteins containing one or two extracellular immunoglobulin (Ig) domains. The majority of NITRs can be classified as inhibitory receptors that possess cytoplasmic immunoreceptor tyrosine-based inhibition motifs (ITIMs). A much smaller number of NITRs can be classified as activating receptors by the lack of cytoplasmic ITIMs and presence of a positively charged residue within their transmembrane domain, which permits partnering with an activating adaptor protein.
Forty-four NITR genes in medaka (Oryzias latipes) are located in three gene clusters on chromosomes 10, 18 and 21 and can be organized into 24 families including inhibitory and activating forms. The particularly large dataset acquired in medaka makes direct comparison possible to another complete dataset acquired in zebrafish in which NITRs are localized in two clusters on different chromosomes. The two largest medaka NITR gene clusters share conserved synteny with the two zebrafish NITR gene clusters. Shared synteny between NITRs and CD8A/CD8B is limited but consistent with a potential common ancestry.
Comprehensive phylogenetic analyses between the complete datasets of NITRs from medaka and zebrafish indicate multiple species-specific expansions of different families of NITRs. The patterns of sequence variation among gene family members are consistent with recent birth-and-death events. Similar effects have been observed with mammalian immunoglobulin (Ig), T cell antigen receptor (TCR) and killer cell immunoglobulin-like receptor (KIR) genes. NITRs likely diverged along an independent pathway from that of the somatically rearranging antigen binding receptors but have undergone parallel evolution of V family diversity.
The two major forms of novel immune-type receptors (NITRs) possess either one N-terminal ectodomain of the Ig variable (V) type, which may contain a joining (J) or J-like sequence or an N-terminal V domain plus a second C-terminal Ig domain of the intermediate (I) type. Most NITRs are predicted to encode type I transmembrane cell surface proteins that possess extracellular immunoglobulin (Ig) domains and can be classified as inhibitory or activating. A smaller number of NITR genes lack a transmembrane domain and are predicted to be secreted proteins [1, 2]. NITRs are expressed in several different hematopoietic lineages [3, 4] and their overall structure and cytoplasmic signaling resemble the mammalian natural killer cell receptors (e.g. KIRs) [1, 5, 6]. NITRs have been reported to bind allogeneic cell surface targets, although their precise ligand specificity is not clear [Cannon et al., unpublished observations and ]. Other observations, including upregulation of transcription in response to pathogens, support a role for NITRs in immunity .
NITRs were identified initially in the compact genome of the Southern pufferfish (Spheroides nephelus); more than 26 NITR genes were shown to be encoded in a single gene cluster [9, 10]. Subsequent analysis of the zebrafish (Danio rerio) genome revealed one large NITR gene cluster encoding 36 NITRs on chromosome 7 and a second smaller NITR gene cluster encoding 3 NITRs on chromosome 14 [2, 11]. NITR cDNAs also have been characterized from channel catfish (Ictalurus punctatus), rainbow trout (Oncorhynchus mykiss) and Japanese flounder (Paralichthys olivaceus) [3, 4, 12, 13] and it is assumed that they are present in all bony fish. The majority of NITRs possess cytoplasmic immunoreceptor tyrosine-based inhibition motifs (ITIMs); inhibitory functions have been confirmed in vitro . A smaller number of NITRs possess a positively charged residue within their transmembrane domain and have been shown to interact with Dap12 through which they mediate an activating function .
We herein report the identification of 44 medaka (Oryzias latipes) NITR genes, which can be classified into 24 families including inhibitory and activating forms; this represents nearly twice the number of NITR families that have been identified in any other species. The resolution of this dataset of NITRs along with the previous genomic resolution of NITRs in zebrafish permits a comprehensive phylogenetic analysis as well as the examination of shared synteny. The NITRs emerge as independent families of immune-type receptors in which the evolution of the V regions exhibit considerable similarity to that seen in multigene families encoding Ig, TCRs and KIRs.
Medaka encode 44 NITRs on three chromosomes
The reduced levels of protein sequence identity between NITRs found in different species make it difficult to identify genetic orthologs between fish species outside of conserved core residues. In order to identify NITR genes from medaka, we data mined the recently released medaka genome database [16, 17]. By using zebrafish NITR sequences as queries, it was possible to identify a number of candidate genes that encode Ig domains. By then restricting the queries only to genes that encode an authentic NITR I-type domain , NITR gene clusters were identified on medaka chromosomes 10, 18 and 21. The search was then directed to include all genes predicted to encode V domains in these genomic regions. In total, 10 NITR genes were identified on chromosome 10; 30 NITR genes and three pseudogenes (defined as possessing a stop codon within the V or I domain) were identified on chromosome 18; and one NITR gene and one pseudogene were identified on chromosome 21 (Figure 1 and Additional File 1). Three additional NITR genes (NITR8b, NITR10b and NITR12b) cannot be assigned to a specific chromosome. It is possible that certain genes described here may represent allelic variants and that additional NITR genes may be identified in future releases of the medaka genome.
The 44 NITR genes that were identified can be grouped into 24 V families, based on >70% predicted amino acid sequence identity (Figure 2a and Additional File 2). Twelve families consist of multiple (up to five) members and the remaining families are singletons. Each NITR family was assigned a number and for families with multiple members, each member was assigned a letter (e.g. NITR1a, NITR1b, etc.) as described . With the exception of the NITR9 family (see below) in which activating potential (rather than V structure) served as the basis for assignment to a specific gene family, unequivocal orthologs between the NITRs of medaka and any other species were not identified and the numbering of these genes was arbitrary.
Two major classes of NITR gene organization
All of the NITR genes currently described from the Southern pufferfish and 11 NITR genes from a second pufferfish species, Tetraodon nigroviridis (Additional file 3), encode two Ig domains (V and I domains) in a single exon . In contrast, all of the two-Ig domain NITR genes described from zebrafish encode the V and I domains in two exons . Both types of NITR gene organization are present in medaka (Figure 2). The medaka NITR genes encoding two Ig domains in a single exon are all located on chromosomes 18 and 21; the medaka NITR genes encoding two Ig domains in two exons are located on chromosome 10 with the exception of the NITR9 genes, which are on chromosome 18.
Twenty-one full-length NITR cDNAs were cloned from medaka kidney and spleen by 3' RACE using nested primers that include the predicted translational start codon (Additional Files 4 and 5) in order to verify the exon-intron boundaries for the NITRs. The predicted exon-intron boundaries most often were verified; however, in several cases, predicted mRNA splice sites were inexact, reflected slight differences in cDNA sequence. It is likely that the genome from inbred Hd-rR medaka and cDNAs from the out-bred orange-red medaka line possess different alleles for certain NITRs. For example, the identification of a NITR6 cDNA from the orange-red line confirms that the NITR6 genes possess a cytoplasmic tail; however, the cDNA cytoplasmic sequence is too divergent to confidently predict these exons from the Hd-rR genome. In addition the orange-red NITR12 cDNA encodes additional nucleotide sequences that are not present in the Hd-rR genome resulting in a stop codon preventing translation of an ITIM-like sequence (itim) in the cDNA (Additional File 6).
Structural features of medaka NITRs
The initial search criteria for identifying medaka NITRs was directed to the I domain, which possesses six highly conserved cysteines. NITR22 and NITR23 lack the sixth and third cysteines, respectively (Figure 3). A joining (J; FGXGTXLXV) or J-like peptide sequence at the carboxyl-terminus of either the V and/or the I domain is a conserved feature of nearly all medaka NITRs (Figures 3 and 4 and Additional File 7). The majority of NITRs in medaka and other species can be classified as inhibitory based on the presence of one or two cytoplasmic ITIM or itim sequences (Figure 4 and Additional File 8) as well as by the absence of a charged transmembrane residue, a conserved feature of activating receptors. In medaka, three closely related NITR genes (NITR9a, NITR9b and NITR9c) possess a charged residue (Arg) within their transmembrane domain corresponding to activating NITRs in other species such as Nitr9, the single activating NITR gene in zebrafish (Additional File 9) [2, 15].
NITR gene clusters in medaka and zebrafish share conserved synteny
The complete annotation of NITRs in zebrafish and medaka permits the investigation of potential shared synteny. The larger zebrafish NITR gene cluster on chromosome 7 is adjacent to the map4k2l, men1 and rtn3 genes  and the larger medaka NITR gene cluster is linked to these genes on chromosome 18 (Table 1). The second smaller zebrafish NITR gene cluster on chromosome 14 is adjacent to the fgfrl1b gene  and the medaka NITR gene cluster on chromosome 10 is tightly linked to FGFRL1(2of2) within the same contig (Figure 5). (Medaka and zebrafish possess two FGFRL1 orthologs referred to as fgfrl1a and fgfrl1b in zebrafish and fgfrl1(1of2) and fgfrl1(2of2) in medaka). Medaka orthologs of additional genes adjacent to the NITR gene cluster on zebrafish chromosome 14 are also encoded on medaka chromosome 10 (Table 1). These observations indicate that the NITR gene clusters on medaka chromosome 18 and zebrafish chromosome 7 share conserved synteny whereas the NITR gene cluster on medaka chromosome 10 and zebrafish chromosome 14 share conserved synteny. The origin of the third, smaller NITR gene cluster on medaka chromosome 21 is not evident, as zebrafish only possess two NITR gene clusters.
A bony fish SPON2/FGFRL1/NITR locus shares conserved synteny with the SPON2/FGFRL1/CD8 locus in Xenopusand chicken
NITRs have not been identified outside of bony fish and their evolutionary origins remain unresolved. The genomes of human, chicken, Xenopus tropicalis and a third fish species stickleback (Gasterosteus aculeatus) were queried in an effort to identify genomic regions, which share conserved synteny with the NITR gene clusters. In medaka, zebrafish and stickleback, one NITR gene cluster is linked to SPON2 and FGFRL1 orthologs. In chicken and Xenopus the SPON2 and FGFRL1 genes are physically linked to the CD8A/CD8B locus (Figure 5). SPON2/FGFRL1 and CD8A/CD8B genes are on different chromosomes in human. In medaka, zebrafish and stickleback the CD8A/CD8B locus does not display conserved synteny to human (Additional File 10). The observation that SPON2 and FGFRL are linked to the CD8 locus in chicken and Xenopus and that SPON2 and FGFRL orthologs are linked to a NITR gene cluster in medaka, zebrafish and stickleback is consistent with a common ancestry for these immune-type molecules.
Medaka NITR24 is tightly linked to a CTLA4ortholog
The NITR24 locus on medaka chromosome 21 is separated from the V-domain receptor CTLA4 by ~0.1 Mb. A small number of predicted olfactory receptor genes lie between CTLA4 and NITR24 (Additional File 11). Although mammalian CTLA4 is tightly linked to CD28, a structurally similar receptor, these genes do not appear to be closely linked in bony fish . Neither NITR gene cluster in zebrafish is linked to CTLA4 or CD28.
Species-specific expansions of NITR families
In general, most NITRs from a single species lack clear orthologs in other species. The absence of such comparison sets skews multi-species phylogenetic comparisons as not all NITR members from each species (and family) have been identified. Three major inferences can be made from the data described here (Additional File 12) and other observations predicted previously: 1) the activating receptors, medaka NITR9, zebrafish Nitr9, and catfish IpNITR2, group together suggesting a common ancestry; 2) the majority of the NITRs encoding both V and I domains within a single exon (medaka, Southern pufferfish and Tetraodon) group together, suggesting that they are likely derived from a common ancestral gene; and 3) the majority of the NITRs, which lack an I domain (from medaka, zebrafish, catfish, and Tetraodon), group together, suggesting that they too may be derived from a common ancestral gene.
Multiple species-specific expansions of individual NITR families are revealed by comparing the two complete NITR datasets from medaka and zebrafish (Figure 6). For example, the medaka NITR gene cluster on chromosome 18 includes expansions of the NITR1, NITR2, NITR3 and NITR5 families and the corresponding zebrafish NITR gene cluster on chromosome 7 includes expansions of the divergent nitr1, nitr2 and nitr3 families. Three NITRs are encoded by the zebrafish NITR gene cluster on chromosome 14 whereas 10 NITRs are encoded by the corresponding gene cluster in medaka suggesting either gene loss in zebrafish or gene duplications in medaka.
Interpreting the evolution of NITRs within teleost species
Forty-four different NITR genes, the majority of which are inhibitory, are present in the medaka genome. Perhaps the most striking observation relates to the high degree of variation between these V domains within medaka and their diversification from NITRs in other fish species. This finding is particularly striking given the exceptionally large number of V families in medaka, i.e., the greater the number of V families present, the higher the expectations of identifying orthologous sequences. The medaka/zebrafish chromosomes 18/7 and 10/14, which encode the larger NITR gene clusters, exhibit a high degree of conserved synteny . Medaka chromosomes 10 and 18, and zebrafish chromosomes 7 and 14 likely derive from a hypothetical ancestral chromosome termed "f" which was duplicated prior to the divergence of medaka and zebrafish . Although the nucleotide identity is low between medaka and zebrafish NITR genes (confounding the confident assignment of genetic orthologs), conserved synteny between the two primary medaka and zebrafish NITR gene clusters is impressive.
Comparisons of predicted amino acid sequence structures within and between different species indicate that NITRs are undergoing family-specific expansions and individual genes are undergoing birth and death events. As the NITR gene products are highly diverse within a single species (Figures 3 and 4) and approximately 10% of the NITRs described in medaka can be classified as pseudogenes (Figure 1), it is less likely that concerted evolution played a primary role in establishing the current NITR sequences . Birth and death events of NITR genes likely have been occurring continuously after the divergence of medaka and zebrafish ~323 MYA . Similar observations have been made concerning the V gene segments of the Ig and TCR genes  and the KIR genes within primate species [reviewed by ]. The KIR gene family is of particular value in comparative interpretations as they also represent type I, Ig domain-containing membrane proteins with potential for inhibitory and activating signalling. KIRs are highly polymorphic and polygenic within single species and drastically diverged in gene sequence and number between different primates . The nature of diversity in KIRs is markedly similar to that observed in the NITRs of medaka and zebrafish. The studies reported here suggest that particularly extensive expansion and diversification of the NITRs have occurred in bony fish. Deducing the status of the NITR gene family from the genomes of more basal fish species (such as gar or sturgeon) will be highly informative regarding the origins of diversity within the NITR family.
Possible links between NITRs and mammalian receptors
All evidence to date strongly suggests that NITRs are restricted to bony fish. Although NITRs possess authentic V (and frequently V-J) domains reminiscent of TCR, their overall structure, genomic organization and the presence of both activating and inhibitory forms, has led to the hypothesis that they may be functional equivalents to mammalian NK receptors, resembling in some way KIRs. We previously reported that the zebrafish NITR gene cluster on chromosome 7 shares a limited syntenic relationship with the mammalian leukocyte receptor cluster (LRC) that includes the KIRs . The new data presented here introduces two alternative models in which the NITRs share a common ancestry with CD8A/CD8B or with CD28/CTLA4. An expansion of CD8A-like genes was recently described within the chicken CD8A/CD8B locus; the reported structural similarity between the V-like domains of the CD8A-like receptors and NITRs supports the former model . Although intriguing, the latter model is not as well supported. Specifically, in mammals CD28 and CTLA4 effect opposing roles in T-cell stimulation and are encoded by adjacent genes. However, the CD28 and CTLA4 genes are not physically linked in teleosts. Linkage of an NITR gene to CTLA4 has not been identified outside of medaka; in zebrafish the NITR gene clusters are on chromosomes 7 and 14 [2, 11]. CTLA4 and CD28 genes in zebrafish are located on chromosomes 9 and 1, respectively .
KIRs, CD8, CD28 and CTLA4 are expressed by cytotoxic lymphocytes and NITRs are similarly expressed in lymphoid populations including NK-like cytotoxic cells [unpublished observations and [3, 4]]. In general, inhibitory natural killer cells receptors (NKRs) repress cytotoxicity by recognition of MHC or other markers of "self" on candidate target cells whereas activating NKRs induce cytotoxicity by detecting a loss of "self" or a gain of "non-self" (e.g. viral or stress-related peptides). In at least one case, a NITR specifically engages a non-self determinant (Cannon et al., unpublished observations), suggesting the possibility that they may be involved in cytotoxic interactions. Ongoing efforts to identify NITR ligands and define cell populations that express NITRs in vivo will prove crucial in deciphering their role in immunity. In particular, such studies will determine whether their function is analogous to mammalian processes (NK recognition) or represents a unique adaptation within bony fish.
The medaka genome encodes 44 NITR genes in 3 gene clusters. NITR genes in this species are divergent in sequence but structurally similar to those characterized in other fish species and include both inhibitory and activating forms. Comparative genomics reveals three possible scenarios in which the NITRs share a common genetic ancestry with CD8, CD28/CTLA4 or KIRs. A comparison of the complete sets of NITRs from medaka and zebrafish indicate that: 1) the two largest NITR gene clusters in medaka and zebrafish are derived from a whole chromosome duplication event, 2) as a gene family, the NITRs have experienced recent, species-specific "birth and death" events and 3) NITR family expansions have led to high levels of sequence variation between fish species.
Using zebrafish nitr1c and nitr3a as queries, putative medaka NITR genes were identified by searching version 1.0 of the medaka genome . Manual annotation of these medaka genes was performed using FSPLICE to predict splice sites (GT/AG) and FGENESH+ to predict structures based on homology with zebrafish and other medaka NITR genes: "fish" specific parameters were applied to both programs which are included in the MolQuest software suite . SEQtools was used to visualize the predicted transcripts and their three-frame translations . In addition, V and I domain predictions were validated by searching the NCBI non-redundant database with BLAST , while transmembrane predictions were validated using the TMHMM 2.0 server . Lastly, groups of medaka NITRs were defined through BLASTCLUST (single-linkage) based on sharing 70% identity within the peptide sequence of the V-domain. Orthologous genes and their genomic positions were identified in medaka, zebrafish, stickleback, Xenopus tropicalis, chicken, mouse and human genomes using ContigView software within the Ensembl genome browser .
Published NITR V or I domains were aligned by ClustalW  and neighbor-joining trees  were constructed from pairwise Poisson correction distances with 2000 bootstrap replications by MEGA2.1 software . NITR sequences described from medaka, zebrafish, Southern pufferfish, channel catfish, rainbow trout, Japanese flounder and Tetraodon (Additional file 3) were included in alignments [2–4, 10–14].
Rapid Amplification of cDNA ends
Tissues from young adult (6–9 month), out-bred, orange-red medaka were dissected directly into RNA-Bee and total RNA purified as recommended by the manufacturer (Tel-Test, Friendswood, TX). Two μg of RNA (1 μg from kidney and 1 μg from spleen) were reverse transcribed using the oligo dT primer and Superscript™ III Reverse Transcriptase supplied with the GeneRacer™ kit (Invitrogen, Carlsbad, CA). Nested, forward primers were designed against the exons encoding predicted leader sequences of medaka NITRs and 3' RACE was completed as recommended in the GeneRacer kit. Primer sequences are listed in Additional File 4. PCR products were cloned into pCRII-TOPO (Invitrogen) or pGEM-T Easy (Promega, Madison, WI) and sequenced.
Sequence data described in this article have been deposited with the GenBank database under accession numbers: NITR1, transcript variant 1 [GenBank:EU419354]; NITR1, transcript variant 2 [GenBank:EU419355]; NITR1, transcript variant 3 [GenBank:EU419356]; NITR2, transcript variant 1 [GenBank:EU419357]; NITR3, transcript variant 1 [GenBank:EU419358]; NITR4, transcript variant 1 [GenBank:EU419359]; NITR5, transcript variant 1 [GenBank:EU419360]; NITR5, transcript variant 2 [GenBank:EU419361]; NITR6, transcript variant 1 [GenBank:EU419362]; NITR8, transcript variant 1 [GenBank:EU419363]; NITR9, transcript variant 1 [GenBank:EU419364]; NITR11, transcript variant 1 [GenBank:EU419365]; NITR12, transcript variant 1 [GenBank:EU419366]; NITR14, transcript variant 1 [GenBank:EU419367]; NITR17, transcript variant 1 [GenBank:EU419368]; NITR18, transcript variant 1 [GenBank:EU419369]; NITR20, transcript variant 1 [GenBank:EU419370]; NITR22, transcript variant 1 [GenBank:EU419371]; NITR24, transcript variant 1 [GenBank:EU419372]; NITR24, transcript variant 2 [GenBank:EU419373]; and NITR24, transcript variant 3 [GenBank:EU419374].
Litman GW, Hawke NA, Yoder JA: Novel immune-type receptor genes. Immunol Rev. 2001, 181: 250-259. 10.1034/j.1600-065X.2001.1810121.x.
Yoder JA, Litman RT, Mueller MG, Desai S, Dobrinski KP, Montgomery JS, Buzzeo MP, Ota T, Amemiya CT, Trede NS, Wei S, Djeu JY, Humphray S, Jekosch K, Hernandez Prada JA, Ostrov DA, Litman GW: Resolution of the novel immune-type receptor gene cluster in zebrafish. Proc Natl Acad Sci U S A. 2004, 101: 15706-15711. 10.1073/pnas.0405242101.
Evenhuis J, Bengten E, Snell C, Quiniou SM, Miller NW, Wilson M: Characterization of additional novel immune type receptors in channel catfish, Ictalurus punctatus. Immunogenetics. 2007
Hawke NA, Yoder JA, Haire RN, Mueller MG, Litman RT, Miracle AL, Stuge T, Shen L, Miller N, Litman GW: Extraordinary variation in a diversified family of immune-type receptor genes. Proc Natl Acad Sci U S A. 2001, 98: 13832-13837. 10.1073/pnas.231418598.
Sambrook JG, Beck S: Evolutionary vignettes of natural killer cell receptors. Curr Opin Immunol. 2007, 19: 553-560. 10.1016/j.coi.2007.08.002.
Yoder JA: Investigating the morphology, function and genetics of cytotoxic cells in bony fish. Comp Biochem Physiol C Toxicol Pharmacol. 2004, 138 (3): 271-280. 10.1016/j.cca.2004.03.008.
Litman GW, Cannon JP, Dishaw LJ: Reconstructing immune phylogeny: new perspectives. Nat Rev Immunol. 2005, 5: 866-879. 10.1038/nri1712.
Rojo I, de Ilarduya OM, Estonba A, Pardo MA: Innate immune gene expression in individual zebrafish after Listonella anguillarum inoculation. Fish Shellfish Immunol. 2007, 23: 1285-1293. 10.1016/j.fsi.2007.07.002.
Rast JP, Haire RN, Litman RT, Pross S, Litman GW: Identification and characterization of T-cell antigen receptor related genes in phylogenetically diverse vertebrate species. Immunogenetics. 1995, 42: 204-212. 10.1007/BF00191226.
Strong SJ, Mueller MG, Litman RT, Hawke NA, Haire RN, Miracle AL, Rast JP, Amemiya CT, Litman GW: A novel multigene family encodes diversified variable regions. Proc Natl Acad Sci U S A. 1999, 96: 15080-15085. 10.1073/pnas.96.26.15080.
Yoder JA, Cannon JP, Litman RT, Murphy C, Freeman JL, Litman GW: Evidence for a transposition event in a second NITR gene cluster in zebrafish. Immunogenetics. 2008, 60: 257-265. 10.1007/s00251-008-0285-3.
Piyaviriyakul P, Kondo H, Hirono I, Aoki T: A novel immune-type receptor of Japanese flounder (Paralichthys olivaceus) is expressed in both T and B lymphocytes. Fish Shellfish Immunol. 2007, 22: 467-476. 10.1016/j.fsi.2006.05.007.
Yoder JA, Mueller MG, Nichols KM, Ristow SS, Thorgaard GH, Ota T, Litman GW: Cloning novel immune-type inhibitory receptors from the rainbow trout, Oncorhynchus mykiss. Immunogenetics. 2002, 54: 662-670. 10.1007/s00251-002-0511-3.
Yoder JA, Mueller MG, Wei S, Corliss BC, Prather DM, Willis T, Litman RT, Djeu JY, Litman GW: Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian lymphocyte receptor cluster. Proc Natl Acad Sci USA. 2001, 98: 6771-6776. 10.1073/pnas.121101598.
Wei S, Zhuo JM, Chen X, Shah R, Liu J, Orcutt TM, Traver D, Djeu JY, Litman GW, Yoder JA: The zebrafish activating immune receptor Nitr9 signals via Dap12. Immunogenetics. 2007, 59: 813-821. 10.1007/s00251-007-0250-6.
Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, Jindo T, Kobayashi D, Shimada A, Toyoda A, Kuroki Y, Fujiyama A, Sasaki T, Shimizu A, Asakawa S, Shimizu N, Hashimoto S, Yang J, Lee Y, Matsushima K, Sugano S, Sakaizumi M, Narita T, Ohishi K, Haga S, Ohta F, Nomoto H, Nogata K, Morishita T, Endo T, Shin I, Takeda H, Morishita S, Kohara Y: The medaka draft genome and insights into vertebrate genome evolution. Nature. 2007, 447: 714-719. 10.1038/nature05846.
UTGB - Medaka. [http://medaka.utgenome.org/]
Bernard D, Hansen JD, Du PL, Lefranc MP, Benmansour A, Boudinot P: Costimulatory receptors in jawed vertebrates: conserved CD28, odd CTLA4 and multiple BTLAs. Dev Comp Immunol. 2007, 31: 255-271. 10.1016/j.dci.2006.06.003.
Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.
Ota T, Sitnikova T, Nei M: Evolution of vertebrate immunoglobulin variable gene segments. Curr Top Microbiol Immunol. 2000, 248: 221-245.
Martinez-Borra J, Khakoo SI: Speed and selection in the evolution of killer-cell immunoglobulin-like receptors. Int J Immunogenet. 2008, 35: 89-96. 10.1111/j.1744-313X.2008.00756.x.
Liaw HJ, Chen WR, Huang YC, Tsai CW, Chang KC, Kuo CL: Genomic organization of the chicken CD8 locus reveals a novel family of immunoreceptor genes. J Immunol. 2007, 178: 3023-3030.
ENSEMBL - Medaka Genome. [http://www.ensembl.org/Oryzias_latipes/]
MolQuest - Bioinformatics Toolbox. [http://www.molquest.com/molquest.phtml]
SEQtools - a program suite for sequence analysis. [http://www.seqtools.dk/]
Ye J, McGinnis S, Madden TL: BLAST: improvements for better sequence analysis. Nucleic Acids Res. 2006, 34: W6-W9. 10.1093/nar/gkl164.
Center for Biological Sequence Analysis - TMHMM Server. [http://www.cbs.dtu.dk/services/TMHMM/]
Flicek P, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Eyre T, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Holland R, Howe KL, Howe K, Johnson N, Jenkinson A, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Megy K, Meidl P, Overduin B, Parker A, Pritchard B, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Slater G, Smedley D, Spudich G, Trevanion S, Vilella AJ, Vogel J, White S, Wood M, Birney E, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Herrero J, Hubbard TJ, Kasprzyk A, Proctor G, Smith J, Ureta-Vidal A, Searle S: Ensembl 2008. Nucleic Acids Res. 2008, 36: D707-D714. 10.1093/nar/gkm988.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.
Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: molecular evolutionary genetics analysis software. Bioinformatics. 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244.
Giudicelli V, Duroux P, Ginestoux C, Folch G, Jabado-Michaloud J, Chaume D, Lefranc MP: IMGT/LIGM-DB, the IMGT comprehensive database of immunoglobulin and T cell receptor nucleotide sequences. Nucleic Acids Res. 2006, 34: D781-D784. 10.1093/nar/gkj088.
Ravetch JV, Lanier LL: Immune inhibitory receptors. Science. 2000, 290: 84-89. 10.1126/science.290.5489.84.
Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks. Nucleic Acids Res. 2006, 34: D257-D260. 10.1093/nar/gkj079.
Cossins AR, Crawford DL: Fish as models for environmental genomics. Nat Rev Genet. 2005, 6: 324-333. 10.1038/nrg1590.
We are indebted to Dave Hinton, Seth Kullman and Sheran Law for medaka tissues. We thank Chris Amemiya for very helpful discussions and advice on phylogenetic relationships, Jacques Robert and Alexander Taranin for sharing the chromosomal location of Xenopus tropicalis CD8B and Barbara Pryor for editorial assistance. This work was supported by grant R01 AI057559 from NIH.
SD completed the gene prediction analyses. AKH, TMO and JAY cloned and sequenced the medaka NITR cDNAs. JAY conducted the phylogenetic analyses and comparative genomics. GWL and JAY interpreted data and wrote the manuscript. All authors have read and approved the final version of the manuscript.
Electronic supplementary material
Additional File 1: Predicted nucleotide sequence of NITRs based on the 1st draft of the medaka Hd-rR genome. Sequences are in FASTA format and include only full-length open reading frames. Predicted exons are color-coded (leader = green; V domain = light blue; I domain = red; transmembrane domain = pink; cytoplasmic exons are dark blue and yellow). (DOC 76 KB)
Additional File 2: Protein sequences encoded by predicted NITR transcripts based on the 1st draft of the medaka Hd-rR genome. Sequences are in FASTA format and include translated sequences from Additional file 1. Predicted protein domains are color-coded (leader = green; V domain = light blue; I domain = red; transmembrane domain = pink; cytoplasmic domain = yellow). (DOC 55 KB)
Additional file 3: Tetraodon nigroviridis encode multiple NITRs. Putative NITR genes and gene transcripts were identified from Tetraodon nigroviridis by a tBLASTn search using a Spheroides nephelus NITR1 protein sequence [GenBank:AAF28781] as a query. This strategy identified six candidate Tetraodon cDNAs [GenBank:CR683680, GenBank:CR675622, GenBank:CR664265, GenBank:CR725053, GenBank:CR724466, and GenBank:CR707240] and a single genomic BAC clone B22K12 [GenBank:BX629355]. Only one of these cDNAs [GenBank:CR683680] appears to include an entire ORF. The partial protein sequences encoded by [GenBank:CR675622] and [GenBank:CR664265] are identical. The partial protein sequences encoded by [GenBank:CR725053] and [GenBank:CR724466] differ at three predicted amino acids residues; [GenBank:CR725053] is 10 amino acids longer at the amino-terminus. BAC B22K12 encodes 12 different NITR genes (designated by their nucleotide location in the reverse complement position within BAC B22K12). The predicted protein sequences of four representative NITR cDNAs are aligned relative to the predicted protein sequences of the NITRs encoded by BAC B22K12. Note that NITR-116851 (BAC B22K12) is predicted to encode a V domain and lacks an I domain. (PDF 89 KB)
Additional File 4: Oligonucleotide primer sequences used for nested 3' RACE. Predicted translational start sites are underlined. Primers designed to amplify NITR10 amplified NITR24 transcripts. Primers designed to amplify NITR13 amplified NITR4 transcripts. Primers designed to amplify NITR16 amplified NITR1 transcripts. (PDF 112 KB)
Additional File 5: Classification of medaka NITR cDNAs. Neighbor-joining tree of V domains encoded by out-bred, orange-red medaka NITR cDNAs (highlighted in yellow) and predicted NITR V domains from the genome of the inbred Hd-rR medaka. The number assigned to each interior branch corresponds to the bootstrap value. Branch lengths correspond to the number of amino acid substitutions estimated by Poisson correction, scale indicated below. (PDF 23 KB)
Additional File 6: Variation between the medaka NITR12 cDNA and genomic sequence. A partial alignment of NITR12 cDNA from out-bred, orange-red medaka to the predicted coding sequence of the NITR12a gene from inbred Hd-rR medaka. Predicted peptide sequence is shown above (cDNA) and below (genomic DNA) the alignment. The sequences shown include a portion of the transmembrane domain exon and the two cytoplasmic domain exons. Dashes indicate gaps in the nucleotide alignment. Asterisks indicate nucleotide identity. Stop codons are indicated by red highlighting in the nucleotide sequence and by dashes in the protein sequence. Downward arrowheads define predicted RNA splice sites. The NITR12 cDNA encodes additional nucleotide sequences that are not present in the genomic sequence (nucleotides highlighted in pink). Although the NITR12a genomic sequence is predicted to encode an ITIM (highlighted in yellow) and a second ITIM-like sequence (itim; highlighted in orange), differences in the NITR12 cDNA alter the predicted ITIM sequence (LNYAAL to RYYAAL) and introduce a stop codon, preventing translation of the itim. The differences in predicted protein sequence likely reflect allelic variations between the out-bred, orange-red transcripts and Hd-rR genomic sequence. (PDF 138 KB)
Additional File 7: Joining (J) and J-like sequences in medaka NITRs. The consensus joining (J) sequence is FGXGTXLX(V/L). A J-like sequence possesses either motif FGXG or TXLX(V/L/I) or partial sequence of both motifs. Sequences were compiled from version 1.0 of the Hd-rR genomic sequence. (PDF 180 KB)
Additional File 10: The CD8 locus does not display conserved synteny between medaka, zebrafish, stickleback and human. Although the CD8A and CD8B genes are tightly linked in (a) medaka, (b) zebrafish, (c) stickleback and (d) human, flanking genes are not well conserved between fish and human. Data were acquired using release NCBI 36 of the human genome, release 1.0 of the medaka genome, release Zv7 of the zebrafish genome and release 1.0 of the stickleback genome. (PDF 23 KB)
Additional File 11: Medaka NITR24 is linked to CTLA4. (a) Human CTLA4 and CD28 are tightly linked. (b) In medaka CTLA4 is closely linked to NITR24. Gene names and chromosomal order were identified via Contig View within the Ensembl Genome Browser . Three predicted olfactory receptor (OR) genes are indicated on medaka chromosome 21, although the exact number of OR genes has not been determined. Data were acquired using release NCBI 36 of the human genome and release 1.0 of the medaka genome. CTLA4 and CD28 are not linked to NITR genes in zebrafish (not shown). Note that Ensembl misidentified medaka CTLA4  as CD28. Predicted protein structures are shown and structural features are as in Figure 4 plus cytoplasmic activating motifs (YxxM); and a functionally unresolved motif (YxxH). (PDF 44 KB)
Additional File 12: Phylogenetic comparisons of published NITR V domains. Cladogram depicts conventional phylogenetic radiations of bony fish [Figure adapted from 35]. Colors denote species in which NITR genes and/or gene transcripts have been characterized. (b) Neighbor-joining tree of V domains encoded by 73 representative NITR genes. Ol, medaka (Oryzias latipes); Dr, zebrafish (Danio rerio); Ip, channel catfish (Ictalurus punctatus); Om, rainbow trout (Oncorhynchus mykiss); Sn, Southern pufferfish (Spheroides nephelus); Po Japanese flounder (Paralichthys olivaceus) and Tn, pufferfish (Tetraodon nigroviridis). Specific NITR gene numbers are indicated (e.g. medaka NITR8a = Ol 8a), except for Tetraodon sequences (see Additional file 3). Only one NITR gene family member is incorporated; additional family members are indicated in parenthesis. Activating NITR genes are indicated by a "+"; NITR families that possess both inhibitory and activating forms are indicated by a "+/-". NITR genes which encode a V and I domain in a single exon ("V and I exons fused") and "No I domain" are indicated. These labels are presented in parentheses for NITR families, which do not consistently conform to these definitions. Number assigned to each interior branch corresponds to the bootstrap value; bootstrap values less than 60 are not shown. Branch lengths correspond to number of amino acid substitutions estimated by Poisson correction, indicated below. (PDF 51 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Desai, S., Heffelfinger, A.K., Orcutt, T.M. et al. The medaka novel immune-type receptor (NITR) gene clusters reveal an extraordinary degree of divergence in variable domains. BMC Evol Biol 8, 177 (2008). https://doi.org/10.1186/1471-2148-8-177
- Bony Fish
- Natural Killer Cell Receptor
- NITR Gene
- Zebrafish Chromosome
- Medaka Genome