|
|
||||||||



*
Department of Microbiology and Immunology, University of Miami School of Medicine, Miami, FL 33136;
Department of Microbiology and Immunology, University of Maryland, Baltimore, MD 21201;
Basel Institute for Immunology, Basel, Switzerland; and
§
Department of Biosystems Science, Graduate University for Advanced Studies, Hayama, Japan
| Abstract |
|---|
|
|
|---|
2
peptide-binding domains, and there are highly diverse alleles (defined
by the
1 peptide-binding domain) in each lineage. Surprisingly, the
3 domains are homogenized in both lineages, suggesting that
interallelic gene conversion/recombination maintains the high sequence
similarity. | Introduction |
|---|
|
|
|---|
The Xenopus MHC, however, has two characteristics that make the model unique. First, MHC genes are silenced in polyploid individuals to a diploid number, even in higher order polyploids that potentially harbor eight copies of the MHC (9, 10). This limiting of MHC gene number may place unique pressures on the frog immune system, reflected in the number and diversity of all types of MHC genes. Second, immunocompetent tadpoles express class Ia proteins and the immune proteasome element lmp7 only at low levels until metamorphosis when the genes become expressed ubiquitously (11, 12), while class II molecules are expressed in high amounts by B cells and macrophages (13, 14). After metamorphosis, class II is expressed by all lymphocytes and APC. This differential expression of MHC during development may be important for limiting autoaggressive responses when new adult-type proteins are expressed at metamorphosis (12, 15).
Concerning the general evolution of the class I gene family, there are several other unusual features of Xenopus class Ia genes and proteins: 1) unique among all vertebrates studied to date, the MHC proper encodes only one class I gene, of the class Ia type, as detected by Southern blotting under low stringency conditions (1); 2) the single serologically detected class I molecule has unusual biochemical properties, including an association with non-ß2-microglobulin proteins on erythrocytes (16), allelic size heterogeneity, and differential affinity of association with ß2-microglobulin among allelic products (17); and 3) the class I gene appears to be sandwiched between the putative class II and class III regions, probably in close linkage to the TAP2 and lmp7 genes (6, 7). To gain more insight into the Xenopus model, and to better understand the evolutionary history of the MHC, we embarked on further analyses of the frog class I genes. Specifically, we asked whether the sequences of some class Ia alleles might reveal the basis for their unusual biochemical and genetic characteristics.
| Materials and Methods |
|---|
|
|
|---|
Xenopus laevis with the MHC haplotypes f, g (18), j (19), or r (17), and X. laevis-Xenopus gilli (LG) interspecies hybrid frogs with the haplotypes a/c (LG15) and b/d (LG3) (20) were used in all of these studies. X. laevis and X. gilli are both pseudotetraploid frogs (2n = 36) that are capable of cross-fertilization and that are believed to have diverged from a common ancestor between 1 and 10 million years ago (20). Outbred Rana pipiens purchased from Carolina Biological Supply (Burlington, NC) were used to prepare the cDNA library for isolation of class I cDNA clones.
Immunoprecipitation
The cell surface labeling of erythrocyte proteins, preparation
of lysates, and immunoprecipitation procedures have been described
previously in detail (17). The gel in Fig. 1
displays
radioiodinated cell surface class I proteins immunoprecipitated with
Xenopus alloantisera. The proteins were first
electrophoresed under nonreducing conditions, and then excised from the
gel and run on another gel under reducing conditions (see legend).
|
The original class Ia clone, encoded from the f
haplotype, was published (1). cDNA libraries were prepared
from LG3 (b/d) liver and spleen RNA and LG15
(a/c) intestine and spleen RNA, as was done
previously for f/f (21). The libraries were
screened under high stringency conditions (1) with the
presumed exon encoding the
3 domain. Only one unique clone was
isolated from the LG3 library (b/d-2 in Fig. 3
). a/c-1, a/c-2, and
a/c-2' were isolated from the LG-15 library (Figs. 3
and 4
). The other
class I cDNAs were cloned by PCR using spleen first strand cDNA as
template (2). r/r was amplified with primers in the 5' and
3' untranslated regions ending on the start and stop codons,
respectively. j/j and g/g were cloned by 3' RACE (2, 22);
the 5' untranslated primer ending on the start codon and an adapter
primer (2, 22) at the 3' end were used in the
amplifications. The 5'UT primer was not successful with b/d-1 (LG3)
cDNA; thus, an internal primer at the beginning of the
1 domain
(1) and the 3' primer ending on the stop codon had to be
used to obtain a PCR product. Thus, the a/c-2, b/d-2, g/g, and j/j
class I cDNA clones are
1.8 kb, a/c-2' is 2.7 kb (see text), a/c-1
is 3 kb, and the b/d-1 and r/r clones encompass only the coding regions
(the mRNA encoding the r/r clone is 3 kb, see Fig. 2
; b/d-1 was not
selected from the cDNA library because a relatively low number of
primary recombinants was obtained from this library and (also) because
of the large size of the mRNA). The PCR products were cloned into the
EcoRV site of pBluescript, and the coding regions were
sequenced on both strands with universal (T3 and T7) and gene-specific
primers (23). Sequences were manually aligned;
dS and
dN were calculated by the method of
Nei and Gojobori (24).
|
|
|
A cDNA library was prepared from the poly(A)+ RNA from pooled spleens of eight outbred R. pipiens adults (21). Screening was performed under low stringency conditions with a probe encoding all three external domains of the X. laevis f/f class I cDNA (1). The coding sequences of two clones were sequenced in their entirety.
Northern blotting
Fifteen micrograms of total splenic RNA from animals homozygous
for the f, g, j, and r
haplotypes and from LG-15 (a/c) were
electrophoresed and blotted, as described (12).
Hybridizations were performed under high stringency conditions with the
exon encoding the
3 domains as probe, as it was assumed (correctly,
in retrospect) to be the most conserved region of class I
mRNA.
| Results |
|---|
|
|
|---|
Previously, we have shown that class I proteins encoded by the
a, b, c, d, f, g, j, and r haplotypes
migrate to different positions by SDS-PAGE (radioiodinated
r, f, g, and j class I molecules in Fig. 1
; also see Refs. 16, 17). f, r, and one allele from the LG15 and LG3 hybrids (a and
b) migrate to a relatively higher position under reducing conditions
than g, j, and the other allele from the LG hydrids (c and d). In
previous experiments with defined crosses of Xenopus, all of
these alleles have been shown to segregate as alleles at a single locus
(1, 7, 16, 17). Such differences in the migration of
allelic sets on SDS gels are attributable to some other biochemical
property(ies) besides glycosylation since the proteins, after
endoglycosidase-F digestion, still display the same m.w. differences
(17). After the cloning of the class Ia gene from the
f haplotype (1), we embarked on sequence
analysis of all eight alleles to examine: 1) the biochemical basis for
these two sizes of class I proteins; 2) the previously described
unusual serological properties; and 3) their level of diversity.
Northern blot analysis
The original class Ia cDNA clone encoded by the f
haploytpe is
3 kb (1), about the size of f class Ia
mRNA, as detected by Northern blotting (Fig. 2
A). Surprisingly, while r
class I mRNA is also of this high m.w., g and j display smaller-sized
mRNAs (
1.8 kb). Note that the very weak band in the f and r lanes
that is approximately the same size as the g and j bands is not a class
I transcript, but the leading edge of 18S ribosomal RNA; no low m.w.
transcript for class I has been cloned from either the f or
r haplotypes, and a repeat of the northern with
f/f RNA reveals no band at
1.8 kb (Fig. 2
B). Surprisingly, for the LG15 (a/c)
interspecies hybrid, three bands are always detected by Northern
blotting, one coincident with the f (3-kb) band and one with the g
(1.8-kb) band, and a third band at
2.7 kb (upper band marked with an
arrow in the class I-2 bracket in Fig. 2
B). This upper class
I-2 band is also found in RNA isolated from g/g
erythrocytes, but not from intestine (compare Fig. 2
, A and
B). Sequence of the 3'UT of the LG-15 2.7-kb band shows it
to be an alternatively spliced form of the lower m.w. transcript
(a/c-2' in Fig. 4
); this species is found
in high amounts in erythrocyte RNA for all class I-2 haplotypes (e.g.,
j and g) and in all tissues for one LG15
haplotype (Fig. 2
B and data not shown). In total, the
Northern blotting data suggest that there are distinct classes of mRNA
for the Xenopus alleles (3 kb for f and r, and 1.8 (and 2.7)
kb for g and j), and further that the X. gilli and X.
laevis alleles also conform to these two classes; this is perhaps
not surprising considering the relatively recent divergence time for
X. laevis and X. gilli. The groupings detected by
Northern blotting correspond precisely to the f, g, j, and r (and LG)
protein data described above (without protein sequences or segregation
analyses of the LG alleles, however, we have not shown definitively
that all of the lower m.w. mRNAs encode the lower m.w. proteins
detected under reducing conditions).
Sequencing of alleles
The coding regions of the eight alleles were sequenced, and the
deduced amino acid sequences were aligned. There is great diversity in
the PBD3 and,
conversely, near invariance in the Ig-like
3 domain (Fig. 3
). Surprisingly, the cytoplasmic, and to
a lesser extent the transmembrane regions are divergent, and the
alleles segregate into the same two groups resolved by the Northern
blotting: f/f, r/r, a/c-1, and b/d-1 fall into one group and g/g, j/j,
a/c-2, and b/d-2 into the other. This segregation into two classes is
also apparent in the majority of the
2 PBD, but is not as noticeable
in the
1 PBD; indeed, the
1 domain is much more diverse than the
2 domain, and even within allelic sets there is a high degree of
diversity (Figs. 3
and 5
).
|
1 and
2 PBD
trees were made with and without amino acids that bind peptides and are
polymorphic in other species; because these residues are very diverse
even within a species, it was logical to remove them from the analysis
when comparing sequences over extremely long distances. In this way we
hoped to reveal true phylogenetic distances in the absence of positive
selection (25, 26). The distances between the
Xenopus alleles in the
1 domain are similar to the
distance between human and mouse genes. In contrast, the distances
between the mouse and human sequences in the
3 domains are nearly
the same as that in the other two domains, while the Xenopus
3 sequences cluster very tightly.
The frequency of mutations in the PBD, at least for the
1 domains,
is higher than any set of alleles characterized to date (Table I
and Fig. 6
; see Ref. 25). However,
while the number of nonsynonymous (dN)
changes is very high, their frequency is the same as (
1) or not
significantly greater than (
2) the synonymous substitutions
(dS). Surprisingly, given the great
distance between alleles in the
1 domain, and especially the
cytoplasmic domains, the
3 domains are extremely well conserved,
suffering neither synonymous nor nonsynonymous changes.
|
|
R. pipiens sequences
Because of the interesting features of the Xenopus MHC,
we thought it appropriate to begin to analyze MHC of other amphibians
to determine whether the same basic phenomena occur in evolution. In
addition, we hoped to provide stability to our phylogenetic trees by
including sequences from other anuran amphibians. R. pipiens
MHC class I cDNAs were obtained by cross-hybridization with the
Xenopus class Ia probe under low stringency conditions.
Sequences of the two clones show them to be allelic and most likely
class Ia clones: these sequences are more similar to the
Xenopus class Ia alleles in the PBD than the Ig-like
3
domain (Figs. 3
and 5
), suggesting, as was previously observed in
chickens and Xenopus, that the PBD are more conserved over
large phylogenetic distances (1, 26). In addition, amino
acid differences between the two presumed Rana alleles are
primarily concentrated in the amino acids that interact with peptide,
and the invariant class Ia amino acids that anchor the ends of peptides
in the groove are generally conserved or substituted for by amino acids
of a similar composition. The genetic distance between the
Rana and the Xenopus class Ia sequences is large,
especially in the
3 domain, consistent with an old divergence
between the two anuran families (at least 150 million years)
(27).
| Discussion |
|---|
|
|
|---|
It is clear from our work that different functional regions of the
Xenopus class Ia alleles evolve at different rates, a result
consistent with other evolutionary studies (e.g., see Ref.
26). Three findings particularly deserve attention: 1) the
class I domains that are most dissimilar are the cytoplasmic regions,
followed by
1,
2, and
3; 2) while there is great diversity in
the PBD among alleles, the number of nonsynonymous substitutions in the
codons that encode peptide-binding residues is not significantly higher
than synonymous ones; and 3) despite the high rate of divergence, there
are very few substitutions of any type in the
3 domain.
Although it is difficult to interpret all of these data, we can
understand some facets of the work. Comparisons among the
Xenopus alleles are more similar to data obtained when
sequences are compared between mouse and human class I genes; in fact,
the distances between the
1 domains (and the CYT regions) of the
Xenopus alleles, as mentioned, are similar to those between
mouse and human class I isotypes (Figs. 3
and 5
). Considering that mice
and humans diverged at least 6065 million years ago with some
estimates of up to 100 million years (28), the
Xenopus alleles may have emerged at around the same time as
the common ancestor of most MHC-linked mammalian class I genes.
If the alleles and the allelic lineages are so old, then why are the
3 domains of the Xenopus alleles so well conserved, to
the point of homogenization? One might argue that since all of the
class Ia allelic forms must associate with
ß2-microglobulin and perhaps other
frog-specific proteins, that the
3 domain must remain essentially
invariant. However, the distances between the mouse and human isotypes
are not that much smaller in
3 as compared with
1 and
2
domains, and class I proteins from these species presumably also have
similar requirements in their associations with other proteins.
Furthermore, this interpretation cannot account for the extremely low
number of synonymous changes detected, which are expected to be similar
to the frequency found in residues in the
1 and
2 domains. Thus,
it is difficult to even speculate at this time as to why the
Xenopus
3 domains are so remarkably similar. It appears
that the sequences are homogenized by interallelic
conversion/recombination, as has been shown for intronic regions of
mammalian class I genes (29). Such homogenization can also
be seen between the allelic lineages in some areas of the 3'UT regions
(Fig. 4
). Because Ig-like
3 domains are often used to detect
phylogenetic distances among class I genes, one obviously must be wary
of the significance of such trees.
The general divergence trend of CYT >
1 >
2 >
3 is also true of the Xenopus class Ib (XNC)
genes, suggesting that the pressures placed on class Ia alleles also
act on different class I genes (7). The differences in the
sequence diversity in the
1 and
2 PBD suggest that the former
domain is selected primarily for its diversity and the latter for
another function, perhaps for interactions with other components of the
class I-processing machinery. Furthermore, it should be pointed out
that although the sequences of CYT regions of the class Ia alleles are
very different between the two allelic classes, within each class the
transmembrane/CYT, like the
3 domains, are very conserved. The one
exception to this rule is that the g and j genes
have accumulated a few changes not found in b/d-2 and a/c-2; the most
likely explanation for this finding is that the two alleles from the LG
hybrids are derived from X. gilli and have been separated
from their counterparts in X. laevis, g and
j, for a longer time than any of the alleles from the other
(f, r) group, which are probably all derived from
X. laevis. Perhaps the cytoplasmic tails of the two lineages
are important for a particular function, e.g., interaction with
cytosolic proteins.
In mammalian class I molecules, there are significantly more
nonsynonymous changes over synonymous ones in the codons specifiying
residues that interact with peptides, supporting the idea of
overdominant selection on MHC proteins (25, 30). In the
Xenopus
1 domains, the frequency of mutation is 23
times that of the human and mouse, and it is likely that the synonymous
sites have become saturated. Thus, while
dN is not significantly greater than
dS, we think this is most likely due
to the age of the alleles, not to a lack of selection. More recently
derived Xenopus alleles must be examined to test this
notion.
Genes vs alleles
The cytoplasmic domains and 3'UT are quite different between the two sets of alleles, suggesting that these sets were originally different genes rather than alleles, although they can segregate as alleles at a single locus. The same lineages (f, r on the one hand, and g, j on the other) have been revealed for lmp7 (31) and TAP2 genes (7). The simplest interpretation is that there was a long separation of these lineages, followed by an introgression sometime in the history of Xenopus as a genus (32). Because the genus Xenopus has been estimated to be more than 100 x 106 yr old, such a finding might not be that surprising and introgression indeed has been documented in this taxon (33); speciation via allopolyploidization, postulated to be common in Xenopus, would favor introgression. However, our preliminary work has shown that the two lmp7 allelic lineages are maintained in approximately equal numbers in a large number of species examined (M. Nonaka, M. Flajnik, and L. Du Pasquier, unpublished), suggesting rather that the lineages may be under balancing selection. We hypothesize that the two allelic sets may function in presentation of very different sets of peptides. The closer genetic linkage of lmp7/TAP2/class Ia in Xenopus (and other vertebrates such as birds (34) and zebrafish (35)) as compared with mammals may have allowed coselection of certain allelic lineages that are perpetuated and maintained in the population (6). Future studies will be aimed at examining the presence of the allelic lineages among diverse Xenopus species, and further to analyze whether there is a functional pressure to maintain the class I processing and presenting genes in a tight linkage group.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Address correspondence and reprint requests to Dr. Martin F. Flajnik, Department of Microbiology and Immunology, University of Maryland, 655 West Baltimore Street, Room 3-056, Baltimore, MD 21201-5159. E-mail address: ![]()
3 Abbreviations used in this paper: PBD, peptide-binding domain; CYT, cytoplasmic region; UT, untranslated. ![]()
Received for publication July 10, 1998. Accepted for publication July 19, 1999.
| References |
|---|
|
|
|---|
and LMP2. J. Immunol. 159:73.
This article has been cited by other articles:
![]() |
Y. Ohta, W. Goetz, M. Z. Hossain, M. Nonaka, and M. F. Flajnik Ancestral Organization of the MHC Revealed in the Amphibian Xenopus J. Immunol., March 15, 2006; 176(6): 3674 - 3685. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. H. Bos and B. Waldman Evolution by Recombination and Transspecies Polymorphism in the MHC Class I Gene of Xenopus laevis Mol. Biol. Evol., January 1, 2006; 23(1): 137 - 143. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Moon, S. M. Veniamin, J. A. Parks-Dely, and K. E. Magor The MHC of the Duck (Anas platyrhynchos) Contains Five Differentially Expressed Class I Genes J. Immunol., November 15, 2005; 175(10): 6702 - 6712. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Ohta, E. C. McKinney, M. F. Criscitiello, and M. F. Flajnik Proteasome, Transporter Associated with Antigen Processing, and Class I Genes in the Nurse Shark Ginglymostoma cirratum: Evidence for a Stable Class I Region and MHC Haplotype Lineages J. Immunol., January 15, 2002; 168(2): 771 - 781. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Aoyagi, J. M. Dijkstra, C. Xia, I. Denda, M. Ototake, K. Hashimoto, and T. Nakanishi Classical MHC Class I Genes Composed of Highly Divergent Sequence Lineages Share a Single Locus in Rainbow Trout (Oncorhynchus mykiss) J. Immunol., January 1, 2002; 168(1): 260 - 273. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. P. Shum, L. Guethlein, L. R. Flodin, M. A. Adkison, R. P. Hedrick, R. B. Nehring, R. J. M. Stet, C. Secombes, and P. Parham Modes of Salmonid MHC Class I and II Evolution Differ from the Primate Paradigm J. Immunol., March 1, 2001; 166(5): 3297 - 3308. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. R. Kirby An Ancient Transpecific Polymorphism Shows Extreme Divergence in a Multitrait Cline in an Intertidal Snail (Nucella lapillus (L.)) Mol. Biol. Evol., December 1, 2000; 17(12): 1816 - 1825. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |