The JI
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     
 


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Aoyagi, K.
Right arrow Articles by Nakanishi, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Aoyagi, K.
Right arrow Articles by Nakanishi, T.
The Journal of Immunology, 2002, 168: 260-273.
Copyright © 2002 by The American Association of Immunologists

Classical MHC Class I Genes Composed of Highly Divergent Sequence Lineages Share a Single Locus in Rainbow Trout (Oncorhynchus mykiss)1 ,2

Kazuhiko Aoyagi3,*, Johannes M. Dijkstra3,*, Chun Xia{dagger}, Ikuo Denda{ddagger}, Mitsuru Ototake4,*, Keiichiro Hashimoto§ and Teruyuki Nakanishi5,*

* National Research Institute of Aquaculture, Fisheries Research Agency, Tamaki, Mie, Japan; {dagger} Department of Microbiology and Infection, College of Veterinary Medicine, China Agricultural University, Beijing, People’s Republic of China; {ddagger} Nagano Prefectural Fisheries Experimental Station, Akashina, Nagano, Japan; and § Institute for Comprehensive Medical Science, Fujita Health University, Toyoake, Aichi, Japan


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 
The classical MHC class I genes have been known to be highly polymorphic in various vertebrates. To date, putative allelic sequences of the classical MHC class I genes in teleost fish have been reported in several studies. However, the establishment of their allelic status has been hampered in most cases by the lack of appropriate genomic information. In the present study, using heterozygous and homozygous fish, we obtained classical-type MHC class I sequences of rainbow trout (Oncorhynchus mykiss) and investigated their allelic relationship by gene amplification and Southern and Northern hybridization analyses. The results indicated that all MHC class I sequences we obtained were derived from a single locus. Based on this, a unique polymorphic nature of the MHC class I locus of rainbow trout has been revealed. The mosaic combination of highly divergent ancient sequences in the peptide-binding domains is notable, and the variable nature around the boundary between the {alpha}3 and transmembrane domains is unprecedented.


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 
Polymorphic MHC class I and II molecules play a crucial role in the acquired immune system by presenting antigenic peptides to T cells (1, 2). The classical MHC class I molecules are expressed in most cell types and present endogenous peptides at the cell surface for recognition by the TCR/CD8 complex.

Although T cells have not been clearly identified in fish, several studies showed the presence of lymphocytes with T cell-like functions, as exhibited in the mixed leukocyte reactions, allograft rejection, and cell-mediated cytotoxicity (3). Similarity between the fish and mammalian immune system has been also revealed by sequence analyses on principal molecules in the immune system. To date, genes for MHC molecules such as class IA, {beta}2-microglobulin ({beta}2m),6 class IIA, and class IIB genes have been reported from more than 25 species of teleosts and elasmobranchii (3, 4, 5). The identification of these polymorphic MHC genes and also of TCR and Ig genes (6) indicates that the acquired immune system is already complex at the level of fish, although current data indicate that jawed fish form the most primitive group with an MHC/TCR system (7).

High sequence variability in MHC class I genes has been reported for such fish species as zebrafish (8), cichlids (9), guppy (10), Atlantic salmon (11), the coelacanth Latimeria chalumnae (12), and shark (13), and apparent sequence lineages could be noticed in some cases. The classical nature of class I molecules includes high allelic polymorphism, the conservation of important amino acid residues revealed in the mammalian classical MHC class I molecules, and expression in most tissues. Due to the lack of appropriate assay systems in fish, classical MHC class I function has not been clearly shown to date. In addition to plausible classical MHC class I genes, apparent nonclassical genes and related pseudogenes were also described for fish (14, 15). Besides the unresolved question about the classical or nonclassical functions of the fish genes, a major impediment in understanding the fish MHC class I system is that most studies do not reveal whether the observed sequence variability is attributed to allelic polymorphism or to variability among sequences derived from different loci. A good exception is the study of Okamura et al. (13), which showed high sequence variability connected to a single locus in shark. Many other studies on the variability of fish classical MHC class I have been limited to the {alpha}1 and {alpha}2 domains without extensive analysis on allelic relationship.

In rainbow trout, classical MHC class I sequences have been described with several full open reading frame sequences (16, 17, 18) and a number of sequences encompassing one or both of the {alpha}1 and {alpha}2 regions (14, 17). Furthermore, nonclassical MHC class I sequences of rainbow trout have also been reported (18, 19). However, the number of the classical MHC class I loci has not been clearly addressed in the previous studies. Although the existence of at least two loci was expected based on the observation that three classical-type class I genes could be amplified by RT-PCR and also could be isolated from a cDNA library from a single fish, more analyses should be performed to reach the conclusion (17). The sequence information focusing on the {alpha}1 and {alpha}2 regions did not address the number of loci clearly, although the existence of two classical loci, A and UA, has been proposed based on the sequence variability (14). Our present study defined a single expressed classical MHC class I locus in rainbow trout that exhibits a unique polymorphic nature.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 
Animals used

Nine outbred rainbow trout (Oncorhynchus mykiss) were used in this analysis: fish W1 and W2 of strain Westralian (20), fish N1 and N2 of strain Nagano (from the Nagano Prefectural Experimental Station of Fisheries, Nagano, Japan), and fish T1, T2, T3, T4, and T5 of strain Tamaru (from the National Research Institute of Aquaculture, Tamaki, Japan). Five different homozygous clones were obtained from strain Nagano: AA1, AB2, AB5, AB7, and C25. The homozygous isogeneic trout were produced by gynogenesis in two generations by suppression of mitosis and meiosis in the first and second generations, respectively (21). Clonality was confirmed by DNA fingerprinting (data not shown). All fish analyzed were adults.

PCR primers and probes

The positions of primers and amplified fragments are indicated in Fig. 1Go. For initial amplification of the {alpha}2 domain, degenerate primers were derived from conserved regions for intramolecular disulfide bonds in MHC class I genes in chicken (GenBank accession no. X12780), frog (L20733), and human (K02883): pMix-{alpha}2f, 5'-ACTGTCGACCARHNGATGTAYGGNTGT (D = A + T + G, H = A + C + T, n = A + C + G + T, R = A + G, Y = C + T) (forward) and pMix-{alpha}2r, 5'-ACTGTCGACYTTNARCCAYTCDATRCA (reverse).



View larger version (70K):
[in this window]
[in a new window]
 
FIGURE 1. (continues)

 
The 5'-RACE PCR was performed using primer p4-{alpha}2r1 derived from the Onmy-UBA*401 {alpha}2 region, 5'-CTCATCATCCCACTCACAACCGTACAT (reverse), and anchor primer pVsq, 5'-TGGCGGCCGCTCTAGAACTAGTGGATC derived from vector {lambda}ZAP-II (Stratagene, La Jolla, CA).

For full-length cDNA amplification, 3'-RACE PCR was performed with forward primers derived from the 5' untranslated region (UTR) of Onmy-UBA*401 or Onmy-UCA*C32 and anchor primer NotI dT18 (first-strand cDNA synthesis kit; Pharmacia Biotech, Uppsala, Sweden). The forward primers were pIV-5'f, 5'-ATAACACTTCCACATTGACGGGCC from Onmy-UBA*401, and pI-5'f, 5'-CCTATTGGAGATAACATACTTCGTC from Onmy-UCA*C32. PCR products were cloned, and for every amplified fragment 18 clones were sequenced to screen for possible gene or allele variants. The sequences could be divided into lineages and sublineages. Minor differences within a sublineage were not due to sequencing mistakes, as confirmed by an independent PCR analysis performed for all fish.

For detection of sublineage-specific expression by PCR, primer sets derived from the 5' and 3' UTRs were used for amplification from single-stranded cDNA. For Onmy-UBA*101, primer pI-5'f was used in conjunction with primer p1–3'r, 5'-GTACTATTCCACTCAAGCAATTAC (reverse). For Onmy-UBA*201 and *301, pI-5'f was used in conjunction with p234-3'r, 5'-CATTCCACTCATTCCATTCTAGCA (reverse). Primers for specific amplification of Onmy-UBA*401 were pIV-5'f and p234-3'r. For detection of Onmy-UBA*501, primer pI-5'f and p5–3'r, 5'-AACATACCTGTAAGAATGAGAAGT (reverse) were used. Primers for specific detection of Onmy-UBA*601 or *701 expression were pI-5'f (forward) and p67–3'r1, 5'-AAGATCAAATGGCATTAAACAAATGG (reverse). Most sublineages could be distinguished, but the primers could not differentiate between alleles of Onmy-UBA*6 and *7 or between alleles of Onmy-UBA*2 and *3.

Various digoxigenin-labeled probes for Northern and Southern blot analysis were prepared by PCR Dig probe synthesis kit (Boehringer, Tokyo, Japan) using primers and templates as follows: the 170-bp {alpha}2 probe of Onmy-UBA*101, pr1-{alpha}2, primers, p1-{alpha}2f, 5'-GGATGATGACACTG GAGCC (forward) and p1-{alpha}2r, 5'-AATGTTTTCTCCGCTCCTCAT (reverse), and template, Onmy-UBA*101 cDNA clone; the 181-bp {alpha}2 probe of Onmy-UBA*201, pr2-{alpha}2, primers, p2-{alpha}2f, 5'-GGATGATGAGGCTG GAGCC (forward) and p2-{alpha}2r, 5'-CTGGGTGCAGTAGTTTTTTC (reverse), and template, Onmy-UBA*201 cDNA clone; the 181-bp {alpha}2 probe of Onmy-UBA*301, pr3-{alpha}2, primers, p3-{alpha}2f, 5'-GGATGATGAGACTG GAGTC (forward) and p3-{alpha}2r, 5'- CTGGGTGATGTAGTTTTTC (reverse), and template, Onmy-UBA*301 cDNA clone; the 220-bp {alpha}2 probe of Onmy-UBA*401, pr4-{alpha}2, primers, p4-{alpha}2f, 5'-GAGGATGTACGGTT GTGAGT (forward) and p4-{alpha}2r2, 5'- CAGCCACTCAATGCACGTC (reverse), and template, Onmy-UBA*401 cDNA clone; the 170-bp {alpha}2 probe of Onmy-UBA*501, pr5-{alpha}2, primers, p5-{alpha}2f, 5'-GAATGATGAGACTG GTGCC (forward) and p5-{alpha}2r, 5'-AATGTTTATCCCGCTCAGTCAT (reverse), and template, Onmy-UBA*501cDNA clone; the 189-bp {alpha}2 probe of Onmy-UBA*601, pr6-{alpha}2, primers, p67-{alpha}2f,5'-GCTGATGTATGGCT GTG (forward) and p67-{alpha}2r, 5'-TTCTTCTGAAAGTTAGCCT (reverse), and template, Onmy-UBA*601 cDNA clone; the 189-bp {alpha}2 probe of Onmy-UBA*701, pr7-{alpha}2, primers, p67-{alpha}2f and p67-{alpha}2r, and template, Onmy-UBA*701 cDNA clone; the 241-bp {alpha}3 probe of Onmy-UBA*401, pr4-{alpha}3, primers, p4-{alpha}3f, 5'-GTCCCTCCGTCAGTGTCTCTG (forward) and p4-{alpha}3r, 5'-TGACTTGAACCACACACTGATACT (reverse), and template, Onmy-UBA*401 cDNA clone; the 377-bp 3' UTR probe of Onmy-UBA*601, pr6–3' UTR, primers, p67–3'f, 5'-GTATACTTACACA CACATCAG (forward) and p67–3'r2, 5'-GATAATTCAACACAATTA CAGC (reverse), and template, Onmy-UBA*601 cDNA clone.

Construction of cDNA and isolation of genomic DNA

Total RNA was isolated from the kidney (including head kidney) using a RNA extraction kit (Pharmacia Biotech), followed by mRNA purification with a mRNA purification kit (Pharmacia Biotech). For the RT-PCR used in initial amplification of the {alpha}2 region, the 3'-RACE PCR, and the PCR for detection of lineage-specific expression, first single-strand cDNA was synthesized from purified mRNA by using the NotI dT18 primer (Pharmacia Biotech) and a single-strand cDNA synthesis kit (Pharmacia Biotech). For 5'-RACE PCR, first a cDNA library was constructed from isolated mRNA by use of pd(N)6 primer (Pharmacia Biotech), the Time Saver kit (Pharmacia Biotech), and the {lambda}ZAP II cloning system (Stratagene). In all cases, manufacturers’ recommendations were followed.

Genomic DNA of rainbow trout was purified from the pectoral fin tissue by phenol extraction (22).

Conditions of PCR from cDNA or genomic DNA

The general conditions for all PCR amplifications using cDNA were 100-µl volume, 2.5 U of KOD DNA polymerase (TOYOBO, Tokyo, Japan), KOD buffer, 0.2 mM dNTPs, and 125 nM of each primer.

Amplification of the {alpha}2 domain of MHC class I from an AA1 fish with degenerate primers was performed on single-stranded cDNA from 50 ng of mRNA. PCR conditions were 94°C for 3 min, then 30 cycles of 94°C for 30 s, 50°C for 30 s, 72°C for 1 min, and finally 72°C for 10 min.

The 5'-RACE PCR was performed on double-stranded cDNA with {lambda}ZAP II arms obtained from 50 ng of mRNA. PCR conditions were 94°C for 3 min, then 30 cycles of 94°C for 1 min, 50°C for 1 min, 72°C for 2 min, and finally 72°C for 10 min.

The 3'-RACE PCR was performed on single-stranded cDNA from 50 ng of mRNA. PCR conditions were 94°C for 3 min, then 30 cycles of 94°C for 30 s, 50°C for 30 s, 72°C for 2 min, and finally 72°C for 10 min.

For expression analysis with specific primers, PCR was performed on single-stranded cDNA from 50 ng of mRNA. PCR conditions were 94°C for 3 min, then 30 cycles of 94°C for 30 s, 60°C for 30 s, 72°C for 1 min, and finally 72°C for 10 min.

PCR with primers pMix-{alpha}2f and pMix-{alpha}2r for amplification of Onmy-UFA was performed using 50 ng of genomic DNA. The 50-µl reaction mixtures contained 1.25 U of AmpliTaq Gold polymerase (PerkinElmer, Tokyo, Japan), AmpliTaq Gold buffer, 0.2 mM dNTPs, and 1 µM of each primer. PCR conditions were 94°C for 10 min, then 43 cycles of 94°C for 1 min, 55°C for 1 min, 72°C for 2 min, and finally 72°C for 10 min.

Sequence analysis

PCR products were cloned into PCR-script Amp SK+ vector (Stratagene). The nucleotide sequences were determined by the dideoxy chain termination method using a dRhodamine terminator cycle sequencing kit (Applied Biosystems, Foster City, CA) and suitable primers; subsequent analysis was performed with an automated sequencer (ABI PRISM 377 DNA sequencer; Applied Biosystems). Comparison of nucleotide sequences and deduced amino acid sequences was performed using the programs Search Homology, Multiple Alignment, Create Text Document, and unweighted pair-group method with arithmetic mean (UPGMA) of GENETYX version 9.0 (Software Development, Tokyo, Japan) computer software.

Southern and Northern blot analysis

For both Southern and Northern blot analysis, digoxigenin-labeled probes were made, and hybridization, washing, and detection were performed by using the Dig Lumin Detection kit (Boehringer). The manufacturer’s recommendations were followed.

For Southern blot analysis, genomic DNA was digested by the restriction enzyme PstI, subjected to electrophoresis in 0.8% agarose gels, and blotted onto a nylon membrane. After hybridization with the probes, the stringency of washing was twice for 15 min with 1x SSC/0.1% SDS at 68°C. For Northern blot analysis, total RNA from specific tissues was isolated by TRIzol (Life Technologies, Grand Island, NY) and electrophoresed through a 1% agarose/10% formalin gel, and then blotted onto a nylon membrane. After hybridization with the probes, the stringency of washing was twice for 15 min with 2x SSC/0.1% SDS at 68°C.


    Results and Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 
Isolation of MHC class I genes from rainbow trout

When we started this study several years ago, the isolation of MHC class I genes from the rainbow trout was not reported. Thus, we tried to isolate MHC class I genes from this species and initially succeeded to amplify a 240-bp fragment from a kidney cDNA library of the clonal homozygous strain AA1 by use of degenerate primers based on the conserved {alpha}2 domain motifs in the classical MHC class I genes of other vertebrates. The amplified fragment, denoted Onmy-UBA*401-{alpha}2, showed homology with the {alpha}2 domain of the MHC class I genes of other species, and this was regarded as a part of a candidate for a classical MHC class I gene as it encodes, in this domain, the conserved amino acid residues critical for the interaction with a bound peptide. Then we conducted 5'-RACE PCR to obtain the 5' region, Onmy-UBA*401-5' end and found that the {alpha}1 domain sequence of this gene has some characteristic features that are somewhat deviated from those of the typical classical MHC class I sequences (Figs. 1Go and 2Go, discussed later). More investigation was necessary to obtain a clear picture of the classical MHC class I genes of the rainbow trout.



View larger version (63K):
[in this window]
[in a new window]
 
FIGURE 2. (continues)

 
In the present study, nine outbred rainbow trout belonging to three different strains plus five different homozygous clonal fish derived from one of these strains by artificial gynogenesis were analyzed for the MHC class I sequences. Because the comparison of Onmy-UBA*401-5' end (5' UTR, leader, {alpha}1, and {alpha}2 domain sequences) with the subsequently published rainbow trout MHC class I sequence Onmy-UCA*C32 (16) showed sequence differences in all domains (Fig. 1Go), we prepared two kinds of 5' primers, pIV-5'f and pI-5'f, based on our gene and C32, respectively, to conduct 3'-RACE PCR. Using these two primers, we obtained 2- to 2.3-kb fragments from all fish investigated. Unexpectedly, in all cases, only PCR with either pI-5'f or pIV-5'f was successful. For homozygous fish, all 18 clones isolated from each fish appeared to be identical, except for single nucleotide changes in individual clones probably caused by PCR error. In most outbred fish, two truly different sequences were detected. A total of 22 independent PCR products was analyzed. The sequences contain open reading frames of approximately 1.1 kb with homology with known classical MHC class I genes in all regions: the leader peptide; the {alpha}1, {alpha}2, and {alpha}3 domains; and the transmembrane (TM) and the cytoplasmic (CY) regions. The 22 PCR products revealed 10 different sequences. Based on >99% nucleotide identity, they were divided into seven sublineages that were designated Onmy-UBA*1, *2, *3, *4, *5, *6, and *7. The individual sequences of a sublineage were numbered Onmy-UBA*101, *102, and so on. Table IGo shows the sequences detected in the fish investigated. Fig. 1Go compares the nucleotide sequences of single representatives of the sublineages. The homology between sublineages varies between 73 and 98% nucleotide identity. Within one sublineage, the full-length cDNAs show only one to two nucleotide changes in the {alpha}2 or 3' UTR (Fig. 1Go).


View this table:
[in this window]
[in a new window]
 
Table I. The MHC class I cDNAs detected

 
Variability between sublineages seems to occur in all regions and appears to be caused by several mechanisms, including replacement, deletion, and insertion. The most noticeable differences between sublineages in the {alpha}3 region, TM/CY region, and 3' UTR are frequent insertions or deletions. Small direct repeats such as the AC repeats in the 3' UTR (Fig. 1Go) seem to be hot spots for mutations.

To confirm that in all fish analyzed only one or two of the different sublineages were expressed, five sublineage-specific primer sets were derived from the 5' and the 3' UTR, as described in Materials and Methods, and fish mentioned in Table IGo were tested for fragment amplification. Both positive and negative PCR data obtained (not shown) fully matched the data derived from sequencing full-length cDNA clones, as described above. This means that probably none of the seven sublineages has escaped detection when expressed.

The pink salmon sequence Ongo-UA-(92H) and the Atlantic salmon sequence Sasa-p30 show high homology with the Onmy-UBA sequences in all their domains, including the 3' UTRs (Fig. 1Go).

Analysis of deduced amino acid sequences

Fig. 2Go shows a comparison of the deduced amino acid sequences of representatives of the seven sublineages described in this study with classical MHC class I molecules from different phylogenetic groups. Atlantic and pink salmon are other salmonid species, carp and zebrafish are cyprinids, and sharks belong to the elasmobranchii; the function of their indicated MHC class I molecules has not been shown, but their classical nature is indicated by sequence analysis. It is noticeable, though, that the depicted zebrafish sequence Brre-UBA*01 lacks the highly conserved Y7 residue (see below) and therefore possibly represents a nonclassical MHC class I molecule.

In all the molecules of Fig. 2Go, some cysteine residues are conserved, and there is a N-glycosylation site at position 86 (residue numbers mentioned in this section are in agreement with Fig. 2Go). A similar structure of the depicted proteins is also expected by conservation of some salt bridges described for human MHC class I (23) in the Onmy-UBA molecules: H3-D28, R41-E61 (except in Onmy-UBA*401), H92-D118, K143-D147, K157-E161 (only in Onmy-UBA*601 and *701), and E166-K169 (D166-K169 in Onmy-UBA*201 and *301). Conserved for mammalian classical MHC class I molecules are the amino acids Y7, Y59, Y84, T142, K145, W146, Y159, and Y171, which are involved in binding the termini of the MHC class I-bound peptide (24). The nonmammalian molecules depicted in Fig. 2Go show almost complete conservation of these amino acids, except for Y84, which is changed into R84 in fish, amphibian, and bird MHC class I molecules. Other exceptions are the lack of Y59 in Onmy-UBA*401 and the presence of L146 instead of W146 in Onmy-UBA*701. Variability in the residues at position 59 or 146 has also been described for some alleles of human HLA-B and HLA-C (25). Concerning the binding of human CD8 to MHC class I, several investigators described the importance of an acidic stretch in the {alpha}3 domain (26, 27). Although not identical to those in humans, the Onmy-UBA sequences all contain an acidic stretch at this position (Fig. 2Go). Of the three residues (Q114, D121, and E127) in the {alpha}2 domain shown to affect CD8 binding in human MHC class I (28), Q114 and D121 have been completely conserved in the trout sublineages. Also, residues that may interact with {beta}2m as F8, T10, V24, Q95, Q114, D118, G119, D121, E229, D235, and Q239 (23) are conserved in mouse (not shown), human, and rainbow trout (Fig. 2Go).

An additional indication that the Onmy-UBA sequences represent classical MHC class I molecules is the distribution of variability. Whereas the {alpha}3 regions of the different sublineages share >73% amino acid identity, within the {alpha}1 and {alpha}2 regions this can be as low as 49 and 53%, respectively. The variability is especially high at those positions involved in binding antigenic peptides and determining their specificity in humans (indicated as P in Fig. 2Go) (29). The average number of detected amino acids in the {alpha}1 and {alpha}2 domains of the Onmy-UBA sequences is 1.6 at non-P positions and 2.6 at P positions.

Most interestingly, Onmy-UBA*601 and *701 have 14-aa extensions around the carboxyl terminus of the {alpha}3 domain, which has already been described for the trout sequence Onmy-UCA*C32 (16). For Onmy-UBA*501, a similar small extension is found as described for the Atlantic salmon sequence Sasa-p30. Onmy-UBA*101, *201, *301, and *401 have no extensions as has pink salmon Ongo-UA-(92H). Thus, not only does the length of this region vary between salmonid species, but also probably within one locus in one species. Because the residues of the extensions are largely hydrophilic, they are probably located extracellularly and may influence the distance of the {alpha}3 domain to the membrane.

Whereas the {alpha}2 region of the sublineages Onmy-UBA*1, *2, *3, *4, and *5 are more similar to those detected in other salmonids, the {alpha}2 regions of Onmy-UBA*6 and *7 are more closely related to those detected in the cyprinid sequences depicted in Fig. 2Go. They will be discussed below as type I and type II {alpha}2 regions, respectively.

Several full-length rainbow trout classical-type MHC class I sequences were reported previously, and they are similar to our Onmy-UBA sequences. Those are Onmy-UCA*C32 (16), Onmy-UAA*KD6, *HC-01, *OSU-01, Onmy-UBA*SP3 (17), and Onmy-UBA*0101 (18), and the similarity can be extended to the 3' UTR sequences. The variability in the {alpha}1 and {alpha}2 sequences of rainbow trout class I sequences reported by Miller et al. (14) and by Hansen et al. (17) is also within the range of the variability we detected in our Onmy-UBA sequences. The classical nature of the class I sequences of rainbow trout was examined in the previous studies (16, 17, 18) and was suggested based on their expression in various tissues and the presence of predicted key amino acid residues known to be conserved in the classical MHC class I molecules of other vertebrates. The classical nature of our Onmy-UBA sequences is also apparent because of their expression profiles (see below), high variability in the {alpha}1 and {alpha}2 domains, and the conservation of some critical amino acids.

Ongoing studies of the Onmy-UBA locus, including those of MHC class I-restricted cytotoxicity, up-regulation of mRNA expression by viral infection, and the class I protein expression at the cellular level using a mAb, further support the conclusion that the UBA locus is a classical MHC class I locus.

Phylogenetic comparison

Fig. 3Go shows a UPGMA analysis of the extracellular domains of the seven sublineages compared with known probable classical MHC class I sequences of trout and a few other fish species. Comparison with the other fish is done to show the existence of trans-specific lineages. In the figure, the sequences are not shown individually, but as homology groups; no rearrangement of the computer software UPGMA analysis was made, the picture was only simplified. Members within each group are identified in Table IIGo and Fig. 3Go. Within the groups I-VII of {alpha}1, and I-II of {alpha}2, the members share >67% amino acid identity with each other. Within the groups I-V of the {alpha}3 region, the members share >88% amino acid identity. Most trout {alpha}1 sequences belong to group {alpha}1-I. Group {alpha}1-IV, containing Onmy-UBA*401, is very different from group {alpha}1-I (±50% amino acid identity). Group {alpha}1-V contains both a cyprinid, Brre-UAA*01, and a salmonid, Onmy-UA*A4.11, sequence. For the {alpha}2 region, only two groups were detected. As mentioned before (16), some trout sequences show high similarity to cyprinid {alpha}2 domains. The division into groups I-III of the {alpha}3 domain shows the existence of trans-specific lineages between salmonids. Group {alpha}3-II contains a small extension at the end of the {alpha}3 domain, and group {alpha}3-III a large extension (see also Fig. 2Go). Analysis of the TM/CY region is not shown, as the present data do not allow a natural division of the sequences of this region into groups. Table IIGo schematically depicts the organization of the sequences of Fig. 3Go that encompass more than one region. The salmonid sequences containing a full open reading frame are divided into lineages Sal-MHCIa*A-D by classifying their extracellular regions into different groups. The sequences Onmy-UBA*101, *201, and *301 can be classified into the same lineage as pink salmon Ongo-UA-(92H), and Onmy-UBA*501 and Atlantic salmon Sasa-p30 group together. This definition of lineage is chosen, as recombination events creating different combinations of {alpha}1 and {alpha}2 groups seem to have been an important cause of polymorphism. For example, the sequences Onmy-UBA*101 and *701 have a completely identical {alpha}1 region, but a very different {alpha}2 region, with differences beginning at the domain border (Fig. 1Go). And Onmy-UBA*601 and Onmy-UA*A5.1 (17) have an identical {alpha}2 domain, but a very different {alpha}1 domain.



View larger version (23K):
[in this window]
[in a new window]
 
FIGURE 3. The amino acid sequences of the {alpha}1, {alpha}2, and {alpha}3 regions of published rainbow trout (rt) MHC class Ia molecules and those of several other fish (Atlantic salmon (as), pink salmon (ps), common carp (ca), zebrafish (zb)) are compared by UPGMA. Numbers indicate amino acid substitutions per site. Division into regions is in agreement with Fig. 2Go. The clustering of the individual sequences is shown in a simplified way, by only depicting the several homology groups, indicated with Roman numbers. The number of sequences included in each group, as well as the amino acid substitutions per site between group members, is indicated in parentheses. Compared are the sequences as described in Fig. 2Go; furthermore, the rainbow trout sequences starting with Onmy-A*, Onmy-UA*, Onmy-UAA*, Onmy-UBA*, or Onmy-UCA* designation reported by others (14 16 17 ) (GenBank accession nos. U55380, AF002171–2179, AF091785, AF104523–4529, AF104579-4582, and AF115518–5528), and the zebrafish sequence Brre-UAA*01 (8 ) (GenBank accession no. Z46776).

 

View this table:
[in this window]
[in a new window]
 
Table II. Lineage organization of salmonid MHC class I sequences1

 
Southern blot analysis

Even for a homozygous clonal fish, more than nine bands can be detected in Southern blot analysis when using a {alpha}3 probe (data not shown), indicating a number of similar genes or pseudogenes. Therefore, the {alpha}3 probe was not suited to identify the Onmy-UBA locus, and locus identification was performed with {alpha}2 probes.

Genomic DNA of all the nine outbred and the five homozygous fish used in this study was digested with PstI and used for Southern blot analysis. For all sublineages, probes derived from the {alpha}2 region were made. The hybridization pattern of all type I {alpha}2 probes (from Onmy-UBA*101, *201, *301, *401, or *501) was identical; Fig. 4GoI was obtained with probe pr5-{alpha}2 derived from Onmy-UBA*501. The hybridization pattern of all type II {alpha}2 probes (from Onmy-UBA*601 or *701) was also identical; Fig. 4GoII was obtained with probe pr6-{alpha}2 from Onmy-UBA*601. The patterns indicate that under the washing conditions used, all type I {alpha}2 probes can recognize all genomic fragments with a type I {alpha}2 region, and all type II {alpha}2 probes can recognize all type II {alpha}2 regions in the genome; there seemed to be no cross-reaction between the two types of {alpha}2 regions.



View larger version (38K):
[in this window]
[in a new window]
 
FIGURE 4. Southern blot analysis with {alpha}2 probes. Genomic DNA of all fish mentioned in Table IGo was digested by restriction enzyme PstI and electrophoretically separated in 0.8% agarose gels, 10 µg/lane. After transfer, nylon membranes were either hybridized with a type I {alpha}2 probe, pr5-a2 (I), or a type II {alpha}2 probe, pr6-{alpha}2 (II). Identical fragments were found for fish expressing identical sublineages. Therefore, the type I {alpha}2 regions are thought to be located on a 2.4-kb PstI fragment for Onmy-UBA*1 and *2 (1 2 ), on a 2.3-kb fragment for Onmy-UBA*3 (3 ), on a 4-kb fragment for Onmy-UBA*4 (4 ), and on a 1.9-kb fragment for Onmy-UBA*5 (5 ). The type II {alpha}2 regions of Onmy-UBA*6 and *7 (6 7 ) are located on a 2.8-kb PstI fragment. The numbers of the sublineages are indicated at the right of the figure. The probe pr6-{alpha}2 also hybridized to a 2.4-kb Onmy-UFA fragment (F).

 
Type I {alpha}2 probes only hybridized to one restriction fragment for homozygous fish; in outbred fish, either none, one, or two fragments were detected (Fig. 4GoI). Hybridization fragments appeared sublineage specific. For fish expressing Onmy-UBA*1 or *2, a 2.4-kb; Onmy-UBA*3, a 2.3-kb; Onmy-UBA*4, a 4-kb; and Onmy-UBA*5, a 1.9-kb fragment was detected, respectively. In cDNA sequences of fish T4, no type I {alpha}2 region could be detected, and also Southern blot analysis could not detect such a region. The type II {alpha}2 probes detected one fragment of 2.4 kb in every fish (Fig. 4GoII). For the fish T1, T4, and T5, which were shown by cDNA analysis to express sublineages Onmy-UBA*6 and/or *7, a second fragment of 2.8 kb was found. Because the combined type I and II {alpha}2 probes detected only one polymorphic fragment per haploid genome set, the most probable explanation is that all the seven sublineages detected by cDNA analysis belong to the same locus, designated Onmy-UBA.

To clarify the nature of the 2.4-kb band detectable in all fish as shown in Fig. 4GoII, we tried to amplify novel {alpha}2 sequences from genomic DNA by the use of degenerate primers. We succeeded in isolation of a new {alpha}2 DNA fragment and named its gene Onmy-UFA. As the new {alpha}2 sequence can be classified as a type II {alpha}2, it should correspond to the 2.4-kb fragment in Fig. 4GoII. It turned out to be closely related to reported sequences such as Onts-HN-1 and Sasa-UA*8 (Fig. 5Go). However, this Onmy-UFA fragment has defects in coding an intact MHC class I protein; a deletion of a DNA segment resulting in a frameshift and the presence of a stop codon could be noticed in the sequence (Fig. 5Go). We amplified this Onmy-UFA fragment from various fish, and the sequence analyses revealed its nonpolymorphic nature. Thus, every fish possesses this Onmy-UFA, as revealed in the Southern hybridization (Fig. 4IGoI), and its expression cannot be detected (see below, Fig. 6Go), in agreement with the mutational defects in the sequence.



View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 5. Comparison of the Onmy-UFA fragment with other {alpha}2 type II sequences. Dots indicate identical nucleotides. Dashes indicate gaps in the sequences. An identical fragment could be amplified from each fish investigated (W2, AA1, and C25). Compared with the other {alpha}2 type II sequences, Onmy-UFA*101 has a 38-bp deletion, resulting in a frame shift, and a stop codon that interrupts the homologous open reading frame (TAG at position 86). An even higher homology than with the reported rainbow trout {alpha}2 type II sequences was found with the {alpha}2 domains of the chinook salmon (Oncorhynchus tshawytscha) sequence Onts-HN-1 and the Atlantic salmon (Salmo salar) sequence Sasa-UA*8. These sequences were amplified by Miller et al. (14 40 ), and neither expression nor the presence of a full-length open reading frame has been shown. However, those sequences are probably derived from a polymorphic locus, because closely related sequences have been detected in the same species (14 40 ).

 


View larger version (79K):
[in this window]
[in a new window]
 
FIGURE 6. Northern blot analysis of total RNA from several tissues of clonal trout AA1. A total of 5 µg of RNA/lane was electrophoretically separated in a 1% agarose/10% formalin gel. After transfer, the nylon membrane was hybridized with the {alpha}2 probe pr4-{alpha}2. In all tissues, a 2.4-kb transcript was detected. As size markers, ribosomal RNAs are indicated.

 
Tissue-specific expression

Ubiquitous tissue expression is an indication for the classical nature of MHC class I genes (30). Fig. 6Go shows the hybridization pattern of probe pr4-{alpha}2 ({alpha}2 domain of Onmy-UBA*401) with total RNA isolated from kidney (including head kidney), brain, gill, heart, liver, spleen, intestine, gonad, and muscle of an AA1 fish. For all tissues, expression of a 2.4-kb transcript was detected, but whereas for tissue from gills or intestine a strong signal was found, expression in liver or muscle was very low. The signals around 1.8 kb observed in intestine and gill are likely to be produced by the presence of ribosomal RNA. To date, we could not amplify shorter transcripts by 3' or 5'-RACE with intestine (data not shown). With a probe derived from the {alpha}3 region of Onmy-UBA*401, pr4-{alpha}3, the same expression pattern for different AA1 tissues as with pr4-{alpha}2 was detected (data not shown). With the pr6-{alpha}2 and pr7-{alpha}2 probes ({alpha}2 type II), no hybridization signal could be detected (data not shown). As the {alpha}2 type II sequences published to date for rainbow trout are very similar (>97% nucleotide identity), it suggests that these sequences are not or are hardly expressed in the different tissues of the AA1 fish.

The UBA*401 expression pattern is consistent with the findings of Hansen et al. for Onmy-UAA*KD6 (17) and Onmy-UCA*C32 (16), probable UBA alleles of lineages Sal-MHCIa*A and *D, respectively.

Sublineage-specific transcripts

Because sequence analysis indicated length differences between transcripts of the several sublineages, this was investigated further by Northern blot analysis for total RNA of the kidney (including head kidney) of 10 of the fish described. The {alpha}3 probe pr4-{alpha}3 indeed detected major transcripts of 2.1–2.9 kb (Fig. 7GoA). Except for the 2.9-kb transcripts for fish T1 and T4, the sizes correlated well with the obtained Onmy-UBA sequence data for each fish. Interestingly, in several fish determined by cDNA analysis to express only one Onmy-UBA sublineage, a minor transcript somewhat larger than the major transcript was detected (in fish N1, N2, AB2, AB7, and C25). Presence and length of the minor transcript also correlated with the Onmy-UBA sublineages detected by cDNA sequence analysis, and therefore might be derived from the same locus. The type I {alpha}2 probe pr4-{alpha}2 reacted more specifically than in Southern blot analysis, and recognized transcripts of the sublineages Onmy-UBA*2, *3, and *4 (Fig. 7GoB). The minor and major transcripts of fish AB2 (Fig. 7GoA) are both recognized by probe pr4-{alpha}2 (Fig. 7GoB), indicating their similarity. The type II {alpha}2 probe pr6-{alpha}2 recognized two transcripts for sublineages Onmy-UBA*6 and probably Onmy-UBA*7 (Fig. 7GoC; comparison with fish N1 and N2 indicates there is no hybridization with Onmy-UBA*101 transcripts). The similarity of the two transcripts recognized by the pr6-{alpha}2 probe was confirmed by probe pr6–3' UTR derived from the 3' UTR of Onmy-UBA*601 (Fig. 7GoD). The recognized bands in Fig. 7Go, B–D, were identical to those recognized by the {alpha}3 probe (Fig. 7GoA).



View larger version (45K):
[in this window]
[in a new window]
 
FIGURE 7. Sublineage-specific transcripts. Northern blot analysis: 15 µg of total RNA isolated from the kidney (including head kidney) of fish N1, N2, W2, AB2, AA1, AB5, AB7, C25, T1, and T4 was electrophoretically separated in 1% agarose/10% formalin gels. After transfer, the nylon membranes were hybridized with the Onmy-UBA*401-derived {alpha}3 probe pr4-{alpha}3 (A), the Onmy-UBA*401-derived {alpha}2 probe pr4-{alpha}2 (B), the Onmy-UBA*601-derived {alpha}2 probe pr6-{alpha}2 (C), and the Onmy-UBA*601-derived 3' UTR probe pr6–3' UTR (D). The same membrane with intermediate washing was used subsequently for C, B, and A. D, Obtained from an independent blot. As size markers, ribosomal RNAs are indicated.

 
The Northern blot analysis in the present study seems to agree with the sequence data, except that in an allele-specific way the locus may give rise to two transcripts. A possible explanation for the two transcripts is the use of different polyadenylation signals, as suggested by comparison of Onmy-UCA*C32 and Onmy-UBA*701 (Fig. 1Go).

The observed variability in transcripts is not uncommon for MHC class I; in the frog Xenopus laevis transcripts of divergent lineages differ in size, and alternative use of polyadenylation signals has been found (31).

Allelic relationship of the MHC class I genes isolated from rainbow trout

Despite a great number of MHC class I sequences reported for teleost fish, interpretation of the data with respect to allelic polymorphism has been limited in most cases by lack of locus identification. To date, a few groups reported the MHC class I sequences from rainbow trout (14, 16, 17, 18). However, the MHC class I locus (or loci) of this species has not been clearly identified; therefore, allelic relationships of the obtained sequences have not been clarified. Our results defined a single classical MHC class I locus of rainbow trout and its unique allelic polymorphism.

In the present study, we concluded that all our MHC class I sequences are derived from the same locus, namely, Onmy-UBA, based on the following observations: 1) the PCR amplification of cDNAs with appropriate primers revealed a single sequence from a homozygous clonal fish, and essentially two sequences from an outbred fish (Table IGo); 2) the Southern blot data are basically consistent with the presence of a single expressed classical MHC class I sequence per haploid genome. Namely, the {alpha}2 domains of all the MHC class I cDNA sequences isolated by RT-PCR could be assigned to the sublineage-specific polymorphic bands in the Southern blot analysis (Fig. 4Go); and 3) the Northern blot data revealed that the isolated MHC class I cDNAs could represent the expressed MHC class I genes (that can hybridize with the {alpha}3 probe) whose transcripts exhibited sublineage-specific variations in length (Fig. 7Go). Thus, these results led us to a conclusion that the MHC class I cDNA sequences we obtained belong to a single expressed MHC class I locus, although they include sequences highly divergent to each other.

It should be mentioned that we could not draw our conclusion only through the amplification studies, because we could not strictly exclude a possibility of some fortuitous amplifications for single-locus interpretation. With the additional results of Southern and Northern blot analyses, it became clear that the expression of the MHC class I genes with the relevant {alpha}3 domain could be explained by our isolated genes. While intralocus homogenization outside the {alpha}1 and {alpha}2 regions has been reported in various vertebrate MHC class I genes (e.g., Refs. 13, 31 , and 32) and some of those regions can be used as locus markers, Onmy-UBA sequences showed significant variability at the {alpha}3, TM/CY, and 3' UTR (Fig. 1Go). The genes with distinct lineages that used to exist at different loci may have come together through genetic interactions. In Xenopus MHC class I locus in which two ancient lineages could be recognized, highly divergent CY and 3' UTR were observed, and those genes of different lineages were speculated to be derived from different loci (31).

The number of classical MHC class I loci in rainbow trout

To date, two classical MHC class I loci have been considered to exist in the rainbow trout genome by some other investigators. Mainly based on the sequence characteristics of the {alpha}2 domain, one group expected two classical type class I loci, A and UA in salmonid fish genome (14), which have been originally found by Grimholt et al. (11) in Atlantic salmon and by Hansen et al. (16) in rainbow trout, respectively. However, our present results revealed that MHC class I sequences with highly distinct lineages can be observed at a single locus of rainbow trout. By Hansen et al. (17), three MHC class I sequences could be amplified from a single fish in RT-PCR experiments, suggesting the existence of at least two loci. Furthermore, it was reported that three MHC class I sequences were also isolated from a cDNA library made from a single fish (Onmy-UCA*C32, -UAA*KD6, and -UBA*SP3) (17). As it has been speculated that the ancestor of salmonid fish has experienced the tetraploidization, and multivalents have been observed at meiosis (33), its possible effect on the number of MHC class I genes was mentioned by these authors, and the allelic relationship of the three sequences remains to be determined. It should be mentioned that the classical-type MHC class I sequences reported by other investigators, including those three just mentioned above, should be detected with the methods used in our study. Thus, it would be interesting to compare more details of our results with those of others to examine the discrepancy. It also should be mentioned that our study does not exclude a possibility of the existence of some other expressed classical MHC class I genes that possess {alpha}3 domain sequences distinct from (namely, not cross-hybridized with) those of Onmy-UBA genes.

Polymorphic nature and ancient lineages at the MHC class I locus of rainbow trout

With respect to the mode of the variability in the peptide-binding domains, basic mutational mechanisms such as point nucleotide substitutions may be observed by comparing the allelic sequences possessing the same sequence lineage (Fig. 1Go). We could observe polymorphism at the amino acid positions that could interact with a bound peptide (Figs. 1Go and 2Go). Thus, the basic mutation and selection mechanism found in the MHC class I genes of the other vertebrates also seem to be operating at this locus.

Importantly, however, mosaic combinations of the sequence lineages (discussed later) associated with the {alpha}1 or {alpha}2 domains were impressively observed in the peptide-binding domains (Table IIGo). We detected in our rainbow trout three types of combinations of sequence groups, namely I/I, I/II, and IV/I, for the combination of the {alpha}1/{alpha}2 domain lineages. In Table IIGo, we showed mosaic structures of the peptide-binding domains not only in our sequences, but also in the rainbow trout sequences reported by others, that may well belong to the Onmy-UBA locus described in the present study. The apparent domain shuffle has plausibly been created by intra- or interlocus recombination between sequences with divergent lineages that had been established by various mutational mechanisms over a long period of time.

In Table IIGo, some ancient lineages can be discerned in the peptide-binding domains. The Atlantic salmon and the Pacific salmon (including rainbow trout and pink salmon) were estimated to diverge from each other approximately 20 million years ago, and the ancestors of rainbow trout and of pink salmon approximately 10 million years ago (14). As discussed by Hansen et al. (16), the group II lineage of the {alpha}2 domain could be found both in rainbow trout and in fishes of cypriniformes (carp and zebrafish). Therefore, it was suggested that the formation of the traits of this lineage predated the separation of these two orders of fish groups (16), which may have occurred more than 130 million years ago (34). Subsequently, the existence of this lineage was also confirmed in other salmonid fishes, namely, coho, chinook, and Atlantic salmon (14). Furthermore, it was noticed that this lineage partly shares its characteristics with the coelacanth sequence Lach-UA-01 (14).

Hansen’s group (17) also reported two sequences with the group V lineage of the {alpha}1 domain (Onmy-UA*A4.11 and Onmy-UAA*KD1.5) from rainbow trout. We noticed that this lineage corresponds to that of the {alpha}1 domain of zebrafish Brre-UAA*01 (Fig. 3Go and Table IIGo), although not mentioned previously (17). We expected that the sequences with the group V lineage in the {alpha}1 domain would also belong to the Onmy-UBA locus reported in the present study, as the reported rainbow trout sequences with this lineage possess the group I {alpha}2 sequences that are exclusively associated with the Onmy-UBA locus in our study. Whereas these rainbow trout sequences possess a V/I combination for the {alpha}1/{alpha}2 domain lineages though, the sequence Brre-UAA*01 exhibits a V/II combination. Considering the mosaic pattern of the {alpha}1 and {alpha}2 lineages in UBA, we anticipated that a V/II combination could be found in rainbow trout as well, which was confirmed by our recent data (our unpublished observation).

Close relationship between the group I of the {alpha}1 domain in the salmonid fish and those in the neoteleost fish, such as guppy (the order cyprinidontiformes), cichlid (perciformes), rockfish (scorpaeniformes), and stickleback (Gasterosteiformes), has been noticed (reviewed in Ref. 14). The group I of the {alpha}2 domain also shares the sequence characteristics with the {alpha}2 domain sequences from the neoteleost fish (also reviewed in Ref. 14). Therefore, it is supposed that these lineage groups were present in the common ancestor of the salmonid and the neoteleost fish.

Although the group IV of the rainbow trout {alpha}1 domain has not been observed in cyprinid fishes or in fishes of other orders, this group may be ancient because of its divergent characteristics compared with the other lineage groups in this domain. Actually, the group IV sequence shares some characteristic amino acid residues, namely, G14, I15, E16, I32, D33, and Y34 with Cyca-UA1*01 sequence (Fig. 2Go), and it also shares these residues, except D33, with the sequence from a wild guppy (Pore-UA-W4.4) (10). Some of these residues could also be found in the group II, III, and V, but not in the group I sequences of the {alpha}1 domain.

Because the study of the variability of MHC class I genes in cyprinid fishes has still been limited, it is possible that groups apparently related to the group I or the group IV in the {alpha}1 domain and those related to the group I in the {alpha}2 domain are present in cyprinid fishes.

The observation that highly divergent lineages exist at a MHC class I locus has precedents in other vertebrates. Namely, in banded houndshark (Triakis scyllia), two divergent lineages of {alpha}2 sequences are present at an MHC class I locus (13). One of them can be observed at another locus, and the interaction between the two loci has been suggested. At the Xenopus MHC class I locus, the existence of two ancient lineages has been reported (31), and the divergence of those lineages, supposedly originally belonging to different loci, was estimated to occur more than 100 million years ago comparing with the divergence between the human and mouse MHC class I sequences. It is also conceivable that, in rainbow trout or in its ancestors, highly distinct sequence lineages had been produced at different loci, then those lineages came to exist at a single locus, namely Onmy-UBA, through the genetic interaction between those loci. Thus, at a single MHC class I locus of various vertebrates, the appearance of ancient distinct lineages that possibly had been originally formed at distinct loci has been observed. And in some cases, mosaic structures of the peptide-binding regions have been produced with distinct domain lineages or sequence motifs derived from them.

Comparison with a recently reported study

After the present study was submitted, a study by Shum et al. (41) reporting the diversity of the rainbow and brown trout MHC class I genes appeared. The most important claim of their study relevant to the present one is that a single trout MHC class I locus exhibits high polymorphism with some ancient divergent lineages and the apparent shuffling of the {alpha}1 and {alpha}2 domains, as observed in our study. However, the most critical point, that the MHC class I sequences from the respective fish species obtained in their study are derived from a single locus, appears not to be substantiated. The result that they did not obtain by PCR more than two sequences from an individual fish (actually only a single sequence in 12 of 25 rainbow or brown trout) is compatible with, but does not demonstrate, the allelic relationship of the sequences. In the present study, we stressed the importance of the allelic relationship of the isolated sequences on which the sequence comparison was based.

Concluding remarks

We established the allelic relationship of the isolated rainbow trout MHC class I sequences by PCR and Southern and Northern blot analyses, and also by using homozygous fish. Because highly divergent sequence lineages turned out to exist at this locus, such various extensive analyses should be necessary to draw a conclusion concerning allelic relationship of the isolated genes. The genes derived from this single locus, Onmy-UBA, essentially exhibit characteristics of classical MHC class I genes, although some members show divergent features. These rainbow trout MHC class I genes exhibit a broad range of variability. In addition to the polymorphism caused by nucleotide substitutions and insertions/deletions of nucleotides, the polymorphism based on mosaic combinations of the {alpha}1 and the {alpha}2 domain sequences with ancient lineages has been revealed. Thus, the present study clarified interesting combinatorial variations at a single rainbow trout MHC class I locus that should be an important characteristic of polymorphism in the peptide-binding domains and presumably broaden peptide-binding repertoire in this species. The example of rainbow trout presents a notable unique case of extensive domain shuffle at the MHC class I peptide-binding region in the vertebrates. The antiquity of some of the sequence lineages has been directly indicated by the presence of those sequence lineages in distantly related orders of fish groups. Finally, unprecedented variability at the {alpha}3/TM boundary region of rainbow trout MHC class I genes is noteworthy.



View larger version (73K):
[in this window]
[in a new window]
 
FIGURE 1B. (continues)

 


View larger version (64K):
[in this window]
[in a new window]
 
FIGURE 1C. Nucleotide sequences of the trout MHC class I cDNAs Onmy-UBA*101, *201, *301, *401, *501, *601, *701, and Onmy-UCA*C32 (16 ); the pink salmon Ongo-UA-(92H) (35 ); and the Atlantic salmon Sasa-p30 (11 ). Alignment is performed by computer software and modified by hand. Dots indicate identical nucleotides. Dashes indicate gaps in the sequences. The amino acid sequence of the Onmy-UBA*101-encoded protein is indicated above the second nucleotides of codons. Not indicated are the sequences Onmy-UBA*102 (differs from Onmy-UBA*101 by having an AC insertion at position 1333), Onmy-UBA*502 (differs from Onmy-UBA*501 by substitutions of C and A at positions 509 and 1172 by T and C, respectively), and Onmy-UBA*503 (differs from Onmy-UBA*501 by having the A at position 1172 substituted by C). The locations of various primers are shown, the name and direction of the primers are indicated with arrows, and the sequences of the primers are thinly underlined. Onmy-UBA*401-{alpha}2 contains the positions 400–625 of Onmy-UBA*401, and Onmy-UBA*401 5' end contains the positions 1–433 of Onmy-UBA*401. AC repeats in the 3' UTRs are underlined and described as such. The trout MHC class I sequence Onmy-UCA*C32 is identical to Onmy-UBA*701, except for extensions at the 3' and 5' UTRs. The polyadenylation signals probably used by Onmy-UBA*701 and Onmy-UCA*C32 are boxed. The open reading frames of Onmy-UBA*101 and Onmy-UAA-OSU*01 (17 ) are identical, but there is a single nucleotide difference in the partly sequenced (65-bp) 3' UTR of Onmy-UAA-OSU*01 (data not shown). Onmy-UBA*401 appears to be identical to Onmy-UBA*SP3 (17 ), except that Onmy-UBA*SP3 has a 274-bp insertion in its 3' UTR and the poly(A) tail is connected at a slightly different (4 nt apart) position (data not shown). The border between the {alpha}3 and TM/CY domains as depicted in Fig. 2Go is based on the intron-exon organization in this region, as determined for the sequences Onmy-UBA*401, *501, and *701 (data not shown). Near the end of the {alpha}3 domain sequence of Onmy-UBA*501, a sequence composed of eight bases (AGACCAAC) was noticed, and the region containing this sequence has been deleted in Onmy-UBA*101, and it appears to be duplicated in Onmy-UBA*701. Thus, the DNA region containing this sequence might have been involved in the formation of some structures to induce a deletion or a duplication of the DNA sequence, resulting in the shortened {alpha}3 domain and the elongated one, respectively.

 


View larger version (56K):
[in this window]
[in a new window]
 
FIGURE 2B. Predicted amino acid sequences of MHC class I molecules of rainbow trout and other vertebrates. Representatives of the Onmy-UBA sublineages, Sasa-p30 (11 ), Ongo-UA-(92H) (35 ), Cyca-UA1*01 (36 ), Brre-UBA*01 (8 ), Trsc-UAA*101 (13 ), Xela-UAA1f (37 ), Gaga-B-F12 (38 ), and HLA-A2 (39 ). Alignment is performed by computer software and modified by hand. Numbering of amino acids is performed for the Onmy-UBA*101 sequence, starting from the {alpha}1 domain. Dashes indicate identical amino acids. Asterisks indicate gaps in the sequence. Downward-pointing arrowheads indicate the conserved positions believed to interact with antigenic peptide termini in mammals (24 ), and upward-pointing arrowheads indicate conserved positions believed to form salt bridges (23 ) or to interact with CD8 (28 ) or {beta}2m (23 ). P, The positions believed to interact with the antigenic peptide and determine its specificity in mammals (29 ); C, conserved cysteine residues; the box at position 86, a conserved glycosylation site; the hooked line denoted CD8, the acidic domain believed to interact with CD8 in humans (26 ). The UBA*401 sequence has a deletion of two amino acid residues, including the conserved Y59 compared with the corresponding region of the mammalian molecules. Some flexibility and compensation may exist in MHC class I molecules in terms of both the structure around this region and the hydrogen bonding to a bound peptide through a conserved residue.

 

    Footnotes
 
1 This study was supported by the promotion of basic research activities for innovative biosciences funded by Bio-Oriented Technology Research Advancement Institution, Japan, and supported by the postdoctoral fellowships for K.A., C.X., and J.M.D. from the Science and Technology Agency of Japan. Back

2 The sequences presented in this article have been submitted to GenBank under accession numbers AF287483-AF287492 and AY044428. Back

3 K.A. and J.M.D. contributed equally to this work. Back

4 Address correspondence and reprint requests to Dr. Mitsuru Ototake, National Research Institute of Aquaculture, Fisheries Research Agency, Tamaki, Mie 519-0423, Japan. E-mail address: ototake{at}affrc.go.jp Back

5 Current address: Department of Veterinary Medicine, College of Bioresource Sciences, Nihon University, Fujisawa, Kanagawa, Japan. Back

6 Abbreviations used in this paper: {beta}2m, {beta}2-microglobulin; TM, transmembrane; UPGMA, unweighted pair-group method with arithmatic mean; CY, cytoplasmic; UTR, untranslated region. Back

Received for publication January 29, 2001. Accepted for publication October 29, 2001.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 References
 

  1. Townsend, A., H. Bodmer. 1989. Antigen recognition by class I-restricted T lymphocytes. Annu. Rev. Immunol. 7:601.[Medline]
  2. Germain, R. N.. 1994. MHC-dependent antigen processing and peptide presentation: providing ligands for T lymphocyte activation. Cell 76:287.[Medline]
  3. Manning, M. J., T. Nakanishi. 1996. The specific immune system: cellular defenses. G. Iwama, and T. Nakanishi, eds. The Fish Immune System 159. Academic Press, San Diego.
  4. Dixon, B., S. H. van Erp, P. N. Rodrigues, E. Egberts, R. J. Stet. 1995. Fish major histocompatibility complex genes: an expansion. Dev. Comp. Immunol. 19:109.[Medline]
  5. Stet, R. J. M., B. Dixon, S. H. M. van Erp, M. C. van Lierop, P. N. S. Rodrigues, E. Egberts. 1996. Inference of structure and function of fish major histocompatibility complex (MHC) molecules from expressed genes. Fish Shellfish Immunol. 6:305.
  6. Rast, J. P., G. W. Litman. 1998. Towards understanding the evolutionary origins and early diversification of rearranging antigen receptors. Immunol. Rev. 166:79.[Medline]
  7. Klein, J., A. Sato. 1998. Birth of the major histocompatibility complex. Scand. J. Immunol. 47:199.[Medline]
  8. Takeuchi, H., F. Figueroa, C. O’hUigin, J. Klein. 1995. Cloning and characterization of class I MHC genes of the zebrafish Brachydanio rerio. Immunogenetics 42:77.[Medline]
  9. Sato, A., D. Klein, H. Sultmann, F. Figueroa, C. O’hUigin, J. Klein. 1997. Class I MHC genes of cichlid fishes: identification, expression, and polymorphism. Immunogenetics 46:63.[Medline]
  10. Sato, A., F. Figueroa, C. O’hUigin, D. N. Reznick, J. Klein. 1996. Identification of major histocompatibility complex genes in the guppy Poecilia reticulata. Immunogenetics 43:38.[Medline]
  11. Grimholt, U., I. Hordvik, V. M. Fosse, I. Olsaker, C. Endresen, O. Lie. 1993. Molecular cloning of major histocompatibility complex class I cDNAs from Atlantic salmon (Salmo salar). Immunogenetics 37:469.[Medline]
  12. Betz, U. A., W. E. Mayer, J. Klein. 1994. Major histocompatibility complex class I genes of the coelacanth Latimeria chalumnae. Proc. Natl. Acad. Sci. USA 91:11065.[Abstract/Free Full Text]
  13. Okamura, K., M. Ototake, T. Nakanishi, Y. Kurosawa, K. Hashimoto. 1997. The most primitive vertebrates with jaws possess highly polymorphic MHC class I genes comparable to those of humans. Immunity 7:777.[Medline]
  14. Miller, K. M., R. E. Withler. 1998. The salmonid class I MHC: limited diversity in a primitive teleost. Immunol. Rev. 166:279.[Medline]
  15. Hashimoto, K., K. Okamura, H. Yamaguchi, M. Ototake, T. Nakanishi, Y. Kurosawa. 1999. Conservation and diversification of MHC class I and its related molecules in vertebrates. Immunol. Rev. 167:81.[Medline]
  16. Hansen, J. D., P. Strassburger, L. Du Pasquir. 1996. Conservation of an {alpha}2 domain within the teleostean world, MHC class I from the rainbow trout Oncorhynchus mykiss. Dev. Comp. Immunol. 20:417.[Medline]
  17. Hansen, J. D., P. Strassburger, G. H. Thorgaard, W. P. Young, L. Du Pasquier. 1999. Expression, linkage, and polymorphism of MHC-related genes in rainbow trout Oncorhynchus mykiss. J. Immunol. 163:774.[Abstract/Free Full Text]
  18. Shum, B. P., R. Rajalingam, K. E. Magor, K. Azumi, W. H. Carr, B. Dixon, R. J. Stet, M. A. Adkison, R. P. Hedrick, P. Parham. 1999. A divergent non-classical class I gene conserved in salmonids. Immunogenetics 49:479.[Medline]
  19. Shum, B. P., K. Azumi, S. Zhang, S. R. Kehrer, R. L. Raison, H. W. Detrich, P. Parham. 1996. Unexpected {beta}2-microglobulin sequence diversity in individual rainbow trout. Proc. Natl. Acad. Sci. USA 93:2779.[Abstract/Free Full Text]
  20. Scott, D., J. Hewittson, J. C. Fraser. 1978. The origin of rainbow trout Salmo gairdneri, Richardson, in New Zealand. Calif. Fish and Game 64:210.
  21. Onozato, H.. 1984. Diploidization of gynogenetically activated salmonid eggs using hydrostatic pressure. Aquaculture 43:91.
  22. J. Sambrook, and E. F. Fritsch, and T. Maniatis, eds. Molecular Cloning: A Laboratory Manual 2nd Ed.1989 Cold Spring Harbor Laboratory Press, Cold Spring Harbor.
  23. Saper, M. A., P. J. Bjorkman, D. C. Wiley. 1991. Refined structure of the human histocompatibility antigen HLA-A2 at 2.6 A resolution. J. Mol. Biol. 219:277.[Medline]
  24. Madden, D. R.. 1995. The three-dimensional structure of peptide-MHC complexes. Annu. Rev. Immunol. 13:587.[Medline]
  25. Parham, P., E. J. Adams, K. L. Arnett. 1995. The origins of HLA-A,B,C polymorphism. Immunol. Rev. 143:141.[Medline]
  26. Potter, T. A., T. V. Rajan, II R. F. Dick, J. A. Bluestone. 1989. Substitution at residue 227 of H-2 class I molecules abrogates recognition by CD8-dependent, but not CD8-independent, cytotoxic T lymphocytes. Nature 337:73.[Medline]
  27. Salter, R. D., R. J. Benjamin, P. K. Wesley, S. E. Buxton, T. P. Garrett, C. Clayberger, A. M. Krensky, A. M. Norment, D. R. Littman, P. Parham. 1990. A binding site for the T-cell co-receptor CD8 on the {alpha}3 domain of HLA-A2. Nature 345:41.[Medline]
  28. Sun, J., D. J. Leahy, P. B. Kavathas. 1995. Interaction between CD8 and major histocompatibility complex (MHC) class I mediated by multiple contact surfaces that include the {alpha}2 and {alpha}3 domains of MHC class I. J. Exp. Med. 182:1275.[Abstract/Free Full Text]
  29. Matsumura, M., D. H. Fremont, P. A. Peterson, I. A. Wilson. 1992. Emerging principles for the recognition of peptide antigens by MHC class I molecules. Science 257:927.[Abstract/Free Full Text]
  30. Klein, J., C. O’hUigin. 1994. The conundrum of nonclassical major histocompatibility complex genes. Proc. Natl. Acad. Sci. USA 91:6251.[Free Full Text]
  31. Flajnik, M. F., Y. Ohta, A. S. Greenberg, L. Salter-Cid, A. Carrizosa, L. Du Pasquier, M. Kasahara. 1999. Two ancient allelic lineages at the single classical class I locus in the Xenopus MHC. J. Immunol. 163:3826.[Abstract/Free Full Text]
  32. Cereb, N., A. L. Hughes, S. Y. Yang. 1997. Locus-specific conservation of the HLA class I introns by intra-locus homogenization. Immunogenetics 47:30.[Medline]
  33. Allendorf, F. W., G. H. Thorgaard. 1984. Tetraploidy and the evolution of salmonid fishes. B. J. Turner, ed. Evolutionary Genetics of Fishes 1. Plenum Press, New York.
  34. Carroll, R. L.. 1988. Vertebrate Paleontology and Evolution Freeman, New York.
  35. Katagiri, T., I. Hirono, T. Aoki, M. Sakai. 1996. Isolation of major histocompatibility complex class I cDNA from pink salmon (Oncorhynchus gorbuscha). Dev. Comp. Immunol. 20:217.[Medline]
  36. Van Erp, S. H., B. Dixon, F. Figueroa, E. Egberts, R. J. Stet. 1996. Identification and characterization of a new major histocompatibility complex class I gene in carp (Cyprinus carpio L.). Immunogenetics 44:49.[Medline]
  37. Shum, B. P., D. Avila, L. Du Pasquier, M. Kasahara, M. F. Flajnik. 1993. Isolation of a classical MHC class I cDNA from an amphibian: evidence for only one class I locus in the Xenopus MHC. J. Immunol. 151:5376.[Abstract]
  38. Guillemot, F., A. Billault, O. Pourquie, G. Behar, A. M. Chausse, R. Zoorob, G. Kreibich, C. Auffray. 1988. A molecular map of the chicken major histocompatibility complex: the class II {beta} genes are closely linked to the class I genes and the nucleolar organizer. EMBO J. 7:2775.[Medline]
  39. Koller, B. H., H. T. Orr. 1985. Cloning and complete sequence of an HLA-A2 gene: analysis of two HLA-A alleles at the nucleotide level. J. Immunol. 134:2727.[Abstract]
  40. Miller, K. M., R. E. Withler, T. D. Beacham. 1997. Molecular evolution at MHC genes in two populations of chinook salmon Oncorhynchus tshawytscha. Mol. Ecol. 6:937.[Medline]
  41. Shum, B. P., L. Guethlein, L. R. Flodin, M. A. Adkison, R. P. Hedrick, R. B. Nehring, R. J. Stet, C. Secombes, P. Parham. 2001. Modes of salmonid MHC class I and II evolution differ from the primate paradigm. J. Immunol. 166:3297.[Abstract/Free Full Text]



This article has been cited by other articles:


Home page
Lab AnimHome page
U Grimholt, R Johansen, and A J Smith
A review of the need and possible uses for genetically standardized Atlantic salmon (Salmo salar) in research
Lab Anim, April 1, 2009; 43(2): 121 - 126.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
B. Woelfing, A. Traulsen, M. Milinski, and T. Boehm
Does intra-individual major histocompatibility complex diversity keep a golden mean?
Phil Trans R Soc B, January 12, 2009; 364(1513): 117 - 128.
[Abstract] [Full Text] [PDF]


Home page
J. Immunol.Home page
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]


Home page
J. Immunol.Home page
K. E. Magor, B. P. Shum, and P. Parham
The {beta}2-Microglobulin Locus of Rainbow Trout (Oncorhynchus mykiss) Contains Three Polymorphic Genes
J. Immunol., March 15, 2004; 172(6): 3635 - 3643.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow Request Permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Aoyagi, K.
Right arrow Articles by Nakanishi, T.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Aoyagi, K.
Right arrow Articles by Nakanishi, T.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS