|
|
||||||||



*
National Research Institute of Aquaculture, Fisheries Research Agency, Tamaki, Mie, Japan;
Department of Microbiology and Infection, College of Veterinary Medicine, China Agricultural University, Beijing, Peoples Republic of China;
Nagano Prefectural Fisheries Experimental Station, Akashina, Nagano, Japan; and
Institute for Comprehensive Medical Science, Fujita Health University, Toyoake, Aichi, Japan
| Abstract |
|---|
|
|
|---|
3 and transmembrane domains
is unprecedented. | Introduction |
|---|
|
|
|---|
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,
2-microglobulin
(
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
1 and
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
1 and
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
1 and
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 |
|---|
|
|
|---|
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. 1
. For initial amplification of the
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-
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-
2r, 5'-ACTGTCGACYTTNARCCAYTCDATRCA
(reverse).
|
2r1 derived
from the Onmy-UBA*401
2 region,
5'-CTCATCATCCCACTCACAACCGTACAT (reverse), and anchor primer
pVsq, 5'-TGGCGGCCGCTCTAGAACTAGTGGATC derived from vector
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 p13'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 p53'r, 5'-AACATACCTGTAAGAATGAGAAGT (reverse) were used. Primers for specific detection of Onmy-UBA*601 or *701 expression were pI-5'f (forward) and p673'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
2 probe of Onmy-UBA*101, pr1-
2,
primers, p1-
2f, 5'-GGATGATGACACTG GAGCC (forward) and
p1-
2r, 5'-AATGTTTTCTCCGCTCCTCAT (reverse), and template,
Onmy-UBA*101 cDNA clone; the 181-bp
2 probe of
Onmy-UBA*201, pr2-
2, primers,
p2-
2f, 5'-GGATGATGAGGCTG GAGCC (forward) and
p2-
2r, 5'-CTGGGTGCAGTAGTTTTTTC (reverse), and template,
Onmy-UBA*201 cDNA clone; the 181-bp
2 probe of
Onmy-UBA*301, pr3-
2, primers,
p3-
2f, 5'-GGATGATGAGACTG GAGTC (forward) and
p3-
2r, 5'- CTGGGTGATGTAGTTTTTC (reverse), and template,
Onmy-UBA*301 cDNA clone; the 220-bp
2 probe of
Onmy-UBA*401, pr4-
2, primers,
p4-
2f, 5'-GAGGATGTACGGTT GTGAGT (forward) and
p4-
2r2, 5'- CAGCCACTCAATGCACGTC (reverse), and template,
Onmy-UBA*401 cDNA clone; the 170-bp
2 probe of
Onmy-UBA*501, pr5-
2, primers,
p5-
2f, 5'-GAATGATGAGACTG GTGCC (forward) and
p5-
2r, 5'-AATGTTTATCCCGCTCAGTCAT (reverse), and template,
Onmy-UBA*501cDNA clone; the 189-bp
2 probe of
Onmy-UBA*601, pr6-
2, primers,
p67-
2f,5'-GCTGATGTATGGCT GTG (forward) and
p67-
2r, 5'-TTCTTCTGAAAGTTAGCCT (reverse), and template,
Onmy-UBA*601 cDNA clone; the 189-bp
2 probe of
Onmy-UBA*701, pr7-
2, primers,
p67-
2f and p67-
2r, and template,
Onmy-UBA*701 cDNA clone; the 241-bp
3 probe of
Onmy-UBA*401, pr4-
3, primers,
p4-
3f, 5'-GTCCCTCCGTCAGTGTCTCTG (forward) and
p4-
3r, 5'-TGACTTGAACCACACACTGATACT (reverse), and
template, Onmy-UBA*401 cDNA clone; the 377-bp 3' UTR probe
of Onmy-UBA*601, pr63' UTR, primers,
p673'f, 5'-GTATACTTACACA CACATCAG (forward) and
p673'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
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
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
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
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-
2f and pMix-
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 manufacturers 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 |
|---|
|
|
|---|
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
2 domain motifs in the classical MHC class I genes of
other vertebrates. The amplified fragment, denoted
Onmy-UBA*401-
2, showed homology with the
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
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. 1
and 2
, discussed later).
More investigation was necessary to obtain a clear picture of the
classical MHC class I genes of the rainbow trout.
|
1, and
2 domain sequences) with the subsequently published rainbow trout
MHC class I sequence Onmy-UCA*C32 (16) showed
sequence differences in all domains (Fig. 1
1,
2, and
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 I
2 or 3' UTR (Fig. 1
|
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. 1
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 I
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. 1
).
Analysis of deduced amino acid sequences
Fig. 2
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. 2
, 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. 2
). 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. 2
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
3 domain (26, 27). Although not identical to those in humans, the Onmy-UBA
sequences all contain an acidic stretch at this position (Fig. 2
). Of
the three residues (Q114, D121, and E127) in the
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
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. 2
).
An additional indication that the Onmy-UBA sequences represent
classical MHC class I molecules is the distribution of variability.
Whereas the
3 regions of the different sublineages share >73%
amino acid identity, within the
1 and
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. 2
)
(29). The average number of detected amino acids in the
1 and
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
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
3 domain to
the membrane.
Whereas the
2 region of the sublineages Onmy-UBA*1, *2, *3, *4, and
*5 are more similar to those detected in other salmonids, the
2
regions of Onmy-UBA*6 and *7 are more closely related to those detected
in the cyprinid sequences depicted in Fig. 2
. They will be discussed
below as type I and type II
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
1 and
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
1 and
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. 3
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 II
and Fig. 3
. Within the groups I-VII of
1, and I-II of
2, the members share >67% amino acid identity
with each other. Within the groups I-V of the
3 region, the members
share >88% amino acid identity. Most trout
1 sequences belong to
group
1-I. Group
1-IV, containing Onmy-UBA*401, is very different
from group
1-I (±50% amino acid identity). Group
1-V contains
both a cyprinid, Brre-UAA*01, and a salmonid, Onmy-UA*A4.11, sequence.
For the
2 region, only two groups were detected. As mentioned before
(16), some trout sequences show high similarity to
cyprinid
2 domains. The division into groups I-III of the
3
domain shows the existence of trans-specific lineages
between salmonids. Group
3-II contains a small extension at the end
of the
3 domain, and group
3-III a large extension (see also Fig. 2
). 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 II
schematically depicts the organization of the
sequences of Fig. 3
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
1 and
2 groups seem to have been an
important cause of polymorphism. For example, the sequences
Onmy-UBA*101 and *701 have a completely identical
1 region, but a
very different
2 region, with differences beginning at the domain
border (Fig. 1
). And Onmy-UBA*601 and Onmy-UA*A5.1 (17)
have an identical
2 domain, but a very different
1 domain.
|
|
Even for a homozygous clonal fish, more than nine bands can be
detected in Southern blot analysis when using a
3 probe (data not
shown), indicating a number of similar genes or pseudogenes. Therefore,
the
3 probe was not suited to identify the Onmy-UBA
locus, and locus identification was performed with
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
2 region
were made. The hybridization pattern of all type I
2 probes (from
Onmy-UBA*101, *201, *301, *401, or *501) was identical; Fig. 4
I was obtained with probe
pr5-
2 derived from Onmy-UBA*501. The
hybridization pattern of all type II
2 probes (from
Onmy-UBA*601 or *701) was also identical; Fig. 4
II was obtained with probe pr6-
2 from
Onmy-UBA*601. The patterns indicate that under the washing
conditions used, all type I
2 probes can recognize all genomic
fragments with a type I
2 region, and all type II
2 probes can
recognize all type II
2 regions in the genome; there seemed to be no
cross-reaction between the two types of
2 regions.
|
2 probes only hybridized to one restriction fragment for
homozygous fish; in outbred fish, either none, one, or two fragments
were detected (Fig. 4
2 region could
be detected, and also Southern blot analysis could not detect such a
region. The type II
2 probes detected one fragment of 2.4 kb in
every fish (Fig. 4
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. 4
II, we tried to amplify novel
2 sequences
from genomic DNA by the use of degenerate primers. We succeeded in
isolation of a new
2 DNA fragment and named its gene
Onmy-UFA. As the new
2 sequence can be classified as a
type II
2, it should correspond to the 2.4-kb fragment in Fig. 4
II. It turned out to be closely related to reported
sequences such as Onts-HN-1 and Sasa-UA*8 (Fig. 5
). 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. 5
). 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. 4I
I), and its expression cannot be detected (see below, Fig. 6
), in agreement with the mutational
defects in the sequence.
|
|
Ubiquitous tissue expression is an indication for the classical
nature of MHC class I genes (30). Fig. 6
shows the
hybridization pattern of probe pr4-
2 (
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
3 region of
Onmy-UBA*401, pr4-
3, the same expression
pattern for different AA1 tissues as with pr4-
2 was
detected (data not shown). With the pr6-
2 and
pr7-
2 probes (
2 type II), no hybridization signal
could be detected (data not shown). As the
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
3 probe pr4-
3 indeed
detected major transcripts of 2.12.9 kb (Fig. 7
A). 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
2 probe pr4-
2 reacted more specifically
than in Southern blot analysis, and recognized transcripts of the
sublineages Onmy-UBA*2, *3, and *4 (Fig. 7
B). The
minor and major transcripts of fish AB2 (Fig. 7
A) are both
recognized by probe pr4-
2 (Fig. 7
B),
indicating their similarity. The type II
2 probe pr6-
2
recognized two transcripts for sublineages Onmy-UBA*6 and
probably Onmy-UBA*7 (Fig. 7
C; 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-
2 probe was confirmed
by probe pr63' UTR derived from the 3' UTR of
Onmy-UBA*601 (Fig. 7
D). The recognized bands in
Fig. 7
, BD, were identical to those recognized by the
3
probe (Fig. 7
A).
|
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 I
); 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
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. 4
); 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
3 probe)
whose transcripts exhibited sublineage-specific variations in length
(Fig. 7
). 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
3 domain could be explained by
our isolated genes. While intralocus homogenization outside the
1
and
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
3, TM/CY, and 3' UTR
(Fig. 1
). 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
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
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. 1
). We could observe
polymorphism at the amino acid positions that could interact with a
bound peptide (Figs. 1
and 2
). 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
1 or
2 domains were
impressively observed in the peptide-binding domains (Table II
). 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
1/
2 domain lineages. In Table II
, 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 II
, 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
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).
Hansens group (17) also reported two sequences with the
group V lineage of the
1 domain (Onmy-UA*A4.11 and Onmy-UAA*KD1.5)
from rainbow trout. We noticed that this lineage corresponds to that of
the
1 domain of zebrafish Brre-UAA*01 (Fig. 3
and Table II
),
although not mentioned previously (17). We expected that
the sequences with the group V lineage in the
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
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
1/
2 domain lineages though, the sequence
Brre-UAA*01 exhibits a V/II combination. Considering the mosaic pattern
of the
1 and
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
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
2
domain also shares the sequence characteristics with the
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
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. 2
), 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
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
1 domain and those
related to the group I in the
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
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
1 and
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
1 and the
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
3/TM boundary
region of rainbow trout MHC class I genes is noteworthy.
|
|
|
| Footnotes |
|---|
2 The sequences presented in this article have been submitted to GenBank under accession numbers AF287483-AF287492 and AY044428. ![]()
3 K.A. and J.M.D. contributed equally to this work. ![]()
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 ![]()
5 Current address: Department of Veterinary Medicine, College of Bioresource Sciences, Nihon University, Fujisawa, Kanagawa, Japan. ![]()
6 Abbreviations used in this paper:
2m,
2-microglobulin; TM, transmembrane; UPGMA, unweighted pair-group method with arithmatic mean; CY, cytoplasmic; UTR, untranslated region. ![]()
Received for publication January 29, 2001. Accepted for publication October 29, 2001.
| References |
|---|
|
|
|---|
2 domain within the teleostean world, MHC class I from the rainbow trout Oncorhynchus mykiss. Dev. Comp. Immunol. 20:417.[Medline]
2-microglobulin sequence diversity in individual rainbow trout. Proc. Natl. Acad. Sci. USA 93:2779.
3 domain of HLA-A2. Nature 345:41.[Medline]
2 and
3 domains of MHC class I. J. Exp. Med. 182:1275.
genes are closely linked to the class I genes and the nucleolar organizer. EMBO J. 7:2775.[Medline]
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
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] |
||||
![]() |
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] |
||||
![]() |
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] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |