Open Access

A genome-wide survey of Major Histocompatibility Complex (MHC) genes and their paralogues in zebrafish

BMC Genomics20056:152

DOI: 10.1186/1471-2164-6-152

Received: 05 June 2005

Accepted: 04 November 2005

Published: 04 November 2005

Abstract

Background

The genomic organisation of the Major Histocompatibility Complex (MHC) varies greatly between different vertebrates. In mammals, the classical MHC consists of a large number of linked genes (e.g. greater than 200 in humans) with predominantly immune function. In some birds, it consists of only a small number of linked MHC core genes (e.g. smaller than 20 in chickens) forming a minimal essential MHC and, in fish, the MHC consists of a so far unknown number of genes including non-linked MHC core genes. Here we report a survey of MHC genes and their paralogues in the zebrafish genome.

Results

Using sequence similarity searches against the zebrafish draft genome assembly (Zv4, September 2004), 149 putative MHC gene loci and their paralogues have been identified. Of these, 41 map to chromosome 19 while the remaining loci are spread across essentially all chromosomes. Despite the fragmentation, a set of MHC core genes involved in peptide transport, loading and presentation are still found in a single linkage group.

Conclusion

The results extend the linkage information of MHC core genes on zebrafish chromosome 19 and show the distribution of the remaining MHC genes and their paralogues to be genome-wide. Although based on a draft genome assembly, this survey demonstrates an essentially fragmented MHC in zebrafish.

Background

The human Major Histocompatibility Complex (MHC) is a gene-dense region on chromosome 6p21.3, and comprises a group of genes that are involved functionally with the adaptive and innate immune system. From centromere to telomere, it is divided into five regions: extended class II, classical class II, class III, classical class I and extended class I. The classical MHC contains 224 genes, many of which are pseudogenes, and with every two in five expressed genes having a potential immune function and a role in disease resistance [1], this region has become a focal locus in comparative genomics. Antibody and T cell mediated immune responses against invading pathogens are initiated through MHC class I and class II molecules [2]. These main components are not only missing from invertebrates, but are also not present in primitive jawless fish, such as hagfish and lamprey [3]. MHC class I and II molecules do, however, exist in all jawed vertebrates, including the cartilaginous fish. The gene loci in the class III region encode a variety of proteins with both immune and non-immune functions [4].

Genomic sequences encompassing the MHC are currently available for human, chimpanzee, macaque, rat, mouse, cat, pig, horse, quail, chicken, frog, teleost fish and shark [5]. While the genomic architecture of the mammalian MHC is conserved, the number of genes between species can differ greatly. In birds, for example, the minimal essential MHC in chicken consists of 19 genes [6] while gene duplications have expanded this region in quail [7] and sparrows [8]. The linkage of class I, II, and III region genes can be traced back to cartilaginous fish, which are the earliest jawed vertebrates known to have diverged from a common ancestor with humans [9]. MHC loci, however, do not always exist in a single tightly linked cluster as generally observed in mammals. A large scale inversion has separated the class IIB cluster from the MHC-linked class IIA cluster in cattle [10]. The class I and class II loci in zebrafish (Brachydanio rerio) are found on different chromosomes [11, 12]. A similar organisation is present in trout, stickleback, common guppy and cichlid fish [13, 14]. This observation demonstrates that the separation of the MHC class I and class II loci is characteristic of teleost fish, which represent half of all vertebrates. Since the genes of the immune system were present in the common ancestor of tetrapods and teleosts, the differences in their genomic organisation may be the result of lineage-specific chromosomal events such as duplications, inversions, deletions and translocations.

The genome of the zebrafish is currently being sequenced at the Wellcome Trust Sanger Institute. The availability of sequence data will allow an insight into the understanding and evolution of the immune system in fish. To date, small regions containing major histocompatibility genes have been mapped by radiation hybrid mapping and sequencing of select genomic clones [15, 16]. With the ongoing whole-genome sequencing efforts and the compilation of physical maps, it is now possible to examine larger genomic regions and assess the degree of shared synteny between mammals and fish on a genome-wide scale. Here we have utilised a comprehensive whole-genome shotgun assembly data set to fully analyse the zebrafish loci that are related to the human MHC. The human genome contains at least three regions that are paralogous to the MHC [17]. These are thought to be the result of two rounds of duplications that occurred early during vertebrate evolution. In this study, we examine the genome-wide distribution of paralogous genes in zebrafish.

Results and Discussion

Analysis of the MHC class I region in zebrafish

The MHC class I region in human embodies gene clusters coding for HLA class I molecules, histones, solute carriers, vomeronasal receptors, olfactory receptors and zinc fingers [18]. These clusters have undergone a large-scale expansion and each has paralogues throughout the genome. Ten representative members of the extended class I region, which are expressed and have paralogues in the human genome, and more than 30 genes that reside in the classical class I region were chosen for this analysis (Table 1). Genes that were excluded from this study included pseudogenes, zinc fingers with the tripartite motif (TRIM), and RING finger proteins. Approximately half of the human genes were identified in the zebrafish whole-genome assembly, predominantly on chromosome 19. Because multigene families evolve by a birth-and-death process [19], orthologous genes are usually difficult to identify through sequence similarity searches.
Table 1

List of human MHC class I and extended class I (xI) region genes used against the zebrafish whole-genome assembly. Genes with the suffix "like" are either gene fragments and/or highly similar to their human counterparts. Genes not identified by sequence similarity searches are marked as not found (NF). The official gene nomenclature for zebrafish is shown in brackets.

MHC CLASS I REGION GENES

HUMAN

ZEBRAFISH

 

name

accession

location

location

location

location

xI

HFE

NM_000410

NF

   

xI

HMGN4

NM_006353

NF

   

xI

PRSS16

NM_005865

NF

   

xI

POM121L2

NM_033482

NF

   

xI

GPX5

NM_001509

NF

   

xI

RFP

NM_006510

NF

   

xI

MAS1L

NM_052967

NF

   

xI

UBD

NM_006398

NF

   

xI

GABBR1

NM_021905

15 [gabbr1]

19 [gabbr1]

  

xI

MOG

NM_002433

NF

   

I

HLA*

 

19 [mhc1uea]

19 [mhc1uaa]

19 [mhc1ufa]

19 [mhc1uda]

   

s_1723 [mhc1ze]

1 [mhc1-like] (X3)

1 [zgc:64115]

 

I

HCG9

NM_005844

NF

   

I

PPP1R11

NM_021959

NF

   

I

RPP21

NM_024839

NF

   

I

GNL1

NM_005275

16 [gnl1]

   

I

PRR3

NM_025263

NF

   

I

ABCF1

NM_001090

19 [zgc:85667]

   

I

PPP1R10

NM_002714

19 [ppp1r10]

   

I

MRPS18B

NM_014046

19 [mrps18b]

   

I

C6orf134

NM_024909

19 [c6orf134]

   

I

C6orf136

NM_145029

19 [c6orf136]

   

I

DHX16

NM_003587

NA5816 [dhx16]

15 [dhx16]

  

I

NRM

NM_007243

NA1826 [nrm]

   

I

MDC1

NM_014641

NF

   

I

TUBB

NM_178014

20 [tubb5]

14 [tubb2]

23 [tubb2]

 

I

FLOT1

NM_005803

19 [flot1]

16 [flot1b]

  

I

IER3

NM_003897

NF

22 [ier5]

  

I

DDR1

NM_013994

16 [ddr1-like]

4 [ddr2]

  

I

GTF2H4

NM_001517

19 [zgc:77721]

   

I

VARS2L

NM_020442

19 [vars2-like]

   

I

DPCR1

NM_080870

NF

   

I

C6orf15

NM_014070

NF

   

I

CDSN

NM_001264

NF

   

I

PSORS1C1

NM_014068

NF

   

I

PSORS1C2

NM_014069

NF

   

I

C6orf18

NM_019052

NF

   

I

TCF19

NM_007109

1 [tcf19-like]

   

I

POU5FI

NM_002701

16 [pou5f1]

   

I

MICA

NM_000247

NF

   

I

HCP5

NM_006674

NF

   

I

MICB

NM_005931

NF

   
Chromosome 19 harbours the largest number of MHC genes in zebrafish (Figure 1). The core region, comprising genes encoding the classical peptide presenting class I molecules (mhc1uea, mhc1uaa, mhc1ufa, mhc1uda), proteasome subunits (psmb8, psmb9a, psmb11, psmb10), ATP-binding cassettes (abcb3, abcb3l1, abcb3l2) and tapasin (tpsn), cluster tightly together and span approximately 175 kb of DNA. Over 5 Mb downstream of these loci are several other genes that span the human MHC region, although these are not all found adjacent to each other. In humans, MHC class I molecules are expressed on the surface of most cells and are recognized by CD8+ cytotoxic T cells. They are heterodimers composed of an alpha subunit, encoded by the MHC, and a beta subunit, or b-2 microglobulin (B2M), which maps to chromosome 15q21–q22.2. In zebrafish too, the b2m gene is found located in a different part of the genome on chromosome 13. Previously unmapped class I ZE lineage genes [20] and novel MHC class I antigen fragments have been identified on chromosome 1 and scaffold 1723. The additional functional class I lineage (mhc1ze) is also present in two other cyprinids (carp and barbus), and not linked to other class I and II genes.
Figure 1

Map of major histocompatibility genes and their paralogues in zebrafish (not to scale). Only chromosomes and two unmapped contigs (NA17761 and NA17767) that harbour MHC-related genes are shown. Orthologues or paralogues of human MHC class I region genes are shown in red, of class II region genes in blue, and class III region genes in green. Genes with the suffix "L" for like are either gene fragments and/or highly similar to their human counterparts. Genes that have more than one copy in the genome are shown by a greater than (>) symbol. Similarities and differences between the whole-genome assembly and mapping information present in the ZFIN database are shown by plus (+) and minus (-) signs, respectively. The core MHC has been compiled from five genomic clones spanning this region [EMBL:AL672216, EMBL:AL672151, EMBL:AL672164, EMBL:AL672176].

Analysis of the MHC class II region in zebrafish

MHC class II molecules are heterodimers comprising A (alpha) and B (beta) chains that present peptides to CD4+ T cells via the endosomal pathway. One class IIA (mhc2daa) and six class IIB (mhc2dab, mhc2dbb, mhc2dcb, mhc2ddb, mhc2deb, mhc2dfb) genes have previously been identified by screening a zebrafish genomic BAC library [21]. Only the mhc2dab and mhc2daa genes are known to be expressed [22]. They are closely linked and were identified on chromosome 8 (Table 2). Analysis of approximately 1 Mb of contiguous DNA surrounding the functional class II region in zebrafish demonstrates the presence of 23 flanking genes mapping to various human chromosomal locations, including two or more genes mapping to human 6q (C6orf117, LOC557721, ppil4), 12q (slc15a4, KIAA1944), 20q (rnpc2, slc12a5), 22q (slc7a4, sf3a1) and Xp (pqbp1l, t541, jarid1c). The lack of genes mapping to the human MHC, in addition to the low gene density of this region, indicates that the functional zebrafish class II region is the result of a translocation event [12].
Table 2

List of human MHC class II and extended class II (xII) region genes used against the zebrafish whole-genome assembly. Genes with the suffix "like" are either gene fragments and/or highly similar to their human counterparts. Genes not identified by sequence similarity searches are marked as not found (NF). The official gene nomenclature for zebrafish is shown in brackets.

MHC CLASS II REGION GENES

HUMAN

ZEBRAFISH

 

name

accession

location

location

location

location

II

C6orf10

NM_006781

NF

   

II

BTNL2

NM_019602

NF

   

II

HLA-D*

 

8 [mhc2dab]

8 [mhc2daa]

19 [mhc2dcb-rs]

s_2399 [mhc2dfb]

   

5 [mhc2-like]

NA14232 [mhc2dbb]

NA6696 [mhc2ddb]

NA6696 [mhc2-like]

   

20 [mhc2-like]

15 [mhc2-like]

  

II

TAP2

NM_018833

19 [abcb3]

19 [abcb3l1]

19 [abcb3l2]

 

II

PSMB8

NM_004159

19 [psmb8]

19 [psmb10]

19 [psmb11]

19 [psmb9a]

   

2 [psmb5]

s_1793 [psmb7]

  

II

TAP1

NM_000593

9 [tap1]

   

II

PSMB9

NM_002800

19 [psmb8]

19 [psmb10]

19 [psmb11]

19 [psmb9a]

   

2 [psmb5]

s_1793 [psmb7]

  

II

BRD2

NM_005104

19 [brd2]

9 [brd2-like]

s_549 [brd3]

5 [zgc:77289]

   

14 [brd8]

14 [brdt-like]

s_277 [brd4]

 

xII

COL11A2

NM_080679

19 [col11a2]

24 [col11a1]

  

xII

RXRB

NM_021976

17 [rxrb]

NA16779 [rxrb]

21 [rxra]

5 [rxra]

   

9 [rxrg]

   

xII

SLC39A7

NM_006979

19 [ke4]

   

xII

HSD17B8

NM_014234

19 [fabgl]

   

xII

RING1

NM_002931

NF

2 [rnf2]

  

xII

VPS52

NM_022553

19 [vps52]

   

xII

RPS18

NM_022551

NA303 [rps18]

   

xII

B3GALT4

NM_003782

NF

22 [zgc:76904]

s_1722 [b3galt1]

 

xII

C6orf11

NM_005452

s_2182 [bing4]

   

xII

HKE2

NM_014260

14 [zgc:66282]

NA8827 [zgc:66282]

  

xII

RGL2

NM_004761

23 [rgl2]

NA14862 [rgl2-like]

NA7773 [zgc:77299]

NA12682 [ralgds]

xII

TAPBP

NM_003190

19 [tpsn]

   

xII

ZNF297

NM_005453

19 [znf297]

   

xII

DAXX

NM_001350

19 [daxx]

   

xII

LYPLA2L

CAB63783

19 [lypla2l-like]

   

xII

KIFC1

NM_002263

19 [kifc1]

   

Mapping to chromosome 19, approximately 22 Mb telomeric of the class I region, is the class IIB gene mhc2dcb-rs (Q95HJ7), with a similar sequence on contig NA4006. Also within this segment of DNA there are two predicted class IIA chain-encoding gene fragments, consisting of only the alpha2 domain and cytoplasmic tail-encoding parts. It is therefore unlikely that this locus is functional. Assuming that these gene fragments are not due to errors in the whole-genome assembly, it is evidence for a linkage of remnants of MHC-related class II genes with the core MHC region containing the class I peptide presenting genes.

Previously identified mhc2dbb and mhc2dfb genes are predicted to be located in currently unmapped contigs, NA14232 and scaffold 2399. However, mhc2deb (found in clone U08874) was not identified in the current assembly, and may be attributed to gaps in the sequence data or allelic or haplotypic differences. Several fragmented sequences resembling class IIA and/or B chain-encoding genes were also identified in contigs NA17244 and NA15003, and chromosomes 5, 15 and 20. Only contig NA6696 harbours putative complete class IIA and class IIB genes. Another gene found in this 15 kb contig resembles mhc2ddb, consisting of exons 3 and 4 only.

The two TAP transporters in human, namely TAP1 and TAP2, are both located in the MHC class II region, and are closely linked to genes encoding the proteasome subunits and the class II molecules. In zebrafish, however, a TAP1-like sequence is found on chromosomes 9 and three TAP2-like genes are found in the class I region [15] on chromosome 19. The latter, named abcb3, abcb3l1, abcb3l2, form part of the zebrafish core MHC region. The trout tap1 gene is not linked to the major class IA region either [23]. Likewise, in Fugu, the tap1 gene is found on an isolated scaffold that is not linked to the main class I region [24].

Analysis of the MHC class III region in zebrafish

The human class III is the most gene dense region within the MHC containing few or no pseudogenes. Unlike the class I and class II loci that are evolutionary and functionally related [25], the class III region genes are not. The class III, however, include immune-related genes such as those encoding complement components, tumour necrosis factors and heat shock proteins. The search for the MHC class III region in zebrafish was first initiated by the identification of the BF/C2 gene using degenerate PCR [26]. This was then followed by the identification of several, but not all, of the zebrafish homologues of the human MHC class III region genes [16].

In the zebrafish, the MHC class III genes are found distributed throughout their genome (Table 3). Several zebrafish genes (zgc:63773, tnf, bat2-like, ck2b, ddah1, skiv2l, C6orf31-like, ppt2, bat3-like) map to chromosome 19, which also encompasses the largest stretch of MHC-related genes. Two other unmapped scaffolds also harbour several class III genes: neu1, zgc:64108, bat8, zbtb12 are found in NA17767; and bf/c2, rdbp, hsp70, skiv2l are in scaffold NA17761. These are syntenic to the human MHC class III region with the conservation of both gene order and content (Figure 1).
Table 3

List of human MHC class III region genes used against the zebrafish whole-genome assembly. Genes with the suffix "like" are either gene fragments and/or highly similar to their human counterparts. Genes not identified by sequence similarity searches are marked as not found (NF). The official gene nomenclature for zebrafish is shown in brackets.

MHC CLASS III REGION GENES

HUMAN

ZEBRAFISH

name

accession

location

location

location

location

BAT1

NM_004640

19 [zgc:63773]

16 [ddx39]

  

ATP6V1G2

NM_130463

NF

5 [atp6v1g1]

  

NFKBIL1

NM_005007

NF

   

TNF

NM_000594

19 [tnf]

NA11949 [tnfa-like]

  

LST1

NM_007161

NF

   

NCR3

NM_147130

NF

   

AIF1

NM_001623

NF

5 [zgc:73179]

  

BAT2

NM_004638

19 [bat2-like]

   

BAT3

NM_004639

1 [bat3]

19 [bat3-like]

  

APOM

NM_019101

NF

   

C6orf47

NM_021184

NF

   

BAT4

NM_033177

s_339 [bat4-like]

   

CSNK2B

NM_001320

19 [ck2b] (X2)

23 [ck2b]

  

LY6G5B

NM_021221

NF

   

LY6G5C

NM_025262

NF

   

BAT5

NM_021160

NA9087 [bat5-like]

   

LY6G6D

NM_021246

NF

   

LY6G6E

NM_024123

NF

   

LY6G6C

NM_025261

NF

   

C6orf25

NM_025260

NF

   

DDAH2

NM_013974

NF

19 [ddah1]

  

CLIC1

NM_001288

18 [zgc:77538]

14 [zgc:92762]

16 [clica]

17 [clic4]

MSH5

NM_002441

15 [msh5]

   

C6orf26

NM_025259

NF

   

C6orf27

NM_025258

1 [c6orf27-like]

   

VARS2

NM_006295

11 [vars2]

   

LSM2

NM_021177

6 [smx5]

   

HSPA1A

NM_005345

NA17761 [hsp70]

3 [hsp70] (X3)

8 [hsp70]

 

C6orf48

NM_016947

NF

   

NEU1

NM_000434

NA17767 [neu1]

   

C6orf29

NM_032794

NA17767 [zgc:64108]

   

BAT8

NM_025256

NA17767 [bat8]

   

ZBTB12

NM_181842

NA17767 [zbtb12]

   

BF/C2

NM_001710

21 [bf]

NA17761 [bf]

  

RDBP

NM_002904

NA17761 [rdbp]

   

SKIV2L

NM_006929

19 [skiv2l]

NA17761 [skiv2l]

  

DOM3Z

NM_032419

NA17766 [dom3z]

   

STK19

NM_004197

18 [stk19-like]

   

C4B

NM_000592

15 [c4]

15 [a2m]

22 [c3]

NA17328 [a2m]

  

NA14479 [c3]

   

CYP21A2

NM_000500

6 [cyp21a2]

   

TNXB

NM_032470

16 [tnx-like]

2 [tnw]

2 [tnr]

5 [tnc]

CREBL1

NM_004381

NF

20 [atf6-like]

  

FKBPL

NM_022110

NF

   

C6orf31

NM_030651

19 [c6orf31-like]

   

PPT2

NM_005155

19 [ppt2]

s_1796 [ppt2]

NA5816 [zgc:55621]

 

EGFL8

NM_030652

NF

   

AGPAT1

NM_032741

NF

1 [agpat3]

13 [agpat4]

 

RNF5

NM_006913

NF

   

AGER

NM_001136

NF

   

PBX2

NM_002586

NA14559 [pbx1a]

11 [pbx1a]

18 [pbx3b]

NA11844 [pbxy]

GPSM3

NM_022107

NF

   

NOTCH4

NM_004557

s_285 [notch3]

5 [notch1]

NA15389 [notch2]

s_1523 [notch3]

Approximately half of the human MHC class III region genes were not identified in the zebrafish assembly. Among these were the Ly6 family members, which may therefore be mammalian-specific. Alternatively, being involved in the immune response, Ly6 genes evolved more rapidly than others, and might have diverged sufficiently to not be recognised by sequence similarity searches. For eight human genes, ATP6V1G2, AIF1, CLIC1, CREBL1, DDAH2, AGPAT1, PBX2 and NOTCH4, the paralogues but not the orthologues of the genes in the human MHC class III region have been identified by sequence similarity in zebrafish. This observation extends to IER3, RING1 and B3GALT4 found in the class I and extended class II regions. The presence or absence of other genes may be attributed to lineage-specific evolution. For example, lamprey [27], zebrafish [26] and medaka [28] possess genes equally similar to both mammalian BF and C2, while a Xenopus clone has clearly been identified as encompassing BF. This BF/C2 ancestral gene has further duplicated in zebrafish [29], copies of which were identified on chromosome 21 and contig NA17761. Similarly, a recent survey of hsp70 genes in Xiphophorus maculates (platyfish) [30] has revealed that a single HSP70 gene gave rise to four distinct groups of genes: mammalian testis-specific HST70, mammalian MHC-linked HSP70, mammalian HSP70B' and the fish HSP70. Human class III HSPA1A, HSPA1B and HSPA1L genes are intronless and encode identical or near-identical proteins. Intronless zebrafish hsp70 genes were identified on chromosome 8 and scaffold NA17761. Three additional copies of HSP70-like genes were identified on chromosome 3, although these contained one, two or three introns in the 3' end of their sequence. The identified zebrafish genes cluster phylogenetically within the fish subgroup, apart from the sequence on scaffold NA17761, which appears to be similar to the mammalian MHC-linked heat shock protein sequences (data not shown).

Construction of in-silico gene maps

The zebrafish MHC gene map (Figure 1) was constructed using primarily mapping information from the whole-genome assembly displayed in the Ensembl platform. The position of at least 49 genes mapped in this survey could be verified by comparison to experimental data stored in the ZFIN database, which also yielded information for genes residing in unmapped contigs. Mapping data from these two sources did not coincide for several genes, including bat3, clica, vars2, c3, pou5f1, rgl2, tubb5, ddx39, col11a1, pbx1a, and stk19-like. These discrepencies may be attributed to the high levels of polymorphisms and regions of misassembly caused by the source DNA used for the whole-genome assembly being collected from over 1000 embryos. Until the genome sequence is complete, it will not be possible to accurately predict the position and number of all human MHC orthologues in zebrafish. There are also a number of genes, notably CD1, MOG and HFE, that could not be identified in the draft assembly used here (these are listed as 'not found' in Tables 1, 2, 3).

There were a number of difficulties associated with assessing the degree of synteny between the zebrafish whole-genome assembly and the human MHC. Although BLAST is a heavily used analysis tool for identifying related sequences, it does not discriminate between large gene family members. To maximise the identification of MHC-related loci in zebrafish, only the highest BLAST scoring sequences were chosen for further analysis (many of which were unique BLAST reciprocal hits), and searches were conducted using blastp or tblastx. Orthology can be confirmed by obtaining mapping data of surrounding genes, with the assumption that groups of syntenic genes would remain in close proximity through evolution, and therefore be maintained in similar segments of DNA. The physical linkage between numerous MHC-related loci is apparent on zebrafish chromosome 19. Genes mapping to zebrafish chromosome 5 have also been observed to be syntenic with human chromosome 9 [31]. Gene blocks, in particular present in contig NA17761 and NA17767, map closely together in the human MHC class III region. This criterion is more difficult to apply when genes are dispersed, as seen in Figure 1. The Ensembl platform [32] used to map the MHC in zebrafish provides a combination of alignment data, genomic location, detailed transcript structures to compare functional domains of orthologous proteins, in conjunction with multi-species comparisons. In combination with further phylogenetic analyses when duplicates of one gene were identified, the distinction between orthologues and paralogues was ascertained. Nevertheless, it is very difficult to assign orthologous comparative relationships with multigene families until all members of the family have been sequenced in both organisms. This might be the case in human, but the zebrafish genome is yet to be completed.

Duplicated genes and MHC paralagous regions

In human, it was observed that four regions, the MHC on 6p21.3 as well as 9q33–q34, 1q21–q25/1p11–p32 and 19p13.1–p13.4 are paralogous regions [33] that share members of the same gene family. Further analysis has shown that paralogues of human MHC genes are also scattered essentially over all chromosomes [18]. Likewise, genome-wide duplications have recently been examined in zebrafish [31], confirming that an extra round of a whole-genome duplication event occurred early in the teleost lineage after it split from the tetrapod lineage. Evidence that paralogous genes exist in zebrafish became apparent when many duplicates associated with the MHC, and mainly the class III region, were identified (Figure 1). Here we discuss five MHC-encoded gene families in more detail.

The Notch gene family members encode evolutionary transmembrane receptors that regulate cell fate determination. Four Notch paralogues (NOTCH1-4) have been identified in human, and only three of their orthologues are found in zebrafish (Table 3). The orthologue of human NOTCH4, found in the MHC, is absent in the fish whole-genome assembly, indicative that gene loss has followed the block duplication events in teleosts. This is also seen in Fugu [24] and may be due to it being the most divergent member of the Notch family [34]. Two duplicates of notch3 are found in scaffolds 1523 and 285. Until the assembly is complete, it is not clear whether there are two copies of this gene or if this is due to the status of the current assembly. Only one surrounding rdh8-like sequence is found in common between the two scaffolds. In addition, notch1 and notch2 have been mapped to chromosome 5 and contig NA15389.

Retinoid receptors are soluble nuclear proteins belonging to the steroid/thyroid hormone receptor superfamily of transcriptional regulators. The RXR subfamily consists of three polypeptide chains, namely alpha, beta and gamma, encoded by separate loci. The human RXRB gene is found within the MHC, while the RXRA and RXRG paralogues are located on chromosomes 9 and 1, respectively. Five rxr loci have been identified in the zebrafish assembly in contig NA16779 and chromosomes 5, 9, 17 and 21. Two semi-orthologues of human rxrb appear to have arisen from a fish-specific duplication, and duplicates of rxra are also present in the zebrafish genome. Mapping data from the current assembly in comparison to the ZFIN linkage maps are contradictory and may highlight problems in both approaches to mapping genes. Interestingly, neither of the rxrb sequences map to chromosome 19 in the Zv4 whole-genome assembly as originally thought [35].

The human PBX2 gene encodes a homeodomain-containing protein. Three paralogues are located within chromosomes 1q23, 9q33 and 19p12, named PBX1, PBX3 and PBX4, respectively. Four PBX-like sequences were identified in the zebrafish: pbx1a on scaffold NA14559 and chromosome 11; pbx3b on chromosome 18; and pbxy on scaffold NA11844. These are related to human PBX1, PBX3 and PBX4, respectively [36]. The orthologous sequence to human PBX2, which is present in the MHC, has not been identified in the current assembly.

The complement component C4 gene encodes a protein that plays a central role in the innate immune response. Structurally, C4, C5 and C3 belong to the α2 macroglobulin (A2M) protein family and are derived from a single common ancestor. In human, they are found on four different chromosomes: 6p21 (C4), 9q33 (C5), 19p13 (C3) and 12p13 (A2M), and a similar situation exists in all tetrapods studied to date. In zebrafish, the a2m and c4 genes are both found on chromosome 15 as previously shown [37]. Multiple copies of a2m have also been observed in zebrafish with a duplicate being present in NA17328. Two copies of the c3 gene were identified on chromosome 22 and NA14479.

Six members of the chloride intracellular channel (CLIC) gene family (CLIC1–CLIC6) have been described in humans. They are involved in chloride ion transport within various subcellular compartments [38]. Four clic genes were identified in zebrafish on chromosomes 14, 16, 17 and 18. On phylogenetic analysis, they cluster with mammalian CLIC2, CLIC3, CLIC4 and CLIC5 genes, respectively. The orthologue of the human MHC-embedded CLIC1 was not found.

Conclusion

Comparative genomics reveals that the organisation of the MHC in more distantly-related organisms varies from the human model [39]. Teleost fish are particularly unusual in their organisation in comparison to mammals, chicken, Xenopus and shark in that their class I and class II loci are found on different chromosomes [14]. The core class I region in zebrafish, medaka, Fugu and rainbow trout [13, 4042] comprises genes involved in class I peptide presentation and processing: the classical class I molecules, the immunoproteasome subunits, ATP-binding cassette transporters and tapasin. The whole-genome shotgun data for zebrafish has allowed a genome-wide analysis of major histocompatibility genes and their paralogues, and highlights that one or more copies exist of MHC-related genes in fish as in humans. The results obtained thus far extend the linkage information regarding major histocompatibility genes on chromosome 19 in zebrafish and the classical mammalian MHC, and further supports previous findings that the functional class II loci are found on a different chromosome. The distribution of the remaining MHC-related genes and their paralogues is genome-wide, confirming a fragmented MHC in fish.

Methods

The genome of the zebrafish is currently being sequenced at the Sanger Institute [43]. The analysis was carried out on the fourth Ensembl whole-genome assembly Zv4, released on September 2004 [32]. This assembly integrates the whole-genome shotgun assembly with data from the physical map [44]. Protein sequences encoded in 117 human MHC genes were chosen for BLAST sequence similarity searches against the zebrafish Ensembl assembly. Previously identified MHC-related zebrafish cDNAs [45] were also used in the analysis. Gene annotations were verified using VEGA [46]. A number of potentially duplicated genes were identified. Their predicted amino acid sequences were then aligned with other vertebrate gene family members present in the UniProt database [47] using CLUSTAL [48], and the multiple sequence alignments were used for phylogenetic analysis by the neighbour joining method [49] via the PHYLO_WIN interface [50]. Official gene symbols assigned by the Zebrafish Nomenclature Committee (ZNC) and HUGO Nomenclature Committee (HGNC) have been used in the annotation. Zebrafish and human gene names are given in lower and upper case, respectively.

The zebrafish source DNA for the whole-genome assembly was acquired from approximately 1000 five day old embryos. This has resulted in possible misassemblies due to polymorphisms between sequences, different haplotypes, and difficulties with the assembly of duplicated genes and regions, in particular with many highly homologous MHC genes that have arisen by duplication. Where possible, finished clones [e.g. EMBL:AL672216, EMBL:AL672151, EMBL:AL672185, EMBL:AL672164, EMBL:AL672176] were used to study the short-range linkage of genes in preference to the Ensembl whole-genome assembly. Alternative mapping information, for comparison and validity of the in-silico approach used during this study, was obtained from the Zebrafish Information Network (ZFIN) [45] and published scientific literature.

Note added in proof

At least one of the genes stated to be missing in assembly Zv4 is present in assembly Zv5 which was released on October 2005.

Declarations

Acknowledgements

We like to thank all members of the sequencing department and zebrafish project group at the Sanger Institute [43] and K. Belov for comments. We would like to acknowledge J. Klein for providing genomic clones used in this study. JS and SB were supported by the Wellcome Trust.

Authors’ Affiliations

(1)
Wellcome Trust Sanger Institute
(2)
Max-Planck-Institut für Biologie, Abteilung Immunogenetik

References

  1. MHC Sequencing Consortium: Complete sequence and gene map of a human major histocompatibility complex. The MHC sequencing consortium. Nature. 1999, 401: 921-923. 10.1038/44853.View ArticleGoogle Scholar
  2. Parham P: Virtual reality in the MHC. Immunol Rev. 1999, 167: 5-15.PubMedView ArticleGoogle Scholar
  3. Marchalonis JJ, Schluter SF: On the relevance of invertebrate recognition and defence mechanisms to the emergence of the immune response of vertebrates. Scand J Immunol. 1990, 32: 13-20.PubMedView ArticleGoogle Scholar
  4. Milner CM, Campbell RD: Genetic organization of the human MHC class III region. Front Biosci. 2001, 6: D914-26.PubMedView ArticleGoogle Scholar
  5. Kelley J, Walter L, Trowsdale J: Comparative genomics of major histocompatibility complexes. Immunogenetics. 2005, 56: 683-695. 10.1007/s00251-004-0717-7.PubMedView ArticleGoogle Scholar
  6. Kaufman J, Milne S, Gobel TW, Walker BA, Jacob JP, Auffray C, Zoorob R, Beck S: The chicken B locus is a minimal essential major histocompatibility complex. Nature. 1999, 401: 923-925. 10.1038/44856.PubMedView ArticleGoogle Scholar
  7. Shiina T, Shimizu S, Hosomichi K, Kohara S, Watanabe S, Hanzawa K, Beck S, Kulski JK, Inoko H: Comparative genomic analysis of two avian (quail and chicken) MHC regions. J Immunol. 2004, 172: 6751-6763.PubMedView ArticleGoogle Scholar
  8. Freeman-Gallant CR, Johnson EM, Saponara F, Stanger M: Variation at the major histocompatibility complex in Savannah sparrows. Mol Ecol. 2002, 11: 1125-1130. 10.1046/j.1365-294X.2002.01508.x.PubMedView ArticleGoogle Scholar
  9. Terado T, Okamura K, Ohta Y, Shin DH, Smith SL, Hashimoto K, Takemoto T, Nonaka MI, Kimura H, Flajnik MF, Nonaka M: Molecular cloning of C4 gene and identification of the class III complement region in the shark MHC. J Immunol. 2003, 171: 2461-2466.PubMedView ArticleGoogle Scholar
  10. Band M, Larson JH, Womack JE, Lewin HA: A radiation hybrid map of BTA23: identification of a chromosomal rearrangement leading to separation of the cattle MHC class II subregions. Genomics. 1998, 53: 269-275. 10.1006/geno.1998.5507.PubMedView ArticleGoogle Scholar
  11. Bingulac-Popovic J, Figueroa F, Sato A, Talbot WS, Johnson SL, Gates M, Postlethwait JH, Klein J: Mapping of mhc class I and class II regions to different linkage groups in the zebrafish, Danio rerio. Immunogenetics. 1997, 46: 129-134. 10.1007/s002510050251.PubMedView ArticleGoogle Scholar
  12. Kuroda N, Figueroa F, O'HUigin C, Klein J: Evidence that the separation of Mhc class II from class I loci in the zebrafish, Danio rerio, occurred by translocation. Immunogenetics. 2002, 54: 418-430. 10.1007/s00251-002-0473-5.PubMedView ArticleGoogle Scholar
  13. Hansen JD, Strassburger P, Thorgaard GH, Young WP, Du Pasquier L: Expression, linkage, and polymorphism of MHC-related genes in rainbow trout, Oncorhynchus mykiss. J Immunol. 1999, 163: 774-786.PubMedGoogle Scholar
  14. Sato A, Figueroa F, Murray BW, Malaga-Trillo E, Zaleska-Rutczynska Z, Sultmann H, Toyosawa S, Wedekind C, Steck N, Klein J: Nonlinkage of major histocompatibility complex class I and class II loci in bony fishes. Immunogenetics. 2000, 51: 108-116. 10.1007/s002510050019.PubMedView ArticleGoogle Scholar
  15. Michalova V, Murray BW, Sultmann H, Klein J: A contig map of the Mhc class I genomic region in the zebrafish reveals ancient synteny. J Immunol. 2000, 164: 5296-5305.PubMedView ArticleGoogle Scholar
  16. Sultmann H, Sato A, Murray BW, Takezaki N, Geisler R, Rauch GJ, Klein J: Conservation of Mhc class III region synteny between zebrafish and human as determined by radiation hybrid mapping. J Immunol. 2000, 165: 6984-6993.PubMedView ArticleGoogle Scholar
  17. Kasahara M: Genome dynamics of the major histocompatibility complex: insights from genome paralogy. Immunogenetics. 1999, 50: 134-145. 10.1007/s002510050589.PubMedView ArticleGoogle Scholar
  18. Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, Khodiyar VK, Lush MJ, Povey S, Talbot CCJ, Wright MW, Wain HM, Trowsdale J, Ziegler A, Beck S: Gene map of the extended human MHC. Nat Rev Genet. 2004, 5: 889-899. 10.1038/nrg1489.PubMedView ArticleGoogle Scholar
  19. Nei M, Gu X, Sitnikova T: Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci U S A. 1997, 94: 7799-7806. 10.1073/pnas.94.15.7799.PubMedPubMed CentralView ArticleGoogle Scholar
  20. Kruiswijk CP, Hermsen TT, Westphal AH, Savelkoul HF, Stet RJ: A novel functional class I lineage in zebrafish (Danio rerio), carp (Cyprinus carpio), and large barbus (Barbus intermedius) showing an unusual conservation of the peptide binding domains. J Immunol. 2002, 169: 1936-1947.PubMedView ArticleGoogle Scholar
  21. Graser R, Vincek V, Takami K, Klein J: Analysis of zebrafish Mhc using BAC clones. Immunogenetics. 1998, 47: 318-325. 10.1007/s002510050364.PubMedView ArticleGoogle Scholar
  22. Sultmann H, Mayer WE, Figueroa F, O'Huigin C, Klein J: Organization of Mhc class II B genes in the zebrafish (Brachydanio rerio). Genomics. 1994, 23: 1-14. 10.1006/geno.1994.1452.PubMedView ArticleGoogle Scholar
  23. Phillips RB, Zimmerman A, Noakes MA, Palti Y, Morasch MR, Eiben L, Ristow SS, Thorgaard GH, Hansen JD: Physical and genetic mapping of the rainbow trout major histocompatibility regions: evidence for duplication of the class I region. Immunogenetics. 2003, 55: 561-569. 10.1007/s00251-003-0615-4.PubMedView ArticleGoogle Scholar
  24. Sambrook JG, Russell R, Umrania Y, Edwards YJ, Campbell RD, Elgar G, Clark MS: Fugu orthologues of human major histocompatibility complex genes: a genome survey. Immunogenetics. 2002, 54: 367-380. 10.1007/s00251-002-0478-0.PubMedView ArticleGoogle Scholar
  25. Klein J, Sato A: Birth of the major histocompatibility complex. Scand J Immunol. 1998, 47: 199-209. 10.1046/j.1365-3083.1998.00292.x.PubMedView ArticleGoogle Scholar
  26. Seeger A, Mayer WE, Klein J: A complement factor B-like cDNA clone from the zebrafish (Brachydanio rerio). Mol Immunol. 1996, 33: 511-520. 10.1016/0161-5890(96)00002-8.PubMedView ArticleGoogle Scholar
  27. Nonaka M, Takahashi M, Sasaki M: Molecular cloning of a lamprey homologue of the mammalian MHC class III gene, complement factor B. J Immunol. 1994, 152: 2263-2269.PubMedGoogle Scholar
  28. Kuroda N, Wada H, Naruse K, Simada A, Shima A, Sasaki M, Nonaka M: Molecular cloning and linkage analysis of the Japanese medaka fish complement Bf/C2 gene. Immunogenetics. 1996, 44: 459-467. 10.1007/s002510050153.PubMedView ArticleGoogle Scholar
  29. Gongora R, Figueroa F, Klein J: Independent duplications of Bf and C3 complement genes in the zebrafish. Scand J Immunol. 1998, 48: 651-658. 10.1046/j.1365-3083.1998.00457.x.PubMedView ArticleGoogle Scholar
  30. Yamashita M, Hirayoshi K, Nagata K: Characterization of multiple members of the HSP70 family in platyfish culture cells: molecular evolution of stress protein HSP70 in vertebrates. Gene. 2004, 336: 207-218. 10.1016/j.gene.2004.04.023.PubMedView ArticleGoogle Scholar
  31. Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS: The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res. 2005, 15: 1307-1314. 10.1101/gr.4134305.PubMedPubMed CentralView ArticleGoogle Scholar
  32. Ensembl: [http://www.ensembl.org/]
  33. Flajnik MF, Kasahara M: Comparative genomics of the MHC: glimpses into the evolution of the adaptive immune system. Immunity. 2001, 15: 351-362. 10.1016/S1074-7613(01)00198-4.PubMedView ArticleGoogle Scholar
  34. Li L, Huang GM, Banta AB, Deng Y, Smith T, Dong P, Friedman C, Chen L, Trask BJ, Spies T, Rowen L, Hood L: Cloning, characterization, and the complete 56.8-kilobase DNA sequence of the human NOTCH4 gene. Genomics. 1998, 51: 45-58. 10.1006/geno.1998.5330.PubMedView ArticleGoogle Scholar
  35. Gongora R, Zaleska-Rutczynska Z, Takami K, Figueroa F, Klein J: Linkage of RXRB-like genes to class I and not to class II Mhc genes in the zebrafish. Immunogenetics. 1998, 48: 141-143. 10.1007/s002510050414.PubMedView ArticleGoogle Scholar
  36. Abi-Rached L, Gilles A, Shiina T, Pontarotti P, Inoko H: Evidence of en bloc duplication in vertebrate genomes. Nat Genet. 2002, 31: 100-105. 10.1038/ng855.PubMedView ArticleGoogle Scholar
  37. Samonte IE, Sato A, Mayer WE, Shintani S, Klein J: Linkage relationships of genes coding for alpha2-macroglobulin, C3 and C4 in the zebrafish: implications for the evolution of the complement and Mhc systems. Scand J Immunol. 2002, 56: 344-352. 10.1046/j.1365-3083.2002.01154.x.PubMedView ArticleGoogle Scholar
  38. Jentsch TJ, Stein V, Weinreich F, Zdebik AA: Molecular structure and physiological function of chloride channels. Physiol Rev. 2002, 82: 503-568.PubMedView ArticleGoogle Scholar
  39. Kulski JK, Shiina T, Anzai T, Kohara S, Inoko H: Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man. Immunol Rev. 2002, 190: 95-122. 10.1034/j.1600-065X.2002.19008.x.PubMedView ArticleGoogle Scholar
  40. Takami K, Zaleska-Rutczynska Z, Figueroa F, Klein J: Linkage of LMP, TAP, and RING3 with Mhc class I rather than class II genes in the zebrafish. J Immunol. 1997, 159: 6052-6060.PubMedGoogle Scholar
  41. Clark MS, Shaw L, Kelly A, Snell P, Elgar G: Characterization of the MHC class I region of the Japanese pufferfish (Fugu rubripes). Immunogenetics. 2001, 52: 174-185. 10.1007/s002510000285.PubMedView ArticleGoogle Scholar
  42. Matsuo MY, Asakawa S, Shimizu N, Kimura H, Nonaka M: Nucleotide sequence of the MHC class I genomic region of a teleost, the medaka (Oryzias latipes). Immunogenetics. 2002, 53: 930-940. 10.1007/s00251-001-0427-3.PubMedView ArticleGoogle Scholar
  43. The Danio rerio Sequencing Project: [http://www.sanger.ac.uk/Projects/D_rerio/]
  44. Zebrafish Genome Fingerprinting Project: [http://www.sanger.ac.uk/Projects/D_rerio/WebFPC/zebrafish/]
  45. Sprague J, Clements D, Conlin T, Edwards P, Frazer K, Schaper K, Segerdell E, Song P, Sprunger B, Westerfield M: The Zebrafish Information Network (ZFIN): the zebrafish model organism database. Nucleic Acids Res. 2003, 31: 241-243. 10.1093/nar/gkg027.PubMedPubMed CentralView ArticleGoogle Scholar
  46. VEGA: [http://vega.sanger.ac.uk/]
  47. Bairoch A, Apweiler R, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, Martin MJ, Natale DA, O'Donovan C, Redaschi N, Yeh LS: The Universal Protein Resource (UniProt). Nucleic Acids Res. 2005, 33 (): D154-9. 10.1093/nar/gki070.PubMedPubMed CentralView ArticleGoogle Scholar
  48. Higgins DG, Thompson JD, Gibson TJ: Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996, 266: 383-402.PubMedView ArticleGoogle Scholar
  49. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
  50. Galtier N, Gouy M, Gautier C: SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996, 12: 543-548.PubMedGoogle Scholar

Copyright

© Sambrook et al; licensee BioMed Central Ltd. 2005

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.