Comprehensive analysis of MHC class II genes in teleost fish genomes reveals dispensability of the peptide-loading DM system in a large part of vertebrates
- Johannes M Dijkstra†1,
- Unni Grimholt†2Email author,
- Jong Leong3,
- Ben F Koop3 and
- Keiichiro Hashimoto1Email author
© Dijkstra et al.; licensee BioMed Central Ltd. 2013
Received: 21 August 2013
Accepted: 4 November 2013
Published: 26 November 2013
Classical major histocompatibility complex (MHC) class II molecules play an essential role in presenting peptide antigens to CD4+ T lymphocytes in the acquired immune system. The non-classical class II DM molecule, HLA-DM in the case of humans, possesses critical function in assisting the classical MHC class II molecules for proper peptide loading and is highly conserved in tetrapod species. Although the absence of DM-like genes in teleost fish has been speculated based on the results of homology searches, it has not been definitively clear whether the DM system is truly specific for tetrapods or not. To obtain a clear answer, we comprehensively searched class II genes in representative teleost fish genomes and analyzed those genes regarding the critical functional features required for the DM system.
We discovered a novel ancient class II group (DE) in teleost fish and classified teleost fish class II genes into three major groups (DA, DB and DE). Based on several criteria, we investigated the classical/non-classical nature of various class II genes and showed that only one of three groups (DA) exhibits classical-type characteristics. Analyses of predicted class II molecules revealed that the critical tryptophan residue required for a classical class II molecule in the DM system could be found only in some non-classical but not in classical-type class II molecules of teleost fish.
Teleost fish, a major group of vertebrates, do not possess the DM system for the classical class II peptide-loading and this sophisticated system has specially evolved in the tetrapod lineage.
KeywordsMHC Class II Classical Non-classical DM Peptide-loading Teleost fish Genomics Evolution
The highly polymorphic classical MHC class II molecules can present exogenous antigenic peptides including those derived from proteins of many pathogens to CD4+ T lymphocytes in the acquired immune system . CD4+ T lymphocytes then can exert various functions such as helper activities toward other immune cells, e.g. B lymphocytes, macrophages and dendritic cells for their activation .
In the mammalian acquired immune system, a non-classical MHC class II molecule, HLA-DM in humans (H2-DM in mice), plays an important role in proper peptide presentation by the classical MHC class II molecules [1, 2]. A newly synthesized classical MHC class II molecule, which is a heterodimer composed of the α and β chains, is transported from the endoplasmic reticulum to endosomal compartments including the late endosomal MIIC (MHC class II compartment) by binding to a protein called the invariant chain. The invariant chain blocks the peptide-binding groove of a classical MHC class II molecule by its “CLIP (class II-associated invariant chain peptide)” region so that other endogenous peptides cannot bind to the groove in the endoplasmic reticulum [1, 2]. After digestion of the invariant chain by endosomal proteases, CLIP is dissociated from the groove by the non-classical class II DM molecule in the MIIC, and then other peptides, including those derived from exogenous antigens, can bind to the groove of the classical MHC class II molecule [1–6]. The DM molecule also can induce the dissociation of relatively weakly bound peptides thus showing peptide-editing function. The DM molecule itself neither exhibits classical-type polymorphism nor shows binding capacity for peptide ligands .
DM molecules are critical for the classical MHC class II function, as exemplified by the observation that human mutant cell lines deficient in DM molecules, and also antigen-presenting cells from the DM-knockout mice, exhibited failure in proper MHC class II peptide presentation [7–10]. Similar to a classical MHC class II molecule, a DM molecule is a heterodimer composed of α and β chains, which are encoded by DMA and DMB genes, respectively, and possesses an overall structure similar to that of a classical MHC class II molecule, but with a unique narrow groove [11, 12]. Orthologous DM genes, DMA and DMB, have been identified not only in many mammals but also in chickens  and frogs , indicating phylogenetic conservation throughout tetrapods. In all investigated tetrapods, the DM genes reside in the Mhc region along with the classical MHC class I and class II genes .
DM-lineage genes, however, have not been reported from the largest group of vertebrates, teleost fish, which include more than 26000 species, about 40% of all the species of vertebrates . Teleost fish appear to possess effective acquired immune functions including presumable T lymphocyte-dependent responses to exogenous antigens . Like in tetrapods, various important genes of the MHC class II system have been found in their genomes, which include genes for MHC class II [17, 18], CD4 [19, 20] and invariant chain [21–23]. The MHC class II loci in teleost fish display some unique features, namely their non-linkage with MHC class I loci and the lack of synteny of the class II loci between several teleost fish species [24–27]. Regarding non-classical class II genes in teleost fish, studies have been rather limited so far e.g. [28–30].
Previously, the absence of DM-like genes was speculated based on extensive homology searches using teleost fish genomes e.g. . However, based on the negative results, one could not obtain a definitive answer, e.g., primitive DM-like molecules might be highly divergent from the known DM sequences. Very recently, the structure of a complex between a classical class II molecule and a DM molecule has been determined and the study revealed the interacting amino acid residues between the two molecules . This allows us, for the first time, to analyze teleost fish class II molecules regarding the possession of the functional residues critically important for “the DM system” (in this paper defined as the peptide-loading system with the DM molecule-involved special mechanism). Taking advantage of recent progress of genome research projects, we searched and analyzed MHC class II genes in teleost fish comprehensively. We typified teleost fish class II genes and their protein products based on various classical/non-classical characteristics and then examined the possession of functional residues critical for the DM system. Based on our results, we could draw a clear conclusion about the lack of the DM system in teleost fish. We also discuss possible functions of the intriguing non-classical class II molecules revealed in teleost fish including ones newly identified in the present study.
Results and discussion
Comprehensive search for teleost fish MHCclass II genes
Novel expressed MHC class II genes in an Mhcregion of the Atlantic salmon genome and their ancient characteristics
In a duplicated Mhc region on chromosome 2 (linkage group 1 (LG1); Additional file 8: Text S2) of the genome of Atlantic salmon (Salmo salar), our model fish, we could find two novel class II genes, named DEA and DEB, the former encoding an α chain and the latter a β chain of an MHC class II molecule (Figure 2B and Additional file 2: Figure S2). These genes are located closely to each other like a typical pair of MHC class II genes of mammals. A teleost fish specific whole genome duplication (Figure 1 TGD) resulted in duplication of the Mhc region to two different chromosomes, each experiencing subsequent rearrangements and gene losses [32, 33]. Teleost classical MHC class I genes were found in only one of those two Mhc regions e.g. . Our mapping of the Atlantic salmon DEA/DEB to chromosome 2 is the first identification of intact MHC class II loci within a teleost Mhc region (Figure 2B and Additional file 2: Figure S2), representing the TGD duplicate Mhc region without class I [32, 33]. Early in the salmonid line an additional whole genome duplication occurred (Figure 1 SGD), which presumably resulted in duplicate Mhc regions on chromosomes 2 and 5 (Figure 2B and Additional file 2: Figure S2). In the Atlantic salmon Mhc region of chromosome 5, DE group genes also can be found, namely the pseudogenes DFAψ and DFBψ (Figure 2B and Additional file 2: Figure S2). Both DEA/DEB on chromosome 2 and DFAψ/DFBψ on chromosome 5 are closely linked with typical Mhc scaffold genes such as BRD2, TAP1, PBX2 and RGL2 genes (Additional file 8: Text S2 and Additional file 2: Figure S2). A possible orthologue of DFAψ, equally partial and inactivated, has previously been described in the salmonid fish rainbow trout .
DE-related MHC class II sequences could be identified in the databases for a few other teleost fish including red seabream (a neoteleost fish), fathead minnow and dojo loach (both Cypriniformes), and importantly several primitive non-teleost ray-finned fish like spotted gar, sturgeon and paddlefish (Figures 3 and 4, Additional file 3: Figure S3, Additional file 10: Text S4). In the genome of the spotted gar (Lepisosteus oculatus; PreEnsembl Lepocu1) we could identify four DE-related sequences (Figure 2A).
In Figure 3, the predicted amino acid sequences of the membrane-distal domains of the Atlantic salmon DE molecule are compared with those of other MHC class II molecules. The DE and related sequences significantly differ from all the teleost fish class II sequences previously published. Particularly, they share several residues with class II molecules of the other jawed vertebrates (not only the classical class II molecules of tetrapods and cartilaginous fish but also the DM molecules of tetrapods) which are not observed in previously reported teleost fish class II. Those include F12, Q14, P18, R44 of class II α1 domain and D41 and R93 of class II β1 domain (shaded with grey in Figure 3). Regarding specific residues of DE and related sequences, those consist of the following (shaded with lime green in Figure 3): R8, L10, N16, W51 of α1 domain and E16 of β1 domain.
The Mhc regions in tetrapods and cartilaginous fish comprise a similar set of MHC class I and II genes, as well as other Mhc scaffold genes [35, 36]. In tetrapods, the DM genes also reside in this region. The teleost fish genomes, in contrast, represent a derived situation in that the classical MHC class I and II genes are not linked [24, 26]. In the primitive ray-finned fish spotted gar genome, we found that a classical-type class I gene (for its molecular features, see Additional file 2: Figure S2) is linked with class II genes (Figure 2A), which suggests that the non-linkage between classical class I and II was established only within teleost fish.
Clarification of the three major MHC class II groups in teleost fish
With the addition of a newly identified group, teleost fish class II genes can be organized into three major groups, namely, DA, DB (named after a zebrafish DB sequence) and DE group, based on several sequence features including specific insertions and/or deletions (Figure 3 and Additional file 10: Text S4). Previously, reported DA and DB group sequences were classified as classical and non-classical class II, respectively, assignments which were based only on polymorphism and expression analyses and not on comparison of characteristic amino acid residues (Additional file 9: Text S3). The DE and DA groups are well supported as distinct lineages by phylogenetic tree analyses. And the tree shows that early in ray-finned fish evolution there was a separation between DE lineage and DA/DB lineage, and that from the latter the DA lineage sprouted in teleost fish (Figure 4 and Additional file 3: Figure S3). How the DA lineage relates to extant lineages within the teleost DB group is not clear at present. Three discernible lineages within the DB group are represented by (i) genes found in zebrafish, (ii) the Atlantic salmon DCA/DCB plus genes found in the neoteleost genomic synteny region “S1”, and (iii) the Atlantic salmon DBA/DBB plus genes found in the neoteleost genomic synteny region “S2” (Figures 2, 3,4 and Additional file 4: Table S1).
Previous studies on teleost fish DA and DB group genes reported a lack of synteny between teleost fish class II genomic regions, suggesting the occurrence of multiple translocations and locus turnovers e.g. [25, 27]. Our study largely agrees with and extends those previous observations in these respects (Figure 2 and Additional file 4: Table S1), although within neoteleosts, some syntenic MHC class II regions could be found (Figure 2: Additional file 1: Figure S1, Additional file 2: Figure S2 and Additional file 3: Figure S3). The observation that the classical DA lineage shows little locus conservation is unexpected since the longevity of a locus should be helpful for the maintenance of allelic variations. For example, a high rate of loci turnover was claimed to be responsible for the limited diversity of the MHC class I genes in the cotton-top tamarin  and, further, polymorphic MHC class I genes tend to map to more ancient genomic regions than non-polymorphic ones throughout the jawed vertebrates e.g. .
Classical and non-classical features of teleost fish class II genes
Comparison of various features of class II genes and their protein products in jawed vertebrates
Cartilaginous fish classical
Teleost fish DA
Teleost fish DB
Teleost fish DE
Located at Mhc
Linked with classical-type class I
Peptide-binding residues c
CD4-binding βS144, βE162
+ /― e
Endosomal sorting motif f
+ /― e
+ /― e
Amino acid at α125
Expression and polymorphism
The DA group includes all the highly expressed teleost class II genes (Figure 2, Additional file 5: Table S2) and previously reported polymorphic ones (Additional file 9: Text S3). There have been relatively few reports on teleost non-classical class II genes e.g. [28–30] (Additional file 9: Text S3), and our present study added identification, gene-specific expression and polymorphism analyses of the Atlantic salmon DCB (DB group), DEA and DEB (DE group) (Additional file 8: Text S2). Together with our previous analyses of Atlantic salmon DB group genes, DBA, DBB, DCA and DDA, and non-specific transcriptome analysis in the present study (Additional file 8: Text S2), the results show that all these genes are essentially non-polymorphic and expressed at much lower levels than the classical DAA and DAB genes of the DA group in the same species. Especially the expression of Atlantic salmon DEA and DEB is very low. The Atlantic salmon genes DCA, DCB and DEA have rather tissue-specific expression patterns unknown for mammalian class II genes (Additional file 8: Text S2).
Peptide binding capacity
Only DA molecules, but not DB and DE, display a high degree of conservation of the α1 domain residues αN62 and αN69 and the β1 domain residues βH81 and βN82, which in mammalian classical molecules make important hydrogen bonds with the backbone of peptide ligands [38, 39] (Figure 3, Table 1, Additional file 10: Text S4 and Additional file 6: Table S3). Best conserved is βN82, known to be of particular importance  (Additional file 6: Table S3). At the other three positions, teleost fish DA molecules exhibit some variations comparable to some species in mammals (Additional file 6: Table S3). Although the presence of such peptide-binding residues in DA molecules has been reported before , the detailed analyses of teleost fish DB molecules have not been reported. Among the DA group, a few DA molecules lack these peptide-binding residues and might exert nonclassical functions akin to DO in mammals, which diverged from classical molecules in relatively recent times.
CD4 binding capacity
Recently, several amino acid residues important for the interaction between CD4 and class II molecules were revealed . Two of those residues at the interface of the two molecules appear highly conserved throughout tetrapod species, namely βS144 and βE162 (Additional file 10: Text S4). We found that these two residues are also highly conserved in teleost fish DA, whereas cartilaginous fish classical molecules, tetrapod DM molecules and some teleost fish DB molecules lack these residues (Table 1 and Additional file 10: Text S4).
Endosomal sorting motif
A characteristic feature of the non-classical DM molecules is the possession of an endosomal sorting motif in the β chain cytoplasmic tail [43, 44] (Table 1). We found potential tyrosine-based endosomal sorting motifs in the β chain of a few teleost DB group molecules although their location differs from those in the DM molecules (Table 1 and Additional file 10: Text S4E).
Preservation of classical genes throughout teleost fish
Based on the current database information, all the teleost fish species that we investigated have DA, most of them have DB, and a few of them have DE group genes. Conservation of gene copies preserving classical features accompanied by seemingly random loss of older non-classical gene duplicates is somewhat reminiscent of MHC class I evolution described for higher vertebrates . However, the mode of the teleost fish MHC class II evolution highly contrasts with that of the tetrapod species in which not only the classical class II but also the non-classical DM genes are highly conserved. Recently, it was reported that Atlantic cod, a teleost fish, does not possess various genes in the MHC class II system such as those for MHC class II molecules, CD4 and invariant chain . Although this apparent loss of the MHC class II system should have various disadvantages and actually the cods are known to have a poor adaptive antibody response, they do survive and thrive. The cod situation reflects the plasticity of the teleost fish immune system in which other factors may adapt to a large variation in the MHC class II system .
Phylogeny of the DM system through the window of the critical functional residues
A critical residue of the DM system is not found in the teleost fish classical class II molecules
Recently, the structure of HLA-DM/HLA-DR1 complex was clarified . The overall structure of the complex is largely consistent with the previous independent estimation of the interface of the two molecules based on experiments using mutagenesis [46–49] and tethered complexes . In the side-by-side structure, the interface is mainly formed by the α chains of the two molecules, and a lateral surface of the DRα1 domain, close to the N-terminus of the peptide-binding groove, interacts with DM α chain and additionally DMβ1 . The structural study revealed two key amino acid residues (αN125 of DM and αW43 of DR1) in the interaction between the DM molecule and the classical class II DR1 molecule at pH 5.5, which is within the range of the physiological late endosomal pH suitable for the DM activity . In the structure of the DM/DR1 complex, a tryptophan residue of the DR1 molecule (αW43) flips from the original location and its indole ring nitrogen atom interacts with an asparagine of the DM molecule (αN125) through a hydrogen bond . This was elucidated by comparison between the structures of the DM-bound  and -unbound DR molecule . The change of αW43 position is accompanied by conformational alterations in the P1 pocket peptide-binding region of the DR1 molecule, which include the novel formation of a long α-helical segment with a short break and the repositioning of the hydrophobic αF51 into the P1 pocket . These changes explain the dissociation of CLIP, the stabilization of empty class II molecules, and further the selection of high affinity peptide ligands . The previous study indicated that αW43F mutation of the DR molecule greatly reduced both the DM function and the binding to DM molecule . The effects of αN125A and αN125R mutations in the DM molecule were also examined and these mutations caused a loss of both the DM activity and the binding to the DR molecule . HLA-DO, a human non-classical class II molecule, which can bind tightly to the DM molecule and is known to be an inhibitory modulator of the DM molecule, also possesses αW43. The structure of HLA-DM/HLA-DO complex was independently reported very recently and it also revealed the important participation of αW43 of HLA-DO and αN125 of HLA-DM in the complex formation in which HLA-DO behaves as a mimic of the classical class II molecule . Very importantly, the classical class II α chains of the DA group of teleost fish do not possess the αW43 critical for the interaction with the DM molecule (Table 1, Figure 4, Additional file 10: Text S4A).
As αW43 constitutes a part of the β-strand 4 and there are a few highly conserved amino acid residues near this position, we did not have any difficulties with the alignment of the position 43 between teleost fish and tetrapod class II α chain sequences. At the amino acid position corresponding to αW43, the teleost fish classical-type class II molecules, namely, those belonging to the DA group, exhibit variability (Additional file 10: Text S4A), but no tryptophan residue was observed. Retrospectively, the absence of a tryptophan and also some variability at the position 43 of the class II α chain of teleost fish can be recognized in a previous study using a few teleost fish sequences  and also in the other studies when the sequence alignments are adjusted e.g. [29, 53]. However, the meaning of these observations could not be understood in relation to the DM system before. In the present study, we comprehensively investigated various teleost fish genomes and examined many DA group members. We did not find any sub-lineages in the DA group in which class II α chains specifically possess αW43 (Additional file 10: Text S4A). The variable nature of the position 43 of the class II α chain without specific conserved residues suggests that the teleost fish classical class II molecules do not use this amino acid position for the interaction with some regulatory molecules. Rather, teleost fish appear to use this position to further increase the variation of the pocket.
Teleost fish αW43-containing class II molecules are classified as non-classical
In the DB group of teleost fish, we could find that six sequences (M5A of medaka, O97A, O9A, O57A and O33A of tilapia and DCA of Atlantic salmon) possess αW43 (Additional file 10: Text S4A). Although our previous study reported that the Atlantic salmon DCA gene for the α chain was not polymorphic based on the EST information , the identification and analysis of DCB gene for the β chain in the present study was important, as we know examples in which the extent of polymorphism is highly different between the genes coding for the α and β chains. In the case of human classical class II DR genes, DRA shows only limited polymorphism while DRB shows very high polymorphism . In the present study, the low polymorphism of DCB gene was clarified, and therefore we could classify the plausible pair of the Atlantic salmon DCA/DCB molecule as non-classical with the other supporting observations described in the following. Consistent with the non-classical feature of lacking polymorphism, DCA lacks both peptide-binding asparagine residues (αN62 and αN69) although DCB retains βN82. For the position 69, DCA possesses a hydrophobic residue like the DM sequences, a feature shared with most of the other teleost fish DB sequences (Additional file 6: Table S3). Further, both DCA and DCB genes showed unique expression patterns. They are expressed predominantly in the digestive tract, as observed both in the RT-PCR study and in the transcriptome analysis (Additional file 8: Text S2). These expression patterns are quite different from those of the classical class II genes of both teleost fish DA group and tetrapods. The classical MHC class II genes of human and Atlantic salmon show their highest expression in various immunologically important tissues when investigated with transcriptome analyses, and also human DM genes, HLA-DMA and HLA-DMB, show similar expression patterns (Additional file 8: Text S2). Another characteristic of DCB is that, as already mentioned briefly, DCB interestingly possesses a putative tyrosine-based endosomal sorting motif in its cytoplasmic tail like the DM molecules, although the position of this motif is different from those observed in the DM molecules (Additional file 10: Text S4E). This endosomal sorting motif can also be observed in molecules relatively closely related to DCB (Additional file 10: Text S4E). Non-classical identity of DCB could not be deduced from the analysis of the CD4-binding residues as DCB possesses both βS144 and βE162 (Additional file 10: Text S4D), but the retention of these residues in non-classical molecules can also be observed in other cases such as HLA-DO (Additional file 10: Text S4D).
Similar to DCA, the other five teleost fish class II α chain sequences that possess αW43 could be classified as non-classical based on the analysis of the residues important for peptide-binding although we do not have clear information about their polymorphism and expression patterns (Additional file 10: Text S4). Their apparently intact, presumable β chain partners (three out of five) could also be classified as non-classical as they lack the important peptide-binding residues (both βH81 and βN82 in two cases and βH81 in one case) and also lack at least one of the CD4-binding residues (Additional file 10: Text S4). All the identified teleost fish class II α chains that contain αW43 belong to the DB group of the synteny region “S1” or “S2” of neoteleost fish or their closely related molecules including the Atlantic salmon DCA and together they seem to have descended from a common ancestral molecule in the early phase of the evolution of the DB group (Figure 4).
In addition to αW43, the residues αK38 and αE40 of the DR molecule also participate in the formation of the hydrogen bonding network between DR and DM molecules  and previous mutagenesis experiments supported the importance of αE40 . Among the six teleost fish molecules which possess αW43, half of them also possess both αK38 and αE40. Although the chicken and frog classical class II molecules possess different amino acids, the teleost fish DB group molecules of the synteny region “S1” or “S2” group and their closely related ones and also the classical-type DA group molecules possess these residues at the relevant positions in relatively high frequency (Additional file 10: Text S4). The coelacanth class II molecule also possesses αE40 . In the hydrogen bonding network, αE40 of the DR molecule interacts with αR98 of the DM molecule. The conservation of αK98 can be observed in chicken  and frog DM molecules (Additional file 10: Text S4), and also often in the teleost fish DB group molecules and in some DA molecules (Additional file 10: Text S4). Therefore, various preconditions for the critical hydrogen bonding network between DM and DR molecules appear to have been already established in the common ancestor between teleost fish and tetrapods.
A teleost fish αN125-containing class II molecule is not a DM-equivalent
αN125 is conserved in mammalian and chicken DM molecules as previously noted  and it is also conserved in frog DM molecules (Additional file 10: Text S4). αN125 is not specific to the DM molecule, but is also observed in many classical class II α chains of the tetrapod lineage including frog, coelacanth and human DQA-related mammalian molecules, as well as in some closely related non-classical ones like human DOA (Additional file 10: TextS4). We could find, thus far, a single αN125-containing class II sequence in teleost fish, namely, DDA of Atlantic salmon (Additional file 10: Text S4). In our previous paper, we reported that DDA has little polymorphism based on available EST information. DDA does not possess two conserved peptide-binding asparagines (αN62 and αN69) and it possesses a hydrophobic residue at the position of 69 as found in most of the other teleost fish DB and also tetrapod DM molecules. DDA is predominantly expressed in spleen although in much lower amounts compared to the classical DAA (Additional file 8: Text S2). Based on various observations described above, DDA could be classified as a non-classical class II molecule, although we do not know about DDB at present. When we conducted homology searches with Atlantic salmon DDA, we could not retrieve DDA-like sequences from genome sequence databases of other teleost fish species. Therefore, although DDA is a non-classical class II α chain possessing αN125, it is not like DM that is highly conserved throughout tetrapod species.
Novel teleost fish class II group DE
In the present study, we identified a new teleost fish class II group called DE. As described above, DEA and DEB genes of Atlantic salmon do not show classical-type polymorphism and their expression levels are very low (Additional file 8: Text S2). Although the predicted Atlantic salmon DE molecule possesses both βS144 and βE162 for CD4-binding, it lacks αN62, βH81 and βN82 for peptide-binding. Based on these observations, we could exclude the Atlantic salmon DE molecule from the classical class II group. As also described in the previous section, the Atlantic salmon DE molecule shares several amino acid residues with the cartilaginous fish and tetrapod classical, and also tetrapod DM molecules (Figure 3, Additional file 10: Text S4). From the standpoint of the conservation of these residues, the DE molecule is closest to the DM molecules among all the known teleost fish class II molecules. All the available DEA sequences found in teleost fish and also in the primitive ray-finned fish possesses an aspartic acid residue at the position of 125 instead of an asparagine (Additional file 10: Text S4C). Intriguingly, an aspartic acid also can participate in a hydrogen-bond interaction with an indole ring nitrogen atom of a tryptophan like an asparagine e.g., .
The observation that the gar class II sequence of the conserved DE lineage possesses αD125 and additionally another gar class II sequence (Additional file 10: Text S4C) possesses αN125 suggests that αN125/D125 already appeared in the class II molecules of an ancestor of Osteichthyes, whereas available cartilaginous fish sequences possess a glycine at this position that appears ancestral to the MHC family  (Table 1, Additional file 10: Text S4C). It should be noted that a glycine cannot provide hydrogen-bond capacity necessary for the DM function. Together with the other residues already discussed above, αN125/D125 in ray-finned fish molecules further suggests the establishment of the preconditions for the critical hydrogen bonding network between DM and DR molecules in the common ancestor of teleost fish and tetrapods.
With the identification of DE genes in primitive ray-finned fish like the spotted gar, the coexistence of classical and non-classical class II lineages for a long evolutionary time was demonstrated in the present study (Figure 1). DE genes are found in some teleost fish, while they could not be identified in the genome databases and the EST databases of the other teleost fish. Thus, although DE molecules share some non-classical features with DM molecules, they seem to have quite different characteristics regarding stable inheritance.
Teleost fish do not possess the DM system
Without a complete coverage of teleost fish genomes, it is logically not possible to deny the existence of DM-equivalent genes. With this limitation, we clarified that some non-classical class II genes exist in the teleost fish genomes whose protein products partially share some characteristics with the DM molecule. Actually there exist some teleost fish non-classical class II molecules that possess αN125 or αD125. Mainly based on their conservation profiles in the teleost fish genomes, they appear not to be DM-equivalents of teleost fish. In the middle of the β2 domain, the tetrapod DM molecules have a unique insertion of several amino acid residues and near or in this region there are a few residues (e.g. βE47 and βL51) influential for the DM/DR interaction supported by the mutation studies  and also by the structural study . So far, teleost fish genes whose protein products possess these features have not been identified.
To obtain a clear conclusion about a possible DM system in teleost fish, it was necessary to investigate the other side of the DM mechanism, namely, the possession of the critical tryptophan residue in the classical class II molecules. Based on the analyses regarding αW43 in the teleost fish classical class II molecules, we could conclude that teleost fish classical class II molecules do not possess αW43 and therefore do not possess the DM system in which the interaction between αW43 of the classical class II molecule and αN125 of the DM molecule is critically important. Although we could observe some teleost fish class II molecules which possess αW43, we could clearly classify them as non-classical. In contrast, coelacanth classical-type class II molecule has this tryptophan , and this might be consistent with the presence of the DM system in the primitive stage of the tetrapod lineage. Although available data are very limited, known classical polymorphic class II α chains of cartilaginous fish do not possess αW43 (Figure 3, Additional file 10: Text S4A), e.g. the nurse shark molecule has an alanine at this position which does not have hydrogen-bonding capacity. Therefore, αW43 appears fixed in the classical class II α chain only from the level of coelacanth in the lobe-finned fish/tetrapod line  (Figure 3, Additional file 10: Text S4A). As the coelacanth class II molecule also possesses αN125, both the important αW43 and αN125 prerequisites for the DM function may have been fixed at the early phase of the tetrapod lineage (Figure 1). All these observations support that the DM system has specifically evolved in the tetrapod lineage.
Possible MHC class II peptide-presentation system of teleost fish without DM
Studies on the teleost fish MHC class II peptide-presentation system have thus far been very limited. If we assume that the teleost fish invariant chain [21–23] possesses functions similar to those of mammals [1, 2], some basic issues need to be considered for the peptide-loading pathway without the DM system. Those include the dissociation of CLIP-equivalent fragments from classical class II molecules and the stabilization of empty classical class II molecules during peptide-exchange reactions. First, as some mammalian classical MHC class II molecules bind CLIP with low affinity, the rapid dissociation of CLIP has been observed at an endosomal low pH [2, 57]. Therefore, if the binding of teleost fish CLIP-equivalent fragments to classical class II molecules is not so strong, the fragments may dissociate without help of a DM-equivalent. Second, classical class II gene duplications can produce an evolutionary reservoir of nonclassical class II genes. Some classical MHC class II molecules possess intrinsic affinity for each other , and similar interactive forces may have been the evolutionary basis of the tetrapod DM molecules for establishing the specific interaction wtih classical class II molecules. Some teleost fish non-classical class II molecules may also have acquired stabilizing activity toward classical class II molecules in evolution. Thus, although mechanisms cannot be identical, there might be some overlap between the tetrapod DM and the teleost non-classical class II functions in aiding classical class II molecules.
For distinct subgroups of teleost fish, we may find different strategies for the MHC class II peptide-loading. With four kinds of non-classical class II molecules, Atlantic salmon might have a peptide-loading system uniquely evolved with these molecules. Other than peptide presentation, some reports indicated that tetrapod MHC class II molecules can be involved in signaling pathways , and some teleost fish class II molecules may participate in similar or yet unknown functions.
The accelerated progress on whole genome sequence and also expressed sequence information of various species certainly is valuable to gain an evolutionary bird’s-eye view of important biological systems. The observation that all the authentic polymorphic classical class II-type molecules of teleost fish do not possess the critical residue αW43 led us to conclude that teleost fish do not possess the DM system. The DM molecule appears to have acquired highly sophisticated and efficient mechanisms for peptide editing and stabilization of the classical class II molecules, ensuring its preservation throughout tetrapod evolution. As teleost fish comprise a significant part of the jawed vertebrates (more than 40% of all the species), the present study revealed that both DM-dependent and -independent systems are present as major fractions in the jawed vertebrates. Our study also suggests that preconditions necessary for the important hydrogen-bonding network in the DM system appeared in the common ancestor of teleost fish and tetrapods. Exploring the teleost fish class II peptide-loading system would constitute an important part for the comprehensive understanding of the MHC class II antigen-presentation systems in the jawed vertebrates. Future studies on the non-classical class II molecules of teleost fish should reveal whether they have functions to support the classical class II molecules like the DM molecules or other and yet unknown functions.
Data mining and bioinformatics
A mixture of annotated and un-annotated MHC class II α and β sequences were identified using Ensembl’s Biomart and the GO term for class II (GO: 0042613) supplemented with various blastN and TblastN searches of Ensembl and NCBI databases using evolutionary diverged as well as species-specific sequences. For Atlantic salmon, we supplemented the six known Atlantic salmon MHC class II genes  with blastN and TblastN searches using available salmon genome sequences including and mostly from our ongoing genome sequencing project (in part available at either cGRASP (http://web.uvic.ca/grasp/) or NCBI (http://www.ncbi.nlm.nih.gov/)). Open reading frames were predicted using GenScan  and Fgenesh  and by comparison with known MHC sequences. Some small pseudogene remnants which did not contribute to evolutionary understanding were neglected.
Expression pattern of Atlantic salmon MHC class II gene transcripts
Expression of Atlantic salmon class II genes was estimated by transcriptome analysis and by gene-specific RT-PCR analysis. The transcriptome data, using tissues of a one-year old salmon, comprised >70.000 non-redundant contigs and >50 million reads per tissue (Additional file 8: Text S2 (Table TS2-2A)), and data on MHC class II expression agree fully with our previous RT-PCR analysis results on DAA, DAB, DBA, DBB, DCA and DDA and our present RT-PCR analysis results on DCB, DEA and DEB for three adult salmon individuals (Additional file 8: Text S2-2B). For details of methods and results see Additional file 8: Text S2-2A and Additional file 8: Text S2-2B. The protocol for the ethics and use of the animals was in accordance with the Animal Care at the University of Victoria, and all animal experiments comply with the current laws of Norway.
The alignments of the MHC class II sequences shown in Additional file 10: Text S4 were done manually, based on structural and evolutionary considerations. The evolutionary history of these manually aligned sequences was inferred using the Neighbor-Joining method  using MEGA5 software . See Additional file 3: Figure S3 for details.
Class II-associated invariant chain peptide
Major histocompatibility complex
MHC class II compartment
Million years ago
Teleost-specific whole genome duplication
Salmonid-specific whole genome duplication.
This work was supported by JSPS KAKENHI Grant Number 22580213 and 25450310 (to J. M. D.), a Norwegian Research Council grant (Havbruk # NFR 185217/S40) (to U. G.), an NSERC grant (to B. F. K.) and the MEXT-Supported Program for the Strategic Research Foundation at Private Universities, 2008–2012.
- Neefjes J, Jongsma MLM, Paul P, Bakke O: Towards a systems understanding of MHC class I and MHC class II antigen presentation. Nature Rev Immunol. 2011, 11: 823-836.
- Schulze MS, Wucherpfennig KW: The mechanism of HLA-DM induced peptide exchange in the MHC class II antigen presentation pathway. Curr Opin Immunol. 2012, 24: 105-111. 10.1016/j.coi.2011.11.004.PubMed CentralPubMedView Article
- Sloan VS, Cameron P, Porter G, Gammon M, Amaya M, Mellins E, Zaller DM: Mediation by HLA-DM of dissociation of peptides from HLA-DR. Nature. 1995, 375: 802-806. 10.1038/375802a0.PubMedView Article
- Denzin LK, Cresswell P: HLA-DM induces CLIP dissociation from MHC classII αβ dimers and facilitates peptide loading. Cell. 1995, 82: 155-165. 10.1016/0092-8674(95)90061-6.PubMedView Article
- Sherman MA, Weber DA, Jensen PE: DM enhances peptide binding to class II MHC by release of invariant chain-derived peptide. Immunity. 1995, 3: 197-205. 10.1016/1074-7613(95)90089-6.PubMedView Article
- Weber DA, Evavold BD, Jensen PE: Enhanced dissociation of HLA-DR-bound peptides in the presence of HLA-DM. Science. 1996, 274: 618-620. 10.1126/science.274.5287.618.PubMedView Article
- Morris P, Shaman J, Attaya M, Amaya M, Goodman S, Bergman C, Monaco JJ, Mellins E: An essential role for HLA-DM in antigen presentation by class II major histocompatibility molecules. Nature. 1994, 368: 551-554. 10.1038/368551a0.PubMedView Article
- Fling SP, Arp B, Pious D: HLA-DMA and –DMB genes are both required for MHC class II/peptide complex formation in antigen-presenting cells. Nature. 1994, 368: 554-558. 10.1038/368554a0.PubMedView Article
- Miyazaki T, Wolf P, Tourne S, Waltzinger C, Dierich A, Barois N, Ploegh H, Benoist C, Mathis D: Mice lacking H2-M complexes, enigmatic elements of the MHC class II peptide-loading pathway. Cell. 1996, 84: 531-541. 10.1016/S0092-8674(00)81029-6.PubMedView Article
- Fung-Leung W-P, Surh CD, Liljedahl M, Pang J, Leturcq D, Peterson PA, Webb SR, Karlsson L: Antigen presentation and T cell development in H2-M-deficient Mice. Science. 1996, 271: 1278-1281. 10.1126/science.271.5253.1278.PubMedView Article
- Mosyak L, Zaller DM, Wiley DC: The structure of HLA-DM, the peptide exchange catalyst that loads antigen onto class II MHC molecules during antigen presentation. Immunity. 1998, 9: 377-383. 10.1016/S1074-7613(00)80620-2.PubMedView Article
- Fremont DH, Crawford F, Marrack P, Hendrickson WA, Kappler J: Crystal structure of mouse H2-M. Immunity. 1998, 9: 385-393. 10.1016/S1074-7613(00)80621-4.PubMedView Article
- Kaufman J, Milne S, Göbel 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 Article
- Ohta Y, Goetz W, Hossain MZ, Nonaka M, Flajnik MF: Ancestral organization of the MHC revealed in the amphibian Xenopus. J Immunol. 2006, 176: 3674-3685.PubMedView Article
- Nelson JS: Fishes of the world. 2006, New York: Wiley, 4
- Iwama G, Nakanishi T: The fish immune system. 1996, San Diego: Academic Press
- Hashimoto K, Nakanishi T, Kurosawa Y: Isolation of carp genes encoding major histocompatibility complex antigens. Proc Natl Acad Sci USA. 1990, 87: 6863-6867. 10.1073/pnas.87.17.6863.PubMed CentralPubMedView Article
- Sultmann H, Mayer WE, Figueroa F, O’hUigin C, Klein J: Zebrafish Mhc class II alpha chain-encoding genes: polymorphism, expression, and function. Immunogenetics. 1993, 38: 408-420.PubMedView Article
- Suetake H, Araki K, Suzuki Y: Cloning, expression, and characterization of fugu CD4, the first ectothermic animal CD4. Immunogenetics. 2004, 56: 368-374.PubMedView Article
- Dijkstra JM, Somamoto T, Moore L, Hordvik I, Ototake M, Fischer U: Identification and characterization of a second CD4-like gene in teleost fish. Mol Immunol. 2006, 43: 410-419. 10.1016/j.molimm.2005.03.005.PubMedView Article
- Yoder JA, Haire RN, Litman GW: Cloning of two zebrafish cDNAs that share domains with MHC class II-associated invariant chain. Immunogenetics. 1999, 50: 84-88. 10.1007/s002510050691.PubMedView Article
- Fujiki K, Smith CM, Liu L, Sundick RS, Dixon B: Alternate forms of MHC class II-associated invariant chain are not produced by alternative splicing in rainbow trout (Oncorhynchus mykiss) but are encoded by separate genes. Dev Comp Immunol. 2003, 27: 377-391. 10.1016/S0145-305X(02)00119-2.PubMedView Article
- Dijkstra JM, Kiryu I, Kollner B, Yoshiura Y, Ototake M: MHC class II invariant chain homologues in rainbow trout (Oncorhynchus mykiss). Fish Shellfish Immunol. 2003, 15: 91-105. 10.1016/S1050-4648(02)00141-9.PubMedView Article
- 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 Article
- 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 Article
- Stet RJ, Kruiswijk CP, Dixon B: Major histocompatibility lineages and immune gene function in teleost fishes: the road not taken. Crit Rev Immunol. 2003, 23: 441-471.PubMedView Article
- Star B, Nederbragt AJ, Jentoft S, Grimholt U, Malmstrøm M, Gregers TF, Rounge TB, Paulsen J, Solbakken MH, Sharma A, Wetten OF, Lanzén A, Winer R, Knight J, Vogel JH, Aken B, Andersen O, Lagesen K, Tooming-Klunderud A, Edvardsen RB, Tina KG, Espelund M, Nepal C, Previti C, Karlsen BO, Moum T, Skage M, Berg PR, Gjøen T, Kuhl H, et al: The genome sequence of Atlantic cod reveals a unique immune system. Nature. 2011, 477: 207-210. 10.1038/nature10342.PubMed CentralPubMedView Article
- McConnell TJ, Godwin UB, Cuthbertson BJ: Expressed major histocompatibility complex class II loci in fishes. Immunol Rev. 1998, 166: 294-300. 10.1111/j.1600-065X.1998.tb01270.x.PubMedView Article
- Harstad H, Lukacs MF, Bakke HG, Grimholt U: Multiple expressed MHC class II loci in salmonids; details of one non-classical region in Atlantic salmon (Salmo salar). BMC Genomics. 2008, 9: 193-10.1186/1471-2164-9-193.PubMed CentralPubMedView Article
- Sato A, Dongak R, Hao L, Shintani S, Sato T: Organization of Mhc class II A and B genes in the tilapiine fish Oreochromis. Immunogenetics. 2012, 64: 679-690. 10.1007/s00251-012-0618-0.PubMedView Article
- Pos W, Sethi DK, Call MJ, Schulze MS, Anders AK, Pyrdol J, Wucherpfennig KW: Crystal structure of the HLA-DM-HLA-DR1 complex defines mechanisms for rapid peptide selection. Cell. 2012, 151: 1557-1568. 10.1016/j.cell.2012.11.025.PubMed CentralPubMedView Article
- Palti Y, Rodriguez MF, Gahr SA, Hansen JD: Evolutionary history of the ABCB2 genomic region in teleosts. Dev Comp Immunol. 2007, 31: 483-498. 10.1016/j.dci.2006.07.010.PubMedView Article
- Ohashi K, Takizawa F, Tokumaru N, Nakayasu C, Toda H, Fischer U, Moritomo T, Hashimoto K, Nakanishi T, Dijkstra JM: A molecule in teleost fish, related with human MHC-encoded G6F, has a cytoplasmic tail with ITAM and marks the surface of thrombocytes and in some fishes also of erythrocytes. Immunogenetics. 2010, 62: 543-559. 10.1007/s00251-010-0460-1.PubMedView Article
- Lukacs MF, Harstad H, Bakke HG, Beetz-Sargent M, McKinnel L, Lubieniecki KP, Koop BF, Grimholt U: Comprehensive analysis of MHC class I genes from the U-, S-, and Z-lineages in Atlantic salmon. BMC Genomics. 2010, 11: 154-10.1186/1471-2164-11-154.PubMed CentralPubMedView Article
- Ohta Y, Okamura K, McKinney EC, Bartl S, Hashimoto K, Flajnik MF: Primitive synteny of vertebrate major histocompatibility complex class I and class II genes. Proc Natl Acad Sci USA. 2000, 97: 4712-4717. 10.1073/pnas.97.9.4712.PubMed CentralPubMedView Article
- 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 Article
- Cadavid LF, Mejia BE, Watkins DI: MHC class I genes in a New World primate, the cotton-top tamarin (Saguinus oedipus), have evolved by an active process of loci turnover. Immunogenetics. 1999, 49: 196-205. 10.1007/s002510050480.PubMedView Article
- Fremont DH, Monnaie D, Nelson CA, Hendrickson WA, Unanue ER: Crystal structure of I-Ak in complex with a dominant epitope of lysozyme. Immunity. 1998, 8: 305-317. 10.1016/S1074-7613(00)80536-1.PubMedView Article
- Painter CA, Stern LJ: Conformational variation in structures of classical and non-classical MHCII proteins and functional implications. Immunol Rev. 2012, 250: 144-157. 10.1111/imr.12003.PubMed CentralPubMedView Article
- Zhou Z, Callaway KA, Weber DA, Jensen PE: Cutting Edge: HLA-DM functions through a mechanism that does not require specific conserved hydrogen bonds in class II MHC-peptide complexes. J Immunol. 2009, 183: 4187-4191. 10.4049/jimmunol.0901663.PubMed CentralPubMedView Article
- Salomonsen J, Marston D, Avila D, Bumstead N, Johansson B, Juul-Madsen H, Olesen GD, Riegert P, Skjødt K, Vainio O, Wiles MV, Kaufman J: The properties of the single chicken MHC classical class II α chain (B-LA) gene indicate an ancient origin for the DR/E-like isotype of class II molecules. Immunogenetics. 2003, 55: 605-614. 10.1007/s00251-003-0620-7.PubMedView Article
- Wang XX, Li Y, Yin Y, Mo M, Wang Q, Gao W, Wang L, Mariuzza RA: Affinity maturation of human CD4 by yeast surface display and crystal structure of a CD4-HLA-DR1 complex. Proc Natl Acad Sci USA. 2011, 108: 15960-15965. 10.1073/pnas.1109438108.PubMed CentralPubMedView Article
- Lindstedt R, Liljedahl M, Peleraux A, Peterson PA, Karlsson L: The MHC class II molecule H2-M is targeted to an endosomal compartment by a tyrosine-based targeting motif. Immunity. 1995, 3: 561-572. 10.1016/1074-7613(95)90127-2.PubMedView Article
- Marks MS, Roche PA, van Donselaar E, Woodruff L, Peters PJ, Bonifacino JS: A lysosomal targeting signal in the cytoplasmic tail of the β chain directs HLA-DM to MHC class II compartments. J Cell Biol. 1995, 131: 351-369. 10.1083/jcb.131.2.351.PubMedView Article
- Hughes AL, Nei M: Evolution of the major histocompatibility complex: independent origin of nonclassical class I genes in different groups of mammals. Mol Biol Evol. 1989, 6: 559-579.PubMed
- Doebele RC, Busch R, Scott HM, Pashine A, Mellins ED: Determination of the HLA-DM interaction site on HLA-DR molecules. Immunity. 2000, 13: 517-527. 10.1016/S1074-7613(00)00051-0.PubMedView Article
- Pashine A, Busch R, Belmares MP, Munning JN, Doebele RC, Buckingham M, Nolan GP, Mellins ED: Interaction of HLA-DR with an acidic face of HLA-DM disrupts sequence-dependent interactions with peptides. Immunity. 2003, 19: 183-192. 10.1016/S1074-7613(03)00200-0.PubMedView Article
- Anders A-K, Call MJ, Schulze M-S ED, Fowler KD, Schubert DA, Seth NP, Sundberg EJ, Wucherpfennig KW: HLA-DM captures partially empty HLA-DR molecules for catalyzed removal of peptide. Nat Immunol. 2011, 12: 54-61. 10.1038/ni.1967.PubMed CentralPubMedView Article
- Painter CA, Negroni MP, Kellersberger KA, Zavala-Ruiz Z, Evans JE, Stern LJ: Conformational lability in the class II MHC 310 helix and adjacent extended strand dictate HLA-DM susceptibility and peptide exchange. Proc Natl Acad Sci USA. 2011, 108: 19329-19334. 10.1073/pnas.1108074108.PubMed CentralPubMedView Article
- Stratikos E, Mosyak L, Zaller DM, Wiley DC: Identification of the lateral interaction surfaces of human histocompatibility leukocyte antigen (HLA)-DM with HLA-DR1 by formation of tethered complexes that present enhanced HLA-DM catalysis. J Exp Med. 2002, 196: 173-183. 10.1084/jem.20020117.PubMed CentralPubMedView Article
- Stern LJ, Brown JH, Jardetzky TS, Gorga JC, Urban RG, Strominger JL, Wiley DC: Crystal structure of the human class II MHC protein HLA-DR1 complexed with an influenza virus peptide. Nature. 1994, 368: 215-221. 10.1038/368215a0.PubMedView Article
- Guce AI, Mortimer SE, Yoon T, Painter CA, Jiang W, Mellins ED, Stern LJ: HLA-DO acts as a substrate mimic to inhibit HLA-DM by a competitive mechanism. Nature Struct Mol Biol. 2013, 20: 90-98.View Article
- Stet RJM, De Vries B, Mudde K, Hermsen T, Van Heerwaarden J, Shum BP, Grimholt U: Unique haplotypes of co-segregating major histocompatibility class II A and class II B alleles in Atlantic salmon (Salmo salar) give rise to diverse class II genotypes. Immunogenetics. 2002, 54: 320-331. 10.1007/s00251-002-0477-1.PubMedView Article
- Murphy K: Janeway’s Immunobiology. 2012, New York: Garland Science, 8
- Cordero-Morales JF, Jogini V, Chakrapani S, Perozo E: A multipoint Hydrogen-bond network underlying KcsA C-type inactivation. Biophys J. 2011, 100: 2387-2393. 10.1016/j.bpj.2011.01.073.PubMed CentralPubMedView Article
- Dijkstra JM, Katagiri T, Hosomichi K, Yanagiya K, Inoko H, Ototake M, Aoki T, Hashimoto K, Shiina T: A third broad lineage of major histocompatibility complex (MHC) class I in teleost fish; MHC class II linkage and processed genes. Immunogenetics. 2007, 59: 305-321. 10.1007/s00251-007-0198-6.PubMedView Article
- Hausmann DHF, Yu B, Hausmann S, Wucherpfennig KW: pH dependent peptide binding properties of the type I diabetes-associated I-Ag7 molecule: rapid release of CLIP at an endosomal pH. J Exp Med. 1999, 189: 1723-1733. 10.1084/jem.189.11.1723.PubMed CentralPubMedView Article
- Schafer PH, Pierce SK, Jardetzky TS: The structure of MHC class II: a role for dimer of dimers. Seminars Immunol. 1995, 7: 389-398. 10.1006/smim.1995.0043.View Article
- Burge C, Karlin S: Prediction of complete gene structures in human genomic DNA. J Mol Biol. 1997, 268: 78-94. 10.1006/jmbi.1997.0951.PubMedView Article
- Salamov AA, Solovyev VV: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10: 516-522. 10.1101/gr.10.4.516.PubMed CentralPubMedView Article
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstucting phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMed
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView Article
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.