Search

The leukocyte receptor complex (LRC) encodes numerous immunoglobulin (Ig)-like receptors involved in innate immunity. These include the killer-cell Ig-like receptors (KIR) and the leukocyte Ig-like receptors (LILR) which can be polymorphic and vary greatly in number between species. Using the recent long-read genome assembly, Sscrofa11.1, we have characterized the porcine LRC on chromosome 6. We identified a ~ 197-kb region containing numerous LILR genes that were missing in previous assemblies. Out of 17 such LILR genes and fragments, six encode functional proteins, of which three are inhibitory and three are activating, while the majority of pseudogenes had the potential to encode activating receptors. Elsewhere in the LRC, between FCAR and GP6, we identified a novel gene that encodes two Ig-like domains and a long inhibitory intracellular tail. Comparison with two other porcine assemblies revealed a second, nearly identical, non-functional gene encoding a short intracellular tail with ambiguous function. These novel genes were found in a diverse range of mammalian species, including a pseudogene in humans, and typically consist of a single long-tailed receptor and a variable number of short-tailed receptors. Using porcine transcriptome data, both the novel inhibitory gene and the LILR were highly expressed in peripheral blood, while the single KIR gene, KIR2DL1, was either very poorly expressed or not at all. These observations are a prerequisite for improved understanding of immune cell functions in the pig and other species. Electronic supplementary material The online version of this article (10.1007/s00251-018-1067-1) contains supplementary material, which is available to authorized users.

In this study, we investigated the immunoglobulin heavy (IGH) gene locus of torafugu (Takifugu rubripes) from publicly available assembly sequences and presented an annotated locus map, including the IGHV genes, pseudogenes, and IGHC genes. Three new IGHV gene families (IGHV3-IGHV5) were discovered. We observed the interspersion of IGHV1 and IGHV2 family members and that they often intermingled with each other, while other family members were further interspersed. Conservation of the promoter and recombination signal sequences (RSS) was observed in a family-specific manner. In addition to known variable region genes present on chromosome 5 (current torafugu genome assembly), we found 34 additional IGHV genes on scaffold 287 and three novel potentially functional IGHD genes on scaffold 483. In total, the variable region of the torafugu IGH locus consists of at least 48 IGHV genes, seven IGHD genes, and six IGHJ genes. IGHC genes have also been mapped in this study, with three genes encoding immunoglobulin classes: IgT, IgM, and IgD. We confirmed the expression of newly identified IGHV3 family sequences in the spleen and kidney of adult torafugu and found a favorable IGHV segment usage by IgM and IgT. Possible structural variation in the IGHδ locus was observed based on the current torafugu assembly. The complete characterization of the torafugu IGH locus will facilitate detailed studies of large-scale mechanisms associated with the recombination of the variable region genes and will offer insights into the genetic basis of the potential diversity in the antibody response observed in torafugu.

The activating natural cytotoxicity receptors on natural killer (NK) cells play a fundamental role in immunosurveillance of infections and cancer. Phylogenetic analyses showed that NKp30 is highly conserved in almost all jawed vertebrates and thus, represents one of the most ancient NK cell receptors. However, in contrast to other higher vertebrates, NKp30 is only a pseudogene in mouse, which contains two premature stop codons. To decipher the evolutionary role and biological function of NKp30 in mouse, we removed these premature stop codons and expressed the putative mouse NKp30 (mNKp30) protein as soluble Fc fusion construct and as full-length receptor on A5-GFP reporter cells. Interestingly, even though both NKp30 variants were expressed, maturation and targeting to the plasma membrane were impaired. Previous studies implicated that N-linked glycosylation is crucial for plasma membrane targeting and ligand binding of human NKp30. However, even though present in all other jawed vertebrates analyzed so far, these three N-linked glycosylation sites are missing in mouse NKp30. Interestingly, reconstitution of N-linked glycosylation enabled secretion of a mNKp30-Fc fusion protein which recognized a yet unknown ligand on the plasma membrane of mastocytoma cells. Based on these data, our study is the first to show expression and functional analysis of a mNKp30 protein suggesting that the mouse NKp30 pseudogene is the result of a species-specific loss of function.

The tripartite motif protein (TRIM)5α/CypA fusion protein TRIMCyp in Old World monkeys is generally considered unable to restrict HIV-1 replication. Monkeys with TRIMCyp can serve as a unique animal model for studies of HIV-1 infection. The present study investigated the distribution and expression status of TRIMCyp in four species of macaques originating from China and its borderlands: pigtail macaques (Macaca nemestrina), rhesus macaques (Macaca mulatta), long-tailed macaques (Macaca fascicularis), and Tibetan macaques (Macaca thibetana). The results revealed that the frequencies of the TRIMCyp genotype were significantly different among different species and even within different populations of the same species. Interestingly, the TRIMCyp genotype was more prevalent among macaques originating from Yunnan and surrounding regions than those from other regions of China. Importantly, TRIMCyp individuals were first identified in Chinese M. mulatta originating from Yunnan, although multiple earlier studies failed to find CypA retrotransposition in this subspecies. Furthermore, TRIMe7-CypA, one of the splicing isoforms of the TRIMCyp transcript was expressed in M. nemestrina and M. mulatta but not M. fascicularis. The intra- and interspecies polymorphisms in the deduced TRIMCyp amino acid sequences of these macaques were also analyzed. Taken together, the data in this study provide important information about the genomic background of TRIMCyp among major species of Chinese macaques.

T-cell immunoglobulin domain and mucin domain containing protein 1 (TIM1), also known as a cellular receptor for hepatitis A virus (HAVCR1) or a molecule induced by ischemic injury in the kidney (KIM1), is involved in the regulation of immune responses. We investigated a natural selection history of TIM1 by comparative sequencing analysis in 24 different primates. It was found that TIM1 had become a pseudogene in multiple lineages of the New World monkey. We also investigated T cell lines originated from four different New World monkey species and confirmed that TIM1 was not expressed at the mRNA level. On the other hand, there were ten amino acid sites in the Ig domain of TIM1 in the other primates, which were suggested to be under positive natural selection. In addition, mucin domain of TIM1 was highly polymorphic in the Old World monkeys, which might be under balanced selection. These data suggested that TIM1 underwent a lineage-specific evolutionary pathway in the primates.

The largest remaining carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is currently under threat of extinction due to a fatal contagious cancer-devil facial tumour disease. Low major histocompatibility complex (MHC) class I diversity is believed to have contributed to the transmission of the tumour allograft through devil populations. Here, we report low MHC class II variability in this species, with DA β chain genes (Saha-DAB1, 2 and 3) exhibiting very limited diversity and the sole α chain gene (Saha-DAA) monomorphic. Three, six and three alleles were found at Saha-DAB1, 2 and 3, respectively, with a predominant allele found at each locus. Heterozygosity at these three loci is low in the eastern population and modestly higher in northwestern individuals. The results are indicative of a selective sweep likely due to an infectious disease resulting in the fixation of selectively favoured alleles and depletion of genetic diversity at devil class II loci. Several attempts were made to isolate the other marsupial classical class II gene family, namely, DB, resulting in only one DBB pseudogene being found. These findings further support the view that this species has a compromised capacity to respond to pathogen evolution, emerging infectious diseases and environmental changes.

Immunoglobulin (IG) gene rearrangement and expression are central to disease resistance and health maintenance in animals. The IG kappa (IGK) locus in swine (Sus scrofa domestica) contributes to approximately half of all antibody molecules, in contrast to many other Cetartiodactyla, whose members provide the majority of human dietary protein and in which kappa locus utilization is limited. The porcine IGK variable locus is 27.9 kb upstream of five IG kappa J genes (IGKJ) which are separated from a single constant gene (IGKC) by 2.8 kb. Fourteen variable genes (IGKV) were identified, of which nine are functional and two are open reading frame (ORF). Of the three pseudogenes, IGKV3-1 contains a frameshift and multiple stop codons, IGKV7-2 contains multiple stop codons, and IGKV2-5 is missing exon 2. The nine functional IGKV genes are phylogenetically related to either the human IGKV1 or IGKV2 subgroups. IGKV2 subgroup genes were found to be dominantly expressed. Polymorphisms were identified on overlapping BACs derived from the same individual such that 11 genes contain amino acid differences. The most striking allelic differences are present in IGKV2 genes, which contain as many as 16 amino acid changes between alleles, the majority of which are in complementarity determining region (CDR) 1. In addition, many IGKV2 CDR1 are shared between genes but not between alleles, suggesting extensive diversification of this locus through gene conversion.

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens. Electronic supplementary material The online version of this article (doi:10.1007/s00251-010-0422-7) contains supplementary material, which is available to authorized users.

Killer cell immunoglobulin-like receptors (KIR) gene frequencies have been shown to be distinctly different between populations and contribute to functional variation in the immune response. We have investigated KIR gene frequencies in 370 individuals representing three Asian populations in Singapore and report here the distribution of 14 KIR genes (2DL1, 2DL2, 2DL3, 2DL4, 2DL5, 2DS1, 2DS2, 2DS3, 2DS4, 2DS5, 3DL1, 3DL2, 3DL3, 3DS1) with two pseudogenes (2DP1, 3DP1) among Singapore Chinese (n = 210); Singapore Malay (n = 80), and Singapore Indian (n = 80). Four framework genes (KIR3DL3, 3DP1, 2DL4, 3DL2) and a nonframework pseudogene 2DP1 were detected in all samples while KIR2DS2, 2DL2, 2DL5, and 2DS5 had the greatest significant variation across the three populations. Fifteen significant linkage patterns, consistent with associations between genes of A and B haplotypes, were observed. Eighty-four distinct KIR profiles were determined in our populations, 38 of which had not been described in other populations. KIR haplotype studies were performed using nine Singapore Chinese families comprising 34 individuals. All genotypes could be resolved into corresponding pairs of existing haplotypes with eight distinct KIR genotypes and eight different haplotypes. The haplotype A2 with frequency of 63.9% was dominant in Singapore Chinese, comparable to that reported in Korean and Chinese Han. The A haplotypes predominate in Singapore Chinese, with ratio of A to B haplotypes of approximately 3:1. Comparison with KIR frequencies in other populations showed that Singapore Chinese shared similar distributions with Chinese Han, Japanese, and Korean; Singapore Indian was found to be comparable with North Indian Hindus while Singapore Malay resembled the Thai.

Most of the previously studied teleost MHC class I molecules can be classified into two broad lineages: "U" and "Z/ZE." However, database reports on genes in cyprinid and salmonid fishes show that there is a third major lineage, which lacks detailed analysis so far. We designated this lineage "L" because of an intriguing linkage characteristic. Namely, one zebrafish L locus is closely linked with MHC class II loci, despite the extensively documented nonlinkage of teleost class I with class II. The L lineage consists of highly variable, nonclassical MHC class I genes, and has no apparent orthologues outside teleost fishes. Characteristics that distinguish the L lineage from most other MHC class I are (1) absence of two otherwise highly conserved tryptophan residues W51 and W60 in the alpha1 domain, (2) a low GC content of the alpha1 and alpha2 exons, and (3) an HINLTL motif including a possible glycosylation site in the alpha3 domain. In rainbow trout (Oncorhynchus mykiss) we analyzed several intact L genes in detail, including their genomic organization and transcription pattern. The gene Onmy-LAA is quite different from the genes Onmy-LBA, Onmy-LCA, Onmy-LDA, and Onmy-LEA, while the latter four are similar and categorized as "Onmy-LBA-like." Whereas the Onmy-LAA gene is organized like a canonical MHC class I gene, the Onmy-LBA-like genes are processed and lack all introns except intron 1. Onmy-LAA is predominantly expressed in the intestine, while the Onmy-LBA-like transcripts display a rather homogeneous tissue distribution. To our knowledge, this is the first description of an MHC class I lineage with multiple copies of processed genes, which are intact and transcribed. The present study significantly improves the knowledge of MHC class I variation in teleosts.

Many genes and proteins are required to carry out the processes of innate and adaptive immunity. For many studies, including systems biology, it is necessary to have a clear and comprehensive definition of the immune system, including the genes and proteins that take part in immunological processes. We have identified and cataloged a large portion of the human immunology-related genes, which we call the essential immunome. The 847 identified genes and proteins were annotated, and their chromosomal localizations were compared to the mouse genome. Relation to disease was also taken into account. We identified numerous pseudogenes, many of which are expressed, and found two putative new genes. We also carried out an evolutionary analysis of immune processes based on gene orthologs to gain an overview of the evolutionary past and molecular present of the human immune system. A list of genes and proteins were compiled. A comprehensive characterization of the member genes and proteins, including the corresponding pseudogenes is presented. Immunome genes were found to have three types of emergence in independent studies of their ontologies, domains, and functions.

In cattle, there are six classical MHC class I genes that are variably present between different haplotypes. Almost all known haplotypes contain between one and three genes, with an allele of Gene 2 present on the vast majority. However, very little is known about the sequence and therefore structure and evolutionary history of this genomic region. To address this, we have refined the MHC class I region in the Hereford cattle genome assembly and sequenced a complete A14 haplotype from a homozygous Holstein. Comparison of the two haplotypes revealed extensive variation within the MHC class Ia region, but not within the flanking regions, with each gene contained within a conserved 63- to 68-kb sequence block. This variable region appears to have undergone block gene duplication and likely deletion at regular breakpoints, suggestive of a site-specific mechanism. Phylogenetic analysis using complete gene sequences provided evidence of allelic diversification via gene conversion, with breakpoints between each of the extracellular domains that were associated with high guanine-cytosine (GC) content. Advancing our knowledge of cattle MHC class I evolution will help inform investigations of cattle genetic diversity and disease resistance.

The catfish IGH locus is large ( approximately 1 Mb) and complex, having undergone multiple internal duplications and transpositions. To define the structure of the locus that contains the single expressed IGHM gene, two overlapping bacterial-artificial-chromosome (BAC) clones spanning the most 3' end of the channel catfish immunoglobulin heavy (IGH) chain locus have been completely sequenced. The analyses created a contig of 257,153 bp containing 55 VH, 6 D, 12 JH genes and the IGH constant region genes encoding the functional secreted and membrane forms of IgM and the membrane form of IgD. This analysis revealed three major features. First, no C-region genes were found aside from the previously described IGHM1 and IGHD1, with the latter gene being the most 3' C-region gene of the catfish IGH locus. There was no evidence in the region sequenced for genes that could encode an Ig class similar to the IgZ/IgT described in zebrafish, trout and pufferfish. Second, there are a high number of VH pseudogenes, 28 out of 55 (51%). In contrast, the entire zebrafish IGH locus has 40 functional VH genes and eight pseudogenes (17%). Third, an internal duplication of a 52.4-kb block of VH genes has occurred. These observations suggest that the IGH locus of teleost fish varies significantly from species to species in the diversity of C-region genes as well as the numbers of genes encoding V regions.

As part of an ongoing elucidation of rainbow trout major histocompatibility complex (MHC) class I, the polymorphism of two MHC class Ib loci was analyzed. These loci, Onmy-UCA and Onmy-UDA, are situated head-to-tail and share more than 89% nucleotide identity in their open reading frames. They share 80% identity with some trout Ia alleles. The deduced amino acid sequences suggest that the UCA and UDA molecules are transported to endosomal compartments and may bind peptides in their binding groove. Our survey revealed seven UCA and eight UDA alleles. Similarity indices overlap when comparing within and between UCA and UDA alleles and some cross-locus motif variation is observed. In most trout both UCA and UDA transcripts were found. However, there probably is functional redundancy, because some trout lacked transcription of one of the two loci. Furthermore, for some UCA and UDA alleles, splicing deficiencies, early stop codons, and upstream start codons were found, which may interfere with efficient protein expression. The present study is the first extensive report on MHC class Ib polymorphism assigned to locus in ectotherm species.

CD45 is a transmembrane protein tyrosine phosphatase, which in mammals plays an important role in T and B cell receptor and cytokine signaling. Recently, a catfish cDNA was shown to contain all characteristic CD45 features: an alternatively spliced amino-terminus, a cysteine-rich region, three fibronectin domains, a transmembrane region, and two phosphotyrosine phosphatase domains. However, analyses of CD45 cDNAs from various catfish lymphoid cell lines demonstrated that catfish CD45 is unique in that it contains a large number of alternatively spliced exons. Sequence analyses of cDNAs derived from the catfish clonal B cell line 3B11 indicated that this cell line expresses up to 13 alternatively spliced exons. Furthermore, sequence similarity among the alternatively spliced exons suggested duplication events. To establish the exact number and organization of alternatively spliced exons, a bacterial artificial chromosome library was screened, and the catfish functional CD45 gene plus six CD45 pseudogenes were sequenced. The catfish functional CD45 gene spans 37 kb and contains 49 exons. In comparison, the human and pufferfish CD45 genes consist of 34 and 30 exons, respectively. This difference in the otherwise structurally conserved catfish gene is due to the presence of 18 alternatively spliced exons that were likely derived through several duplication events. In addition, duplication events were also likely involved in generating the six pseudogenes, truncated at the 3' ends. A similarly 3' truncated CD45 pseudogene is also present in the pufferfish genome, suggesting that this specific CD45 gene duplication occurred before catfish and pufferfish diverged (approximately 400 million years ago).

Major histocompatibility complex (MHC) class I genes have complicated and profound evolutionary histories. To reconstruct and better understand their histories, partial class I genes (exon 2-intron 2-exon 3) were sequenced in a sampling of prosimians (Strepsirhini, Primates). In total, we detected 117 different sequences from 36 Malagasy prosimians (lemurs) and 1 non-Malagasy prosimian (galago) representing 4 families, 7 genera, and 13 species. Unlike the MHC class II genes ( MHC-DRB), MHC class I genes show a generally genus-specific mode of evolution in lemurs. Additionally, no prosimian class I loci were found to be orthologous to HLA genes, even at highly conserved loci (such as HLA-E, HLA-F). Phylogenetic analysis indicates that nucleotide diversity among loci was very small and the persistence time of the polymorphisms was short, suggesting that the origin of the lemur MHC class I genes detected in this study was relatively recent. The evolutionary mode of these genes is similar to that of classic HLA genes, HLA-A, HLA-B, and HLA-C, in terms of their recent origin and rarity of pseudogenes, and differs from them with respect to the degree of gene duplications. From the viewpoint of MHC genes evolution, some interlocus sequence exchanges were apparently observed in the lemur lineage upon phylogenetic and amino acid motif analyses. This is also in contrast to the evolutionary mode of HLA genes, where intralocus exchanges have certainly occurred but few interlocus exchanges have taken place. Consequently, the gene conversion model for explaining the generation of the MHC diversity among different loci can be thought to play more important roles in the evolution of lemur MHC class I genes than in that of HLA genes.

Conventional phylogenetic trees for the human leukocyte antigen (HLA)-DRB1 alleles constructed by the neighbor-joining (Saitou and Nei 1987) and UPGMA (Sneath and Sokal 1973) methods using nucleotide sequences of the DRB1 alleles suggest that DRB1*0701 may have diverged from other DRB1 alleles before the separation of the human and chimpanzee species, because of a large number of nucleotide changes in DRB1*0701 compared with any of the other DRB1 alleles. Here we show new evidence that the haplotypes centering on DRB1*0701 and DRB1*04 alleles are the most homologous. This suggests that these haplotypes have derived from the common ancestral haplotype, and that they have likely retained complete linkage disequilibrium even after the divergence of the DRB1*0701 and DRB1*04 allelic lineages. Together with the corresponding haplotype carrying chimpanzee DRB1*0701, which has a high sequence homology to HLA-DRB1*0701, these haplotypes reveal that: (1) the DRB1*04 allelic lineage may have been generated from the DRB1*0701 lineage after the separation of the human and chimpanzee species; (2) the DRB1*04 allelic lineage possibly has a higher substitution rate of DRB1 compared with pseudogene and neutral region; (3) there could be a significant difference in the substitution rate of DRB1 between the DRB1*0701 and DRB1*04 allelic lineages. Based on the difference between the present and previous results, we would like to propose that phylogenetic studies using not only nucleotide sequences of the DRB1 alleles but also haplotypes centering on the alleles should be conducted for understanding detailed phylogenetic relationships of the DRB1 alleles.

To date, there are no published MHC sequences from the California sea lion (Zalophus californianus), a thriving species that, by feeding high on the marine food web, could be a sentinel for disturbances in marine and coastal ecosystems. In this study, degenerate primers and RACE technology were used to amplify near-full-length (MhcZaca- DQB) and full-length (MhcZaca- DQA) expressed class II MHC gene products from the peripheral blood mononuclear cells of two California sea lions in rehabilitation. Five unique Zaca- DQA sequences and eight unique Zaca- DQB sequences, all encoding functional proteins, were identified in the two animals, indicating the presence of multiple DQ- loci in this species. An additional three Zaca- DQB sequences containing features compatible with pseudogenes or null alleles were also identified. Despite the identification of multiple DQA and DQB sequences, the degree of heterogeneity between them was extremely low. To confirm the limited degree of Zaca-DQ nucleotide variation between individuals, we used denaturing gradient gel electrophoresis to examine putative peptide binding region sequences from the peripheral blood leukocyte-derived RNAs of 19 wild-caught California sea lions from physically distinct populations. The pattern of Zaca-DQ sequence migration was identical between individuals and independent of geographical region. This apparent Zaca-DQ sequence identity between sea lions was confirmed by direct sequencing of individual bands. In combination, these findings raise important questions regarding immunogenetic diversity within this thriving species, and should prompt further research into the existence of a highly polymorphic sea lion class II MHC molecule with sequence features that support traditional peptide binding functions.

The murine Nkrp1 gene family encodes three previously identified activation and inhibitory receptors expressed on natural killer (NK) cells. This family includes the gene for NKR-P1C (NK1.1), the most specific serologic marker on C57BL/6-derived NK cells and is localized in a gene cluster in the NK gene complex (NKC). To further analyze the Nkrp1 family, we constructed and analyzed a bacterial artificial chromosome contig. A genomic organization of the Nkrp1 family was obtained and three new Nkrp1 genes were isolated from interleukin-2-activated NK cells. Thus, the Nkrp1 family adds to the repertoire of receptors expressed by NK cells.

The kinesin-related protein (HSET) gene belongs to the kinesin superfamily, the members of which are involved in cellular transport processes. The HSET gene product was previously characterized by partial cDNA sequencing. The gene is located on the short arm of human Chromosome 6 (6p21.3), at the centromeric end of the major histocompatibility complex. Here, we report the genomic structure of the complete HSET gene together with its flanking loci. Sequence analysis of the 40 kilobase (kb) cosmid clone containing the HSET gene also revealed the presence of several new genes not related to the kinesin superfamily. These include a 60S ribosomal protein L35A-like pseudogene (rPL35A-like) on the telomeric side and a polycomb-like gene (PHF1), a copper tolerance-like gene (CUTA1) and the 5' part of the synaptic ras-GTPase-activating protein (SynGAP) gene centromeric of HSET. In addition, a complete 60S ribosomal protein L12-like (rPL12L) gene in intron 3 of the HSET gene was identified which appears to have an open reading frame. The possible involvement of the HSET gene and a beta-tubulin gene (TUBB) in the pathogenesis of immotile cilia syndrome (ICS) was studied by screening two unrelated ICS families with microtubular defects and suspected HLA linkage for mutations within the HSET gene and the TUBB gene. Four single base substitutions were detected in the HSET gene, and none in the TUBB gene. On the basis of these data, a role of the HSET and TUBB products in the pathogenesis of ICS in the two families is unlikely.

The V10 variable gene of the human T-cell receptor gamma locus (TCRG-V10), the only member of the subgroup III, has a structural defect which inhibits the splicing of the leader intron. We show that there is a single point mutation in the V10 leader donor splice site responsible for this situation and that this mutation is found in the different populations tested, indicating that V10 corresponds to a pseudogene in humans. We restored the splice site by mutagenesis and obtained correct splicing in vitro. Analysis of the V10 germline gene in different primates reveals functional splice sites in the closest human apes, the chimpanzee and the gorilla. The splice competence of TCRG-V10 in higher primates was addressed in peripheral blood lymphocytes from chimpanzee by specific cDNA amplification, and correct splicing of the TCRG-V10 leader intron was found as well as a majority of in frame rearrangements involving only the TCRG-J1 or J2 segments. These results suggest that V10(+)gamma /delta T cells may represent an important subset in the non-human higher primates, contrary to the situation observed in the human.

The major histocompatibility complex (MHC) is integral to the vertebrate adaptive immune system. Characterizing diversity at functional MHC genes is invaluable for elucidating patterns of adaptive variation in wild populations, and is particularly interesting in species of conservation concern, which may suffer from reduced genetic diversity and compromised disease resilience. Here, we use next generation sequencing to investigate MHC class II B (MHCIIB) diversity in two sister taxa of New Zealand birds: South Island saddleback (SIS), Philesturnus carunculatus, and North Island saddleback (NIS), Philesturnus rufusater. These two species represent a passerine family outside the more extensively studied Passerida infraorder, and both have experienced historic bottlenecks. We examined exon 2 sequence data from populations that represent the majority of genetic diversity remaining in each species. A high level of locus co-amplification was detected, with from 1 to 4 and 3 to 12 putative alleles per individual for South and North Island birds, respectively. We found strong evidence for historic balancing selection in peptide-binding regions of putative alleles, and we identified a cluster combining non-classical loci and pseudogene sequences from both species, although no sequences were shared between the species. Fewer total alleles and fewer alleles per bird in SIS may be a consequence of their more severe bottleneck history; however, overall nucleotide diversity was similar between the species. Our characterization of MHCIIB diversity in two closely related species of New Zealand saddlebacks provides an important step in understanding the mechanisms shaping MHC diversity in wild, bottlenecked populations.

The human HLA-DR3 haplotype consists of two functional genes (DRB1*03 and DRB3*01) and one pseudogene (DRB2), arranged in the order DRB1... DRB2... DRB3 on the chromosome. To shed light on the origin of the haplotype, we sequenced 1480 nucleotides of the HLA-DRB2 gene and aong stretches of two other genes, Gogo-DRB2 from a gorilla, “Sylvia” and Patr-DRB2 from a chimpanzee, “Hugo”. All three sequences (HLA-DRB2, Gogo-DRB2, Patr-DRB2) are pseudogenes. The HLA-DRB2 and Gogo-DRB2 pseudogenes lack exon 2 and contain a twenty-nucleotide deletion in exon 3, which destroys the correct translational reading frame and obliterates the highly conserved cysteine residue at position 173. The Patr-DRB2 pseudogene lacks exons 1 and 2; it does not contain the twenty-nucleotide deletion, but does contain a characteristic duplication of that part of exon 6 which codes for the last four amino acid residues of the cytoplasmic region. When the nucleotide sequences of these three genes are compared to those of all other known DRB genes, the HLA-DRB2 is seen as most closely related to Gogo-DRB2, indicating orthologous relationship between the two sequences. The Patr-DRB2 gene is more distantly related to these two DRB2 genes and whether it is orthologous to them is uncertain. The three genes are in turn most closely related to HLA-DRBVI (the pseudogene of the DR2 haplotype) and Patr-DRB6 (another pseudogene of the Hugo haplotype), followed by HLA-DRB4 (the functional but nonpolymorphic gene of the DR4 haplotype). These relationships suggest that these six genes evolved from a common ancestor which existed before the separation of the human, gorilla, and chimpanzee lineages. The DRB2 and DRB6 have apparently been pseudogenes for at least six million years (myr). In the human and the gorilla haplotype, the DRB2 pseudogene is flanked on each side by what appear to be related genes. Apparently, the DR3 haplotype has existed in its present form for more than six myr.

The HLA-DR2 haplotype contains three beta-chain encoding DRB genes and one alpha-chain encoding DRA gene. Of the three DRB genes, two are presumably functional (HLA-DRB1 and HLA-DRB5), whereas the third (HLA-DRBVI) is a pseudogene. A pseudogene closely related to HLA-DRBVI is present in the chimpanzee (Patr-DRB6) and in the gorilla (Gogo-DRB6). We sequenced the HLA-DRBVI and Patr-DRB6 pseudogenes (all exons and most of the introns), and compared the sequence to that of the Gogo-DRB6 gene (of which only the exon sequence is available). All three pseudogenes seem to lack exon 1 and contain other deletions responsible for shifts in the translational reading frame. At least the HLA-DRBVI pseudogene, however, seems to be transcribed nevertheless. The chimpanzee pseudogene contains two inserts in intron 2, one of which is an Alu repeat belonging to the Sb subfamily, while the other remains unidentified. These inserts are lacking in the human gene. A comparison with sequences published by other investigators revealed the presence of the HLA-DRBVI pseudogene also in the DR1 and DRw10 haplotypes. Measurements of genetic distances indicate DRB6 to be closely related to the DRB2 pseudogene and to the HLA-DRB4 functional gene. In humans, gorillas, and chimpanzees, the DRB6 pseudogene is associated with the same functional gene (DRB5) indicating that this linkage disequilibrium is at least six million years old and that DR2 is one of the oldest DR haplotypes in higher primates.

The CD33-related sialic acid binding Ig-like lectins (CD33rSiglecs) are predominantly inhibitory receptors expressed on leukocytes. They are distinguishable from conserved Siglecs, such as Sialoadhesin and MAG, by their rapid evolution. A comparison of the CD33rSiglec gene cluster in different mammalian species showed that it can be divided into subclusters, A and B. The two subclusters, inverted in relation to each other, each encode a set of CD33rSiglec genes arranged head-to-tail. Two regions of strong correspondence provided evidence for a large-scale inverse duplication, encompassing the framework CEACAM-18 (CE18) and ATPBD3 (ATB3) genes that seeded the mammalian CD33rSiglec cluster. Phylogenetic analysis was consistent with the predicted inversion. Rodents appear to have undergone wholesale loss of CD33rSiglec genes after the inverse duplication. In contrast, CD33rSiglecs expanded in primates and many are now pseudogenes with features consistent with activating receptors. In contrast to mammals, the fish CD33rSiglecs clusters show no evidence of an inverse duplication. They display greater variation in cluster size and structure than mammals. The close arrangement of other Siglecs and CD33rSiglecs in fish is consistent with a common ancestral region for Siglecs. Expansion of mammalian CD33rSiglecs appears to have followed a large inverse duplication of a smaller primordial cluster over 180 million years ago, prior to eutherian/marsupial divergence. Inverse duplications in general could potentially have a stabilizing effect in maintaining the size and structure of large gene clusters, facilitating the rapid evolution of immune gene families.

The major histocompatibility complex (MHC) is an immunological gene-dense region of high diversity in mammalian species. Sus scrofa was domesticated by at least six independent events over Eurasia during the Holocene period. It has been hypothesized that the level and distribution of MHC variation in pig populations reflect genetic selection and environmental influences. In an effort to define the complexity of MHC polymorphisms and the role of selection in the generation of class II gene diversity (DQB, DRB1, and pseudogene PsiDRB3), DNA from globally distributed unrelated domestic pigs of European and Asian origins and a Suidae out-group was analyzed. The number of pseudogene alleles identified (PsiDRB3 33) was greater than those found in the expressed genes (DQB 20 and DRB1 23) but the level of observed heterozygosity (PsiDRB3 0.452, DQB 0.732, and DRB1 0.767) and sequence diversity (PsiDRB3 0.029, DQB 0.062, and DRB1 0.074) were significantly lower in the pseudogene, respectively. The substitution ratios reflected an excess of d (N) (DQB 1.476, DRB1 1.724, and PsiDRB3 0.508) and the persistence of expressed gene alleles suggesting the influence of balancing selection, while the pseudogene was undergoing purifying selection. The lack of a clear MHC phylogeographic tree, coupled with close genetic distances observed between the European and Asian populations (DQB 0.047 and DRB1 0.063) suggested that unlike observations using mtDNA, the MHC diversity lacks phylogeographic structure and appears to be globally uniform. Taken together, these results suggest that, despite regional differences in selective breeding and environments, no skewing of MHC diversity has occurred.

A number of statistical methods are widely used to describe allelic variation at specific genetic loci and its implication on the evolutionary history of these loci. Although the methods were developed primarily to study allelic variation at loci that are virtually always present in the genome, they are often applied to data of gene content variation (i.e., presence/absence of multiple homologous genes) at the killer cell immunoglobulin-like receptor (KIR) gene cluster. In this paper, we discuss methodological issues involved in the analysis of gene content variation data in the KIR region and also its covariation with polymorphism at the human leukocyte antigen class I loci, which encode ligands for KIR. A comparison of several statistical methods and measures (gene frequency, haplotype frequency, and linkage disequilibrium estimation) using the Centre d’Etude du Polymorphisme Humain data will be provided using KIR haplotypes that have been determined by segregation analysis, noting the strengths and weaknesses of the methods when only the presence/absence data is considered. Finally, application of these methods to a set of globally distributed populations is described (see Single et al., Nat Genet 39:1114–1119, 2007) in order to illustrate the challenges faced when inferring the joint effects of natural selection and demographic history on these immune-related genes.

The human T-cell receptor (TCR) alpha/delta variable loci are interspersed on the chromosome 14q11 and consist of 57 intergenic spaces ranging from 4 to 100 kb in length. To elucidate the evolutionary history of this locus, we searched the intergenic spaces of all TCR alpha/delta variable (TRAV/DV) genes for pseudogenes and potential protein-coding genes. We applied direct open reading frame (ORF) searches, an exon-finding algorithm and comparative genomics. Two TRAV/DV pseudogenes were discovered bearing 80 and 65% sequence similarity to TRAV14DV4 and TRAV9-1/9-2 genes, respectively. A gene bearing 85% sequence identity to B lymphocyte activation-related protein, BC-1514, upstream of TRAV26-2 was also discovered. This ORF (BC-1514tcra) is a member of a gene family whose evolutionary history and function are not known. In total, 36 analogs of this gene exist in the human, the chimpanzee, the Rhesus monkey, the frog and the zebrafish. Phylogenetic analyses show convergent evolution of these genes. Assays for the expression of BC-1514tcra revealed transcripts in the bone marrow, thymus, spleen, and small intestine. These assays also showed the expression of another analog to BC-1514, found on chromosome 5 in the bone marrow and thymus RNA. The existence of at least 17 analogs at various locations in the human genome and in nonsyntenic chromosomes of the chimpanzee suggest that BC-1514tcra, along with its analogs may be transposable elements with evolved function(s). The identification of conserved putative serine phosphorylation sites provide evidence of their possible role(s) in signal transduction events involved in B cell development and differentiation.

The availability of a contig of bacterial artificial chromosome (BAC) clones spanning the equine major histocompatibility complex (MHC) made possible a detailed analysis of horse MHC class I genes. Prior to this study, only a single horse MHC class I gene had been sequenced at the genomic level. Although many ( approximately 60) MHC class I cDNA sequences had been determined and published, from this information, it was not possible to determine how many class I loci are expressed in horses or to assign individual sequences to allelic series. In this study, 15 MHC class I genes were identified in BAC subclones and fully sequenced. Because the BAC library donor horse had been bred for homozygosity at the MHC, these 15 genomic clones represent distinct MHC class I genes and pseudogenes and not alleles at a smaller number of loci. For five of the genes, cDNA sequences from these loci had previously been identified. Two additional expressed class I genes were discovered, bringing the known total of different equine MHC class I genes (loci) expressed as mRNA to seven. Expression of all seven loci was detected by reverse transcriptase-polymerase chain reaction in adult, fetal, and placental tissues. The remaining eight genes were designated as pseudogenes. This work resulted in moderate expansion of the horse MHC BAC contig length, and the remaining gap was shortened. The information contained in these equine MHC class I sequences will permit comparison of MHC class I genes expressed across different horse MHC haplotypes and between horses and other mammalian species.

More than 100 variable (V), 27 diversity (D), and six joining (J) genes are encoded in the human heavy chain locus, and many of these genes exists in different allelic forms. The number of genes and the allelic differences help to create diversity in the immunoglobulin receptors, a key feature of the adaptive immune system. We here report the identification of two novel and seemingly functional alleles of human heavy chain genes. The variable IGHV3-23*04 allele is found with an allele frequency of 0.21 amongst Danish Caucasians, whereas the novel joining IGHJ6*04 allele is rare (allele frequency 0.02). We also report the full sequence of IGHV3-h. The gene exists in two allelic forms but is only found in 58% of the Danish Caucasians studied. The methionine translation initiation codon is mutated, ATG--AAG, and we therefore propose that the gene is a pseudogene incapable of being translated.

The most telomeric class I region of the MHC in rat and mouse is the M region, which contains about 20 class I genes or gene fragments. The central part carries three class I genes- M4, M5, and M6-which are orthologous between the two species. M4 and M6 are pseudogenes in the mouse but transcribed, intact genes in the rat. To analyze the pseudogene status for the mouse genes in more detail, we have sequenced the respective exons in multiple representative haplotypes. The stop codons are conserved in all mouse strains analyzed, and, consistent with the pseudogene status, all strains show additional insertions and deletions, taking the genes further away from functionality. Thus, M4 and M6 indeed have a split status. They are silent in the mouse but intact in the closely related rodent, the rat.

Genome analysis of the swine leukocyte antigen ( SLA) region is needed to obtain information on the MHC genomic sequence similarities and differences between the swine and human, given the possible use of swine organs for xenotransplantation. Here, the genomic sequences of a 433-kb segment located between the non-classical and classical SLA class I gene clusters were determined and analyzed for gene organization and contents of repetitive sequences. The genomic organization and diversity of this swine non-class I gene region was compared with the orthologous region of the human leukocyte antigen ( HLA) complex. The length of the fully sequenced SLA genomic segment was 433 kb compared with 595 kb in the corresponding HLA class I region. This 162-kb difference in size between the swine and human genomic segments can be explained by indel activity, and the greater variety and density of repetitive sequences within the human MHC. Twenty-one swine genes with strong sequence similarity to the corresponding human genes were identified, with the gene order from the centromere to telomere of HCR - SPR1 - SEEK1 - CDSN - STG - DPCR1 - KIAA1885 - TFIIH - DDR - IER3 - FLOT1 - TUBB - KIAA0170 - NRM - KIAA1949 - DDX16 - FLJ13158 - MRPS18B - FB19 - ABCFI - CAT56. The human SEEK1 and DPCR1 genes are pseudogenes in swine. We conclude that the swine non-class I gene region that we have sequenced is highly conserved and therefore homologous to the corresponding region located between the HLA-C and HLA-E genes in the human.

Genomic DNA from a panel of 41 healthy unrelated Japanese individuals was typed for the presence or absence of 16 KIR genes and pseudogenes. Only eight different KIR genotypes were found. Group A haplotypes outnumbered group B haplotypes in frequency by approximately 3:1, with individuals having two group A haplotypes accounting for 56% of the panel. The frequency of A haplotypes in the Japanese is higher than that observed in other populations. Flow cytometric comparison of KIR expression in 19 panel members showed considerable diversity in NK cell repertoire, which was also seen within the group of individuals having two A haplotypes. This diversity is likely due to allelic polymorphism in expressed genes of the A haplotype. In comparison to other populations, the Japanese appear less heterogeneous in KIR genotype as assessed by gene content.

We report the isolation and characterization of cDNA clones of expressed, functional major histocompatibility complex class-I ( Mhc-I) genes from two species of monotremes: the duck-billed platypus and the short-beaked echidna. The cDNA clones were isolated from libraries constructed from spleen RNA, clearly establishing their expression in at least this one peripheral lymphoid organ. From the presence of conserved amino acid residues, it appears the expressed sequences encode molecules that likely function as classical Mhc-I. These clones were isolated using monotreme Mhc-I processed pseudogenes as probes. These processed pseudogenes were isolated from genomic DNA and, based on their structure, are likely independently derived in the platypus and echidna. When all the monotreme sequences were included in phylogenetic analyses, we found no apparent orthologous relationships between the platypus and echidna Mhc-I. Analyses that included a large number of Mhc-I sequences from other taxa support a separate monotreme Mhc-I clade, basal to a therian Mhc-I clade that is comprised of sequences from marsupial and placental mammals. The phylogenies also support the hypothesis that Mhc-I genes of placental mammals, marsupials, and monotremes are derived from three separate lineages of Mhc-I genes, best explained by two rounds of duplications and deletions. The first round would have occurred prior to the divergence of monotremes and therians, and the second prior to the divergence of marsupials and placental mammals. The sequences described here represent the first reported functional monotreme Mhc-I, as well as the first processed pseudogenes of any type from monotremes.

Human killer-cell immunoglobulin-like receptors (KIR) show three types of organization of their extracellular domains: D0-D1-D2 in KIR3D, D1-D2 in the majority of KIR2D, and D0-D2 in KIR2DL4 and the novel KIR2DL5. The gene for a KIR2DL3 variant, which has a D1-D2 structure, has been shown previously to have a nonexpressed region (pseudoexon 3) that is paralogous to the exon encoding the D0 domain of other KIR. This pseudoexon is not expressed because it is skipped during splicing of pre-mRNA. In this study, we demonstrate that all eight genes encoding human KIR with D1-D2 configuration (KIR2DL1-KIR2DL3, KIR2DS1-KIR2DS5) have similarly untranslated pseudoexons. Whereas the pseudoexons of four of these KIR genes bear nonsense mutations and/or altered splicing sites, the pseudoexons in the other four KIR genes have no major structural abnormalities, indicating that other mechanisms are responsible for inactivation of their exons 3. A comparison of the sequences on pseudoexons 3 with the paralogous expressed exons suggests that an exonic splicing enhancer may be necessary for the expression of exon 3 in KIR genes.

Three rat BV13S1 alleles (T-cell receptor beta-chain variable gene 13) were characterized by new BV13S1-allele specific monoclonal antibodies (18B1 and 17D5) and sequence analysis of expressed and genomic BV13S1. Two alleles were functional and designated BV13S1A1 present in strains LEW, BUF, PVG, and BV13S1A2 present in BN and WF. Their products differed by six amino acids, two of them in complementarity-determing region (CDR)1 and one in CDR2. A third nonfunctional allele, BV13S1A3P, was found in strains F344 and DA. Apart from a single nucleotide insertion, it was identical to BV13S1A2. All 12 rat strains tested showed association of TCRBC1 with BV8S2/4 alleles but not with the BV13S1 alleles, which may reflect a different gene order of the rat BV compared to mouse. BV13S1A1-encoded T-cell receptors (TCRs) which bind both monoclonal antibody (mAb) 18B1 and mAb 17D5 are over-represented in the CD4 lymphocyte subset. BV13S1A2-encoded TCRs which are stained by mAb 18B1 but not by mAb 17D5 show a slight CD8-biased expression. Preferential usage of BV13S1A1-positive TCRs by CD4 but not by CD8 cells in (LEW x WF)F1 hybrids and cosegregation of BV13SA1 and increased frequency of BV13S1 TCR-positive CD4 cells in a (LEW x BN)x BN backcross suggest structural differences of the two allelic products as the reason for their contrasting CD4/CD8 subset bias.

Lymphocytes of a New World primate, the cotton-top tamarin (Saguinus oedipus), express classical G–related major histocompatibility complex (MHC) class I molecules with unusually limited polymorphism and variability. Three G-related loci, an F locus, an E locus, and two pseudogenes (So-N1 and So-N3) have been identified by cDNA library screening and extensive PCR analysis of both cDNA and genomic DNA from the cotton-top tamarin. Furthermore, each genus of the subfamily Callitrichinae (tamarins and marmosets) appears to express its own unique set of MHC class I genes, likely due to a rapid turnover of loci. The rapid emergence of unique MHC class I genes in the Callitrichinae genera, resulting from an active process of duplication and inactivation of loci, may account for the limited diversity of the MHC class I genes in the cotton-top tamarin. To determine the nature of the entire complement of MHC class I genes in the cotton-top tamarin, we synthesized a genomic DNA library and screened it with MHC class I-specific probes. We isolated nine new MHC class I pseudogenes from this library. These newly isolated tamarin G–related MHC class I pseudogenes are not closely related to any of their functional counterparts in the tamarin, suggesting that they do not share a recent common ancestral gene with the tamarin's currently expressed MHC class I loci. In addition, these tamarin sequences display a high rate of nonsynonymous substitutions in their putative peptide binding region. This indicates that the genes from which they have derived were likely subject to positive selection and, therefore, were once functional. Our data support the notion that an extremely high rate of loci turnover is largely responsible for the limited diversity of the MHC class I genes in the cotton-top tamarin.