|
|
||||||||
Departments of Structural Biology and Microbiology and Immunology, Stanford University School of Medicine, Stanford, CA 94305
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The CD94 and NKG2A polypeptides form a heterodimeric receptor (13, 14) that recognizes the MHC nonclassical class I HLA-E ligand complexed with peptides derived from the leader sequences of certain classical HLA-A, B, C, and nonclassical HLA-G H chains (15, 16, 17, 18, 19). Ligand recognition transduces inhibitory signals in NK cells via the immunomodulatory tyrosine-based inhibition motifs on the cytoplasmic tail of NKG2A (20, 21, 22). The CD94:NKG2C heterodimer is a receptor with ligand specificity similar to CD94:NKG2A, but which generates activating signals within NK cells (22, 23). Human NKG2D is a homodimeric activating receptor that interacts with MHC class I-like chain A (24, 25), a MHC class I-like protein with considerable polymorphism, and with unique long binding proteins (ULBPs), newly identified molecules distantly related to MHC class I (26). In NK cells, activating signals are conveyed via positively charged amino acid residues in the transmembrane domains of activating receptors (NKG2C, D) that interact with signal-transducing adaptor molecules such as DAP-10 and DAP-12 (27, 28). The exact roles of human NKG2E and NKG2F are uncertain, but both encode a positively charged lysine residue in the transmembrane domain, as does NKG2C.
The KIR genes of the human LRC encode receptors that are specific for polymorphic determinants of HLA-A, B, and C, and also for HLA-G (1, 29). Previous work from our laboratory revealed diversity of KIR haplotype in the human population, with respect to the number and content of KIR genes (30) and also allelic polymorphism at individual genes (31, 32, 33, 34). Furthermore, comparison of human and chimpanzee KIR revealed many species-specific features to the KIR gene systems (35, 36). Given the diversity and rapid evolution of the KIR genes it became of particular interest to investigate intraspecies and interspecies diversity of NKC genes that encode NK cell receptors for MHC class I molecules. Therefore, we have investigated the diversity and variation in the CD94 and NKG2 genes of humans and apes.
| Materials and Methods |
|---|
|
|
|---|
Healthy human donors recruited from our laboratory represent a
variety of ethnic and racial backgrounds, as shown in Fig. 1
. Tissue
samples from non-human primates were purchased from the Yerkes Regional
Primate Center at Emory University (Atlanta, GA) and from the
Laboratory for Experimental Medicine and Surgery in Primates (New York
UniversityMedical Center, Tuxedo, NY). Blood was drawn from
animal subjects as part of routine health examinations. Total
cytoplasmic RNA was prepared from PMBC isolated from human and
chimpanzee donors. In addition, RNA was extracted from NK cell clones
NK1.3 and M1.1 established from chimpanzees Cathy and Mason,
respectively (35); a splenic sample from the chimpanzee
Edwina was also used as source material. Genomic DNA of higher primates
used for PCR typing analysis was isolated from previously established
EBV-transformed B cell lines (36, 37). These include
samples from the family Hominidae: humans (Homo sapiens,
n = 11), common chimpanzees (Pan
troglodytes, n = 48), pygmy chimpanzees or bonobos
(P. paniscus, n = 11), gorillas
(Gorilla gorilla, n = 2), and orangutans
(Pongo pygmaeus, n = 5); and from the family
Hylobatidae: common gibbons (Hylobates lar,
n = 3).
|
Total cytoplasmic RNA was isolated with the reagent RNAzol following manufacturers protocol (Tel-Test, Friendswood, TX). First-strand cDNA was prepared from total RNA using murine Moloney leukemia virus-reverse transcriptase (Life Technologies, Rockville, MD) by standard procedures (38).
cDNA was amplified by PCR in 1x AmpliTaq buffer (Applied Biosystems, Foster City, CA), 0.25 mM of each dNTP, 0.8 pM each of a sense (forward) and antisense (reverse) primer, and 1 U of AmpliTaq DNA polymerase (Applied Biosystems). The following PCR conditions were used: initial denaturation at 96°C for 1 min; then 30 cycles of denaturation at 94°C for 30 s, primer annealing at 5260°C for 30 s, and extension at 72°C for 1 min; followed by 72°C for 10 min to promote complete extension.
Primers were designed from available human cDNA sequences in the
GenBank database. Initially, certain primer sets were designed to
amplify more than one NKG2 gene; gene-specific primer sets
were used subsequently. The following oligonucleotide primers were used
for our cDNA amplifications: CD94-PCRF and CD94-PCRR (for
CD94 and Pt-CD94), NKG2ABC-PCRF and NKG2ABC-PCRR
(for NKG2A/B, C, and Pt-NKG2A/B, C); hNKG2A-PCRF
and hNKG2AC-R1 (for NKG2A/B and Pt-NKG2A/B);
NKG2C-PCRF and hNKG2AC-PR1 (for NKG2C); NKG2C-F1 and
NKG2C1-R1 (for Pt-NKG2CI); NKG2C-F1 and NKG2CN-R1 (for
Pt-NKG2CII); DFORa and DREV (for NKG2D and
Pt-NKG2D); ECF-F1 and E-R1 (for NKG2E); ECF-F4
and E-R4 (for Pt-NKG2E); ECF-E1F1 and E-E7R1 (for
Pt-NKG2E partial); and NKG2F-F1 and NKG2F-R1 (for
NKG2F and Pt-NKG2F) (Table I
).
|
DNA sequences were assembled and analyzed using the computer program AutoAssembler (version 2.1; Applied Biosystems) and the Wisconsin Package sequence analysis software (version 10.1; Genetics Computer Group, Madison, WI).
Nucleotide sequences
All common chimpanzee CD94 and NKG2 genes
were assigned the prefix "Pt- " for P.
troglodytes. New nucleotide sequences determined here have been
deposited into the GenBank database under the accession numbers
AF25906163, AF2601356, and AF35000519 (Table II
). Sequences previously reported by us
have the accession numbers AF25905460 and AF260134 (35)
(Table II
). Sequences isolated from this study and reported by others
are also shown in Table II
; the "NM" prefix in the accession number
denotes GenBank reference sequences. Numerous alternatively spliced
variants were also identified in our analysis; these have not been
deposited into the GenBank database but are available from B. P.
Shum by request.
|
To examine whether alternative splicing could generate the
chimpanzee CD94 and CD94B variants, which differ
by insertion of three nucleotides (encoding a glutamine) between exon 4
and exon 5 in the cDNA, we determined the DNA sequence at the intron
4/exon 5 boundary of the Pt-CD94 gene. Genomic DNA from the
common chimpanzee Cathy was amplified by PCR with oligonucleotide
primers CD94-E4F2 and CD94-E5R2 (Table I
). The 50-µl reaction
contained 1x AmpliTaq buffer (Applied Biosystems), 0.25 mM of each
dNTP, 0.1 µM of each primer, 2.5 U of AmpliTaq DNA polymerase
(Applied Biosystems), and 200 ng of genomic DNA. PCR was performed with
an initial denaturation at 96°C for 2 min; 30 cycles of denaturation
at 94°C for 30 s, annealing at 60°C for 30 s, extension
at 72°C for 2 min; and a final extension of 72°C for 10 min. The
PCR product (
1969 bp) was cloned into the pCR-TOPO vector
(Invitrogen). Fifteen plasmid clones were isolated and their nucleotide
sequence was determined using T3 and T7 universal primers.
Transcription of the common chimpanzee Pt-NKG2CI and CII genes
To examine the transcription of Pt-NKG2CI and
Pt-NKG2CII genes in NK cells, we amplified cDNA by RT-PCR
with gene-specific primers from individual NK cell clones; these NK
clones established from the common chimpanzee Cathy have been described
previously (35). Each cDNA typing reaction of 25-µl
volume contained 1x AmpliTaq buffer (Applied Biosystems),
0.25 mM of each dNTP, 0.8 µM of each primer, and 1 U of AmpliTaq DNA
polymerase (Applied Biosystems). PCRs were performed with an initial
denaturation at 96°C for 1 min; followed by 35 cycles of denaturation
at 94°C for 15 s, annealing at 52°C for 15 s, and
extension at 72°C for 30 s; then a final extension at 72°C for
10 min. Amplified product was visualized under UV light after
separation by electrophoresis on a 1.2% agarose E-gel following
manufacturers instructions (Invitrogen). Gene-specific
oligonucleotide primers used for cDNA typing were as follows: for
CD94 and Pt-CD94, 94-CyF1 and 94-E6R1; for
NKG2A and Pt-NKG2A, A-E2F1 and A-E5R1; for
NKG2C and Pt-NKG2CI, CI-E1F1 and CI-E6R1; for
Pt-NKG2CII, CII-E1F1 and CII-E6R1; and for NKG2D
and Pt-NKG2D, D-CyF1 and D-E10R1 (Table I
). Selected
amplification products were sequenced to confirm specificity of
amplifications.
Genomic DNA typing
The presence of CD94 and NKG2 genes was
typed by PCR amplification of genomic DNA isolated from a panel of
higher primates using gene-specific oligonucleotide primers. Because
these amplification primers were based on human and chimpanzee
sequences, negative results in other species should not be interpreted
to mean the absence of a gene. All typing experiments were repeated at
least once for confirmation. The genomic typing primers used were as
follows: for CD94 and Pt-CD94, 94g-E5F1 and
94g-E6R1; for NKG2A and Pt-NKG2A, Ag-E1F1a and
Ag-E2R1; for NKG2C, Pt-NKG2CI, and
CII, Cg-E1F1 and Cg-E4R1; for Pt-NKG2CI, cG2C1-F1
and cG2C1-R1; for Pt-NKG2CII, cG2C2-F1 and cG2C2-R1; for
NKG2D and Pt-NKG2D, Dg-E4F1 and Dg-E5R1; and for
2-microglobulin
(
2m) control,
2m-F1 and
2m-R1
(Table I
).
PCRs were performed in a volume of 25 µl containing 1x AmpliTaq
buffer (Applied Biosystems), 0.1 mM of each dNTP, 0.4 µM of each
primer, 2 U of AmpliTaq DNA polymerase (Applied Biosystems), and
100
ng of genomic DNA. The following conditions were used for
amplification: initial denaturation at 96°C for 2 min; 10 cycles of a
first PCR (denaturation at 94°C for 15 s, annealing at 65°C
for 30 s, extension at 72°C for 1.5 min); 27 cycles of a second
PCR (denaturation at 94°C for 15 s, annealing at 60°C for
30 s, extension at 72°C for 1 min); and a final extension at
72°C for 10 min. The amplified products were separated by
electrophoresis on a 1.5% agarose gel in 0.5x Tris buffer with EDTA
buffer and were visualized by ethidium bromide staining afterward.
Included as positive controls for specificity were
1 pg of plasmids
containing cDNAs of each human and chimpanzee gene. Selected
amplification products were sequenced to confirm specificity.
Phylogenetic analysis
Sequences were aligned with the Pileup program of the Wisconsin Package (version 10.1; Genetics Computer Group) and manually adjusted afterward. Phylogenetic analysis was performed using the PAUP* 4.0 software package (Phylogenetic Analysis Using Parsimony (*and Other Methods), version 4.0b4a; Sinauer Associates, Sutherland, MA). The neighbor-joining method was used to construct dendrograms. Confidence of individual nodes was evaluated by 1000 bootstrap replications, and majority-rule consensus trees were generated.
| Results |
|---|
|
|
|---|
To investigate the polymorphism and interspecies divergence of
CD94 and NKG2 genes, we determined their coding
region cDNA sequences from four common chimpanzees and from healthy
human donors. The panel of 14 human donors represented major ethnic
groups, but not all of the genes were characterized for all donors
(Fig. 1
). No allelic polymorphism was
found for CD94 in either human or common chimpanzee, whereas
polymorphism was seen for all NKG2 genes in either one or
both species. Many alternatively spliced variants emerged from our
analysis, including the previously described CD94/B and
NKG2A/B (GenBank accession no. AJ000001) (7).
Alternatively spliced variants of several types were identified: ones
lacking exons, ones including introns, ones generated from cryptic
splice sites, and one having combinations of the above. Assignments for
cryptic splice sites were made from comparison of the variant cDNA
sequence with the sequence of full-length cDNA and/or corresponding
gene sequences. The presence of highly conserved "GT" or "AG"
dinucleotides at most of the alternatively spliced junctions revealed
the use of cryptic splice sites (data not shown) (39).
CD94
Data from six humans and four chimpanzees reveals no
CD94 polymorphism in either species (Fig. 1
). Human and
chimpanzee CD94 differ at only two nucleotide positions in
the coding sequence, one of which is nonsynonymous and the other
synonymous (Fig. 2
).
|
|
|
In summary, CD94 is encoded by a nonpolymorphic gene that is highly conserved in human and chimpanzee, having a nucleotide sequence similarity of 99.6%. Alternatively spliced CD94 variants were found in both species.
NKG2D
Three NKG2D alleles were characterized from the
analysis of 14 human donors (Fig. 1
). Each allele was identified from
more than one donor and at least five donors were NKG2D
heterozygotes (Fig. 1
A); NKG2D01 corresponds to
the sequence originally described by Houchins et al. (7).
The three NKG2D alleles are very similar, differing by
substitutions at only two nucleotide positions, of which one is
nonsynonymous and the other synonymous (Figs. 2
and 3
). A single
Pt-NKG2D sequence was encountered in the analysis of the
four chimpanzees, and it differs from the human NKG2D
alleles by 7- to 8-nt substitutions which produce 1- to 2-aa
differences (Figs. 1
B, 2, and 3).
Several alternatively spliced NKG2D variants were
characterized, two of which (variants 1 and 2) were found in both human
and common chimpanzee (Fig. 4
). Variants were found that included 106
bp of the 3' end of intron 3, intron 6, 21 bp at the 3' end of intron
8, or 23 bp within intron 9, or were missing exon 8. The additional
intron sequence in variant 1 is in the 5' untranslated region (3' end
of intron 3) and should not affect the protein product. The intron 6
inclusion (in variants 2, 3, and 4) results in a stop codon almost
immediately after the start of the intron 6 sequence and the encoded
proteins therefore lack an extracellular domain. Variant 5 has seven
additional amino acid residues in the carbohydrate recognition domain
(CRD) due to inclusion in the mRNA 21 bp from the 3' end of intron 8.
The additional 23 bp of intron 9 in variant 6 introduces a stop codon
that prevents the translation of exon 10.
In summary, we have defined three very similar human NKG2D
alleles and a single chimpanzee Pt-NKG2D allele. The
sequences have a nucleotide identity of
98.9%, demonstrating that
NKG2D, like CD94, is a highly conserved gene. A
number of alternatively spliced NKG2D variants were
identified in both human and chimpanzee.
NKG2A
Whereas a single NKG2A sequence was identified from
analysis of six human donors, three NKG2A alleles were found
in the four chimpanzees (Fig. 1
). The human sequence was identical to
that previously reported (7, 41). The three
Pt-NKG2A alleles differ at four nucleotide substitutions of
which two are nonsynonymous and two are synonymous (
Figs. 13![]()
![]()
). Human
and chimpanzee NKG2A differ by 710 nt substitutions in the
coding region and the encoded proteins differ by 46 amino acid
substitutions (Figs. 2
and 3
).
Previous analysis indicated that NKG2A and NKG2B
are alternatively spliced forms of the same gene, with NKG2A
including all exons in the message, whereas NKG2B message
lacks exon 4 (Fig. 4
and Refs. 7 and 41).
Clones corresponding to NKG2B were obtained from all six
human donors analyzed and from one of the four chimpanzees (Fig. 1
). An
additional splice variant found in chimpanzee Pt-NKG2A03
(variant 1) lacking exons 4 and 6 causes a frameshift and premature
termination in exon 7.
In summary, NKG2A is a highly conserved gene with
98.8%
nucleotide sequence identity in human and chimpanzee and share a
pattern of alternative mRNA splicing to give the NKG2B form.
Human NKG2A appears to have no polymorphism, whereas
chimpanzee Pt-NKG2A has modest polymorphism.
NKG2C
Analysis of nine human donors identified two NKG2C
alleles, both of which were carried by several individuals, two of whom
were shown to be NKG2C heterozygotes (Fig. 1
A).
The NKG2C01 allele has a sequence identical to that
deposited by Cantoni et al. (23) into GenBank
(accession no. Y13055). NKG2C02 differs from it by two
nonsynonymous substitutions at nucleotide positions 5 and 305 (Figs. 2
and 3
). The NKG2C sequence described by Houchins et al.
(7) (GenBank accession no. NM_002260) was not encountered
in our study and differs from NKG2C01 by one nonsynonymous
nucleotide substitution.
Six different Pt-NKG2C sequences were seen in the cDNA
clones isolated from the four common chimpanzees (Fig. 1
B).
Four of the sequences encode proteins that are similar in length and
sequence to human NKG2C (Pt-NKG2CI0104; Figs. 2
and 3
). Of
these the Pt-NKG2CI01, 02, and 04 proteins have two extra amino acid
residues in the carboxyl terminus when compared with the human NKG2C;
Pt-NKG2CI02 is five residues shorter due to an additional 21-bp
deletion in exon 1 that maintains the proper reading frame (Figs. 2
and 3
). The other two chimpanzee sequences (Pt-NKG2CII01 and
02) share a nucleotide deletion at position 358 that alters
the reading frame near the start of the region encoding the CRD.
Starting with position 140 (peptide alignment in Fig. 3
) the proteins
encoded by Pt-NKG2CII01 and 02 become totally
different from those encoded by Pt-NKG2CI and
NKG2C, and they terminate 35 amino acid residues following
the frameshift. Thus the predicted Pt-NKG2CII proteins have
cytoplasmic, transmembrane, and extracellular stalk regions like NKG2C,
but the CRD is replaced by a shorter region of unknown structure.
Pt-NKG2CII terminates at the same position as human and
chimpanzee NKG2F, but the two genes have different
frameshift mechanisms (Fig. 3
).
In comparison to the CD94, NKG2D, and
NKG2A genes, NKG2C was found to be most diverged
(Fig. 2
). Human NKG2C and the chimpanzee alleles that encode
similar protein products (Pt-NKG2CI) differ by 2224
nucleotide substitutions (96.696.9% sequence identities) that
produce 1417 amino acid differences (Figs. 3
and 5
A). Furthermore, that three
different Pt-NKG2C sequences were isolated from two
chimpanzees (Edwina and Elwood; Fig. 1
B) indicates that
there can be two Pt-NKG2C genes in common chimpanzee
genomes. In addition, clones corresponding to both Pt-NKG2CI
and CII sequences were isolated from all four chimpanzees
examined. Pairwise comparison of the coding region (Fig. 5
A)
and the 3' untranslated region (Fig. 5
B) of six
Pt-NKG2C sequences clearly shows that they divide into two
groups: one comprising the Pt-NKG2CI sequences, the other
comprising the Pt-NKG2CII sequences. Comparison between the
coding regions of the two common chimpanzee genes yields differences of
1116 nucleotide substitutions, not including gaps (Fig. 5
A). The range of nucleotide comparisons within each group
(NKG2C, Pt-NK2CI, and CII) is lower
and does not overlap with the range of differences for intergroup
comparisons (Fig. 5
). That both chimpanzee Pt-NKG2C genes
are equally divergent in nucleotide sequence from human
NKG2C is consistent with a history in which the
Pt-NKG2CI and CII genes arose by gene duplication
within the chimpanzee lineage after its split from the human
lineage.
|
Transcription of the Pt-NKG2CI and CII genes in a
panel of NK cell clones was examined. Seventeen NK cell clones derived
from the PBMC of chimpanzee Cathy (35) were typed by
RT-PCR for the presence of Pt-NKG2CI and
Pt-NKG2CII mRNA (Fig. 6
A). Six of the clones (35%)
typed positively for Pt-NKG2CI and three (18%) typed
positively for Pt-NKG2CII. One of the NK cell clones typed
positively for both Pt-NKG2CI and Pt-NKG2CII. The
negative typing reactions were not due to poor quality mRNA, as all the
NK cell clones typed positively for expression of Pt-NKG2D.
In conclusion, there is differential expression of the two
Pt-NKG2C genes in chimpanzee NK cell populations.
|
3.23.3% nucleotide differences). Within-species variation is
more striking for chimpanzee than for human, the former species having
two distinctive, transcribed NKG2C genes. NKG2E
Human NKG2E is closely related to NKG2C, but
its mRNA splices differentially into an additional exon 7 not used by
any other NKG2 gene (Fig. 4
). An alternatively spliced
NKG2H form (more similar in its splicing pattern to
NKG2C) has also been described (Figs. 3
and 4
and Ref.
42). The oligonucleotide primers we used for amplification
were designed to detect the NKG2E form, but not the
NKG2H form. Three NKG2E alleles were identified
from nine human donors analyzed, with at least two donors being
heterozygotes (Fig. 1
A). The human NKG2E alleles
vary at only two nucleotide positions, of which one is nonsynonymous
and the other synonymous (Figs. 2
and 3
). NKG2E01 is
identical to the GenBank reference sequence NM 002261
(43).
The oligonucleotide primers used to amplify human NKG2E
failed to amplify from chimpanzee cDNA. Successful amplification of
chimpanzee Pt-NKG2E was eventually achieved with two other
primer sets, one that amplified a segment of 552 nt corresponding to
about three quarters of the coding region, and another that amplified
the entire coding region. Clones corresponding to Pt-NKG2E
were obtained from three of the four chimpanzees tested. Most of the
clones corresponded to alternatively spliced variants or partial
Pt-NKG2E sequences; only one full-length clone was obtained.
The divergence between the human and chimpanzee NKG2E
nucleotide sequences (
3.4%) is comparable to that seen for
NKG2C. The nucleotide sequences of Pt-NKG2E
clones analyzed provided evidence for five alleles, distinguished by
the patterns of nucleotide substitution at six positions. Four of these
positions involve nonsynonymous substitutions and two involve
synonymous substitutions. All three chimpanzees from whom
Pt-NKG2E sequences were obtained appeared to be
heterozygotes; repeated attempts to amplify Pt-NKG2E from
chimpanzee Elwood failed.
Two alternatively spliced variants were identified from the human
NKG2E (Fig. 4
, variants 1 and 2) which were distinct from
the eight variants of Pt-NKG2E (variants 310), but all
contain a change in the normal reading frame resulting in premature
termination. Whereas all variants involved "conventional" splicing
differences within the exons and introns of Pt-NKG2E,
variant 10 (Fig. 4
) consists of four exons from Pt-NKG2CII02
(though missing the 5' third of exon 4), followed by 50 bp identical to
a part of the NKG2E promoter (also identical to part of the
NKG2C promoter), and then continues with exons 3 through 6
of Pt-NKG2E03 (but missing the 3' part of exon 4). This
suggests that the duplicated Pt-NKG2CII gene is directly
upstream of Pt-NKG2E.
Of the NKC genes compared in chimpanzee and human,
NKG2E and NKG2C show the highest interspecies
divergence (
3.23.4% difference in nucleotide sequence) and they
are both polymorphic within each species. Alternatively spliced
variants of NKG2E were frequently encountered, all of them
containing frameshifts that cause premature termination.
NKG2F
Four NKG2F alleles were obtained from the analysis of
five human donors and five alleles from the four chimpanzees (Fig. 1
),
resulting in three NKG2F and five Pt-NKG2F allotypes (Figs. 2
and 3
).
The human NKG2F01 allele is identical to the original
sequence from Plougastel and Trowsdale (44) (GenBank
accession no. U96845) and NKG2F02 is identical to the
GenBank reference sequence NM 013431 originally described by Glienke
et al. (10). At least two humans and three chimpanzees
appear to be NKG2F heterozygotes (Fig. 1
). Clones
corresponding to three NKG2F alleles were obtained from
Elwood, the chimpanzee from whom NKG2E could not be
amplified. The alleles Pt-NKG2F04 and 05 differ
from the others by a 3-bp deletion in the 5' untranslated region and a
4-bp deletion in the 3' untranslated region (data not shown).
Human and chimpanzee NKG2F have
98.6% sequence
similarity, differing by 5- to 9-bp differences between the two
species. Striking is that only 2 of the 16 polymorphic positions are
species specific (Fig. 2
).
Three alternatively spliced variants of Pt-NKG2F01 were encountered, all involving exon 4 and the use of cryptic splice sites. Variant 1 lacked 183 bp within exon 4, but the deletion was to the 3' side of the NKG2F stop codon. Variant 2 lacked 68 bp at the 5' end of exon 4, which changes the reading frame to the one used by NKG2A, C, and E in this exon. Variant 3 lacked the same internal exon 4 sequence as variant 1 but also lacked 49 bp from the 5' end of the exon and is prematurely terminated.
In summary, the NKG2F gene has a pattern of variation that distinguishes it from the other NKC genes examined. It combines relative conservation between the two species with greater polymorphism within each species.
Haplotypic variation in the chimpanzee NKC
Based upon the sequences of the human and chimpanzee
CD94 and NKG2A, C, and D
genes, a system of PCR sequence-specific primers typing of genomic DNA
was developed and used to analyze samples obtained from 11 human donors
and 69 apes representing six taxonomic species in total. Typing for the
2m gene was included as a
positive control (Fig. 6
B). Typing for NKG2C
involved three reactions: one aimed at a generic typing for
(Pt-)NKG2C, one targeted at
Pt-NKG2CI, and one targeted at Pt-NKG2CII.
All 80 samples typed positively for
2m and NKG2A, and all
but the three gibbons typed positively for CD94. The 11
humans gave identical typing reactions in which the only negatives were
for the two Pt-NKG2C genes. In contrast, a heterogeneity of
typing pattern was seen within the panel of 48 common chimpanzees. The
most frequent pattern, seen in 30 individuals, comprised positive
typing for all the genes tested. Six other patterns, all seen in at
least two individuals, were distinguished by negative reactions for
NKG2C and NKG2D genes. Most divergent from the
dominant pattern 1 was pattern 2, in which typing for NKG2D
and all three NKG2C reactions was negative (Fig. 6
B). This pattern of reactivity was also given by the two
gorillas tested.
The 11 pygmy chimpanzees tested gave an identical pattern which distinguished them from all individuals of the other five species. This pattern was comprised of positive reactions for all the genes tested, with the exception of Pt-NKG2CI. Species-specific patterns of typing were also seen for the five orangutans and three gibbons tested. Whereas the gibbons were negative for NKG2C and NKG2D (as well as CD94), the orangutans were negative for NKG2C but positive for NKG2D.
It must be emphasized that negative typing reactions do not show that a gene is absent, merely that at least one of the target sequences used for PCR is not present. Thus the results of the typing provide a measure of divergence from the human and chimpanzee sequences used to design the typing system. Most striking is the phylogenetic conservation of NKG2A and CD94 compared with NKG2C and NKG2D, a result that is consistent with and complementary to the comparison of cDNA sequences in human and common chimpanzee. Also impressive is the variability of the NKC genes within the common chimpanzee. The actual heterogeneity and polymorphism of the chimpanzee NKC is probably much greater than is apparent from this analysis because only individuals who are negative on both haplotypes will give the negative reactions needed to be distinguished from pattern 1.
| Discussion |
|---|
|
|
|---|
|
In parallel with the analysis of human diversity in NKC
genes we performed a similar analysis in common chimpanzee. For each of
the human genes we identified chimpanzee counterparts that are >96%
identical in nucleotide sequence. Thus, as a group these NKC
lectin-like receptors of NK cells are more conserved than the
KIR genes encoding Ig-like receptors (Fig. 8
). With the exception of
NKG2C, each of the human NKC genes corresponds to
a single chimpanzee NKC gene which we interpret to be its
ortholog. In contrast, the single human NKG2C gene is
equally related to two chimpanzee paralogs, one of which
(Pt-NKG2CI) encodes a protein like human NKG2C, whereas the
other (Pt-NKG2CII) has a frameshift that produces a shorter
protein which lacks an extracellular CRD. Both Pt-NKG2C
genes were demonstrated to be independently and clonally expressed in
chimpanzee NK cell clones (Fig. 6
A). More polymorphism was
observed in the NKC genes of four chimpanzees than was
observed in the larger panel of human donors. However, the numbers of
allotypes were few, as were the number of amino acid residues that
distinguished them. Thus in both human and chimpanzee the
NKC genes are conserved.
|
Vance et al. (53) proposed two models for the evolution of the genes named NKG2A, C, and E in humans and mice: in the first model they are not orthologs and derived from lineage-specific pathways of gene duplication, deletion, and divergence; in the second model they are orthologs but have acquired strong species-specific characteristics by concerted evolution through gene conversions. In either case the process led to order-specific and species-specific divergence. Distinguishing the two alternatives is difficult, but the latter model was favored by Vance et al. (53) for the following reasons. First, they considered it unlikely to be coincidental for both activating and inhibitory genes to be inherited in both mammalian lineages; second, the overall gene orders of NKG2A, C, and E are preserved in the human and mouse NKC (12, 53); and third, the 5' ends of NKG2A (inhibitory form) and NKG2C (activating form) are divergent (52, 53), indicating that the two genes have been separated for a long time. Providing some contrast to this view is our comparison of human and chimpanzee NKC genes showing how genes within the third lineage can evolve by species-specific duplications over short periods of evolutionary time. Thus the two chimpanzee Pt-NKG2C genes are both paralogs of the human NKG2C and neither is an ortholog. Moreover, the apparent haplotypic diversity of the Pt-NKG2C genes seen in common chimpanzees indicates that this subregion of the NKC has repeatedly been a target for evolutionary change. Thus a model in which the human and mouse NKG2A, C, and E genes are paralogs rather than orthologs cannot be ruled out.
Further evidence for species-specific evolution within the
NKC is seen in a tree made from just the primate
NKG2A, C, E, F, and related
sequences (Fig. 7
B). This tree contains three
distinguishable groups of sequences, corresponding to the
NKG2A sequences, the NKG2F sequences, and the
combination of NKG2C and NKG2E sequences. The two
rhesus monkey (Macaca mulatta) sequences in the
NKG2A group (48) are considerably diverged,
both from human and chimpanzee NKG2A and from themselves (11
amino acid differences). Either there is a polymorphism greater than is
seen in human and chimpanzee at a single Mm-NKG2A locus or
there are two Mm-NKG2A genes. It is also intriguing that the
rhesus monkey MHC-E is oligomorphic or polymorphic
(59). Similarly, within the NKG2F group there
are two divergent rhesus monkey NKG2F-like sequences
(48) of which one is closer to human and chimpanzee
NKG2F. Within the NKG2C group, the eight
Mm-NKG2C sequences (48) form a species-specific
subgroup that does not interleaf with the human and chimpanzee
NKG2C sequences. All these data point to the NKC
of the rhesus monkey having significant differences in gene content and
haplotypic variation from the NKCs of human and
chimpanzee.
Although the CD94 and NKG2 genes of the
NKC are generally older, less variable, and slower evolving
than the KIR gene family of the LRC (Fig. 8
),
phylogenetic comparisons suggest that they have nontrivially diverged
during the course of mammalian evolution. That is clearly the case for
the ligands recognized by these lectin-like receptors. Thus the Qa-1
class I ligand for mouse CD94:NKG2A and CD94:NKG2C receptors does not
in its sequence have particular affinity with the HLA-E class I ligand
for human CD94:NKG2A and CD94:NKG2C receptors (52, 53, 60, 61). Even more different are the divergent MHC class I-related
ligands for NKG2D, which comprise MIC (24, 25) and ULBPs
(47) in humans, and the retinoic acid early inducible
RAE-1 and H60 in mice (62, 63).
Comparison of the human and chimpanzee NKC genes indicates
that they have evolved under different selection pressures:
CD94 is most conserved, followed by NKG2D,
NKG2A, and NKG2F, with NKG2C and
NKG2E being most divergent (Fig. 8
) (35). The
human and common chimpanzee CD94 have only 0.4% nucleotide
difference, well below the average
1.24% variation as suggested for
intergenic regions of human and chimpanzee genes (64).
NKG2A and NKG2D also appear to be under purifying
selection. For CD94 and NKG2A, the identical
inhibitory MHC class I specificity of the human and chimpanzee
CD94:NKG2A receptors supports this mode of selection (35).
The human and common chimpanzee NKG2C and NKG2E
genes differ by 3.23.4% in nucleotide sequences, values higher than
expected for the divergence of human and chimpanzee genes under neutral
selection, and similar to those KIR genes which are
orthologous in chimpanzee and human (Fig. 8
). Most nucleotide
substitutions that distinguish NKG2C and
Pt-NKG2CI change residues in the CRD loops (Figs. 2
and 3
)
predicted to interact with ligand (65). These differences
might affect the ligand-binding specificity and be a consequence of
selection upon activating receptors by pathogens. Such selection might
also be the cause of the greater haplotypic diversity of
NKG2C compared with other NKC genes, as
seen both within species (as in the common chimpanzee) and between
species.
The rodent NKC contains a diverse family of
Ly49 genes that encode lectin-like receptors specific for
polymorphic determinants of classical MHC class I molecules
(66). In humans the Ly49 family is represented by a single
nonfunctional gene, Ly49L (67), while KIR
provide functions analogous to those of Ly49 in mice. In our
investigation we did not examine the structure and diversity of the
Ly49L gene. However, Mager et al. (68) have
recently shown that a gene related to human Ly49L is present
in common chimpanzee as well as in other primates (gorilla, orangutan,
gibbon, baboon, and African green monkey) and other mammalian orders.
In contrast, genes related to the expanded rodent Ly49
family were not detectable by Southern blot analysis in non-rodent
species. Partial sequence analysis showed that common chimpanzee
Ly49L has the same inactivating mutations as the human gene,
whereas that was not so for the other primates. The extent to which
these differences are fixed in humans, chimpanzees, and other species
was not addressed, as analysis was confined to individuals of a
species. For baboon, a complete cDNA sequence for Ly49L was
determined (
95% nucleotide identity with the human Ly49L
in the coding region) and it was shown to be expressed in baboon
lymphocytes along with transcripts from several KIR genes.
These observations raise the possibility that Ly49L encodes
a functional receptor in some primate species and that its inactivation
in humans and chimpanzee may be exceptional. Also worthy of
consideration is that inactivation of the Ly49L gene may not
be a property of all human and chimpanzee NKC
haplotypes.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Current address: Department of Cell and Molecular Medicine, Southampton General Hospital, Southampton, U.K. ![]()
3 Current address: Institute for Transplantation Diagnostics and Cell Therapeutics, Heinrich Heine University Medical Center, Dusseldorf, Germany. ![]()
4 Address correspondence to Dr. Peter Parham, Department of Structural Biology, Stanford University School of Medicine, 299 Campus Drive West, Sherman Fairchild Building, D-157, Stanford, CA 94305-5126. E-mail address: peropa{at}leland.stanford.edu ![]()
5 Abbreviations used in this paper: KIR, killer cell Ig-like receptor;
2m,
2-microglobulin; CRD, carbohydrate recognition domain; LRC, leukocyte receptor complex or cluster; NKC, NK complex; ULBP, unique long binding protein. ![]()
Received for publication May 4, 2001. Accepted for publication October 22, 2001.
| References |
|---|
|
|
|---|

T cells by NKG2D via engagement by MIC induced on virus-infected cells. Nat. Immunol. 2:255.[Medline]
This article has been cited by other articles:
![]() |
E. J. Petrie, C. S. Clements, J. Lin, L. C. Sullivan, D. Johnson, T. Huyton, A. Heroux, H. L. Hoare, T. Beddoe, H. H. Reid, et al. CD94-NKG2A recognition of human leukocyte antigen (HLA)-E bound to an HLA class I leader sequence J. Exp. Med., March 17, 2008; 205(3): 725 - 735. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Yu, G. Heller, J. Chewning, S. Kim, W. M. Yokoyama, and K. C. Hsu Hierarchy of the Human Natural Killer Cell Response Is Determined by Class and Quantity of Inhibitory Receptors for Self-HLA-B and HLA-C Ligands J. Immunol., November 1, 2007; 179(9): 5977 - 5989. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Averdam, H. Kuhl, M. Sontag, T. Becker, A. L. Hughes, R. Reinhardt, and L. Walter Genomics and Diversity of the Common Marmoset Monkey NK Complex J. Immunol., June 1, 2007; 178(11): 7151 - 7161. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Karimi, T. M. Cao, J. A. Baker, M. R. Verneris, L. Soares, and R. S. Negrin Silencing Human NKG2D, DAP10, and DAP12 Reduces Cytotoxicity of Activated CD8+ T Cells and NK Cells J. Immunol., December 15, 2005; 175(12): 7819 - 7828. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. P. Williams, A. R. Bateman, and S. I. Khakoo HANGING IN THE BALANCE: KIR and Their Role in Disease Mol. Interv., August 1, 2005; 5(4): 226 - 240. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. G. de Groot, C. A. Garcia, E. J. Verschoor, G. G. M. Doxiadis, S. G. E. Marsh, N. Otting, and R. E. Bontrop Reduced MIC Gene Repertoire Variation in West African Chimpanzees as Compared to Humans Mol. Biol. Evol., June 1, 2005; 22(6): 1375 - 1385. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Han, M. Zhang, N. Li, T. Chen, Y. Zhang, T. Wan, and X. Cao KLRL1, a novel killer cell lectinlike receptor, inhibits natural killer cell cytotoxicity Blood, November 1, 2004; 104(9): 2858 - 2866. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Miyashita, N. Tsuchiya, K. Hikami, K. Kuroki, T. Fukazawa, M. Bijl, C. G. M. Kallenberg, H. Hashimoto, T. Yabe, and K. Tokunaga Molecular genetic analyses of human NKG2C (KLRC2) gene deletion Int. Immunol., January 1, 2004; 16(1): 163 - 168. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A. Rabinovich, J. Li, J. Shannon, R. Hurren, J. Chalupny, D. Cosman, and R. G. Miller Activated, But Not Resting, T Cells Can Be Recognized and Killed by Syngeneic NK Cells J. Immunol., April 1, 2003; 170(7): 3572 - 3576. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. A. Guethlein, L. R. Flodin, E. J. Adams, and P. Parham NK Cell Receptors of the Orangutan (Pongo pygmaeus): A Pivotal Species for Tracking the Coevolution of Killer Cell Ig-Like Receptors with MHC-C J. Immunol., July 1, 2002; 169(1): 220 - 229. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. G. Shilling, N. Young, L. A. Guethlein, N. W. Cheng, C. M. Gardiner, D. Tyan, and P. Parham Genetic Control of Human NK Cell Repertoire J. Immunol., July 1, 2002; 169(1): 239 - 247. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |